Search

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., 8 pages, 4 figures, Submitted to ApJ, this article supersedes arXiv:1709.09260

A striking shape was recently observed for the endoplasmic reticulum, a cellular organelle consisting of stacked sheets connected by helical ramps [Terasaki et al., Cell 154, 285 (2013)]. This shape is interesting both for its biological function, to synthesize proteins using an increased surface area for ribosome factories, and its geometric properties that may be insensitive to details of the microscopic interactions. In the present work, we find very similar shapes in our molecular dynamics simulations of the nuclear pasta phases of dense nuclear matter that are expected deep in the crust of neutron stars. There are dramatic differences between nuclear pasta and terrestrial cell biology. Nuclear pasta is 14 orders of magnitude denser than the aqueous environs of the cell nucleus and involves strong interactions between protons and neutrons, while cellular-scale biology is dominated by the entropy of water and complex assemblies of biomolecules. Nonetheless, the very similar geometry suggests both systems may have similar coarse-grained dynamics and that the shapes are indeed determined by geometrical considerations, independent of microscopic details. Many of our simulations self-assemble into flat sheets connected by helical ramps. These ramps may impact the thermal and electrical conductivities, viscosity, shear modulus, and breaking strain of neutron star crust. The interaction we use, with Coulomb frustration, may provide a simple model system that reproduces many biologically important shapes.

Background: Exotic non-spherical nuclear pasta shapes are expected in nuclear matter at just below saturation density because of competition between short range nuclear attraction and long range Coulomb repulsion. Purpose: We explore the impact of nuclear pasta on nucleosynthesis, during neutron star mergers, as cold dense nuclear matter is ejected and decompressed. Methods: We perform classical molecular dynamics simulations with 51200 and 409600 nucleons, that are run on GPUs. We expand our simulation region to decompress systems from an initial density of 0.080 fm^{-3} down to 0.00125 fm^{-3}. We study proton fractions of Y_P=0.05, 0.10, 0.20, 0.30, and 0.40 at T =0.5, 0.75, and 1.0 MeV. We calculate the composition of the resulting systems using a cluster algorithm. Results: We find final compositions that are in good agreement with nuclear statistical equilibrium models for temperatures of 0.75 and 1 MeV. However, for proton fractions greater than Y_P=0.2 at a temperature of T = 0.5 MeV, the MD simulations produce non-equilibrium results with large rod-like nuclei. Conclusions: Our MD model is valid at higher densities than simple nuclear statistical equilibrium models and may help determine the initial temperatures and proton fractions of matter ejected in mergers., 13 pages

Nuclear pasta, with non-spherical shapes, is expected near the base of the crust in neutron stars. Large scale molecular dynamics simulations of pasta show long lived topological defects that could increase electron scattering and reduce both the thermal and electrical conductivities. We model a possible low conductivity pasta layer by increasing an impurity parameter Q_{imp}. Predictions of light curves for the low mass X-ray binary MXB 1659-29, assuming a large Q_{imp}, find continued late time cooling that is consistent with Chandra observations. The electrical and thermal conductivities are likely related. Therefore observations of late time crust cooling can provide insight on the electrical conductivity and the possible decay of neutron star magnetic fields (assuming these are supported by currents in the crust)., 5 pages, 2 figures

The dense neutron-rich matter found in supernovae and neutron stars is expected to form complex nonuniform phases referred to as nuclear pasta. The pasta shapes depend on density, temperature and proton fraction and determine many transport properties in supernovae and neutron star crusts. We use two recently developed hybrid CPU/GPU codes to perform large scale molecular dynamics (MD) simulations with $51200$ and $409600$ nucleons of nuclear pasta. From the output of the MD simulations we characterize the topology and compute two observables, the radial distribution function $g(r)$ and the structure factor $S(q)$, for systems with proton fractions $Y_p=0.10, 0.20, 0.30$ and $0.40$ at about one third of nuclear saturation density and temperatures near $1.0$ MeV. We observe that the two lowest proton fraction systems simulated, $Y_p=0.10$ and $0.20$, equilibrate quickly and form liquid-like structures. Meanwhile, the two higher proton fraction systems, $Y_p=0.30$ and $0.40$, take a longer time to equilibrate and organize themselves in solid-like periodic structures. Furthermore, the $Y_p=0.40$ system is made up of slabs, lasagna phase, interconnected by defects while the $Y_p=0.30$ systems consist of a stack of perforated plates, the nuclear waffle phase. The periodic configurations observed in our MD simulations for proton fractions $Y_p\ge0.30$ have important consequences for the structure factors $S(q)$ of protons and neutrons, which relate to many transport properties of supernovae and neutron star crust. A detailed study of the waffle phase and how its structure depends on temperature, size of the simulation and the screening length showed that finite-size effects appear to be under control and, also, that the plates in the waffle phase merge at temperatures slightly above $1.0$ MeV and the holes in the plates form an hexagonal lattice at temperatures slightly lower than $1.0$ MeV., Comment: 16 pages, 12 figires, 6 tables, small changes with respect to previous version, Phys Rev C in press

Using molecular dynamics simulations, we determine the structure of neutron star crust made of rapid proton capture nucleosynthesis material. We find a regular body centered cubic lattice, even with the large number of impurities that are present. Low charge $Z$ impurities tend to occupy interstitial positions, while high $Z$ impurities tend to occupy substitutional lattice sites. We find strong attractive correlations between low $Z$ impurities that could significantly increase the rate of pycnonuclear (density driven) nuclear reactions. The thermal conductivity is significantly reduced by electron impurity scattering. Our results will be used in future work to study the effects of impurities on mechanical properties such as the shear modulus and breaking strain., 7 pages, 9 figures, added A values to table I, version as published

Fusion reactions in the crust of an accreting neutron star are an important source of heat, and the depth at which these reactions occur is important for determining the temperature profile of the star. Fusion reactions depend strongly on the nuclear charge $Z$. Nuclei with $Z\le 6$ can fuse at low densities in a liquid ocean. However, nuclei with Z=8 or 10 may not burn until higher densities where the crust is solid and electron capture has made the nuclei neutron rich. We calculate the $S$ factor for fusion reactions of neutron rich nuclei including $^{24}$O + $^{24}$O and $^{28}$Ne + $^{28}$Ne. We use a simple barrier penetration model. The $S$ factor could be further enhanced by dynamical effects involving the neutron rich skin. This possible enhancement in $S$ should be studied in the laboratory with neutron rich radioactive beams. We model the structure of the crust with molecular dynamics simulations. We find that the crust of accreting neutron stars may contain micro-crystals or regions of phase separation. Nevertheless, the screening factors that we determine for the enhancement of the rate of thermonuclear reactions are insensitive to these features. Finally, we calculate the rate of thermonuclear $^{24}$O + $^{24}$O fusion and find that $^{24}$O should burn at densities near $10^{11}$ g/cm$^3$. The energy released from this and similar reactions may be important for the temperature profile of the star., Comment: 7 pages, 4 figs, minor changes, to be published in Phys. Rev. C

Recently, crust cooling times have been measured for neutron stars after extended outbursts. These observations are very sensitive to the thermal conductivity $\kappa$ of the crust and strongly suggest that $\kappa$ is large. We perform molecular dynamics simulations of the structure of the crust of an accreting neutron star using a complex composition that includes many impurities. The composition comes from simulations of rapid proton capture nucleosynthesys followed by electron captures. We find that the thermal conductivity is reduced by impurity scattering. In addition, we find phase separation. Some impurities with low atomic number $Z$ are concentrated in a subregion of the simulation volume. For our composition, the solid crust must separate into regions of different compositions. This could lead to an asymmetric star with a quadrupole deformation that radiates gravitational waves. Observations of crust cooling can constrain impurity concentrations., Comment: 8 pages, 8 figures, minor changes, PRE in press

We calculate the shear viscosity $\eta$ and thermal conductivity $\kappa$ of a nuclear pasta phase in neutron star crusts. This involves complex non-spherical shapes. We use semiclassical molecular dynamics simulations involving 40,000 to 100,000 nucleons. The viscosity $\eta$ can be simply expressed in terms of the height $Z^*$ and width $\Delta q$ of the peak in the static structure factor $S_p(q)$. We find that $\eta$ increases somewhat, compared to a lower density phase involving spherical nuclei, because $Z^*$ decreases from form factor and ion screening effects. However, we do not find a dramatic increase in $\eta$ from non-spherical shapes, as may occur in conventional complex fluids., Comment: 8 pages, 7 figures, minor changes, Phys. Rev. C in press

Nucleosynthesis, on the surface of accreting neutron stars, produces a range of chemical elements. We perform molecular dynamics simulations of crystallization to see how this complex composition forms new neutron star crust. We find chemical separation, with the liquid ocean phase greatly enriched in low atomic number elements compared to the solid crust. This phase separation should change many crust properties such as the thermal conductivity and shear modulus. The concentration of carbon, if present, is enriched in the ocean. This may allow unstable thermonuclear burning of the carbon and help explain the ignition of the very energetic explosions known as superbursts., Comment: 8 pages, 6 figures, minor changes, Physical Review E in press

The nuclear pasta -- a novel state of matter having nucleons arranged in a variety of complex shapes -- is expected to be found in the crust of neutron stars and in core-collapse supernovae at subnuclear densities of about $10^{14}$ g/cm$^3$. Due to frustration, a phenomenon that emerges from the competition between short-range nuclear attraction and long-range Coulomb repulsion, the nuclear pasta displays a preponderance of unique low-energy excitations. These excitations could have a strong impact on many transport properties, such as neutrino propagation through stellar environments. The excitation spectrum of the nuclear pasta is computed via a molecular-dynamics simulation involving up to 100,000 nucleons. The dynamic response of the pasta displays a classical plasma oscillation in the 1-2 MeV region. In addition, substantial strength is found at low energies. Yet this low-energy strength is missing from a simple ion model containing a single-representative heavy nucleus. The low-energy strength observed in the dynamic response of the pasta is likely to be a density wave involving the internal degrees of freedom of the clusters., Comment: 4 pages, 3 figures, Phys Rev C in press