JEJonathan Engel

Jonathan Engel

*This is an algorithmically generated author record based on publicly available information.

Current Organizations / Institutes

The organizations/institutes from the scholar's latest article

Previous Organizations / Institutes

The organizations / institutes from the scholar's other articles in the past ten years


Jonathan Engel has a total of 235 co-authors. Here are the top co-authors by number of shared publications:


Latest Publications

  • Open Access

    Lanczos Method for Quasiparticle Random-Phase Approximation Strength Functions in Atomic Nuclei

    by Dong Min Roh, Chao Yang, Jonathan Engel, Matthew Dai

    We present a Lanczos method for computing charge-changing quasiparticle random-phase approximation (QRPA) strength functions in atomic nuclei. Starting from the finite amplitude method (FAM) formulation of the QRPA linear-response problem, we derive equivalent spectral representations and, in the real case, a reduced eigenvalue problem involving the matrix products $MK$ and $KM$, where $M\equiv A+B$ and $K\equiv A-B$ are formed from the usual QRPA matrices $A$ and $B$. The resulting formulation enables a matrix-free Lanczos approximation of the Lorentzian-smeared strength function over a broad energy interval from a single run, in contrast to conventional frequency-by-frequency response calculations. Numerical tests for $^{112}$Sn and $^{150}$Nd first show that the generalized minimal residual method (GMRES) reproduces the converged iterative FAM strength profiles while requiring fewer iterations. Using GMRES as the frequency-by-frequency reference, we then show that the Lanczos approximation reproduces the same strength profiles with reduced overall cost. These results indicate that Lanczos projection provides an efficient and accurate approach for QRPA strength-function calculations when spectral information is required over an extended frequency range.

  • Open Access

    Detecting the QCD Axion via the Ferroaxionic Force with Piezoelectric Materials

    by Asimina Arvanitaki, Jonathan Engel, Andrew A. Geraci, Alexander Hepburn, Amalia Madden, Ken Van Tilburg

    We show that piezoelectric materials can be used to source virtual QCD axions, generating a new axion-mediated force. Spontaneous parity violation within the piezoelectric crystal combined with time-reversal violation from aligned spins provide the necessary symmetry breaking to produce an effective in-medium scalar coupling of the axion to nucleons up to seven orders of magnitude larger than that in vacuum. We propose a detection scheme based on nuclear spin precession caused by the axion’s pseudoscalar coupling to nuclear spins. This signal is resonantly enhanced when the distance between the source crystal and the spin sample is modulated at the spin precession frequency. Using this effect, future experimental setups can be sensitive to the QCD axion in the unexplored mass range from 10 − 5 eV to 10 − 2 eV .


  • Open Access

    Effects of beyond-mean-field correlations on nuclear Schiff moments

    by E. F. Zhou, J. M. Yao, J. Engel, J. Meng

    We compute the nuclear Schiff moments of the diamagnetic atoms $^{129}$Xe, $^{199}$Hg, and $^{225}$Ra in multireference covariant density functional theory. Beyond-mean-field correlations, arising from symmetry restoration and shape mixing, are incorporated via the generator coordinate method with projection onto states with well-defined parity, particle number, and angular momentum. Our results reveal a correlation between the contributions of nuclear intermediate states to Schiff moments and the electric dipole transition strengths from these states to the ground state. The new beyond-mean-field effects can either enhance or suppress the Schiff moments. In $^{225}$Ra, they do the latter, reducing the enhancement from octupole deformation somewhat.


  • Open Access

    Nuclear Schiff Moments and CP Violation

    by Jonathan Engel

    This paper reviews the calculation of nuclear Schiff moments, which one must know in order to interpret experiments that search for time-reversal-violating electric dipole moments in certain atoms and molecules. After briefly reviewing the connection between dipole moments and CP violation in and beyond the Standard Model of particle physics, Schiff's theorem, which concerns the screening of nuclear electric dipole moments by electrons, Schiff moments, and experiments to measure dipole moments in atoms and molecules, the paper examines attempts to compute Schiff moments in nuclei such as $^{199}$Hg and octupole-deformed isotopes such as $^{225}$Ra, which are particularly useful in experiments. It then turns to ab initio nuclear-structure theory, describing ways in which both the In-Medium Similarity Renormalization Group and coupled-cluster theory can be used to compute important Schiff moments more accurately than the less controlled methods that have been applied so far.


  • Top Cited Papers

    Ab Initio Uncertainty Quantification of Neutrinoless Double-Beta Decay in Ge76

    by A. Belley, J. M. Yao, B. Bally, J. Pitcher, J. Engel, H. Hergert, J. D. Holt, T. Miyagi, T. R. Rodríguez, A. M. Romero, et al.

    The observation of neutrinoless double-beta (0νββ) decay would offer proof of lepton number violation, demonstrating that neutrinos are Majorana particles, while also helping us understand why there is more matter than antimatter in the Universe. If the decay is driven by the exchange of the three known light neutrinos, a discovery would, in addition, link the observed decay rate to the neutrino mass scale through a theoretical quantity known as the nuclear matrix element (NME). Accurate values of the NMEs for all nuclei considered for use in 0νββ experiments are therefore crucial for designing and interpreting those experiments. Here, we report the first comprehensive ab initio uncertainty quantification of the 0νββ-decay NME, in the key nucleus Ge76. Our method employs nuclear strong and weak interactions derived within chiral effective field theory and recently developed many-body emulators. Our result, with a conservative treatment of uncertainty, is an NME of 2.60−1.36+1.28, which, together with the best-existing half-life sensitivity and phase-space factor, sets an upper limit for effective neutrino mass of 187−62+205  meV. The result is important for designing next-generation germanium detectors aiming to cover the entire inverted hierarchy region of neutrino masses.


  • Effects of quasiparticle-vibration coupling on Gamow-Teller strength and β decay with the Skyrme proton-neutron finite-amplitude method

    by Qunqun Liu, Jonathan Engel, Nobuo Hinohara, Markus Kortelainen

    We adapt the proton-neutron finite-amplitude method, which in its original form is an efficient implementation of the Skyrme quasiparticle random phase approximation, to include the coupling of quasiparticles to like-particle phonons. The approach allows us to add beyond-quasiparticle random-phase approximation correlations to computations of Gamow-Teller strength and β-decay rates in deformed nuclei for the first time. We test the approach in several deformed isotopes for which measured strength distributions are available. The additional correlations dramatically improve agreement with the data, and will lead to improved global β-decay rates.


  • Open Access

    The Influence of β-decay Rates on r-process Observables

    by Kelsey A. Lund, J. Engel, G. C. McLaughlin, M. R. Mumpower, E. M. Ney, R. Surman

    The rapid neutron capture process (r-process) is one of the main mechanisms whereby elements heavier than iron are synthesized, and is entirely responsible for the natural production of the actinides. Kilonova emissions are modeled as being largely powered by the radioactive decay of species synthesized via the r-process. Given that the r-process occurs far from nuclear stability, unmeasured beta-decay rates play an essential role in setting the timescale for the r-process. In an effort to better understand the sensitivity of kilonova modeling to different theoretical global beta-decay descriptions, we incorporate these into nucleosynthesis calculations. We compare the results of these calculations and highlight differences in kilonova nuclear energy generation and light-curve predictions, as well as final abundances and their implications for nuclear cosmochronometry. We investigate scenarios where differences in beta-decay rates are responsible for increased nuclear heating on timescales of days that propagates into a significantly increased average bolometric luminosity between 1 and 10 days post-merger. We identify key nuclei, both measured and unmeasured, whose decay rates directly impact nuclear heating generation on timescales responsible for light-curve evolution. We also find that uncertainties in beta-decay rates significantly impact age estimates from cosmochronometry.


  • Ab initio studies of the double–Gamow-Teller transition and its correlation with neutrinoless double- β decay

    by J. M. Yao, I. Ginnett, A. Belley, T. Miyagi, R. Wirth, S. Bogner, J. Engel, H. Hergert, J. D. Holt, S. R. Stroberg

    We use chiral interactions and several ab initio methods to compute the nuclear matrix elements (NMEs) for ground-state-to-ground-state double Gamow-Teller transitions in a range of isotopes and explore the correlation of these NMEs with those for neutrinoless double beta decay produced by the exchange of a light Majorana neutrino. When all the NMEs of both isospin-conserving and isospin-changing transitions from the ab initio calculations are considered, the correlation is strong. For the experimentally relevant isospin-changing transitions by themselves, however, the correlation is weaker and may not be helpful for reducing the uncertainty in the NMEs for neutrinoless double beta decay.


  • Open Access

    Solving Nuclear Structure Problems with the Adaptive Variational Quantum Algorithm

    by A. M. Romero, J. Engel, Ho Lun Tang, Sophia E. Economou

    We use the Lipkin-Meshkov-Glick (LMG) model and the valence-space nuclear shell model to examine the likely performance of variational quantum eigensolvers in nuclear-structure theory. The LMG model exhibits both a phase transition and spontaneous symmetry breaking at the mean-field level in one of the phases, features that characterize collective dynamics in medium-mass and heavy nuclei. We show that with appropriate modifications, the ADAPT-VQE algorithm, a particularly flexible and accurate variational approach, is not troubled by these complications. We treat up to 12 particles and show that the number of quantum operations needed to approach the ground-state energy scales linearly with the number of qubits. We find similar scaling when the algorithm is applied to the nuclear shell model with realistic interactions in the $sd$ and $pf$ shells. Although most of these simulations contain no noise, we use a noise model from real IBM hardware to show that for the LMG model with four particles, weak noise has no effect on the efficiency of the algorithm.


  • Top Cited Papers

    Solving nuclear structure problems with the adaptive variational quantum algorithm

    by A. M. Romero, J. Engel, Ho Lun Tang, Sophia E. Economou

    We use the Lipkin-Meshkov-Glick (LMG) model and the valence-space nuclear shell model to examine the likely performance of variational quantum eigensolvers in nuclear-structure theory. The LMG model exhibits both a phase transition and spontaneous symmetry breaking at the mean-field level in one of the phases, features that characterize collective dynamics in medium-mass and heavy nuclei. We show that with appropriate modifications, the ADAPT-VQE algorithm, a particularly flexible and accurate variational approach, is not troubled by these complications. We treat up to 12 particles and show that the number of quantum operations needed to approach the ground-state energy scales linearly with the number of qubits. We find similar scaling when the algorithm is applied to the nuclear shell model with realistic interactions in the sd and pf shells. Although most of these simulations contain no noise, we use a noise model from real IBM hardware to show that for the LMG model with four particles, weak noise has no effect on the efficiency of the algorithm.