Showing 1–50 of 138 results for author: Kennedy, D

ProACT: An Augmented Reality Testbed for Intelligent Prosthetic Arms

Authors: Shivani Guptasarma, Monroe D. Kennedy III

Abstract: Upper-limb amputees face tremendous difficulty in operating dexterous powered prostheses. Previous work has shown that aspects of prosthetic hand, wrist, or elbow control can be improved through "intelligent" control, by combining movement-based or gaze-based intent estimation with low-level robotic autonomy. However, no such solutions exist for whole-arm control. Moreover, hardware platforms for… ▽ More Upper-limb amputees face tremendous difficulty in operating dexterous powered prostheses. Previous work has shown that aspects of prosthetic hand, wrist, or elbow control can be improved through "intelligent" control, by combining movement-based or gaze-based intent estimation with low-level robotic autonomy. However, no such solutions exist for whole-arm control. Moreover, hardware platforms for advanced prosthetic control are expensive, and existing simulation platforms are not well-designed for integration with robotics software frameworks. We present the Prosthetic Arm Control Testbed (ProACT), a platform for evaluating intelligent control methods for prosthetic arms in an immersive (Augmented Reality) simulation setting. Using ProACT with non-amputee participants, we compare performance in a Box-and-Blocks Task using a virtual myoelectric prosthetic arm, with and without intent estimation. Our results show that methods using intent estimation improve both user satisfaction and the degree of success in the task. To the best of our knowledge, this constitutes the first study of semi-autonomous control for complex whole-arm prostheses, the first study including sequential task modeling in the context of wearable prosthetic arms, and the first testbed of its kind. Towards the goal of supporting future research in intelligent prosthetics, the system is built upon on existing open-source frameworks for robotics. △ Less

Submitted 6 July, 2024; originally announced July 2024.

Comments: 10 pages, 8 figures. Under review. Code and data will be available at https://armlabstanford.github.io/proact after publication

Suppression of temperature-gradient-driven turbulence by sheared flows in fusion plasmas

Authors: P. G. Ivanov, T. Adkins, D. Kennedy, M. Giacomin, M. Barnes, A. A. Schekochihin

Abstract: Starting from the assumption that saturation of plasma turbulence driven by temperature-gradient instabilities in fusion plasmas is achieved by a local energy cascade between a long-wavelength outer scale, where energy is injected into the fluctuations, and a small-wavelength dissipation scale, where fluctuation energy is thermalized by particle collisions, we formulate a detailed phenomenological… ▽ More Starting from the assumption that saturation of plasma turbulence driven by temperature-gradient instabilities in fusion plasmas is achieved by a local energy cascade between a long-wavelength outer scale, where energy is injected into the fluctuations, and a small-wavelength dissipation scale, where fluctuation energy is thermalized by particle collisions, we formulate a detailed phenomenological theory for the influence of perpendicular flow shear on magnetized-plasma turbulence. Our theory introduces two distinct regimes, called the weak-shear and strong-shear regimes, each with its own set of scaling laws for the scale and amplitude of the fluctuations and for the level of turbulent heat transport. We discover that the ratio of the typical radial and poloidal wavenumbers of the fluctuations (i.e., their aspect ratio) at the outer scale plays a central role in determining the dependence of the turbulent transport on the imposed flow shear. Our theoretical predictions are found to be in excellent agreement with numerical simulations of two paradigmatic models of fusion-relevant plasma turbulence: (i) an electrostatic fluid model of slab electron-scale turbulence, and (ii) Cyclone-base-case gyrokinetic ion-scale turbulence. Additionally, our theory envisions a potential mechanism for the suppression of electron-scale turbulence by perpendicular ion-scale flows based on the role of the aforementioned aspect ratio of the electron-scale fluctuations. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 21 pages, 13 figures, 1 table

A quasi-linear model of electromagnetic turbulent transport and its application to flux-driven transport predictions for STEP

Authors: M. Giacomin, D. Dickinson, W. Dorland, N. R. Mandell, A. Bokshi, F. J. Casson, H. G. Dudding, D. Kennedy, B. S. Patel, C. M. Roach

Abstract: A quasi-linear reduced transport model is developed from a database of high-$β$ electromagnetic nonlinear gyrokinetic simulations performed with Spherical Tokamak for Energy Production (STEP) relevant parameters. The quasi-linear model is fully electromagnetic and accounts for the effect of equilibrium flow shear using a novel approach. Its flux predictions are shown to agree quantitatively with p… ▽ More A quasi-linear reduced transport model is developed from a database of high-$β$ electromagnetic nonlinear gyrokinetic simulations performed with Spherical Tokamak for Energy Production (STEP) relevant parameters. The quasi-linear model is fully electromagnetic and accounts for the effect of equilibrium flow shear using a novel approach. Its flux predictions are shown to agree quantitatively with predictions from local nonlinear gyrokinetic simulations across a broad range of STEP-relevant local equilibria. This reduced transport model is implemented in the T3D transport solver that is used to perform the first flux-driven simulations for STEP to account for transport from hybrid-KBM turbulence, which dominates over a wide region of the core plasma. Nonlinear gyrokinetic simulations of the final transport steady state from T3D return turbulent fluxes that are consistent with the reduced model, indicating that the quasi-linear model may also be appropriate for describing the transport steady state. Within the assumption considered here, our simulations support the existence of a transport steady state in STEP with a fusion power comparable to that in the burning flat-top of the conceptual design, but do not demonstrate how this state can be accessed. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Hybrid data assimilation techniques using the adjoint method in a coupled Lorenz system

Authors: Philip David Kennedy, Abhirup Banerjee, Armin Köhl, Detlef Stammer

Abstract: A hybrid 4D-variational data assimilation method for numerical climate models is introduced using the Lorenz '63 model. This new approach has the potential to optimise a high complexity Earth system model (ESM) by utilising the adjoint equations of an intermediate complexity ESM. The method is conceptually demonstrated by consecutively synchronising two Lorenz '63 systems to observations before op… ▽ More A hybrid 4D-variational data assimilation method for numerical climate models is introduced using the Lorenz '63 model. This new approach has the potential to optimise a high complexity Earth system model (ESM) by utilising the adjoint equations of an intermediate complexity ESM. The method is conceptually demonstrated by consecutively synchronising two Lorenz '63 systems to observations before optimisation. The first represents a 'high complexity' model and the second an 'intermediate complexity' model which has adjoint equations. This method will save computational power for a full ESM and has negligible error and uncertainty change compared to the optimisation of a single model with adjoint equations. A similar setup can be applied to sparse observations. An alternative assimilation setup, with two identical models, is used to filter noisy data. This reduces optimised parametric model uncertainty by approximately one third. Such a precision gain could prove valuable for seasonal, annual, and decadal predictions. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 13 pages, 9 figures

On the importance of parallel magnetic-field fluctuations for electromagnetic instabilities in STEP

Authors: D. Kennedy, C. M. Roach, M. Giacomin, P. Ivanov, T. Adkins, F. Sheffield, T. G örler, A. Bokshi, D. Dickinson, H. G. Dudding, B. S. Patel

Abstract: [ABRIDGED] This paper discusses the importance of parallel perturbations of the magnetic-field in gyrokinetic simulations of electromagnetic instabilities and turbulence at mid-radius in the burning plasma phase of the conceptual high-$β$, reactor-scale, tight-aspect-ratio tokamak STEP. Previous studies have revealed the presence of unstable hybrid kinetic ballooning modes (hKBMs) at binormal scal… ▽ More [ABRIDGED] This paper discusses the importance of parallel perturbations of the magnetic-field in gyrokinetic simulations of electromagnetic instabilities and turbulence at mid-radius in the burning plasma phase of the conceptual high-$β$, reactor-scale, tight-aspect-ratio tokamak STEP. Previous studies have revealed the presence of unstable hybrid kinetic ballooning modes (hKBMs) at binormal scales approaching the ion Larmor radius. In this STEP plasma it was found that the hKBM requires the inclusion of parallel magnetic-field perturbations to be linearly unstable. Here, the extent to which the inclusion of fluctuations in the parallel magnetic-field can be relaxed is explored through gyrokinetic simulations. In particular, the frequently used MHD approximation (dropping $δ\! B_{\parallel}$ and setting the $\nabla B$ drift frequency equal to the curvature drift frequency) is discussed and simulations explore whether this approximation is useful for modelling STEP plasmas. It is shown that the MHD approximation can reproduce some of the linear properties of the full STEP gyrokinetic system, but is too stable at low $k_y$ and nonlinear simulations using the MHD approximation result in very different transport states. It is demonstrated that the MHD approximation is challenged by the high $β^{\prime}$ values in STEP, and that the approximation improves considerably at lower $β^{\prime}$. Furthermore, it is shown that the sensitivity of STEP to $δ\! B_{\parallel}$ fluctuations is primarily because the plasma sits close to marginality and it is shown that in slightly more strongly driven conditions the hKBM is unstable without $δ\! B_{\parallel}.$ Crucially, it is demonstrated that the state of large transport typically predicted by local electromagnetic gyrokinetic simulations of STEP plasmas is not solely due to $δ\! B_{\parallel}$ physics. △ Less

Submitted 20 May, 2024; v1 submitted 16 February, 2024; originally announced February 2024.

Comments: 29th IAEA FEC conference

Bayesian Optimization Algorithms for Accelerator Physics

Authors: Ryan Roussel, Auralee L. Edelen, Tobias Boltz, Dylan Kennedy, Zhe Zhang, Fuhao Ji, Xiaobiao Huang, Daniel Ratner, Andrea Santamaria Garcia, Chenran Xu, Jan Kaiser, Angel Ferran Pousa, Annika Eichler, Jannis O. Lubsen, Natalie M. Isenberg, Yuan Gao, Nikita Kuklev, Jose Martinez, Brahim Mustapha, Verena Kain, Weijian Lin, Simone Maria Liuzzo, Jason St. John, Matthew J. V. Streeter, Remi Lehe , et al. (1 additional authors not shown)

Abstract: Accelerator physics relies on numerical algorithms to solve optimization problems in online accelerator control and tasks such as experimental design and model calibration in simulations. The effectiveness of optimization algorithms in discovering ideal solutions for complex challenges with limited resources often determines the problem complexity these methods can address. The accelerator physics… ▽ More Accelerator physics relies on numerical algorithms to solve optimization problems in online accelerator control and tasks such as experimental design and model calibration in simulations. The effectiveness of optimization algorithms in discovering ideal solutions for complex challenges with limited resources often determines the problem complexity these methods can address. The accelerator physics community has recognized the advantages of Bayesian optimization algorithms, which leverage statistical surrogate models of objective functions to effectively address complex optimization challenges, especially in the presence of noise during accelerator operation and in resource-intensive physics simulations. In this review article, we offer a conceptual overview of applying Bayesian optimization techniques towards solving optimization problems in accelerator physics. We begin by providing a straightforward explanation of the essential components that make up Bayesian optimization techniques. We then give an overview of current and previous work applying and modifying these techniques to solve accelerator physics challenges. Finally, we explore practical implementation strategies for Bayesian optimization algorithms to maximize their performance, enabling users to effectively address complex optimization challenges in real-time beam control and accelerator design. △ Less

Submitted 5 April, 2024; v1 submitted 9 December, 2023; originally announced December 2023.

Data downloaded via parachute from a NASA super-pressure balloon

Authors: Ellen L. Sirks, Richard Massey, Ajay S. Gill, Jason Anderson, Steven J. Benton, Anthony M. Brown, Paul Clark, Joshua English, Spencer W. Everett, Aurelien A. Fraisse, Hugo Franco, John W. Hartley, David Harvey, Bradley Holder, Andrew Hunter, Eric M. Huff, Andrew Hynous, Mathilde Jauzac, William C. Jones, Nikky Joyce, Duncan Kennedy, David Lagattuta, Jason S. -Y. Leung, Lun Li, Stephen Lishman , et al. (18 additional authors not shown)

Abstract: In April to May 2023, the superBIT telescope was lifted to the Earth's stratosphere by a helium-filled super-pressure balloon, to acquire astronomical imaging from above (99.5% of) the Earth's atmosphere. It was launched from New Zealand then, for 40 days, circumnavigated the globe five times at a latitude 40 to 50 degrees South. Attached to the telescope were four 'DRS' (Data Recovery System) cap… ▽ More In April to May 2023, the superBIT telescope was lifted to the Earth's stratosphere by a helium-filled super-pressure balloon, to acquire astronomical imaging from above (99.5% of) the Earth's atmosphere. It was launched from New Zealand then, for 40 days, circumnavigated the globe five times at a latitude 40 to 50 degrees South. Attached to the telescope were four 'DRS' (Data Recovery System) capsules containing 5 TB solid state data storage, plus a GNSS receiver, Iridium transmitter, and parachute. Data from the telescope were copied to these, and two were dropped over Argentina. They drifted 61 km horizontally while they descended 32 km, but we predicted their descent vectors within 2.4 km: in this location, the discrepancy appears irreducible below 2 km because of high speed, gusty winds and local topography. The capsules then reported their own locations to within a few metres. We recovered the capsules and successfully retrieved all of superBIT's data - despite the telescope itself being later destroyed on landing. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: 12 pages

Journal ref: Aerospace 2023, 10, 960

Gaussian Process Regression models for the properties of micro-tearing modes in spherical tokamak

Authors: William Hornsby, Ander Gray, James Buchanan, Bhavin Patel, Daniel Kennedy, Francis Casson, Colin Roach, Mikkel Lykkegaard, Huy Nguyen, Nikolaos Papadimas, Ben Fourcin, Jordan Hart

Abstract: Spherical tokamaks (STs) have many desirable features that make them an attractive choice for a future fusion power plant. Power plant viability is intrinsically related to plasma heat and particle confinement and this is often determined by the level of micro-instability driven turbulence. Accurate calculation of the properties of turbulent micro-instabilities is therefore critical for tokamak de… ▽ More Spherical tokamaks (STs) have many desirable features that make them an attractive choice for a future fusion power plant. Power plant viability is intrinsically related to plasma heat and particle confinement and this is often determined by the level of micro-instability driven turbulence. Accurate calculation of the properties of turbulent micro-instabilities is therefore critical for tokamak design, however, the evaluation of these properties is computationally expensive. The considerable number of geometric and thermodynamic parameters and the high resolutions required to accurately resolve these instabilities makes repeated use of direct numerical simulations in integrated modelling workflows extremely computationally challenging and creates the need for fast, accurate, reduced-order models. This paper outlines the development of a data-driven reduced-order model, often termed a {\it surrogate model} for the properties of micro-tearing modes (MTMs) across a spherical tokamak reactor-relevant parameter space utilising Gaussian Process Regression (GPR) and classification; techniques from machine learning. These two components are used in an active learning loop to maximise the efficiency of data acquisition thus minimising computational cost. The high-fidelity gyrokinetic code GS2 is used to calculate the linear properties of the MTMs: the mode growth rate, frequency and normalised electron heat flux; core components of a quasi-linear transport model. Five-fold cross-validation and direct validation on unseen data is used to ascertain the performance of the resulting surrogate models. △ Less

Submitted 14 November, 2023; v1 submitted 18 September, 2023; originally announced September 2023.

Demonstration of Autonomous Emittance Characterization at the Argonne Wakefield Accelerator

Authors: Ryan Roussel, Auralee Edelen, Dylan Kennedy, Seongyeol Kim, Eric Wisniewski, John Power

Abstract: Transverse beam emittance plays a key role in the performance of high brightness accelerators. Characterizing beam emittance is often done using a quadrupole scan, which fits beam matrix elements to experimental measurements using first order optics. Despite its simplicity at face value, this procedure is difficult to automate due to practical limitations. Key issues that must be addressed include… ▽ More Transverse beam emittance plays a key role in the performance of high brightness accelerators. Characterizing beam emittance is often done using a quadrupole scan, which fits beam matrix elements to experimental measurements using first order optics. Despite its simplicity at face value, this procedure is difficult to automate due to practical limitations. Key issues that must be addressed include maintaining beam size measurement validity by keeping beams within the radius to diagnostic screens, ensuring that measurement fitting produces physically valid results, and accurately characterizing emittance uncertainty. We describe a demonstration of the Bayesian Exploration technique towards solving this problem at the Argonne Wakefield Accelerator, enabling a turnkey, autonomous quadrupole scan tool that can be used to quickly measure beam emittances at various locations in accelerators with limited operator input. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Electromagnetic gyrokinetic instabilities in STEP

Authors: Daniel Kennedy, Maurizio Giacomin, Francis J Casson, David Dickinson, William A Hornsby, Bhavin S Patel, Colin M Roach

Abstract: We present herein the results of a linear gyrokinetic analysis of electromagnetic microinstabilites in the conceptual high-$β$, reactor-scale, tight-aspect-ratio tokamak STEP (Spherical Tokamak for Energy Production, https://step.ukaea.uk). We examine a range of flux surfaces between the deep core and the pedestal top for two candidate flat-top operating points of the prototype device. Local linea… ▽ More We present herein the results of a linear gyrokinetic analysis of electromagnetic microinstabilites in the conceptual high-$β$, reactor-scale, tight-aspect-ratio tokamak STEP (Spherical Tokamak for Energy Production, https://step.ukaea.uk). We examine a range of flux surfaces between the deep core and the pedestal top for two candidate flat-top operating points of the prototype device. Local linear gyrokinetic analysis is performed to determine the type of microinstabilities that arise under these reactor-relevant conditions. We find that the equilibria are dominated at ion binormal scales by a hybrid version of the Kinetic Ballooning Mode (KBM) instability that has significant linear drive contributions from the ion temperature gradient and from trapped electrons, while collisional Microtearing Modes (MTMs) are sub-dominantly also unstable at similar binormal scales. The hybrid-KBM and MTM exhibit very different radial scales. We study the sensitivity of these instabilities to physics parameters, and discuss potential mechanisms for mitigating them. The results of this investigation are compared to a small set of similar conceptual reactor designs in the literature. A detailed benchmark of the linear results is performed using three gyrokinetic codes; alongside extensive resolution testing and sensitivity to numerical parameters providing confidence in the results of our calculations, and paving the way for detailed nonlinear studies in a companion article. △ Less

Submitted 5 October, 2023; v1 submitted 4 July, 2023; originally announced July 2023.

Journal ref: Nucl. Fusion 63 126061 (2023)

On electromagnetic turbulence and transport in STEP

Authors: Maurizio Giacomin, Daniel Kennedy, Francis J Casson, Ajay C. J., David Dickinson, Bhavin S. Patel, Colin M. Roach

Abstract: In this work, we present first-of-their-kind nonlinear local gyrokinetic simulations of electromagnetic turbulence at mid-radius in the burning plasma phase of the conceptual high-$β$, reactor-scale, tight-aspect-ratio tokamak STEP (Spherical Tokamak for Energy Production). A prior linear analysis in D. Kennedy et al. 2023 Nucl. Fusion 63 126061 reveals the presence of unstable hybrid kinetic ball… ▽ More In this work, we present first-of-their-kind nonlinear local gyrokinetic simulations of electromagnetic turbulence at mid-radius in the burning plasma phase of the conceptual high-$β$, reactor-scale, tight-aspect-ratio tokamak STEP (Spherical Tokamak for Energy Production). A prior linear analysis in D. Kennedy et al. 2023 Nucl. Fusion 63 126061 reveals the presence of unstable hybrid kinetic ballooning modes, where inclusion of the compressional magnetic field fluctuation, $δB_{\parallel}$, is crucial, and subdominant microtearing modes are found at binormal scales approaching the ion-Larmor radius. Local nonlinear gyrokinetic simulations on the selected surface in the central core region suggest that hybrid kinetic ballooning modes can drive large turbulent transport, and that there is negligible turbulent transport from subdominant microtearing modes when hybrid kinetic ballooning modes are artificially suppressed (through the omission of $δB_{\parallel}$). Nonlinear simulations that include perpendicular equilibrium flow shear can saturate at lower fluxes that are more consistent with the available sources in STEP. This analysis suggests that hybrid kinetic ballooning modes could play an important role in setting the turbulent transport in STEP, and possible mechanisms to mitigate turbulent transport are discussed. Increasing the safety factor or the pressure gradient strongly reduces turbulent transport from hybrid kinetic ballooning modes in the cases considered here. Challenges of simulating electromagnetic turbulence in this high-$β$ regime are highlighted. In particular the observation of radially extended turbulent structures in the absence of equilibrium flow shear motivates future advanced global gyrokinetic simulations that include $δB_\parallel$. △ Less

Submitted 24 March, 2024; v1 submitted 4 July, 2023; originally announced July 2023.

Journal ref: Plasma Phys. Control. Fusion 2024

Nonlinear microtearing modes in MAST and their stochastic layer formation

Authors: M. Giacomin, D. Dickinson, D. Kennedy, B. S. Patel, C. M. Roach

Abstract: First nonlinear gyrokinetic simulations of microtearing modes in the core of a MAST case are performed on two surfaces of the high-collisionality discharge used in Valovič et al. Nucl. Fusion 51.7 (2011) to obtain the favorable energy confinement scaling with collisionality, $τ_E\propto\,ν_*^{-1}$. On the considered surfaces microtearing modes dominate linearly at binormal length scales of the ord… ▽ More First nonlinear gyrokinetic simulations of microtearing modes in the core of a MAST case are performed on two surfaces of the high-collisionality discharge used in Valovič et al. Nucl. Fusion 51.7 (2011) to obtain the favorable energy confinement scaling with collisionality, $τ_E\propto\,ν_*^{-1}$. On the considered surfaces microtearing modes dominate linearly at binormal length scales of the order of the ion Larmor radius. While the effect of electron collision frequency is moderate in linear simulations, a strong dependence on this parameter is found in nonlinear simulations at $r/a=0.5$, where $r$ and $a$ are the surface and tokamak minor radius, respectively. The dynamics of magnetic islands generated by microtearing modes is analysed, showing that the radial extent of the stochastic region caused by islands overlapping plays an important role in determining the saturation level of the microtearing mode driven heat flux. Local nonlinear gyrokinetic simulations show that the microtearing mode driven heat flux, $Q_e^\mathrm{MTM}$, is largely dominated by magnetic flutter and depends strongly on the magnetic shear, $\hat{s}$. Comparing two surfaces, $r/a=0.5$ and $r/a=0.6$, reveals that $Q_e^\mathrm{MTM}$ is negligible at $r/a=0.5$ ($\hat{s}=0.34$), with the electron temperature gradient driven heat flux, $Q_e^\mathrm{ETG}$, comparable to the experimental electron heat flux, $Q_e^\mathrm{exp}$, while $Q_e^\mathrm{MTM}$ is significantly larger and comparable to $Q_e^\mathrm{ETG}$ and $Q_e^\mathrm{exp}$ at $r/a=0.6$ ($\hat{s}=1.1$). Microtearing modes cause more experimentally significant transport in higher $\hat{s}$ regions and may influence (together with electron temperature gradient modes) the observed scaling of energy confinement time with collisionality (Valovič et al. Nucl. Fusion 51.7 (2011)). △ Less

Submitted 30 August, 2023; v1 submitted 4 March, 2023; originally announced March 2023.

Journal ref: Plasma Phys. Control. Fusion 65 095019 (2023)

Nebular-Phase Spectra of Type Ia Supernovae from the Las Cumbres Observatory Global Supernova Project

Authors: M. L. Graham, T. D. Kennedy, S. Kumar, R. C. Amaro, D. J. Sand, S. W. Jha, L. Galbany, J. Vinko, J. C. Wheeler, E. Y. Hsiao, K. A. Bostroem, J. Burke, D. Hiramatsu, G. Hosseinzadeh, C. McCully, D. A. Howell, T. Diamond, P. Hoeflich, X. Wang, W. Li

Abstract: The observed diversity in Type Ia supernovae (SNe Ia) -- the thermonuclear explosions of carbon-oxygen white dwarf stars used as cosmological standard candles -- is currently met with a variety of explosion models and progenitor scenarios. To help improve our understanding of whether and how often different models contribute to the occurrence of SNe Ia and their assorted properties, we present a c… ▽ More The observed diversity in Type Ia supernovae (SNe Ia) -- the thermonuclear explosions of carbon-oxygen white dwarf stars used as cosmological standard candles -- is currently met with a variety of explosion models and progenitor scenarios. To help improve our understanding of whether and how often different models contribute to the occurrence of SNe Ia and their assorted properties, we present a comprehensive analysis of seven nearby SNe Ia. We obtained one to two epochs of optical spectra with Gemini Observatory during the nebular phase ($>$200 days past peak) for each of these events, all of which had time-series of photometry and spectroscopy at early times (the first $\sim$8 weeks after explosion). We use the combination of early- and late-time observations to assess the predictions of various models for the explosion (e.g., double-detonation, off-center detonation, stellar collisions), progenitor star (e.g., ejecta mass, metallicity), and binary companion (e.g., another white dwarf or a non-degenerate star). Overall, we find general consistency in our observations with spherically-symmetric models for SN Ia explosions, and with scenarios in which the binary companion is another degenerate star. We also present an in-depth analysis of SN 2017fzw, a member of the sub-group of SNe Ia which appear to be transitional between the subluminous "91bg-like" events and normal SNe Ia, and for which nebular-phase spectra are rare. △ Less

Submitted 19 January, 2022; originally announced January 2022.

Comments: 26 pages, 14 figures, 8 tables, accepted to MNRAS

Recommendations for repositories and scientific gateways from a neuroscience perspective

Authors: Malin Sandström, Mathew Abrams, Jan Bjaalie, Mona Hicks, David Kennedy, Arvind Kumar, JB Poline, Prasun Roy, Paul Tiesinga, Thomas Wachtler, Wojtek Goscinski

Abstract: Digital services such as repositories and science gateways have become key resources for the neuroscience community, but users often have a hard time orienting themselves in the service landscape to find the best fit for their particular needs. INCF (International Neuroinformatics Coordinating Facility) has developed a set of recommendations and associated criteria for choosing or setting up and r… ▽ More Digital services such as repositories and science gateways have become key resources for the neuroscience community, but users often have a hard time orienting themselves in the service landscape to find the best fit for their particular needs. INCF (International Neuroinformatics Coordinating Facility) has developed a set of recommendations and associated criteria for choosing or setting up and running a repository or scientific gateway, intended for the neuroscience community, with a FAIR neuroscience perspective. These recommendations have neurosciences as their primary use case but are often general. Considering the perspectives of researchers and providers of repositories as well as scientific gateways, the recommendations harmonize and complement existing work on criteria for repositories and best practices. The recommendations cover a range of important areas including accessibility, licensing, community responsibility and technical and financial sustainability of a service. △ Less

Submitted 3 January, 2022; originally announced January 2022.

Comments: 10 pages, submitted to Scientific Data

Perturbative heavy quark contributions to the anomalous magnetic moment of the muon

Authors: P. D. Kennedy, J. Erler, H. Spiesberger

Abstract: We discuss a method for calculating the heavy quark vacuum polarisation contribution to the muon anomalous magnetic moment, $a_μ$, using perturbative QCD. This approach is independent of $e^+e^-$ cross-section data allowing a fully theoretical evaluation of these contributions. We confirm an existing result at lower orders in $α_s$ and state a new explicit analytic formula which includes terms up… ▽ More We discuss a method for calculating the heavy quark vacuum polarisation contribution to the muon anomalous magnetic moment, $a_μ$, using perturbative QCD. This approach is independent of $e^+e^-$ cross-section data allowing a fully theoretical evaluation of these contributions. We confirm an existing result at lower orders in $α_s$ and state a new explicit analytic formula which includes terms up to $\mathcal{O}\left(α_s^3\right)$. Numerically the charm quark contribution to $a_μ$ is found to be $a_μ^c = (14.5 \pm 0.2)\times 10^{-10}$ and the bottom contributes $a_μ^b = (0.302 \pm 0.002) \times 10^{-10} $. Our uncertainty estimates include both parametric uncertainties from $\hat{m}_q(\hat{m}_q)$ and $α_s(\hat{m}_q)$, and theoretical uncertainties in the perturbative expansion. Comparison is made between these results and alternative approaches such as lattice QCD or those based on a dispersion relation and cross-section data. △ Less

Submitted 25 November, 2021; originally announced November 2021.

Comments: 8 pages, 2 figures, contribution to 16th International Workshop on Tau Lepton Physics (TAU2021), September 27 - October 1, 2021

Report number: MITP/21-064

The evolutionary epidemiology of pathogens during vaccination campaigns

Authors: Troy Day, David A. Kennedy, Andrew F. Read, Sylvain Gandon

Abstract: With the unprecedented global vaccination campaign against SARS-CoV-2 attention has now turned to the potential impact of this large-scale intervention on the evolution of the virus. In this perspective we summarize what is currently known about evolution in the context of vaccination from research on other pathogen species, with an eye towards the future evolution of SARS-CoV-2. With the unprecedented global vaccination campaign against SARS-CoV-2 attention has now turned to the potential impact of this large-scale intervention on the evolution of the virus. In this perspective we summarize what is currently known about evolution in the context of vaccination from research on other pathogen species, with an eye towards the future evolution of SARS-CoV-2. △ Less

Submitted 31 January, 2022; v1 submitted 28 September, 2021; originally announced September 2021.

Comments: 12 pages, 5 figures, 3 boxes, 2 appendices, 1 supplemental figure

Empirical Decision Rules for Improving the Uncertainty Reporting of Small Sample System Usability Scale Scores

Authors: Nicholas Clark, Matthew Dabkowski, Patrick Driscoll, Dereck Kennedy, Ian Kloo, Heidy Shi

Abstract: The System Usability Scale (SUS) is a short, survey-based approach used to determine the usability of a system from an end user perspective once a prototype is available for assessment. Individual scores are gathered using a 10-question survey with the survey results reported in terms of central tendency (sample mean) as an estimate of the system's usability (the SUS study score), and confidence i… ▽ More The System Usability Scale (SUS) is a short, survey-based approach used to determine the usability of a system from an end user perspective once a prototype is available for assessment. Individual scores are gathered using a 10-question survey with the survey results reported in terms of central tendency (sample mean) as an estimate of the system's usability (the SUS study score), and confidence intervals on the sample mean are used to communicate uncertainty levels associated with this point estimate. When the number of individuals surveyed is large, the SUS study scores and accompanying confidence intervals relying upon the central limit theorem for support are appropriate. However, when only a small number of users are surveyed, reliance on the central limit theorem falls short, resulting in confidence intervals that suffer from parameter bound violations and interval widths that confound mappings to adjective and other constructed scales. These shortcomings are especially pronounced when the underlying SUS score data is skewed, as it is in many instances. This paper introduces an empirically-based remedy for such small-sample circumstances, proposing a set of decision rules that leverage either an extended bias-corrected accelerated (BCa) bootstrap confidence interval or an empirical Bayesian credibility interval about the sample mean to restore and bolster subsequent confidence interval accuracy. Data from historical SUS assessments are used to highlight shortfalls in current practices and to demonstrate the improvements these alternate approaches offer while remaining statistically defensible. A freely available, online application is introduced and discussed that automates SUS analysis under these decision rules, thereby assisting usability practitioners in adopting the advocated approaches. △ Less

Submitted 16 January, 2021; v1 submitted 2 January, 2021; originally announced January 2021.

A statistical machine learning approach for benchmarking in the presence of complex contextual factors and peer groups

Authors: Daniel W. Kennedy, Jessica Cameron, Paul P. -Y. Wu, Kerrie Mengersen

Abstract: The ability to compare between individuals or organisations fairly is important for the development of robust and meaningful quantitative benchmarks. To make fair comparisons, contextual factors must be taken into account, and comparisons should only be made between similar organisations such as peer groups. Previous benchmarking methods have used linear regression to adjust for contextual factors… ▽ More The ability to compare between individuals or organisations fairly is important for the development of robust and meaningful quantitative benchmarks. To make fair comparisons, contextual factors must be taken into account, and comparisons should only be made between similar organisations such as peer groups. Previous benchmarking methods have used linear regression to adjust for contextual factors, however linear regression is known to be sub-optimal when nonlinear relationships exist between the comparative measure and covariates. In this paper we propose a random forest model for benchmarking that can adjust for these potential nonlinear relationships, and validate the approach in a case-study of high noise data. We provide new visualisations and numerical summaries of the fitted models and comparative measures to facilitate interpretation by both analysts and non-technical audiences. Comparisons can be made across the cohort or within peer groups, and bootstrapping provides a means of estimating uncertainty in both adjusted measures and rankings. We conclude that random forest models can facilitate fair comparisons between organisations for quantitative measures including in cases on complex contextual factor relationships, and that the models and outputs are readily interpreted by stakeholders. △ Less

Submitted 16 November, 2020; originally announced November 2020.

Comments: 18 pages, 8 figures

Peer groups for organisational learning: clustering with practical constraints

Authors: Daniel William Kennedy, Jessica Cameron, Paul Pao-Yen Wu, Kerrie Mengersen

Abstract: Peer-grouping is used in many sectors for organisational learning, policy implementation, and benchmarking. Clustering provides a statistical, data-driven method for constructing meaningful peer groups, but peer groups must be compatible with business constraints such as size and stability considerations. Additionally, statistical peer groups are constructed from many different variables, and can… ▽ More Peer-grouping is used in many sectors for organisational learning, policy implementation, and benchmarking. Clustering provides a statistical, data-driven method for constructing meaningful peer groups, but peer groups must be compatible with business constraints such as size and stability considerations. Additionally, statistical peer groups are constructed from many different variables, and can be difficult to understand, especially for non-statistical audiences. We developed methodology to apply business constraints to clustering solutions and allow the decision-maker to choose the balance between statistical goodness-of-fit and conformity to business constraints. Several tools were utilised to identify complex distinguishing features in peer groups, and a number of visualisations are developed to explain high-dimensional clusters for non-statistical audiences. In a case study where peer group size was required to be small ($\leq 100$ members), we applied constrained clustering to a noisy high-dimensional data-set over two subsequent years, ensuring that the clusters were sufficiently stable between years. Our approach not only satisfied clustering constraints on the test data, but maintained an almost monotonic negative relationship between goodness-of-fit and stability between subsequent years. We demonstrated in the context of the case study how distinguishing features between clusters can be communicated clearly to different stakeholders with substantial and limited statistical knowledge. △ Less

Submitted 16 November, 2020; originally announced November 2020.

Comments: 22 pages, 4 figures

Deep Generative Modeling in Network Science with Applications to Public Policy Research

Authors: Gavin S. Hartnett, Raffaele Vardavas, Lawrence Baker, Michael Chaykowsky, C. Ben Gibson, Federico Girosi, David P. Kennedy, Osonde A. Osoba

Abstract: Network data is increasingly being used in quantitative, data-driven public policy research. These are typically very rich datasets that contain complex correlations and inter-dependencies. This richness both promises to be quite useful for policy research, while at the same time posing a challenge for the useful extraction of information from these datasets - a challenge which calls for new data… ▽ More Network data is increasingly being used in quantitative, data-driven public policy research. These are typically very rich datasets that contain complex correlations and inter-dependencies. This richness both promises to be quite useful for policy research, while at the same time posing a challenge for the useful extraction of information from these datasets - a challenge which calls for new data analysis methods. In this report, we formulate a research agenda of key methodological problems whose solutions would enable new advances across many areas of policy research. We then review recent advances in applying deep learning to network data, and show how these methods may be used to address many of the methodological problems we identified. We particularly emphasize deep generative methods, which can be used to generate realistic synthetic networks useful for microsimulation and agent-based models capable of informing key public policy questions. We extend these recent advances by developing a new generative framework which applies to large social contact networks commonly used in epidemiological modeling. For context, we also compare and contrast these recent neural network-based approaches with the more traditional Exponential Random Graph Models. Lastly, we discuss some open problems where more progress is needed. △ Less

Submitted 16 October, 2020; v1 submitted 15 October, 2020; originally announced October 2020.

Comments: 77 pages, 38 figures. This article is a RAND Working Report. v2: corrected typo in metadata

Report number: WR-A843-1

Physics-informed Gaussian Process for Online Optimization of Particle Accelerators

Authors: Adi Hanuka, X. Huang, J. Shtalenkova, D. Kennedy, A. Edelen, V. R. Lalchand, D. Ratner, J. Duris

Abstract: High-dimensional optimization is a critical challenge for operating large-scale scientific facilities. We apply a physics-informed Gaussian process (GP) optimizer to tune a complex system by conducting efficient global search. Typical GP models learn from past observations to make predictions, but this reduces their applicability to new systems where archive data is not available. Instead, here we… ▽ More High-dimensional optimization is a critical challenge for operating large-scale scientific facilities. We apply a physics-informed Gaussian process (GP) optimizer to tune a complex system by conducting efficient global search. Typical GP models learn from past observations to make predictions, but this reduces their applicability to new systems where archive data is not available. Instead, here we use a fast approximate model from physics simulations to design the GP model. The GP is then employed to make inferences from sequential online observations in order to optimize the system. Simulation and experimental studies were carried out to demonstrate the method for online control of a storage ring. We show that the physics-informed GP outperforms current routinely used online optimizers in terms of convergence speed, and robustness on this task. The ability to inform the machine-learning model with physics may have wide applications in science. △ Less

Submitted 8 September, 2020; originally announced September 2020.

Journal ref: Phys. Rev. Accel. Beams 24, 072802 (2021)

What the reproductive number R_0 can and cannot tell us about COVID-19 dynamics

Authors: Clara L. Shaw, David A. Kennedy

Abstract: The reproductive number R_0 (and its value after initial disease emergence R) has long been used to predict the likelihood of pathogen invasion, to gauge the potential severity of an epidemic, and to set policy around interventions. However, often ignored complexities have generated confusion around use of the metric. This is particularly apparent with the emergent pandemic virus SARS-CoV-2, the c… ▽ More The reproductive number R_0 (and its value after initial disease emergence R) has long been used to predict the likelihood of pathogen invasion, to gauge the potential severity of an epidemic, and to set policy around interventions. However, often ignored complexities have generated confusion around use of the metric. This is particularly apparent with the emergent pandemic virus SARS-CoV-2, the causative agent of COVID-19. We address some of these misconceptions, namely, how R changes over time, varies over space, and relates to epidemic size by referencing the mathematical definition of R and examples from the current pandemic. We hope that a better appreciation of the uses, nuances, and limitations of R facilitates a better understanding of epidemic spread, epidemic severity, and the effects of interventions in the context of SARS-CoV-2. △ Less

Submitted 25 June, 2020; originally announced June 2020.

Comments: 25 pages, 2 figures

Coulomb collisions in strongly anisotropic plasmas II. Cyclotron cooling in laboratory pair plasmas

Authors: Daniel Kennedy, Per Helander

Abstract: The behaviour of a strongly-magnetized collisional electron-positron plasma which is optically thin to cyclotron radiation is considered, and the distribution functions accessible to it on the various timescales in the system are calculated. Particular attention is paid to the limit in which the collision time exceeds the radiation emission time, making the electron distribution function strongly… ▽ More The behaviour of a strongly-magnetized collisional electron-positron plasma which is optically thin to cyclotron radiation is considered, and the distribution functions accessible to it on the various timescales in the system are calculated. Particular attention is paid to the limit in which the collision time exceeds the radiation emission time, making the electron distribution function strongly anisotropic. Indeed, these are the exact conditions likely to be attained in the first laboratory electron-positron plasma experiments currently being developed, which will typically have very low densities and be confined in very strong magnetic fields. The constraint of strong-magnetization adds an additional complication in that long-range Coulomb collisions, which are usually negligible, must now be considered. A rigorous collision operator for these long-range collisions has never been written down. Nevertheless, we show that the collisional scattering can be accounted for without knowing the explicit form of this collision operator. The rate of radiation emission is calculated and it is found that the loss of energy from the plasma is proportional to the parallel collision frequency multiplied by a factor that only depends logarithmically on plasma parameters. That is, this is a self-accelerating process, meaning that the bulk of the energy will be lost in a few collision times. We show that in a simple case, that of straight field-line geometry, there are no unstable drift waves in such plasmas, despite being far from Maxwellian. △ Less

Submitted 8 September, 2020; v1 submitted 16 June, 2020; originally announced June 2020.

Coulomb collisions in strongly anisotropic plasmas I. Cyclotron cooling in electron-ion plasmas

Authors: Daniel Kennedy, Per Helander

Abstract: The behaviour of a collisional plasma which is optically thin to cyclotron radiation is considered, and the distribution functions accessible to it on the various timescales in the system are calculated. Particular attention is paid to the limit in which the collision time exceeds the radiation emission time, making the electron distribution function strongly anisotropic. Unusually for plasma phys… ▽ More The behaviour of a collisional plasma which is optically thin to cyclotron radiation is considered, and the distribution functions accessible to it on the various timescales in the system are calculated. Particular attention is paid to the limit in which the collision time exceeds the radiation emission time, making the electron distribution function strongly anisotropic. Unusually for plasma physics, the collision operator can nevertheless be calculated analytically although the plasma is far from Maxwellian. The rate of radiation emission is calculated and found to be governed by the collision frequency multiplied by a factor that only depends logarithmically on plasma parameters. △ Less

Submitted 28 October, 2020; v1 submitted 16 June, 2020; originally announced June 2020.

Mathematical model of COVID-19 intervention scenarios for Sao Paulo- Brazil

Authors: Osmar Pinto Neto, Jose Clark Reis, Ana Carolina Brisola Brizzi, Gustavo Jose Zambrano, Joabe Marcos de Souza, Wellington Amorim Pedroso, Rodrigo Cunha de Mello Pedreiro, Bruno de Matos Brizzi, Ellysson Oliveira Abinader, Deanna M. Kennedy, Renato Amaro Zangaro

Abstract: An epidemiological compartmental model was used to simulate social distancing strategies to contain the COVID-19 pandemic and prevent a second wave in Sao Paulo, Brazil. Optimization using genetic algorithm was used to determine the optimal solutions. Our results suggest the best-case strategy for Sao Paulo is to maintain or increase the current magnitude of social distancing for at least 60 more… ▽ More An epidemiological compartmental model was used to simulate social distancing strategies to contain the COVID-19 pandemic and prevent a second wave in Sao Paulo, Brazil. Optimization using genetic algorithm was used to determine the optimal solutions. Our results suggest the best-case strategy for Sao Paulo is to maintain or increase the current magnitude of social distancing for at least 60 more days and increase the current levels of personal protection behaviors by a minimum of 10% (e.g., wearing facemasks, proper hand hygiene and avoid agglomeration). Followed by a long-term oscillatory level of social distancing with a stepping-down approach every 80 days over a period of two years with continued protective behavior. △ Less

Submitted 18 May, 2020; originally announced May 2020.

Comments: 24 pages, 4 figures, 1 table

The turbulent dynamics of Jupiter's and Saturn's weather layers: order out of chaos?

Authors: Peter L Read, Roland M B Young, Daniel Kennedy

Abstract: The weather layers of the gas giant planets, Jupiter and Saturn, comprise the shallow atmospheric layers that are influenced energetically by a combination of incoming solar radiation and localised latent heating of condensates, as well as by upwelling heat from their planetary interiors. They are also the most accessible regions of those planets to direct observations. Recent analyses in Oxford o… ▽ More The weather layers of the gas giant planets, Jupiter and Saturn, comprise the shallow atmospheric layers that are influenced energetically by a combination of incoming solar radiation and localised latent heating of condensates, as well as by upwelling heat from their planetary interiors. They are also the most accessible regions of those planets to direct observations. Recent analyses in Oxford of cloud-tracked winds on Jupiter have demonstrated that kinetic energy is injected into the weather layer at scales comparable to the Rossby radius of deformation and cascades both upscale, mostly into the extra-tropical zonal jets, and downscale to the smallest resolvable scales in Cassini images. The large-scale flow on both Jupiter and Saturn appears to equilibrate towards a state which is close to marginal instability according to Arnol'd's 2nd stability theorem. This scenario is largely reproduced in a hierarchy of numerical models of giant planet weather layers, including relatively realistic models which seek to predict thermal and dynamical structures using a full set of parameterisations of radiative transfer, interior heat sources and even moist convection. Such models include the Jason GCM, developed in Oxford, which also represents the formation of (energetically passive) clouds of NH3, NH4SH and H2O condensates and the transport of condensable tracers. Recent results show some promise in comparison with observations from the Cassini and Juno missions, but some observed features (such as Jupiter's Great Red Spot and other compact ovals) are not yet captured spontaneously by any weather layer model. We review recent work in this vein and discuss a number of open questions for future study. △ Less

Submitted 6 April, 2020; originally announced April 2020.

Comments: 27 pages, including 7 figures. Submitted to Geoscience Letters as invited review from 2019 assembly of the Asia-Oceania Geosciences Society

MSC Class: 85-02

Online tuning and light source control using a physics-informed Gaussian process Adi

Authors: A. Hanuka, J. Duris, J. Shtalenkova, D. Kennedy, A. Edelen, D. Ratner, X. Huang

Abstract: Operating large-scale scientific facilities often requires fast tuning and robust control in a high dimensional space. In this paper we introduce a new physics-informed optimization algorithm based on Gaussian process regression. Our method takes advantage of the existing domain knowledge in the form of realizations of a physics model of the observed system. We have applied a physics-informed Gaus… ▽ More Operating large-scale scientific facilities often requires fast tuning and robust control in a high dimensional space. In this paper we introduce a new physics-informed optimization algorithm based on Gaussian process regression. Our method takes advantage of the existing domain knowledge in the form of realizations of a physics model of the observed system. We have applied a physics-informed Gaussian Process method experimentally at the SPEAR3 storage ring to demonstrate online accelerator optimization. This method outperforms Gaussian Process trained on data as well as the standard approach routinely used for operation, in terms of convergence speed and optimal point. The proposed method could be applicable to automatic tuning and control of other complex systems, without a prerequisite for any observed data. △ Less

Submitted 4 November, 2019; originally announced November 2019.

Journal ref: https://ml4physicalsciences.github.io/2019/files/NeurIPS_ML4PS_2019_85.pdf

Bayesian optimization of a free-electron laser

Authors: Joseph Duris, Dylan Kennedy, Adi Hanuka, Jane Shtalenkova, Auralee Edelen, Adam Egger, Tyler Cope, Daniel Ratner

Abstract: The Linac Coherent Light Source changes configurations multiple times per day, necessitating fast tuning strategies to reduce setup time for successive experiments. To this end, we employ a Bayesian approach to transport optics tuning to optimize groups of quadrupole magnets. We use a Gaussian process to provide a probabilistic model of the machine response with respect to control parameters from… ▽ More The Linac Coherent Light Source changes configurations multiple times per day, necessitating fast tuning strategies to reduce setup time for successive experiments. To this end, we employ a Bayesian approach to transport optics tuning to optimize groups of quadrupole magnets. We use a Gaussian process to provide a probabilistic model of the machine response with respect to control parameters from a modest number of samples. Subsequent samples are selected during optimization using a statistical test combining the model prediction and uncertainty. The model parameters are fit from archived scans, and correlations between devices are added from a simple beam transport model. The result is a sample-efficient optimization routine, which we show significantly outperforms existing optimizers. △ Less

Submitted 12 September, 2019; originally announced September 2019.

Comments: 6 pages, 3 figures

Journal ref: Phys. Rev. Lett. 124, 124801 (2020)

Software Citation Implementation Challenges

Authors: Daniel S. Katz, Daina Bouquin, Neil P. Chue Hong, Jessica Hausman, Catherine Jones, Daniel Chivvis, Tim Clark, Mercè Crosas, Stephan Druskat, Martin Fenner, Tom Gillespie, Alejandra Gonzalez-Beltran, Morane Gruenpeter, Ted Habermann, Robert Haines, Melissa Harrison, Edwin Henneken, Lorraine Hwang, Matthew B. Jones, Alastair A. Kelly, David N. Kennedy, Katrin Leinweber, Fernando Rios, Carly B. Robinson, Ilian Todorov , et al. (2 additional authors not shown)

Abstract: The main output of the FORCE11 Software Citation working group (https://www.force11.org/group/software-citation-working-group) was a paper on software citation principles (https://doi.org/10.7717/peerj-cs.86) published in September 2016. This paper laid out a set of six high-level principles for software citation (importance, credit and attribution, unique identification, persistence, accessibilit… ▽ More The main output of the FORCE11 Software Citation working group (https://www.force11.org/group/software-citation-working-group) was a paper on software citation principles (https://doi.org/10.7717/peerj-cs.86) published in September 2016. This paper laid out a set of six high-level principles for software citation (importance, credit and attribution, unique identification, persistence, accessibility, and specificity) and discussed how they could be used to implement software citation in the scholarly community. In a series of talks and other activities, we have promoted software citation using these increasingly accepted principles. At the time the initial paper was published, we also provided guidance and examples on how to make software citable, though we now realize there are unresolved problems with that guidance. The purpose of this document is to provide an explanation of current issues impacting scholarly attribution of research software, organize updated implementation guidance, and identify where best practices and solutions are still needed. △ Less

Submitted 21 May, 2019; originally announced May 2019.

Hamiltonian Monte Carlo on Symmetric and Homogeneous Spaces via Symplectic Reduction

Authors: Alessandro Barp, Anthony Kennedy, Mark Girolami

Abstract: The Hamiltonian Monte Carlo method generates samples by introducing a mechanical system that explores the target density. For distributions on manifolds it is not always simple to perform the mechanics as a result of the lack of global coordinates, the constraints of the manifold, and the requirement to compute the geodesic flow. In this paper we explain how to construct the Hamiltonian system on… ▽ More The Hamiltonian Monte Carlo method generates samples by introducing a mechanical system that explores the target density. For distributions on manifolds it is not always simple to perform the mechanics as a result of the lack of global coordinates, the constraints of the manifold, and the requirement to compute the geodesic flow. In this paper we explain how to construct the Hamiltonian system on naturally reductive homogeneous spaces using symplectic reduction, which lifts the HMC scheme to a matrix Lie group with global coordinates and constant metric. This provides a general framework that is applicable to many manifolds that arise in applications, such as hyperspheres, hyperbolic spaces, symmetric positive-definite matrices, Grassmannian, and Stiefel manifolds. △ Less

Submitted 18 April, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Evolving embodied intelligence from materials to machines

Authors: David Howard, Agoston E. Eiben, Danielle Frances Kennedy, Jean-Baptiste Mouret, Philip Valencia, Dave Winkler

Abstract: Natural lifeforms specialise to their environmental niches across many levels; from low-level features such as DNA and proteins, through to higher-level artefacts including eyes, limbs, and overarching body plans. We propose Multi-Level Evolution (MLE), a bottom-up automatic process that designs robots across multiple levels and niches them to tasks and environmental conditions. MLE concurrently e… ▽ More Natural lifeforms specialise to their environmental niches across many levels; from low-level features such as DNA and proteins, through to higher-level artefacts including eyes, limbs, and overarching body plans. We propose Multi-Level Evolution (MLE), a bottom-up automatic process that designs robots across multiple levels and niches them to tasks and environmental conditions. MLE concurrently explores constituent molecular and material 'building blocks', as well as their possible assemblies into specialised morphological and sensorimotor configurations. MLE provides a route to fully harness a recent explosion in available candidate materials and ongoing advances in rapid manufacturing processes. We outline a feasible MLE architecture that realises this vision, highlight the main roadblocks and how they may be overcome, and show robotic applications to which MLE is particularly suited. By forming a research agenda to stimulate discussion between researchers in related fields, we hope to inspire the pursuit of multi-level robotic design all the way from material to machine. △ Less

Submitted 17 January, 2019; originally announced January 2019.

Journal ref: Nature Machine Intelligence. Vol. 1, Number 1, pages 12--19. 2019

NeuroStorm: Accelerating Brain Science Discovery in the Cloud

Authors: Gregory Kiar, Robert J. Anderson, Alex Baden, Alexandra Badea, Eric W. Bridgeford, Andrew Champion, Vikram Chandrashekhar, Forrest Collman, Brandon Duderstadt, Alan C. Evans, Florian Engert, Benjamin Falk, Tristan Glatard, William R. Gray Roncal, David N. Kennedy, Jeremy Maitin-Shepard, Ryan A. Marren, Onyeka Nnaemeka, Eric Perlman, Sharmishtaas Seshamani, Eric T. Trautman, Daniel J. Tward, Pedro Antonio Valdés-Sosa, Qing Wang, Michael I. Miller , et al. (2 additional authors not shown)

Abstract: Neuroscientists are now able to acquire data at staggering rates across spatiotemporal scales. However, our ability to capitalize on existing datasets, tools, and intellectual capacities is hampered by technical challenges. The key barriers to accelerating scientific discovery correspond to the FAIR data principles: findability, global access to data, software interoperability, and reproducibility… ▽ More Neuroscientists are now able to acquire data at staggering rates across spatiotemporal scales. However, our ability to capitalize on existing datasets, tools, and intellectual capacities is hampered by technical challenges. The key barriers to accelerating scientific discovery correspond to the FAIR data principles: findability, global access to data, software interoperability, and reproducibility/re-usability. We conducted a hackathon dedicated to making strides in those steps. This manuscript is a technical report summarizing these achievements, and we hope serves as an example of the effectiveness of focused, deliberate hackathons towards the advancement of our quickly-evolving field. △ Less

Submitted 20 March, 2018; v1 submitted 8 March, 2018; originally announced March 2018.

Comments: 10 pages, 4 figures, hackathon report

Geometry and Dynamics for Markov Chain Monte Carlo

Authors: Alessandro Barp, Francois-Xavier Briol, Anthony D. Kennedy, Mark Girolami

Abstract: Markov Chain Monte Carlo methods have revolutionised mathematical computation and enabled statistical inference within many previously intractable models. In this context, Hamiltonian dynamics have been proposed as an efficient way of building chains which can explore probability densities efficiently. The method emerges from physics and geometry and these links have been extensively studied by a… ▽ More Markov Chain Monte Carlo methods have revolutionised mathematical computation and enabled statistical inference within many previously intractable models. In this context, Hamiltonian dynamics have been proposed as an efficient way of building chains which can explore probability densities efficiently. The method emerges from physics and geometry and these links have been extensively studied by a series of authors through the last thirty years. However, there is currently a gap between the intuitions and knowledge of users of the methodology and our deep understanding of these theoretical foundations. The aim of this review is to provide a comprehensive introduction to the geometric tools used in Hamiltonian Monte Carlo at a level accessible to statisticians, machine learners and other users of the methodology with only a basic understanding of Monte Carlo methods. This will be complemented with some discussion of the most recent advances in the field which we believe will become increasingly relevant to applied scientists. △ Less

Submitted 8 May, 2017; originally announced May 2017.

Comments: Submitted to "Annual Review of Statistics and Its Applications"

Investigation of New Methods for Numerical Stochastic Perturbation Theory in $\varphi^4$ Theory

Authors: Mattia Dalla Brida, Marco Garofalo, A. D. Kennedy

Abstract: Numerical stochastic perturbation theory is a powerful tool for estimating high-order perturbative expansions in lattice field theory. The standard algorithms based on the Langevin equation, however, suffer from several limitations which in practice restrict the potential of this technique. In this work we investigate some alternative methods which could in principle improve on the standard approa… ▽ More Numerical stochastic perturbation theory is a powerful tool for estimating high-order perturbative expansions in lattice field theory. The standard algorithms based on the Langevin equation, however, suffer from several limitations which in practice restrict the potential of this technique. In this work we investigate some alternative methods which could in principle improve on the standard approach. In particular, we present a study of the recently proposed Instantaneous Stochastic Perturbation Theory, as well as a formulation of numerical stochastic perturbation theory based on Generalized Hybrid Molecular Dynamics algorithms. The viability of these methods is investigated in $\varphi^4$ theory. △ Less

Submitted 8 September, 2017; v1 submitted 13 March, 2017; originally announced March 2017.

Comments: 45 pages, 12 figures. Added new section on cost comparison with Langevin NSPT. Matches published version

Journal ref: Phys. Rev. D 96, 054502 (2017)

Numerical Stochastic Perturbation Theory and Gradient Flow in φ^4 Theory

Authors: Mattia Dalla Brida, Marco Garofalo, Anthony D. Kennedy

Abstract: In this contribution we present an exploratory study of several novel methods for numerical stochastic perturbation theory. For the investigation we consider observables defined through the gradient flow in the simple φ^4 theory. In this contribution we present an exploratory study of several novel methods for numerical stochastic perturbation theory. For the investigation we consider observables defined through the gradient flow in the simple φ^4 theory. △ Less

Submitted 27 December, 2015; originally announced December 2015.

Comments: Talk given at the 33rd International Symposium on Lattice Field Theory, 14 -18 July, 2015, Kobe, Japan; LaTeX source, 7 pages, 7 figures

Renormalizability of the Schrödinger Functional

Authors: A. D. Kennedy, Stefan Sint

Abstract: Symanzik showed that quantum field theory can be formulated on a space with boundaries by including suitable surface interactions in the action to implement boundary conditions. We show that to all orders in perturbation theory all the divergences induced by these surface interactions can be absorbed by a renormalization of their coefficients. Symanzik showed that quantum field theory can be formulated on a space with boundaries by including suitable surface interactions in the action to implement boundary conditions. We show that to all orders in perturbation theory all the divergences induced by these surface interactions can be absorbed by a renormalization of their coefficients. △ Less

Submitted 26 November, 2015; originally announced November 2015.

Comments: Parallel talk presented at Lattice 2015, Kobe, Japan

Performance of two Askaryan Radio Array stations and first results in the search for ultra-high energy neutrinos

Authors: ARA Collaboration, P. Allison, R. Bard, J. J. Beatty, D. Z. Besson, C. Bora, C. -C. Chen, C. -H. Chen, P. Chen, A. Christenson, A. Connolly, J. Davies, M. Duvernois, B. Fox, R. Gaior, P. W. Gorham, K. Hanson, J. Haugen, B. Hill, K. D. Hoffman, E. Hong, S. -Y. Hsu, L. Hu, J. -J. Huang, M. -H. A. Huang , et al. (42 additional authors not shown)

Abstract: Ultra-high energy neutrinos are interesting messenger particles since, if detected, they can transmit exclusive information about ultra-high energy processes in the Universe. These particles, with energies above $10^{16}\mathrm{eV}$, interact very rarely. Therefore, detectors that instrument several gigatons of matter are needed to discover them. The ARA detector is currently being constructed at… ▽ More Ultra-high energy neutrinos are interesting messenger particles since, if detected, they can transmit exclusive information about ultra-high energy processes in the Universe. These particles, with energies above $10^{16}\mathrm{eV}$, interact very rarely. Therefore, detectors that instrument several gigatons of matter are needed to discover them. The ARA detector is currently being constructed at South Pole. It is designed to use the Askaryan effect, the emission of radio waves from neutrino-induced cascades in the South Pole ice, to detect neutrino interactions at very high energies. With antennas distributed among 37 widely-separated stations in the ice, such interactions can be observed in a volume of several hundred cubic kilometers. Currently 3 deep ARA stations are deployed in the ice of which two have been taking data since the beginning of the year 2013. In this publication, the ARA detector "as-built" and calibrations are described. Furthermore, the data reduction methods used to distinguish the rare radio signals from overwhelming backgrounds of thermal and anthropogenic origin are presented. Using data from only two stations over a short exposure time of 10 months, a neutrino flux limit of $3 \cdot 10^{-6} \mathrm{GeV} / (\mathrm{cm^2 \ s \ sr})$ is calculated for a particle energy of 10^{18}eV, which offers promise for the full ARA detector. △ Less

Submitted 16 May, 2016; v1 submitted 31 July, 2015; originally announced July 2015.

Comments: 21 pages, 34 figures, 1 table, includes supplementary material

Journal ref: Phys. Rev. D 93, 082003 (2016)

Design, Modeling and Testing of the Askaryan Radio Array South Pole Autonomous Renewable Power Stations

Authors: D. Z. Besson, D. M. Kennedy, K. Ratzlaff, R. Young

Abstract: We describe the design, construction and operation of the Askaryan Radio Array (ARA) Autonomous Renewable Power Stations, initially installed at the South Pole in December, 2010 with the goal of providing an independently operating 100 W power source capable of year-round operation in extreme environments. In addition to particle astrophysics applications at the South Pole, such a station can easi… ▽ More We describe the design, construction and operation of the Askaryan Radio Array (ARA) Autonomous Renewable Power Stations, initially installed at the South Pole in December, 2010 with the goal of providing an independently operating 100 W power source capable of year-round operation in extreme environments. In addition to particle astrophysics applications at the South Pole, such a station can easily be, and has since been, extended to operation elsewhere, as described herein. △ Less

Submitted 6 March, 2014; originally announced March 2014.

Comments: 33 pages, 28 figures

Optical hyperpolarization and NMR detection of $^{129}$Xe on a microfluidic chip

Authors: Ricardo Jimenez-Martinez, Daniel J. Kennedy, Michael Rosenbluh, Elizabeth A. Donley, Svenja Knappe, Scott J. Seltzer, Hattie L. Ring, Vikram S. Bajaj, John Kitching

Abstract: Optically hyperpolarized $^{129}$Xe gas has become a powerful contrast agent in nuclear magnetic resonance (NMR) spectroscopy and imaging, with applications ranging from studies of the human lung to the targeted detection of biomolecules. Equally attractive is its potential use to enhance the sensitivity of microfluidic NMR experiments, in which small sample volumes yield poor sensitivity. Unfortu… ▽ More Optically hyperpolarized $^{129}$Xe gas has become a powerful contrast agent in nuclear magnetic resonance (NMR) spectroscopy and imaging, with applications ranging from studies of the human lung to the targeted detection of biomolecules. Equally attractive is its potential use to enhance the sensitivity of microfluidic NMR experiments, in which small sample volumes yield poor sensitivity. Unfortunately, most $^{129}$Xe polarization systems are large and non-portable. Here we present a microfabricated chip that optically polarizes $^{129}$Xe gas. We have achieved $^{129}$Xe polarizations greater than 0.5$\%$ at flow rates of several microliters per second, compatible with typical microfluidic applications. We employ in situ optical magnetometry to sensitively detect and characterize the $^{129}$Xe polarization at magnetic fields of 1 $μ$T. We construct the device using standard microfabrication techniques, which will facilitate its integration with existing microfluidic platforms. This device may enable the implementation of highly sensitive $^{129}$Xe NMR in compact, low-cost, portable devices. △ Less

Submitted 30 December, 2013; originally announced January 2014.

How Scaling Symmetry Solves a Second-Order Differential Equation

Authors: Sidney Bludman, Andres Guzman, Dallas C. Kennedy

Abstract: While not generally a conservation law, any symmetry of the equations of motion implies a useful reduction of any second-order equationto a first-order equation between invariants, whose solutions (first integrals) can then be integrated by quadrature (Lie's Theorem on the solvability of differential equations). We illustrate this theorem by applying scale invariance to the equations for the hydro… ▽ More While not generally a conservation law, any symmetry of the equations of motion implies a useful reduction of any second-order equationto a first-order equation between invariants, whose solutions (first integrals) can then be integrated by quadrature (Lie's Theorem on the solvability of differential equations). We illustrate this theorem by applying scale invariance to the equations for the hydrostatic equilibrium of stars in local thermodynamic equilibrium: Scaling symmetry reduces the Lane-Emden equation to a first-order equation between scale invariants un; vn, whose phase diagram encapsulates all the properties of index-n polytropes. From this reduced equation, we obtain the regular (Emden) solutions and demonstrate graphically how they transform under scale transformations. △ Less

Submitted 1 December, 2012; originally announced December 2012.

Comments: 7 pages, 6 figures

Numerical determination of partial spectrum of Hermitian matrices using a Lanczos method with selective reorthogonalization

Authors: Chris Johnson, A. D. Kennedy

Abstract: We introduce a new algorithm for finding the eigenvalues and eigenvectors of Hermitian matrices within a specified region, based upon the LANSO algorithm of Parlett and Scott. It uses selective reorthogonalization to avoid the duplication of eigenpairs in finite-precision arithmetic, but uses a new bound to decide when such reorthogonalization is required, and only reorthogonalizes with respect to… ▽ More We introduce a new algorithm for finding the eigenvalues and eigenvectors of Hermitian matrices within a specified region, based upon the LANSO algorithm of Parlett and Scott. It uses selective reorthogonalization to avoid the duplication of eigenpairs in finite-precision arithmetic, but uses a new bound to decide when such reorthogonalization is required, and only reorthogonalizes with respect to eigenpairs within the region of interest. We investigate its performance for the Hermitian Wilson--Dirac operator (γ_5D) in lattice quantum chromodynamics, and compare it with previous methods. △ Less

Submitted 6 November, 2012; originally announced November 2012.

Shadow Hamiltonians, Poisson Brackets, and Gauge Theories

Authors: A. D. Kennedy, P. J. Silva, M. A. Clark

Abstract: Numerical lattice gauge theory computations to generate gauge field configurations including the effects of dynamical fermions are usually carried out using algorithms that require the molecular dynamics evolution of gauge fields using symplectic integrators. Sophisticated integrators are in common use but are hard to optimise, and force-gradient integrators show promise especially for large latti… ▽ More Numerical lattice gauge theory computations to generate gauge field configurations including the effects of dynamical fermions are usually carried out using algorithms that require the molecular dynamics evolution of gauge fields using symplectic integrators. Sophisticated integrators are in common use but are hard to optimise, and force-gradient integrators show promise especially for large lattice volumes. We explain why symplectic integrators lead to very efficient Monte Carlo algorithms because they exactly conserve a shadow Hamiltonian. The shadow Hamiltonian may be expanded in terms of Poisson brackets, and can be used to optimize the integrators. We show how this may be done for gauge theories by extending the formulation of Hamiltonian mechanics on Lie groups to include Poisson brackets and shadows, and by giving a general method for the practical computation of forces, force-gradients, and Poisson brackets for gauge theories. △ Less

Submitted 24 October, 2012; originally announced October 2012.

The Electroweak Fit of the Standard Model after the Discovery of a New Boson at the LHC

Authors: M. Baak, M. Goebel, J. Haller, A. Hoecker, D. Kennedy, R. Kogler, K. Moenig, M. Schott, J. Stelzer

Abstract: In view of the discovery of a new boson by the ATLAS and CMS Collaborations at the LHC, we present an update of the global Standard Model (SM) fit to electroweak precision data. Assuming the new particle to be the SM Higgs boson, all fundamental parameters of the SM are known allowing, for the first time, to overconstrain the SM at the electroweak scale and assert its validity. Including the effec… ▽ More In view of the discovery of a new boson by the ATLAS and CMS Collaborations at the LHC, we present an update of the global Standard Model (SM) fit to electroweak precision data. Assuming the new particle to be the SM Higgs boson, all fundamental parameters of the SM are known allowing, for the first time, to overconstrain the SM at the electroweak scale and assert its validity. Including the effects of radiative corrections and the experimental and theoretical uncertainties, the global fit exhibits a p-value of 0.07. The mass measurements by ATLAS and CMS agree within 1.3sigma with the indirect determination M_H=(94 +25 -22) GeV. Within the SM the W boson mass and the effective weak mixing angle can be accurately predicted to be M_W=(80.359 +- 0.011) GeV and sin^2(theta_eff^ell)=(0.23150 +- 0.00010) from the global fit. These results are compatible with, and exceed in precision, the direct measurements. For the indirect determination of the top quark mass we find m_t=(175.8 +2.7 -2.4) GeV, in agreement with the kinematic and cross-section based measurements. △ Less

Submitted 25 September, 2012; v1 submitted 12 September, 2012; originally announced September 2012.

Comments: 11 pages, 5 figures, version submitted to EPJ-C

Report number: DESY-12-154

Journal ref: Eur. Phys. J. C72 (2012) 2205

Monte Carlo Integration with Subtraction

Authors: Rudy Arthur, A. D. Kennedy

Abstract: This paper investigates a class of algorithms for numerical integration of a function in d dimensions over a compact domain by Monte Carlo methods. We construct a histogram approximation to the function using a partition of the integration domain into a set of bins specified by some parameters. We then consider two adaptations; the first is to subtract the histogram approximation, whose integral w… ▽ More This paper investigates a class of algorithms for numerical integration of a function in d dimensions over a compact domain by Monte Carlo methods. We construct a histogram approximation to the function using a partition of the integration domain into a set of bins specified by some parameters. We then consider two adaptations; the first is to subtract the histogram approximation, whose integral we may easily evaluate explicitly, from the function and integrate the difference using Monte Carlo; the second is to modify the bin parameters in order to make the variance of the Monte Carlo estimate of the integral the same for all bins. This allows us to use Student's t-test as a trigger for rebinning, which we claim is more stable than the χ-squared test that is commonly used for this purpose. We provide a program that we have used to study the algorithm for the case where the histogram is represented as a product of one-dimensional histograms. We discuss the assumptions and approximations made, as well as giving a pedagogical discussion of the myriad ways in which the results of any such Monte Carlo integration program can be misleading. △ Less

Submitted 4 September, 2012; originally announced September 2012.

Comments: Code PANIC included as a set of ancillary files

Room Temperature Operation of a Radiofrequency Diamond Magnetometer near the Shot Noise Limit

Authors: Chang S. Shin, Claudia E. Avalos, Mark C. Butler, David R. Trease, Scott J. Seltzer, J. Peter Mustonen, Daniel J. Kennedy, Victor M. Acosta, Dmitry Budker, Alexander Pines, Vikram S. Bajaj

Abstract: We operate a nitrogen vacancy (NV-) diamond magnetometer at ambient temperatures and study the dependence of its bandwidth on experimental parameters including optical and microwave excitation powers. We introduce an analytical theory that yields an explicit formula for the response of an ensemble of NV- spins to an oscillating magnetic field, such as in NMR applications. We measure a detection ba… ▽ More We operate a nitrogen vacancy (NV-) diamond magnetometer at ambient temperatures and study the dependence of its bandwidth on experimental parameters including optical and microwave excitation powers. We introduce an analytical theory that yields an explicit formula for the response of an ensemble of NV- spins to an oscillating magnetic field, such as in NMR applications. We measure a detection bandwidth of 1.6 MHz and a sensitivity of 4.6 nT/Hz^(1/2), unprecedented in a detector with this active volume and close to the photon shot noise limit of our experiment. △ Less

Submitted 15 January, 2012; originally announced January 2012.

Comments: 4 pages; 4 figures; supporting information

Scaling Symmetry and Integrable Spherical Hydrostatics

Authors: Sidney Bludman, Dallas C. Kennedy

Abstract: Any symmetry reduces a second-order differential equation to a first integral: variational symmetries of the action (exemplified by central field dynamics) lead to conservation laws, but symmetries of only the equations of motion (exemplified by scale-invariant hydrostatics) yield first-order {\em non-conservation laws} between invariants. We obtain these non-conservation laws by extending Noether… ▽ More Any symmetry reduces a second-order differential equation to a first integral: variational symmetries of the action (exemplified by central field dynamics) lead to conservation laws, but symmetries of only the equations of motion (exemplified by scale-invariant hydrostatics) yield first-order {\em non-conservation laws} between invariants. We obtain these non-conservation laws by extending Noether's Theorem to non-variational symmetries and present an innovative variational formulation of spherical adiabatic hydrostatics. For the scale-invariant case, this novel synthesis of group theory, hydrostatics, and astrophysics allows us to recover all the known properties of polytropes and define a {\em core radius}, inside which polytropes of index $n$ share a common core mass density structure, and outside of which their envelopes differ. The Emden solutions (regular solutions of the Lane-Emden equation) are obtained, along with useful approximations. An appendix discusses the $n=3$ polytrope in order to emphasize how the same mechanical structure allows different thermal structures in relativistic degenerate white dwarfs and zero age main sequence stars. △ Less

Submitted 6 November, 2012; v1 submitted 18 December, 2011; originally announced December 2011.

Comments: 10 pages, 4 figures, 2 tables. arXiv admin note: substantial text overlap with arXiv:1106.1222

Journal ref: Journal of Modern Physics, Vol. 4 No. 4 (April 2013) 486-494

Bounds on the convergence of Ritz values from Krylov subspaces to interior eigenvalues of Hermitean matrices

Authors: Chris Johnson, A. D. Kennedy

Abstract: We consider bounds on the convergence of Ritz values from a sequence of Krylov subspaces to interior eigenvalues of Hermitean matrices. These bounds are useful in regions of low spectral density, for example near voids in the spectrum, as is required in many applications. Our bounds are obtained by considering the usual Kaniel-Paige-Saad formalism applied to the shifted and squared matrix. We consider bounds on the convergence of Ritz values from a sequence of Krylov subspaces to interior eigenvalues of Hermitean matrices. These bounds are useful in regions of low spectral density, for example near voids in the spectrum, as is required in many applications. Our bounds are obtained by considering the usual Kaniel-Paige-Saad formalism applied to the shifted and squared matrix. △ Less

Submitted 13 October, 2011; originally announced October 2011.

Comments: 1 figure

MSC Class: 65F15 (Primary) 15A18; 15B57; 65Z05; 65B99 (Secondary)

Improving dynamical lattice QCD simulations through integrator tuning using Poisson brackets and a force-gradient integrator

Authors: M. A. Clark, Bálint Joó, A. D. Kennedy, P. J. Silva

Abstract: We show how the integrators used for the molecular dynamics step of the Hybrid Monte Carlo algorithm can be further improved. These integrators not only approximately conserve some Hamiltonian $H$ but conserve exactly a nearby shadow Hamiltonian $\tilde{H}$. This property allows for a new tuning method of the molecular dynamics integrator and also allows for a new class of integrators (force-gradi… ▽ More We show how the integrators used for the molecular dynamics step of the Hybrid Monte Carlo algorithm can be further improved. These integrators not only approximately conserve some Hamiltonian $H$ but conserve exactly a nearby shadow Hamiltonian $\tilde{H}$. This property allows for a new tuning method of the molecular dynamics integrator and also allows for a new class of integrators (force-gradient integrators) which is expected to reduce significantly the computational cost of future large-scale gauge field ensemble generation. △ Less

Submitted 20 September, 2011; v1 submitted 8 August, 2011; originally announced August 2011.

Comments: 5 pages, 1 figure; minor changes

Updated Status of the Global Electroweak Fit and Constraints on New Physics

Authors: M. Baak, M. Goebel, J. Haller, A. Hoecker, D. Kennedy, K. Moenig, M. Schott, J. Stelzer

Abstract: We present an update of the Standard Model fit to electroweak precision data. We include newest experimental results on the top quark mass, the W mass and width, and the Higgs boson mass bounds from LEP, Tevatron and the LHC. We also include a new determination of the electromagnetic coupling strength at the Z pole. We find for the Higgs boson mass (91 +30 -23) GeV and (120 +12 -5) GeV when not in… ▽ More We present an update of the Standard Model fit to electroweak precision data. We include newest experimental results on the top quark mass, the W mass and width, and the Higgs boson mass bounds from LEP, Tevatron and the LHC. We also include a new determination of the electromagnetic coupling strength at the Z pole. We find for the Higgs boson mass (91 +30 -23) GeV and (120 +12 -5) GeV when not including and including the direct Higgs searches, respectively. From the latter fit we indirectly determine the W mass to be (80.360 +0.014 -0.013) GeV. We exploit the data to determine experimental constraints on the oblique vacuum polarisation parameters, and confront these with predictions from the Standard Model (SM) and selected SM extensions. By fitting the oblique parameters to the electroweak data we derive allowed regions in the BSM parameter spaces. We revisit and consistently update these constraints for a fourth fourth fermion generation, two Higgs doublet, inert Higgs and littlest Higgs models, models with large, universal or warped extra dimensions and technicolour. In most of the models studied a heavy Higgs boson can be made compatible with the electroweak precision data. △ Less

Submitted 11 September, 2012; v1 submitted 5 July, 2011; originally announced July 2011.

Comments: 58 pages, 27 figures, published version

Report number: CERN-OPEN-2011-033, DESY-11-107

Journal ref: Eur.Phys.J. C72 (2012) 2003

Design and Initial Performance of the Askaryan Radio Array Prototype EeV Neutrino Detector at the South Pole

Authors: P. Allison, J. Auffenberg, R. Bard, J. J. Beatty, D. Z. Besson, S. Boeser, C. Chen, P. Chen, A. Connolly, J. Davies, M. DuVernois, B. Fox, P. W. Gorham, E. W. Grashorn, K. Hanson, J. Haugen, K. Helbing, B. Hill, K. D. Hoffman, M. Huang, M. H. A. Huang, A. Ishihara, A. Karle, D. Kennedy, H. Landsman , et al. (23 additional authors not shown)

Abstract: We report on studies of the viability and sensitivity of the Askaryan Radio Array (ARA), a new initiative to develop a Teraton-scale ultra-high energy neutrino detector in deep, radio-transparent ice near Amundsen-Scott station at the South Pole. An initial prototype ARA detector system was installed in January 2011, and has been operating continuously since then. We report on studies of the backg… ▽ More We report on studies of the viability and sensitivity of the Askaryan Radio Array (ARA), a new initiative to develop a Teraton-scale ultra-high energy neutrino detector in deep, radio-transparent ice near Amundsen-Scott station at the South Pole. An initial prototype ARA detector system was installed in January 2011, and has been operating continuously since then. We report on studies of the background radio noise levels, the radio clarity of the ice, and the estimated sensitivity of the planned ARA array given these results, based on the first five months of operation. Anthropogenic radio interference in the vicinity of the South Pole currently leads to a few-percent loss of data, but no overall effect on the background noise levels, which are dominated by the thermal noise floor of the cold polar ice, and galactic noise at lower frequencies. We have also successfully detected signals originating from a 2.5 km deep impulse generator at a distance of over 3 km from our prototype detector, confirming prior estimates of kilometer-scale attenuation lengths for cold polar ice. These are also the first such measurements for propagation over such large slant distances in ice. Based on these data, ARA-37, the 200 km^2 array now under construction, will achieve the highest sensitivity of any planned or existing neutrino detector in the 10^{16}-10^{19} eV energy range. △ Less

Submitted 9 August, 2011; v1 submitted 13 May, 2011; originally announced May 2011.

Comments: 25 pages, 37 figures, this version with improved ice attenuation length analysis; for submission to Astroparticle Physics