We investigate the impact of the QCD vacuum at nonzero $\theta$ on the properties of light nuclei, Big Bang nucleosynthesis, and stellar nucleosynthesis. Our analysis starts with a calculation of the $\theta$-dependence of the neutron-proton mass difference and neutron decay using chiral perturbation theory. We then discuss the $\theta$-dependence of the nucleon-nucleon interaction using a one-boson-exchange model and compute the properties of the two-nucleon system. Using the universal properties of four-component fermions at large scattering length, we then deduce the binding energies of the three-nucleon and four-nucleon systems. Based on these results, we discuss the implications for primordial abundances of light nuclei, the production of nuclei in stellar environments, and implications for an anthropic view of the universe.

We derive the axion-nucleon interaction Lagrangian in heavy baryon chiral perturbation theory up to next-to-next-to-leading order. The effective axion-nucleon coupling is calculated to a few percent accuracy.

The triple-alpha process, whereby evolved stars create carbon and oxygen, is believed to be fine-tuned to a high degree. Such fine-tuning is suggested by the unusually strong temperature dependence of the triple-alpha reaction rate at stellar temperatures. This sensitivity is due to the resonant character of the triple-alpha process, which proceeds through the so-called "Hoyle state" of $^{12}$C with spin-parity $0^+$. The question of fine-tuning can be studied within the {\it ab initio} framework of nuclear lattice effective field theory, which makes it possible to relate {\it ad hoc} changes in the energy of the Hoyle state to changes in the fundamental parameters of the nuclear Hamiltonian, which are the light quark mass $m_q$ and the electromagnetic fine-structure constant. Here, we update the effective field theory calculation of the sensitivity of the triple-alpha process to small changes in the fundamental parameters. In particular, we consider recent high-precision lattice QCD calculations of the nucleon axial coupling $g_A$, as well as new and more comprehensive results from stellar simulations of the production of carbon and oxygen. While the updated stellar simulations allow for much larger {\it ad hoc} shifts in the Hoyle state energy than previously thought, recent lattice QCD results for the nucleon S-wave singlet and triplet scattering lengths now disfavor the scenario of no fine-tuning in the light quark mass $m_q$.

Nuclear clustering describes the appearance of structures resembling smaller nuclei such as alpha particles (4He nuclei) within the interior of a larger nucleus. While clustering is important for several well-known examples, much remains to be discovered about the general nature of clustering in nuclei. In this letter we present lattice Monte Carlo calculations based on chiral effective field theory for the ground states of helium, beryllium, carbon, and oxygen isotopes. By computing model-independent measures that probe three- and four-nucleon correlations at short distances, we determine the shape of the alpha clusters and the entanglement of nucleons comprising each alpha cluster with the outside medium. We also introduce a new computational approach called the pinhole algorithm, which solves a long-standing deficiency of auxiliary-field Monte Carlo simulations in computing density correlations relative to the center of mass. We use the pinhole algorithm to determine the proton and neutron density distributions and the geometry of cluster correlations in 12C, 14C, and 16C. The structural similarities among the carbon isotopes suggest that 14C and 16C have excitations analogous to the well-known Hoyle state resonance in 12C.

Brueckner theory is used to investigate the in-medium properties of a $\Lambda$-hyperon in nuclear and neutron matter, based on hyperon-nucleon interactions derived within SU(3) chiral effective field theory (EFT). It is shown that the resulting $\Lambda$ single-particle potential $U_\Lambda(p_\Lambda =0,\rho)$ becomes strongly repulsive for densities $\rho$ of two-to-three times that of normal nuclear matter. Adding a density-dependent effective $\Lambda N$-interaction constructed from chiral $\Lambda NN$ three-body forces increases the repulsion further. Consequences of these findings for neutron stars are discussed. It is argued that for hyperon-nuclear interactions with properties such as those deduced from the SU(3) EFT potentials, the onset for hyperon formation in the core of neutron stars is expected to be shifted to extremely high baryon density, thus potentially resolving the so-called hyperon puzzle.

We present a determination of the pion-nucleon ($\pi N$) $\sigma$-term $\sigma_{\pi N}$ based on the Cheng-Dashen low-energy theorem (LET), taking advantage of the recent high-precision data from pionic atoms to pin down the $\pi N$ scattering lengths as well as of constraints from analyticity, unitarity, and crossing symmetry in the form of Roy-Steiner equations to perform the extrapolation to the Cheng-Dashen point in a reliable manner. With isospin-violating corrections included both in the scattering lengths and the LET, we obtain $\sigma_{\pi N}=(59.1\pm 1.9\pm 3.0)$ MeV $=(59.1\pm 3.5)$ MeV, where the first error refers to uncertainties in the $\pi N$ amplitude and the second to the LET. Consequences for the scalar nucleon couplings relevant for the direct detection of dark matter are discussed.

I discuss the fine-tuning of the nuclear forces and in the formation of nuclei in the production of the elements in the Big Bang and in stars.

The quark mass dependences of light element binding energies and nuclear scattering lengths are derived using chiral perturbation theory in combination with non-perturbative methods. In particular, we present new, improved values for the quark mass dependence of meson resonances that enter the nuclear force. A detailed analysis of the theoretical uncertainties arising in this determination is presented. As an application we derive from a comparison of observed and calculated primordial deuterium and helium abundances a stringent limit on the variation of the light quark mass, $\delta m_q/m_q = 0.02 \pm 0.04$. Inclusion of the neutron lifetime modification under the assumption of a variation of the Higgs vacuum expectation value that translates into changing quark, electron, and weak gauge boson masses, leads to a stronger limit, $|\delta m_q/m_q| < 0.009$.

The excited state of the 12C nucleus known as the "Hoyle state" constitutes one of the most interesting, difficult and timely challenges in nuclear physics, as it plays a key role in the production of carbon via fusion of three alpha particles in red giant stars. In this letter, we present ab initio lattice calculations which unravel the structure of the Hoyle state, along with evidence for a low-lying spin-2 rotational excitation. For the 12C ground state and the first excited spin-2 state, we find a compact triangular configuration of alpha clusters. For the Hoyle state and the second excited spin-2 state, we find a "bent-arm" or obtuse triangular configuration of alpha clusters. We also calculate the electromagnetic transition rates between the low-lying states of 12C.

The Hoyle state plays a crucial role in the hydrogen burning of stars heavier than our sun and in the production of carbon and other elements necessary for life. This excited state of the carbon-12 nucleus was postulated by Hoyle [1] as a necessary ingredient for the fusion of three alpha particles to produce carbon at stellar temperatures. Although the Hoyle state was seen experimentally more than a half century ago [2,3], nuclear theorists have not yet uncovered the nature of this state from first principles. In this letter we report the first ab initio calculation of the low-lying states of carbon-12 using supercomputer lattice simulations and a theoretical framework known as effective field theory. In addition to the ground state and excited spin-2 state, we find a resonance at -85(3) MeV with all of the properties of the Hoyle state and in agreement with the experimentally observed energy. These lattice simulations provide insight into the structure of this unique state and new clues as to the amount of fine-tuning needed in nature for the production of carbon in stars.

We discuss lattice simulations of the ground state of dilute neutron matter at next-to-leading order in chiral effective field theory. In a previous paper the coefficients of the next-to-leading-order lattice action were determined by matching nucleon-nucleon scattering data for momenta up to the pion mass. Here the same lattice action is used to simulate the ground state of up to 12 neutrons in a periodic cube using Monte Carlo. We explore the density range from 2% to 8% of normal nuclear density and analyze the ground state energy as an expansion about the unitarity limit with corrections due to finite scattering length, effective range, and P-wave interactions.

We calculate the two-pion continuum contribution to the nucleon isovector spectral functions drawing upon the new high statistics measurements of the pion form factor by the CMD-2, KLOE, and SND collaborations. The general structure of the spectral functions remains unchanged, but the magnitude increases by about 10%. Using the updated spectral functions, we calculate the contribution of the two-pion continuum to the nucleon isovector form factors and radii. We compare the isovector radii with simple rho-pole models and illustrate their strong underestimation in such approaches. Moreover, we give a convenient parametrization of the result for use in future form factor analyses.