BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Calvin W. Johnson

Calvin W. Johnson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
nucl-thastro-ph.SRphysics.comp-phcond-mat.quant-gasmath-phmath.MPphysics.atom-phquant-phhep-phnucl-exphysics.data-an

Catalog footprint

What is connected

27works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

27 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Operator evolution for ab initio electric dipole transitions of 4He

A goal of nuclear theory is to make quantitative predictions of low-energy nuclear observables starting from accurate microscopic internucleon forces. A major element of such an effort is applying unitary transformations to soften the nuclear Hamiltonian and hence accelerate the convergence of ab initio calculations as a function of the model space size. The consistent simultaneous transformation of external operators, however, has been overlooked in applications of the theory, particularly for nonscalar transitions. We study the evolution of the electric dipole operator in the framework of the similarity renormalization group method and apply the renormalized matrix elements to the calculation of the 4He total photoabsorption cross section and electric dipole polarizability. All observables are calculated within the ab initio no-core shell model. We find that, although seemingly small, the effects of evolved operators on the photoabsorption cross section are comparable in magnitude to the correction produced by including the chiral three-nucleon force and cannot be neglected.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spin-orbit decomposition of ab initio wavefunctions

Although the modern shell-model picture of atomic nuclei is built from single-particle orbits with good total angular momentum $j$, leading to $j$-$j$ coupling, phenomenological models suggested decades ago that for $0p$-shell nuclides a simpler picture can be realized via coupling of total spin $S$ and total orbital angular momentum $L$. I revisit this idea with large-basis, no-core shell model (NCSM) calculations using modern \textit{ab initio} two-body interactions, and dissect the resulting wavefunctions into their component $L$- and $S$-components. Remarkably, there is broad agreement with calculations using the phenomenological Cohen-Kurath forces, despite a gap of nearly fifty years and six orders of magnitude in basis dimensions. I suggest $L$-$S$ may be a useful tool for analyzing \textit{ab initio} wavefunctions of light nuclei, for example in the case of rotational bands.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Systematics of strength function sum rules

Sum rules provide useful insights into transition strength functions and are often expressed as expectation values of an operator. In this letter I demonstrate that non-energy-weighted transition sum rules have strong secular dependences on the energy of the initial state. Such non-trivial systematics have consequences: the simplification suggested by the generalized Brink-Axel hypothesis, for example, does not hold for most cases, though it weakly holds in at least some cases for electric dipole transitions. Furthermore, I show the systematics can be understood through spectral distribution theory, calculated via traces of operators and of products of operators. Seen through this lens, violation of the generalized Brink-Axel hypothesis is unsurprising: one \textit{expects} sum rules to evolve with excitation energy. Furthermore, to lowest order the slope of the secular evolution can be traced to a component of the Hamiltonian being positive (repulsive) or negative (attractive).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Operator evolution for ab initio nuclear theory

The past two decades have seen a revolution in ab initio calculations of nuclear properties. One key element has been the development of a rigorous effective interaction theory, applying unitary transformations to soften the nuclear Hamiltonian and hence accelerate the convergence as a function of the model space size. For consistency, however, one ought to apply the same transformation to other operators when calculating transitions and mean values from the eigenstates of the renormalized Hamiltonian. Working in a translationally-invariant harmonic oscillator basis for the two- and three-nucleon systems, we evolve the Hamiltonian, square-radius and total dipole strength operators by the similarity renormalization group (SRG). The inclusion of up to three-body matrix elements in the 4He nucleus all but completely restores the invariance of the expectation values under the transformation. We also consider a Gaussian operator with adjustable range and find at short ranges an increased contribution from such induced three-body terms.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Testing the spin-cutoff parameterization with shell-model calculations

The nuclear level density, an important input to Hauser-Feshbach calculations, depends not only on excitation energy but also on angular momentum J. The J-dependence of the level density at fixed excitation energy E_x is usually parameterized via the spin-cutoff factor sigma. We carefully test the statistical accuracy of this parameterization for a large number of spectra computed using semi-realistic interactions in the interacting shell model, with a nonlinear least-squares fit of sigma and finding the error bar in sigma. The spin-cutoff parameterization works well as long as there are enough states to be statistical. In turn, the spin-cutoff factor can be related to the average value of J^2 at a fixed excitation energy, and we briefly investigate extracting <J^2 (E_x)> from a thermal calculation such as one might do via Monte Carlo.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Comparison of Gamow-Teller strengths in the random phase approximation

The Gamow-Teller response is astrophysically important for a number of nuclides, particularly around iron. The random phase approximation (RPA) is an efficient way to generate strength distributions. In order to better understand both theoretical systematics and uncertainties, we compare the Gamow-Teller strength distributions for a suite of nuclides and for a suite of interactions, including semi-realistic interactions in the $1p$-$0f$ space with the RPA and a separable multi-shell interaction in the quasi-particle RPA. We also compare with experimental results for GT$_-$ on ${54}$Fe.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Computational Nuclear Quantum Many-Body Problem: The UNEDF Project

The UNEDF project was a large-scale collaborative effort that applied high-performance computing to the nuclear quantum many-body problem. UNEDF demonstrated that close associations among nuclear physicists, mathematicians, and computer scientists can lead to novel physics outcomes built on algorithmic innovations and computational developments. This review showcases a wide range of UNEDF science results to illustrate this interplay.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Excitation spectra from angular momentum projection of Hartree-Fock states and the configuration-interaction shell-model

We make numerical comparison of spectra from angular-momentum projection on Hartree-Fock states with spectra from configuration-interaction nuclear shell-model calculations, all carried out in the same model spaces (in this case the sd, lower pf, and p-sd_5/2 shells) and using the same input Hamiltonians. We find, unsurprisingly, that the low-lying excitation spectra for rotational nuclides are well reproduced, but the spectra for vibrational nuclides, and more generally the complex specta for odd-A and odd-odd nuclides are less well reproduced in detail.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Factorization in large-scale many-body calculations

One approach for solving interacting many-fermion systems is the configuration-interaction method, also sometimes called the interacting shell model, where one finds eigenvalues of the Hamiltonian in a many-body basis of Slater determinants (antisymmeterized products of single-particle wavefunctions). The resulting Hamiltonian matrix is typically very sparse, but for large systems the nonzero matrix elements can nonetheless require terabytes or more of storage. An alternate algorithm, applicable to a broad class of systems with symmetry, in our case rotational invariance, is to exactly factorize both the basis and the interaction using additive/multiplicative quantum numbers; such an algorithm can reduce the storage requirements by an order of magnitude or more. We discuss factorization in general as well as in the context of a specific configuration-interaction code, BIGSTICK, which runs both on serial and parallel machines.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Many-body approximations for atomic binding energies

We benchmark three standard approximations for the many-body problem -- the Hartree-Fock, projected Hartree-Fock, and random phase approximations -- against full numerical configuration-interaction calculations of the electronic structure of atoms, from Li through to Ne. These configuration-interaction calculations used up to $2 \times 10^8$ uncoupled basis states, equivalent to $ 10^7$ coupled basis states (configuration state functions.) Each method uses exactly the same input, i.e., the same single-particle basis and Coulomb matrix elements, so any differences are strictly due to the approximation itself. Although it consistently overestimates the ground state binding energy, the random phase approximation has the smallest overall errors; furthermore, we suggest it may be useful as a method for efficient optimization of single-particle basis functions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The origin of order in random matrices with symmetries

From Noether's theorem we know symmetries lead to conservation laws. What is left to nature is the ordering of conserved quantities; for example, the quantum numbers of the ground state. In physical systems the ground state is generally associated with `low' quantum numbers and symmetric, low-dimensional irreps, but there is no \textit{a priori} reason to expect this. By constructing random matrices with nontrivial point-group symmetries, I find the ground state is always dominated by extremal low-dimensional irreps. Going further, I suggest this explains the dominance of J=0 g.s. even for random two-body interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Random matrices, symmetries, and many-body states

All nuclei with even numbers of protons and of neutrons have ground states with zero angular momentum. This is ascribed to the pairing force between nucleons, but simulations with random interactions suggest a much broader many-body phenomenon. In this Letter I project out random Hermitian matrices that have good quantum numbers and, computing the width of the Hamiltonian in subspaces, find ground states dominated by low quantum numbers, e.g. J=0. Furthermore I find odd-$Z$, odd-$N$ systems with isospin conservation have relatively fewer J=0 ground states.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Many-body fits of phase-equivalent effective interactions

In many-body theory it is often useful to renormalize short-distance, high-momentum components of an interaction via unitary transformations. Such transformations preserve the on-shell physical observables of the two-body system (mostly phase-shifts, hence unitarily-connected effective interactions are often called phase-equivalent), while modifying off-shell T-matrix elements influential in many-body systems. In this paper I lay out a general and systematic approach for controlling the off-shell behavior of an effective interaction, which can be adjusted to many-body properties, and present an application to trapped fermions at the unitary

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Sensitivity analysis of random two-body interactions

The input to the configuration-interaction shell model includes many dozens or hundreds of independent two-body matrix elements. Previous studies have shown that when fitting to experimental low-lying spectra, the greatest sensitivity is to only a few linear combinations of matrix elements. Here we consider interactions drawn from the two-body random ensemble, or TBRE, and find that the low-lying spectra are also most sensitive to only a few linear combinations of two-body matrix elements, in a fashion nearly indistinguishable from an interaction empirically fit to data. We find in particular the spectra for both the random and empirical interactions are sensitive to similar matrix elements, which we analyze using monopole and contact interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

A new puzzle for random interaction

We continue a series of numerical experiments on many-body systems with random two-body interactions, by examining correlations in ratios in excitation energies of yrast $J$ = 0, 2, 4, 6, 8 states. Previous studies, limited only to $J$ = 0,2,4 states, had shown strong correlations in boson systems but not fermion systems. By including $J \ge 6$ states and considering different scatter plots, strong and realistic correlations appear in both boson and fermion systems. Such correlations are a challenge to explanations of random interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Behavior of shell-model configuration moments

An important input into reaction theory is the density of states or the level density. Spectral distribution theory (also known as nuclear statistical spectroscopy) characterizes the secular behavior of the density of states through moments of the Hamiltonian. One particular approach is to partition the model space into subspaces and find the moments in those subspaces; a convenient choice of subspaces are spherical shell-model configurations. We revisit these configuration moments and find, for complete $0\hbarω$ many-body spaces, the following behaviors: (a) the configuration width is nearly constant for all configurations; (b) the configuration asymmetry or third moment is strongly correlated with the configuration centroid; (c) the configuration fourth moment, or excess is linearly related to the square to the configuration asymmetry. Such universal behavior may allow for more efficient modeling of the density of states in a shell-model framework.

0 citations0 reviewsNo code linkedOpen access