Search

Neutron stars are not observed to spin faster than about half their breakup rate. This limiting rotational frequency may be related to the strength of their crusts. As a star spins up from accretion, centrifugal forces stress the crust. We perform finite-element simulations of rotating neutron stars and find that the crust fails at rotation rates about half the breakup rate. Given uncertainties in microphysics, we have not determined the crust configuration after this failure. Instead, we argue that the crust may fail in an asymmetric way and could produce a configuration with a significant ellipticity (fractional difference in moments of inertia). If the ellipticity is large, a rotating star will radiate gravitational waves that may limit further spin up. These stars may be promising sources for LIGO / VIRGO and next generation gravitational wave detectors., Comment: 7 pages, 3 figures, 1 table. Additional discussion added including dependence on equation of state. ApJL in press

``Mountains'', or non-axisymmetrical deformations in the elastic crust of rotating neutron stars are efficient radiators of continuous gravitational waves. Recently, small anisotropies were observed in the solid innermost inner core of the Earth. We implement three-dimensional finite-element simulations to study mountains sourced by modest anisotropies in the solid crust of rotating neutron stars. We find that anisotropic mountains may be detectable by current ground-based gravitational-wave detectors and might explain several observed phenomena, confirming the results of a previous work on less realistic neutron star models. In particular, we find that a slightly anisotropic neutron star crust that changes its rotation rate modestly can support an ellipticity of a few $\times$ 10$^{-9}$, which is equivalent to the upper bounds on the ellipticity of some nearby and rapidlly spinning pulsars., Comment: 13 pages, 10 figures

"Mountains" or non-axisymmetric deformations of rotating neutron stars (NS) efficiently radiate gravitational waves (GW). We consider analogies between NS mountains and surface features of solar system bodies. Both NS and moons such as Europa or Enceladus have thin crusts over deep oceans while Mercury has a thin crust over a large metallic core. Thin sheets may wrinkle in universal ways. Europa has linear features, Enceladus has "Tiger" stripes, and Mercury has lobate scarps. NS may have analogous features. The innermost inner core of the Earth is anisotropic with a shear modulus that depends on direction. If NS crust material is also anisotropic this will produce an ellipticity, when the crust is stressed, that grows with spin frequency. This yields a braking index (log derivative of spin down rate assuming only GW spin down) very different from $n=5$ and could explain the maximum spin observed for neutron stars and a possible minimum ellipticity of millisecond pulsars., Comment: 8 pages, 2 figures, Added finite element simulations that confirm results, PRD in press

The recently published CREX results suggest a rather peculiar picture for the density dependence of the symmetry energy. Whereas PREX favors a large neutron skin thickness in $^{208}$Pb, thereby suggesting a stiff equation of state, CREX suggests instead a much softer equation of state. This discrepancy has caused a large spur in the theoretical community since no model has been able to simultaneously reproduce within $1\sigma$ the PREX and CREX results. Motivated by a novel correlation between a CREX observable and a combination of bulk symmetry energy parameters, we calibrate three new covariant energy density functionals that reproduce binding energies and charge radii of spherical nuclei - and also accommodate the constraints imposed by PREX and CREX. Given that these models suggest a stiff equation of state at high densities, predictions for neutron star properties are also discussed., Comment: 10 pages, resubmitted to PRC

We study the detectability of gravitational-wave signals from sub-solar mass binary neutron star systems by the current generation of ground-based gravitational-wave detectors. We find that finite size effects from large tidal deformabilities of the neutron stars and lower merger frequencies can significantly impact the sensitivity of the detectors to these sources. By simulating a matched-filter based search using injected binary neutron star signals with tidal deformabilities derived from physically motivated equations of state, we calculate the reduction in sensitivity of the detectors. We conclude that the loss in sensitive volume can be as high as $78.4 \%$ for an equal mass binary system of chirp mass $0.17 \, \textrm{M}_{\odot}$, in a search conducted using binary black hole template banks. We use this loss in sensitive volume, in combination with the results from the search for sub-solar mass binaries conducted on data collected by the LIGO-Virgo observatories during their first three observing runs, to obtain a conservative upper limit on the merger rate of sub-solar mass binary neutron stars. Since the discovery of a low-mass neutron star would provide new insight into formation mechanisms of neutron stars and further constrain the equation of state of dense nuclear matter, our result merits a dedicated search for sub-solar mass binary neutron star signals., Comment: 12 pages, 7 figures, supplemental materials at https://github.com/sugwg/sub-solar-ns-detectability

Non-axisymmetrical deformations of the crust on rapidly rotating neutron stars are one of the main targets of searches for continuous gravitational waves. The maximum ellipticity, or fractional difference in moments of inertia, that can be supported by deformations of the crust (known as "mountains") provides an important upper limit on the strength of these continuous gravitational wave sources. We use the formalism of Gittins et al 2021, along with a deforming force that acts mainly in the transverse direction, to obtain a maximum ellipticity of 7.4$\times$10$^{-6}$. This is larger than the original results of Gittins et al 2021 but consistent with earlier calculations by Ushomirsky et al 2000. This suggests that rotating neutron stars could be strong sources of continuous gravitational waves., Comment: 9 pages, 5 figures, minor revisions MNRAS in press

Type-Ia supernovae (SN Ia) are powerful stellar explosions that provide important distance indicators in cosmology. Recently, we proposed a new SN Ia mechanism that involves a nuclear fission chain reaction in an isolated white dwarf (WD) [PRL 126, 1311010]. The first solids that form as a WD starts to freeze are actinide rich and potentially support a fission chain reaction. In this letter we explore thermonuclear ignition from fission heating. We perform thermal diffusion simulations and find at high densities, above about 7x10^8 g/cm^3, that the fission heating can ignite carbon burning. This could produce a SN Ia or another kind of astrophysical transient., Comment: 6 pages, 3 figures, Accepted ApJ Letters

We report a precise measurement of the parity-violating asymmetry $A_{\rm PV}$ in the elastic scattering of longitudinally polarized electrons from $^{48}{\rm Ca}$. We measure $A_{\rm PV} =2668\pm 106\ {\rm (stat)}\pm 40\ {\rm (syst)}$ parts per billion, leading to an extraction of the neutral weak form factor $F_{\rm W} (q=0.8733$ fm$^{-1}) = 0.1304 \pm 0.0052 \ {\rm (stat)}\pm 0.0020\ {\rm (syst)}$ and the charge minus the weak form factor $F_{\rm ch} - F_{\rm W} = 0.0277\pm 0.0055$. The resulting neutron skin thickness $R_n-R_p=0.121 \pm 0.026\ {\rm (exp)} \pm 0.024\ {\rm (model)}$~fm is relatively thin yet consistent with many model calculations. The combined CREX and PREX results will have implications for future energy density functional calculations and on the density dependence of the symmetry energy of nuclear matter., Comment: 8 pages, 5 figures Replace 6-16-22: included ancillary files, corrected errors in references, author affiliations. Small text changes for clarity

A thorough understanding of neutrino-nucleus scattering physics is crucial for the successful execution of the entire US neutrino physics program. Neutrino-nucleus interaction constitutes one of the biggest systematic uncertainties in neutrino experiments - both at intermediate energies affecting long-baseline Deep Underground Neutrino Experiment (DUNE), as well as at low energies affecting coherent scattering neutrino program - and could well be the difference between achieving or missing discovery level precision. To this end, electron-nucleus scattering experiments provide vital information to test, assess and validate different nuclear models and event generators intended to be used in neutrino experiments. In this white paper, we highlight connections between electron- and neutrino-nucleus scattering physics at energies ranging from 10s of MeV to a few GeV, review the status of ongoing and planned electron scattering experiments, identify gaps, and layout a path forward that benefits the neutrino community. We also highlight the systemic challenges with respect to the divide between the nuclear and high-energy physics communities and funding that presents additional hurdle in mobilizing these connections to the benefit of neutrino programs., Comment: 37 pages, contribution to Snowmass 2021

Several semilocal exchange potentials usually employed in the framework of density-functional theory (DFT) are tested and compared with their exact counterpart, the exchange Optimized Effective Potential (OEP), as applied to the jellium-slab model of a metal-vacuum interface. Driven by their explicit dependence on the ground-state density, its gradient, and its kinetic-energy density, the three analyzed semilocal exchange potentials approach their respective asymptotic limits faster than in the case of the OEP, all of them having an asymptotic scaling of the form $-\alpha\,e^2/z + V_{\infty}$, with $\alpha < 1$. Here we provide the leading analytic asymptotics of the three model potentials under study, and we find that none of them exhibits the exact OEP slab asymptotics $-\;e^2/z$. While the so-called Becke-Roussel potential's leading asymptote is close to its exact OEP counterpart, the other two model potentials under study approach a material-dependent positive constant value far into the vacuum, resulting in considerably overestimated ionization potentials., Comment: To appear in Phys. Rev. B

Primordial black holes, with masses comparable to asteroids, are an attractive possibility for dark matter. In addition, other forms of dark matter could form compact dark objects (CDO). We search for small tidal accelerations from low mass black holes or CDOs orbiting near the Earth, and find none. Using about 10 years of data from the superconducting gravimeters in the Black Forest Observatory in South-Western Germany and at Djougou, Northern Benin in Western Africa we set an upper limit on the maximum mass of any dark object orbiting the Earth as a function of orbital radius. For semi-major axis less than two earth radii we exclude all black holes or CDOs with masses larger than 6.7x10^{13} kg. Lower mass primordial black holes may be strongly constrained by Hawking radiation. We conclude that near Earth black holes are extremely unlikely., Comment: six pages, 6 figures

We report the first measurement of the parity-violating elastic electron scattering asymmetry on 27Al. The 27Al elastic asymmetry is A_PV = 2.16 +- 0.11 (stat) +- 0.16 (syst) ppm, and was measured at =0.02357 +- 0.0001 GeV^2, = 7.61 +- 0.02 degrees, and = 1.157 GeV with the Qweak apparatus at Jefferson Lab. Predictions using a simple Born approximation as well as more sophisticated distorted-wave calculations are in good agreement with this result. From this asymmetry the 27Al neutron radius R_n = 2.89 +- 0.12 fm was determined using a many-models correlation technique. The corresponding neutron skin thickness R_n-R_p = -0.04 +- 0.12 fm is small, as expected for a light nucleus with a neutron excess of only 1. This result thus serves as a successful benchmark for electroweak determinations of neutron radii on heavier nuclei. A tree-level approach was used to extract the 27Al weak radius R_w = 3.00 +- 0.15 fm, and the weak skin thickness R_wk - R_ch = -0.04 +- 0.15 fm. The weak form factor at this Q^2 is F_wk = 0.39 +- 0.04., Comment: Revised version after referee's comments and suggestions. Some improvements and clarification to the text, no changes to the figures, tables, results or conclusions. 7 pages, 3 figures, 3 tables

The workshop Dawn VI: Next Generation Observatories took place online over three days, 5-7 October, 2021. More than 200 physicists and astronomers attended to contribute to, and learn from, a discussion of next-generation ground-based gravitational-wave detectors. The program was centered on the next generation of ground-based gravitational-wave observatories and their synergy with the greater landscape of scientific observatories of the 2030s. Cosmic Explorer (CE), a concept developed with US National Science Foundation support, was a particular focus; Einstein Telescope (ET), the European next generation concept, is an important complement and partner in forming a network. The concluding summary of the meeting expressed the sentiment that the observational science accessible to CE and ET, also in combination with data from other non-GW observatories, will stimulate a very broad community of analysts and yield insights which are exciting given the access to GWs from the entire universe. The need, and desire, for closer collaboration between ET and CE was expressed; a three-detector network is optimal for delivering much of the science. The science opportunities afforded by CE and ET are broad and compelling, impacting a wide range of disciplines in physics and high energy astrophysics. There was a consensus that CE is a concept that can deliver the promised science. A strong endorsement of Cosmic Explorer, as described in the CE Horizon Study, is a primary outcome of DAWN VI., Comment: Proceedings of the Dawn VI workshop

Type-Ia supernovae (SN Ia) are powerful stellar explosions that provide important distance indicators in cosmology. There is significant tension between values of the Hubble constant (expansion rate of the universe) determined from SN Ia and from other data. Recently, we proposed a new SN Ia mechanism that involves a nuclear fission chain reaction in an isolated white dwarf [PRL 126, 1311010]. We find the average mass of an exploding star decreases with increasing enrichment f_5 -- the fraction of uranium that is the isotope U-235. As a result, the average SN Ia luminosity decreases with increasing f_5. Furthermore, f_5 is likely higher in the host galaxies of SN Ia observed at large redshift $z$ because of younger galaxy ages. This change of f_5 leads to the evolution of SN Ia properties with redshift. If f_5 increases with redshift this results in an increased SN Ia rate, but a lower average SN Ia luminosity., Comment: 9 pages, 7 figures, submitted to ApJ

We report precision determinations of the beam normal single spin asymmetries ($A_n$) in the elastic scattering of 0.95 and 2.18~GeV electrons off $^{12}$C, $^{40}$Ca, $^{48}$Ca, and $^{208}$Pb at very forward angles where the most detailed theoretical calculations have been performed. The first measurements of $A_n$ for $^{40}$Ca and $^{48}$Ca are found to be similar to that of $^{12}$C, consistent with expectations thus demonstrating the validity of theoretical calculations for nuclei with Z~$\leq20$. We also report $A_n$ for $^{208}$Pb at two new momentum transfers (Q$^2$) extending the previous measurement. Our new data confirm the surprising result previously reported, with all three data points showing significant disagreement with the results from the $Z\leq 20$ nuclei. These data confirm our basic understanding of the underlying dynamics that govern $A_n$ for nuclei containing $\lesssim 50$ nucleons, but point to the need for further investigation to understand the unusual $A_n$ behaviour discovered for scattering off $^{208}$Pb., Comment: 7 pages, 2 figures

The asymptotics of the Kohn-Sham (KS) exact exchange potential $V_x(z)$ of a jelliumlike semi-infinite metal is investigated, in the framework of the optimized-effective-potential formalism of density-functional theory. Our numerical calculations clearly show that deep into the vacuum side of the surface $V_x(z) \propto e^2 \ln(az) / z$, with $a$ being a system-dependent constant, thus confirming the analytical calculations reported in Phys. Rev. B {\bf 81}, 121106(R) (2010). A criticism of this work published in Phys. Rev. B {\bf 85}, 115124 (2012) is shown to be incorrect. Our rigorous exchange-only results provide strong constraints both for the building of approximate exchange functionals and for the determination of the still unknown KS correlation potential., Comment: To appear in Phys. Rev. B

Type-Ia supernovae (SN Ia) are powerful stellar explosions that provide important distance indicators in cosmology. Recently, we proposed a new SN Ia mechanism that involves a nuclear fission chain-reaction in an isolated white dwarf [PRL 126, 1311010]. Here we perform novel reaction network simulations of the actinide-rich first solids in a cooling white dwarf. The network includes neutron-capture and fission reactions on a range of U and Th isotopes with various possible values for U-235 enrichment. We find, for modest U-235 enrichments, neutron-capture on U-238 and Th-232 can breed additional fissile nuclei so that a significant fraction of all U and Th nuclei may fission during the chain-reaction. The resulting large energy release could ignite thermonuclear carbon burning and possibly trigger a SN Ia., Comment: 8 pages, 5 figures, revised introduction, accepted Phys. Rev. C

Do white dwarfs have inner cores made of iron? Neutron rich nuclei like $^{56}$Fe experience a net gravitational force and sediment toward the core. Using new phase diagrams and molecular dynamics simulations, we show that $^{56}$Fe should separate into mesoscopic Fe-rich crystallites due to its large charge relative to the background. At solar abundances, these crystallites rapidly precipitate and form an inner core of order 100 km and $10^{-3} M_\odot$ that may be detectable with asteroseismology. Associated cooling delays could be up to a Gyr for low mass white dwarfs but are only $\sim$0.1 Gyr for massive white dwarfs, so while this mechanism may contribute to the Q-branch the heating is insufficient to fully explain it., Comment: 8 pages, 5 figures, accepted for publication in ApJ Letters

The first solids that form as a white dwarf (WD) starts to crystallize are expected to be greatly enriched in actinides. Previously [PRL 126, 1311010] we found that these solids might support a nuclear fission chain reaction that could ignite carbon burning and provide a new Type Ia supernova (SN Ia) mechanism involving an {\it isolated} WD. Here we explore this fission mechanism in more detail and calculate the final temperature and density after the chain reaction and discuss a number of open physics questions., Comment: 8 pages, 4 figures

Searches for continuous gravitational waves from \textit{unknown} Galactic neutron stars provide limits on the shapes of neutron stars. A rotating neutron star will produce gravitational waves if asymmetric deformations exist in its structure that are characterized by the star's ellipticity. In this study, we use a simple model of the spatial and spin distribution of Galactic neutron stars to estimate the total number of neutron stars probed, using gravitational waves, to a given upper limit on the ellipticity. This may help optimize future searches with improved sensitivity. The improved sensitivity of third-generation gravitational wave detectors may increase the number of neutron stars probed, to a given ellipticity, by factors of 100 to 1000., Comment: Accepted by ApJ

The first solids that form as a cooling white dwarf (WD) starts to crystallize are expected to be greatly enriched in actinides. This is because the melting points of WD matter scale as $Z^{5/3}$ and actinides have the largest charge $Z$. We estimate that the solids may be so enriched in actinides that they could support a fission chain reaction. This reaction could ignite carbon burning and lead to the explosion of an isolated WD in a thermonuclear supernova (SN Ia). Our mechanism could potentially explain SN Ia with sub-Chandrasekhar ejecta masses and short delay times., Comment: 8 pages, 6 figures total including Appendix, Phys. Rev. Let. in press

We report a precision measurement of the parity-violating asymmetry $A_{PV}$ in the elastic scattering of longitudinally polarized electrons from $^{208}$Pb. We measure $A_{PV}=550\pm 16 {\rm (stat)}\pm 8\ {\rm (syst)}$ parts per billion, leading to an extraction of the neutral weak form factor $F_W(Q^2 = 0.00616\ {\rm GeV}^2) = 0.368 \pm 0.013$. Combined with our previous measurement, the extracted neutron skin thickness is $R_n-R_p=0.283 \pm 0.071$~fm. The result also yields the first significant direct measurement of the interior weak density of $^{208}$Pb: $\rho^0_W = -0.0796\pm0.0036\ {\rm (exp.)}\pm0.0013\ {\rm (theo.)}\ {\rm fm}^{-3}$ leading to the interior baryon density $\rho^0_b = 0.1480\pm0.0036\ {\rm (exp.)}\pm0.0013\ {\rm (theo.)}\ {\rm fm}^{-3}$. The measurement accurately constrains the density dependence of the symmetry energy of nuclear matter near saturation density, with implications for the size and composition of neutron stars.

Laboratory experiments sensitive to the equation of state of neutron rich matter in the vicinity of nuclear saturation density provide the first rung in a "density ladder" that connects terrestrial experiments to astronomical observations. In this context, the neutron skin thickness of 208Pb (Rskin) provides a stringent laboratory constraint on the density dependence of the symmetry energy. In turn, an improved value of Rskin has been reported recently by the PREX collaboration. Exploiting the strong correlation between Rskin and the slope of the symmetry energy L within a specific class of relativistic energy density functionals, we report a value of L=(106 +/- 37)MeV -- that systematically overestimates current limits based on both theoretical approaches and experimental measurements. The impact of such a stiff symmetry energy on some critical neutron-star observables is also examined., Comment: 5 pages, 4 figures. Manuscript accepted for publication in the Physical Review Letters

Recent observations of Galactic white dwarfs (WDs) with Gaia suggest there is a population of massive crystallizing WDs exhibiting anomalous cooling -- the Q branch. While single-particle $^{22}$Ne sedimentation has long been considered a possible heat source, recent work suggests that $^{22}$Ne must separate into clusters, enhancing diffusion, in order for sedimentation to provide heating on the observed timescale. We show definitively that $^{22}$Ne cannot separate to form clusters in C/O WDs using molecular dynamics simulations, and we further present a general C/O/Ne phase diagram showing that strong $^{22}$Ne enrichment is not achievable for $^{22}$Ne abundance $\lesssim 30\%$. We conclude that the anomalous heating cannot be due to $^{22}$Ne cluster sedimentation and that Q branch WDs may have an unusual composition, possibly rich with heavier elements., Comment: 7 pages, 4 figures, accepted for publication in ApJ Letters

The present PREX-II and CREX experiments are measuring the rms radius of the weak charge density of $^{208}$Pb and $^{48}$Ca. We discuss the feasibility of a new parity violating electron scattering experiment to measure the surface thickness of the weak charge density of a heavy nucleus. Once PREX-II and CREX have constrained weak radii, an additional parity violating measurement at a momentum transfer near 0.76 fm$^{-1}$ for $^{208}$Pb or 1.28 fm$^{-1}$ for $^{48}$Ca can determine the surface thickness., Comment: 8 pages, 10 figures

The energy radiated in supernova neutrinos is a fundamental quantity that is closely related to the gravitational binding energy of a neutron star. Recently the tidal deformability of neutron stars was constrained by gravitational wave observations. By considering several equations of state, we find a strong correlation between the tidal deformability and neutron star binding energy. We use this correlation to sharpen predictions of the binding energy of neutron stars and the total neutrino energy in supernovae. We find a minimum binding energy for a neutron star formed in a supernova of $\sim1.5\times 10^{53}$ ergs. Should the neutrino energy in a supernova be significantly below this value, it would strongly suggest new unobserved particles are carrying away some of the supernova energy. Alternatively, if the neutrino energy is observed above $\sim 6\times 10^{53}$ ergs, it would strongly imply the formation of a (perhaps surprisingly) massive neutron star., Comment: 8 pages, 6 figures

Recently, Cheng et al. identified a number of massive white dwarfs (WD) that appear to have an additional heat source providing a luminosity near $\approx 10^{-3}L_\odot$ for multiple Gyr. In this paper we explore heating from electron capture and pycnonuclear reactions. We also explore heating from dark matter annihilation. WD stars appear to be too small to capture enough dark matter for this to be important. Finally, if dark matter condenses to very high densities inside a WD this could ignite nuclear reactions. We calculate the enhanced central density of a WD in the gravitational potential of a very dense dark matter core. While this might start a supernova, it seems unlikely to provide modest heating for a long time. We conclude that electron capture, pycnonuclear, and dark matter reactions are unlikely to provide significant heating in the massive WD that Cheng considers., Comment: 5 pages, 3 figures

The dependence of fusion dynamics on neutron excess for light nuclei is extracted. This is accomplished by comparing the average fusion cross-section at energies just above the fusion barrier for $^{12-15}$C + $^{12}$C with measurements of the interaction cross-section from high evergy collisions. The experimental results indicate that the fusion cross-section associated with dynamics increases with increasing neutron excess. Calculations with a time-dependent Hartree-Fock model fail to describe the observed trend., Comment: 5 pqages, 4 figures. arXiv admin note: substantial text overlap with arXiv:2005.03553

The saturation density of nuclear matter $\rho_0$ is a fundamental nuclear physics property that is difficult to predict from fundamental principles. The saturation density is closely related to the interior density of a heavy nucleus, such as $^{208}$Pb. We use parity violating electron scattering to determine the average interior weak charge and baryon densities in $^{208}$Pb. This requires not only measuring the weak radius $R_{\rm wk}$ but also determining the surface thickness of the weak charge density $a$. We obtain $\rho_0=0.150\pm0.010$ fm$^{-3}$, where the 7\% error has contributions form the PREX error on the weak radius, an assumed 10\% uncertainty in the surface thickness $a$, and from the extrapolation to infinite nuclear matter. These errors can be improved with the upcoming PREX II results and with a new parity violating electron scattering experiment, at a somewhat higher momentum transfer, to determine $a$., Comment: 6 pages, 3 figures, minor changes, Phys. Rev. C in press

Is the secondary component of GW190814 the lightest black hole or the heaviest neutron star ever discovered in a double compact-object system [R. Abbott et al., ApJ Lett., 896, L44 (2020)]? This is the central question animating this letter. Covariant density functional theory provides a unique framework to investigate both the properties of finite nuclei and neutron stars, while enforcing causality at all densities. By tuning existing energy density functionals we were able to: (a) account for a 2.6 Msun neutron star, (b) satisfy the original constraint on the tidal deformability of a 1.4 Msun neutron star, and (c) reproduce ground-state properties of finite nuclei. Yet, for the class of models explored in this work, we find that the stiffening of the equation of state required to support super-massive neutron stars is inconsistent with either constraints obtained from energetic heavy-ion collisions or from the low deformability of medium-mass stars. Thus, we speculate that the maximum neutron star mass can not be significantly higher than the existing observational limit and that the 2.6 Msun compact object is likely to be the lightest black hole ever discovered., Comment: 6 pages and 3 figures

The dynamics present in the fusion of neutron-rich nuclei is explored through the comparison of experimental cross-sections at above-barrier energies with measurements of the interaction cross-section at relativistic energies. The increase of fusion dynamics with increasing neutron excess is clearly demonstrated. Experimental cross-sections are compared with the predictions of a Sao Paulo model using relativistic mean field density distributions and the impact of different interactions is explored., Comment: 5 pages, 5 figures

Above-barrier fusion cross-sections for an isotopic chain of oxygen isotopes with A=16-19 incident on a $^{12}$C target are presented. Experimental data are compared with both static and dynamical microscopic calculations. These calculations are unable to explain the $\sim$37\% increase in the average above-barrier fusion cross-section observed for $^{19}$O as compared to $\beta$-stable oxygen isotopes. This result suggests that for neutron-rich nuclei existing time-dependent Hartree-Fock calculations underpredict the role of dynamics at near-barrier energies. High-quality measurement of above-barrier fusion for an isotopic chain of increasingly neutron-rich nuclei provides an effective means to probe this fusion dynamics., Comment: 6 pages, 4 figures

Dark matter could be composed of compact dark objects (CDOs). These objects may interact very weakly with normal matter and could move freely {\it inside} the Earth. A CDO moving in the inner core of the Earth will have an orbital period near 55 min and produce a time dependent signal in a gravimeter. Data from superconducting gravimeters rule out such objects moving inside the Earth unless their mass $m_D$ and or orbital radius $a$ are very small so that $m_D\, a < 1.2\times 10^{-13}M_\oplus R_\oplus$. Here $M_\oplus$ and $R_\oplus$ are the mass and radius of the Earth., Comment: 8 pages, 7 figures

The historic observations of the neutron star merger GW170817 advanced our understanding of r-process nucleosynthesis and the equation of state (EOS) of neutron rich matter. Simple neutrino physics suggests that supernovae are not the site of the main r-process. Instead, the very red color of the kilonova associated with GW170817 shows that neutron star (NS) mergers are an important r-process site. We now need to measure the masses and beta decay half-lives of very neutron rich heavy nuclei so that we can more accurately predict the abundances of heavy elements that are produced. This can be done with new radioactive beam accelerators such as the Facility for Rare Isotope Beams (FRIB). GW170817 provided information on the deformability of NS and the equation of state of dense matter. The PREX II experiment will measure the neutron skin of ${}^{208}$Pb and help constrain the low density EOS. As the sensitivity of gravitational wave detectors improve, we expect to observe many more events. We look forward to exciting advances and surprises!, Comment: 12 pages, 4 figures

The historic first detection of the binary neutron star merger GW170817 by the LIGO-Virgo collaboration has set a limit on the gravitational deformability of neutron stars. In contrast, radio observations of PSR J0740+6620 find a very massive neutron star. Tension between the small deformability and the large maximum mass may suggest that the pressure rises rapidly with density and thus the speed of sound in dense matter is likely a large fraction of the speed of light. We use these observations and simple constant sound-speed model equations of state to set a lower bound on the maximum speed of sound in neutron stars. If the tidal deformability of a 1.4$M_\odot$ neutron star is less than 600, as is suggested by subsequent analyses of GW170817, we find that the sound speed in the cores of neutron stars is likely larger than the conformal limit of $c/\sqrt{3}$. Implications of this for our understanding of both hadronic and quark-gluon descriptions of dense matter are discussed., Comment: 5 pages, 3 figures

The complex nature of the nuclear forces generates a broad range and diversity of observational phenomena. Heavy nuclei, though orders of magnitude less massive than neutron stars, are governed by the same underlying physics, which is enshrined in the nuclear equation of state. Heavy nuclei are expected to develop a neutron-rich skin where many neutrons collect near the surface. Such a skin thickness is strongly sensitive to the poorly-known density dependence of the symmetry energy near saturation density. An accurate and model-independent determination of the neutron-skin thickness of heavy nuclei would provide a significant first constraint on the density dependence of the nuclear symmetry energy. The determination of the neutron-skin thickness of heavy nuclei has far reaching consequences in many areas of physics as diverse as heavy-ion collisions, polarized electron and proton scattering off nuclei, precision tests of the standard model using atomic parity violation, and nuclear astrophysics. While a systematic and concerted experimental effort has been made to measure the neutron-skin thickness of heavy nuclei, a precise and model-independent determination remains elusive. How to move forward at a time when many new facilities are being commissioned and how to strengthen the synergy with other areas of physics are primary goals of this review., Comment: Submitted for publication to Journal of Physics G

Dark matter could be composed of compact dark objects (CDOs). A close binary of CDOs orbiting in the interior of solar system bodies can be a loud source of gravitational waves (GWs) for the LIGO and VIRGO detectors. We perform the first search ever for this type of signal and rule out close binaries, with separations of order 300 m, orbiting near the center of the Sun with GW frequencies (twice the orbital frequency) between 50 and 550 Hz and CDO masses above $\approx 10^{-9} M_\odot$. This mass limit is eight orders of magnitude lower than the mass probed in a LIGO search at extra galactic distances., Comment: 5 pages, 3 figures, Phys. Lett. B in press

Dark matter could be composed of compact dark objects (CDOs). We find that the oscillation of CDOs inside neutron stars can be a detectable source of gravitational waves (GWs). The GW strain amplitude depends on the mass of the CDO, and its frequency is typically in the range 3-5 kHz as determined by the central density of the star. In the best cases, LIGO may be sensitive to CDO masses greater than or of order $10^{-8}$ solar masses., Comment: 5 pages, 1 figure, Phys. Rev. Lett. in press

The weak charge of the proton determines its coupling to the $Z^0$ boson. The distribution of weak charge is found to be dramatically different from the distribution of electric charge. The proton's weak radius $R_W= 1.545\pm 0.017$ fm is 80\% larger than its charge radius $R_{ch}\approx 0.84$ fm because of a very large pion cloud contribution. This large weak radius can be measured with parity violating electron scattering and may provide insight into the structure of the proton, various radiative corrections, and possible strange quark contributions., Comment: 4 pages, 2 figures, Error estimate improved, Phys. Lett. B in press

In July 2018 an FRIB Theory Alliance program was held on the implications of GW170817 and its associated kilonova for r-process nucleosynthesis. Topics of discussion included the astrophysical and nuclear physics uncertainties in the interpretation of the GW170817 kilonova, what we can learn about the astrophysical site or sites of the r process from this event, and the advances in nuclear experiment and theory most crucial to pursue in light of the new data. Here we compile a selection of scientific contributions to the workshop, broadly representative of progress in r-process studies since the GW170817 event., Comment: Proceedings for the FRIB Theory Alliance workshop "FRIB and the GW170817 kilonova", held 16-27 July 2018 at the Facility for Rare Isotope Beams, Michigan State University, East Lansing, MI USA; 44 pages, 14 figures

The elastic properties of neutron star crusts are relevant for a variety of currently observable or near-future electromagnetic and gravitational wave phenomena. These phenomena may depend on the elastic properties of nuclear pasta found in the inner crust. We present large scale classical molecular dynamics simulations where we deform nuclear pasta. We simulate idealized samples of nuclear pasta and describe their breaking mechanism. We also deform nuclear pasta that is arranged into many domains, similar to what is known for the ions in neutron star crusts. Our results show that nuclear pasta may be the strongest known material, perhaps with a shear modulus of $10^{30}\,\text{ergs/cm}^3$ and breaking strain greater than 0.1., Comment: 5 pages, 2 figures. Submitted to Physical Review Letters

This is an exciting time for the study of r-process nucleosynthesis. Recently, a neutron star merger GW170817 was observed in extraordinary detail with gravitational waves and electromagnetic radiation from radio to gamma rays. The very red color of the associated kilonova suggests that neutron star mergers are an important r-process site. Astrophysical simulations of neutron star mergers and core collapse supernovae are making rapid progress. Detection of both, electron neutrinos and antineutrinos from the next galactic supernova will constrain the composition of neutrino-driven winds and provide unique nucleosynthesis information. Finally FRIB and other rare-isotope beam facilities will soon have dramatic new capabilities to synthesize many neutron-rich nuclei that are involved in the r-process. The new capabilities can significantly improve our understanding of the r-process and likely resolve one of the main outstanding problems in classical nuclear astrophysics. However, to make best use of the new experimental capabilities and to fully interpret the results, a great deal of infrastructure is needed in many related areas of astrophysics, astronomy, and nuclear theory. We will place these experiments in context by discussing astrophysical simulations and observations of r-process sites, observations of stellar abundances, galactic chemical evolution, and nuclear theory for the structure and reactions of very neutron-rich nuclei. This review paper was initiated at a three-week International Collaborations in Nuclear Theory program in June 2016 where we explored promising r-process experiments and discussed their likely impact, and their astrophysical, astronomical, and nuclear theory context., Comment: 132 pages, 18 figures, submitted to Journal of Physics G

The crust of accreting neutron stars plays a central role in many different observational phenomena. In these stars, heavy elements produced by H-He burning in the rapid proton capture (rp-) process continually freeze to form new crust. In this paper, we explore the expected composition of the solid phase. We first demonstrate using molecular dynamics that two distinct types of chemical separation occur, depending on the composition of the rp-process ashes. We then calculate phase diagrams for three-component mixtures and use them to determine the allowed crust compositions. We show that, for the large range of atomic numbers produced in the rp-process ($Z\sim 10$--$50$), the solid that forms has only a small number of available compositions. We conclude that accreting neutron star crusts should be polycrystalline, with domains of distinct composition. Our results motivate further work on the size of the compositional domains, and have implications for crust physics and accreting neutron star phenomenology., Comment: 8 pages, 4 figures, Submitted to ApJ, this article supersedes arXiv:1709.09260

The limiting rotational frequency of neutron stars may be determined by the strength of their crusts. As a star spins up from accretion, centrifugal forces will cause the crust to fail. If the crust breaks unevenly, a rotating mass quadrupole moment will radiate gravitational waves (GW). This radiation can prevent further spin up and may be a promising source for continuous GW searches. We calculate that for a breaking strain (strength) of neutron star crust that is consistent with molecular dynamics simulations, the crust may fail at rotational frequencies in agreement with observations., Comment: 5 pages, 4 figures, 2 tables

We show that the neutron star in the transient system MXB~1659-29 has a core neutrino luminosity that substantially exceeds that of the modified Urca reactions (i.e., $n+n\to n+p+e^{-}+\bar{\nu}_{e}$ and inverse) and is consistent with the direct Urca reactions ($n\to p+e^{-}+\bar{\nu}_{e}$ and inverse) occurring in a small fraction of the core. Observations of the thermal relaxation of the neutron star crust following 2.5 years of accretion allow us to measure the energy deposited into the core during accretion, which is then reradiated as neutrinos, and infer the core temperature. For a nucleonic core, this requires that the nucleons are unpaired and that the proton fraction exceed a critical value to allow the direct Urca reaction to proceed. The neutron star in MXB~1659-29 is the first with a firmly detected thermal component in its X-ray spectrum that needs a fast neutrino cooling process. Measurements of the temperature variation of the neutron star core during quiescence would place an upper limit on the core specific heat and serve as a check on the fraction of the neutron star core in which nucleons are unpaired., Comment: 5 pages, 3 figures

Core-collapse supernova explosions are driven by a central engine that converts a small fraction of the gravitational binding energy released during core collapse to outgoing kinetic energy. The suspected mode for this energy conversion is the neutrino mechanism, where a fraction of the neutrinos emitted from the newly formed protoneutron star are absorbed by and heat the matter behind the supernova shock. Accurate neutrino-matter interaction terms are crucial for simulating these explosions. In this proceedings for IAUS 331, SN 1987A, 30 years later, we explore several corrections to the neutrino-nucleon scattering opacity and demonstrate the effect on the dynamics of the core-collapse supernova central engine via two dimensional neutrino-radiation-hydrodynamics simulations. Our results reveal that the explosion properties are sensitive to corrections to the neutral-current scattering cross section at the 10-20% level, but only for densities at or above $\sim 10^{12}$ g cm$^{-3}$, Comment: 6 pages, 3 figures, appears in Proc. IAU Symposium 331, SN 1987A, 30 years later - Cosmic Rays and Nuclei from Supernovae and Their Aftermaths

The historical first detection of a binary neutron star merger by the LIGO-Virgo collaboration [B. P. Abbott et al. Phys. Rev. Lett. 119, 161101 (2017)] is providing fundamental new insights into the astrophysical site for the $r$-process and on the nature of dense matter. A set of realistic models of the equation of state (EOS) that yield an accurate description of the properties of finite nuclei, support neutron stars of two solar masses, and provide a Lorentz covariant extrapolation to dense matter are used to confront its predictions against tidal polarizabilities extracted from the gravitational-wave data. Given the sensitivity of the gravitational-wave signal to the underlying EOS, limits on the tidal polarizability inferred from the observation translate into constraints on the neutron-star radius. Based on these constraints, models that predict a stiff symmetry energy, and thus large stellar radii, can be ruled out. Indeed, we deduce an upper limit on the radius of a $1.4\,M_{\odot}$ neutron star of $R_{\star}^{1.4}\!<\!13.76\,{\rm km}$. Given the sensitivity of the neutron-skin thickness of ${}^{208}$Pb to the symmetry energy, albeit at a lower density, we infer a corresponding upper limit of about $R_{\rm skin}^{208}\!\lesssim\!0.25\,{\rm fm}$. However, if the upcoming PREX-II experiment measures a significantly thicker skin, this may be evidence of a softening of the symmetry energy at high densities---likely indicative of a phase transition in the interior of neutron stars., Comment: 6 pages, 4 figures, edited figures and the text with a more correct interpretation of the LIGO-Virgo results

We present a new mechanism for deep crustal heating in accreting neutron stars. Charged pions ($\pi^+$) are produced in nuclear collisions on the neutron star surface during active accretion and upon decay they provide a flux of neutrinos into the neutron star crust. For massive and/or compact neutron stars, neutrinos deposit $\approx 1\textrm{--} 2 \, \mathrm{MeV}$ of heat per accreted nucleon into the inner crust. The strength of neutrino heating is comparable to the previously known sources of deep crustal heating, such as from pycnonuclear fusion reactions, and is relevant for studies of cooling neutron stars. We model the thermal evolution of a transient neutron star in a low-mass X-ray binary, and in the particular case of the neutron star MXB~1659-29 we show that additional deep crustal heating requires a higher thermal conductivity for the neutron star inner crust. A better knowledge of pion production cross sections near threshold would improve the accuracy of our predictions., Comment: 12 pages, 9 figures, 3 tables; [Added a new figure and edited the text in response to Referee's remarks and suggestions]

Nucleosynthesis in the oceans of accreting neutron stars can produce novel mixtures of nuclides, whose composition is dependent on the exact astrophysical conditions. Many simulations have now been done to determine the nucleosynthesis products in the ocean, but the phase separation at the base of the ocean, which determines the composition of the crust, has not been as well studied. In this work, we simulate the phase separation of a composition, which was predicted to produce a crust enriched in light nuclei, in contrast with past work which predicts that crust is enriched in heavy nuclei. We perform molecular dynamics simulations of the phase separation of this mixture using the methods of Horowitz $\textit{et. al.}$ (2007). We find good agreement with the predictions of Mckinven $\textit{et al.}$ (2016) for the phase separation of this mixture. Moreover, this supports their method as a computationally efficient alternative to molecular dynamics for calculating phase separation for a wider regime of astrophysical conditions., Comment: 8 pages, 5 figures, Submitted to Phys. Rev. E