Showing 1–50 of 122 results for author: Shaw, L

Universal fluctuations and noise learning from Hilbert-space ergodicity

Authors: Adam L. Shaw, Daniel K. Mark, Joonhee Choi, Ran Finkelstein, Pascal Scholl, Soonwon Choi, Manuel Endres

Abstract: Systems reaching thermal equilibrium are ubiquitous. For classical systems, this phenomenon is typically understood statistically through ergodicity in phase space, but translating this to quantum systems is a long-standing problem of interest. Recently a quantum notion of ergodicity has been proposed, namely that isolated, global quantum states uniformly explore their available state space, dubbe… ▽ More Systems reaching thermal equilibrium are ubiquitous. For classical systems, this phenomenon is typically understood statistically through ergodicity in phase space, but translating this to quantum systems is a long-standing problem of interest. Recently a quantum notion of ergodicity has been proposed, namely that isolated, global quantum states uniformly explore their available state space, dubbed Hilbert-space ergodicity. Here we observe signatures of this process with an experimental Rydberg quantum simulator and various numerical models, before generalizing to the case of a local quantum system interacting with its environment. For a closed system, where the environment is a complementary subsystem, we predict and observe a smooth quantum-to-classical transition in that observables progress from large, quantum fluctuations to small, Gaussian fluctuations as the bath size grows. This transition is universal on a quantitative level amongst a wide range of systems, including those at finite temperature, those with itinerant particles, and random circuits. Then, we consider the case of an open system interacting noisily with an external environment. We predict the statistics of observables under largely arbitrary noise channels including those with correlated errors, allowing us to discriminate candidate error models both for continuous Hamiltonian time evolution and for digital random circuits. Ultimately our results clarify the role of ergodicity in quantum dynamics, with fundamental and practical consequences. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: ALS and DKM contributed equally

Report number: MIT-CTP/5693

A Maximum Entropy Principle in Deep Thermalization and in Hilbert-Space Ergodicity

Authors: Daniel K. Mark, Federica Surace, Andreas Elben, Adam L. Shaw, Joonhee Choi, Gil Refael, Manuel Endres, Soonwon Choi

Abstract: We report universal statistical properties displayed by ensembles of pure states that naturally emerge in quantum many-body systems. Specifically, two classes of state ensembles are considered: those formed by i) the temporal trajectory of a quantum state under unitary evolution or ii) the quantum states of small subsystems obtained by partial, local projective measurements performed on their comp… ▽ More We report universal statistical properties displayed by ensembles of pure states that naturally emerge in quantum many-body systems. Specifically, two classes of state ensembles are considered: those formed by i) the temporal trajectory of a quantum state under unitary evolution or ii) the quantum states of small subsystems obtained by partial, local projective measurements performed on their complements. These cases respectively exemplify the phenomena of "Hilbert-space ergodicity" and "deep thermalization." In both cases, the resultant ensembles are defined by a simple principle: the distributions of pure states have maximum entropy, subject to constraints such as energy conservation, and effective constraints imposed by thermalization. We present and numerically verify quantifiable signatures of this principle by deriving explicit formulae for all statistical moments of the ensembles; proving the necessary and sufficient conditions for such universality under widely-accepted assumptions; and describing their measurable consequences in experiments. We further discuss information-theoretic implications of the universality: our ensembles have maximal information content while being maximally difficult to interrogate, establishing that generic quantum state ensembles that occur in nature hide (scramble) information as strongly as possible. Our results generalize the notions of Hilbert-space ergodicity to time-independent Hamiltonian dynamics and deep thermalization from infinite to finite effective temperature. Our work presents new perspectives to characterize and understand universal behaviors of quantum dynamics using statistical and information theoretic tools. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Report number: MIT-CTP/5692

Universal quantum operations and ancilla-based readout for tweezer clocks

Authors: Ran Finkelstein, Richard Bing-Shiun Tsai, Xiangkai Sun, Pascal Scholl, Su Direkci, Tuvia Gefen, Joonhee Choi, Adam L. Shaw, Manuel Endres

Abstract: Enhancing the precision of measurements by harnessing entanglement is a long-sought goal in the field of quantum metrology. Yet attaining the best sensitivity allowed by quantum theory in the presence of noise is an outstanding challenge, requiring optimal probe-state generation and readout strategies. Neutral atom optical clocks, leading systems for measuring time, have shown recent progress in t… ▽ More Enhancing the precision of measurements by harnessing entanglement is a long-sought goal in the field of quantum metrology. Yet attaining the best sensitivity allowed by quantum theory in the presence of noise is an outstanding challenge, requiring optimal probe-state generation and readout strategies. Neutral atom optical clocks, leading systems for measuring time, have shown recent progress in terms of entanglement generation, but currently lack the control capabilities to realize such schemes. Here we show universal quantum operations and ancilla-based readout for ultranarrow optical transitions of neutral atoms. Our demonstration in a tweezer clock platform enables a circuit-based approach to quantum metrology with neutral atom optical clocks. To this end, we demonstrate two-qubit entangling gates with 99.62% fidelity via Rydberg interactions and dynamical connectivity for optical clock qubits, which we combine with local addressing to implement universally programmable quantum circuits. Using this approach, we generate a near-optimal entangled probe state, a cascade of Greenberger-Horne-Zeilinger (GHZ) states of different sizes, and perform dual-quadrature GHZ readout. We also show repeated fast phase detection with non-destructive conditional reset of clock qubits and minimal dead time between repetitions by implementing ancilla-based quantum logic spectroscopy (QLS) for neutral atoms. Finally, we extend this to multi-qubit parity checks and measurement-based, heralded, Bell state preparation. Our work lays the foundation for hybrid processor-clock devices with neutral atoms and more generally points to a future of practical applications for quantum processors linked with quantum sensors. △ Less

Submitted 8 July, 2024; v1 submitted 25 February, 2024; originally announced February 2024.

Comments: RF, RBST, and XS contributed equally to this work

The Influence of Laser Focusing Conditions on the Direct Laser Acceleration of Electrons

Authors: H. Tang, K. Tangtartharakul, R. Babjak, I-L. Yeh, F. Albert, H. Chen, P. T. Campbell, Y. Ma, P. M. Nilson, B. K. Russell, J. L. Shaw, A. G. R. Thomas, M. Vranic, A. V. Arefiev, L. Willingale

Abstract: Direct Laser Acceleration (DLA) of electrons during a high-energy, picosecond laser interaction with an underdense plasma has been demonstrated to be substantially enhanced by controlling the laser focusing geometry. Experiments using the OMEGA EP facility measured electrons accelerated to maximum energies exceeding 120 times the ponderomotive energy under certain laser focusing, pulse energy, and… ▽ More Direct Laser Acceleration (DLA) of electrons during a high-energy, picosecond laser interaction with an underdense plasma has been demonstrated to be substantially enhanced by controlling the laser focusing geometry. Experiments using the OMEGA EP facility measured electrons accelerated to maximum energies exceeding 120 times the ponderomotive energy under certain laser focusing, pulse energy, and plasma density conditions. Two-dimensional particle-in-cell simulations show that the laser focusing conditions alter the laser field evolution, channel fields generation, and electron oscillation, all of which contribute to the final electron energies. The optimal laser focusing condition occurs when the transverse oscillation amplitude of the accelerated electron in the channel fields matches the laser beam width, resulting in efficient energy gain. Through this observation, a simple model was developed to calculate the optimal laser focal spot size in more general conditions and is validated by experimental data. △ Less

Submitted 12 February, 2024; originally announced February 2024.

Comments: 14 pages, 6 figures

Erasure-cooling, control, and hyper-entanglement of motion in optical tweezers

Authors: Pascal Scholl, Adam L. Shaw, Ran Finkelstein, Richard Bing-Shiun Tsai, Joonhee Choi, Manuel Endres

Abstract: We demonstrate how motional degrees of freedom in optical tweezers can be used as quantum information carriers. To this end, we first implement a species-agnostic cooling mechanism via conversion of motional excitations into erasures - errors with a known location - reminiscent of Maxwell's demon thought experiment. We find that this cooling mechanism fundamentally outperforms idealized traditiona… ▽ More We demonstrate how motional degrees of freedom in optical tweezers can be used as quantum information carriers. To this end, we first implement a species-agnostic cooling mechanism via conversion of motional excitations into erasures - errors with a known location - reminiscent of Maxwell's demon thought experiment. We find that this cooling mechanism fundamentally outperforms idealized traditional sideband cooling, which we experimentally demonstrate in specific scenarios. By coherently manipulating the motional state, we perform mid-circuit readout and mid-circuit erasure detection of an optical qubit via local shelving into motional superposition states. We finally entangle the motion of two atoms in separate tweezers, and utilize this to generate hyper-entanglement by preparing a simultaneous Bell state of motional and optical qubits. This work shows how controlling motion enriches the toolbox of quantum information processing with neutral atoms, and opens unique prospects for metrology enhanced by mid-circuit readout and a large class of quantum operations enabled via hyper-entanglement. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: PS, ALS and RF contributed equally to this work

Generalized extrapolation methods based on compositions of a basic 2nd-order scheme

Authors: Sergio Blanes, Fernando Casas, Luke Shaw

Abstract: We propose new linear combinations of compositions of a basic second-order scheme with appropriately chosen coefficients to construct higher order numerical integrators for differential equations. They can be considered as a generalization of extrapolation methods and multi-product expansions. A general analysis is provided and new methods up to order 8 are built and tested. The new approach is sh… ▽ More We propose new linear combinations of compositions of a basic second-order scheme with appropriately chosen coefficients to construct higher order numerical integrators for differential equations. They can be considered as a generalization of extrapolation methods and multi-product expansions. A general analysis is provided and new methods up to order 8 are built and tested. The new approach is shown to reduce the latency problem when implemented in a parallel environment and leads to schemes that are significantly more efficient than standard extrapolation when the linear combination is delayed by a number of steps. △ Less

Submitted 23 April, 2024; v1 submitted 20 November, 2023; originally announced November 2023.

Comments: 17 figures

Keyword-optimized Template Insertion for Clinical Information Extraction via Prompt-based Learning

Authors: Eugenia Alleva, Isotta Landi, Leslee J Shaw, Erwin Böttinger, Thomas J Fuchs, Ipek Ensari

Abstract: Clinical note classification is a common clinical NLP task. However, annotated data-sets are scarse. Prompt-based learning has recently emerged as an effective method to adapt pre-trained models for text classification using only few training examples. A critical component of prompt design is the definition of the template (i.e. prompt text). The effect of template position, however, has been insu… ▽ More Clinical note classification is a common clinical NLP task. However, annotated data-sets are scarse. Prompt-based learning has recently emerged as an effective method to adapt pre-trained models for text classification using only few training examples. A critical component of prompt design is the definition of the template (i.e. prompt text). The effect of template position, however, has been insufficiently investigated. This seems particularly important in the clinical setting, where task-relevant information is usually sparse in clinical notes. In this study we develop a keyword-optimized template insertion method (KOTI) and show how optimizing position can improve performance on several clinical tasks in a zero-shot and few-shot training setting. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Security in Cryptocurrency

Authors: Chelsea Medina, Lily Shaw, Dissy Vargas, Sundar Krishnan

Abstract: This paper discusses the mechanisms of cryptocurrency, the idea of using security in the system, and the popularity of it. To begin, the authors provide a background on cryptocurrency and how it works. The authors understand that while most people may be familiar with the concept, they may not know how it works. Next, the authors discuss the security of cryptocurrency in-depth within the paper. Th… ▽ More This paper discusses the mechanisms of cryptocurrency, the idea of using security in the system, and the popularity of it. To begin, the authors provide a background on cryptocurrency and how it works. The authors understand that while most people may be familiar with the concept, they may not know how it works. Next, the authors discuss the security of cryptocurrency in-depth within the paper. The authors also provide examples of attacks on cryptocurrency systems to show the vulnerabilities within the system. Lastly, the authors discuss the popularity of the system to further express the need for security in cryptocurrency. △ Less

Submitted 16 October, 2023; originally announced October 2023.

Dephasingless laser wakefield acceleration in the bubble regime

Authors: Kyle G. Miller, Jacob R. Pierce, Manfred V. Ambat, Jessica L. Shaw, Kale Weichman, Warren B. Mori, Dustin H. Froula, John P. Palastro

Abstract: Laser wakefield accelerators (LWFAs) have electric fields that are orders of magnitude larger than those of conventional accelerators, promising an attractive, small-scale alternative for next-generation light sources and lepton colliders. The maximum energy gain in a single-stage LWFA is limited by dephasing, which occurs when the trapped particles outrun the accelerating phase of the wakefield.… ▽ More Laser wakefield accelerators (LWFAs) have electric fields that are orders of magnitude larger than those of conventional accelerators, promising an attractive, small-scale alternative for next-generation light sources and lepton colliders. The maximum energy gain in a single-stage LWFA is limited by dephasing, which occurs when the trapped particles outrun the accelerating phase of the wakefield. Here, we demonstrate that a single space-time structured laser pulse can be used for ionization injection and electron acceleration over many dephasing lengths in the bubble regime. Simulations of a dephasingless laser wakefield accelerator driven by a 6.2-J laser pulse show 25 pC of injected charge accelerated over 20 dephasing lengths (1.3 cm) to a maximum energy of 2.1 GeV. The space-time structured laser pulse features an ultrashort, programmable-trajectory focus. Accelerating the focus, reducing the focused spot-size variation, and mitigating unwanted self-focusing stabilize the electron acceleration, which improves beam quality and leads to projected energy gains of 125 GeV in a single, sub-meter stage driven by a 500-J pulse. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Comments: 18 pages, 4 figures

Benchmarking highly entangled states on a 60-atom analog quantum simulator

Authors: Adam L. Shaw, Zhuo Chen, Joonhee Choi, Daniel K. Mark, Pascal Scholl, Ran Finkelstein, Andreas Elben, Soonwon Choi, Manuel Endres

Abstract: Quantum systems have entered a competitive regime where classical computers must make approximations to represent highly entangled quantum states. However, in this beyond-classically-exact regime, fidelity comparisons between quantum and classical systems have so far been limited to digital quantum devices, and it remains unsolved how to estimate the actual entanglement content of experiments. Her… ▽ More Quantum systems have entered a competitive regime where classical computers must make approximations to represent highly entangled quantum states. However, in this beyond-classically-exact regime, fidelity comparisons between quantum and classical systems have so far been limited to digital quantum devices, and it remains unsolved how to estimate the actual entanglement content of experiments. Here we perform fidelity benchmarking and mixed-state entanglement estimation with a 60-atom analog Rydberg quantum simulator, reaching a high entanglement entropy regime where exact classical simulation becomes impractical. Our benchmarking protocol involves extrapolation from comparisons against an approximate classical algorithm, introduced here, with varying entanglement limits. We then develop and demonstrate an estimator of the experimental mixed-state entanglement, finding our experiment is competitive with state-of-the-art digital quantum devices performing random circuit evolution. Finally, we compare the experimental fidelity against that achieved by various approximate classical algorithms, and find that only the algorithm we introduce is able to keep pace with the experiment on the classical hardware we employ. Our results enable a new paradigm for evaluating the ability of both analog and digital quantum devices to generate entanglement in the beyond-classically-exact regime, and highlight the evolving divide between quantum and classical systems. △ Less

Submitted 4 December, 2023; v1 submitted 15 August, 2023; originally announced August 2023.

Comments: ALS, ZC, JC, and DKM contributed equally

Journal ref: Nature 628, 71-77 (2024)

Programmable and arbitrary-trajectory ultrafast flying focus pulses

Authors: M. V. Ambat, J. L. Shaw, J. J. Pigeon, K. G. Miller, T. T. Simpson, D. H. Froula, J. P. Palastro

Abstract: "Flying focus" techniques produce laser pulses with dynamic focal points that travels distances much greater than a Rayleigh length. The implementation of these techniques in laser-based applications requires the design of optical configurations that can both extend the focal range and structure the radial group delay. This article describes a method for designing optical configurations that produ… ▽ More "Flying focus" techniques produce laser pulses with dynamic focal points that travels distances much greater than a Rayleigh length. The implementation of these techniques in laser-based applications requires the design of optical configurations that can both extend the focal range and structure the radial group delay. This article describes a method for designing optical configurations that produce ultrashort flying focus pulses with arbitrary-trajectory focal points. The method is illustrated by several examples that employ an axiparabola for extending the focal range and either a reflective echelon or a deformable mirror-spatial light modulator pair for structuring the radial group delay. The latter configuration enables rapid exploration and optimization of flying foci, which could be ideal for experiments. △ Less

Submitted 11 July, 2023; originally announced July 2023.

Optimal dimensions of cone and pyramid moth-eye structures for $\mathrm{SiO}_{2}$ windows

Authors: Chaoran Tu, Jonathan Hu, Curtis R. Menyuk, Thomas F. Carruthers, L. Brandon Shaw, Lynda E. Busse, Jasbinder S. Sanghera

Abstract: We computationally investigate the transmission efficiency through moth-eye nanostructures that are fabricated on $\mathrm{SiO}_{2}$ windows in the wavelength range from 0.4 to 2 $\mathrm{μm}$. We investigated both truncated cones and truncated pyramids, and we varied the height, bottom width, and top width of these shapes in order to maximize the transmission efficiency. We found that there is no… ▽ More We computationally investigate the transmission efficiency through moth-eye nanostructures that are fabricated on $\mathrm{SiO}_{2}$ windows in the wavelength range from 0.4 to 2 $\mathrm{μm}$. We investigated both truncated cones and truncated pyramids, and we varied the height, bottom width, and top width of these shapes in order to maximize the transmission efficiency. We found that there is no substantial difference in transmission between truncated cone and pyramid structures. Using the constraints from the current achievable experimental limits, a relatively uniform transmission coefficient of larger than 98.8% can be obtained from 0.4 $\mathrm{μm}$ to 2 $\mathrm{μm}$. These transmission results are only 0.4% in absolute value lower than the transmission of a structure that is not constrained by current experimental limits. △ Less

Submitted 24 June, 2023; originally announced June 2023.

Stability and Fluctuations in Complex Ecological Systems

Authors: Eric Forgoston, Sarah Day, Peter C. de Ruiter, Arjen Doelman, Nienke Hartemink, Alan Hastings, Lia Hemerik, Alexandru Hening, Josef Hofbauer, Sonia Kefi, David A. Kessler, Toni Klauschies, Christian Kuehn, Xiaoxiao Li, John C. Moore, Elly Morrien, Anje-Margriet Neutel, Jelena Pantel, Sebastian J. Schreiber, Leah B. Shaw, Nadav Shnerb, Eric Siero, Laura S. Storch, Michael A. S. Thorne, Ingrid van de Leemput , et al. (2 additional authors not shown)

Abstract: From 08-12 August, 2022, 32 individuals participated in a workshop, Stability and Fluctuations in Complex Ecological Systems, at the Lorentz Center, located in Leiden, The Netherlands. An interdisciplinary dialogue between ecologists, mathematicians, and physicists provided a foundation of important problems to consider over the next 5-10 years. This paper outlines eight areas including (1) improv… ▽ More From 08-12 August, 2022, 32 individuals participated in a workshop, Stability and Fluctuations in Complex Ecological Systems, at the Lorentz Center, located in Leiden, The Netherlands. An interdisciplinary dialogue between ecologists, mathematicians, and physicists provided a foundation of important problems to consider over the next 5-10 years. This paper outlines eight areas including (1) improving our understanding of the effect of scale, both temporal and spatial, for both deterministic and stochastic problems; (2) clarifying the different terminologies and definitions used in different scientific fields; (3) developing a comprehensive set of data analysis techniques arising from different fields but which can be used together to improve our understanding of existing data sets; (4) having theoreticians/computational scientists collaborate closely with empirical ecologists to determine what new data should be collected; (5) improving our knowledge of how to protect and/or restore ecosystems; (6) incorporating socio-economic effects into models of ecosystems; (7) improving our understanding of the role of deterministic and stochastic fluctuations; (8) studying the current state of biodiversity at the functional level, taxa level and genome level. △ Less

Submitted 12 June, 2023; originally announced June 2023.

Comments: 22 pages

Erasure conversion in a high-fidelity Rydberg quantum simulator

Authors: Pascal Scholl, Adam L. Shaw, Richard Bing-Shiun Tsai, Ran Finkelstein, Joonhee Choi, Manuel Endres

Abstract: Minimizing and understanding errors is critical for quantum science, both in noisy intermediate scale quantum (NISQ) devices and for the quest towards fault-tolerant quantum computation. Rydberg arrays have emerged as a prominent platform in this context with impressive system sizes and proposals suggesting how error-correction thresholds could be significantly improved by detecting leakage errors… ▽ More Minimizing and understanding errors is critical for quantum science, both in noisy intermediate scale quantum (NISQ) devices and for the quest towards fault-tolerant quantum computation. Rydberg arrays have emerged as a prominent platform in this context with impressive system sizes and proposals suggesting how error-correction thresholds could be significantly improved by detecting leakage errors with single-atom resolution, a form of erasure error conversion. However, two-qubit entanglement fidelities in Rydberg atom arrays have lagged behind competitors and this type of erasure conversion is yet to be realized for matter-based qubits in general. Here we demonstrate both erasure conversion and high-fidelity Bell state generation using a Rydberg quantum simulator. We implement erasure conversion via fast imaging of alkaline-earth atoms, which leaves atoms in a metastable state unperturbed and yields additional information independent of the final qubit readout. When excising data with observed erasure errors, we achieve a lower-bound for the Bell state generation fidelity of ${\geq} 0.9971^{+10}_{-13}$, which improves to ${\geq}0.9985^{+7}_{-12}$ when correcting for remaining state preparation errors. We further demonstrate erasure conversion in a quantum simulation experiment for quasi-adiabatic preparation of long-range order across a quantum phase transition, where we explicitly differentiate erasure conversion of preparation and Rydberg decay errors. We unveil the otherwise hidden impact of these errors on the simulation outcome by evaluating correlations between erasures and the final readout as well as between erasures themselves. Our work demonstrates the capability for Rydberg-based entanglement to reach fidelities in the ${\sim} 0.999$ regime, with higher fidelities a question of technical improvements, and shows how erasure conversion can be utilized in NISQ devices. △ Less

Submitted 12 October, 2023; v1 submitted 5 May, 2023; originally announced May 2023.

Comments: PS and ALS contributed equally to this work

Journal ref: Nature 622, 273 (2023)

Multi-ensemble metrology by programming local rotations with atom movements

Authors: Adam L. Shaw, Ran Finkelstein, Richard Bing-Shiun Tsai, Pascal Scholl, Tai Hyun Yoon, Joonhee Choi, Manuel Endres

Abstract: Current optical atomic clocks do not utilize their resources optimally. In particular, an exponential gain in sensitivity could be achieved if multiple atomic ensembles were to be controlled or read-out individually, even without entanglement. However, controlling optical transitions locally remains an outstanding challenge for neutral atom based clocks and quantum computing platforms. Here we sho… ▽ More Current optical atomic clocks do not utilize their resources optimally. In particular, an exponential gain in sensitivity could be achieved if multiple atomic ensembles were to be controlled or read-out individually, even without entanglement. However, controlling optical transitions locally remains an outstanding challenge for neutral atom based clocks and quantum computing platforms. Here we show arbitrary, single-site addressing for an optical transition via sub-wavelength controlled moves of tweezer-trapped atoms, which we perform with $99.84(5)\%$ fidelity and with $0.1(2)\%$ crosstalk to non-addressed atoms. The scheme is highly robust as it relies only on relative position changes of tweezers and requires no additional addressing beams. Using this technique, we implement single-shot, dual-quadrature readout of Ramsey interferometry using two atomic ensembles simultaneously, and show an enhancement of the usable interrogation time at a given phase-slip error probability. Finally, we program a sequence which performs local dynamical decoupling during Ramsey evolution to evolve three ensembles with variable phase sensitivities, a key ingredient of optimal clock interrogation. Our results demonstrate the potential of fully programmable quantum optical clocks even without entanglement and could be combined with metrologically useful entangled states in the future. △ Less

Submitted 6 December, 2023; v1 submitted 29 March, 2023; originally announced March 2023.

Comments: ALS and RF contributed equally

Journal ref: Nature Physics (2024)

High energy and high power primary Li-CF$_x$ batteries enabled by the combined effects of the binder and the electrolyte

Authors: Haobin Huo, Leon L. Shaw, Karoly Nemeth

Abstract: Several effective methods have been developed recently to demonstrate simultaneous high energy and high power density in Li - carbon fluoride (CF$_x$) batteries. These methods can achieve as high as 1000 Wh/kg energy density at 60-70 kW/kg power density (40-50 C rate) in coin cells and 750 Wh/kg energy density at 12.5 kW/kg power density (20 C rate) in pouch cells. This performance is made possibl… ▽ More Several effective methods have been developed recently to demonstrate simultaneous high energy and high power density in Li - carbon fluoride (CF$_x$) batteries. These methods can achieve as high as 1000 Wh/kg energy density at 60-70 kW/kg power density (40-50 C rate) in coin cells and 750 Wh/kg energy density at 12.5 kW/kg power density (20 C rate) in pouch cells. This performance is made possible by ingenious nano-architecture design, controlled porosity, boron doping and electrolyte additives. In the present study, we show that a similarly great performance, 931 Wh/kg energy density at 59 kW/kg power density, can be achieved by using a polyacrylonitrile binder and a LiBF4 electrolyte in Li - graphite fluoride coin cells. We also demonstrate that the observed effect is the result of the right combination of the binder and the electrolyte. We propose that the mechanistic origin of the observed phenomena is an electro-catalytic effect by the polyacrylonitrile binder. While our proposed method has a competitive performance, it also offers a simple implementation and a scalable production of high energy and high power primary Li-CF$_x$ cells. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 13 pages, 8 figures, 1 table

Accurate simulation of direct laser acceleration in a laser wakefield accelerator

Authors: Kyle G. Miller, John P. Palastro, Jessica L. Shaw, Fei Li, Frank S. Tsung, Viktor K. Decyk, Warren B. Mori

Abstract: In a laser wakefield accelerator (LWFA), an intense laser pulse excites a plasma wave that traps and accelerates electrons to relativistic energies. When the pulse overlaps the accelerated electrons, it can enhance the energy gain through direct laser acceleration (DLA) by resonantly driving the betatron oscillations of the electrons in the plasma wave. The particle-in-cell (PIC) algorithm, althou… ▽ More In a laser wakefield accelerator (LWFA), an intense laser pulse excites a plasma wave that traps and accelerates electrons to relativistic energies. When the pulse overlaps the accelerated electrons, it can enhance the energy gain through direct laser acceleration (DLA) by resonantly driving the betatron oscillations of the electrons in the plasma wave. The particle-in-cell (PIC) algorithm, although often the tool of choice to study DLA, contains inherent errors due to numerical dispersion and the time staggering of the electric and magnetic fields. Further, conventional PIC implementations cannot reliably disentangle the fields of the plasma wave and laser pulse, which obscures interpretation of the dominant acceleration mechanism. Here, a customized field solver that reduces errors from both numerical dispersion and time staggering is used in conjunction with a field decomposition into azimuthal modes to perform PIC simulations of DLA in an LWFA. Comparisons with traditional PIC methods, model equations, and experimental data show improved accuracy with the customized solver and convergence with an order-of-magnitude fewer cells. The azimuthal-mode decomposition reveals that the most energetic electrons receive comparable energy from DLA and LWFA. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Comments: 10 pages, 5 figures, to submit to Physics of Plasmas

Dark-state enhanced loading of an optical tweezer array

Authors: Adam L. Shaw, Pascal Scholl, Ran Finklestein, Ivaylo S. Madjarov, Brandon Grinkemeyer, Manuel Endres

Abstract: Neutral atoms and molecules trapped in optical tweezers have become a prevalent resource for quantum simulation, computation, and metrology. However, the maximum achievable system sizes of such arrays are often limited by the stochastic nature of loading into optical tweezers, with a typical loading probability of only 50%. Here we present a species-agnostic method for dark-state enhanced loading… ▽ More Neutral atoms and molecules trapped in optical tweezers have become a prevalent resource for quantum simulation, computation, and metrology. However, the maximum achievable system sizes of such arrays are often limited by the stochastic nature of loading into optical tweezers, with a typical loading probability of only 50%. Here we present a species-agnostic method for dark-state enhanced loading (DSEL) based on real-time feedback, long-lived shelving states, and iterated array reloading. We demonstrate this technique with a 95-tweezer array of $^{88}$Sr atoms, achieving a maximum loading probability of 84.02(4)% and a maximum array size of 91 atoms in one dimension. Our protocol is complementary to, and compatible with, existing schemes for enhanced loading based on direct control over light-assisted collisions, and we predict it can enable close-to-unity filling for arrays of atoms or molecules. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Journal ref: Phys. Rev. Lett. 130, 193402 (2023)

A New Optimality Property of Strang's Splitting

Authors: Fernando Casas, Jesús María Sanz-Serna, Luke Shaw

Abstract: For systems of the form $\dot q = M^{-1} p$, $\dot p = -Aq+f(q)$, common in many applications, we analyze splitting integrators based on the (linear/nonlinear) split systems $\dot q = M^{-1} p$, $\dot p = -Aq$ and $\dot q = 0$, $\dot p = f(q)$. We show that the well-known Strang splitting is optimally stable in the sense that, when applied to a relevant model problem, it has a larger stability reg… ▽ More For systems of the form $\dot q = M^{-1} p$, $\dot p = -Aq+f(q)$, common in many applications, we analyze splitting integrators based on the (linear/nonlinear) split systems $\dot q = M^{-1} p$, $\dot p = -Aq$ and $\dot q = 0$, $\dot p = f(q)$. We show that the well-known Strang splitting is optimally stable in the sense that, when applied to a relevant model problem, it has a larger stability region than alternative integrators. This generalizes a well-known property of the common Störmer/Verlet/leapfrog algorithm, which of course arises from Strang splitting based on the (kinetic/potential) split systems $\dot q = M^{-1} p$, $\dot p = 0$ and $\dot q = 0$, $\dot p = -Aq+f(q)$. △ Less

Submitted 15 February, 2023; v1 submitted 13 October, 2022; originally announced October 2022.

Comments: 2 figures

Single-Shot Electron Radiography Using a Laser-Plasma Accelerator

Authors: G. Bruhaug, M. S. Freeman, H. G. Rinderknecht, L. P. Neukirch, C. H. Wilde, F. E Merrill, J. R. Rygg, M. S. Wei, G. W. Collins, J. L. Shaw

Abstract: Contact and projection electron radiography of static targets was demonstrated using a laser plasma accelerator driven by a kilojoule, picosecond class laser as a source of relativistic electrons with an average energy of 20 MeV. Objects with areal densities as high as 7.7 g/cm^2 were probed in materials ranging from plastic to tungsten, and radiographs with resolution as good as 90 micrometers we… ▽ More Contact and projection electron radiography of static targets was demonstrated using a laser plasma accelerator driven by a kilojoule, picosecond class laser as a source of relativistic electrons with an average energy of 20 MeV. Objects with areal densities as high as 7.7 g/cm^2 were probed in materials ranging from plastic to tungsten, and radiographs with resolution as good as 90 micrometers were produced. The effects of electric fields produced by the laser ablation of the radiography objects were observed and are well described by an analytic expression relating imaging magnification change to electric field strength. △ Less

Submitted 18 January, 2023; v1 submitted 29 September, 2022; originally announced September 2022.

Split Hamiltonian Monte Carlo revisited

Authors: Fernando Casas, Jesús María Sanz-Serna, Luke Shaw

Abstract: We study Hamiltonian Monte Carlo (HMC) samplers based on splitting the Hamiltonian $H$ as $H_0(θ,p)+U_1(θ)$, where $H_0$ is quadratic and $U_1$ small. We show that, in general, such samplers suffer from stepsize stability restrictions similar to those of algorithms based on the standard leapfrog integrator. The restrictions may be circumvented by preconditioning the dynamics. Numerical experiments… ▽ More We study Hamiltonian Monte Carlo (HMC) samplers based on splitting the Hamiltonian $H$ as $H_0(θ,p)+U_1(θ)$, where $H_0$ is quadratic and $U_1$ small. We show that, in general, such samplers suffer from stepsize stability restrictions similar to those of algorithms based on the standard leapfrog integrator. The restrictions may be circumvented by preconditioning the dynamics. Numerical experiments show that, when the $H_0(θ,p)+U_1(θ)$ splitting is combined with preconditioning, it is possible to construct samplers far more efficient than standard leapfrog HMC. △ Less

Submitted 15 July, 2022; originally announced July 2022.

Comments: 25 pages, 6 figures

Benchmarking Quantum Simulators using Ergodic Quantum Dynamics

Authors: Daniel K. Mark, Joonhee Choi, Adam L. Shaw, Manuel Endres, Soonwon Choi

Abstract: We propose and analyze a sample-efficient protocol to estimate the fidelity between an experimentally prepared state and an ideal target state, applicable to a wide class of analog quantum simulators without advanced sophisticated spatiotemporal control. Our approach utilizes newly discovered universal fluctuations emerging from generic Hamiltonian dynamics, and it does not require any fine-tuned… ▽ More We propose and analyze a sample-efficient protocol to estimate the fidelity between an experimentally prepared state and an ideal target state, applicable to a wide class of analog quantum simulators without advanced sophisticated spatiotemporal control. Our approach utilizes newly discovered universal fluctuations emerging from generic Hamiltonian dynamics, and it does not require any fine-tuned control over state preparation, quantum evolution, or readout capability. It only needs a small number of experimental measurements, achieving near optimal sample complexity: in ideal cases, a percent-level precision is obtained with $\sim 10^3$ measurements independent of system size. Furthermore, the accuracy of our fidelity estimation improves with increasing system size. We numerically demonstrate our protocol for a variety of quantum simulator platforms such as itinerant particles on optical lattices, trapped ions, and Rydberg atoms. We discuss further applications of our method for advanced tasks such as multi-parameter estimation of quantum states and processes. △ Less

Submitted 16 August, 2023; v1 submitted 24 May, 2022; originally announced May 2022.

Comments: 7+27 pages, 3+9 figures v2: version accepted in Physical Review Letters

Report number: MIT-CTP 5438

Journal ref: Phys. Rev. Lett. 131, 110601 (2023)

The logic of planetary combination in Vettius Valens

Authors: Claire Hall, Liam P. Shaw

Abstract: The Anthologies of the second-century astrologer Vettius Valens (120-c.175 CE) is the most extensive surviving practical astrological text from the period. Despite this, the theoretical underpinnings of the Anthologies have been understudied; in general, the work has been overshadowed by Ptolemy's contemporaneous Tetrabiblos. While the Tetrabiblos explicitly aims to present a systematic account of… ▽ More The Anthologies of the second-century astrologer Vettius Valens (120-c.175 CE) is the most extensive surviving practical astrological text from the period. Despite this, the theoretical underpinnings of the Anthologies have been understudied; in general, the work has been overshadowed by Ptolemy's contemporaneous Tetrabiblos. While the Tetrabiblos explicitly aims to present a systematic account of astrology, Valens' work is often characterised as a miscellaneous collection, of interest to historians only for the evidence it preserves about the practical methods used in casting horoscopes. In this article, we argue that the Anthologies is also an invaluable resource for engagement with the conceptual basis of astrology. As a case study, we take a section of Anthologies Book 1 which lists the possible astrological effects of planets, both alone and in 'combinations' of two and three. We demonstrate that analysing Valens' descriptions quantitatively with textual analysis reveals a consistent internal logic of planetary combination. By classifying descriptive terms as positive or negative, we show that the resulting 'sentiment' of planetary combinations is well-correlated with their component parts. Furthermore, we find that the sentiment of three-planet combinations is more strongly correlated with the average sentiment of their three possible component pairs than with the average sentiment of individual planets, suggesting an iterative combinatorial logic. Recognition of this feature of astrological practice has been neglected compared to the mathematical methods for calculating horoscopes. We argue that this analysis not only provides evidence that the astrological lore detailed in Valens is more consistent than is often assumed, but is also indicative of a wider methodological technique in practical astrology: combinatorial reasoning from existing astrological lore. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 30 pages, 7 figures

Emergent quantum state designs from individual many-body wavefunctions

Authors: Jordan S. Cotler, Daniel K. Mark, Hsin-Yuan Huang, Felipe Hernandez, Joonhee Choi, Adam L. Shaw, Manuel Endres, Soonwon Choi

Abstract: Quantum chaos in many-body systems provides a bridge between statistical and quantum physics with strong predictive power. This framework is valuable for analyzing properties of complex quantum systems such as energy spectra and the dynamics of thermalization. While contemporary methods in quantum chaos often rely on random ensembles of quantum states and Hamiltonians, this is not reflective of mo… ▽ More Quantum chaos in many-body systems provides a bridge between statistical and quantum physics with strong predictive power. This framework is valuable for analyzing properties of complex quantum systems such as energy spectra and the dynamics of thermalization. While contemporary methods in quantum chaos often rely on random ensembles of quantum states and Hamiltonians, this is not reflective of most real-world systems. In this paper, we introduce a new perspective: across a wide range of examples, a single non-random quantum state is shown to encode universal and highly random quantum state ensembles. We characterize these ensembles using the notion of quantum state $k$-designs from quantum information theory and investigate their universality using a combination of analytic and numerical techniques. In particular, we establish that $k$-designs arise naturally from generic states as well as individual states associated with strongly interacting, time-independent Hamiltonian dynamics. Our results offer a new approach for studying quantum chaos and provide a practical method for sampling approximately uniformly random states; the latter has wide-ranging applications in quantum information science from tomography to benchmarking. △ Less

Submitted 5 March, 2021; originally announced March 2021.

Comments: 7+19 pages, 6 figures

Journal ref: PRX Quantum 4, 010311 (2023)

Preparing random states and benchmarking with many-body quantum chaos

Authors: Joonhee Choi, Adam L. Shaw, Ivaylo S. Madjarov, Xin Xie, Ran Finkelstein, Jacob P. Covey, Jordan S. Cotler, Daniel K. Mark, Hsin-Yuan Huang, Anant Kale, Hannes Pichler, Fernando G. S. L. Brandão, Soonwon Choi, Manuel Endres

Abstract: Producing quantum states at random has become increasingly important in modern quantum science, with applications both theoretical and practical. In particular, ensembles of such randomly-distributed, but pure, quantum states underly our understanding of complexity in quantum circuits and black holes, and have been used for benchmarking quantum devices in tests of quantum advantage. However, creat… ▽ More Producing quantum states at random has become increasingly important in modern quantum science, with applications both theoretical and practical. In particular, ensembles of such randomly-distributed, but pure, quantum states underly our understanding of complexity in quantum circuits and black holes, and have been used for benchmarking quantum devices in tests of quantum advantage. However, creating random ensembles has necessitated a high degree of spatio-temporal control, placing such studies out of reach for a wide class of quantum systems. Here we solve this problem by predicting and experimentally observing the emergence of random state ensembles naturally under time-independent Hamiltonian dynamics, which we use to implement an efficient, widely applicable benchmarking protocol. The observed random ensembles emerge from projective measurements and are intimately linked to universal correlations built up between subsystems of a larger quantum system, offering new insights into quantum thermalization. Predicated on this discovery, we develop a fidelity estimation scheme, which we demonstrate for a Rydberg quantum simulator with up to 25 atoms using fewer than 10^4 experimental samples. This method has broad applicability, as we show for Hamiltonian parameter estimation, target-state generation benchmarking, and comparison of analog and digital quantum devices. Our work has implications for understanding randomness in quantum dynamics, and enables applications of this concept in a much wider context. △ Less

Submitted 16 May, 2023; v1 submitted 5 March, 2021; originally announced March 2021.

Comments: JC and ALS contributed equally to this work

Journal ref: Nature 613, 468 (2023)

Predominant Contribution of Direct Laser Acceleration to High-Energy Electron Spectra in a Low-Density Self-Modulated Laser Wakefield Accelerator

Authors: P. M. King, K. Miller, N. Lemos, J. L. Shaw, B. F. Kraus, M. Thibodeau, B. M. Hegelich, J. Hinojosa, P. Michel, C. Joshi, K. A. Marsh, W. Mori, A. Pak, A. G. R. Thomas, F. Albert

Abstract: The two-temperature relativistic electron spectrum from a low-density ($3\times10^{17}$~cm$^{-3}$) self-modulated laser wakefield accelerator (SM-LWFA) is observed to transition between temperatures of $19\pm0.65$ and $46\pm2.45$ MeV at an electron energy of about 100 MeV. When the electrons are dispersed orthogonally to the laser polarization, their spectrum above 60 MeV shows a forking structure… ▽ More The two-temperature relativistic electron spectrum from a low-density ($3\times10^{17}$~cm$^{-3}$) self-modulated laser wakefield accelerator (SM-LWFA) is observed to transition between temperatures of $19\pm0.65$ and $46\pm2.45$ MeV at an electron energy of about 100 MeV. When the electrons are dispersed orthogonally to the laser polarization, their spectrum above 60 MeV shows a forking structure characteristic of direct laser acceleration (DLA). Both the two-temperature distribution and the forking structure are reproduced in a quasi-3D \textsc{Osiris} simulation of the interaction of the 1-ps, moderate-amplitude ($a_{0}=2.7$) laser pulse with the low-density plasma. Particle tracking shows that while the SM-LWFA mechanism dominates below 40 MeV, the highest-energy ($>60$ MeV) electrons gain most of their energy through DLA. By separating the simulated electric fields into modes, the DLA-dominated electrons are shown to lose significant energy to the longitudinal laser field from the tight focusing geometry, resulting in a more accurate measure of net DLA energy gain than previously possible. △ Less

Submitted 2 December, 2020; originally announced December 2020.

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

Structural study of TATB under detonation-induced shock conditions

Authors: Elissaios Stavrou, Michael Bagge-Hansen, Joshua A. Hammons, Michael H. Nielsen, William L. Shaw, Will Bassett, Thomas W. Myers, Lisa M. Lauderbach, Ralph L. Hodgin, Nicholas A. Perez-Marty, Kamel Fezzaa, Alex Deriy, Sorin Bastea, Laurence E. Fried, Lara D. Leininger, Trevor M. Willey

Abstract: We explore the response of the insensitive high explosive (IHE) 1,3,5-Triamino-2,4,6- trinitrobenzene (TATB) under detonation-induced shock conditions using in-situ synchrotron X-ray diffraction in the 100 ns time scale using either a conventional or a colliding detonation drive. In all of the detonation experiments on various sizes and morphologies of TATB, we observe an extended stability of the… ▽ More We explore the response of the insensitive high explosive (IHE) 1,3,5-Triamino-2,4,6- trinitrobenzene (TATB) under detonation-induced shock conditions using in-situ synchrotron X-ray diffraction in the 100 ns time scale using either a conventional or a colliding detonation drive. In all of the detonation experiments on various sizes and morphologies of TATB, we observe an extended stability of the TATB triclinic crystal structure. As the detonation front passes through the TATB, X-ray diffraction indicates a portion of the TATB exhibits a compression up to 30+ GPa, followed subsequently by a pressure release and continued decomposition over a few hundred nanoseconds. Likewise, for colliding detonation-driven shock compression of single crystals of TATB, a significant portion of the triclinic crystal structure appears to be stable up to 60+ GPa. Conversely, in similar detonations of a LLM-105 PBX, X-ray diffraction is simply indicative of decomposition without the apparent compression and slow decomposition seen in TATB. The results indicate the surprising resilience of TATB under these high-pressure, temperature and shock conditions, providing a baseline for understanding the insensitivity of TATB that is considered the industry standard for an insensitive high explosive. The results also provide intriguing information for the extended reaction zone in TATB, and the hot-spot mechanisms for initiating and propagating detonation in this uniquely insensitive explosive. △ Less

Submitted 17 September, 2021; v1 submitted 22 November, 2020; originally announced November 2020.

Comments: Accepted in Physical Review Materials, https://journals.aps.org/prmaterials/accepted/a6076KafM7c1c00414b62063da29797895775b4b6

Impacts of density-dependent predation, cannibalism and fishing in a stage-structured population model of the blue crab in Chesapeake Bay

Authors: Fangming Xu, Leah B. Shaw, Junping Shi, Romuald N. Lipcius

Abstract: The blue crab (Callinectes sapidus) is a dominant ecological species of high commercial value. Spawning stock and recruitment of the Chesapeake Bay population declined by 80% in the 1990s. After severe management actions were implemented in 2008, female abundance rebounded to pre-1994 levels and stabilized. The stepwise decline in the early 1990s, followed by a consistently low level of abundance… ▽ More The blue crab (Callinectes sapidus) is a dominant ecological species of high commercial value. Spawning stock and recruitment of the Chesapeake Bay population declined by 80% in the 1990s. After severe management actions were implemented in 2008, female abundance rebounded to pre-1994 levels and stabilized. The stepwise decline in the early 1990s, followed by a consistently low level of abundance for 15 y and a jump to high abundance after 2008, suggested the existence of alternative stable states. Alternatively, high fishing pressure combined with low recruitment in 1992 could have triggered a proportional decline in the population, followed by a population increase in 2008 due to rigorous management actions that reduced fishing. We evaluated these alternatives with a stage-structured dynamic population model using ordinary differential equations. In addition, stock assessment models assume that fishing and mortality are independent of density. Hence, we also investigated the role of density-dependent predation, cannibalism and fishing in blue crab population dynamics. We conclude that for the blue crab population in Chesapeake Bay: (1) bistable positive states are not likely with biologically realistic parameter values; (2) hyperbolic (depensatory) fishing will not produce extinction at the range of population densities observed in the bay; and (3) crabs can survive a higher fishing rate under the more realistic assumption of sigmoidal (density-dependent) predation and cannibalism than under constant (density-independent) predation and cannibalism. These collectively indicate that the blue crab population in Chesapeake Bay is resilient to a range of biotic and abiotic disturbances. △ Less

Submitted 16 November, 2020; originally announced November 2020.

Generating ultra-dense pair beams using 400 GeV/c protons

Authors: C. D. Arrowsmith, N. Shukla, N. Charitonidis, R. Boni, H. Chen, T. Davenne, D. H. Froula, B. T. Huffman, Y. Kadi, B. Reville, S. Richardson, S. Sarkar, J. L. Shaw, L. O. Silva, R. M. G. M. Trines, R. Bingham, G. Gregori

Abstract: A previously unexplored experimental scheme is presented for generating low-divergence, ultra-dense, relativistic, electron-positron beams using 400 GeV/c protons available at facilities such as HiRadMat and AWAKE at CERN. Preliminary Monte-Carlo and Particle-in-cell simulations demonstrate the possibility of generating beams containing $10^{13}-10^{14}$ electron-positron pairs at sufficiently hig… ▽ More A previously unexplored experimental scheme is presented for generating low-divergence, ultra-dense, relativistic, electron-positron beams using 400 GeV/c protons available at facilities such as HiRadMat and AWAKE at CERN. Preliminary Monte-Carlo and Particle-in-cell simulations demonstrate the possibility of generating beams containing $10^{13}-10^{14}$ electron-positron pairs at sufficiently high densities to drive collisionless beam-plasma instabilities, which are expected to play an important role in magnetic field generation and the related radiation signatures of relativistic astrophysical phenomena. The pair beams are quasi-neutral, with size exceeding several skin-depths in all dimensions, allowing for the first time the examination of the effect of competition between transverse and longitudinal instability modes on the growth of magnetic fields. Furthermore, the presented scheme allows for the possibility of controlling the relative density of hadrons to electron-positron pairs in the beam, making it possible to explore the parameter spaces for different astrophysical environments. △ Less

Submitted 9 November, 2020; originally announced November 2020.

Journal ref: Phys. Rev. Research 3, 023103 (2021)

Laser-plasma acceleration beyond wave breaking

Authors: J. P. Palastro, B. Malaca, J. Vieira, D. Ramsey, T. T. Simpson, P. Franke, J. L. Shaw, D. H. Froula

Abstract: Laser wakefield accelerators rely on the extremely high electric fields of nonlinear plasma waves to trap and accelerate electrons to relativistic energies over short distances. When driven strongly enough, plasma waves break, trapping a large population of the background electrons that support their motion. This limits the maximum electric field. Here we introduce a novel regime of plasma wave ex… ▽ More Laser wakefield accelerators rely on the extremely high electric fields of nonlinear plasma waves to trap and accelerate electrons to relativistic energies over short distances. When driven strongly enough, plasma waves break, trapping a large population of the background electrons that support their motion. This limits the maximum electric field. Here we introduce a novel regime of plasma wave excitation and wakefield acceleration that removes this limit, allowing for arbitrarily high electric fields. The regime, enabled by spatiotemporal shaping of laser pulses, exploits the property that nonlinear plasma waves with superluminal phase velocities cannot trap charged particles and are therefore immune to wave breaking. A laser wakefield accelerator operating in this regime provides energy tunability independent of the plasma density and can accommodate the large laser amplitudes delivered by modern and planned high-power, short pulse laser systems. △ Less

Submitted 20 August, 2020; originally announced August 2020.

Detonation-induced transformation of graphite to hexagonal diamond

Authors: Elissaios Stavrou, Michael Bagge-Hansen, Joshua A. Hammons, Michael H. Nielsen, Bradley A. Steele, Penghao Xiao, Matthew P. Kroonblawd, Matthew D. Nelms, William L. Shaw, Will Bassett, Sorin Bastea, Lisa M. Lauderbach, Ralph L. Hodgin, Nicholas A. Perez-Marty, Saransh Singh, Pinaki Das, Yuelin Li, Adam Schuman, Nicholas Sinclair, Kamel Fezzaa, Alex Deriy, Lara D. Leininger, Trevor M. Willey

Abstract: We explore the structural evolution of highly oriented pyrolytic graphite (HOPG) under detonation-induced shock conditions using in-situ synchrotron X-ray diffraction in the ns time scale. We observe the formation of hexagonal diamond (lonsdaleite) at pressures above 50 GPa, in qualitative agreement with recent gas gun experiments. First-principles density functional calculations reveal that under… ▽ More We explore the structural evolution of highly oriented pyrolytic graphite (HOPG) under detonation-induced shock conditions using in-situ synchrotron X-ray diffraction in the ns time scale. We observe the formation of hexagonal diamond (lonsdaleite) at pressures above 50 GPa, in qualitative agreement with recent gas gun experiments. First-principles density functional calculations reveal that under uniaxial compression the energy barrier for the transition towards hexagonal diamond is lower than cubic diamond. Finally, no indication of cubic diamond formation was observed up to >70 GPa. △ Less

Submitted 22 July, 2020; originally announced July 2020.

Journal ref: Phys. Rev. B 102, 104116 (2020)

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

Authors: Clara L. Shaw, David A. Kennedy

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

Submitted 25 June, 2020; originally announced June 2020.

Comments: 25 pages, 2 figures

Efficient extraction of high pulse energy from partly quenched highly Er3+-doped fiber amplifiers

Authors: Pablo G. Rojas Hernandez, Mohammad Belal, Colin Baker, Shankar Pidishety, Yutong Feng, E. Joseph Friebele, L. Brandon Shaw, Daniel Rhonehouse, Jasbinder Sanghera, Johan Nilsson

Abstract: We demonstrate efficient pulse-energy extraction from a partly-quenched erbium-doped aluminosilicate fiber amplifier. This has a high erbium-concentration, which allows for short devices with reduced nonlinear distortions, but which also results in partial quenching and thus significant unsaturable absorption, even though the fiber is still able to amplify. Although the quenching degrades the aver… ▽ More We demonstrate efficient pulse-energy extraction from a partly-quenched erbium-doped aluminosilicate fiber amplifier. This has a high erbium-concentration, which allows for short devices with reduced nonlinear distortions, but which also results in partial quenching and thus significant unsaturable absorption, even though the fiber is still able to amplify. Although the quenching degrades the average-power efficiency, the pulse energy remains high, and our results point to an increasingly promising outcome for short pulses. Furthermore, unlike unquenched fibers, the conversion efficiency improves at low repetition rates, which we attribute to smaller relative energy loss to quenched ions at higher pulse energy. A short (2.6 m) cladding-pumped partly-quenched Er-doped-fiber with 95-dB/m 1530-nm peak absorption and saturation energy estimated to 85 μJ, reached 0.8 mJ of output energy when seeded by 0.2-μs, 23-μJ pulses. Thus, according to our results, pulses can be amplified to high energy in short highly-Er-doped fibers designed to reduce nonlinear distortions, at the expense of average-power efficiency △ Less

Submitted 11 February, 2020; originally announced February 2020.

Comments: Manuscript submitted to Optics Express. Peer-Review extended up to March, 31

High-Fidelity Entanglement and Detection of Alkaline-Earth Rydberg Atoms

Authors: Ivaylo S. Madjarov, Jacob P. Covey, Adam L. Shaw, Joonhee Choi, Anant Kale, Alexandre Cooper, Hannes Pichler, Vladimir Schkolnik, Jason R. Williams, Manuel Endres

Abstract: Trapped neutral atoms have become a prominent platform for quantum science, where entanglement fidelity records have been set using highly-excited Rydberg states. However, controlled two-qubit entanglement generation has so far been limited to alkali species, leaving the exploitation of more complex electronic structures as an open frontier that could lead to improved fidelities and fundamentally… ▽ More Trapped neutral atoms have become a prominent platform for quantum science, where entanglement fidelity records have been set using highly-excited Rydberg states. However, controlled two-qubit entanglement generation has so far been limited to alkali species, leaving the exploitation of more complex electronic structures as an open frontier that could lead to improved fidelities and fundamentally different applications such as quantum-enhanced optical clocks. Here we demonstrate a novel approach utilizing the two-valence electron structure of individual alkaline-earth Rydberg atoms. We find fidelities for Rydberg state detection, single-atom Rabi operations, and two-atom entanglement surpassing previously published values. Our results pave the way for novel applications, including programmable quantum metrology and hybrid atom-ion systems, and set the stage for alkaline-earth based quantum computing architectures. △ Less

Submitted 28 November, 2020; v1 submitted 13 January, 2020; originally announced January 2020.

Comments: ISM and JPC contributed equally to this work

Journal ref: Nature Physics 16, 857-861 (2020)

Deep learning within a priori temporal feature spaces for large-scale dynamic MR image reconstruction: Application to 5-D cardiac MR Multitasking

Authors: Yuhua Chen, Jaime L. Shaw, Yibin Xie, Debiao Li, Anthony G. Christodoulou

Abstract: High spatiotemporal resolution dynamic magnetic resonance imaging (MRI) is a powerful clinical tool for imaging moving structures as well as to reveal and quantify other physical and physiological dynamics. The low speed of MRI necessitates acceleration methods such as deep learning reconstruction from under-sampled data. However, the massive size of many dynamic MRI problems prevents deep learnin… ▽ More High spatiotemporal resolution dynamic magnetic resonance imaging (MRI) is a powerful clinical tool for imaging moving structures as well as to reveal and quantify other physical and physiological dynamics. The low speed of MRI necessitates acceleration methods such as deep learning reconstruction from under-sampled data. However, the massive size of many dynamic MRI problems prevents deep learning networks from directly exploiting global temporal relationships. In this work, we show that by applying deep neural networks inside a priori calculated temporal feature spaces, we enable deep learning reconstruction with global temporal modeling even for image sequences with >40,000 frames. One proposed variation of our approach using dilated multi-level Densely Connected Network (mDCN) speeds up feature space coordinate calculation by 3000x compared to conventional iterative methods, from 20 minutes to 0.39 seconds. Thus, the combination of low-rank tensor and deep learning models not only makes large-scale dynamic MRI feasible but also practical for routine clinical application. △ Less

Submitted 2 October, 2019; originally announced October 2019.

Comments: Early accepted by MICCAI 2019

Journal ref: Medical Image Computing and Computer Assisted Intervention 2019 pp 495-504

An atomic array optical clock with single-atom readout

Authors: Ivaylo S. Madjarov, Alexandre Cooper, Adam L. Shaw, Jacob P. Covey, Vladimir Schkolnik, Tai Hyun Yoon, Jason R. Williams, Manuel Endres

Abstract: Currently, the most accurate and stable clocks use optical interrogation of either a single ion or an ensemble of neutral atoms confined in an optical lattice. Here, we demonstrate a new optical clock system based on an array of individually trapped neutral atoms with single-atom readout, merging many of the benefits of ion and lattice clocks as well as creating a bridge to recently developed tech… ▽ More Currently, the most accurate and stable clocks use optical interrogation of either a single ion or an ensemble of neutral atoms confined in an optical lattice. Here, we demonstrate a new optical clock system based on an array of individually trapped neutral atoms with single-atom readout, merging many of the benefits of ion and lattice clocks as well as creating a bridge to recently developed techniques in quantum simulation and computing with neutral atoms. We evaluate single-site resolved frequency shifts and short-term stability via self-comparison. Atom-by-atom feedback control enables direct experimental estimation of laser noise contributions. Results agree well with an ab initio Monte Carlo simulation that incorporates finite temperature, projective read-out, laser noise, and feedback dynamics. Our approach, based on a tweezer array, also suppresses interaction shifts while retaining a short dead time, all in a comparatively simple experimental setup suited for transportable operation. These results establish the foundations for a third optical clock platform and provide a novel starting point for entanglement-enhanced metrology, quantum clock networks, and applications in quantum computing and communication with individual neutral atoms that require optical clock state control. △ Less

Submitted 29 October, 2019; v1 submitted 15 August, 2019; originally announced August 2019.

Comments: 14 pages, 8 figures, 1 table; accepted in PRX on October 25th, 2019

Journal ref: Phys. Rev. X 9, 041052 (2019)

Model of pattern formation in marsh ecosystems with nonlocal interactions

Authors: Sofya Zaytseva, Junping Shi, Leah B Shaw

Abstract: Smooth cordgrass Spartina alterniflora is a grass species commonly found in tidal marshes. It is an ecosystem engineer, capable of modifying the structure of its surrounding environment through various feedbacks. The scale-dependent feedback between marsh grass and sediment volume is particularly of interest. Locally, the marsh vegetation attenuates hydrodynamic energy, enhancing sediment accretio… ▽ More Smooth cordgrass Spartina alterniflora is a grass species commonly found in tidal marshes. It is an ecosystem engineer, capable of modifying the structure of its surrounding environment through various feedbacks. The scale-dependent feedback between marsh grass and sediment volume is particularly of interest. Locally, the marsh vegetation attenuates hydrodynamic energy, enhancing sediment accretion and promoting further vegetation growth. In turn, the diverted water flow promotes the formation of erosion troughs over longer distances. This scale-dependent feedback may explain the characteristic spatially varying marsh shoreline, commonly observed in nature. We propose a mathematical framework to model grass-sediment dynamics as a system of reaction-diffusion equations with an additional nonlocal term quantifying the short-range positive and long-range negative grass-sediment interactions. We use a Mexican-hat kernel function to model this scale-dependent feedback. We perform a steady state biharmonic approximation of our system and derive conditions for the emergence of spatial patterns, corresponding to a spatially varying marsh shoreline. We find that the emergence of such patterns depends on the spatial scale and strength of the scale-dependent feedback, specified by the width and amplitude of the Mexican-hat kernel function. △ Less

Submitted 3 July, 2019; originally announced July 2019.

Divergence-Free Magnetohydrodynamics on Conformally Moving, Adaptive Meshes Using a Vector Potential Method

Authors: P. Chris Fragile, Daniel Nemergut, Payden L. Shaw, Peter Anninos

Abstract: We present a new method for evolving the equations of magnetohydrodynamics (both Newtonian and relativistic) that is capable of maintaining a divergence-free magnetic field ($\nabla \cdot \mathbf{B} = 0$) on adaptively refined, conformally moving meshes. The method relies on evolving the magnetic vector potential and then using it to reconstruct the magnetic fields. The advantage of this approach… ▽ More We present a new method for evolving the equations of magnetohydrodynamics (both Newtonian and relativistic) that is capable of maintaining a divergence-free magnetic field ($\nabla \cdot \mathbf{B} = 0$) on adaptively refined, conformally moving meshes. The method relies on evolving the magnetic vector potential and then using it to reconstruct the magnetic fields. The advantage of this approach is that the vector potential is not subject to a constraint equation in the same way the magnetic field is, and so can be refined and moved in a straightforward way. We test this new method against a wide array of problems from simple Alfven waves on a uniform grid to general relativistic MHD simulations of black hole accretion on a nested, spherical-polar grid. We find that the code produces accurate results and in all cases maintains a divergence-free magnetic field to machine precision. △ Less

Submitted 4 December, 2018; originally announced December 2018.

Comments: 25 pages, 16 figures, submitted to Journal of Computational Physics

Journal ref: Journal of Computational Physics: X 2 (2019) 100020

Intervention Strategies for Epidemics: Does Ignoring Time Delay Lead to Incorrect Predictions?

Authors: Adrienna Bingham, Leah B. Shaw

Abstract: Our paper investigates distributions of exposed and infectious time periods in an epidemic model and how applying a disease control strategy affects the model's accuracy. While ordinary differential equations are widely used for their simplicity, they incorporate an exponential distribution for time spent exposed or infectious. This allows for a high probability of unrealistically short exposed an… ▽ More Our paper investigates distributions of exposed and infectious time periods in an epidemic model and how applying a disease control strategy affects the model's accuracy. While ordinary differential equations are widely used for their simplicity, they incorporate an exponential distribution for time spent exposed or infectious. This allows for a high probability of unrealistically short exposed and infectious time periods. We propose that caution must be taken when applying intervention methods to basic models in order to avoid inaccurate predictions. Delay differential equations, which use a delta distribution for exposed and infectious periods, can provide better realism but are more difficult to use and analyze. We introduce a multi-infected compartment model to interpolate between an ODE model with exponential distributions and a DDE model with delta distributions in order to investigate the effect these distributions have on the dynamics of the system when an intervention method is also included. Using steady state stability and bifurcation analysis, this paper considers when simpler infectious disease models can be used versus when more realistic time periods must be incorporated. We find that the placement of control measures on subpopulations and the length of the time delay impacts the accuracy of the simpler models. △ Less

Submitted 27 September, 2018; originally announced September 2018.

MSC Class: 92D25; 92D30

Ionization waves of arbitrary velocity driven by a flying focus

Authors: J. P. Palastro, D. Turnbull, S. -W. Bahk, R. K. Follett, J. L. Shaw, D. Haberberger, J. Bromage, D. H. Froula

Abstract: A chirped laser pulse focused by a chromatic lens exhibits a dynamic, or "flying," focus in which the trajectory of the peak intensity decouples from the group velocity. In a medium, the flying focus can trigger an ionization front that follows this trajectory. By adjusting the chirp, the ionization front can be made to travel at an arbitrary velocity along the optical axis. We present analytical… ▽ More A chirped laser pulse focused by a chromatic lens exhibits a dynamic, or "flying," focus in which the trajectory of the peak intensity decouples from the group velocity. In a medium, the flying focus can trigger an ionization front that follows this trajectory. By adjusting the chirp, the ionization front can be made to travel at an arbitrary velocity along the optical axis. We present analytical calculations and simulations describing the propagation of the flying focus pulse, the self-similar form of its intensity profile, and ionization wave formation. The ability to control the speed of the ionization wave and, in conjunction, mitigate plasma refraction has the potential to advance several laser-based applications, including Raman amplification, photon acceleration, high harmonic generation, and THz generation. △ Less

Submitted 20 December, 2017; originally announced December 2017.

Journal ref: Phys. Rev. A 97, 033835 (2018)

Enhancement of large fluctuations to extinction in adaptive networks

Authors: Jason Hindes, Ira B. Schwartz, Leah B. Shaw

Abstract: During an epidemic, individual nodes in a network may adapt their connections to reduce the chance of infection. A common form of adaption is avoidance rewiring, where a noninfected node breaks a connection to an infected neighbor and forms a new connection to another noninfected node. Here we explore the effects of such adaptivity on stochastic fluctuations in the susceptible-infected-susceptible… ▽ More During an epidemic, individual nodes in a network may adapt their connections to reduce the chance of infection. A common form of adaption is avoidance rewiring, where a noninfected node breaks a connection to an infected neighbor and forms a new connection to another noninfected node. Here we explore the effects of such adaptivity on stochastic fluctuations in the susceptible-infected-susceptible model, focusing on the largest fluctuations that result in extinction of infection. Using techniques from large-deviation theory, combined with a measurement of heterogeneity in the susceptible degree distribution at the endemic state, we are able to predict and analyze large fluctuations and extinction in adaptive networks. We find that in the limit of small rewiring there is a sharp exponential reduction in mean extinction times compared to the case of zero adaption. Furthermore, we find an exponential enhancement in the probability of large fluctuations with increased rewiring rate, even when holding the average number of infected nodes constant. △ Less

Submitted 16 January, 2018; v1 submitted 10 August, 2017; originally announced August 2017.

Journal ref: Phys. Rev. E 97, 012308 (2018)

Self-modulated laser wakefield accelerators as x-ray sources

Authors: N. Lemos, J. L. Martins, F. S. Tsung, J. L. Shaw, K. A. Marsh, F. Albert, B. B. Pollock, C. Joshi

Abstract: The development of a directional, small-divergence, and short-duration picosecond x-ray probe beam with an energy greater than 50 keV is desirable for high energy density science experiments. We therefore explore through particle-in-cell (PIC) computer simulations the possibility of using x-rays radiated by betatron-like motion of electrons from a self-modulated laser wakefield accelerator as a po… ▽ More The development of a directional, small-divergence, and short-duration picosecond x-ray probe beam with an energy greater than 50 keV is desirable for high energy density science experiments. We therefore explore through particle-in-cell (PIC) computer simulations the possibility of using x-rays radiated by betatron-like motion of electrons from a self-modulated laser wakefield accelerator as a possible candidate to meet this need. Two OSIRIS 2D PIC simulations with mobile ions are presented, one with a normalized vector potential a0 = 1.5 and the other with an a0 = 3. We find that in both cases direct laser acceleration (DLA) is an important additional acceleration mechanism in addition to the longitudinal electric field of the plasma wave. Together these mechanisms produce electrons with a continuous energy spectrum with a maximum energy of 300 MeV for a0 = 3 case and 180 MeV in the a0 = 1.5 case. Forward-directed x-ray radiation with a photon energy up to 100 keV was calculated for the a0 = 3 case and up to 12 keV for the a0 = 1.5 case. The x-ray spectrum can be fitted with a sum of two synchrotron spectra with critical photon energy of 13 and 45 keV for the a0 of 3 and critical photon energy of 0.3 and 1.4 keV for a0 of 1.5 in the plane of polarization of the laser. The full width at half maximum divergence angle of the x-rays was 62 x 1.9 mrad for a0 = 3 and 77 x 3.8 mrad for a0 = 1.5. △ Less

Submitted 5 October, 2015; originally announced October 2015.

Estimation of direct laser acceleration in laser wakefield accelerators using particle-in-cell simulations

Authors: J. L. Shaw, N. Lemos, K. A. Marsh, F. S. Tsung, W. B. Mori, C. Joshi

Abstract: Many current laser wakefield acceleration (LWFA) experiments are carried out in a regime where the laser pulse length is on the order of or longer than the wake wavelength and where ionization injection is employed to inject electrons into the wake. In these experiments, the trapped electrons will co-propagate with the longitudinal wakefield and the transverse laser field. In this scenario, the el… ▽ More Many current laser wakefield acceleration (LWFA) experiments are carried out in a regime where the laser pulse length is on the order of or longer than the wake wavelength and where ionization injection is employed to inject electrons into the wake. In these experiments, the trapped electrons will co-propagate with the longitudinal wakefield and the transverse laser field. In this scenario, the electrons can gain a significant amount of energy from both the direct laser acceleration (DLA) mechanism as well as the usual LWFA mechanism. Particle-in-cell (PIC) codes are frequently used to discern the relative contribution of these two mechanisms. However, if the longitudinal resolution used in the PIC simulations is inadequate, it can produce numerical heating that can overestimate the transverse motion, which is important in determining the energy gain due to DLA. We have therefore carried out a systematic study of this LWFA regime by varying the longitudinal resolution of PIC simulations from the standard, best-practice resolution of 30 points per laser wavelength to four times that value and then examining the energy gain characteristics of both the highest-energy electrons and the bulk electrons. By calculating the contribution of DLA to the final energies of the electrons produced from the LWFA, we find that although the transverse momentum and oscillation radii are over-estimated in the lower-resolution simulations, this over-estimation does not lead to artificial energy gain by DLA. Rather, the DLA contribution to the highest-energy electrons is larger in the higher-resolution cases because the DLA resonance is better maintained. Thus, even at the highest longitudinal resolutions, DLA contributes a significant portion of the energy gained by the highest-energy electrons and also contributes to accelerating the bulk of the charge in the electron beam produced by the LWFA. △ Less

Submitted 25 September, 2015; originally announced September 2015.

Comments: 12 pages. 5 figures

Sub-femtosecond electron bunches created by direct laser acceleration in a laser wakefield accelerator with ionization injection

Authors: N. Lemos, J. L. Shaw, K. A. Marsh, C. Joshi

Abstract: In this work, we will show through three-dimensional particle-in-cell simulations that direct laser acceleration in laser a wakefield accelerator can generate sub-femtosecond electron bunches. Two simulations were done with two laser pulse durations, such that the shortest laser pulse occupies only a fraction of the first bubble, whereas the longer pulse fills the entire first bubble. In the latte… ▽ More In this work, we will show through three-dimensional particle-in-cell simulations that direct laser acceleration in laser a wakefield accelerator can generate sub-femtosecond electron bunches. Two simulations were done with two laser pulse durations, such that the shortest laser pulse occupies only a fraction of the first bubble, whereas the longer pulse fills the entire first bubble. In the latter case, as the trapped electrons moved forward and interacted with the high intensity region of the laser pulse, micro-bunching occurred naturally, producing 0.5 fs electron bunches. This is not observed in the short pulse simulation. △ Less

Submitted 26 February, 2015; originally announced February 2015.

Comments: AAC 2014

Satisfying the Direct Laser Acceleration Resonance Condition in a Laser Wakefield Accelerator

Authors: J. L. Shaw, N. Vafaei-Najafabadi, K. A. Marsh, N. Lemos, F. S. Tsung, W. B Mori, C. Joshi

Abstract: In this proceeding, we show that when the drive laser pulse overlaps the trapped electrons in a laser wakefield accelerator (LWFA), those electrons can gain energy from direct laser acceleration (DLA) over extended distances despite the evolution of both the laser and the wake. Through simulations, the evolution of the properties of both the laser and the electron beam is quantified, and then the… ▽ More In this proceeding, we show that when the drive laser pulse overlaps the trapped electrons in a laser wakefield accelerator (LWFA), those electrons can gain energy from direct laser acceleration (DLA) over extended distances despite the evolution of both the laser and the wake. Through simulations, the evolution of the properties of both the laser and the electron beam is quantified, and then the resonance condition for DLA is examined in the context of this change. We find that although the electrons produced from the LWFA cannot continuously satisfy the DLA resonance condition, they nevertheless can gain a significant amount of energy from DLA. △ Less

Submitted 26 February, 2015; originally announced February 2015.

Comments: 2014 Advanced Accelerator Concepts Workshop

Rare Event Extinction on Stochastic Networks

Authors: Brandon S. Lindley, Leah B. Shaw, Ira B. Schwartz

Abstract: We consider the problem of extinction processes on random networks with a given structure. For sufficiently large well-mixed populations, the process of extinction of one or more state variable components occurs in the tail of the quasi-stationary probability distribution, thereby making it a rare event. Here we show how to extend the theory of large deviations to random networks to predict extinc… ▽ More We consider the problem of extinction processes on random networks with a given structure. For sufficiently large well-mixed populations, the process of extinction of one or more state variable components occurs in the tail of the quasi-stationary probability distribution, thereby making it a rare event. Here we show how to extend the theory of large deviations to random networks to predict extinction times. In particular, we use the theory to find the most probable path leading to extinction. We apply the methodology to epidemic models and discover how mean extinction times scale with epidemiological and network parameters in Erdos-Renyi networks. The results are shown to compare quite well with Monte Carlo simulations of the network in predicting both the most probable paths to extinction and mean extinction times. △ Less

Submitted 8 December, 2014; v1 submitted 30 October, 2014; originally announced November 2014.

Comments: 5 pages 6 figures

Journal ref: EPL (Europhysics Letters) Volume 108 Number 5 58008 (2014)

Inferring unobserved multistrain epidemic sub-populations using synchronization dynamics

Authors: Eric Forgoston, Leah B. Shaw, Ira B. Schwartz

Abstract: A new method is proposed to infer unobserved epidemic sub-populations by exploiting the synchronization properties of multistrain epidemic models. A model for dengue fever is driven by simulated data from secondary infective populations. Primary infective populations in the driven system synchronize to the correct values from the driver system. Most hospital cases of dengue are secondary infection… ▽ More A new method is proposed to infer unobserved epidemic sub-populations by exploiting the synchronization properties of multistrain epidemic models. A model for dengue fever is driven by simulated data from secondary infective populations. Primary infective populations in the driven system synchronize to the correct values from the driver system. Most hospital cases of dengue are secondary infections, so this method provides a way to deduce unobserved primary infection levels. We derive center manifold equations that relate the driven system to the driver system and thus motivate the use of synchronization to predict unobserved primary infectives. Synchronization stability between primary and secondary infections is demonstrated through numerical measurements of conditional Lyapunov exponents and through time series simulations. △ Less

Submitted 30 October, 2014; originally announced October 2014.

Comments: 20 pages, 5 figures. arXiv admin note: text overlap with arXiv:1309.2600

DNA hairpins primarily promote duplex melting rather than inhibiting hybridization

Authors: John S. Schreck, Thomas E. Ouldridge, Flavio Romano, Petr Sulc, Liam Shaw, Ard A. Louis, Jonathan P. K. Doye

Abstract: The effect of secondary structure on DNA duplex formation is poorly understood. We use a coarse-grained model of DNA to show that specific 3- and 4-base pair hairpins reduce hybridization rates by factors of 2 and 10 respectively, in good agreement with experiment. By contrast, melting rates are accelerated by factors of ~100 and ~2000. This surprisingly large speed-up occurs because hairpins form… ▽ More The effect of secondary structure on DNA duplex formation is poorly understood. We use a coarse-grained model of DNA to show that specific 3- and 4-base pair hairpins reduce hybridization rates by factors of 2 and 10 respectively, in good agreement with experiment. By contrast, melting rates are accelerated by factors of ~100 and ~2000. This surprisingly large speed-up occurs because hairpins form during the melting process, stabilizing partially melted states, and facilitating dissociation. These results may help guide the design of DNA devices that use hairpins to modulate hybridization and dissociation pathways and rates. △ Less

Submitted 19 August, 2014; originally announced August 2014.

Comments: 5 pages + 14 pages of appendices

Journal ref: Nucl. Acids Res. 43, 6181-6190 (2015)

A measurement of secondary cosmic microwave background anisotropies from the 2500-square-degree SPT-SZ survey

Authors: E. M. George, C. L. Reichardt, K. A. Aird, B. A. Benson, L. E. Bleem, J. E. Carlstrom, C. L. Chang, H-M. Cho, T. M. Crawford, A. T. Crites, T. de Haan, M. A. Dobbs, J. Dudley, N. W. Halverson, N. L. Harrington, G. P. Holder, W. L. Holzapfel, Z. Hou, J. D. Hrubes, R. Keisler, L. Knox, A. T. Lee, E. M. Leitch, M. Lueker, D. Luong-Van , et al. (23 additional authors not shown)

Abstract: We present measurements of secondary cosmic microwave background (CMB) anisotropies and cosmic infrared background (CIB) fluctuations using data from the South Pole Telescope (SPT) covering the complete 2540 sq.deg. SPT-SZ survey area. Data in the three SPT-SZ frequency bands centered at 95, 150, and 220 GHz, are used to produce six angular power spectra (three single-frequency auto-spectra and th… ▽ More We present measurements of secondary cosmic microwave background (CMB) anisotropies and cosmic infrared background (CIB) fluctuations using data from the South Pole Telescope (SPT) covering the complete 2540 sq.deg. SPT-SZ survey area. Data in the three SPT-SZ frequency bands centered at 95, 150, and 220 GHz, are used to produce six angular power spectra (three single-frequency auto-spectra and three cross-spectra) covering the multipole range 2000 < ell < 11000 (angular scales 5' > θ> 1'). These are the most precise measurements of the angular power spectra at ell > 2500 at these frequencies. The main contributors to the power spectra at these angular scales and frequencies are the primary CMB, CIB, thermal and kinematic Sunyaev-Zel'dovich effects (tSZ and kSZ), and radio galaxies. We include a constraint on the tSZ power from a measurement of the tSZ bispectrum from 800 sq.deg. of the SPT-SZ survey. We measure the tSZ power at 143 GHz to be DtSZ = 4.08 +0.58 -0.67 μK^2 and the kSZ power to be DkSZ = 2.9 +- 1.3 μK^2. The data prefer positive kSZ power at 98.1% CL. We measure a correlation coefficient of ξ= 0.113 +0.057 -0.054 between sources of tSZ and CIB power, with ξ< 0 disfavored at a confidence level of 99.0%. The constraint on kSZ power can be interpreted as an upper limit on the duration of reionization. When the post-reionization homogeneous kSZ signal is accounted for, we find an upper limit on the duration Δz < 5.4 at 95% CL. △ Less

Submitted 1 December, 2014; v1 submitted 13 August, 2014; originally announced August 2014.

Comments: Accepted by ApJ (v2 is accepted version); 12 figures; 24 pages

Journal ref: 2015, ApJ, 799, 177

Protein Folding in the Hexagonal Prism Lattice with Diagonals

Authors: Dipan Lal Shaw, M. Sohel Rahman, A. S. M. Sohidull Islam, Shuvasish Karmaker

Abstract: Predicting protein secondary structure using lattice model is one of the most studied computational problem in bioinformatics. Here secondary structure or three dimensional structure of protein is predicted from its amino acid sequence. Secondary structure refers to local sub-structures of protein. Mostly founded secondary structures are alpha helix and beta sheets. Since, it is a problem of great… ▽ More Predicting protein secondary structure using lattice model is one of the most studied computational problem in bioinformatics. Here secondary structure or three dimensional structure of protein is predicted from its amino acid sequence. Secondary structure refers to local sub-structures of protein. Mostly founded secondary structures are alpha helix and beta sheets. Since, it is a problem of great potential complexity many simplified energy model have been proposed in literature on basis of interaction of amino acid residue in protein. Here we use well researched Hydrophobic-Polar (HP) energy model. In this paper, we proposed hexagonal prism lattice with diagonal that can overcome the problems of other lattice structure, e.g., parity problem. We give two approximation algorithm for protein folding on this lattice. Our first algorithm leads us to similar structure of helix structure which is commonly found in protein structure. This motivated us to find next algorithm which improves the algorithm ratio of 9/7. △ Less

Submitted 17 July, 2014; originally announced July 2014.

Comments: 12 page, 8 figure, ISSAC