We calculate the $p$-wave contribution to the proton-proton fusion $S$ factor and its energy derivative in pionless effective field theory (EFT) up to next-to-leading order. The leading contributions are given by a recoil piece from the Gamow-Teller and Fermi operators, and from relativistic $1/m$ suppressed weak interaction operators. We obtain the value of $(2.5\pm0.3 )\times 10^{-28}~\mathrm{MeV\ fm^2}$ for the $S$ factor and $(2.2\pm0.2) \times 10^{-26}~\mathrm{fm^2}$ for its energy derivative at threshold. These are smaller than the results of a prior study that employed chiral EFT by several orders of magnitude. We conclude that, contrary to what has been previously reported, the $p$-wave contribution does not need to be considered in a high-precision determination of the $S$ factor at astrophysical energies. Combined with the chiral EFT calculation of Acharya {\it et al.} [Phys. Lett. B \bf{760}, 584 (2016)] for the $s$-wave channel, this gives a total threshold $S$ factor of $S(0) = (4.047^{+0.024}_{-0.032}) \times 10^{-23}~{\rm MeV~fm}^2$.

Bayesian analysis of the radiative capture reactions $^3\mathrm{He}(\alpha,\gamma)^7\mathrm{Be}$ and $^3\mathrm{H}(\alpha,\gamma)^7\mathrm{Li}$ are performed to draw inferences about the cross sections at threshold. We do a model comparison of two competing effective field theory power countings for the capture reactions. The two power countings differ in the contribution of two-body electromagnetic currents. In one power counting, two-body currents contribute at leading order, and in the other they contribute at higher orders. The former is favored for $^3\mathrm{He}(\alpha,\gamma)^7\mathrm{Be}$ if elastic scattering data in the incoming channel is considered in the analysis. Without constraints from elastic scattering data, both the power countings are equally favored. For $^3\mathrm{H}(\alpha,\gamma)^7\mathrm{Li}$, the first power counting with two-body current contributions at leading order is favored with or without constraints from elastic scattering data.

In this work we study the radiative capture of ${\rm {}^3He}$ on ${\rm {}^4He}$ within the halo effective field theory (EFT) framework. At leading order the capture amplitude comprises the initial state $s$-wave strong and Coulomb interactions summed to all orders. At the same order in the expansion, leading two-body currents contribute as well. We find delicate cancelations between the various contributions, and the two-body current contributions can be replaced by appropriately enhancing the asymptotic normalizations of the $^7$Be ground and first excited state wave functions. The next-to-leading order corrections come from the $s$-wave shape parameter and the pure Coulomb $d$-wave initial state interactions. We fit the EFT parameters to available scattering data and most recent capture data. Our zero-energy astrophysical $S$-factor estimate, $S_{34}\sim 0.55$ keV b, is consistent within error bars with the average in the literature.

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.

The r-mode instability, believed to limit the rotation speed of compact stars, can provide empirical confirmation for the existence of stable deconfined phases of quark matter that are predicted by weak coupling calculations in Quantum Chromodynamics. We construct a model for strange quark stars as heavy as 2 solar masses that are made of superconducting quark matter in the bulk and a thin crystalline quark matter crust. This crystalline quark crust is sufficiently robust to withstand r-mode heating and viscous rubbing for realistic mode amplitudes O(10^{-2}), unlike a crust made of neutron-rich nuclei. The dissipation provided by viscous rubbing at the core-crust boundary is both necessary and sufficient to obtain stable rotation speeds that are consistent with the majority of rapidly spinning pulsars in low mass X-ray binaries. Our analysis implies that while bare strange stars are ruled out by the existence of rapidly spinning pulsars, a strange star with a quark matter crust is a distinct possibility.

The cross section for radiative capture of neutron on carbon-14 is calculated using the model-independent formalism of halo effective field theory. The dominant contribution from E1 transition is considered, and the cross section is expressed in terms of elastic scattering parameters of the effective range expansion. Contributions from both resonant and non-resonant interaction are calculated. Significant interference between these leads to a capture contribution that deviates from simple Breit-Wigner resonance form.

Using halo effective field theory, we provide a model-independent calculation of the radiative neutron capture on lithium-7 over an energy range where the contribution from the 3+ resonance becomes important. One finds that a satisfactory description of the capture reaction, in the present single-particle approximation, suggests the use of a resonance width about three times larger than the experimental value. We also present power counting arguments that establish a hierarchy for electromagnetic one- and two-body currents.

The radiative neutron capture on lithium-7 is calculated model independently using a low energy halo effective field theory. The cross section is expressed in terms of scattering parameters directly related to the S-matrix element. The cross section depends on the poorly known p-wave effective range parameter r. This constitutes the leading order uncertainty in traditional model calculations. It is explicitly demonstrated by comparing with potential model calculations. A single parameter fit describes the low energy data extremely well and yields r = -1.47 fm^{-1}.

The r-mode instability in rotating compact stars is used to constrain the phase of matter at high density. The color-flavor-locked phase with kaon condensation (CFL-K0) and without (CFL) is considered in the temperature range 10^8K < T <10^{11} K. While the bulk viscosity in either phase is only effective at damping the r-mode at temperatures T > 10^{11} K, the shear viscosity in the CFL-K0 phase is the only effective damping agent all the way down to temperatures T > 10^8 K characteristic of cooling neutron stars. However, it cannot keep the star from becoming unstable to gravitational wave emission for rotation frequencies f ~ 56-11 Hz at T ~ 10^8-10^9 K. Stars composed almost entirely of CFL or CFL-K0 matter are ruled out by observation of rapidly rotating neutron stars, indicating that dissipation at the quark-hadron interface or nuclear crust interface must play a key role in damping the instability.

We determine characteristic timescales for the viscous damping of r-mode oscillations in rapidly rotating compact stars that contain quark matter. We present results for the color-flavor-locked (CFL) phase of dense quark matter, in which the up, down and strange quarks are gapped, as well as the normal (ungapped) quark phase. While the ungapped quark phase supports a temperature window between 10^8 K and 5x10^9 K where the r-mode is damped even for rapid rotation, the r-mode in a rapidly rotating pure CFL star is not damped in the temperature range 10^10 K - 10^11 K. Rotating hybrid stars with quark matter cores display an instability window whose width is determined by the amount of quark matter present, and they can have large spin frequencies outside this window. Except at high temperatures T > 10^10 K, the presence of a quark phase allows for larger critical frequencies and smaller spin-periods compared to rotating neutron stars. If low-mass X-ray binaries contain a large amount of ungapped or CFL quark matter, then our estimates of the r-mode instability suggest that there should be a population of rapidly rotating binaries at frequencies greater than 1000 Hz which have not yet been observed.

We analyze the free energy of charge and color neutral 2-flavor quark matter within the BCS approximation. We consider both the homogeneous gapless superconducting phase and the heterogeneous mixed phase where normal and BCS superconducting phases coexist. We calculate the surface tension between normal and superconducting phases and use it to compare the free energies of the gapless and mixed phases. Our calculation, which retains only the leading order gradient contribution to the free energy, indicates that the mixed phase is energetically favored over an interesting range of densities of relevance to 2 flavor quark matter in neutron stars.

At the high densities present in the interior of neutron stars, the neutrons are condensed into the 3P2 superfluid phase. While this condensation has little impact on the equation of state, it can have an important role in determining the low-temperature energy-momentum transport properties. The spontaneous breaking of baryon number by the condensate gives rise to the familiar Goldstone boson, but in addition, the spontaneous breaking of rotational invariance by the condensate gives rise to three Goldstone bosons, in general, one for each broken generator of rotations. These Goldstone bosons, which couple to the Z0, provide a new mechanism for neutrino emission. Using a low-energy effective field theory to describe the dynamics of these Goldstone bosons we estimate the neutrino emissivity of dense neutron matter and show that their annihilation is the dominant energy-loss mechanism over a range of temperatures.

An effective field theory calculation of the n+p->d+gamma cross section accurate to 1% for center of mass energy E up to ~1 MeV is presented. At these energies, which are relevant for big-bang nucleosynthesis, isovector magnetic transitions M1 and isovector electric transitions E1 give the dominant contributions. The M1 amplitude is calculated up to next-to-next-to-leading order (NNLO) and the contribution from the associated four-nucleon-one-photon operator is determined from the cold neutron capture rate. The E1 amplitude is calculated up to NNNNLO. The four-nucleon-one-photon operator contribution to E1 is determined from the related deuteron photodisintegration reaction gamma+d->n+p.