Showing 1–50 of 211 results for author: Harris, R

Far from Perfect: Quantum Error Correction with (Hyperinvariant) Evenbly Codes

Authors: Matthew Steinberg, Junyu Fan, Robert J. Harris, David Elkouss, Sebastian Feld, Alexander Jahn

Abstract: We introduce a new class of qubit codes that we call Evenbly codes, building on a previous proposal of hyperinvariant tensor networks. Its tensor network description consists of local, non-perfect tensors describing CSS codes interspersed with Hadamard gates, placed on a hyperbolic $\{p,q\}$ geometry with even $q\geq 4$, yielding an infinitely large class of subsystem codes. We construct an exampl… ▽ More We introduce a new class of qubit codes that we call Evenbly codes, building on a previous proposal of hyperinvariant tensor networks. Its tensor network description consists of local, non-perfect tensors describing CSS codes interspersed with Hadamard gates, placed on a hyperbolic $\{p,q\}$ geometry with even $q\geq 4$, yielding an infinitely large class of subsystem codes. We construct an example for a $\{5,4\}$ manifold and describe strategies of logical gauge fixing that lead to different rates $k/n$ and distances $d$, which we calculate analytically, finding distances which range from $d=2$ to $d \sim n^{2/3}$ in the ungauged case. Investigating threshold performance under erasure, depolarizing, and pure Pauli noise channels, we find that the code exhibits a depolarizing noise threshold of about $19.1\%$ in the code-capacity model and $50\%$ for pure Pauli and erasure channels under suitable gauges. We also test a constant-rate version with $k/n = 0.125$, finding excellent error resilience (about $40\%$) under the erasure channel. Recovery rates for these and other settings are studied both under an optimal decoder as well as a more efficient but non-optimal greedy decoder. We also consider generalizations beyond the CSS tensor construction, compute error rates and thresholds for other hyperbolic geometries, and discuss the relationship to holographic bulk/boundary dualities. Our work indicates that Evenbly codes may show promise for practical quantum computing applications. △ Less

Submitted 16 July, 2024; originally announced July 2024.

Comments: 22 pages, 13 figures

The Blue Multi Unit Spectroscopic Explorer (BlueMUSE) on the VLT: science drivers and overview of instrument design

Authors: Johan Richard, Rémi Giroud, Florence Laurent, Davor Krajnović, Alexandre Jeanneau, Roland Bacon, Manuel Abreu, Angela Adamo, Ricardo Araujo, Nicolas Bouché, Jarle Brinchmann, Zhemin Cai, Norberto Castro, Ariadna Calcines, Diane Chapuis, Adélaïde Claeyssens, Luca Cortese, Emanuele Daddi, Christopher Davison, Michael Goodwin, Robert Harris, Matthew Hayes, Mathilde Jauzac, Andreas Kelz, Jean-Paul Kneib , et al. (24 additional authors not shown)

Abstract: BlueMUSE is a blue-optimised, medium spectral resolution, panoramic integral field spectrograph under development for the Very Large Telescope (VLT). With an optimised transmission down to 350 nm, spectral resolution of R$\sim$3500 on average across the wavelength range, and a large FoV (1 arcmin$^2$), BlueMUSE will open up a new range of galactic and extragalactic science cases facilitated by its… ▽ More BlueMUSE is a blue-optimised, medium spectral resolution, panoramic integral field spectrograph under development for the Very Large Telescope (VLT). With an optimised transmission down to 350 nm, spectral resolution of R$\sim$3500 on average across the wavelength range, and a large FoV (1 arcmin$^2$), BlueMUSE will open up a new range of galactic and extragalactic science cases facilitated by its specific capabilities. The BlueMUSE consortium includes 9 institutes located in 7 countries and is led by the Centre de Recherche Astrophysique de Lyon (CRAL). The BlueMUSE project development is currently in Phase A, with an expected first light at the VLT in 2031. We introduce here the Top Level Requirements (TLRs) derived from the main science cases, and then present an overview of the BlueMUSE system and its subsystems fulfilling these TLRs. We specifically emphasize the tradeoffs that are made and the key distinctions compared to the MUSE instrument, upon which the system architecture is built. △ Less

Submitted 19 June, 2024; originally announced June 2024.

Comments: 20 pages, 10 figures, proceedings of the SPIE astronomical telescopes and instrumentation conference, Yokohama, 16-21 June

Exotic definite four-manifolds with non-cyclic fundamental group

Authors: Robert Harris, Patrick Naylor, B. Doug Park

Abstract: We construct infinitely many pairwise non-diffeomorphic smooth structures on a definite $4$-manifold with non-cyclic fundamental group $\mathbb{Z}/2\times \mathbb{Z}/2$. We construct infinitely many pairwise non-diffeomorphic smooth structures on a definite $4$-manifold with non-cyclic fundamental group $\mathbb{Z}/2\times \mathbb{Z}/2$. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 11 pages, comments welcome!

MSC Class: 57R55; 57K41; 57M10; 14E20

Methodological Reconstruction of Historical Landslide Tsunamis Using Bayesian Inference

Authors: Raelynn Wonnacott, Dallin Stewart, Jared P Whitehead, Ronald A Harris

Abstract: Indonesia is one of the world's most densely populated regions and lies among the epicenters of Earth's greatest natural hazards. Effectively reducing the disaster potential of these hazards through resource allocation and preparedness first requires an analysis of the risk factors of the region. Since destructive tsunamis present one of the most eminent dangers to coastal communities, understandi… ▽ More Indonesia is one of the world's most densely populated regions and lies among the epicenters of Earth's greatest natural hazards. Effectively reducing the disaster potential of these hazards through resource allocation and preparedness first requires an analysis of the risk factors of the region. Since destructive tsunamis present one of the most eminent dangers to coastal communities, understanding their sources and geological history is necessary to determine the potential future risk. Inspired by results from Cummins et al. 2020, and previous efforts that identified source parameters for earthquake-generated tsunamis, we consider landslide-generated tsunamis. This is done by constructing a probability distribution of potential landslide sources based on anecdotal observations of the 1852 Banda Sea tsunami, using Bayesian inference and scientific computing. After collecting over 100,000 samples (simulating 100,000 landslide induced tsunamis), we conclude that a landslide event provides a reasonable match to the tsunami reported in the anecdotal accounts. However, the most viable landslides may push the boundaries of geological plausibility. Future work creating a joint landslide-earthquake model may compensate for the weaknesses associated with an individual landslide or earthquake source event. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Computational supremacy in quantum simulation

Authors: Andrew D. King, Alberto Nocera, Marek M. Rams, Jacek Dziarmaga, Roeland Wiersema, William Bernoudy, Jack Raymond, Nitin Kaushal, Niclas Heinsdorf, Richard Harris, Kelly Boothby, Fabio Altomare, Andrew J. Berkley, Martin Boschnak, Kevin Chern, Holly Christiani, Samantha Cibere, Jake Connor, Martin H. Dehn, Rahul Deshpande, Sara Ejtemaee, Pau Farré, Kelsey Hamer, Emile Hoskinson, Shuiyuan Huang , et al. (37 additional authors not shown)

Abstract: Quantum computers hold the promise of solving certain problems that lie beyond the reach of conventional computers. Establishing this capability, especially for impactful and meaningful problems, remains a central challenge. One such problem is the simulation of nonequilibrium dynamics of a magnetic spin system quenched through a quantum phase transition. State-of-the-art classical simulations dem… ▽ More Quantum computers hold the promise of solving certain problems that lie beyond the reach of conventional computers. Establishing this capability, especially for impactful and meaningful problems, remains a central challenge. One such problem is the simulation of nonequilibrium dynamics of a magnetic spin system quenched through a quantum phase transition. State-of-the-art classical simulations demand resources that grow exponentially with system size. Here we show that superconducting quantum annealing processors can rapidly generate samples in close agreement with solutions of the Schrödinger equation. We demonstrate area-law scaling of entanglement in the model quench in two-, three- and infinite-dimensional spin glasses, supporting the observed stretched-exponential scaling of effort for classical approaches. We assess approximate methods based on tensor networks and neural networks and conclude that no known approach can achieve the same accuracy as the quantum annealer within a reasonable timeframe. Thus quantum annealers can answer questions of practical importance that classical computers cannot. △ Less

Submitted 1 March, 2024; originally announced March 2024.

A back-to-back diode model applied to MoS2 van der Waals Schottky diodes

Authors: Jeffrey A. Cloninger, Raine Harris, Kristine L. Haley, Randy M. Sterbentz, Takashi Taniguchi, Kenji Watanabe, Joshua O. Island

Abstract: The use of metal van der Waals contacts and the implicit reduction in Fermi-level pinning in contacted semiconductors has led to remarkable device optimizations. For example, using graphene as an electrical contact allows for tunable Schottky barriers in transistors and barristors. In this study, we present a double Schottky barrier model and apply it to barrier tunable all van der Waals transisto… ▽ More The use of metal van der Waals contacts and the implicit reduction in Fermi-level pinning in contacted semiconductors has led to remarkable device optimizations. For example, using graphene as an electrical contact allows for tunable Schottky barriers in transistors and barristors. In this study, we present a double Schottky barrier model and apply it to barrier tunable all van der Waals transistors. In a molybdenum disulfide (MoS$_2$) transistor with graphene and few-layer graphene contacts, we find that the model can be applied to extract Schottky barrier heights that agree with the Schottky-Mott rule from simple two-terminal current-voltage measurements at room temperature. Furthermore, we show tunability of the Schottky barrier \textit{in-situ} using a regional contact gate. Our results show that a basic back-to-back diode model, applied to two terminal measurements, can capture the diode properties of all-van-der-Waals transistors relatively well. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: 6 pages, 3 figures

Optimally Managing the Impacts of Convergence Tolerance for Distributed Optimal Power Flow

Authors: Rachel Harris, Mohannad Alkhraijah, Daniel K. Molzahn

Abstract: The future power grid may rely on distributed optimization to determine the set-points for huge numbers of distributed energy resources. There has been significant work on applying distributed algorithms to optimal power flow (OPF) problems, which require separate computing agents to agree on shared boundary variable values. Looser tolerances for the mismatches in these shared variables generally… ▽ More The future power grid may rely on distributed optimization to determine the set-points for huge numbers of distributed energy resources. There has been significant work on applying distributed algorithms to optimal power flow (OPF) problems, which require separate computing agents to agree on shared boundary variable values. Looser tolerances for the mismatches in these shared variables generally yield faster convergence at the expense of exacerbating constraint violations, but there is little quantitative understanding of how the convergence tolerance affects solution quality. To address this gap, we first quantify how convergence tolerance impacts constraint violations when the distributed OPF generator dispatch is applied to the power system. Using insights from this analysis, we then develop a bound tightening algorithm which guarantees that operating points from distributed OPF algorithms will not result in violations despite the possibility of shared variable mismatches within the convergence tolerance. We also explore how bounding the cumulative shared variable mismatches can prevent unnecessary conservativeness in the bound tightening. The proposed approach enables control of the trade-off between computational speed, which improves as the convergence tolerance increases, and distributed OPF solution cost, which increases with convergence tolerance due to tightened constraints, while ensuring feasibility. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Quantum error mitigation in quantum annealing

Authors: Mohammad H. Amin, Andrew D. King, Jack Raymond, Richard Harris, William Bernoudy, Andrew J. Berkley, Kelly Boothby, Anatoly Smirnov, Fabio Altomare, Michael Babcock, Catia Baron, Jake Connor, Martin Dehn, Colin Enderud, Emile Hoskinson, Shuiyuan Huang, Mark W. Johnson, Eric Ladizinsky, Trevor Lanting, Allison J. R. MacDonald, Gaelen Marsden, Reza Molavi, Travis Oh, Gabriel Poulin-Lamarre, Hugh Ramp , et al. (10 additional authors not shown)

Abstract: Quantum Error Mitigation (QEM) presents a promising near-term approach to reduce error when estimating expectation values in quantum computing. Here, we introduce QEM techniques tailored for quantum annealing, using Zero-Noise Extrapolation (ZNE). We implement ZNE through zero-temperature extrapolation as well as energy-time rescaling. We conduct experimental investigations into the quantum critic… ▽ More Quantum Error Mitigation (QEM) presents a promising near-term approach to reduce error when estimating expectation values in quantum computing. Here, we introduce QEM techniques tailored for quantum annealing, using Zero-Noise Extrapolation (ZNE). We implement ZNE through zero-temperature extrapolation as well as energy-time rescaling. We conduct experimental investigations into the quantum critical dynamics of a transverse-field Ising spin chain, demonstrating the successful mitigation of thermal noise through both of these techniques. Moreover, we show that energy-time rescaling effectively mitigates control errors in the coherent regime where the effect of thermal noise is minimal. Our ZNE results agree with exact calculations of the coherent evolution over a range of annealing times that exceeds the coherent annealing range by almost an order of magnitude. △ Less

Submitted 2 November, 2023; originally announced November 2023.

Comments: 10 pages, 5 figures

2023 Astrophotonics Roadmap: pathways to realizing multi-functional integrated astrophotonic instruments

Authors: Nemanja Jovanovic, Pradip Gatkine, Narsireddy Anugu, Rodrigo Amezcua-Correa, Ritoban Basu Thakur, Charles Beichman, Chad Bender, Jean-Philippe Berger, Azzurra Bigioli, Joss Bland-Hawthorn, Guillaume Bourdarot, Charles M. Bradford, Ronald Broeke, Julia Bryant, Kevin Bundy, Ross Cheriton, Nick Cvetojevic, Momen Diab, Scott A. Diddams, Aline N. Dinkelaker, Jeroen Duis, Stephen Eikenberry, Simon Ellis, Akira Endo, Donald F. Figer , et al. (55 additional authors not shown)

Abstract: Photonics offer numerous functionalities that can be used to realize astrophotonic instruments. The most spectacular example to date is the ESO Gravity instrument at the Very Large Telescope in Chile. Integrated astrophotonic devices stand to offer critical advantages for instrument development, including extreme miniaturization, as well as integration, superior thermal and mechanical stabilizatio… ▽ More Photonics offer numerous functionalities that can be used to realize astrophotonic instruments. The most spectacular example to date is the ESO Gravity instrument at the Very Large Telescope in Chile. Integrated astrophotonic devices stand to offer critical advantages for instrument development, including extreme miniaturization, as well as integration, superior thermal and mechanical stabilization owing to the small footprint, and high replicability offering cost savings. Numerous astrophotonic technologies have been developed to address shortcomings of conventional instruments to date, including for example the development of photonic lanterns, complex aperiodic fiber Bragg gratings, complex beam combiners to enable long baseline interferometry, and laser frequency combs for high precision spectral calibration of spectrometers. Despite these successes, the facility implementation of photonic solutions in astronomical instrumentation is currently limited because of (1) low throughputs from coupling to fibers, coupling fibers to chips, propagation and bend losses, device losses, etc, (2) difficulties with scaling to large channel count devices needed for large bandwidths and high resolutions, and (3) efficient integration of photonics with detectors, to name a few. In this roadmap, we identify 24 areas that need further development. We outline the challenges and advances needed across those areas covering design tools, simulation capabilities, fabrication processes, the need for entirely new components, integration and hybridization and the characterization of devices. To realize these advances the astrophotonics community will have to work cooperatively with industrial partners who have more advanced manufacturing capabilities. With the advances described herein, multi-functional instruments will be realized leading to novel observing capabilities for both ground and space platforms. △ Less

Submitted 1 November, 2023; originally announced November 2023.

Comments: 191 pages, 47 figures. This is the version of the article before peer review or editing, as submitted by an author to J. Phys. Photonics. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at https://iopscience.iop.org/article/10.1088/2515-7647/ace869/meta

Journal ref: J. Phys. Photonics 5 042501 (2023)

Detecting Shared Data Manipulation in Distributed Optimization Algorithms

Authors: Mohannad Alkhraijah, Rachel Harris, Samuel Litchfield, David Huggins, Daniel K. Molzahn

Abstract: This paper investigates the vulnerability of the Alternating Direction Method of Multipliers (ADMM) algorithm to shared data manipulation, with a focus on solving optimal power flow (OPF) problems. Deliberate data manipulation may cause the ADMM algorithm to converge to suboptimal solutions. We derive two sufficient conditions for detecting data manipulation based on the theoretical convergence tr… ▽ More This paper investigates the vulnerability of the Alternating Direction Method of Multipliers (ADMM) algorithm to shared data manipulation, with a focus on solving optimal power flow (OPF) problems. Deliberate data manipulation may cause the ADMM algorithm to converge to suboptimal solutions. We derive two sufficient conditions for detecting data manipulation based on the theoretical convergence trajectory of the ADMM algorithm. We evaluate the detection conditions' performance on three data manipulation strategies we previously proposed: simple, feedback, and bilevel optimization attacks. We then extend these three data manipulation strategies to avoid detection by considering both the detection conditions and a neural network (NN) detection model in the attacks. We also propose an adversarial NN training framework to detect shared data manipulation. We illustrate the performance of our data manipulation strategy and detection framework on OPF problems. The results show that the proposed detection conditions successfully detect most of the data manipulation attacks. However, a bilevel optimization attack strategy that incorporates the detection methods may avoid being detected. Countering this, our proposed adversarial training framework detects all the instances of the bilevel optimization attack. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Spatio-temporal spread of COVID-19 and its associations with socioeconomic, demographic and environmental factors in England: A Bayesian hierarchical spatio-temporal model

Authors: Xueqing Yin, John M. Aiken, Richard Harris, Jonathan L. Bamber

Abstract: Exploring the spatio-temporal variations of COVID-19 transmission and its potential determinants could provide a deeper understanding of the dynamics of disease spread. This study aims to investigate the spatio-temporal spread of COVID-19 infection rate in England, and examine its associations with socioeconomic, demographic and environmental risk factors. Using weekly reported COVID-19 cases from… ▽ More Exploring the spatio-temporal variations of COVID-19 transmission and its potential determinants could provide a deeper understanding of the dynamics of disease spread. This study aims to investigate the spatio-temporal spread of COVID-19 infection rate in England, and examine its associations with socioeconomic, demographic and environmental risk factors. Using weekly reported COVID-19 cases from 7 March 2020 to 26 March 2022 at Middle Layer Super Output Area (MSOA) level in mainland England, we developed a Bayesian hierarchical spatio-temporal model to predict the COVID-19 infection rates and investigate the influencing factors. The analysis showed that our model outperformed the ordinary least squares (OLS) and geographically weighted regression (GWR) models in terms of prediction accuracy. The results showed that the spread of COVID-19 infection rates over space and time was heterogeneous. Hotspots of infection rate exhibited inconsistent clustered patterns over time. Among the selected risk factors, the annual household income, unemployment rate, population density, percentage of Caribbean population, percentage of adults aged 45-64 years old, and particulate matter concentrations were found to be positively associated with the COVID-19 infection rate. The findings assist policymakers in developing tailored public health interventions for COVID-19 prevention and control. △ Less

Submitted 18 August, 2023; originally announced August 2023.

Comments: 19 pages, 5 figures

Liquid-Crystal-Based Controllable Attenuators Operating in the 1-4 Terahertz Band

Authors: Aniela Dunn, Zhaopeng Zhang, Michael D. Horbury, Eleanor V. Nuttall, Yingjun Han, Mohammed Salih, Lianhe Li, Abigail Bond, Ehab Saleh, Russell Harris, Diego Pardo, Brian N. Ellison, Andrew D. Burnett, Helen F. Gleeson, Alexander Valavanis

Abstract: Liquid-crystal devices (LCDs) offer a potential route toward adaptive optical components for use in the < 2 THz band of the electromagnetic spectrum. We demonstrate LCDs using a commercially available material (E7), with unbiased birefringence values of 0.14-0.18 in the 0.3-4 THz band. We exploit the linear dichroism of the material to modulate the emission from a 3.4-THz quantum cascade laser by… ▽ More Liquid-crystal devices (LCDs) offer a potential route toward adaptive optical components for use in the < 2 THz band of the electromagnetic spectrum. We demonstrate LCDs using a commercially available material (E7), with unbiased birefringence values of 0.14-0.18 in the 0.3-4 THz band. We exploit the linear dichroism of the material to modulate the emission from a 3.4-THz quantum cascade laser by up to 40%, dependent upon both the liquid-crystal layer thickness and the bias voltage applied. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: 10 pages, 6 figures

A tale of two faults: Statistical reconstruction of the 1820 Flores Sea earthquake using tsunami observations alone

Authors: T. Paskett, J. P. Whitehead, R. A. Harris, C. Ashcroft, J. A. Krometis, I. Sorensen, R. Wonnacott

Abstract: Using a Bayesian approach we compare anecdotal tsunami runup observations from the 29 December 1820 Flores Sea earthquake with close to 200,000 tsunami simulations to determine the most probable earthquake parameters causing the tsunami. Using a dual hypothesis of the source earthquake either originating from the Flores Thrust or the Walanae/Selayar Fault, we found that neither source perfectly ma… ▽ More Using a Bayesian approach we compare anecdotal tsunami runup observations from the 29 December 1820 Flores Sea earthquake with close to 200,000 tsunami simulations to determine the most probable earthquake parameters causing the tsunami. Using a dual hypothesis of the source earthquake either originating from the Flores Thrust or the Walanae/Selayar Fault, we found that neither source perfectly matches the observational data, particularly while satisfying seismic constraints of the region. However, there is clear quantitative evidence that a major earthquake on the Walanae/Selayar Fault more closely aligns with historical records of the tsunami, and earthquake shaking. The simulated data available from this study alludes to the potential for a different source in the region or the occurrence of an earthquake near where both faults potentially merge and simultaneously rupture similar to the 2016 Kaikoura, New Zealand event. △ Less

Submitted 2 May, 2023; originally announced May 2023.

Comments: 22 pages, 9 figures

PowerModelsADA: A Framework for Solving Optimal Power Flow using Distributed Algorithms

Authors: Mohannad Alkhraijah, Rachel Harris, Carleton Coffrin, Daniel K. Molzahn

Abstract: This paper presents PowerModelsADA, an open-source framework for solving Optimal Power Flow (OPF) problems using Alternating Distributed Algorithms (ADA). PowerModelsADA provides a framework to test, verify, and benchmark both existing and new ADAs. This paper demonstrates use cases for PowerModelsADA and validates its implementation with multiple OPF formulations. This paper presents PowerModelsADA, an open-source framework for solving Optimal Power Flow (OPF) problems using Alternating Distributed Algorithms (ADA). PowerModelsADA provides a framework to test, verify, and benchmark both existing and new ADAs. This paper demonstrates use cases for PowerModelsADA and validates its implementation with multiple OPF formulations. △ Less

Submitted 5 October, 2023; v1 submitted 2 April, 2023; originally announced April 2023.

Comments: This work has been submitted to the IEEE for possible publication

Experimental verification of the temperature coefficient of resistivity

Authors: Robert D. Polak, Michael R. Harris, Kiet A. Nguyen, Anthony Kearns

Abstract: We have created an experimental procedure for determining the temperature coefficient of resistivity, $α_R$, for introductory physics laboratories. This method examines the relationship between temperature and resistivity to establish $α_R$ within 10% of the accepted value. We have created an experimental procedure for determining the temperature coefficient of resistivity, $α_R$, for introductory physics laboratories. This method examines the relationship between temperature and resistivity to establish $α_R$ within 10% of the accepted value. △ Less

Submitted 21 November, 2022; originally announced January 2023.

Determination of the coefficient of thermal expansion by measuring frequency of a heated music wire

Authors: Robert D. Polak, Michael R. Harris, Kiet A. Nguyen, Anthony Kearns

Abstract: Engaging with physical and material properties through empirical observation is a fundamental part of undergraduate physics and engineering education. Several works have proposed experiments to determine thermal physical constants of materials such as finding the coefficient of linear expansion. As Dajbych and Polak et al. have shown, methods for experimentally verifying physical constants can be… ▽ More Engaging with physical and material properties through empirical observation is a fundamental part of undergraduate physics and engineering education. Several works have proposed experiments to determine thermal physical constants of materials such as finding the coefficient of linear expansion. As Dajbych and Polak et al. have shown, methods for experimentally verifying physical constants can be done by measuring the frequency of a plucked high-carbon steel wire on a guitar. Building upon our previous work, we have extended our method to verify the coefficient of linear thermal expansion through an accessible procedure directed at introductory physics education. △ Less

Submitted 19 November, 2022; originally announced November 2022.

Report of the Topical Group on Physics Beyond the Standard Model at Energy Frontier for Snowmass 2021

Authors: Tulika Bose, Antonio Boveia, Caterina Doglioni, Simone Pagan Griso, James Hirschauer, Elliot Lipeles, Zhen Liu, Nausheen R. Shah, Lian-Tao Wang, Kaustubh Agashe, Juliette Alimena, Sebastian Baum, Mohamed Berkat, Kevin Black, Gwen Gardner, Tony Gherghetta, Josh Greaves, Maxx Haehn, Phil C. Harris, Robert Harris, Julie Hogan, Suneth Jayawardana, Abraham Kahn, Jan Kalinowski, Simon Knapen , et al. (297 additional authors not shown)

Abstract: This is the Snowmass2021 Energy Frontier (EF) Beyond the Standard Model (BSM) report. It combines the EF topical group reports of EF08 (Model-specific explorations), EF09 (More general explorations), and EF10 (Dark Matter at Colliders). The report includes a general introduction to BSM motivations and the comparative prospects for proposed future experiments for a broad range of potential BSM mode… ▽ More This is the Snowmass2021 Energy Frontier (EF) Beyond the Standard Model (BSM) report. It combines the EF topical group reports of EF08 (Model-specific explorations), EF09 (More general explorations), and EF10 (Dark Matter at Colliders). The report includes a general introduction to BSM motivations and the comparative prospects for proposed future experiments for a broad range of potential BSM models and signatures, including compositeness, SUSY, leptoquarks, more general new bosons and fermions, long-lived particles, dark matter, charged-lepton flavor violation, and anomaly detection. △ Less

Submitted 18 October, 2022; v1 submitted 26 September, 2022; originally announced September 2022.

Comments: 108 pages + 38 pages references and appendix, 37 figures, Report of the Topical Group on Beyond the Standard Model Physics at Energy Frontier for Snowmass 2021. The first nine authors are the Conveners, with Contributions from the other authors

Quantum critical dynamics in a 5000-qubit programmable spin glass

Authors: Andrew D. King, Jack Raymond, Trevor Lanting, Richard Harris, Alex Zucca, Fabio Altomare, Andrew J. Berkley, Kelly Boothby, Sara Ejtemaee, Colin Enderud, Emile Hoskinson, Shuiyuan Huang, Eric Ladizinsky, Allison J. R. MacDonald, Gaelen Marsden, Reza Molavi, Travis Oh, Gabriel Poulin-Lamarre, Mauricio Reis, Chris Rich, Yuki Sato, Nicholas Tsai, Mark Volkmann, Jed D. Whittaker, Jason Yao , et al. (2 additional authors not shown)

Abstract: Experiments on disordered alloys suggest that spin glasses can be brought into low-energy states faster by annealing quantum fluctuations than by conventional thermal annealing. Due to the importance of spin glasses as a paradigmatic computational testbed, reproducing this phenomenon in a programmable system has remained a central challenge in quantum optimization. Here we achieve this goal by rea… ▽ More Experiments on disordered alloys suggest that spin glasses can be brought into low-energy states faster by annealing quantum fluctuations than by conventional thermal annealing. Due to the importance of spin glasses as a paradigmatic computational testbed, reproducing this phenomenon in a programmable system has remained a central challenge in quantum optimization. Here we achieve this goal by realizing quantum critical spin-glass dynamics on thousands of qubits with a superconducting quantum annealer. We first demonstrate quantitative agreement between quantum annealing and time-evolution of the Schrödinger equation in small spin glasses. We then measure dynamics in 3D spin glasses on thousands of qubits, where simulation of many-body quantum dynamics is intractable. We extract critical exponents that clearly distinguish quantum annealing from the slower stochastic dynamics of analogous Monte Carlo algorithms, providing both theoretical and experimental support for a scaling advantage in reducing energy as a function of annealing time. △ Less

Submitted 18 April, 2023; v1 submitted 27 July, 2022; originally announced July 2022.

Journal ref: Nature, 2023

Summarizing experimental sensitivities of collider experiments to dark matter models and comparison to other experiments

Authors: Antonio Boveia, Caterina Doglioni, Boyu Gao, Josh Greaves, Philip Harris, Katherine Pachal, Etienne Dreyer, Giuliano Gustavino, Robert Harris, Daniel Hayden, Tetiana Hrynova, Ashutosh Kotwal, Jared Little, Kevin Black, Tulika Bose, Yuze Chen, Sridhara Dasu, Haoyi Jia, Deborah Pinna, Varun Sharma, Nikhilesh Venkatasubramanian, Carl Vuosalo

Abstract: Comparisons of the coverage of current and proposed dark matter searches can help us to understand the context in which a discovery of particle dark matter would be made. In some scenarios, a discovery could be reinforced by information from multiple, complementary types of experiments; in others, only one experiment would see a signal, giving only a partial, more ambiguous picture; in still other… ▽ More Comparisons of the coverage of current and proposed dark matter searches can help us to understand the context in which a discovery of particle dark matter would be made. In some scenarios, a discovery could be reinforced by information from multiple, complementary types of experiments; in others, only one experiment would see a signal, giving only a partial, more ambiguous picture; in still others, no experiment would be sensitive and new approaches would be needed. In this whitepaper, we present an update to a similar study performed for the European Strategy Briefing Book performed within the dark matter at the Energy Frontier (EF10) Snowmass Topical Group We take as a starting point a set of projections for future collider facilities and a method of graphical comparisons routinely performed for LHC DM searches using simplified models recommended by the LHC Dark Matter Working Group and also used for the BSM and dark matter chapters of the European Strategy Briefing Book. These comparisons can also serve as launching point for cross-frontier discussions about dark matter complementarity. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: Submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021)

Mass distributions of dijet resonances from excited quarks at proton-proton colliders

Authors: Emine Gurpinar Guler, Yalcin Guler, Robert M. Harris

Abstract: We study the expected experimental mass distributions of dijet resonances from excited quarks in proton-proton collisions at energies $\sqrt{s}=$ 13, 14, 27, 100, 300, and 500 TeV. We explore in detail the expected shapes at both the generator and experimental levels, and identify within the distributions the effects of the excited quark natural width, parton momentum distributions of the proton,… ▽ More We study the expected experimental mass distributions of dijet resonances from excited quarks in proton-proton collisions at energies $\sqrt{s}=$ 13, 14, 27, 100, 300, and 500 TeV. We explore in detail the expected shapes at both the generator and experimental levels, and identify within the distributions the effects of the excited quark natural width, parton momentum distributions of the proton, radiation, and experimental resolution. We present both differential and cumulative probability distributions as a function of dijet mass, and the signal acceptance of a window in dijet mass centered on each resonance. We find that for a range of resonance masses, between 10\% and 50\% of $\sqrt{s}$, the dijet mass distributions and window acceptance are practically universal, approximately invariant under changes in resonance mass and $\sqrt{s}$. This work supports our Snowmass 2021 study on the sensitivity to dijet resonances at proton-proton colliders. △ Less

Submitted 5 May, 2022; originally announced May 2022.

Comments: 24 pages, 12 figures, to be submitted to JHEP

Report number: FERMILAB-PUB-22-301-PPD-SCD

All-fibre wavefront sensor

Authors: Thomas A. Wright, Stephanos Yerolatsitis, Kerrianne Harrington, Robert J. Harris, Tim A. Birks

Abstract: We report on a tapered three-core optical fibre that can be used as a tip-tilt wavefront sensor. In this device, a coupled region of a few millimetres at the sensing tip of the fibre converts fragile phase information from an incoming wavefront into robust intensity information within each of the cores. The intensity information can be easily converted to linear wavefront error over small ranges,… ▽ More We report on a tapered three-core optical fibre that can be used as a tip-tilt wavefront sensor. In this device, a coupled region of a few millimetres at the sensing tip of the fibre converts fragile phase information from an incoming wavefront into robust intensity information within each of the cores. The intensity information can be easily converted to linear wavefront error over small ranges, making it ideal for closed loop systems. The sensor uses minimal information to infer tip-tilt and is compatible with remote detector arrays. We explore its application within adaptive optics and present a validation case to show its applicability to astronomy. △ Less

Submitted 25 July, 2022; v1 submitted 29 April, 2022; originally announced April 2022.

Comments: 7 pages, 6 figures

Journal ref: MNRAS Volume 514, Issue 4, August 2022, Pages 5422-5428

The Future Circular Collider: a Summary for the US 2021 Snowmass Process

Authors: G. Bernardi, E. Brost, D. Denisov, G. Landsberg, M. Aleksa, D. d'Enterria, P. Janot, M. L. Mangano, M. Selvaggi, F. Zimmermann, J. Alcaraz Maestre, C. Grojean, R. M. Harris, A. Pich, M. Vos, S. Heinemeyer, P. Giacomelli, P. Azzi, F. Bedeschi, M. Klute, A. Blondel, C. Paus, F. Simon, M. Dam, E. Barberis , et al. (19 additional authors not shown)

Abstract: In this white paper for the 2021 Snowmass process, we give a description of the proposed Future Circular Collider (FCC) project and its physics program. The paper summarizes and updates the discussion submitted to the European Strategy on Particle Physics. After construction of an approximately 90 km tunnel, an electron-positron collider based on established technologies allows world-record instan… ▽ More In this white paper for the 2021 Snowmass process, we give a description of the proposed Future Circular Collider (FCC) project and its physics program. The paper summarizes and updates the discussion submitted to the European Strategy on Particle Physics. After construction of an approximately 90 km tunnel, an electron-positron collider based on established technologies allows world-record instantaneous luminosities at center-of-mass energies from the Z resonance up to tt thresholds, enabling a rich set of fundamental measurements including Higgs couplings determinations at the sub percent level, precision tests of the weak and strong forces, and searches for new particles, including dark matter, both directly and via virtual corrections or mixing. Among other possibilities, the FCC-ee will be able to (i) indirectly discover new particles coupling to the Higgs and/or electroweak bosons up to scales around 7 and 50 TeV, respectively; (ii) perform competitive SUSY tests at the loop level in regions not accessible at the LHC; (iii) study heavy-flavor and tau physics in ultra-rare decays beyond the LHC reach, and (iv) achieve the best potential in direct collider searches for dark matter, sterile neutrinos, and axion-like particles with masses up to around 90 GeV. The tunnel can then be reused for a proton-proton collider, establishing record center-of-mass collision energy, allowing unprecedented reach for direct searches for new particles up to the around 50 TeV scale, and a diverse program of measurements of the Standard Model and Higgs boson, including a precision measurement of the Higgs self-coupling, and conclusively testing weakly-interacting massive particle scenarios of thermal relic dark matter. △ Less

Submitted 19 December, 2022; v1 submitted 12 March, 2022; originally announced March 2022.

Comments: 84 pages

MSC Class: for Snowmass 2021

Coherent quantum annealing in a programmable 2000-qubit Ising chain

Authors: Andrew D. King, Sei Suzuki, Jack Raymond, Alex Zucca, Trevor Lanting, Fabio Altomare, Andrew J. Berkley, Sara Ejtemaee, Emile Hoskinson, Shuiyuan Huang, Eric Ladizinsky, Allison MacDonald, Gaelen Marsden, Travis Oh, Gabriel Poulin-Lamarre, Mauricio Reis, Chris Rich, Yuki Sato, Jed D. Whittaker, Jason Yao, Richard Harris, Daniel A. Lidar, Hidetoshi Nishimori, Mohammad H. Amin

Abstract: Quantum simulation has emerged as a valuable arena for demonstrating and understanding the capabilities of near-term quantum computers. Quantum annealing has been used successfully in simulating a range of open quantum systems, both at equilibrium and out of equilibrium. However, in all previous experiments, annealing has been too slow to simulate a closed quantum system coherently, due to the ons… ▽ More Quantum simulation has emerged as a valuable arena for demonstrating and understanding the capabilities of near-term quantum computers. Quantum annealing has been used successfully in simulating a range of open quantum systems, both at equilibrium and out of equilibrium. However, in all previous experiments, annealing has been too slow to simulate a closed quantum system coherently, due to the onset of thermal effects from the environment. Here we demonstrate coherent evolution through a quantum phase transition in the paradigmatic setting of the 1D transverse-field Ising chain, using up to 2000 superconducting flux qubits in a programmable quantum annealer. In large systems we observe the quantum Kibble-Zurek mechanism with theoretically predicted kink statistics, as well as characteristic positive kink-kink correlations, independent of system temperature. In small chains, excitation statistics validate the picture of a Landau-Zener transition at a minimum gap. In both cases, results are in quantitative agreement with analytical solutions to the closed-system quantum model. For slower anneals we observe anti-Kibble-Zurek scaling in a crossover to the open quantum regime. These experiments demonstrate that large-scale quantum annealers can be operated coherently, paving the way to exploiting coherent dynamics in quantum optimization, machine learning, and simulation tasks. △ Less

Submitted 10 March, 2022; v1 submitted 11 February, 2022; originally announced February 2022.

Journal ref: Nature Physics Vol. 18, 1324-1328 (2022)

Sensitivity to Dijet Resonances at Proton-Proton Colliders

Authors: Robert M. Harris, Emine Gurpinar Guler, Yalcin Guler

Abstract: A significant benchmark for discovery at a proton-proton collider is the sensitivity to a dijet resonance, X, the intermediate state of the s-channel process $pp \rightarrow X \rightarrow 2\mbox{ jets}$. To probe the highest resonance masses, hadron collider experiments have used the classic technique of searching for bumps in the mass spectrum of two individually resolved jets. In this Snowmass 2… ▽ More A significant benchmark for discovery at a proton-proton collider is the sensitivity to a dijet resonance, X, the intermediate state of the s-channel process $pp \rightarrow X \rightarrow 2\mbox{ jets}$. To probe the highest resonance masses, hadron collider experiments have used the classic technique of searching for bumps in the mass spectrum of two individually resolved jets. In this Snowmass 2021 study, we explore the search sensitivity to multiple benchmark models of dijet resonances at current and future proton-proton colliders. We present the expected masses for $5σ$ discovery and 95\% confidence level exclusion of diquarks, colorons, excited quarks, $W^{\prime}s$, $Z^{\prime}s$ and Randall-Sundrum gravitons, resulting from accumulation of integrated luminosities between 10 and $10^5$ fb$^{-1}$, at proton-proton colliders operating at energies $\sqrt{s}=$ 13, 14, 27, 75, 100, 150, 300 and 500 TeV. △ Less

Submitted 7 February, 2022; originally announced February 2022.

Comments: Contribution to Snowmass 2021. Containing 22 pages, 11 Figures, 4 Tables

Report number: FERMILAB-FN-1143-PPD-SCD

Expressing and Executing Informed Consent Permissions Using SWRL: The All of Us Use Case

Authors: Muhammad Amith, Marcelline R. Harris, Cooper Stansbury, Kathleen Ford, Frank J. Manion, Cui Tao

Abstract: The informed consent process is a complicated procedure involving permissions as well a variety of entities and actions. In this paper, we discuss the use of Semantic Web Rule Language (SWRL) to further extend the Informed Consent Ontology (ICO) to allow for semantic machine-based reasoning to manage and generate important permission-based information that can later be viewed by stakeholders. We p… ▽ More The informed consent process is a complicated procedure involving permissions as well a variety of entities and actions. In this paper, we discuss the use of Semantic Web Rule Language (SWRL) to further extend the Informed Consent Ontology (ICO) to allow for semantic machine-based reasoning to manage and generate important permission-based information that can later be viewed by stakeholders. We present four use cases of permissions from the All of Us informed consent document and translate these permissions into SWRL expressions to extend and operationalize ICO. Our efforts show how SWRL is able to infer some of the implicit information based on the defined rules, and demonstrate the utility of ICO through the use of SWRL extensions. Future work will include developing formal and generalized rules and expressing permissions from the entire document, as well as working towards integrating ICO into software systems to enhance the semantic representation of informed consent for biomedical research. △ Less

Submitted 23 August, 2021; originally announced August 2021.

Comments: To appear in: Proceedings of the American Medical Informatics Associations (AMIA) 2021 Annual Symposium Oct 30-Nov 03 San Diego, CA, USA. Please visit and cite the canonical version once available. M Amith and M Harris contributed equally to this work

ACM Class: F.4; H.4; E.2; I.7; H.5; A.m

Long-term stability of planets in and around binary stars

Authors: Harry A. Ballantyne, Tore Espaas, Bethan Z. Norgrove, Bethany A. Wootton, Benjamin R. Harris, Isaac L. Pepper, Richard D. Smith, Rosie E. Dommett, Richard J. Parker

Abstract: Planets are observed to orbit the component star(s) of stellar binary systems on so-called circumprimary or circumsecondary orbits, as well as around the entire binary system on so-called circumbinary orbits. Depending on the orbital parameters of the binary system a planet will be dynamically stable if it orbits within some critical separation of the semimajor axis in the circumprimary case, or b… ▽ More Planets are observed to orbit the component star(s) of stellar binary systems on so-called circumprimary or circumsecondary orbits, as well as around the entire binary system on so-called circumbinary orbits. Depending on the orbital parameters of the binary system a planet will be dynamically stable if it orbits within some critical separation of the semimajor axis in the circumprimary case, or beyond some critical separation for the circumbinary case. We present N-body simulations of star-forming regions that contain populations of primordial binaries to determine the fraction of binary systems that can host stable planets at various semimajor axes, and how this fraction of stable systems evolves over time. Dynamical encounters in star-forming regions can alter the orbits of some binary systems, which can induce long-term dynamical instabilities in the planetary system and can even change the size of the habitable zone(s) of the component stars. However, the overall fraction of binaries that can host stable planetary systems is not greatly affected by either the assumed binary population, or the density of the star-forming region. Instead, the critical factor in determining how many stable planetary systems exist in the Galaxy is the stellar binary fraction - the more stars that are born as singles in stellar nurseries, the higher the fraction of stable planetary systems. △ Less

Submitted 17 August, 2021; originally announced August 2021.

Comments: 13 pages, 7 figures, 2 short appendices, accepted for publication in MNRAS. A video abstract created by Helena Gibbon is available here: https://youtu.be/76kANnTK9-s

Decision-making with distorted memory: Escaping the trap of past experience

Authors: Evangelos Mitsokapas, Rosemary J. Harris

Abstract: Snapshots of "best" (or "worst") experience are known to dominate human memory and may thus also have a significant effect on future behaviour. We consider here a model of repeated decision-making where, at every time step, an agent takes one of two choices with probabilities which are functions of the maximum utilities previously experienced. Depending on the utility distributions and the level o… ▽ More Snapshots of "best" (or "worst") experience are known to dominate human memory and may thus also have a significant effect on future behaviour. We consider here a model of repeated decision-making where, at every time step, an agent takes one of two choices with probabilities which are functions of the maximum utilities previously experienced. Depending on the utility distributions and the level of noise in the decision process, it is possible for an agent to become "trapped" in one of the choices on the basis of their early experiences. If the utility distributions for the two choices are different, then the agent may even become trapped in the choice which is objectively worse in the sense of expected long-term returns; crucially we extend earlier work to address this case. Using tools from statistical physics and extreme-value theory, we show that for exponential utilities there is an optimal value of noise which maximizes the expected returns in the long run. We also briefly discuss the behaviour for other utility distributions. △ Less

Submitted 17 February, 2022; v1 submitted 12 August, 2021; originally announced August 2021.

Comments: 28 pages, 12 figures. Minor revisions. This version very similar to that published in Physica A: Statistical Mechanics and its Applications

Journal ref: Physica A: Statistical Mechanics and its Applications, 593 (2022) 126762

Architectural considerations in the design of a third-generation superconducting quantum annealing processor

Authors: Kelly Boothby, Colin Enderud, Trevor Lanting, Reza Molavi, Nicholas Tsai, Mark H. Volkmann, Fabio Altomare, Mohammad H. Amin, Michael Babcock, Andrew J. Berkley, Catia Baron Aznar, Martin Boschnak, Holly Christiani, Sara Ejtemaee, Bram Evert, Matthew Gullen, Markus Hager, Richard Harris, Emile Hoskinson, Jeremy P. Hilton, Kais Jooya, Ann Huang, Mark W. Johnson, Andrew D. King, Eric Ladizinsky , et al. (24 additional authors not shown)

Abstract: Early generations of superconducting quantum annealing processors have provided a valuable platform for studying the performance of a scalable quantum computing technology. These studies have directly informed our approach to the design of the next-generation processor. Our design priorities for this generation include an increase in per-qubit connectivity, a problem Hamiltonian energy scale simil… ▽ More Early generations of superconducting quantum annealing processors have provided a valuable platform for studying the performance of a scalable quantum computing technology. These studies have directly informed our approach to the design of the next-generation processor. Our design priorities for this generation include an increase in per-qubit connectivity, a problem Hamiltonian energy scale similar to previous generations, reduced Hamiltonian specification errors, and an increase in the processor scale that also leaves programming and readout times fixed or reduced. Here we discuss the specific innovations that resulted in a processor architecture that satisfies these design priorities. △ Less

Submitted 4 August, 2021; originally announced August 2021.

On-sky results for the novel integrated micro-lens ring tip-tilt sensor

Authors: Philipp Hottinger, Robert J. Harris, Jonathan Crass, Philipp-Immanuel Dietrich, Matthias Blaicher, Andrew Bechter, Brian Sands, Tim J. Morris, Alastair G. Basden, Nazim Ali Bharmal, Jochen Heidt, Theodoros Anagnos, Philip L. Neureuther, Martin Glück, Jennifer Power, Jörg-Uwe Pott, Christian Koos, Oliver Sawodny, Andreas Quirrenbach

Abstract: We present the first on-sky results of the micro-lens ring tip-tilt (MLR-TT) sensor. This sensor utilizes a 3D printed micro-lens ring feeding six multi-mode fibers to sense misaligned light, allowing centroid reconstruction. A tip-tilt mirror allows the beam to be corrected, increasing the amount of light coupled into a centrally positioned single-mode (science) fiber. The sensor was tested with… ▽ More We present the first on-sky results of the micro-lens ring tip-tilt (MLR-TT) sensor. This sensor utilizes a 3D printed micro-lens ring feeding six multi-mode fibers to sense misaligned light, allowing centroid reconstruction. A tip-tilt mirror allows the beam to be corrected, increasing the amount of light coupled into a centrally positioned single-mode (science) fiber. The sensor was tested with the iLocater acquisition camera at the Large Binocular Telescope in November 2019. The limit on the maximum achieved root mean square reconstruction accuracy was found to be 0.19 $λ$/D in both tip and tilt, of which approximately 50% of the power originates at frequencies below 10 Hz. We show the reconstruction accuracy is highly dependent on the estimated Strehl ratio and simulations support the assumption that residual adaptive optics aberrations are the main limit to the reconstruction accuracy. We conclude that this sensor is ideally suited to remove post-adaptive optics non-common path tip tilt residuals. We discuss the next steps for the concept development, including optimizations of the lens and fiber, tuning of the correction algorithm and selection of optimal science cases. △ Less

Submitted 18 June, 2021; originally announced June 2021.

Comments: Accepted for publication in Journal of the Optical Society of America B (JOSA B)

Embracing Uncertainty in "Small Data" Problems: Estimating Earthquakes from Historical Anecdotes

Authors: Justin A. Krometis, Hayden Ringer, Jared P. Whitehead, Nathan E. Glatt-Holtz, Ronald A. Harris, Andrew J. Holbrook

Abstract: Improving understanding of current seismic risk is often dependent on developing a more complete characterization of earthquakes that have occurred in the past, and in particular before the development of modern sensing equipment in the middle of the twentieth century. However, accounts of such events are typically anecdotal in nature, limiting efforts to model them to more heuristic approaches. T… ▽ More Improving understanding of current seismic risk is often dependent on developing a more complete characterization of earthquakes that have occurred in the past, and in particular before the development of modern sensing equipment in the middle of the twentieth century. However, accounts of such events are typically anecdotal in nature, limiting efforts to model them to more heuristic approaches. To address this shortfall, we develop a framework based in Bayesian inference to provide a more rigorous methodology for estimating pre-instrumental earthquakes. By directly modeling accounts of resultant tsunamis via probability distributions, the framework allows practitioners to make principled estimates of key characteristics (e.g., magnitude and location) of historical earthquakes. To illustrate this idea, we apply the methodology to the estimation of an earthquake in Eastern Indonesia in the mid 19th century, the source of which is currently the subject of considerable debate in the geological community. The approach taken here gives evidence that even "small data" that is limited in scope and extremely uncertain can still be used to yield information on past seismic events. Moreover, sensitivity bounds indicate that the results obtained here are robust despite the inherent uncertainty in the observations. △ Less

Submitted 4 October, 2022; v1 submitted 14 June, 2021; originally announced June 2021.

MSC Class: 62P35; 86A22

3D-M3: High-spatial resolution spectroscopy with extreme AO and 3D printed micro-lenslets

Authors: Theodoros Anagnos, Mareike Trappen, Blaise C. Kuo Tiong, Tobias Feger, Stephanos Yerolatsitis, Robert J. Harris, Julien Lozi, Nemanja Jovanovic, Tim A. Birks, Sébastien Vievard, Olivier Guyon, Itandehui Gris-Sánchez, Sergio G. Leon-Saval, Barnaby Norris, Sebastiaan Y. Haffert, Phillip Hottinger, Matthias Blaicher, Yilin Xu, Christopher H. Betters, Christian Koos, David W. Coutts, Christian Schwab, Andreas Quirrenbach

Abstract: By combining IFS with ExAO we are now able to resolve objects close to the diffraction-limit of large telescopes, exploring new science cases. We introduce an IFU designed to couple light with a minimal platescale from the SCExAO facility at NIR wavelengths to a SM spectrograph. The IFU has a 3D-printed MLA on top of a custom SM MCF, to optimize the coupling of light into the fiber cores. We demon… ▽ More By combining IFS with ExAO we are now able to resolve objects close to the diffraction-limit of large telescopes, exploring new science cases. We introduce an IFU designed to couple light with a minimal platescale from the SCExAO facility at NIR wavelengths to a SM spectrograph. The IFU has a 3D-printed MLA on top of a custom SM MCF, to optimize the coupling of light into the fiber cores. We demonstrate the potential of the instrument via initial results from the first on-sky runs at the 8.2 m Subaru Telescope with a spectrograph using off-the-shelf optics, allowing for rapid development with low cost. △ Less

Submitted 14 June, 2021; v1 submitted 12 May, 2021; originally announced May 2021.

Hospitalisation risk for COVID-19 patients infected with SARS-CoV-2 variant B.1.1.7: cohort analysis

Authors: Tommy Nyberg, Katherine A. Twohig, Ross J. Harris, Shaun R. Seaman, Joe Flannagan, Hester Allen, Andre Charlett, Daniela De Angelis, Gavin Dabrera, Anne M. Presanis

Abstract: Objective: To evaluate the relationship between coronavirus disease 2019 (COVID-19) diagnosis with SARS-CoV-2 variant B.1.1.7 (also known as Variant of Concern 202012/01) and the risk of hospitalisation compared to diagnosis with wildtype SARS-CoV-2 variants. Design: Retrospective cohort, analysed using stratified Cox regression. Setting: Community-based SARS-CoV-2 testing in England, individu… ▽ More Objective: To evaluate the relationship between coronavirus disease 2019 (COVID-19) diagnosis with SARS-CoV-2 variant B.1.1.7 (also known as Variant of Concern 202012/01) and the risk of hospitalisation compared to diagnosis with wildtype SARS-CoV-2 variants. Design: Retrospective cohort, analysed using stratified Cox regression. Setting: Community-based SARS-CoV-2 testing in England, individually linked with hospitalisation data. Participants: 839,278 laboratory-confirmed COVID-19 patients, of whom 36,233 had been hospitalised within 14 days, tested between 23rd November 2020 and 31st January 2021 and analysed at a laboratory with an available TaqPath assay that enables assessment of S-gene target failure (SGTF). SGTF is a proxy test for the B.1.1.7 variant. Patient data were stratified by age, sex, ethnicity, deprivation, region of residence, and date of positive test. Main outcome measures: Hospitalisation between 1 and 14 days after the first positive SARS-CoV-2 test. Results: 27,710 of 592,409 SGTF patients (4.7%) and 8,523 of 246,869 non-SGTF patients (3.5%) had been hospitalised within 1-14 days. The stratum-adjusted hazard ratio (HR) of hospitalisation was 1.52 (95% confidence interval [CI] 1.47 to 1.57) for COVID-19 patients infected with SGTF variants, compared to those infected with non-SGTF variants. The effect was modified by age (P<0.001), with HRs of 0.93-1.21 for SGTF compared to non-SGTF patients below age 20 years, 1.29 in those aged 20-29, and 1.45-1.65 in age groups 30 years or older. Conclusions: The results suggest that the risk of hospitalisation is higher for individuals infected with the B.1.1.7 variant compared to wildtype SARS-CoV-2, likely reflecting a more severe disease. The higher severity may be specific to adults above the age of 30. △ Less

Submitted 29 May, 2021; v1 submitted 12 April, 2021; originally announced April 2021.

Comments: 29 pages, 9 figures

Journal ref: BMJ 2021;373:n1412

Statistical Characterization of Airplane Delays

Authors: Evangelos Mitsokapas, Benjamin Schäfer, Rosemary J. Harris, Christian Beck

Abstract: The aviation industry is of great importance for a globally connected economy. Customer satisfaction with airlines and airport performance is considerably influenced by how much flights are delayed. But how should the delay be quantified with thousands of flights for each airport and airline? Here, we present a statistical analysis of arrival delays at several UK airports between 2018 and 2020. We… ▽ More The aviation industry is of great importance for a globally connected economy. Customer satisfaction with airlines and airport performance is considerably influenced by how much flights are delayed. But how should the delay be quantified with thousands of flights for each airport and airline? Here, we present a statistical analysis of arrival delays at several UK airports between 2018 and 2020. We establish a procedure to compare both mean delay and extreme events among airlines and airports, identifying a power-law decay of large delays. Furthermore, we note drastic changes in plane delay statistics during the COVID-19 pandemic. Finally, we find that delays are described by a superposition of simple distributions, leading to a superstatistics. △ Less

Submitted 31 March, 2021; v1 submitted 21 December, 2020; originally announced January 2021.

Comments: 14 pages, 10 figures

An innovative integral field unit upgrade with 3D-printed micro-lenses for the RHEA at Subaru

Authors: Theodoros Anagnos, Pascal Maier, Philipp Hottinger, Chris Betters, Tobias Feger, Sergio G. Leon-Saval, Itandehui Gris-Sánchez, Stephanos Yerolatsitis, Julien Lozi, Tim A. Birks, Sebastian Vievard, Nemanja Jovanovic, Adam D. Rains, Michael J. Ireland, Robert J. Harris, Blaise C. Kuo Tiong, Olivier Guyon, Barnaby Norris, Sebastiaan Y. Haffert, Matthias Blaicher, Yilin Xu, Moritz Straub, Jörg-Uwe Pott, Oliver Sawodny, Philip L. Neureuther , et al. (4 additional authors not shown)

Abstract: In the new era of Extremely Large Telescopes (ELTs) currently under construction, challenging requirements drive spectrograph designs towards techniques that efficiently use a facility's light collection power. Operating in the single-mode (SM) regime, close to the diffraction limit, reduces the footprint of the instrument compared to a conventional high-resolving power spectrograph. The custom bu… ▽ More In the new era of Extremely Large Telescopes (ELTs) currently under construction, challenging requirements drive spectrograph designs towards techniques that efficiently use a facility's light collection power. Operating in the single-mode (SM) regime, close to the diffraction limit, reduces the footprint of the instrument compared to a conventional high-resolving power spectrograph. The custom built injection fiber system with 3D-printed micro-lenses on top of it for the replicable high-resolution exoplanet and asteroseismology spectrograph at Subaru in combination with extreme adaptive optics of SCExAO, proved its high efficiency in a lab environment, manifesting up to ~77% of the theoretical predicted performance. △ Less

Submitted 24 January, 2021; originally announced January 2021.

Multi-Contrast Computed Tomography Healthy Kidney Atlas

Authors: Ho Hin Lee, Yucheng Tang, Kaiwen Xu, Shunxing Bao, Agnes B. Fogo, Raymond Harris, Mark P. de Caestecker, Mattias Heinrich, Jeffrey M. Spraggins, Yuankai Huo, Bennett A. Landman

Abstract: The construction of three-dimensional multi-modal tissue maps provides an opportunity to spur interdisciplinary innovations across temporal and spatial scales through information integration. While the preponderance of effort is allocated to the cellular level and explore the changes in cell interactions and organizations, contextualizing findings within organs and systems is essential to visualiz… ▽ More The construction of three-dimensional multi-modal tissue maps provides an opportunity to spur interdisciplinary innovations across temporal and spatial scales through information integration. While the preponderance of effort is allocated to the cellular level and explore the changes in cell interactions and organizations, contextualizing findings within organs and systems is essential to visualize and interpret higher resolution linkage across scales. There is a substantial normal variation of kidney morphometry and appearance across body size, sex, and imaging protocols in abdominal computed tomography (CT). A volumetric atlas framework is needed to integrate and visualize the variability across scales. However, there is no abdominal and retroperitoneal organs atlas framework for multi-contrast CT. Hence, we proposed a high-resolution CT retroperitoneal atlas specifically optimized for the kidney across non-contrast CT and early arterial, late arterial, venous and delayed contrast enhanced CT. Briefly, we introduce a deep learning-based volume of interest extraction method and an automated two-stage hierarchal registration pipeline to register abdominal volumes to a high-resolution CT atlas template. To generate and evaluate the atlas, multi-contrast modality CT scans of 500 subjects (without reported history of renal disease, age: 15-50 years, 250 males & 250 females) were processed. We demonstrate a stable generalizability of the atlas template for integrating the normal kidney variation from small to large, across contrast modalities and populations with great variability of demographics. The linkage of atlas and demographics provided a better understanding of the variation of kidney anatomy across populations. △ Less

Submitted 23 December, 2020; v1 submitted 22 December, 2020; originally announced December 2020.

Comments: 8 pages, 5 figures

Parallel decoding of multiple logical qubits in tensor-network codes

Authors: Terry Farrelly, Robert J. Harris, Nathan A. McMahon, Thomas M. Stace

Abstract: We consider tensor-network stabilizer codes and show that their tensor-network decoder has the property that independent logical qubits can be decoded in parallel. As long as the error rate is below threshold, we show that this parallel decoder is essentially optimal. As an application, we verify this for the max-rate holographic Steane (heptagon) code. For holographic codes this tensor-network de… ▽ More We consider tensor-network stabilizer codes and show that their tensor-network decoder has the property that independent logical qubits can be decoded in parallel. As long as the error rate is below threshold, we show that this parallel decoder is essentially optimal. As an application, we verify this for the max-rate holographic Steane (heptagon) code. For holographic codes this tensor-network decoder was shown to be efficient with complexity polynomial in n, the number of physical qubits. Here we show that, by using the parallel decoding scheme, the complexity is also linear in k, the number of logical qubits. Because the tensor-network contraction is computationally efficient, this allows us to exactly contract tensor networks corresponding to codes with up to half a million qubits. Finally, we calculate the bulk threshold (the threshold for logical qubits a fixed distance from the code centre) under depolarizing noise for the max-rate holographic Steane code to be 9.4%. △ Less

Submitted 14 December, 2020; originally announced December 2020.

Comments: 8 + 3 pages, 5 figures

Constraints on Lightly Ionizing Particles from CDMSlite

Authors: SuperCDMS Collaboration, I. Alkhatib, D. W. P. Amaral, T. Aralis, T. Aramaki, I. J. Arnquist, I. Ataee Langroudy, E. Azadbakht, S. Banik, D. Barker, C. Bathurst, D. A. Bauer, L. V. S. Bezerra, R. Bhattacharyya, M. A. Bowles, P. L. Brink, R. Bunker, B. Cabrera, R. Calkins, R. A. Cameron, C. Cartaro, D. G. Cerdeño, Y. -Y. Chang, M. Chaudhuri, R. Chen , et al. (93 additional authors not shown)

Abstract: The Cryogenic Dark Matter Search low ionization threshold experiment (CDMSlite) achieved efficient detection of very small recoil energies in its germanium target, resulting in sensitivity to Lightly Ionizing Particles (LIPs) in a previously unexplored region of charge, mass, and velocity parameter space. We report first direct-detection limits calculated using the optimum interval method on the v… ▽ More The Cryogenic Dark Matter Search low ionization threshold experiment (CDMSlite) achieved efficient detection of very small recoil energies in its germanium target, resulting in sensitivity to Lightly Ionizing Particles (LIPs) in a previously unexplored region of charge, mass, and velocity parameter space. We report first direct-detection limits calculated using the optimum interval method on the vertical intensity of cosmogenically-produced LIPs with an electric charge smaller than $e/(3\times10^5$), as well as the strongest limits for charge $\leq e/160$, with a minimum vertical intensity of $1.36\times10^{-7}$\,cm$^{-2}$s$^{-1}$sr$^{-1}$ at charge $e/160$. These results apply over a wide range of LIP masses (5\,MeV/$c^2$ to 100\,TeV/$c^2$) and cover a wide range of $βγ$ values (0.1 -- $10^6$), thus excluding non-relativistic LIPs with $βγ$ as small as 0.1 for the first time. △ Less

Submitted 19 February, 2022; v1 submitted 18 November, 2020; originally announced November 2020.

Journal ref: Phys. Rev. Lett. 127, 081802 (2021)

Four annular structures in a protostellar disk less than 500,000 years old

Authors: Dominique M. Segura-Cox, Anika Schmiedeke, Jaime E. Pineda, Ian W. Stephens, Manuel Fernández-López, Leslie W. Looney, Paola Caselli, Zhi-Yun Li, Lee G. Mundy, Woojin Kwon, Robert J. Harris

Abstract: Annular structures, or rings and gaps, in disks around pre-main sequence stars have been detected in abundance towards Class II objects ~1,000,000 years in age. These structures are often interpreted as evidence of planet formation, with planet-mass bodies carving rings and gaps in the disk. This implies that planet formation may already be underway in even younger disks in the Class I phase, when… ▽ More Annular structures, or rings and gaps, in disks around pre-main sequence stars have been detected in abundance towards Class II objects ~1,000,000 years in age. These structures are often interpreted as evidence of planet formation, with planet-mass bodies carving rings and gaps in the disk. This implies that planet formation may already be underway in even younger disks in the Class I phase, when the protostar is still embedded in a larger scale dense envelope of gas and dust. While younger disks likely play an important role in the onset of planet formation, only within the past decade have detailed properties of disks in the youngest star-forming phases begun to be observed. Here we present 1.3 mm dust emission observations with 5 au resolution that show four annular substructures in the disk of the young (<500,000 years) protostar IRS 63. IRS 63, a single Class I source located in the nearby Ophiuchus molecular cloud (at a distance of 144 pc), is one of the brightest Class I protostars at (sub)millimeter wavelengths that also has a relatively large disk (>50 au). Multiple annular substructures observed towards disks at young times can act as an early foothold for dust grain growth, a prerequisite of planet formation. Whether planets already exist or not in the disk of IRS 63, it is clear that the planet formation process begins in the young protostellar phases, earlier than predicted by current planet formation theories. △ Less

Submitted 7 October, 2020; originally announced October 2020.

Comments: Accepted for publication in Nature, appearing in the 8 Oct 2020 issue, DOI: 10.1038/s41586-020-2779-6

Methodological reconstruction of historical seismic events from anecdotal accounts of destructive tsunamis: a case study for the great 1852 Banda arc mega-thrust earthquake and tsunami

Authors: Hayden Ringer, Jared P. Whitehead, Justin Krometis, Ronald A. Harris, Nathan Glatt-Holtz, Spencer Giddens, Claire Ashcraft, Garret Carver, Adam Robertson, McKay Harward, Joshua Fullwood, Kameron Lightheart, Ryan Hilton, Ashley Avery, Cody Kesler, Martha Morrise, Michael Hunter Klein

Abstract: We demonstrate the efficacy of a Bayesian statistical inversion framework for reconstructing the likely characteristics of large pre-instrumentation earthquakes from historical records of tsunami observations. Our framework is designed and implemented for the estimation of the location and magnitude of seismic events from anecdotal accounts of tsunamis including shoreline wave arrival times, heigh… ▽ More We demonstrate the efficacy of a Bayesian statistical inversion framework for reconstructing the likely characteristics of large pre-instrumentation earthquakes from historical records of tsunami observations. Our framework is designed and implemented for the estimation of the location and magnitude of seismic events from anecdotal accounts of tsunamis including shoreline wave arrival times, heights, and inundation lengths over a variety of spatially separated observation locations. As an initial test case we use our framework to reconstruct the great 1852 earthquake and tsunami of eastern Indonesia. Relying on the assumption that these observations were produced by a subducting thrust event, the posterior distribution indicates that the observables were the result of a massive mega-thrust event with magnitude near 8.8 Mw and a likely rupture zone in the north-eastern Banda arc. The distribution of predicted epicentral locations overlaps with the largest major seismic gap in the region as indicated by instrumentally recorded seismic events. These results provide a geologic and seismic context for hazard risk assessment in coastal communities experiencing growing population and urbanization in Indonesia. In addition, the methodology demonstrated here highlights the potential for applying a Bayesian approach to enhance understanding of the seismic history of other subduction zones around the world. △ Less

Submitted 14 February, 2021; v1 submitted 29 September, 2020; originally announced September 2020.

Tensor-network codes

Authors: Terry Farrelly, Robert J. Harris, Nathan A. McMahon, Thomas M. Stace

Abstract: Inspired by holographic codes and tensor-network decoders, we introduce tensor-network stabilizer codes which come with a natural tensor-network decoder. These codes can correspond to any geometry, but, as a special case, we generalize holographic codes beyond those constructed from perfect or block-perfect isometries, and we give an example that corresponds to neither. Using the tensor-network de… ▽ More Inspired by holographic codes and tensor-network decoders, we introduce tensor-network stabilizer codes which come with a natural tensor-network decoder. These codes can correspond to any geometry, but, as a special case, we generalize holographic codes beyond those constructed from perfect or block-perfect isometries, and we give an example that corresponds to neither. Using the tensor-network decoder, we find a threshold of 18.8% for this code under depolarizing noise. We also show that for holographic codes the exact tensor-network decoder (with no bond-dimension truncation) is efficient with a complexity that is polynomial in the number of physical qubits, even for locally correlated noise. △ Less

Submitted 22 September, 2020; originally announced September 2020.

Comments: 5 pages main body + 5 pages appendices; 6 figures

Journal ref: Phys. Rev. Lett. 127, 040507 (2021)

Diffraction-limited integral-field spectroscopy for extreme adaptive optics systems with the Multi-Core fiber-fed Integral-Field Unit

Authors: Sebastiaan Y. Haffert, Robert J. Harris, Alessio Zanutta, Fraser A. Pike, Andrea Bianco, Eduardo Redaelli, Aurélien Benoît, David G. MacLachlan, Calum A. Ross, Itandehui Gris-Sánchez, Mareike D. Trappen, Yilin Xu, Matthias Blaicher, Pascal Maier, Giulio Riva, Baptiste Sinquin, Caroline Kulcsár, Nazim Ali Bharmal, Eric Gendron, Lazar Staykov, Tim J. Morris, Santiago Barboza, Norbert Muench, Lisa Bardou, Léonard Prengère , et al. (9 additional authors not shown)

Abstract: Direct imaging instruments have the spatial resolution to resolve exoplanets from their host star. This enables direct characterization of the exoplanets atmosphere, but most direct imaging instruments do not have spectrographs with high enough resolving power for detailed atmospheric characterization. We investigate the use of a single-mode diffraction-limited integral-field unit that is compact… ▽ More Direct imaging instruments have the spatial resolution to resolve exoplanets from their host star. This enables direct characterization of the exoplanets atmosphere, but most direct imaging instruments do not have spectrographs with high enough resolving power for detailed atmospheric characterization. We investigate the use of a single-mode diffraction-limited integral-field unit that is compact and easy to integrate into current and future direct imaging instruments for exoplanet characterization. This achieved by making use of recent progress in photonic manufacturing to create a single-mode fiber-fed image reformatter. The fiber-link is created with 3D printed lenses on top of a single-mode multi-core fiber that feeds an ultrafast laser inscribed photonic chip that reformats the fiber into a pseudo-slit. We then couple it to a first-order spectrograph with a triple stacked volume phase holographic grating for a high efficiency over a large bandwidth. The prototype system has had a successful first-light observing run at the 4.2 meter William Herschel Telescope. The measured on-sky resolving power is between 2500 and 3000, depending on the wavelength. With our observations we show that single-mode integral-field spectroscopy is a viable option for current and future exoplanet imaging instruments. △ Less

Submitted 8 September, 2020; originally announced September 2020.

Comments: submitted to JATIS

Decoding Holographic Codes with an Integer Optimisation Decoder

Authors: Robert J. Harris, Elliot Coupe, Nathan A. McMahon, Gavin K. Brennen, Thomas M. Stace

Abstract: We develop a most likely error Pauli error decoding algorithm for stabiliser codes based on general purpose integer optimisation. Using this decoder we analyse the performance of holographic codes against Pauli errors and find numerical evidence for thresholds against Pauli errors for bulk qubits. We compare the performance of holographic code families of various code rates and find phenomenologic… ▽ More We develop a most likely error Pauli error decoding algorithm for stabiliser codes based on general purpose integer optimisation. Using this decoder we analyse the performance of holographic codes against Pauli errors and find numerical evidence for thresholds against Pauli errors for bulk qubits. We compare the performance of holographic code families of various code rates and find phenomenological Pauli error thresholds ranging from $7\%$ to $16\%$, depending on the code rate. Additionally we give numerical evidence that specific distance measures of the codes we consider scales polynomially with number of physical qubits. △ Less

Submitted 17 December, 2020; v1 submitted 24 August, 2020; originally announced August 2020.

Comments: 7 pages, 6 figures

Journal ref: Phys. Rev. A 102, 062417 (2020)

Light Dark Matter Search with a High-Resolution Athermal Phonon Detector Operated Above Ground

Authors: I. Alkhatib, D. W. P. Amaral, T. Aralis, T. Aramaki, I. J. Arnquist, I. Ataee Langroudy, E. Azadbakht, S. Banik, D. Barker, C. Bathurst, D. A. Bauer, L. V. S. Bezerra, R. Bhattacharyya, T. Binder, M. A. Bowles, P. L. Brink, R. Bunker, B. Cabrera, R. Calkins, R. A. Cameron, C. Cartaro, D. G. Cerdeño, Y. -Y. Chang, M. Chaudhuri, R. Chen , et al. (99 additional authors not shown)

Abstract: We present limits on spin-independent dark matter-nucleon interactions using a $10.6$ $\mathrm{g}$ Si athermal phonon detector with a baseline energy resolution of $σ_E=3.86 \pm 0.04$ $(\mathrm{stat.})^{+0.19}_{-0.00}$ $(\mathrm{syst.})$ $\mathrm{eV}$. This exclusion analysis sets the most stringent dark matter-nucleon scattering cross-section limits achieved by a cryogenic detector for dark matte… ▽ More We present limits on spin-independent dark matter-nucleon interactions using a $10.6$ $\mathrm{g}$ Si athermal phonon detector with a baseline energy resolution of $σ_E=3.86 \pm 0.04$ $(\mathrm{stat.})^{+0.19}_{-0.00}$ $(\mathrm{syst.})$ $\mathrm{eV}$. This exclusion analysis sets the most stringent dark matter-nucleon scattering cross-section limits achieved by a cryogenic detector for dark matter particle masses from $93$ to $140$ $\mathrm{MeV}/c^2$, with a raw exposure of $9.9$ $\mathrm{g}\cdot\mathrm{d}$ acquired at an above-ground facility. This work illustrates the scientific potential of detectors with athermal phonon sensors with eV-scale energy resolution for future dark matter searches. △ Less

Submitted 12 October, 2021; v1 submitted 21 July, 2020; originally announced July 2020.

Comments: 7 pages, 4 figures, this version includes ancillary files from official data release

Journal ref: Phys. Rev. Lett. 127, 061801 (2021)

PrototypeML: A Neural Network Integrated Design and Development Environment

Authors: Daniel Reiss Harris

Abstract: Neural network architectures are most often conceptually designed and described in visual terms, but are implemented by writing error-prone code. PrototypeML is a machine learning development environment that bridges the dichotomy between the design and development processes: it provides a highly intuitive visual neural network design interface that supports (yet abstracts) the full capabilities o… ▽ More Neural network architectures are most often conceptually designed and described in visual terms, but are implemented by writing error-prone code. PrototypeML is a machine learning development environment that bridges the dichotomy between the design and development processes: it provides a highly intuitive visual neural network design interface that supports (yet abstracts) the full capabilities of the PyTorch deep learning framework, reduces model design and development time, makes debugging easier, and automates many framework and code writing idiosyncrasies. In this paper, we detail the deep learning development deficiencies that drove the implementation of PrototypeML, and propose a hybrid approach to resolve these issues without limiting network expressiveness or reducing code quality. We demonstrate the real-world benefits of a visual approach to neural network design for research, industry and teaching. Available at https://PrototypeML.com △ Less

Submitted 1 July, 2020; originally announced July 2020.

Comments: 10 pages, 6 figures. Submitted to NeurIPS 2020. More details available at https://PrototypeML.com

ACM Class: I.2.5; D.1.7

Array Programming with NumPy

Authors: Charles R. Harris, K. Jarrod Millman, Stéfan J. van der Walt, Ralf Gommers, Pauli Virtanen, David Cournapeau, Eric Wieser, Julian Taylor, Sebastian Berg, Nathaniel J. Smith, Robert Kern, Matti Picus, Stephan Hoyer, Marten H. van Kerkwijk, Matthew Brett, Allan Haldane, Jaime Fernández del Río, Mark Wiebe, Pearu Peterson, Pierre Gérard-Marchant, Kevin Sheppard, Tyler Reddy, Warren Weckesser, Hameer Abbasi, Christoph Gohlke , et al. (1 additional authors not shown)

Abstract: Array programming provides a powerful, compact, expressive syntax for accessing, manipulating, and operating on data in vectors, matrices, and higher-dimensional arrays. NumPy is the primary array programming library for the Python language. It plays an essential role in research analysis pipelines in fields as diverse as physics, chemistry, astronomy, geoscience, biology, psychology, material sci… ▽ More Array programming provides a powerful, compact, expressive syntax for accessing, manipulating, and operating on data in vectors, matrices, and higher-dimensional arrays. NumPy is the primary array programming library for the Python language. It plays an essential role in research analysis pipelines in fields as diverse as physics, chemistry, astronomy, geoscience, biology, psychology, material science, engineering, finance, and economics. For example, in astronomy, NumPy was an important part of the software stack used in the discovery of gravitational waves and the first imaging of a black hole. Here we show how a few fundamental array concepts lead to a simple and powerful programming paradigm for organizing, exploring, and analyzing scientific data. NumPy is the foundation upon which the entire scientific Python universe is constructed. It is so pervasive that several projects, targeting audiences with specialized needs, have developed their own NumPy-like interfaces and array objects. Because of its central position in the ecosystem, NumPy increasingly plays the role of an interoperability layer between these new array computation libraries. △ Less

Submitted 17 June, 2020; originally announced June 2020.

Journal ref: Nature 585, 357 (2020)

Constraints on low-mass, relic dark matter candidates from a surface-operated SuperCDMS single-charge sensitive detector

Authors: SuperCDMS Collaboration, D. W. Amaral, T. Aralis, T. Aramaki, I. J. Arnquist, E. Azadbakht, S. Banik, D. Barker, C. Bathurst, D. A. Bauer, L. V. S. Bezerra, R. Bhattacharyya, T. Binder, M. A. Bowles, P. L. Brink, R. Bunker, B. Cabrera, R. Calkins, R. A. Cameron, C. Cartaro, D. G. Cerdeño, Y. -Y. Chang, R. Chen, N. Chott, J. Cooley , et al. (94 additional authors not shown)

Abstract: This article presents an analysis and the resulting limits on light dark matter inelastically scattering off of electrons, and on dark photon and axion-like particle absorption, using a second-generation SuperCDMS high-voltage eV-resolution detector. The 0.93 gram Si detector achieved a 3 eV phonon energy resolution; for a detector bias of 100 V, this corresponds to a charge resolution of 3% of a… ▽ More This article presents an analysis and the resulting limits on light dark matter inelastically scattering off of electrons, and on dark photon and axion-like particle absorption, using a second-generation SuperCDMS high-voltage eV-resolution detector. The 0.93 gram Si detector achieved a 3 eV phonon energy resolution; for a detector bias of 100 V, this corresponds to a charge resolution of 3% of a single electron-hole pair. The energy spectrum is reported from a blind analysis with 1.2 gram-days of exposure acquired in an above-ground laboratory. With charge carrier trapping and impact ionization effects incorporated into the dark matter signal models, the dark matter-electron cross section $\barσ_{e}$ is constrained for dark matter masses from 0.5--$10^{4} $MeV$/c^{2}$; in the mass range from 1.2--50 eV$/c^{2}$ the dark photon kinetic mixing parameter $\varepsilon$ and the axioelectric coupling constant $g_{ae}$ are constrained. The minimum 90% confidence-level upper limits within the above mentioned mass ranges are $\barσ_{e}\,=\,8.7\times10^{-34}$ cm$^{2}$, $\varepsilon\,=\,3.3\times10^{-14}$, and $g_{ae}\,=\,1.0\times10^{-9}$. △ Less

Submitted 29 January, 2021; v1 submitted 28 May, 2020; originally announced May 2020.

Comments: 5 pages + title and references, 3 figures and 1 table

Journal ref: Phys. Rev. D 102, 091101 (2020)

The GWR route map: a guide to the informed application of Geographically Weighted Regression

Authors: Alexis Comber, Chris Brunsdon, Martin Charlton, Guanpeng Dong, Rich Harris, Binbin Lu, Yihe Lü, Daisuke Murakami, Tomoki Nakaya, Yunqiang Wang, Paul Harris

Abstract: Geographically Weighted Regression (GWR) is increasingly used in spatial analyses of social and environmental data. It allows spatial heterogeneities in processes and relationships to be investigated through a series of local regression models rather than a global one. Standard GWR assumes that the relationships between the response and predictor variables operate at the same spatial scale, which… ▽ More Geographically Weighted Regression (GWR) is increasingly used in spatial analyses of social and environmental data. It allows spatial heterogeneities in processes and relationships to be investigated through a series of local regression models rather than a global one. Standard GWR assumes that the relationships between the response and predictor variables operate at the same spatial scale, which is frequently not the case. To address this, several GWR variants have been proposed. This paper describes a route map to inform the choice of whether to use a GWR model or not, and if so which of three core variants to apply: a standard GWR, a mixed GWR or a multiscale GWR (MS-GWR). The route map comprises primary steps: a basic linear regression, a MS-GWR, and investigations of the results of these. The paper provides guidance for deciding whether to use a GWR approach, and if so for determining the appropriate GWR variant. It describes the importance of investigating a number of secondary issues at global and local scales including collinearity, the influence of outliers, and dependent error terms. Code and data for the case study used to illustrate the route map are provided, and further considerations are described in an extensive Appendix. △ Less

Submitted 14 April, 2020; v1 submitted 13 April, 2020; originally announced April 2020.

Comments: 34 pages, 8 figures

MSC Class: 91B72 (Primary); 62M30 (Secondary) ACM Class: G.3

Astrophotonics: astronomy and modern optics

Authors: S. Minardi, R. Harris, L. Labadie

Abstract: Much of the progress in Astronomy has been driven by instrumental developments, from the first telescopes to fiber fed spectrographs. In this review we describe the field of astrophotonics, a combination of photonics and astronomical instrumentation that has the potential to drive the next generation of developements. We begin with the science cases that have been identified as possibly benefiting… ▽ More Much of the progress in Astronomy has been driven by instrumental developments, from the first telescopes to fiber fed spectrographs. In this review we describe the field of astrophotonics, a combination of photonics and astronomical instrumentation that has the potential to drive the next generation of developements. We begin with the science cases that have been identified as possibly benefiting from astrophotonic devices. We then discuss devices, methods and developments in the field along with the advantages they provide. We conclude by describing possible future developments in the field and their influence on astronomy. △ Less

Submitted 8 November, 2020; v1 submitted 27 March, 2020; originally announced March 2020.

Comments: 99 pages, 7 figures - Invited review article for The Astronomy and Astrophysics Review. Submitted on 9 Nov 2018, Resubmission 8 Nov 2020

Giant leaps and long excursions: fluctuation mechanisms in systems with long-range memory

Authors: Robert L. Jack, Rosemary J. Harris

Abstract: We analyse large deviations of time-averaged quantities in stochastic processes with long-range memory, where the dynamics at time t depends itself on the value q_t of the time-averaged quantity. First we consider the elephant random walk and a Gaussian variant of this model, identifying two mechanisms for unusual fluctuation behaviour, which differ from the Markovian case. In particular, the memo… ▽ More We analyse large deviations of time-averaged quantities in stochastic processes with long-range memory, where the dynamics at time t depends itself on the value q_t of the time-averaged quantity. First we consider the elephant random walk and a Gaussian variant of this model, identifying two mechanisms for unusual fluctuation behaviour, which differ from the Markovian case. In particular, the memory can lead to large deviation principles with reduced speeds, and to non-analytic rate functions. We then explain how the mechanisms operating in these two models are generic for memory-dependent dynamics and show other examples including a non-Markovian symmetric exclusion process. △ Less

Submitted 16 June, 2020; v1 submitted 7 March, 2020; originally announced March 2020.

Comments: longer version (16 pages), with more detailed discussion

Journal ref: Phys. Rev. E 102, 012154 (2020)

Astronomical Applications of Multi-Core Fiber Technology

Authors: Nemanja Jovanovic, Robert J. Harris, Nick Cvetojevic

Abstract: Optical fibers have altered astronomical instrument design by allowing for a complex, often large instrument to be mounted in a remote and stable location with respect to the telescope. The fibers also enable the possibility to rearrange the signal from a focal plane to form a psuedo-slit at the entrance to a spectrograph, optimizing the detector usage and enabling the study of hundreds of thousan… ▽ More Optical fibers have altered astronomical instrument design by allowing for a complex, often large instrument to be mounted in a remote and stable location with respect to the telescope. The fibers also enable the possibility to rearrange the signal from a focal plane to form a psuedo-slit at the entrance to a spectrograph, optimizing the detector usage and enabling the study of hundreds of thousands of stars or galaxies simultaneously. Multi-core fibers in particular offer several favorable properties with respect to traditional fibers: 1) the separation between single-mode cores is greatly reduced and highly regular with respect to free standing fibers, 2) they offer a monolithic package with multi-fiber capabilities and 3) they operate at the diffraction limit. These properties have enabled the realization of single component photonic lanterns, highly simplified fiber Bragg gratings, and advanced fiber mode scramblers. In addition, the precise grid of cores has enabled the design of efficient single-mode fiber integral field units for spectroscopy. In this paper, we provide an overview of the broad range of applications enabled by multi-core fiber technology in astronomy and outline future areas of development. △ Less

Submitted 31 March, 2020; v1 submitted 4 February, 2020; originally announced February 2020.

Comments: 10 pages, 5 figures

Journal ref: IEEE Journal of Selected Topics in Quantum Electronics, Vol 26, No. 4, 2020