Search

We systematically investigate the mergers of neutron star-white dwarf binaries from beginning to end, with focus on the properties of the inflows and outflows in accretion disks and their electromagnetic emissions. Using population synthesis models, we determine a subset of these binaries in which the white dwarf companion undergoes unstable mass transfer and complete tidal disruption, forming a large accretion disk around the neutron star. The material evolves according to an advection-dominated accretion disk model with nuclear burning, neutrino-emissions, and disk-surface wind ejection. The extreme dynamics of the entire process has proven difficult for analytic analysis, and thus currently the properties are poorly understood. The outflows from the mergers are iron- and nickel-rich, giving rise to optical and infrared emissions powered from the decay of the radioactive iron-type isotopes, calculated via the SuperNu light-curve code. We find these systems capable of powering bright, yet short-lived, optical transients with the potential to power gamma-ray bursts., Comment: LA-UR-22-29879

Neutron star mergers are cosmic catastrophes that produce some of the most energetic observed phenomena: short gamma-ray bursts, gravitational wave signals, and kilonovae. The latter are optical transients, powered by radioactive nuclides which are synthesized when the neutron-rich ejecta of a disrupted neutron star undergoes decompression. We model this decompression phase using data from simulations of post-merger accretion disk winds. We use smoothed particle hydrodynamics with realistic nuclear heating to model the expansion over multiple scales, from initially several thousand km to billions of km. We then render a realistic image of a kilonova ejecta as it would appear for a nearby observer. This is the first time such a visualization is performed using input from state-of-the-art accretion disk simulations, nuclear physics and atomic physics. The volume rendering in our model computes an opacity transfer function on the basis of the physical opacity, varying significantly with the inhomogeneity of the neutron richness in the ejecta. Other physical quantities such as temperature or electron fraction can be visualized using an independent color transfer function. We discuss several difficulties with the ParaView application that we encountered during the visualization process, and give descriptions of our solutions and workarounds which could be used for future improvements., Comment: Presented at the Scientific Visualization Supercomputing Conference 2021, St. Louis, MO

We investigate the nuclear pasta phases in neutron star crusts by conducting a large number of three-dimensional Hartree-Fock+BCS calculations at densities leading to the crust-core transition. We survey the shape parameter space of pasta at constant pressure. Spaghetti, waffles, lasagna, bi-continuous phases and cylindrical holes occupy local minima in the resulting Gibbs energy surfaces. The bi-continuous phase, in which both the neutron gas and nuclear matter extend continuously in all dimensions and therefore protons are delocalized, appears over a large range of depths. Our results support the idea that nuclear pasta is a glassy system. Multiple pasta configurations coexist in a given layer of the crust. At a characteristic temperature, of order $10^8$-$10^9$K, different phases become frozen into domains whose sizes we estimate to be 1-50 times the lattice spacing and over which the local density and electron fraction can vary. Above this temperature, there is very little long-range order and matter is an amorphous solid. Electron scattering off domain boundaries may contribute to the disorder resistivity of the pasta phases. Annealing of the domains may occur during cooling; repopulating of local minima during crustal heating might lead to temperature dependent transport properties in the deep layers of the crust. We identify 4 distinct regions: (1) nuclear pasta first appears as a local minima, but spherical nuclei are the ground state; (2) nuclear pasta become the absolute minimum, but spherical nuclei are still a local minimum (3) only nuclear pasta appears in local minima, and protons are still localized in at least one dimension (4) only pasta appears, and protons are delocalized. The whole pasta region can occupy up to 70% of the crust by mass and 40% by thickness, and the layer in which protons are delocalized could occupy 45% of the crust mass and 25% of its thickness., Comment: 32 pages, 23 figures, submitted to Physical Review C

FleCSPH is a smoothed particle hydrodynamics simulation tool, based on the compile-time configurable framework FleCSI. The asynchronous distributed tree topology combined with a fast multipole method allows FleCSPH to efficiently compute hydrodynamics and long range particle-particle interactions. FleCSPH provides initial data generators, particle relaxation techniques, and standard evolution drivers, which can be easily modified and extended to user-specific setups. Data input/output uses the H5part format, compatible with modern visualization software.

A two-phase description of the quark-nuclear matter hybrid equation of state that takes into account the effect of excluded volume in both the hadronic and the quark matter phases is introduced. The nuclear phase manifests a reduction of the available volume as density increases, leading to a stiffening of the matter. The quark matter phase displays a reduction of the effective string tension in the confining density functional from available volume contributions. The nuclear equation of state is based upon the relativistic density functional model DD2 with excluded volume. The quark matter is based upon a mean-field modification to the free fermi gas and will be discussed in greater detail. The interactions are decomposed into mean scalar and vector components. The scalar interaction is motivated by a string potential between quarks, whereas the vector interaction potential is motivated by higher order interactions of quarks leading to an increased stiffening at high densities. As an application, we consider matter under compact star constraints of electric neutrality and $\beta$-equilibrium. We obtain mass-radius relations for hybrid stars that form a third family, disconnected from the purely hadronic star branch, and fulfill the $2M_\odot$ constraint., Comment: 12 pages, 18 figures, two figures added and others remastered, text and references added, typos removed