Local Chiral Interactions Research Articles

Path-integral quantum Monte Carlo calculations of light nuclei

We describe a path-integral ground-state quantum Monte Carlo method for light nuclei in continuous space. We show how to efficiently update and sample the paths with spin-isospin dependent and spin-orbit interactions. We apply the method to the triton and alpha particle using both local chiral interactions with next-to-next-to-leading-order %(N$^2$LO) and the Argonne interactions. For operators, like the total energy, that commute with the Hamiltonian, our results agree with Green's function Monte Carlo and auxiliary field diffusion Monte Carlo calculations. For operators that do not commute with the Hamiltonian and for Euclidean response functions, the path-integral formulation allows straightforward calculation without forward walking or the increased variance typical of diffusion methods. We demonstrate this by calculating density distributions, root mean square radii, and Euclidean response functions for single-nucleon couplings.

Implementation of Local Chiral Interactions in the Hyperspherical Harmonics Formalism

With the goal of using chiral interactions at various orders to explore the properties of the few-body nuclear systems, we write the recently developed local chiral interactions as spherical irreducible tensors and implement them in the hyperspherical harmonics expansion method. We devote particular attention to three-body forces at next-to-next-to leading order, which play an important role in reproducing experimental data. We check our implementation by benchmarking the ground-state properties of 3H, 3He, and 4He against the available Monte Carlo calculations. We then confirm their order-by-order truncation error estimates and further investigate uncertainties in the charge radii obtained by using the precise muonic atom data for single-nucleon radii. Having local chiral Hamiltonians at various orders implemented in our hyperspherical harmonics suites of codes opens up the possibility to test such interactions on other light-nuclei properties, such as electromagnetic reactions.

Nuclear and neutron-star matter from local chiral interactions

We report a quantum Monte Carlo calculation of the equation of state of symmetric nuclear matter using local interactions derived from chiral effective field theory up to next-to-next-to-leading order fit to few-body observables only. The empirical saturation density and energy are well reproduced within statistical and systematic uncertainties. We have also derived the symmetry energy as a function of the density, finding good agreement with available experimentally derived constraints at saturation and twice saturation density. We find that the corresponding pressure is also in excellent agreement with recent constraints extracted from gravitational waves of the neutron-star merger GW170817 by the LIGO-Virgo detection.

Local chiral interactions and magnetic structure of few-nucleon systems

The magnetic form factors of $^2$H, $^3$H, and $^3$He, deuteron photodisintegration cross sections at low energies, and deuteron threshold electrodisintegration cross sections at backward angles in a wide range of momentum transfers, are calculated with the chiral two-nucleon (and three-nucleon) interactions including $\Delta$ intermediate states that have recently been constructed in configuration space. The $A\,$=$\,$3 wave functions are obtained from hyperspherical-harmonics solutions of the Schr\"odinger equation. The electromagnetic current includes one- and two-body terms, the latter induced by one- and two-pion exchange (OPE and TPE, respectively) mechanisms and contact interactions. The contributions associated with $\Delta$ intermediate states are only retained at the OPE level, and are neglected in TPE loop (tree-level) corrections to two-body (three-body) current operators. Expressions for these currents are derived and regularized in configuration space for consistency with the interactions. The low-energy constants that enter the contact few-nucleon systems. The predicted form factors and deuteron electrodisintegration cross section are in excellent agreement with experiment for momentum transfers up to 2--3 fm$^{-1}$. However, the experimental values for the deuteron photodisintegration cross section are consistently underestimated by theory, unless use is made of the Siegert form of the electric dipole transition operator. A complete analysis of the results is provided, including the clarification of the origin of the aforementioned discrepancy.

Local chiral interactions, the tritium Gamow-Teller matrix element, and the three-nucleon contact term

The Gamow-Teller (GT) matrix element contributing to tritium $\ensuremath{\beta}$ decay is calculated with trinucleon wave functions obtained from hyperspherical-harmonics solutions of the Schr\"odinger equation with the chiral two- and three-nucleon interactions including $\mathrm{\ensuremath{\Delta}}$ intermediate states that have recently been constructed in configuration space. Predictions up to next-to-next-to-next-to-leading order (N3LO) in the chiral expansion of the axial current (with $\mathrm{\ensuremath{\Delta}}$'s) overestimate the empirical value by 1--4%. By exploiting the relation between the low-energy constant (LEC) in the contact three-nucleon interaction and two-body axial current, we provide new determinations of the LECs ${c}_{D}$ and ${c}_{E}$ that characterize this interaction by fitting the trinucleon binding energy and tritium GT matrix element. Some of the implications that the resulting models of three-nucleon interactions have on the spectra of light nuclei and the equation of state of neutron matter are briefly discussed. We also provide a partial analysis, which ignores $\mathrm{\ensuremath{\Delta}}$'s, of the contributions due to loop corrections in the axial current at next-to-next-to-next-to-next-to-leading order (N4LO). Finally, explicit expressions for the axial current up to N4LO have been derived in configuration space, which other researchers in the field may find useful.

Single- and two-nucleon momentum distributions for local chiral interactions

We present quantum Monte Carlo calculations of the single- and two-nucleon momentum distributions in selected nuclei for $A\le16$. We employ local chiral interactions at next-to-next-to-leading order. We find good agreement at low momentum with the single-nucleon momentum distributions derived for phenomenological potentials. The same agreement is found for the integrated two-nucleon momentum distributions at low relative momentum $q$ and low center-of-mass momentum $Q$. We provide results for the two-nucleon momentum distributions as a function of both $q$ and $Q$. The large ratio of $pn$ to $pp$ pairs around $q=2\,\rm fm^{-1}$ for back-to-back $(Q=0)$ pairs is confirmed up to $^{16}$O, and results are compatible with those extracted from available experimental data.

Auxiliary field diffusion Monte Carlo calculations of light and medium-mass nuclei with local chiral interactions

Quantum Monte Carlo methods have recently been employed to study properties of nuclei and infinite matter using local chiral effective field theory interactions. In this work, we present a detailed description of the auxiliary field diffusion Monte Carlo algorithm for nuclei in combination with local chiral two- and three-nucleon interactions up to next-to-next-to-leading order. We show results for the binding energy, charge radius, charge form factor, and Coulomb sum rule in nuclei with $3\le A\le16$. Particular attention is devoted to the effect of different operator structures in the three-body force for different cutoffs. The outcomes suggest that local chiral interactions fit to few-body observables give a very good description of the ground-state properties of nuclei up to $^{16}$O, with the exception of one fit for the softer cutoff which predicts overbinding in larger nuclei.

Light-Nuclei Spectra from Chiral Dynamics.

In recent years local chiral interactions have been derived and implemented in quantum MonteCarlo methods in order to test to what extent the chiral effective field theory framework impacts our knowledge of few- and many-body systems. In this Letter, we present Green's function MonteCarlo calculations of light nuclei based on the family of local two-body interactions presented by our group in a previous paper in conjunction with chiral three-body interactions fitted to bound- and scattering-state observables in the three-nucleon sector. These interactions include Δ intermediate states in their two-pion-exchange components. We obtain predictions for the energy levels and level ordering of nuclei in the mass range A=4-12, accurate to ≤2% of the binding energy, in very satisfactory agreement with experimental data.

Is a Trineutron Resonance Lower in Energy than a Tetraneutron Resonance?

We present quantum MonteCarlo calculations of few-neutron systems confined in external potentials based on local chiral interactions at next-to-next-to-leading order in chiral effective field theory. The energy and radial densities for these systems are calculated in different external Woods-Saxon potentials. We assume that their extrapolation to zero external-potential depth provides a quantitative estimate of three- and four-neutron resonances. The validity of this assumption is demonstrated by benchmarking with an exact diagonalization in the two-body case. We find that the extrapolated trineutron resonance, as well as the energy for shallow well depths, is lower than the tetraneutron resonance energy. This suggests that a three-neutron resonance exists below a four-neutron resonance in nature and is potentially measurable. To confirm that the relative ordering of three- and four-neutron resonances is not an artifact of the external confinement, we test that the odd-even staggering in the helium isotopic chain is reproduced within this approach. Finally, we discuss similarities between our results and ultracold Fermi gases.

Reconfigurable self-assembly through chiral control of interfacial tension

From determining the optical properties of simple molecular crystals to establishing the preferred handedness in highly complex vertebrates, molecular chirality profoundly influences the structural, mechanical and optical properties of both synthetic and biological matter on macroscopic length scales. In soft materials such as amphiphilic lipids and liquid crystals, the competition between local chiral interactions and global constraints imposed by the geometry of the self-assembled structures leads to frustration and the assembly of unique materials. An example of particular interest is smectic liquid crystals, where the two-dimensional layered geometry cannot support twist and chirality is consequently expelled to the edges in a manner analogous to the expulsion of a magnetic field from superconductors. Here we demonstrate a consequence of this geometric frustration that leads to a new design principle for the assembly of chiral molecules. Using a model system of colloidal membranes, we show that molecular chirality can control the interfacial tension, an important property of multi-component mixtures. This suggests an analogy between chiral twist, which is expelled to the edges of two-dimensional membranes, and amphiphilic surfactants, which are expelled to oil-water interfaces. As with surfactants, chiral control of interfacial tension drives the formation of many polymorphic assemblages such as twisted ribbons with linear and circular topologies, starfish membranes, and double and triple helices. Tuning molecular chirality in situ allows dynamical control of line tension, which powers polymorphic transitions between various chiral structures. These findings outline a general strategy for the assembly of reconfigurable chiral materials that can easily be moved, stretched, attached to one another and transformed between multiple conformational states, thus allowing precise assembly and nanosculpting of highly dynamical and designable materials with complex topologies.