Researcher profile

Calvin W. Johnson

Calvin W. Johnson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

A modified Brink-Axel hypothesis for astrophysical Gamow-Teller transitions

Weak interaction charged current transition strengths from highly excited nuclear states are fundamental ingredients for accurate modeling of compact object composition and dynamics, but are difficult to obtain either from experiment or theory. For lack of alternatives, calculations have often fallen back upon a generalized Brink-Axel hypothesis, that is, assuming the strength function (transition probability) is independent of the initial nuclear state but depends only upon the transition energy and the weak interaction properties of the parent nucleus ground state. Here we present numerical evidence for a modified `local' Brink-Axel hypothesis for Gamow-Teller transitions for $pf$-shell nuclei relevant to astrophysical applications. Specifically, while the original Brink-Axel hypothesis does not hold globally, strength functions from initial states nearby in energy are similar within statistical fluctuations. This agrees with previous work on strength function moments. Using this modified hypothesis, we can tackle strength functions at previously intractable initial energies, using semi-converged initial states at arbitrary excitation energy. Our work provides a well-founded method for computing accurate thermal weak transition rates for medium-mass nuclei at temperatures occurring in stellar cores near collapse. We finish by comparing to previous calculations of astrophysical rates.

preprint2022arXiv

An Alternative Approach to Quantum Imaginary Time Evolution

There is increasing interest in quantum algorithms that are based on the imaginary-time evolution (ITE), a successful classical numerical approach to obtain ground states. However, most of the proposals so far require heavy post-processing computational steps on a classical computer, such as solving linear equations. Here we provide an alternative approach to implement ITE. A key feature in our approach is the use of an orthogonal basis set: the propagated state is efficiently expressed in terms of orthogonal basis states at every step of the evolution. We argue that the number of basis states needed at those steps to achieve an accurate solution can be kept of the order of $n$, the number of qubits, by controlling the precision (number of significant digits) and the imaginary-time increment. The number of quantum gates per imaginary-time step is estimated to be polynomial in $n$. Additionally, while in many QAs the locality of the Hamiltonian is a key assumption, in our algorithm this restriction is not required. This characteristic of our algorithm renders it useful for studying highly nonlocal systems, such as the occupation-representation nuclear shell model. We illustrate our algorithm through numerical implementation on an IBM quantum simulator.

preprint2022arXiv

Collective neutrino oscillations with tensor networks using a time-dependent variational principle

If a system of flavor-oscillating neutrinos is at high enough densities that neutrino-neutrino coherent forward scatterings are non-negligible, the system becomes a time-dependent many-body problem. An important and open question is whether the flavor evolution is sufficiently described by a mean-field approach or can be strongly affected by correlations arising from two-body interactions in the neutrino Hamiltonian, as measured by nontrivial quantum entanglement. Numerical computations of the time evolution of many-body quantum systems are challenging because the size of the Hilbert space scales exponentially with the number of particles N in the system. Thus, it is important to investigate approximate but beyond-mean-field numerical treatments at larger values of N. Here we investigate the efficacy of tensor network methods to calculate the time evolution of interacting neutrinos at larger values of N than are possible with conventional methods. In particular, we introduce the use of time-dependent variational principle methods to address the long-range (in momentum space) interactions of the neutrino Hamiltonian when including many distinct vacuum oscillation frequencies. We also define new error measures based upon the instantaneously conserved charge operators known for this Hamiltonian to determine validity of large-N tensor network calculations.

preprint2022arXiv

New Insights into Backbending in the Symmetry-adapted Framework

We provide new insights into backbending phenomenon within the symmetry-adapted framework which naturally describes the intrinsic deformation of atomic nuclei. For $^{20}\text{Ne}$, the canonical example of backbending in light nuclei, the ab initio symmetry-adapted no-core shell model shows that while the energy spectrum replicates the backbending from experimental energies under the rigid rotor assumption, there is no change in the intrinsic deformation or intrinsic spin of the yrast band around the backbend. For the traditional example of $^{48}\text{Cr}$, computed in the valence shell with empirical interactions, we confirm a high-spin nucleus that is effectively spherical, in agreement with previous models. However, we find that this spherical distribution results, on average, from an almost equal mixing of deformed prolate shapes with deformed oblate and triaxial shapes. Microscopic calculations confirm the importance of spin alignment and configuration mixing, but surprisingly unveil no anomalous increase in moment of inertia. This finding opens the path toward further understanding the rotational behavior and moment of inertia of medium-mass nuclei.

preprint2022arXiv

Nuclear states projected from a pair condensate

Atomic nuclei exhibit deformation, pairing correlations, and rotational symmetries. To meet these competing demands in a computationally tractable formalism, we revisit the use of general pair condensates with good particle number as a trial wave function for even-even nuclei. After minimizing the energy of the condensate, we project out states with good angular momentum with a fast projection technique, allowing for general triaxial deformations. To show applicability, we present example calculations from pair condensates in several model spaces, and compare against projected Hartree-Fock and full configuration-interaction shell model calculations. This approach successfully generates spherical, vibrational and rotational spectra, demonstrating potential for modeling medium- to heavy-mass nuclei.

preprint2022arXiv

Nucleon-pair truncation of the shell model for medium-heavy nuclei

Background: Computationally tractable models of atomic nuclei is a long-time goal of nuclear structure physics. A flexible framework which easily includes excited states and many-body correlations is the configuration-interaction shell model (SM), but the exponential growth of the basis means one needs an efficient truncation scheme, ideally one that includes both deformation and pairing correlations. Purpose: We propose an efficient truncation scheme of the SM: starting from a pair condensate variationally defined by Hartree-Fock single-particle states and the particle-number conserved Bardeen-Cooper-Schrieffer (NBCS) approximation, we carry out projection of states with good angular momentum. Methods: After generating Hartree-Fock single-particle states with Kramers degeneracy in a SM space, we optimize the pair amplitudes in the NBCS by minimizing the energy, and then use linear algebra projection (LAP) of states with good angular momentum. Both NBCS and LAP are computationally fast. Results: Our calculations yield good agreement with full configuration-interaction SM calculations for low-lying states of transitional and rotational nuclei with axially symmetric and triaxial deformation in medium- and heavy-mass regions: $^{44,46,48}$Ti, $^{48,50}$Cr, $^{52}$Fe, $^{60,62,64}$Zn, $^{66,68}$Ge, $^{68}$Se, and $^{108,110}$Xe. We predict low-lying states of $^{112-114}\textrm{Ba}$ and $^{116-120}\textrm{Ce}$, nuclei difficult to reach by large-scale SM calculations. Conclusions: Both pair correlation and the configuration mixing between different intrinsic states play a key role in reproducing collectivity and shape coexistence, demonstrating the utility of this truncation scheme of the SM to study transitional and deformed nuclei.

preprint2021arXiv

Nucleon-pair coupling scheme in Elliott's SU(3) model

Elliott's SU(3) model is at the basis of the shell-model description of rotational motion in atomic nuclei. We demonstrate that SU(3) symmetry can be realized in a truncated shell-model space if constructed in terms of a sufficient number of collective $S$, $D$, $G$, $\dots$ pairs (i.e., with angular momentum zero, two, four, $\dots$) and if the structure of the pairs is optimally determined either by a conjugate-gradient minimization method or from a Hartree-Fock intrinsic state. We illustrate the procedure for 6 protons and 6 neutrons in the $pf$ ($sdg$) shell and exactly reproduce the level energies and electric quadrupole properties of the ground-state rotational band with $SDG$ ($SDGI$) pairs. The $SD$-pair approximation without significant renormalization, on the other hand, cannot describe the full SU(3) collectivity. A mapping from Elliott's fermionic SU(3) model to systems with $s$, $d$, $g$, $\dots$ bosons provides insight into the existence of a decoupled collective subspace in terms of $S$, $D$, $G$, $\dots$ pairs.

preprint2020arXiv

Exact sum rules with approximate ground states

Electromagnetic and weak transitions tell us a great deal about the structure of atomic nuclei. Yet modeling transitions can be difficult: it is often easier to compute the ground state, if only as an approximation, than excited states. One alternative is through transition sum rules, in particular the non-energy-weighted and energy-weighted sum rules, which can be computed as expectation values of operators. We investigate by computing sum rules for a variety of nuclei, comparing the numerically exact full configuration-interaction shell model, as a reference, to Hartree-Fock, projected Hartree-Fock, and the nucleon pair approximation. These approximations yield reasonable agreement, which we explain by prior work on the systematics of transition moments.

preprint2020arXiv

From deformed Hartree-Fock to the nucleon-pair approximation

The nucleon-pair approximation (NPA) can be a compact alternative to full configuration-interaction (FCI) diagonalization in nuclear shell-model spaces, but selecting good pairs is a long-standing problem. While seniority-based pairs work well for near-spherical nuclides, they do not work well for deformed nuclides with strong rotational bands. We propose an alternate approach. We show how one can write any Slater determinant for an even number of particles as a general pair condensate, from which one can project out pairs of good angular momentum. We implement this by generating unconstrained Hartree-Fock states in a shell model basis and extracting $S$, $D$, and $G$ pairs. The subsequent NPA calculations yield good agreement with FCI results using the same effective interactions.

preprint2020arXiv

Uncertainty quantification of an empirical shell-model interaction using principal component analysis

Recent investigations have emphasized the importance of uncertainty quantification (UQ) to describe errors in nuclear theory. We carry out UQ for configuration-interaction shell model calculations in the $1s$-$0d$ valence space, investigating the sensitivity of observables to perturbations in the 66 parameters (matrix elements) of a high-quality empirical interaction. The large parameter space makes computing the corresponding Hessian numerically costly, so we compare a cost-effective approximation, using the Feynman-Hellmann theorem, to the full Hessian and find it works well. Diagonalizing the Hessian yields the principal components of the interaction: linear combinations of parameters ordered by sensitivity. This approximately decoupled distribution of parameters facilitates theoretical error propagation onto structure observables: electromagnetic transitions, Gamow-Teller decays, and dark matter-nucleus scattering matrix elements.

preprint2020arXiv

Unmixing symmetries

The low-lying spectra of atomic nuclei display diverse behaviors, for example rotational bands, which can be described phenomenologically by simple symmetry groups such as spatial SU(3). This leads to the idea of dynamical symmetry, where the Hamiltonian commutes with the Casimir operator(s) of a group, and is block-diagonal in subspaces defined by the group's irreducible representations or irreps. Detailed microscopic calculations, however, show these symmetries are in fact often strongly mixed and the wave function fragmented across many irreps. More commonly the fragmentation across members of a band are similar, or a quasi-dynamical symmetry. In this Letter I explicitly, albeit numerically, construct unitary transformations from a quasi-dynamical symmetry to a dynamical symmetry, adapting the similarity renormalization group, or SRG. The standard SRG produces unsatisfactory results, forcing the induced dynamical symmetry to be dominated by high-weight irreps irrespective of the original decomposition. Using spectral distribution theory to rederive and diagnose standard SRG, I introduce a new form of SRG. The new SRG transforms a quasi-dynamical symmetry to a dynamical symmetry, that is, unmixes the mixed symmetries, with intuitively more appealing results.