Showing 1–50 of 114 results for author: Hamilton, C

Observation of disorder-free localization and efficient disorder averaging on a quantum processor

Authors: Gaurav Gyawali, Tyler Cochran, Yuri Lensky, Eliott Rosenberg, Amir H. Karamlou, Kostyantyn Kechedzhi, Julia Berndtsson, Tom Westerhout, Abraham Asfaw, Dmitry Abanin, Rajeev Acharya, Laleh Aghababaie Beni, Trond I. Andersen, Markus Ansmann, Frank Arute, Kunal Arya, Nikita Astrakhantsev, Juan Atalaya, Ryan Babbush, Brian Ballard, Joseph C. Bardin, Andreas Bengtsson, Alexander Bilmes, Gina Bortoli, Alexandre Bourassa , et al. (195 additional authors not shown)

Abstract: One of the most challenging problems in the computational study of localization in quantum manybody systems is to capture the effects of rare events, which requires sampling over exponentially many disorder realizations. We implement an efficient procedure on a quantum processor, leveraging quantum parallelism, to efficiently sample over all disorder realizations. We observe localization without d… ▽ More One of the most challenging problems in the computational study of localization in quantum manybody systems is to capture the effects of rare events, which requires sampling over exponentially many disorder realizations. We implement an efficient procedure on a quantum processor, leveraging quantum parallelism, to efficiently sample over all disorder realizations. We observe localization without disorder in quantum many-body dynamics in one and two dimensions: perturbations do not diffuse even though both the generator of evolution and the initial states are fully translationally invariant. The disorder strength as well as its density can be readily tuned using the initial state. Furthermore, we demonstrate the versatility of our platform by measuring Renyi entropies. Our method could also be extended to higher moments of the physical observables and disorder learning. △ Less

Submitted 9 October, 2024; originally announced October 2024.

Quantum error correction below the surface code threshold

Authors: Rajeev Acharya, Laleh Aghababaie-Beni, Igor Aleiner, Trond I. Andersen, Markus Ansmann, Frank Arute, Kunal Arya, Abraham Asfaw, Nikita Astrakhantsev, Juan Atalaya, Ryan Babbush, Dave Bacon, Brian Ballard, Joseph C. Bardin, Johannes Bausch, Andreas Bengtsson, Alexander Bilmes, Sam Blackwell, Sergio Boixo, Gina Bortoli, Alexandre Bourassa, Jenna Bovaird, Leon Brill, Michael Broughton, David A. Browne , et al. (224 additional authors not shown)

Abstract: Quantum error correction provides a path to reach practical quantum computing by combining multiple physical qubits into a logical qubit, where the logical error rate is suppressed exponentially as more qubits are added. However, this exponential suppression only occurs if the physical error rate is below a critical threshold. In this work, we present two surface code memories operating below this… ▽ More Quantum error correction provides a path to reach practical quantum computing by combining multiple physical qubits into a logical qubit, where the logical error rate is suppressed exponentially as more qubits are added. However, this exponential suppression only occurs if the physical error rate is below a critical threshold. In this work, we present two surface code memories operating below this threshold: a distance-7 code and a distance-5 code integrated with a real-time decoder. The logical error rate of our larger quantum memory is suppressed by a factor of $Λ$ = 2.14 $\pm$ 0.02 when increasing the code distance by two, culminating in a 101-qubit distance-7 code with 0.143% $\pm$ 0.003% error per cycle of error correction. This logical memory is also beyond break-even, exceeding its best physical qubit's lifetime by a factor of 2.4 $\pm$ 0.3. We maintain below-threshold performance when decoding in real time, achieving an average decoder latency of 63 $μ$s at distance-5 up to a million cycles, with a cycle time of 1.1 $μ$s. To probe the limits of our error-correction performance, we run repetition codes up to distance-29 and find that logical performance is limited by rare correlated error events occurring approximately once every hour, or 3 $\times$ 10$^9$ cycles. Our results present device performance that, if scaled, could realize the operational requirements of large scale fault-tolerant quantum algorithms. △ Less

Submitted 24 August, 2024; originally announced August 2024.

Comments: 10 pages, 4 figures, Supplementary Information

Comparing NASA Discovery and New Frontiers Class Mission Concepts for the Io Volcano Observer (IVO)

Authors: Christopher W. Hamilton, Alfred S. McEwen, Laszlo Keszthelyi, Lynn M. Carter, Ashley G. Davies, Katherine de Kleer, Kandis Lea Jessup, Xianzhe Jia, James T. Keane, Kathleen Mandt, Francis Nimmo, Chris Paranicas, Ryan S. Park, Jason E. Perry, Anne Pommier, Jani Radebaugh, Sarah S. Sutton, Audrey Vorburger, Peter Wurz, Cauê Borlina, Amanda F. Haapala, Daniella N. DellaGiustina, Brett W. Denevi, Sarah M. Hörst, Sascha Kempf , et al. (9 additional authors not shown)

Abstract: Jupiter's moon Io is a highly compelling target for future exploration that offers critical insight into tidal dissipation processes and the geology of high heat flux worlds, including primitive planetary bodies, such as the early Earth, that are shaped by enhanced rates of volcanism. Io is also important for understanding the development of volcanogenic atmospheres and mass-exchange within the Ju… ▽ More Jupiter's moon Io is a highly compelling target for future exploration that offers critical insight into tidal dissipation processes and the geology of high heat flux worlds, including primitive planetary bodies, such as the early Earth, that are shaped by enhanced rates of volcanism. Io is also important for understanding the development of volcanogenic atmospheres and mass-exchange within the Jupiter System. However, fundamental questions remain about the state of Io's interior, surface, and atmosphere, as well as its role in the evolution of the Galilean satellites. The Io Volcano Observer (IVO) would address these questions by achieving the following three key goals: (A) Determine how and where tidal heat is generated inside Io; (B) Understand how tidal heat is transported to the surface of Io; and (C) Understand how Io is evolving. IVO was selected for Phase A study through the NASA Discovery program in 2020 and, in anticipation of a New Frontiers 5 opportunity, an enhanced IVO-NF mission concept was advanced that would increase the Baseline mission from 10 flybys to 20, with an improved radiation design; employ a Ka-band communications to double IVO's total data downlink; add a wide angle camera for color and stereo mapping; add a dust mass spectrometer; and lower the altitude of later flybys to enable new science. This study compares and contrasts the mission architecture, instrument suite, and science objectives for Discovery (IVO) and New Frontiers (IVO-NF) missions to Io, and advocates for continued prioritization of Io as an exploration target for New Frontiers. △ Less

Submitted 14 August, 2024; originally announced August 2024.

Comments: Submitted to The Planetary Science Journal for peer-review on 14 August 2024

Hybrid Magnonics with Localized Spoof Surface Plasmon Polaritons

Authors: Yuzan Xiong, Andrew Christy, Zixin Yan, Amin Pishehvar, Muntasir Mahdi, Junming Wu, James F. Cahoon, Binbin Yang, Michael C. Hamilton, Xufeng Zhang, Wei Zhang

Abstract: Hybrid magnonic systems have emerged as a promising direction for information propagation with preserved coherence. Due to high tunability of magnons, their interactions with microwave photons can be engineered to probe novel phenomena based on strong photon-magnon coupling. Improving the photon-magnon coupling strength can be done by tuning the structure of microwave resonators to better interact… ▽ More Hybrid magnonic systems have emerged as a promising direction for information propagation with preserved coherence. Due to high tunability of magnons, their interactions with microwave photons can be engineered to probe novel phenomena based on strong photon-magnon coupling. Improving the photon-magnon coupling strength can be done by tuning the structure of microwave resonators to better interact with the magnon counterpart. Planar resonators have been explored due to their potential for on-chip integration, but only common modes from stripline-based resonators have been used. Here, we present a microwave spiral resonator supporting the spoof localized surface plasmons (LSPs) and implement it to the investigation of photon-magnon coupling for hybrid magnonic applications. We showcase strong magnon-LSP photon coupling using a ferrimagnetic yttrium iron garnet sphere. We discuss the dependence of the spiral resonator design to the engineering capacity of the photon mode frequency and spatial field distributions, via both experiment and simulation. By the localized photon mode profiles, the resulting magnetic field concentrates near the surface dielectrics, giving rise to an enhanced magnetic filling factor. The strong coupling and large engineering space render the spoof LSPs an interesting contender in developing novel hybrid magnonic systems and functionalities. △ Less

Submitted 2 September, 2024; v1 submitted 13 August, 2024; originally announced August 2024.

Comments: 13 pages, 13 figures

Journal ref: Phys. Rev. Applied (2024)

Galactokinetics

Authors: Chris Hamilton, Shaunak Modak, Scott Tremaine

Abstract: Galactic disks lie at the heart of many of the most pressing astrophysical puzzles. There are sophisticated kinetic theories that describe some aspects of galaxy disk dynamics, but extracting quantitative predictions from those theories has proven very difficult, meaning they have shed little light on observations/simulations of galaxies. Here, we begin to address this issue by developing a tracta… ▽ More Galactic disks lie at the heart of many of the most pressing astrophysical puzzles. There are sophisticated kinetic theories that describe some aspects of galaxy disk dynamics, but extracting quantitative predictions from those theories has proven very difficult, meaning they have shed little light on observations/simulations of galaxies. Here, we begin to address this issue by developing a tractable theory describing fluctuations and transport in thin galactic disks. Our main conceptual advance is to split potential fluctuations into asymptotic wavelength regimes relative to orbital guiding radius and epicyclic amplitude (similar to plasma gyrokinetics), and then to treat separately the dynamics in each regime. As an illustration, we apply our results to quasilinear theory, calculating the angular-momentum transport due to a transient spiral. At each stage we verify our formulae with numerical examples. Our approach should simplify many important calculations in galactic disk dynamics. △ Less

Submitted 6 August, 2024; originally announced August 2024.

Comments: 18 pages, 8 figures. Submitted to ApJ

A Journey in Implementing Computational Physics from the Ground Up

Authors: Maria C. Babiuc Hamilton

Abstract: This chapter narrates the journey of developing and integrating computing into the physics curriculum through three consecutive courses, each tailored to the learners' level. It starts with the entry-level "Physics Playground in Python" for high school and freshman students with no programming experience, designed in the spirit of the "Hello World" approach. At the sophomore and junior level, stud… ▽ More This chapter narrates the journey of developing and integrating computing into the physics curriculum through three consecutive courses, each tailored to the learners' level. It starts with the entry-level "Physics Playground in Python" for high school and freshman students with no programming experience, designed in the spirit of the "Hello World" approach. At the sophomore and junior level, students from all sciences and engineering disciplines learn "Scientific Computing with Python" in an environment based on the "Two Bites at Every Apple" approach. Ultimately, upper undergraduate and entry-level graduate students take "Computational Physics," to develop their skills in solving advanced problems using complex numerical algorithms and computational tools. This journey showcases the increasing complexity and sophistication of computational tools and techniques that can be incorporated into the physical science curriculum, serving as a guide for educators looking to integrate computing into their teaching. It also aims to inspire students by showcasing the impact and potential of computational methods in science education and research. △ Less

Submitted 3 August, 2024; originally announced August 2024.

Comments: contributed chapter for the Computing in Physics Education book

Thermalization and Criticality on an Analog-Digital Quantum Simulator

Authors: Trond I. Andersen, Nikita Astrakhantsev, Amir H. Karamlou, Julia Berndtsson, Johannes Motruk, Aaron Szasz, Jonathan A. Gross, Alexander Schuckert, Tom Westerhout, Yaxing Zhang, Ebrahim Forati, Dario Rossi, Bryce Kobrin, Agustin Di Paolo, Andrey R. Klots, Ilya Drozdov, Vladislav D. Kurilovich, Andre Petukhov, Lev B. Ioffe, Andreas Elben, Aniket Rath, Vittorio Vitale, Benoit Vermersch, Rajeev Acharya, Laleh Aghababaie Beni , et al. (202 additional authors not shown)

Abstract: Understanding how interacting particles approach thermal equilibrium is a major challenge of quantum simulators. Unlocking the full potential of such systems toward this goal requires flexible initial state preparation, precise time evolution, and extensive probes for final state characterization. We present a quantum simulator comprising 69 superconducting qubits which supports both universal qua… ▽ More Understanding how interacting particles approach thermal equilibrium is a major challenge of quantum simulators. Unlocking the full potential of such systems toward this goal requires flexible initial state preparation, precise time evolution, and extensive probes for final state characterization. We present a quantum simulator comprising 69 superconducting qubits which supports both universal quantum gates and high-fidelity analog evolution, with performance beyond the reach of classical simulation in cross-entropy benchmarking experiments. Emulating a two-dimensional (2D) XY quantum magnet, we leverage a wide range of measurement techniques to study quantum states after ramps from an antiferromagnetic initial state. We observe signatures of the classical Kosterlitz-Thouless phase transition, as well as strong deviations from Kibble-Zurek scaling predictions attributed to the interplay between quantum and classical coarsening of the correlated domains. This interpretation is corroborated by injecting variable energy density into the initial state, which enables studying the effects of the eigenstate thermalization hypothesis (ETH) in targeted parts of the eigenspectrum. Finally, we digitally prepare the system in pairwise-entangled dimer states and image the transport of energy and vorticity during thermalization. These results establish the efficacy of superconducting analog-digital quantum processors for preparing states across many-body spectra and unveiling their thermalization dynamics. △ Less

Submitted 8 July, 2024; v1 submitted 27 May, 2024; originally announced May 2024.

Boson Sampling from Non-Gaussian States

Authors: Craig S. Hamilton, Igor Jex

Abstract: Boson sampling has emerged as an important tool to demonstrate the difference between quantum and classical computers and has attracted the interest of experimentalists and theoreticians. In this work we study Boson sampling from general, single-mode states using a scheme that can generate any such state by combining Gaussian states and photon number measurements. We derive a formula that can be u… ▽ More Boson sampling has emerged as an important tool to demonstrate the difference between quantum and classical computers and has attracted the interest of experimentalists and theoreticians. In this work we study Boson sampling from general, single-mode states using a scheme that can generate any such state by combining Gaussian states and photon number measurements. We derive a formula that can be used to calculate the output photon number probabilities of these states after they travel through a linear interferometer. This extends the Boson sampling protocol to the widest array of possible single-mode states and from this we show that the complexity scaling of all such states is similar and hence there is no complexity advantage of using complex input states over simpler ones. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Kinetic Theory of Stellar Systems: A Tutorial

Authors: Chris Hamilton, Jean-Baptiste Fouvry

Abstract: Stellar systems - star clusters, galaxies, dark matter haloes, and so on - are ubiquitous characters in the evolutionary tale of our Universe. This tutorial article is an introduction to the collective dynamical evolution of the very large numbers of stars and/or other self-gravitating objects that comprise such systems, i.e. their kinetic theory. We begin by introducing the basic phenomenology… ▽ More Stellar systems - star clusters, galaxies, dark matter haloes, and so on - are ubiquitous characters in the evolutionary tale of our Universe. This tutorial article is an introduction to the collective dynamical evolution of the very large numbers of stars and/or other self-gravitating objects that comprise such systems, i.e. their kinetic theory. We begin by introducing the basic phenomenology of stellar systems, and explaining why and when we must develop a kinetic theory that transcends the traditional two-body relaxation picture of Chandrasekhar. We study the orbits that comprise stellar systems, how those orbits are modified by perturbations, how a system responds self-consistently to fluctuations in its gravitational potential, and how one can predict the long term fate of a stellar system in various dynamical regimes. Though our treatment is necessarily mathematical, we develop the formalism only to the extent that it facilitates real calculations. We give many examples throughout the text of the equations being applied to topics of major astrophysical importance. Furthermore, in the 1960s and 1970s the kinetic theory of stellar systems was a fledgling subject which developed in tandem with the kinetic theory of plasmas. However, the two fields have long since diverged. Yet once one has become fluent in both Plasmaish and Galacticese, and has a dictionary relating the two, one can pull ideas directly from one field to solve a problem in the other. Therefore, another aim of this tutorial article is to provide our plasma colleagues with a jargon-light understanding of the key properties of stellar systems, to point out the many direct analogies between stellar- and plasma-kinetic calculations, and ultimately to convince them that stellar dynamics and plasma kinetics are, in a deep and beautiful and useful sense, the same thing. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Comments: 66 Pages, 34 Figures. Invited tutorial article for Physics of Plasmas. Comments/criticism strongly encouraged, especially from non-specialists!

How Do Nuclear Isomers Influence the Gamma-Ray Bursts in Binary Neutron Star Mergers?

Authors: Maria C. Babiuc Hamilton, Joseph I. Powell

Abstract: Neutron star mergers are astrophysical `gold mines,' synthesizing over half of the elements heavier than iron through rapid neutron capture nucleosynthesis. The observation of the binary neutron star merger GW170817, detected both in gravitational waves and electromagnetic radiation, marked a breakthrough. One electromagnetic component of this event, the gamma ray burst GRB 170817A, has an unresol… ▽ More Neutron star mergers are astrophysical `gold mines,' synthesizing over half of the elements heavier than iron through rapid neutron capture nucleosynthesis. The observation of the binary neutron star merger GW170817, detected both in gravitational waves and electromagnetic radiation, marked a breakthrough. One electromagnetic component of this event, the gamma ray burst GRB 170817A, has an unresolved aspect: the characteristics of its prompt gamma-ray emission spectrum. In this work, we investigate how gamma-ray spectra in such GRBs may be influenced by de-excitations from isomeric transitions. Our study begins with a review of current knowledge on GRB structure and of r-process nucleosynthesis in neutron star collisions, focusing on the role of nuclear isomers in these settings. We then test our hypothesis by developing criteria to select representative isomers, based on known solar element abundances, for modeling GRB spectral characteristics. We integrate these criteria into an interactive web page, facilitating the construction and analysis of relevant gamma-ray spectra from isomeric transitions. Our analysis reveals that three isomers (zirconium, lead and yttrium) stand out for their potential to impact the prompt GRB spectrum due to their specific properties. This information allows us to incorporate nuclear isomer data into astrophysical simulations and calculate isomeric abundances generated by astrophysical r-processes in neutron star mergers and their imprint on the detected signal. △ Less

Submitted 3 August, 2024; v1 submitted 9 February, 2024; originally announced February 2024.

Comments: 24 pages, 9 figures

Journal ref: Front. Astron. Space Sci., Sec. Cosmology, Volume 11 (2024)

Conditional constrained and unconstrained quantization for a uniform distribution on a hexagon

Authors: Christina Hamilton, Evans Nyanney, Megha Pandey, Mrinal K. Roychowdhury

Abstract: In this paper, we have considered a uniform distribution on a regular hexagon and the set of all its six vertices as a conditional set. For the uniform distribution under the conditional set first, for all positive integers $n\geq 6$, we obtain the conditional optimal sets of $n$-points and the $n$th conditional quantization errors, and then we calculate the conditional quantization dimension and… ▽ More In this paper, we have considered a uniform distribution on a regular hexagon and the set of all its six vertices as a conditional set. For the uniform distribution under the conditional set first, for all positive integers $n\geq 6$, we obtain the conditional optimal sets of $n$-points and the $n$th conditional quantization errors, and then we calculate the conditional quantization dimension and the conditional quantization coefficient in the unconstrained scenario. Then, for the uniform distribution on the hexagon taking the same conditional set, we investigate the conditional constrained optimal sets of n-points and the conditional constrained quantization errors for all $n \geq 6$, taking the constraint as the circumcircle, incircle, and then the different diagonals of the hexagon. △ Less

Submitted 23 January, 2024; v1 submitted 19 January, 2024; originally announced January 2024.

MSC Class: 60Exx; 94A34

Binary mergers in the centers of galaxies: synergy between stellar flybys and tidal fields

Authors: Mila Winter-Granic, Cristobal Petrovich, Valentín Peña-Donaire, Chris Hamilton

Abstract: Galactic centers are very dynamically active environments, often harbouring a nuclear star cluster and supermassive black hole at their cores. Binaries in these environments are subject to strong tidal fields that can efficiently torque its orbit, exciting near unity eccentricities that ultimately lead to their merger. In turn, frequent close interactions with passing stars impulsively perturb the… ▽ More Galactic centers are very dynamically active environments, often harbouring a nuclear star cluster and supermassive black hole at their cores. Binaries in these environments are subject to strong tidal fields that can efficiently torque its orbit, exciting near unity eccentricities that ultimately lead to their merger. In turn, frequent close interactions with passing stars impulsively perturb the orbit of the binary, generally softening their orbit until their evaporation, potentially hindering the role of tides to drive these mergers. In this work, we study the evolution of compact object binaries in the galactic center and their merger rates, focusing for the first time on the combined effect of the cluster's tidal field and flyby interactions. We find a significant synergy between both processes, where merger rates increase by a factor of ~10-30 compared to models in which only flybys or tides are taken into account. This synergy is a consequence of the persistent tides-driven eccentricity excitation that is enhanced by the gradual diffusion of $j_z$ driven by flybys. The merger efficiency peaks when the diffusion rate is ~10-100 slower than the tides-driven torquing. Added to this synergy, we also find that the gradual softening of the binary can lift the relativistic quenching of initially tight binaries, otherwise unable to reach extreme eccentricities, and thus expanding the available phase-space for mergers. Cumulatively, we conclude that despite the gradual softening of binaries due to flybys, these greatly enhance their merger rates in galaxy centers by promoting the tidal field-driven eccentricity excitation. △ Less

Submitted 25 June, 2024; v1 submitted 28 December, 2023; originally announced December 2023.

Comments: 15 pages, 11 figures. Revised version, submitted to the AAS Journals

Eccentricity dynamics of wide binaries -- II. The effect of stellar encounters and constraints on formation channels

Authors: Chris Hamilton, Shaunak Modak

Abstract: GAIA wide stellar binaries (separations $\sim 10^3-10^{4.5}$ AU) are observed to have a superthermal eccentricity distribution function (DF), well-fit by $P(e) \propto e^α$ with $α\sim 1.2$. In Modak \& Hamilton (2023), we proved that this DF cannot have been produced by Galactic tidal torques starting from any realistic DF that was not already superthermal. Here, we consider the other major dynam… ▽ More GAIA wide stellar binaries (separations $\sim 10^3-10^{4.5}$ AU) are observed to have a superthermal eccentricity distribution function (DF), well-fit by $P(e) \propto e^α$ with $α\sim 1.2$. In Modak \& Hamilton (2023), we proved that this DF cannot have been produced by Galactic tidal torques starting from any realistic DF that was not already superthermal. Here, we consider the other major dynamical effect on wide binaries: encounters with passing stars. We derive and solve the Fokker-Planck equation governing the evolution of binaries in semimajor axis and eccentricity under many weak, impulsive, penetrative stellar encounters. We show analytically that these encounters drive the eccentricity DF towards thermal on the same timescale as they drive the semimajor axes $a$ towards disruption, $t_\mathrm{ion} \sim 4\,\mathrm{Gyr}\,(a/10^4\,\mathrm{AU})^{-1}$. We conclude that the observed superthermal DF must derive from an even more superthermal (i.e. higher $α$) birth distribution. This requirement places strong constraints on the dominant binary formation channels. A testable prediction of our theory is that $α$ should be a monotonically decreasing function of binary age. △ Less

Submitted 2 July, 2024; v1 submitted 7 November, 2023; originally announced November 2023.

Comments: Accepted version, to be published in MNRAS

Rising Tides: Analytic Modeling of Tidal Effects in Binary Neutron Star Mergers

Authors: Alexander O'Dell, Maria C. Babiuc Hamilton

Abstract: The gravitational waves produced by binary neutron star mergers offer a unique window into matter behavior under extreme conditions. In this context, we model analytically the effect of matter on the gravitational waves from binary neutron star mergers. We start with a binary black hole system, leveraging the post-Newtonian formalism for the inspiral and the Backwards-one-Body model for the merger… ▽ More The gravitational waves produced by binary neutron star mergers offer a unique window into matter behavior under extreme conditions. In this context, we model analytically the effect of matter on the gravitational waves from binary neutron star mergers. We start with a binary black hole system, leveraging the post-Newtonian formalism for the inspiral and the Backwards-one-Body model for the merger. We combine the two methods to generate a baseline waveform and we validate our results against numerical relativity simulations. Next, we integrate tidal effects in phase and amplitude to account for matter and spacetime interaction, by using the NRTidal model, and test its accuracy against numerical relativity predictions, for two equations of state, finding a mismatch around the merger. Subsequently, we lift the restriction on the coefficients to be independent of the tidal deformability, and recalibrate them using the numerical relativity predictions. We obtain better fits for phase and amplitude around the merger, and are able to extend the phase modeling beyond the merger. We implement our method in a new open-source Python code, steered by a Jupyter Notebook. Our research offers new perspectives on analytically modeling the effect of tides on the gravitational waves from binary neutron star mergers. △ Less

Submitted 22 February, 2024; v1 submitted 29 July, 2023; originally announced July 2023.

Comments: 28 pages, 8 figures

Journal ref: Astronomy 2024, 3(3), 167-188

Relativistic phase space diffusion of compact object binaries in stellar clusters and hierarchical triples

Authors: Chris Hamilton, Roman R. Rafikov

Abstract: The LIGO/Virgo detections of compact object mergers have posed a challenge for theories of binary evolution and coalescence. One promising avenue for producing mergers dynamically is through secular eccentricity oscillations driven by an external perturber, be it a tertiary companion (as in the Lidov-Kozai (LK) mechanism) or the tidal field of the stellar cluster in which the binary orbits. The si… ▽ More The LIGO/Virgo detections of compact object mergers have posed a challenge for theories of binary evolution and coalescence. One promising avenue for producing mergers dynamically is through secular eccentricity oscillations driven by an external perturber, be it a tertiary companion (as in the Lidov-Kozai (LK) mechanism) or the tidal field of the stellar cluster in which the binary orbits. The simplest theoretical models of these oscillations use a 'doubly-averaged' (DA) approximation, averaging both over the binary's internal Keplerian orbit and its 'outer' barycentric orbit relative to the perturber. However, DA theories do not account for fluctuations of the perturbing torque on the outer orbital timescale, which are known to increase a binary's eccentricity beyond the maximum DA value, potentially accelerating mergers. Here we reconsider the impact of these short-timescale fluctuations in the test-particle quadrupolar limit for binaries perturbed by arbitrary spherical cluster potentials (including LK as a special case), {in particular including 1pN} general relativistic (GR) apsidal precession of the internal orbit. Focusing on the behavior of the binary orbital elements around peak eccentricity, we discover a new effect, relativistic phase space diffusion (RPSD), in which a binary can jump to a completely new dynamical trajectory on an outer orbital timescale, violating the approximate conservation of DA integrals of motion. RPSD arises from an interplay between secular behavior at extremely high eccentricity, short-timescale fluctuations, and rapid GR precession, and can change the subsequent secular evolution dramatically. This effect occurs even in hierarchical triples, but has not been uncovered until now. △ Less

Submitted 4 December, 2023; v1 submitted 6 June, 2023; originally announced June 2023.

Comments: Revised version, accepted for publication in ApJ

Stable Quantum-Correlated Many Body States through Engineered Dissipation

Authors: X. Mi, A. A. Michailidis, S. Shabani, K. C. Miao, P. V. Klimov, J. Lloyd, E. Rosenberg, R. Acharya, I. Aleiner, T. I. Andersen, M. Ansmann, F. Arute, K. Arya, A. Asfaw, J. Atalaya, J. C. Bardin, A. Bengtsson, G. Bortoli, A. Bourassa, J. Bovaird, L. Brill, M. Broughton, B. B. Buckley, D. A. Buell, T. Burger , et al. (142 additional authors not shown)

Abstract: Engineered dissipative reservoirs have the potential to steer many-body quantum systems toward correlated steady states useful for quantum simulation of high-temperature superconductivity or quantum magnetism. Using up to 49 superconducting qubits, we prepared low-energy states of the transverse-field Ising model through coupling to dissipative auxiliary qubits. In one dimension, we observed long-… ▽ More Engineered dissipative reservoirs have the potential to steer many-body quantum systems toward correlated steady states useful for quantum simulation of high-temperature superconductivity or quantum magnetism. Using up to 49 superconducting qubits, we prepared low-energy states of the transverse-field Ising model through coupling to dissipative auxiliary qubits. In one dimension, we observed long-range quantum correlations and a ground-state fidelity of 0.86 for 18 qubits at the critical point. In two dimensions, we found mutual information that extends beyond nearest neighbors. Lastly, by coupling the system to auxiliaries emulating reservoirs with different chemical potentials, we explored transport in the quantum Heisenberg model. Our results establish engineered dissipation as a scalable alternative to unitary evolution for preparing entangled many-body states on noisy quantum processors. △ Less

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

Journal ref: Science 383, 1332-1337 (2024)

Phase transition in Random Circuit Sampling

Authors: A. Morvan, B. Villalonga, X. Mi, S. Mandrà, A. Bengtsson, P. V. Klimov, Z. Chen, S. Hong, C. Erickson, I. K. Drozdov, J. Chau, G. Laun, R. Movassagh, A. Asfaw, L. T. A. N. Brandão, R. Peralta, D. Abanin, R. Acharya, R. Allen, T. I. Andersen, K. Anderson, M. Ansmann, F. Arute, K. Arya, J. Atalaya , et al. (160 additional authors not shown)

Abstract: Undesired coupling to the surrounding environment destroys long-range correlations on quantum processors and hinders the coherent evolution in the nominally available computational space. This incoherent noise is an outstanding challenge to fully leverage the computation power of near-term quantum processors. It has been shown that benchmarking Random Circuit Sampling (RCS) with Cross-Entropy Benc… ▽ More Undesired coupling to the surrounding environment destroys long-range correlations on quantum processors and hinders the coherent evolution in the nominally available computational space. This incoherent noise is an outstanding challenge to fully leverage the computation power of near-term quantum processors. It has been shown that benchmarking Random Circuit Sampling (RCS) with Cross-Entropy Benchmarking (XEB) can provide a reliable estimate of the effective size of the Hilbert space coherently available. The extent to which the presence of noise can trivialize the outputs of a given quantum algorithm, i.e. making it spoofable by a classical computation, is an unanswered question. Here, by implementing an RCS algorithm we demonstrate experimentally that there are two phase transitions observable with XEB, which we explain theoretically with a statistical model. The first is a dynamical transition as a function of the number of cycles and is the continuation of the anti-concentration point in the noiseless case. The second is a quantum phase transition controlled by the error per cycle; to identify it analytically and experimentally, we create a weak link model which allows varying the strength of noise versus coherent evolution. Furthermore, by presenting an RCS experiment with 67 qubits at 32 cycles, we demonstrate that the computational cost of our experiment is beyond the capabilities of existing classical supercomputers, even when accounting for the inevitable presence of noise. Our experimental and theoretical work establishes the existence of transitions to a stable computationally complex phase that is reachable with current quantum processors. △ Less

Submitted 21 December, 2023; v1 submitted 21 April, 2023; originally announced April 2023.

The linear response of stellar systems does not diverge at marginal stability

Authors: Chris Hamilton, Tobias Heinemann

Abstract: The linear response of a stellar system's gravitational potential to a perturbing mass comprises two distinct contributions. Most famously, the system will respond by forming a polarization 'wake' around the perturber. At the same time, the perturber may also excite one or more 'Landau modes', i.e. coherent oscillations of the entire stellar system which are either stable or unstable depending on… ▽ More The linear response of a stellar system's gravitational potential to a perturbing mass comprises two distinct contributions. Most famously, the system will respond by forming a polarization 'wake' around the perturber. At the same time, the perturber may also excite one or more 'Landau modes', i.e. coherent oscillations of the entire stellar system which are either stable or unstable depending on the system parameters. The amplitude of the first (wake) contribution is known to diverge as a system approaches marginal stability. In this paper we consider the linear response of a homogeneous stellar system to a point mass moving on a straight line orbit. We prove analytically that the divergence of the wake response is in fact cancelled by a corresponding divergence in the Landau mode response, rendering the total response finite. We demonstrate this cancellation explicitly for a box of stars with Maxwellian velocity distribution. Our results imply that polarization wakes may be much less efficient drivers of secular evolution than previously thought. More generally, any prior calculation that accounted for wakes but ignored modes - such as those based on the Balescu-Lenard equation - may need to be revised. △ Less

Submitted 10 August, 2023; v1 submitted 14 April, 2023; originally announced April 2023.

Comments: Final version, to be published in MNRAS. 5 pages, 2 figures

The James Webb Space Telescope Mission

Authors: Jonathan P. Gardner, John C. Mather, Randy Abbott, James S. Abell, Mark Abernathy, Faith E. Abney, John G. Abraham, Roberto Abraham, Yasin M. Abul-Huda, Scott Acton, Cynthia K. Adams, Evan Adams, David S. Adler, Maarten Adriaensen, Jonathan Albert Aguilar, Mansoor Ahmed, Nasif S. Ahmed, Tanjira Ahmed, Rüdeger Albat, Loïc Albert, Stacey Alberts, David Aldridge, Mary Marsha Allen, Shaune S. Allen, Martin Altenburg , et al. (983 additional authors not shown)

Abstract: Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono… ▽ More Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

Wide binary stars formed in the turbulent interstellar medium

Authors: Siyao Xu, Hsiang-Chih Hwang, Chris Hamilton, Dong Lai

Abstract: The ubiquitous interstellar turbulence regulates star formation and the scaling relations between the initial velocity differences and the initial separations of stars. We propose that the formation of wide binaries with initial separations $r$ in the range $\sim 10^3~\text{AU} \lesssim r \lesssim 10^5$ AU is a natural consequence of star formation in the turbulent interstellar medium. With the de… ▽ More The ubiquitous interstellar turbulence regulates star formation and the scaling relations between the initial velocity differences and the initial separations of stars. We propose that the formation of wide binaries with initial separations $r$ in the range $\sim 10^3~\text{AU} \lesssim r \lesssim 10^5$ AU is a natural consequence of star formation in the turbulent interstellar medium. With the decrease of $r$, the mean turbulent relative velocity $v_\text{tur}$ between a pair of stars decreases, while the largest velocity $v_\text{bon}$ at which they still may be gravitationally bound increases. When $v_\text{tur} < v_\text{bon}$, a wide binary can form. In this formation scenario, we derive the eccentricity distribution $p(e)$ of wide binaries for an arbitrary relative velocity distribution. By adopting a turbulent velocity distribution, we find that wide binaries at a given initial separation generally exhibit a superthermal $p(e)$. This provides a natural explanation for the observed superthermal $p(e)$ of the wide binaries in the Solar neighborhood. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Comments: 7 pages, 4 figures, submitted to ApJL

Eccentricity dynamics of wide binaries -- I. The effect of Galactic tides

Authors: Shaunak Modak, Chris Hamilton

Abstract: A major puzzle concerning the wide stellar binaries (semimajor axes $a\gtrsim 10^3$\,AU) in the Solar neighborhood is the origin of their observed superthermal eccentricity distribution function (DF), which is well-approximated by $P(e)\propto e^α$ with $α\approx 1.3$. This DF evolves under the combined influence of (i) tidal torques from the Galactic disk and (ii) scattering by passing stars, mol… ▽ More A major puzzle concerning the wide stellar binaries (semimajor axes $a\gtrsim 10^3$\,AU) in the Solar neighborhood is the origin of their observed superthermal eccentricity distribution function (DF), which is well-approximated by $P(e)\propto e^α$ with $α\approx 1.3$. This DF evolves under the combined influence of (i) tidal torques from the Galactic disk and (ii) scattering by passing stars, molecular clouds, and substructure. Recently, it was demonstrated that Galactic tides alone cannot produce a superthermal eccentricity DF from an initially isotropic, non-superthermal one, under the restrictive assumptions that the eccentricity DF was initially of power law form and then was rapidly phase-mixed toward a steady state by the tidal perturbation. In this paper we first prove analytically that this conclusion is valid at all times, regardless of these assumptions. We then adopt a thin Galactic disk model and numerically integrate the equations of motion for several ensembles of tidally perturbed wide binaries to study the time evolution in detail. We find that even non-power law DFs can be described by an effective power law index $α_\mathrm{eff}$ which accurately characterizes both their initial and final states, and that a DF with initial (effective or exact) power law index $α_\mathrm{i}$ is transformed by Galactic tides into another power law with index $α_\mathrm{f}\approx (1+α_\mathrm{i})/2$ on a timescale $ \sim 4\,\mathrm{Gyr}\,(a/10^4\mathrm{AU})^{-3/2}$. In a companion paper, we investigate separately the effect of stellar scattering. As the GAIA data continues to improve, these results will place strong constraints on wide binary formation channels. △ Less

Submitted 27 July, 2023; v1 submitted 27 March, 2023; originally announced March 2023.

Comments: Updated to accepted version: 14 pages, 9 figures. Main result in Figure 7. Comments still welcome!

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 524, Issue 2, September 2023, Pages 3102-3115

Measurement-induced entanglement and teleportation on a noisy quantum processor

Authors: Jesse C. Hoke, Matteo Ippoliti, Eliott Rosenberg, Dmitry Abanin, Rajeev Acharya, Trond I. Andersen, Markus Ansmann, Frank Arute, Kunal Arya, Abraham Asfaw, Juan Atalaya, Joseph C. Bardin, Andreas Bengtsson, Gina Bortoli, Alexandre Bourassa, Jenna Bovaird, Leon Brill, Michael Broughton, Bob B. Buckley, David A. Buell, Tim Burger, Brian Burkett, Nicholas Bushnell, Zijun Chen, Ben Chiaro , et al. (138 additional authors not shown)

Abstract: Measurement has a special role in quantum theory: by collapsing the wavefunction it can enable phenomena such as teleportation and thereby alter the "arrow of time" that constrains unitary evolution. When integrated in many-body dynamics, measurements can lead to emergent patterns of quantum information in space-time that go beyond established paradigms for characterizing phases, either in or out… ▽ More Measurement has a special role in quantum theory: by collapsing the wavefunction it can enable phenomena such as teleportation and thereby alter the "arrow of time" that constrains unitary evolution. When integrated in many-body dynamics, measurements can lead to emergent patterns of quantum information in space-time that go beyond established paradigms for characterizing phases, either in or out of equilibrium. On present-day NISQ processors, the experimental realization of this physics is challenging due to noise, hardware limitations, and the stochastic nature of quantum measurement. Here we address each of these experimental challenges and investigate measurement-induced quantum information phases on up to 70 superconducting qubits. By leveraging the interchangeability of space and time, we use a duality mapping, to avoid mid-circuit measurement and access different manifestations of the underlying phases -- from entanglement scaling to measurement-induced teleportation -- in a unified way. We obtain finite-size signatures of a phase transition with a decoding protocol that correlates the experimental measurement record with classical simulation data. The phases display sharply different sensitivity to noise, which we exploit to turn an inherent hardware limitation into a useful diagnostic. Our work demonstrates an approach to realize measurement-induced physics at scales that are at the limits of current NISQ processors. △ Less

Submitted 17 October, 2023; v1 submitted 8 March, 2023; originally announced March 2023.

Journal ref: Nature 622, 481-486 (2023)

Comparative study of magnetocaloric properties for Gd$^{3+}$ compounds with different frustrated lattice geometries

Authors: EliseAnne C. Koskelo, Paromita Mukherjee, Cheng Liu, Alice C. Sackville Hamilton, Harapan S. Ong, M. E. Zhitomirsky, Claudio Castelnovo, Siân E. Dutton

Abstract: As materials with suppressed ordering temperatures and enhanced ground state entropies, frustrated magnetic oxides are ideal candidates for cryogenic magnetocaloric refrigeration. While previous materials design has focused on tuning the magnetic moments, their interactions, and density of moments on the lattice, there has been relatively little attention to frustrated lattices. Prior theoretical… ▽ More As materials with suppressed ordering temperatures and enhanced ground state entropies, frustrated magnetic oxides are ideal candidates for cryogenic magnetocaloric refrigeration. While previous materials design has focused on tuning the magnetic moments, their interactions, and density of moments on the lattice, there has been relatively little attention to frustrated lattices. Prior theoretical work has shown that the magnetocaloric cooling rate at the saturation field is proportional to a macroscopic number of soft mode excitations that arise due to the classical ground state degeneracy. The number of these modes is directly determined by the geometry of the frustrating lattice. For corner-sharing geometries, the pyrochlore has 50\% more modes than the garnet and kagome lattices, whereas the edge-sharing \emph{fcc} has only a subextensive number of soft modes. Here, we study the role of soft modes in the magnetocaloric effect of four large-spin Gd$^{3+}$ ($L=0$, $J=S=7/2$) Heisenberg antiferromagnets on a kagome, garnet, pyrochlore, and \emph{fcc} lattice. By comparing measurements of the magnetic entropy change $ΔS_m$ of these materials at fields up to $9$~T with predictions using mean-field theory and Monte Carlo simulations, we are able to understand the relative importance of spin correlations and quantization effects. We observe that tuning the value of the nearest neighbor coupling has a more dominant contribution to the magnetocaloric entropy change in the liquid-He cooling regime ($2$-$20$~K), rather than tuning the number of soft mode excitations. Our results inform future materials design in terms of dimensionality, degree of magnetic frustration, and lattice geometry. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Comments: 15 pages, 14 figures

Saturation of spiral instabilities in disk galaxies

Authors: Chris Hamilton

Abstract: Spiral density waves can arise in galactic disks as linear instabilities of the underlying stellar distribution function. Such an instability grows exponentially in amplitude at some fixed growth rate $β$ before saturating nonlinearly. However, the mechanisms behind saturation, and the resulting saturated spiral amplitude, have received little attention. Here we argue that one important saturation… ▽ More Spiral density waves can arise in galactic disks as linear instabilities of the underlying stellar distribution function. Such an instability grows exponentially in amplitude at some fixed growth rate $β$ before saturating nonlinearly. However, the mechanisms behind saturation, and the resulting saturated spiral amplitude, have received little attention. Here we argue that one important saturation mechanism is the nonlinear trapping of stars near the spiral's corotation resonance. Under this mechanism, we show analytically that an $m$-armed spiral instability will saturate when the libration frequency of resonantly trapped orbits reaches $ω_\mathrm{lib} \sim \mathrm{a\,\, few}\times m^{1/2} β$. For a galaxy with a flat rotation curve this implies a maximum relative spiral surface density $\vert δΣ/Σ_0\vert \sim \mathrm{a\,\,few} \times (β/Ω_\mathrm{p})^2 \cot α$, where $Ω_\mathrm{p}$ is the spiral pattern speed and $α$ is its pitch angle. This result is in reasonable agreement with recent $N$-body simulations, and suggests that spirals driven by internally-generated instabilities reach relative amplitudes of at most a few tens of percent; higher amplitude spirals, like in M51 and NGC 1300, are likely caused by very strong bars and/or tidal perturbations. △ Less

Submitted 26 January, 2024; v1 submitted 13 February, 2023; originally announced February 2023.

Comments: Revised version, accepted for publication in MNRAS

Stream Fanning and Bifurcations: Observable Signatures of Resonances in Stellar Stream Morphology

Authors: Tomer D. Yavetz, Kathryn V. Johnston, Sarah Pearson, Adrian M. Price-Whelan, Chris Hamilton

Abstract: Recent observations have revealed a trove of unexpected morphological features in many of the Milky Way's stellar streams. Explanations for such features include time-dependent deformations of the Galactic gravitational potential, local disruptions induced by dark matter substructure, and special configurations of the streams' progenitors. In this paper, we study how these morphologies can also ar… ▽ More Recent observations have revealed a trove of unexpected morphological features in many of the Milky Way's stellar streams. Explanations for such features include time-dependent deformations of the Galactic gravitational potential, local disruptions induced by dark matter substructure, and special configurations of the streams' progenitors. In this paper, we study how these morphologies can also arise in certain static, non-spherical gravitational potentials that host a subset of resonantly-trapped orbit families. The transitions, or separatrices, between these orbit families mark abrupt discontinuities in the orbital structure of the potential. We develop a novel numerical approach for measuring the libration frequencies of resonant and near-resonant orbits, and apply it to study the evolution of stellar streams on these orbits. We reveal two distinct morphological features that arise in streams on near-resonant orbits: fans, that come about due to a large spread in the libration frequencies near a separatrix; and bifurcations, that arise when a separatrix splits the orbital distribution of the stellar stream between two (or more) distinct orbit families. We demonstrate that these effects can arise in some Milky Way streams for certain choices of the dark matter halo potential, and discuss how this might be used to probe and constrain the global shape of the Milky Way's gravitational potential. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: 24 pages, 8 figures. Submitted to ApJ; comments welcome!

Overcoming leakage in scalable quantum error correction

Authors: Kevin C. Miao, Matt McEwen, Juan Atalaya, Dvir Kafri, Leonid P. Pryadko, Andreas Bengtsson, Alex Opremcak, Kevin J. Satzinger, Zijun Chen, Paul V. Klimov, Chris Quintana, Rajeev Acharya, Kyle Anderson, Markus Ansmann, Frank Arute, Kunal Arya, Abraham Asfaw, Joseph C. Bardin, Alexandre Bourassa, Jenna Bovaird, Leon Brill, Bob B. Buckley, David A. Buell, Tim Burger, Brian Burkett , et al. (92 additional authors not shown)

Abstract: Leakage of quantum information out of computational states into higher energy states represents a major challenge in the pursuit of quantum error correction (QEC). In a QEC circuit, leakage builds over time and spreads through multi-qubit interactions. This leads to correlated errors that degrade the exponential suppression of logical error with scale, challenging the feasibility of QEC as a path… ▽ More Leakage of quantum information out of computational states into higher energy states represents a major challenge in the pursuit of quantum error correction (QEC). In a QEC circuit, leakage builds over time and spreads through multi-qubit interactions. This leads to correlated errors that degrade the exponential suppression of logical error with scale, challenging the feasibility of QEC as a path towards fault-tolerant quantum computation. Here, we demonstrate the execution of a distance-3 surface code and distance-21 bit-flip code on a Sycamore quantum processor where leakage is removed from all qubits in each cycle. This shortens the lifetime of leakage and curtails its ability to spread and induce correlated errors. We report a ten-fold reduction in steady-state leakage population on the data qubits encoding the logical state and an average leakage population of less than $1 \times 10^{-3}$ throughout the entire device. The leakage removal process itself efficiently returns leakage population back to the computational basis, and adding it to a code circuit prevents leakage from inducing correlated error across cycles, restoring a fundamental assumption of QEC. With this demonstration that leakage can be contained, we resolve a key challenge for practical QEC at scale. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: Main text: 7 pages, 5 figures

Purification-based quantum error mitigation of pair-correlated electron simulations

Authors: T. E. O'Brien, G. Anselmetti, F. Gkritsis, V. E. Elfving, S. Polla, W. J. Huggins, O. Oumarou, K. Kechedzhi, D. Abanin, R. Acharya, I. Aleiner, R. Allen, T. I. Andersen, K. Anderson, M. Ansmann, F. Arute, K. Arya, A. Asfaw, J. Atalaya, D. Bacon, J. C. Bardin, A. Bengtsson, S. Boixo, G. Bortoli, A. Bourassa , et al. (151 additional authors not shown)

Abstract: An important measure of the development of quantum computing platforms has been the simulation of increasingly complex physical systems. Prior to fault-tolerant quantum computing, robust error mitigation strategies are necessary to continue this growth. Here, we study physical simulation within the seniority-zero electron pairing subspace, which affords both a computational stepping stone to a ful… ▽ More An important measure of the development of quantum computing platforms has been the simulation of increasingly complex physical systems. Prior to fault-tolerant quantum computing, robust error mitigation strategies are necessary to continue this growth. Here, we study physical simulation within the seniority-zero electron pairing subspace, which affords both a computational stepping stone to a fully correlated model, and an opportunity to validate recently introduced ``purification-based'' error-mitigation strategies. We compare the performance of error mitigation based on doubling quantum resources in time (echo verification) or in space (virtual distillation), on up to $20$ qubits of a superconducting qubit quantum processor. We observe a reduction of error by one to two orders of magnitude below less sophisticated techniques (e.g. post-selection); the gain from error mitigation is seen to increase with the system size. Employing these error mitigation strategies enables the implementation of the largest variational algorithm for a correlated chemistry system to-date. Extrapolating performance from these results allows us to estimate minimum requirements for a beyond-classical simulation of electronic structure. We find that, despite the impressive gains from purification-based error mitigation, significant hardware improvements will be required for classically intractable variational chemistry simulations. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: 10 pages, 13 page supplementary material, 12 figures. Experimental data available at https://doi.org/10.5281/zenodo.7225821

Journal ref: Nat. Phys. (2023)

Non-Abelian braiding of graph vertices in a superconducting processor

Authors: Trond I. Andersen, Yuri D. Lensky, Kostyantyn Kechedzhi, Ilya Drozdov, Andreas Bengtsson, Sabrina Hong, Alexis Morvan, Xiao Mi, Alex Opremcak, Rajeev Acharya, Richard Allen, Markus Ansmann, Frank Arute, Kunal Arya, Abraham Asfaw, Juan Atalaya, Ryan Babbush, Dave Bacon, Joseph C. Bardin, Gina Bortoli, Alexandre Bourassa, Jenna Bovaird, Leon Brill, Michael Broughton, Bob B. Buckley , et al. (144 additional authors not shown)

Abstract: Indistinguishability of particles is a fundamental principle of quantum mechanics. For all elementary and quasiparticles observed to date - including fermions, bosons, and Abelian anyons - this principle guarantees that the braiding of identical particles leaves the system unchanged. However, in two spatial dimensions, an intriguing possibility exists: braiding of non-Abelian anyons causes rotatio… ▽ More Indistinguishability of particles is a fundamental principle of quantum mechanics. For all elementary and quasiparticles observed to date - including fermions, bosons, and Abelian anyons - this principle guarantees that the braiding of identical particles leaves the system unchanged. However, in two spatial dimensions, an intriguing possibility exists: braiding of non-Abelian anyons causes rotations in a space of topologically degenerate wavefunctions. Hence, it can change the observables of the system without violating the principle of indistinguishability. Despite the well developed mathematical description of non-Abelian anyons and numerous theoretical proposals, the experimental observation of their exchange statistics has remained elusive for decades. Controllable many-body quantum states generated on quantum processors offer another path for exploring these fundamental phenomena. While efforts on conventional solid-state platforms typically involve Hamiltonian dynamics of quasi-particles, superconducting quantum processors allow for directly manipulating the many-body wavefunction via unitary gates. Building on predictions that stabilizer codes can host projective non-Abelian Ising anyons, we implement a generalized stabilizer code and unitary protocol to create and braid them. This allows us to experimentally verify the fusion rules of the anyons and braid them to realize their statistics. We then study the prospect of employing the anyons for quantum computation and utilize braiding to create an entangled state of anyons encoding three logical qubits. Our work provides new insights about non-Abelian braiding and - through the future inclusion of error correction to achieve topological protection - could open a path toward fault-tolerant quantum computing. △ Less

Submitted 31 May, 2023; v1 submitted 18 October, 2022; originally announced October 2022.

Readout of a quantum processor with high dynamic range Josephson parametric amplifiers

Authors: T. C. White, Alex Opremcak, George Sterling, Alexander Korotkov, Daniel Sank, Rajeev Acharya, Markus Ansmann, Frank Arute, Kunal Arya, Joseph C. Bardin, Andreas Bengtsson, Alexandre Bourassa, Jenna Bovaird, Leon Brill, Bob B. Buckley, David A. Buell, Tim Burger, Brian Burkett, Nicholas Bushnell, Zijun Chen, Ben Chiaro, Josh Cogan, Roberto Collins, Alexander L. Crook, Ben Curtin , et al. (69 additional authors not shown)

Abstract: We demonstrate a high dynamic range Josephson parametric amplifier (JPA) in which the active nonlinear element is implemented using an array of rf-SQUIDs. The device is matched to the 50 $Ω$ environment with a Klopfenstein-taper impedance transformer and achieves a bandwidth of 250-300 MHz, with input saturation powers up to -95 dBm at 20 dB gain. A 54-qubit Sycamore processor was used to benchmar… ▽ More We demonstrate a high dynamic range Josephson parametric amplifier (JPA) in which the active nonlinear element is implemented using an array of rf-SQUIDs. The device is matched to the 50 $Ω$ environment with a Klopfenstein-taper impedance transformer and achieves a bandwidth of 250-300 MHz, with input saturation powers up to -95 dBm at 20 dB gain. A 54-qubit Sycamore processor was used to benchmark these devices, providing a calibration for readout power, an estimate of amplifier added noise, and a platform for comparison against standard impedance matched parametric amplifiers with a single dc-SQUID. We find that the high power rf-SQUID array design has no adverse effect on system noise, readout fidelity, or qubit dephasing, and we estimate an upper bound on amplifier added noise at 1.6 times the quantum limit. Lastly, amplifiers with this design show no degradation in readout fidelity due to gain compression, which can occur in multi-tone multiplexed readout with traditional JPAs. △ Less

Submitted 22 November, 2022; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: 10 pages, 10 figures

Journal ref: Appl. Phys. Lett. 122, 014001 (2023)

Galactic bar resonances with diffusion: an analytic model with implications for bar-dark matter halo dynamical friction

Authors: Chris Hamilton, Elizabeth A. Tolman, Lev Arzamasskiy, Vinícius N. Duarte

Abstract: The secular evolution of disk galaxies is largely driven by resonances between the orbits of 'particles' (stars or dark matter) and the rotation of non-axisymmetric features (spiral arms or a bar). Such resonances may also explain kinematic and photometric features observed in the Milky Way and external galaxies. In simplified cases, these resonant interactions are well understood: for instance, t… ▽ More The secular evolution of disk galaxies is largely driven by resonances between the orbits of 'particles' (stars or dark matter) and the rotation of non-axisymmetric features (spiral arms or a bar). Such resonances may also explain kinematic and photometric features observed in the Milky Way and external galaxies. In simplified cases, these resonant interactions are well understood: for instance, the dynamics of a test particle trapped near a resonance of a steadily rotating bar is easily analyzed using the angle-action tools pioneered by Binney, Monari and others. However, such treatments do not address the stochasticity and messiness inherent to real galaxies - effects which have, with few exceptions, been previously explored only with complex N-body simulations. In this paper, we propose a simple kinetic equation describing the distribution function of particles near an orbital resonance with a rigidly rotating bar, allowing for diffusion of the particles' slow actions. We solve this equation for various values of the dimensionless diffusion strength $Δ$, and then apply our theory to the calculation of bar-halo dynamical friction. For $Δ= 0$ we recover the classic result of Tremaine & Weinberg that friction ultimately vanishes, owing to the phase-mixing of resonant orbits. However, for $Δ> 0$ we find that diffusion suppresses phase-mixing, leading to a finite torque. Our results suggest that stochasticity - be it physical or numerical - tends to increase bar-halo friction, and that bars in cosmological simulations might experience significant artificial slowdown, even if the numerical two-body relaxation time is much longer than a Hubble time. △ Less

Submitted 18 May, 2023; v1 submitted 7 August, 2022; originally announced August 2022.

Comments: Accepted version, to appear in ApJ

Suppressing quantum errors by scaling a surface code logical qubit

Authors: Rajeev Acharya, Igor Aleiner, Richard Allen, Trond I. Andersen, Markus Ansmann, Frank Arute, Kunal Arya, Abraham Asfaw, Juan Atalaya, Ryan Babbush, Dave Bacon, Joseph C. Bardin, Joao Basso, Andreas Bengtsson, Sergio Boixo, Gina Bortoli, Alexandre Bourassa, Jenna Bovaird, Leon Brill, Michael Broughton, Bob B. Buckley, David A. Buell, Tim Burger, Brian Burkett, Nicholas Bushnell , et al. (132 additional authors not shown)

Abstract: Practical quantum computing will require error rates that are well below what is achievable with physical qubits. Quantum error correction offers a path to algorithmically-relevant error rates by encoding logical qubits within many physical qubits, where increasing the number of physical qubits enhances protection against physical errors. However, introducing more qubits also increases the number… ▽ More Practical quantum computing will require error rates that are well below what is achievable with physical qubits. Quantum error correction offers a path to algorithmically-relevant error rates by encoding logical qubits within many physical qubits, where increasing the number of physical qubits enhances protection against physical errors. However, introducing more qubits also increases the number of error sources, so the density of errors must be sufficiently low in order for logical performance to improve with increasing code size. Here, we report the measurement of logical qubit performance scaling across multiple code sizes, and demonstrate that our system of superconducting qubits has sufficient performance to overcome the additional errors from increasing qubit number. We find our distance-5 surface code logical qubit modestly outperforms an ensemble of distance-3 logical qubits on average, both in terms of logical error probability over 25 cycles and logical error per cycle ($2.914\%\pm 0.016\%$ compared to $3.028\%\pm 0.023\%$). To investigate damaging, low-probability error sources, we run a distance-25 repetition code and observe a $1.7\times10^{-6}$ logical error per round floor set by a single high-energy event ($1.6\times10^{-7}$ when excluding this event). We are able to accurately model our experiment, and from this model we can extract error budgets that highlight the biggest challenges for future systems. These results mark the first experimental demonstration where quantum error correction begins to improve performance with increasing qubit number, illuminating the path to reaching the logical error rates required for computation. △ Less

Submitted 20 July, 2022; v1 submitted 13 July, 2022; originally announced July 2022.

Comments: Main text: 6 pages, 4 figures. v2: Update author list, references, Fig. S12, Table IV

Anatomy of a slow merger: dissecting secularly-driven inspirals of LIGO/Virgo gravitational wave sources

Authors: Chris Hamilton, Roman R. Rafikov

Abstract: The dozens of compact object mergers detected by LIGO/Virgo raise a key theoretical question: how do initially wide binaries shrink sufficiently quickly that they are able to merge via gravitational wave (GW) radiation within a Hubble time? One promising class of answers involves secular driving of binary eccentricity by some external tidal perturbation. This perturbation can arise due to the pres… ▽ More The dozens of compact object mergers detected by LIGO/Virgo raise a key theoretical question: how do initially wide binaries shrink sufficiently quickly that they are able to merge via gravitational wave (GW) radiation within a Hubble time? One promising class of answers involves secular driving of binary eccentricity by some external tidal perturbation. This perturbation can arise due to the presence of a tertiary point mass, in which case the system exhibits Lidov-Kozai (LK) dynamics, or it can stem from the tidal field of the stellar cluster in which the binary orbits. While these secular tide-driven mechanisms have been studied exhaustively in the case of no GW emission, when GWs are included the dynamical behavior is still incompletely understood. In this paper we consider compact object binaries driven to merger via high eccentricity excitation by (doubly-averaged, test-particle quadrupole level) cluster tides - which includes LK-driven mergers as a special case - and include the effects of both general relativistic precession and GW emission. We provide for the first time an analytical understanding of the different evolutionary stages of the binary's semimajor axis, secular oscillation timescale, and phase space structure all the way to merger. Our results will inform future population synthesis calculations of compact object binary mergers from hierarchical triples and stellar clusters. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Comments: 30 pages, 12 figures. To be submitted to ApJ. Comments welcome

Formation of robust bound states of interacting microwave photons

Authors: Alexis Morvan, Trond I. Andersen, Xiao Mi, Charles Neill, Andre Petukhov, Kostyantyn Kechedzhi, Dmitry Abanin, Rajeev Acharya, Frank Arute, Kunal Arya, Abraham Asfaw, Juan Atalaya, Ryan Babbush, Dave Bacon, Joseph C. Bardin, Joao Basso, Andreas Bengtsson, Gina Bortoli, Alexandre Bourassa, Jenna Bovaird, Leon Brill, Michael Broughton, Bob B. Buckley, David A. Buell, Tim Burger , et al. (125 additional authors not shown)

Abstract: Systems of correlated particles appear in many fields of science and represent some of the most intractable puzzles in nature. The computational challenge in these systems arises when interactions become comparable to other energy scales, which makes the state of each particle depend on all other particles. The lack of general solutions for the 3-body problem and acceptable theory for strongly cor… ▽ More Systems of correlated particles appear in many fields of science and represent some of the most intractable puzzles in nature. The computational challenge in these systems arises when interactions become comparable to other energy scales, which makes the state of each particle depend on all other particles. The lack of general solutions for the 3-body problem and acceptable theory for strongly correlated electrons shows that our understanding of correlated systems fades when the particle number or the interaction strength increases. One of the hallmarks of interacting systems is the formation of multi-particle bound states. In a ring of 24 superconducting qubits, we develop a high fidelity parameterizable fSim gate that we use to implement the periodic quantum circuit of the spin-1/2 XXZ model, an archetypal model of interaction. By placing microwave photons in adjacent qubit sites, we study the propagation of these excitations and observe their bound nature for up to 5 photons. We devise a phase sensitive method for constructing the few-body spectrum of the bound states and extract their pseudo-charge by introducing a synthetic flux. By introducing interactions between the ring and additional qubits, we observe an unexpected resilience of the bound states to integrability breaking. This finding goes against the common wisdom that bound states in non-integrable systems are unstable when their energies overlap with the continuum spectrum. Our work provides experimental evidence for bound states of interacting photons and discovers their stability beyond the integrability limit. △ Less

Submitted 21 December, 2022; v1 submitted 10 June, 2022; originally announced June 2022.

Comments: 7 pages + 15 pages supplements

Journal ref: Nature 612, 240-245 (2022)

Merging Black Holes: Assessing the Performance of Two Analytic Gravitational Waves Models

Authors: Dillon Buskirk, Maria C. Babiuc Hamilton

Abstract: Merging black holes produce the loudest signal in the detectors. However, this is the most difficult signal to accurately predict with analytical techniques. Only computer simulations can account for the nonlinear physics during the collision, but they are inherently complex, costly, and affected by numerical errors. In order to bypass this problem, two analytical models for the merger have been d… ▽ More Merging black holes produce the loudest signal in the detectors. However, this is the most difficult signal to accurately predict with analytical techniques. Only computer simulations can account for the nonlinear physics during the collision, but they are inherently complex, costly, and affected by numerical errors. In order to bypass this problem, two analytical models for the merger have been developed: the Implicit Rotating Source (IRS) and the newer Backwards one Body (BoB). In this work, we assess the performance of the BoB model by comparing it with the older IRS model and with the numerical data, identifying its strengths and weaknesses. Our main finding reveals discrepancies in amplitude, but overall excellent accord in frequency. The BoB model is comparable with the IRS and NR simulations, having the added advantage that it depends only indirectly on numerical data, it accounts for spin, and it offers a seamless fit with the analytical formalisms for the inspiral. By independently evaluating and testing those models, we bring evidence of their reproducibility, thus upholding high scientific standards, and make available our implementation, enabling readers to evaluate our results themselves. △ Less

Submitted 1 March, 2023; v1 submitted 29 May, 2022; originally announced May 2022.

Comments: 28 pages, 8 figures

Journal ref: Journal of Cosmology and Astroparticle Physics (2023) 005

Wide twin binaries are extremely eccentric: evidence of twin binary formation in circumbinary disks

Authors: Hsiang-Chih Hwang, Kareem El-Badry, Hans-Walter Rix, Chris Hamilton, Yuan-Sen Ting, Nadia Zakamska

Abstract: The Gaia mission recently revealed an excess population of equal-mass "twin" wide binaries, with mass ratio $q\gtrsim 0.95$, extending to separations of at least 1000 AU. The origin of this population is an enigma: twin binaries are thought to form via correlated accretion in circumbinary disks, but the typical observed protostellar disks have radii of $\sim100$ AU, far smaller than the separation… ▽ More The Gaia mission recently revealed an excess population of equal-mass "twin" wide binaries, with mass ratio $q\gtrsim 0.95$, extending to separations of at least 1000 AU. The origin of this population is an enigma: twin binaries are thought to form via correlated accretion in circumbinary disks, but the typical observed protostellar disks have radii of $\sim100$ AU, far smaller than the separations of the widest twins. Here, we infer the eccentricity distribution of wide twins from the distribution of their $v$-$r$ angles, i.e., the angle between the components' separation and relative velocity vectors. We find that wide twins must be on extremely eccentric orbits. For the excess-twin population at 400-1000 AU, we infer a near-delta function excess of high-eccentricity system, with eccentricity $0.95 \lesssim e \leq 1$. These high eccentricities for wide twins imply pericenter distances of order $10$ AU and suggest that their orbits were scattered via dynamical interactions in their birth environments, consistent with a scenario in which twins are born in circumbinary disks and subsequently widened. These results further establish twin wide binaries as a distinct population and imply that wide twins can be used as a probe of the dynamical history of stellar populations. △ Less

Submitted 7 July, 2022; v1 submitted 11 May, 2022; originally announced May 2022.

Comments: accepted for publication in The Astrophysical Journal Letter

Eccentric Pairs: Analytic Gravitational Waves from Binary Black Holes in Elliptic Orbits

Authors: Dillon Buskirk, Maria C. Babiuc Hamilton

Abstract: Gravitational waves (GW) from eccentric binaries have intricate signals encoding important features about the location, creation and evolution of the sources. Eccentricity shortens the merger time, making the emitted GW statistically predominant in the observed data once detectors will reach the required sensitivity. We present a novel implementation of fully analytical GW templates from eccentric… ▽ More Gravitational waves (GW) from eccentric binaries have intricate signals encoding important features about the location, creation and evolution of the sources. Eccentricity shortens the merger time, making the emitted GW statistically predominant in the observed data once detectors will reach the required sensitivity. We present a novel implementation of fully analytical GW templates from eccentric binary black hole (BBH) mergers within the \texttt{Wolfram Mathematica} software. We increase the accuracy by identifying and minimizing the possible source of errors. We start with an overview of the physics involved in eccentric mergers, then assemble the strain for the inspiral by employing up to six post-Newtonian (PN) corrections. We complete the eccentric inspiral with the quasi-circular Backwards one Body (BoB) merger model in frequency, amplitude and phase, then we build the hybrid GW strain for the whole evolution of the binary. For low eccentricity we reach coincidence in the overlap, with no ambiguity in the time interval, a remarkable improvement from the usual matching techniques. For high-eccentricity we compensate for the implicit quasi-circular assumption of the BOB approach, by introducing a small rescaling in amplitude. Our streamlined implementation is relevant for the new field of GW astronomy and is straightforward to understand, use and extend, offering researchers in the field a valuable open resource tool. △ Less

Submitted 27 October, 2022; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: 27 pages, 10 figures, accepted for publication in the International Journal of Modern Physics D

Journal ref: International Journal of Modern Physics D, Vol. 32, No. 02, 2250138 (2023)

On the phase-mixed eccentricity and inclination distributions of wide binaries in the Galaxy

Authors: Chris Hamilton

Abstract: Modern observational surveys allow us to probe the phase space distribution function (DF) of wide binaries in the Solar neighbourhood. This DF exhibits non-trivial features, in particular a superthermal distribution of eccentricities for semimajor axes $a\gtrsim 10^3$AU. To interpret such features we must first understand how the binary DF is affected by dynamical perturbations, which typically fa… ▽ More Modern observational surveys allow us to probe the phase space distribution function (DF) of wide binaries in the Solar neighbourhood. This DF exhibits non-trivial features, in particular a superthermal distribution of eccentricities for semimajor axes $a\gtrsim 10^3$AU. To interpret such features we must first understand how the binary DF is affected by dynamical perturbations, which typically fall into two classes: (i) stochastic kicks from passing stars, molecular clouds, etc. and (ii) secular torques from the Galactic tide. Here we isolate effect (ii) and calculate the time-asymptotic, phase-mixed DF for an ensemble of wide binaries under quadrupole-order tides. For binaries wide enough that the phase-mixing assumption is valid, none of our results depend explicitly on semimajor axes, masses, etc. We show that unless the initial DF is both isotropic in binary orientation and thermal in eccentricity, then the final phase-mixed DF is always both anisotropic and non-thermal. However, the only way to produce a superthermal DF under phase mixing is for the initial DF to itself be superthermal. △ Less

Submitted 23 April, 2022; v1 submitted 2 February, 2022; originally announced February 2022.

Comments: Accepted for publication in ApJ Letters

On the Energy Dependence of Galactic Cosmic Ray Anisotropies in the Very Local Interstellar Medium

Authors: Romina Nikoukar, Matthew E. Hill, Lawrence Brown, Stamatios M. Krimigis, Robert B. Decker, Konstantinos Dialynas, Jozsef Kota, Edmond C. Roelof, Scott Lasley, Douglas C. Hamilton, Vladimir Florinski, Joe Giacalone, John Richardson, Merav Opher

Abstract: We report on the energy dependence of galactic cosmic rays (GCRs) in the very local interstellar medium (VLISM) as measured by the Low Energy Charged Particle (LECP) instrument on the Voyager 1 (V1) spacecraft. The LECP instrument includes a dual-ended solid state detector particle telescope mechanically scanning through 360 deg across eight equally-spaced angular sectors. As reported previously,… ▽ More We report on the energy dependence of galactic cosmic rays (GCRs) in the very local interstellar medium (VLISM) as measured by the Low Energy Charged Particle (LECP) instrument on the Voyager 1 (V1) spacecraft. The LECP instrument includes a dual-ended solid state detector particle telescope mechanically scanning through 360 deg across eight equally-spaced angular sectors. As reported previously, LECP measurements showed a dramatic increase in GCR intensities for all sectors of the >=211 MeV count rate (CH31) at the V1 heliopause (HP) crossing in 2012, however, since then the count rate data have demonstrated systematic episodes of intensity decrease for particles around 90° pitch angle. To shed light on the energy dependence of these GCR anisotropies over a wide range of energies, we use V1 LECP count rate and pulse height analyzer (PHA) data from >=211 MeV channel together with lower energy LECP channels. Our analysis shows that while GCR anisotropies are present over a wide range of energies, there is a decreasing trend in the amplitude of second-order anisotropy with increasing energy during anisotropy episodes. A stronger pitch-angle scattering at the higher velocities is argued as a potential cause for this energy dependence. A possible cause for this velocity dependence arising from weak rigidity dependence of the scattering mean free path and resulting velocity-dominated scattering rate is discussed. This interpretation is consistent with a recently reported lack of corresponding GCR electron anisotropies. △ Less

Submitted 19 January, 2022; originally announced January 2022.

Essential Metadata for 3D BRAIN Microscopy

Authors: Alexander J. Ropelewski, Megan A. Rizzo, Jason R. Swedlow, Jan Huisken, Pavel Osten, Neda Khanjani, Kurt Weiss, Vesselina Bakalov, Michelle Engle, Lauren Gridley, Michelle Krzyzanowski, Tom Madden, Deborah Maiese, Justin Waterfield, David Williams, Carol Hamilton, Wayne Huggins

Abstract: Recent advances in fluorescence microscopy techniques and tissue clearing, labeling, and staining provide unprecedented opportunities to investigate brain structure and function. These experiments' images make it possible to catalog brain cell types and define their location, morphology, and connectivity in a native context, leading to a better understanding of normal development and disease etiol… ▽ More Recent advances in fluorescence microscopy techniques and tissue clearing, labeling, and staining provide unprecedented opportunities to investigate brain structure and function. These experiments' images make it possible to catalog brain cell types and define their location, morphology, and connectivity in a native context, leading to a better understanding of normal development and disease etiology. Consistent annotation of metadata is needed to provide the context necessary to understand, reuse, and integrate these data. This report describes an effort to establish metadata standards for 3D microscopy datasets for use by the Brain Research through Advancing Innovative Neurotechnologies (BRAIN) Initiative and the neuroscience research community. These standards were built on existing efforts and developed with input from the brain microscopy community to promote adoption. The resulting Essential Metadata for 3D BRAIN Microscopy includes 91 fields organized into seven categories: Contributors, Funders, Publication, Instrument, Dataset, Specimen, and Image. Adoption of these metadata standards will ensure that investigators receive credit for their work, promote data reuse, facilitate downstream analysis of shared data, and encourage collaboration. △ Less

Submitted 18 May, 2021; originally announced May 2021.

Comments: 10 pages, 1 table

High efficacy of layered controls for reducing transmission of airborne pathogens

Authors: Laura Fierce, Alison Robey, Cathrine Hamilton

Abstract: To optimize strategies for curbing the transmission of airborne pathogens, the efficacy of three key controls -- face masks, ventilation, and physical distancing -- must be well understood. In this study we used the Quadrature-based model of Respiratory Aerosol and Droplets to quantify the reduction in exposure to airborne pathogens from various combinations of controls. For each combination of co… ▽ More To optimize strategies for curbing the transmission of airborne pathogens, the efficacy of three key controls -- face masks, ventilation, and physical distancing -- must be well understood. In this study we used the Quadrature-based model of Respiratory Aerosol and Droplets to quantify the reduction in exposure to airborne pathogens from various combinations of controls. For each combination of controls, we simulated thousands of scenarios that represent the tremendous variability in factors governing airborne transmission and the efficacy of mitigation strategies. While the efficacy of any individual control was highly variable among scenarios, combining universal mask-wearing with distancing of 1~m or more reduced the median exposure by more than 99\% relative to a close, unmasked conversation, with further reductions if ventilation is also enhanced. The large reductions in exposure to airborne pathogens translated to large reductions in the risk of initial infection in a new host. These findings suggest that layering controls is highly effective for reducing transmission of airborne pathogens and will be critical for curbing outbreaks of novel viruses in the future. △ Less

Submitted 21 July, 2022; v1 submitted 6 April, 2021; originally announced April 2021.

Journal ref: Indoor Air, 2022

Simulating near-field enhancement in transmission of airborne viruses with a quadrature-based model

Authors: Laura Fierce, Alison Robey, Cathrine Hamilton

Abstract: Airborne viruses, such as influenza, tuberculosis, and SARS-CoV-2, are transmitted through virus-laden particles expelled when an infectious person sneezes, coughs, talks, or breathes. These virus-laden particles are more highly concentrated in the expiratory jet of an infectious person than in a well-mixed room, but this near-field enhancement in virion exposure has not been well quantified. Tran… ▽ More Airborne viruses, such as influenza, tuberculosis, and SARS-CoV-2, are transmitted through virus-laden particles expelled when an infectious person sneezes, coughs, talks, or breathes. These virus-laden particles are more highly concentrated in the expiratory jet of an infectious person than in a well-mixed room, but this near-field enhancement in virion exposure has not been well quantified. Transmission of airborne viruses depends on factors that are inherently variable and, in many cases, poorly constrained, and quantifying this uncertainty requires large ensembles of model simulations that span the variability in input parameters. However, models that are well-suited to simulate the near-field evolution of respiratory particles are also computationally expensive, which limits the exploration of parametric uncertainty. In order to perform many simulations that span the wide variability in factors governing transmission, we developed the Quadrature-based model of Respiratory Aerosol and Droplets (QuaRAD). QuaRAD is an efficient framework for simulating the evolution of virus-laden particles after they are expelled from an infectious person, their deposition to the nasal cavity of a susceptible person, and the subsequent risk of initial infection. We simulated 10,000 scenarios to quantify the risk of initial infection by a particular virus, SARS-CoV-2. The predicted risk of infection was highly variable among scenarios and, in each scenario, was strongly enhanced near the infectious individual. In more than 50% of scenarios, the physical distancing needed to avoid near-field enhancements in airborne transmission was beyond the recommended safe distance of two meters (six feet) if the infectious person is not wearing a mask, though this distance defining the near-field extent was also highly variable among scenarios. △ Less

Submitted 2 April, 2021; originally announced April 2021.

Journal ref: Indoor Air, 2021

Resonant and non-resonant relaxation of globular clusters

Authors: Jean-Baptiste Fouvry, Chris Hamilton, Simon Rozier, Christophe Pichon

Abstract: Globular clusters contain a finite number of stars. As a result, they inevitably undergo secular evolution (`relaxation') causing their mean distribution function (DF) to evolve on long timescales. On one hand, this long-term evolution may be interpreted as driven by the accumulation of local deflections along each star's mean field trajectory -- so-called `non-resonant relaxation'. On the other h… ▽ More Globular clusters contain a finite number of stars. As a result, they inevitably undergo secular evolution (`relaxation') causing their mean distribution function (DF) to evolve on long timescales. On one hand, this long-term evolution may be interpreted as driven by the accumulation of local deflections along each star's mean field trajectory -- so-called `non-resonant relaxation'. On the other hand, it can be thought of as driven by non-local, collectively dressed and resonant couplings between stellar orbits, a process termed `resonant relaxation'. In this paper we consider a model globular cluster represented by a spherical, isotropic isochrone DF, and compare in detail the predictions of both resonant and non-resonant relaxation theories against tailored direct $N$-body simulations. In the space of orbital actions (namely the radial action and total angular momentum), we find that both resonant and non-resonant theories predict the correct morphology for the secular evolution of the cluster's DF, although non-resonant theory over-estimates the amplitude of the relaxation rate by a factor ${\sim 2}$. We conclude that the secular relaxation of hot isotropic spherical clusters is not dominated by collectively amplified large-scale potential fluctuations, despite the existence of a strong ${\ell = 1}$ damped mode. Instead, collective amplification affects relaxation only marginally even on the largest scales. The predicted contributions to relaxation from smaller scale fluctuations are essentially the same from resonant and non-resonant theories. △ Less

Submitted 18 March, 2021; originally announced March 2021.

Comments: 26 pages, 14 figures, submitted to MNRAS

Secular dynamics of binaries in stellar clusters -- III. Doubly-averaged dynamics in the presence of general relativistic precession

Authors: Chris Hamilton, Roman R. Rafikov

Abstract: Secular evolution of binaries driven by an external (tidal) potential is a classic astrophysical problem. Tidal perturbations can arise due to an external point mass, as in the Lidov-Kozai (LK) theory of hierarchical triples, or due to an extended stellar system (e.g. galaxy or globular cluster) in which the binary resides. For many applications, general-relativistic (GR) apsidal precession is imp… ▽ More Secular evolution of binaries driven by an external (tidal) potential is a classic astrophysical problem. Tidal perturbations can arise due to an external point mass, as in the Lidov-Kozai (LK) theory of hierarchical triples, or due to an extended stellar system (e.g. galaxy or globular cluster) in which the binary resides. For many applications, general-relativistic (GR) apsidal precession is important, and has been accounted for in some LK calculations. Here we generalise and extend these studies by exploring in detail the effect of GR precession on (quadrupole-level) tidal evolution of binaries orbiting in arbitrary axisymmetric potentials (which includes LK theory as a special case). We study the (doubly-averaged) orbital dynamics for arbitrary strengths of GR and binary initial conditions and uncover entirely new phase space morphologies with important implications for the binary orbital evolution. We also explore how GR precession affects secular evolution of binary orbital elements when the binary reaches high eccentricity ($e\to 1$) and delineate several different dynamical regimes. Our results are applicable to a variety of astrophysical systems. In particular, they can be used to understand the high-eccentricity behaviour of (cluster) tide-driven compact object mergers -- i.e. LIGO/Virgo gravitational wave sources -- for which GR effects are crucial. △ Less

Submitted 19 May, 2021; v1 submitted 30 November, 2020; originally announced November 2020.

Comments: Final version, accepted for publication in MNRAS

Noise and waves: a unified kinetic theory for stellar systems

Authors: Chris Hamilton, Tobias Heinemann

Abstract: The traditional Chandrasekhar picture of the slow relaxation of stellar systems assumes that stars' orbits are only modified by occasional, uncorrelated, two-body flyby encounters with other stars. However, the long-range nature of gravity means that in reality large numbers of stars can behave collectively. In stable systems this collective behaviour (i) amplifies the noisy fluctuations in the sy… ▽ More The traditional Chandrasekhar picture of the slow relaxation of stellar systems assumes that stars' orbits are only modified by occasional, uncorrelated, two-body flyby encounters with other stars. However, the long-range nature of gravity means that in reality large numbers of stars can behave collectively. In stable systems this collective behaviour (i) amplifies the noisy fluctuations in the system's gravitational potential, effectively 'dressing' the two-body (star-star) encounters, and (ii) allows the system to support large-scale density waves (a.k.a. normal modes) which decay through resonant wave-star interactions. If the relaxation of the system is dominated by effect (i) then it is described by the Balescu-Lenard (BL) kinetic theory. Meanwhile if (ii) dominates, one must describe relaxation using quasilinear (QL) theory, though in the stellar-dynamical context the full set of QL equations has never been presented. Moreover, in some systems like open clusters and galactic disks, both (i) and (ii) might be important. Here we present for the first time the equations of a unified kinetic theory of stellar systems in angle-action variables that accounts for both effects (i) and (ii) simultaneously. We derive the equations in a heuristic, physically-motivated fashion and work in the simplest possible regime by accounting only for very weakly damped waves. This unified theory is effectively a superposition of BL and QL theories, both of which are recovered in appropriate limits. The theory is a first step towards a comprehensive description of those stellar systems for which neither the QL or BL theory will suffice. △ Less

Submitted 30 November, 2020; originally announced November 2020.

Comments: 9 pages, 2 figures. Comments welcome

Quantum transport in a combined kicked rotor and quantum walk system

Authors: Adrian Ortega, Thomas Gorin, Craig S. Hamilton

Abstract: We present a theoretical and numerical study of the competition between two opposite interference effects, namely interference-induced ballistic transport on one hand, and strong (Anderson) localization on the other. While the former effect allows for resistance free transport, the latter brings the transport to a complete halt. As a model system, we consider the quantum kicked rotor, where strong… ▽ More We present a theoretical and numerical study of the competition between two opposite interference effects, namely interference-induced ballistic transport on one hand, and strong (Anderson) localization on the other. While the former effect allows for resistance free transport, the latter brings the transport to a complete halt. As a model system, we consider the quantum kicked rotor, where strong localization is observed in the discrete momentum coordinate. In this model, we introduce the ballistic transport in the form of a Hadamard quantum walk in that momentum coordinate. The two transport mechanisms are combined by alternating the corresponding Floquet operators. Extending the corresponding calculation for the kicked rotor, we estimate the classical diffusion coefficient for thecombined dynamics. Another argument, based on the introduction of an effective Heisenberg time should then allow to estimate the localization time and the localization length. While this is known to work reasonably well in the kicked rotor case, we find that it fails in our case. While the combined dynamics still shows localization, it takes place at much larger times and shows much larger localization lengths than predicted. Finally, we combine the kicked rotor with other types of quantum walks, namely diffusive and localizing quantum walks. In the diffusive case, the localizing dynamics of the kicked rotor is completely canceled and we get pure diffusion. In the case of the localizing quantum walk, the combined system remains localized, but with a larger localization length. △ Less

Submitted 3 March, 2021; v1 submitted 23 November, 2020; originally announced November 2020.

Evidence for geologically recent explosive volcanism in Elysium Planitia, Mars

Authors: David G. Horvath, Pranabendu Moitra, Christopher W. Hamilton, Robert A. Craddock, Jeffrey C. Andrews-Hanna

Abstract: Volcanic activity on Mars peaked during the Noachian and Hesperian periods but has continued since then in isolated locales. Elysium Planitia hosts numerous young, fissure-fed flood lavas with ages ranging from approximately 500 to 2.5 million years (Ma). We present evidence for what may be the youngest volcanic deposit yet documented on Mars: a low albedo, high thermal inertia, high-calcium pyrox… ▽ More Volcanic activity on Mars peaked during the Noachian and Hesperian periods but has continued since then in isolated locales. Elysium Planitia hosts numerous young, fissure-fed flood lavas with ages ranging from approximately 500 to 2.5 million years (Ma). We present evidence for what may be the youngest volcanic deposit yet documented on Mars: a low albedo, high thermal inertia, high-calcium pyroxene-rich deposit distributed symmetrically around a segment of the Cerberus Fossae fissure system in Elysium Planitia. This deposit is similar to features interpreted as pyroclastic deposits on the Moon and Mercury. However, unlike previously documented lava flows in Elysium Planitia, this feature is morphologically consistent with a fissure-fed pyroclastic deposit, mantling the surrounding lava flows with a thickness on the order of tens of cm over most of the deposit and a volume of 1.1-2.8E7 cubic meters. Thickness and volume estimates are consistent with tephra fall deposits on Earth. Stratigraphic relationships indicate a relative age younger than the surrounding volcanic plains and the Zunil impact crater (~0.1-1 Ma), with crater counting suggesting an absolute model age of 53 to 210 ka. This young age implies that if this deposit is of volcanic origin then the Cerberus Fossae region may not be extinct and Mars may still be volcanically active today. This interpretation is consistent with the identification of seismicity in this region by the Interior Exploration using Seismic Investigations, Geodesy, and Heat Transport (InSight) lander, and has additional implications for astrobiology and the source of transient atmospheric methane. △ Less

Submitted 11 November, 2020; originally announced November 2020.

Comments: 41 pages, 11 figures

Lunar Near-Surface Volatile Sample Return

Authors: Igor Aleinov, Michael J. Way, Christopher W. Hamilton, James W. Head

Abstract: The origin, distribution, depth and volume of lunar volatiles remain open questions. One of the possible sources of Moon's volatiles is their volcanic outgassing during the peak of lunar volcanic activity ~3.5 Ga. This same outgassing would also produce a tenuous transient atmosphere which would promote the delivery of volatiles from the volcanic sources to the polar cold traps. Though such an atm… ▽ More The origin, distribution, depth and volume of lunar volatiles remain open questions. One of the possible sources of Moon's volatiles is their volcanic outgassing during the peak of lunar volcanic activity ~3.5 Ga. This same outgassing would also produce a tenuous transient atmosphere which would promote the delivery of volatiles from the volcanic sources to the polar cold traps. Though such an atmosphere could have played an important role in the evolution of the Moon, little is known about it due to high uncertainty level in the mechanisms involved. The only reliable proxy for the ancient lunar atmosphere are the primordial volatiles deposited by it, which are expected to be preserved in the polar cold traps, and could be studied through sample return. In this white paper we therefore advocate that a volatile sample return from the Moon's polar cold traps should be a fundamental part of the Artemis program. △ Less

Submitted 10 September, 2020; originally announced September 2020.

Comments: 2 pages, A white paper submitted to the Artemis Science Definition Team of NASA's Science Mission Directorate

A simple, heuristic derivation of the Balescu-Lenard kinetic equation for stellar systems

Authors: Chris Hamilton

Abstract: The unshielded nature of gravity means that stellar systems are inherently inhomogeneous. As a result, stars do not move in straight lines. This obvious fact severely complicates the kinetic theory of stellar systems because position and velocity turn out to be poor coordinates with which to describe stellar orbits - instead, one must use angle-action variables. Moreover, the slow relaxation of st… ▽ More The unshielded nature of gravity means that stellar systems are inherently inhomogeneous. As a result, stars do not move in straight lines. This obvious fact severely complicates the kinetic theory of stellar systems because position and velocity turn out to be poor coordinates with which to describe stellar orbits - instead, one must use angle-action variables. Moreover, the slow relaxation of star clusters and galaxies can be enhanced or suppressed by collective interactions ('polarisation' effects) involving many stars simultaneously. These collective effects are also present in plasmas; in that case, they are accounted for by the Balescu-Lenard (BL) equation, which is a kinetic equation in velocity space. Recently several authors have shown how to account for both inhomogeneity and collective effects in the kinetic theory of stellar systems by deriving an angle-action generalisation of the BL equation. Unfortunately their derivations are long and complicated, involving multiple coordinate transforms, contour integrals in the complex plane, and so on. On the other hand, Rostoker's superposition principle allows one to pretend that a long-range interacting $N$-body system, such as a plasma or star cluster, consists merely of uncorrelated particles that are 'dressed' by polarisation clouds. In this paper we use Rostoker's principle to provide a simple, intuitive derivation of the BL equation for stellar systems which is much shorter than others in the literature. It also allows us to straightforwardly connect the BL picture of self-gravitating kinetics to the classical 'two-body relaxation' theory of uncorrelated flybys pioneered by Chandrasekhar. △ Less

Submitted 8 December, 2020; v1 submitted 14 July, 2020; originally announced July 2020.

Comments: Accepted for publication in MNRAS

Compact object binary mergers driven by cluster tides: a new channel for LIGO/Virgo gravitational wave events

Authors: Chris Hamilton, Roman R. Rafikov

Abstract: The detections of gravitational waves produced in mergers of binary black holes (BH) and neutron stars (NS) by LIGO/Virgo have stimulated interest in the origin of the progenitor binaries. Dense stellar systems - globular and nuclear star clusters - are natural sites of compact object binary formation and evolution towards merger. Here we explore a new channel for the production of binary mergers… ▽ More The detections of gravitational waves produced in mergers of binary black holes (BH) and neutron stars (NS) by LIGO/Virgo have stimulated interest in the origin of the progenitor binaries. Dense stellar systems - globular and nuclear star clusters - are natural sites of compact object binary formation and evolution towards merger. Here we explore a new channel for the production of binary mergers in clusters, in which the tidal field of the cluster secularly drives the binary to high eccentricity (even in the absence of a central massive black hole) until gravitational wave emission becomes important. We employ the recently developed secular theory of cluster tide-driven binary evolution to compute present day merger rates for BH-BH, NS-BH and NS-NS binaries, varying cluster potential and central concentration of the binary population (but ignoring cluster evolution and stellar flybys for now). Unlike other mechanisms, this new dynamical channel can produce a significant number of mergers out to cluster-centric distances of several pc. For NS-NS binaries we find merger rates in the range $0.01-0.07$ Gpc$^{-3}$ yr$^{-1}$ from globular clusters and $0.1-0.2$ Gpc$^{-3}$ yr$^{-1}$ from cusped nuclear clusters. For NS-BH and BH-BH binaries we find small merger rates from globular clusters, but a rate of $0.1 - 0.2$ Gpc$^{-3}$ yr$^{-1}$ from cusped nuclear clusters, contributing to the observed LIGO/Virgo rate at the level of several per cent. Therefore, cluster tide-driven mergers constitute a new channel that can be further explored with current and future gravitational wave detectors. △ Less

Submitted 16 August, 2019; v1 submitted 1 July, 2019; originally announced July 2019.

Comments: Final version, accepted for publication in ApJ Letters (9 pages)

Secular dynamics of binaries in stellar clusters II: dynamical evolution

Authors: Chris Hamilton, Roman R. Rafikov

Abstract: Dense stellar clusters are natural sites for the origin and evolution of exotic objects such as relativistic binaries (potential gravitational wave sources), blue stragglers, etc. We investigate the secular dynamics of a binary system driven by the global tidal field of an axisymmetric stellar cluster in which the binary orbits. In a companion paper (Hamilton & Rafikov 2019a) we developed a genera… ▽ More Dense stellar clusters are natural sites for the origin and evolution of exotic objects such as relativistic binaries (potential gravitational wave sources), blue stragglers, etc. We investigate the secular dynamics of a binary system driven by the global tidal field of an axisymmetric stellar cluster in which the binary orbits. In a companion paper (Hamilton & Rafikov 2019a) we developed a general Hamiltonian framework describing such systems. The effective (doubly-averaged) Hamiltonian derived there encapsulates all information about the tidal potential experienced by the binary in its orbit around the cluster in a single parameter $Γ$. Here we provide a thorough exploration of the phase-space of the corresponding secular problem as $Γ$ is varied. We find that for $Γ> 1/5$ the phase-space structure and the evolution of binary orbital element are qualitatively similar to the Lidov-Kozai problem. However, this is only one of four possible regimes, because the dynamics are qualitatively changed by bifurcations at $Γ= 1/5,0,-1/5$. We show how the dynamics are altered in each regime and calculate characteristics such as secular evolution timescale, maximum possible eccentricity, etc. We verify the predictions of our doubly-averaged formalism numerically and find it to be very accurate when its underlying assumptions are fulfilled, typically meaning that the secular timescale should exceed the period of the binary around the cluster by $\gtrsim 10-10^2$ (depending on the cluster potential and binary orbit). Our results may be relevant for understanding the nature of a variety of exotic systems harboured by stellar clusters. △ Less

Submitted 16 August, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

Comments: 25 pages, 14 figures. Final version accepted for publication in MNRAS