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.

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.

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.

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!

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.

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.

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.

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.

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.

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.

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.

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.

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.

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}$

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.

Muons can be created in nascent neutron stars (NSs) due to the high electron chemical potentials and the high temperatures. Because of their relatively lower abundance compared to electrons, their role has so far been ignored in numerical simulations of stellar core collapse and NS formation. However, the appearance of muons softens the NS equation of state, triggers faster NS contraction and thus leads to higher luminosities and mean energies of the emitted neutrinos. This strengthens the postshock heating by neutrinos and can facilitate explosions by the neutrino-driven mechanism.

Core collapse supernova simulations can be sensitive to neutrino interactions near the neutrinosphere. This is the surface of last scattering. We model the neutrinosphere region as a warm unitary gas of neutrons. A unitary gas is a low density system of particles with large scattering lengths. We calculate modifications to neutrino scattering cross sections because of the universal spin and density correlations of a unitary gas. These correlations can be studied in laboratory cold atom experiments. We find significant reductions in cross sections, compared to free space interactions, even at relatively low densities. These reductions could reduce the delay time from core bounce to successful explosion in multidimensional supernova simulations.

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.

We define `astromaterial science' as the study of materials in astronomical objects that are qualitatively denser than materials on earth. Astromaterials can have unique properties related to their large density, though they may be organized in ways similar to more conventional materials. By analogy to terrestrial materials, we divide our study of astromaterials into hard and soft and discuss one example of each. The hard astromaterial discussed here is a crystalline lattice, such as the Coulomb crystals in the interior of cold white dwarfs and in the crust of neutron stars, while the soft astromaterial is nuclear pasta found in the inner crusts of neutron stars. In particular, we discuss how molecular dynamics simulations have been used to calculate the properties of astromaterials to interpret observations of white dwarfs and neutron stars. Coulomb crystals are studied to understand how compact stars freeze. Their incredible strength may make crust "mountains" on rotating neutron stars a source for gravitational waves that the Laser Interferometer Gravitational-Wave Observatory (LIGO) may detect. Nuclear pasta is expected near the base of the neutron star crust at densities of $10^{14}$ g/cm$^3$. Competition between nuclear attraction and Coulomb repulsion rearranges neutrons and protons into complex non-spherical shapes such as sheets (lasagna) or tubes (spaghetti). Semi-classical molecular dynamics simulations of nuclear pasta have been used to study these phases and calculate their transport properties such as neutrino opacity, thermal conductivity, and electrical conductivity. Observations of neutron stars may be sensitive to these properties, and can be be used to interpret observations of supernova neutrinos, magnetic field decay, and crust cooling of accreting neutron stars. We end by comparing nuclear pasta shapes with some similar shapes seen in biological systems.

Complex and exotic nuclear geometries are expected to appear naturally in dense nuclear matter found in the crust of neutron stars and supernovae environment collectively referred to as nuclear pasta. The pasta geometries depend on the average baryon density, proton fraction and temperature and are critically important in the determination of many transport properties of matter in supernovae and the crust of neutron stars. Using a set of self-consistent microscopic nuclear energy density functionals we present the first results of large scale quantum simulations of pasta phases at baryon densities $0.03 \leq \rho \leq 0.10$ fm$^{-3}$, proton fractions $0.05 \leq Y_p \leq 0.40$, and zero temperature. The full quantum simulations, in particular, allow us to thoroughly investigate the role and impact of the nuclear symmetry energy on pasta configurations. We use the Sky3D code that solves the Skyrme Hartree-Fock equations on a three-dimensional Cartesian grid. For the nuclear interaction we use the state of the art UNEDF1 parametrization, which was introduced to study largely deformed nuclei, hence is suitable for studies of the nuclear pasta. Density dependence of the nuclear symmetry energy is simulated by tuning two purely isovector observables that are insensitive to the current available experimental data. We find that a minimum total number of nucleons $A=2000$ is necessary to prevent the results from containing spurious shell effects and to minimize finite size effects. We find that a variety of nuclear pasta geometries are present in the neutron star crust and the result strongly depends on the nuclear symmetry energy. The impact of the nuclear symmetry energy is less pronounced as the proton fractions increase. Quantum nuclear pasta calculations at $T=0$ MeV are shown to get easily trapped in meta-stable states, and possible remedies to avoid meta-stable solutions are discussed.

Nuclear pasta, with nucleons arranged into tubes, sheets, or other complex shapes, is expected in core collapse supernovae (SNe) at just below nuclear density. We calculate the additional opacity from neutrino-pasta coherent scattering using molecular dynamics simulations. We approximately include this opacity in simulations of SNe. We find that pasta slows neutrino diffusion and greatly increases the neutrino signal at late times of 10 or more seconds after stellar core collapse. This signal, for a galactic SN, should be clearly visible in large detectors such as Super-Kamiokande.

We generalize our virial approach to study the neutral current neutrino response of nuclear matter at low densities. In the long-wavelength limit, the virial expansion makes model-independent predictions for neutrino-nucleon scattering rates and the density S_V and spin S_A responses. We find S_A is significantly reduced from one even at low densities. We provide a simple fit S_A^f(n,T,Y_p) of the axial response as a function of density n, temperature T and proton fraction Y_p. This fit reproduces our model independent virial results at low densities and reproduces the Burrows and Sawyer random phase approximation (RPA) results at high densities. Our fit can be incorporated into supernova simulations in a straight forward manner. Preliminary one dimensional supernova simulations suggest that the reduction in the axial response may enhance neutrino heating rates in the gain region during the accretion phase of a core-collapse supernovae.

We show that observations of the core temperature of transiently-accreting neutron stars combined with observations of an accretion outburst give a lower limit to the neutron star core heat capacity. For the neutron stars in the low mass X-ray binaries KS 1731-260, MXB 1659-29, and XTE J1701-462, we show that the lower limit is a factor of a few below the core heat capacity expected if neutrons and protons in the core are paired, so that electrons provide the dominant contribution to the heat capacity. This limit rules out a core dominated by a quark color-flavor-locked (CFL) phase, which would have a much lower heat capacity. Future observations of or limits on cooling during quiescence will further constrain the core heat capacity.

Neutron star and supernova matter at densities just below the nuclear matter saturation density is expected to form a lattice of exotic shapes. These so-called nuclear pasta phases are caused by Coulomb frustration. Their elastic and transport properties are believed to play an important role for thermal and magnetic field evolution, rotation and oscillation of neutron stars. Furthermore, they can impact neutrino opacities in core-collapse supernovae. In this work, we present proof-of-principle 3D Skyrme Hartree-Fock (SHF) simulations of nuclear pasta with the Multi-resolution ADaptive Numerical Environment for Scientific Simulations (MADNESS). We perform benchmark studies of $^{16} \mathrm{O}$, $^{208} \mathrm{Pb}$ and $^{238} \mathrm{U}$ nuclear ground states and calculate binding energies via 3D SHF simulations. Results are compared with experimentally measured binding energies as well as with theoretically predicted values from an established SHF code. The nuclear pasta simulation is initialized in the so-called waffle geometry as obtained by the Indiana University Molecular Dynamics (IUMD) code. The size of the unit cell is 24\:fm with an average density of about $\rho = 0.05 \:\mathrm{fm}^{-3}$, proton fraction of $Y_p = 0.3$ and temperature of $T=0\:$MeV. Our calculations reproduce the binding energies and shapes of light and heavy nuclei with different geometries. For the pasta simulation, we find that the final geometry is very similar to the initial waffle state. In the present pasta calculations spin-orbit forces are not included but will be added in the future. Within the MADNESS framework, we can successfully perform calculations of inhomogeneous nuclear matter. By using pasta configurations from IUMD it is possible to explore different geometries and test the impact of self-consistent calculations on the latter.

[Background] The "pasta" phase of nuclear matter may play an important role in the structure and evolution of neutron stars. Recent works suggest nuclear pasta has a high resistivity which could be explained by the presence of long lived topological defects. The defects act as impurities that decrease thermal and electrical conductivity of the pasta. [Purpose] To quantify how topological defects affect transport properties of nuclear pasta and estimate this effect using an impurity parameter $Q_{\text{imp}}$. [Methods] Contrast molecular dynamics simulations of up to 409\,600 nucleons arranged in parallel nuclear pasta slabs (perfect pasta) with simulations of pasta slabs connected by topological defects (impure pasta). From these simulations compare the viscosity and heat conductivity of perfect and impure pasta to obtain an effective impurity parameter $Q_{\text{imp}}$ due to the presence of defects. [Results] Both the viscosity and thermal conductivity calculated for both perfect and impure pasta are anisotropic, peaking along directions perpendicular to the slabs and reaching a minimum close to zero parallel to them. In our 409\,600 nucleon simulation topological defects connecting slabs of pasta reduce both the thermal conductivity and viscosity on average by about 37\%. We estimate an effective impurity parameter due to the defects of order $Q_{\text{imp}}\sim30$. [Conclusions] Topological defects in the pasta phase of nuclear matter have an effect similar to impurities in a crystal lattice. The irregularities introduced by the defects reduce the thermal and electrical conductivities and the viscosity of the system. This effect can be parameterized by a large impurity parameter $Q_{\text{imp}}\sim30$.

A striking shape was recently observed for the cellular organelle endoplasmic reticulum consisting of stacked sheets connected by helical ramps. 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.

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.

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).

Pulsar glitches-the sudden spin-up in the rotational frequency of a neutron star-suggest the existence of an angular-momentum reservoir confined to the inner crust of the neutron star. Large and regular glitches observed in the Vela pulsar have originally constrained the fraction of the stellar moment of inertia that must reside in the solid crust to about 1.4%. However, crustal entrainment-which until very recently has been ignored-suggests that in order to account for the Vela glitches, the fraction of the moment of inertia residing in the crust must increase to about 7%. This indicates that the required angular momentum reservoir may exceed that which is available in the crust. We explore the possibility that uncertainties in the equation of state provide enough flexibility for the construction of models that predict a large crustal thickness and consequently a large crustal moment of inertia. Given that analytic results suggest that the crustal moment of inertia is sensitive to the transition pressure at the crust-core interface, we tune the parameters of the model to maximize the transition pressure, while still providing an excellent description of nuclear observables. In this manner we are able to obtain fractional moments of inertia as large as 7% for neutron stars with masses below 1.6 solar masses. In particular, we find that if the neutron-skin thickness of 208Pb falls within the (0.20-0.26) fm range, large enough transition pressures can be generated to explain the large Vela glitches without invoking an additional angular-momentum reservoir beyond that confined to the solid crust. Our results suggest that the crust may be enough.

The symmetry energy describes how the energy of nuclear matter rises as one goes away from equal numbers of neutrons and protons. This is very important to describe neutron rich matter in astrophysics. This article reviews our knowledge of the symmetry energy from theoretical calculations, nuclear structure measurements, heavy ion collisions, and astronomical observations. We then present a roadmap to make progress in areas of relevance to the symmetry energy that promotes collaboration between the astrophysics and the nuclear physics communities.

The formation of complex nonuniform phases of nuclear matter, known as nuclear pasta, is studied with molecular dynamics simulations containing 51200 nucleons. A phenomenological nuclear interaction is used that reproduces the saturation binding energy and density of nuclear matter. Systems are prepared at an initial density of 0.10fm$^{-3}$ and then the density is decreased by expanding the simulation volume at different rates to densities of 0.01 fm$^{-3}$ or less. An originally uniform system of nuclear matter is observed to form spherical bubbles ("swiss cheese"), hollow tubes, flat plates ("lasagna"), thin rods ("spaghetti") and, finally, nearly spherical nuclei with decreasing density. We explicitly observe nucleation mechanisms, with decreasing density, for these different pasta phase transitions. Topological quantities known as Minkowski functionals are obtained to characterize the pasta shapes. Different pasta shapes are observed depending on the expansion rate. This indicates non equilibrium effects. We use this to determine the best ways to obtain lower energy states of the pasta system from MD simulations and to place constrains on the equilibration time of the system.

At very high densities, electrons react with protons to form neutron rich matter. This material is central to many fundamental questions in nuclear physics and astrophysics. Moreover, neutron rich matter is being studied with an extraordinary variety of new tools such as the Facility for Rare Isotope Beams (FRIB) and the Laser Interferometer Gravitational Wave Observatory (LIGO). We describe the Lead Radius Experiment (PREX) that uses parity violating electron scattering to measure the neutron radius of 208Pb. This has important implications for neutron stars and their crusts. We discuss X-ray observations of neutron star radii. These also have important implications for neutron rich matter. Gravitational waves (GW) open a new window on neutron rich matter. They come from sources such as neutron star mergers, rotating neutron star mountains, and collective r-mode oscillations. Using large scale molecular dynamics simulations, we find neutron star crust to be very strong. It can support mountains on rotating neutron stars large enough to generate detectable gravitational waves. Finally, neutrinos from core collapse supernovae (SN) provide another, qualitatively different probe of neutron rich matter. Neutrinos escape from the surface of last scattering known as the neutrino-sphere. This is a low density warm gas of neutron rich matter. Neutrino-sphere conditions can be simulated in the laboratory with heavy ion collisions. Observations of neutrinos can probe nucleosyntheses in SN. We believe that combing astronomical observations using photons, GW, and neutrinos, with laboratory experiments on nuclei, heavy ion collisions, and radioactive beams will fundamentally advance our knowledge of compact objects in the heavens, the dense phases of QCD, the origin of the elements, and of neutron rich matter.

Core-collapse supernovae may depend sensitively on charged current neutrino interactions in warm, low density neutron rich matter. A proton in neutron rich matter is more tightly bound than is a neutron. This energy shift \Delta U increases the electron energy in \nu_e + n --> p + e, increasing the available phase space and absorption cross section. Likewise \Delta U decreases the positron energy in \bar \nu_e + p --> n + e^+, decreasing the phase space and cross section. We have calculated \Delta U using a model independent virial expansion and we find \Delta U is much larger, at low densities, than the predictions of many mean field models. Therefore \Delta U could have a significant impact on charged current neutrino interactions in supernovae. Preliminary simulations of the accretion phase of core-collapse supernovae find that \Delta U increases \bar \nu_e energies and decreases the \nu_e luminosity.

The neutron rich isotope 22Ne may be a significant impurity in carbon and oxygen white dwarfs and could impact how the stars freeze. We perform molecular dynamics simulations to determine the influence of 22Ne in carbon-oxygen-neon systems on liquid-solid phase equilibria. Both liquid and solid phases are present simultaneously in our simulation volumes. We identify liquid, solid, and interface regions in our simulations using a bond angle metric. In general we find good agreement for the composition of liquid and solid phases between our MD simulations and the semi analytic model of Medin and Cumming. The trace presence of a third component, neon, does not appear to strongly impact the chemical separation found previously for two component carbon and oxygen systems. This suggests that small amounts of 22Ne may not qualitatively change how the material in white dwarf stars freezes. However, we do find systematically lower melting temperatures (higher Gamma) in our MD simulations compared to the semi analytic model. This difference seems to grow with impurity parameter Q_imp and suggests a problem with simple corrections to the linear mixing rule for the free energy of multicomponent solid mixtures that is used in the semi analytic model.

We determine the liquid-solid phase diagram for carbon-oxygen and oxygen-selenium plasma mixtures using two-phase MD simulations. We identified liquid, solid, and interface regions using a bond angle metric. To study finite size effects, we perform 27648 and 55296 ion simulations. To help monitor non-equilibrium effects, we calculate diffusion constants $D_i$. For the carbon-oxygen system we find that $D_O$ for oxygen ions in the solid is much smaller than $D_C$ for carbon ions and that both diffusion constants are 80 or more times smaller than diffusion constants in the liquid phase. There is excellent agreement between our carbon-oxygen phase diagram and that predicted by Medin and Cumming. This suggests that errors from finite size and non-equilibrium effects are small and that the carbon-oxygen phase diagram is now accurately known. The oxygen-selenium system is a simple two-component model for more complex rapid proton capture nucleosynthesis ash compositions for an accreting neutron star. Diffusion of oxygen, in a predominately selenium crystal, is remarkably fast, comparable to diffusion in the liquid phase. We find a somewhat lower melting temperature for the oxygen-selenium system than that predicted by Medin and Cumming. This is probably because of electron screening effects.

Neutron stars are attractive places to look for dark matter because their high densities allow repeated interactions. Weakly interacting massive particles (WIMPs) may scatter efficiently in the core or in the crust of a neutron star. In this paper we focus on WIMP contributions to transport properties, such as shear viscosity or thermal conductivity, because these can be greatly enhanced by long mean free paths. We speculate that WIMPs increase the shear viscosity of neutron star matter and help stabilize r-mode oscillations. These are collective oscillations where the restoring force is the Coriolis force. At present r-modes are thought to be unstable in many observed rapidly rotating stars. If WIMPs stabilize the r-modes, this would allow neutron stars to spin rapidly. This likely requires WIMP-nucleon cross sections near present experimental limits and an appropriate density of WIMPs in neutron stars.

Fusion of very neutron rich nuclei may be important to determine the composition and heating of the crust of accreting neutron stars. Fusion cross sections are calculated using time-dependent Hartree-Fock theory coupled with density-constrained Hartree-Fock calculations to deduce an effective potential. Systems studied include 16O+16O, 16O+24O, 24O+24O, 12C+16O, and 12C+24O. We find remarkable agreement with experimental cross sections for the fusion of stable nuclei. Our simulations use the SLy4 Skyrme force that has been previously fit to the properties of stable nuclei, and no parameters have been fit to fusion data. We compare our results to the simple S\~{a}o Paulo static barrier penetration model. For the asymmetric systems 12C+24O or 16O+24O we predict an order of magnitude larger cross section than those predicted by the S\~{a}o Paulo model. This is likely due to the transfer of neutrons from the very neutron rich nucleus to the stable nucleus and dynamical rearrangements of the nuclear densities during the collision process. These effects are not included in potential models. This enhancement of fusion cross sections, for very neutron rich nuclei, can be tested in the laboratory with radioactive beams.

The symmetry energy contribution to the nuclear Equation of State (EoS) impacts various phenomena in nuclear astrophysics, nuclear structure, and nuclear reactions. Its determination is a key objective of contemporary nuclear physics with consequences for the understanding of dense matter within neutron stars. We examine the results of laboratory experiments that have provided initial constraints on the nuclear symmetry energy and its density dependence at and somewhat below normal nuclear matter density. Some of these constraints have been derived from properties of nuclei. Others have been derived from the nuclear response to electroweak and hadronic probes. We also examine the most frequently used theoretical models that predict the symmetry energy and its slope. By comparing existing constraints on the symmetry pressure to theories, we demonstrate how the contribution of the three-body force, an essential ingredient in neutron matter models, can be determined.

At very high densities, electrons react with protons to form neutron rich matter. This material is central to many fundamental questions in nuclear physics and astrophysics. Moreover, neutron rich matter is being studied with an extraordinary variety of new tools such as the Facility for Rare Isotope Beams (FRIB) and the Laser Interferometer Gravitational Wave Observatory (LIGO). We describe the Lead Radius Experiment (PREX) that uses parity violating electron scattering to measure the neutron radius in $^{208}$Pb. This has important implications for neutron stars and their crusts. We discuss X-ray observations of neutron star radii. These also have important implications for neutron rich matter. Gravitational waves (GW) open a new window on neutron rich matter. They come from sources such as neutron star mergers, rotating neutron star mountains, and collective r-mode oscillations. Using large scale molecular dynamics simulations, we find neutron star crust to be very strong. It can support mountains on rotating neutron stars large enough to generate detectable gravitational waves. We believe that combing astronomical observations using photons, GW, and neutrinos, with laboratory experiments on nuclei, heavy ion collisions, and radioactive beams will fundamentally advance our knowledge of compact objects in the heavens, the dense phases of QCD, the origin of the elements, and of neutron rich matter.

Neutron rich matter is central to many fundamental questions in nuclear physics and astrophysics. Moreover, this material is being studied with an extraordinary variety of new tools such as the Facility for Rare Isotope Beams (FRIB) and the Laser Interferometer Gravitational Wave Observatory (LIGO). We describe the Lead Radius Experiment (PREX) that uses parity violating electron scattering to measure the neutron radius in $^{208}$Pb. This has important implications for neutron stars and their crusts. We discuss X-ray observations of neutron star radii. These also have important implications for neutron rich matter. Gravitational waves (GW) open a new window on neutron rich matter. They come from sources such as neutron star mergers, rotating neutron star mountains, and collective r-mode oscillations. Using large scale molecular dynamics simulations, we find neutron star crust to be very strong. It can support mountains on rotating neutron stars large enough to generate detectable gravitational waves. Finally, neutrinos from core collapse supernovae (SN) provide another, qualitatively different probe of neutron rich matter. Neutrinos escape from the surface of last scattering known as the neutrino-sphere. This is a low density warm gas of neutron rich matter. Observations of neutrinos can probe nucleosyntheses in SN. Simulations of SN depend on the equation of state (EOS) of neutron rich matter. We discuss a new EOS based on virial and relativistic mean field calculations. We believe that combing astronomical observations using photos, GW, and neutrinos, with laboratory experiments on nuclei, heavy ion collisions, and radioactive beams will fundamentally advance our knowledge of compact objects in the heavens, the dense phases of QCD, the origin of the elements, and of neutron rich matter.

In this book chapter we review plasma crystals in the laboratory, in the interior of white dwarf stars, and in the crust of neutron stars. We describe a molecular dynamics formalism and show results for many neutron star crust properties including phase separation upon freezing, diffusion, breaking strain, shear viscosity and dynamics response of nuclear pasta. We end with a summary and discuss open questions and challenges for the future.

Diffusion in coulomb crystals can be important for the structure of neutron star crusts. We determine diffusion constants $D$ from molecular dynamics simulations. We find that $D$ for coulomb crystals with relatively soft-core $1/r$ interactions may be larger than $D$ for Lennard-Jones or other solids with harder-core interactions. Diffusion, for simulations of nearly perfect body-centered-cubic lattices, involves the exchange of ions in ring-like configurations. Here ions "hop" in unison without the formation of long lived vacancies. Diffusion, for imperfect crystals, involves the motion of defects. Finally, we find that diffusion, for an amorphous system rapidly quenched from coulomb parameter $\Gamma=175$ to coulomb parameters up to $\Gamma=1750$, is fast enough so that the system starts to crystallize during long simulation runs. These results strongly suggest that coulomb solids in cold white dwarf stars, and the crust of neutron stars, will be crystalline and not amorphous.

We generate a second equation of state (EOS) of nuclear matter for a wide range of temperatures, densities, and proton fractions for use in supernovae, neutron star mergers, and black hole formation simulations. We employ full relativistic mean field (RMF) calculations for matter at intermediate density and high density, and the Virial expansion of a non-ideal gas for matter at low density. For this EOS we use the RMF effective interaction FSUGold, whereas our earlier EOS was based on the RMF effective interaction NL3. The FSUGold interaction has a lower pressure at high densities compared to the NL3 interaction. We calculate the resulting EOS at over 100,000 grid points in the temperature range $T$ = 0 to 80 MeV, the density range $n_B$ = 10$^{-8}$ to 1.6 fm$^{-3}$, and the proton fraction range $Y_p$ = 0 to 0.56. We then interpolate these data points using a suitable scheme to generate a thermodynamically consistent equation of state table on a finer grid. We discuss differences between this EOS, our NL3 based EOS, and previous EOSs by Lattimer-Swesty and H. Shen et al for the thermodynamic properties, composition, and neutron star structure. The original FSUGold interaction produces an EOS, that we call FSU1.7, that has a maximum neutron star mass of 1.7 solar masses. A modification in the high density EOS is introduced to increase the maximum neutron star mass to 2.1 solar masses and results in a slightly different EOS that we call FSU2.1. The EOS tables for FSU1.7 and FSU2.1 are available for download.

We generate a new complete equation of state (EOS) of nuclear matter for a wide range of temperatures, densities, and proton fractions ready for use in astrophysical simulations of supernovae and neutron star mergers. Our previous two papers tabulated the EOS at over 180,000 grid points in the temperature range $T$ = 0 to 80 MeV, the density range $n_B$ = 10$^{-8}$ to 1.6 fm$^{-3}$, and the proton fraction range $Y_P$ = 0 to 0.56. In this paper we combine these data points using a suitable interpolation scheme to generate a single equation of state table on a finer grid. This table is thermodynamically consistent and conserves entropy during adiabatic compression tests. We present various thermodynamic quantities and the composition of matter in the new EOS, along with several comparisons with existing EOS tables. Our equation of state table is available for download.

Sedimentation of the neutron rich isotope $^{22}$Ne may be an important source of gravitational energy during the cooling of white dwarf stars. This depends on the diffusion constant for $^{22}$Ne in strongly coupled plasma mixtures. We calculate self-diffusion constants $D_i$ from molecular dynamics simulations of carbon, oxygen, and neon mixtures. We find that $D_i$ in a mixture does not differ greatly from earlier one component plasma results. For strong coupling (coulomb parameter $\Gamma>$ few), $D_i$ has a modest dependence on the charge $Z_i$ of the ion species, $D_i \propto Z_i^{-2/3}$. However $D_i$ depends more strongly on $Z_i$ for weak coupling (smaller $\Gamma$). We conclude that the self-diffusion constant $D_{\rm Ne}$ for $^{22}$Ne in carbon, oxygen, and neon plasma mixtures is accurately known so that uncertainties in $D_{\rm Ne}$ should be unimportant for simulations of white dwarf cooling.

We study the equation of state (EOS) of nuclear matter at subnuclear density in a Virial expansion for a nonideal gas. The gas consists of neutrons, protons, alpha particles, and 8980 species of nuclei with $A \ge 12$ and masses from the finite-range droplet model (FRDM). At very low density, the Virial expansion reduces to nuclear statistical equilibrium. At higher density, the Virial results match smoothly to the relativistic mean field results discussed in our previous paper. We tabulate the resulting EOS at over 73,000 grid points in the temperature range $T$ = 0.158 to 15.8 MeV, the density range $n_B$ = 10$^{-8}$ to 0.1 fm$^{-3}$, and the proton fraction range $Y_P$ = 0.05 to 0.56. In the future we plan to match these low density results to our earlier high density mean field results, and generate a full EOS table for use in supernova and neutron star merger simulations. This Virial EOS is exact in the low density limit.

Nuclear effective interactions are useful tools in astrophysical applications especially if one can guide the extrapolations to the extremes regions of isospin and density that are required to simulate dense, neutron-rich systems. Isospin extrapolations may be constrained in the laboratory by measuring the neutron skin thickness of a heavy nucleus, such as 208Pb. Similarly, future observations of massive neutron stars will constrain the extrapolations to the high-density domain. In this contribution we introduce a new relativistic effective interaction that is simultaneously constrained by the properties of finite nuclei, their collective excitations, and neutron-star properties. By adjusting two of the empirical parameters of the theory, one can efficiently tune the neutron skin thickness of 208Pb and the maximum neutron star mass. We illustrate this procedure in response to the recent interpretation of X-ray observations by Steiner, Lattimer, and Brown that suggests that the FSUGold effective interaction predicts neutron star radii that are too large and a maximum stellar mass that is too small. The new effective interaction is fitted to a neutron skin thickness in 208Pb of only R_n - R_p = 0.16 fm and a moderately large maximum neutron star mass of 1.94 Msun.

Neutron rich matter is at the heart of many fundamental questions in Nuclear Physics and Astrophysics. What are the high density phases of QCD? Where did the chemical elements come from? What is the structure of many compact and energetic objects in the heavens, and what determines their electromagnetic, neutrino, and gravitational-wave radiations? Moreover, neutron rich matter is being studied with an extraordinary variety of new tools such as Facility for Rare Isotope Beams (FRIB) and the Laser Interferometer Gravitational Wave Observatory (LIGO). We describe the Lead Radius Experiment (PREX) that is using parity violation to measure the neutron radius in 208Pb. This has important implications for neutron stars and their crusts. Using large scale molecular dynamics, we model the formation of solids in both white dwarfs and neutron stars. We find neutron star crust to be the strongest material known, some 10 billion times stronger than steel. It can support mountains on rotating neutron stars large enough to generate detectable gravitational waves. Finally, we describe a new equation of state for supernova and neutron star merger simulations based on the Virial expansion at low densities, and large scale relativistic mean field calculations.

The breaking stress (the maximum of the stress-strain curve) of neutron star crust is important for neutron star physics including pulsar glitches, emission of gravitational waves from static mountains, and flares from star quakes. We perform many molecular dynamic simulations of the breaking stress at different coupling parameters (inverse temperatures) and strain rates. We describe our results with the Zhurkov model of strength. We apply this model to estimate the breaking stress for timescales ~1 s - 1 year, which are most important for applications, but much longer than can be directly simulated. At these timescales the breaking stress depends strongly on the temperature. For coupling parameter <200, matter breaks at very small stress, if it is applied for a few years. This viscoelastic creep can limit the lifetime of mountains on neutron stars. We also suggest an alternative model of timescale-independent breaking stress, which can be used to estimate an upper limit on the breaking stress.

We calculate the equation of state (EoS) of dense matter, using a relativistic mean field (RMF) model with a density dependent coupling that is a slightly modified form of the original NL3 interaction. For nonuniform nuclear matter we approximate the unit lattice as a spherical Wigner-Seitz cell, wherein the meson mean fields and nucleon Dirac wave functions are solved fully self-consistently. We also calculate uniform nuclear matter for a wide range of temperatures, densities, and proton fractions, and match them to non-uniform matter as the density decreases. The calculations took over 6,000 CPU days in Indiana University's supercomputer clusters. We tabulate the resulting EoS at over 107,000 grid points in the proton fraction range $Y_P$ = 0 to 0.56. For the temperature range $T$ = 0.16 to 15.8 MeV we cover the density range $n_B$ = 10$^{-4}$ to 1.6 fm$^{-3}$; and for the higher temperature range $T$ = 15.8 to 80 MeV we cover the larger density range $n_B$ = 10$^{-8}$ to 1.6 fm$^{-3}$. In the future we plan to study low density, low temperature (T$<$15.8 MeV), nuclear matter using a Virial expansion, and we will match the low density and high density results to generate a complete EoS table for use in astrophysical simulations of supernova and neutron star mergers.

We determine the phase diagram for dense carbon/ oxygen mixtures in White Dwarf (WD) star interiors using molecular dynamics simulations involving liquid and solid phases. Our phase diagram agrees well with predictions from Ogata et al. and Medin and Cumming and gives lower melting temperatures than Segretain et al. Observations of WD crystallization in the globular cluster NGC 6397 by Winget et al. suggest that the melting temperature of WD cores is close to that for pure carbon. If this is true, our phase diagram implies that the central oxygen abundance in these stars is less than about 60%. This constraint, along with assumptions about convection in stellar evolution models, limits the effective S factor for the $^{12}$C($\alpha,\gamma$)$^{16}$O reaction to S_{300} <= 170 keV barns.

Low mass neutron stars may be uniquely strong sources of gravitational waves (GW). The neutron star crust can support large deformations for low mass stars. This is because of the star's weaker gravity. We find maximum ellipticities $\epsilon$ (fractional difference in moments of inertia) that are 1000 times larger, and maximum quadrupole moments $Q_{22}$ over 100 times larger, for low mass stars than for 1.4 $M_\odot$ neutron stars. Indeed, we calculate that the crust can support an $\epsilon$ as large as 0.01 for a minimum mass neutron star. A 0.12 $M_\odot$ star, that is maximally strained and rotating at 100 Hz, will produce a characteristic gravitational wave strain of $h_0=2.1\times 10^{-24}$ at a distance of 1 kpc. The GW detector Advanced LIGO should be sensitive to such objects through out the Milky Way Galaxy.

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.

Mountains on rapidly rotating neutron stars efficiently radiate gravitational waves. The maximum possible size of these mountains depends on the breaking strain of neutron star crust. With multi-million ion molecular dynamics simulations of Coulomb solids representing the crust, we show that the breaking strain of pure single crystals is very large and that impurities, defects, and grain boundaries only modestly reduce the breaking strain to around 0.1. Due to the collective behavior of the ions during failure found in our simulations, the neutron star crust is likely very strong and can support mountains large enough so that their gravitational wave radiation could limit the spin periods of some stars and might be detectable in large scale interferometers. Furthermore, our microscopic modeling of neutron star crust material can help analyze mechanisms relevant in magnetar giant and micro flares.

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.

Torsional (shear) oscillations of neutron stars may have been observed in quasiperiodic oscillations of Magnetar Giant Flares. The frequencies of these modes depend on the shear modulus of neutron star crust. We calculate the shear modulus of Coulomb crystals from molecular dynamics simulations. We find that electron screening reduces the shear modulus by about 10% compared to previous Ogata et al. results. Our MD simulations can be extended to calculate the effects of impurities and or polycrystalline structures on the shear modulus.

The shear viscosity of the crust might have a damping effect on the amplitude of r-modes of rotating neutron stars. This damping has implications for the emission of gravitational waves. We calculate the contribution to the shear viscosity coming from the ions using both semi-analytical methods, that consider binary collisions, and Molecular Dynamics simulations. We compare these results with the contribution coming from electrons. We study how the shear viscosity depends on density for conditions of interest in neutron star envelopes and outer crusts. In the low density limit, we find good agreement between results of our molecular dynamics simulations and classical semi-analytic calculations.

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.

The equation of state (EOS) for nuclear and neutron rich matter is investigated in a Relativistic Mean Field (RMF) model. New shell states are found that minimize the free energy per baryon, calculated in a spherical Wigner-Seitz (WS) approximation, over a significant range of baryon densities. These shell states, that have both inside and outside surfaces, minimize the Coulomb energy of large proton number configurations at the expense of a larger surface energy. This is related to a possible depression in the central density of super heavy nuclei. As the baryon density increases, we find the system changes from normal nuclei, to shell states, and then to uniform matter.

We study the composition of the outer layers of a protoneutron star and show that light nuclei are present in substantial amounts. The composition is dominated by nucleons, deuterons, tritons and alpha particles; 3He is present in smaller amounts. This composition can be studied in laboratory experiments with new neutron-rich radioactive beams that can reproduce similar densities and temperatures. After including the corresponding neutrino interactions, we demonstrate that light nuclei have a small impact on the average energy of the emitted electron neutrinos, but are significant for the average energy of antineutrinos. During the early post-explosion phase, the average energy of electron antineutrinos is slightly increased, while at later times during the protoneutron star cooling it is reduced by about 1 MeV. The consequences of these changes for nucleosynthesis in neutrino-driven supernova outflows are discussed.

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.

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.

We extend the virial equation of state to include 3H and 3He nuclei, and predict significant mass-three fractions near the neutrinosphere in supernovae. While alpha particles are often more abundant, we demonstrate that energy transfer cross-sections for muon and tau neutrinos at low densities are dominated by breakup of the loosely-bound 3H and 3He nuclei. The virial coefficients involving A=3 nuclei are calculated directly from the corresponding nucleon-3H and nucleon-3He scattering phase shifts. For the neutral-current inelastic cross-sections and the energy transfer cross sections, we perform ab-initio calculations based on microscopic two- and three-nucleon interactions and meson-exchange currents.

We generalize our virial approach to study spin-polarized neutron matter and the consistent neutrino response at low densities. In the long-wavelength limit, the virial expansion makes model-independent predictions for the density and spin response, based only on nucleon-nucleon scattering data. Our results for the neutrino response provide constraints for random-phase approximation or other model calculations, and we compare the virial vector and axial response to response functions used in supernova simulations. The virial expansion is suitable to describe matter near the supernova neutrinosphere, and this work extends the virial equation of state to predict neutrino interactions in neutron matter.

Neutrinos in core collapse supernovae are likely trapped by neutrino-nucleus elastic scattering. Using molecular dynamics simulations, we calculate neutrino mean free paths and ion-ion correlation functions for heterogeneous plasmas. Mean free paths are systematically shorter in plasmas containing a mixture of ions compared to a plasma composed of a single ion species. This is because neutrinos can scatter from concentration fluctuations. The dynamical response function of a heterogeneous plasma is found to have an extra peak at low energies describing the diffusion of concentration fluctuations. Our exact molecular dynamics results for the static structure factor reduce to the Debye Huckel approximation, but only in the limit of very low momentum transfers.

We present the virial equation of state of low-density nuclear matter composed of neutrons, protons and alpha particles. The virial equation of state is model-independent, and therefore sets a benchmark for all nuclear equations of state at low densities. We calculate the second virial coefficients for nucleon-nucleon, nucleon-alpha and alpha-alpha interactions directly from the relevant binding energies and scattering phase shifts. The virial approach systematically takes into account contributions from bound nuclei and the resonant continuum, and consequently provides a framework to include strong-interaction corrections to nuclear statistical equilibrium models. The virial coefficients are used to make model-independent predictions for a variety of properties of nuclear matter over a range of densities, temperatures and compositions. Our results provide constraints on the physics of the neutrinosphere in supernovae. The resulting alpha particle concentration differs from all equations of state currently used in supernova simulations. Finally, the virial equation of state greatly improves our conceptual understanding of low-density nuclear matter.

We present a model-independent description of low-density neutron matter based on the virial expansion. The virial equation of state provides a benchmark for all nuclear equations of state at densities and temperatures where the interparticle separation is large compared to the thermal wavelength. We calculate the second virial coefficient directly from the nucleon-nucleon scattering phase shifts. Our results for the pressure, energy, entropy and the free energy correctly include the physics of the large neutron-neutron scattering length. We find that, as in the universal regime, thermodynamic properties of neutron matter scale over a wide range of temperatures, but with a significantly reduced interaction coefficient compared to the unitary limit.

Heavy ion experiments provide important data to test astrophysical models. The high density equation of state can be probed in HI collisions and applied to the hot protoneutron star formed in core collapse supernovae. The Parity Radius Experiment (PREX) aims to accurately measure the neutron radius of $^{208}$Pb with parity violating electron scattering. This determines the pressure of neutron rich matter and the density dependence of the symmetry energy. Competition between nuclear attraction and coulomb repulsion can form exotic shapes called nuclear pasta in neutron star crusts and supernovae. This competition can be probed with multifragmentation HI reactions. We use large scale semiclassical simulations to study nonuniform neutron rich matter in supernovae. We find that the coulomb interactions in astrophysical systems suppress density fluctuations. As a result, there is no first order liquid vapor phase transition. Finally, the virial expansion for low density matter shows that the nuclear vapor phase is complex with significant concentrations of alpha particles and other light nuclei in addition to free nucleons.

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.

Core collapse supernovae are gigantic explosions of massive stars that radiate 99% of their energy in neutrinos. This provides a unique opportunity for large scale parity or charge conjugation violation. Parity violation in a strong magnetic field could lead to an asymmetry in the neutrino radiation and recoil of the newly formed neutron star. Charge conjugation violation in the neutrino-nucleon interaction reduces the ratio of neutrons to protons in the neutrino driven wind above the neutron star. This is a problem for r-process nucleosynthesis in this wind. On earth, parity violation is an excellent probe of neutrons because the weak charge of a neutron is much larger than that of a proton. The Parity Radius Experiment (PREX) at Jefferson Laboratory aims to precisely measure the neutron radius of $^{208}$Pb with parity violating elastic electron scattering. This has many implications for astrophysics, including the structure of neutron stars, and for atomic parity nonconservation experiments.}

Nonuniform neutron-rich matter present in both core-collapse supernovae and neutron-star crusts is described in terms of a semiclassical model that reproduces nuclear-matter properties and includes long-range Coulomb interactions. The neutron-neutron correlation function and the corresponding static structure factor are calculated from molecular dynamics simulations involving 40,000 to 100,000 nucleons. The static structure factor describes coherent neutrino scattering which is expected to dominate the neutrino opacity. At low momentum transfers the static structure factor is found to be small because of ion screening. In contrast, at intermediate momentum transfers the static structure factor displays a large peak due to coherent scattering from all the neutrons in a cluster. This peak moves to higher momentum transfers and decreases in amplitude as the density increases. A large static structure factor at zero momentum transfer, indicative of large density fluctuations during a first-order phase transition, may increase the neutrino opacity. However, no evidence of such an increase has been found. Therefore, it is unlikely that the system undergoes a simple first-order phase transition. It is found that corrections to the commonly used single heavy nucleus approximation first appear at a density of the order of $10^{13}$ g/cm$^3$ and increase rapidly with increasing density. Thus, neutrino opacities are overestimated in the single heavy nucleus approximation relative to the complete molecular dynamics simulations.

Neutron-rich matter at subnuclear densities may involve complex structures displaying a variety of shapes, such as spherical, slablike, and/or rodlike shapes. These phases of the {\it nuclear pasta} are expected to exist in the crust of neutron stars and in core-collapse supernovae. The dynamics of core-collapse supernovae is very sensitive to the interactions between neutrinos and nucleons/nuclei. Indeed, neutrino excitation of the low-energy modes of the pasta may allow for a significant energy transfer to the nuclear medium, thereby reviving the stalled supernovae shock. The linear response of the nuclear pasta to neutrinos is modeled via a simple semi-classical simulation. The transport mean-free path for $\mu$ and $\tau$ neutrinos (and antineutrinos) is expressed in terms of the static structure factor of the pasta, which is evaluated using Metropolis Monte Carlo simulations.

Advances in neutrino transport allow realistic neutrino interactions to be incorporated into supernova simulations. We add tensor couplings to relativistic RPA calculations of neutrino opacities. Our results reproduce free-space neutrino-nucleon cross sections at low density, including weak magnetism and recoil corrections. In addition, our opacities are thermodynamically consistent with relativistic mean field equations of state. We find antineutrino mean free paths that are considerably larger then those for neutrinos. This difference depends little on density. In a supernova, this difference could lead to an average energy of $\bar\nu_\mu$ that is larger than that for $\nu_\mu$ by an amount that is comparable to the energy difference between $\nu_\mu$ and $\bar\nu_e$

Development of large mass detectors for low-energy neutrinos and dark matter may allow supernova detection via neutrino-nucleus elastic scattering. An elastic-scattering detector could observe a few, or more, events per ton for a galactic supernova at 10 kpc ($3.1 \times 10^{20}$ m). This large yield, a factor of at least 20 greater than that for existing light-water detectors, arises because of the very large coherent cross section and the sensitivity to all flavors of neutrinos and antineutrinos. An elastic scattering detector can provide important information on the flux and spectrum of $\nu_\mu$ and $\nu_\tau$ from supernovae. We consider many detectors and a range of target materials from $^4$He to $^{208}$Pb. Monte Carlo simulations of low-energy backgrounds are presented for the liquid-neon-based Cryogenic Low Energy Astrophysics with Noble gases (CLEAN) detector. The simulated background is much smaller than the expected signal from a galactic supernova.

This paper is being withdrawn. The electic filed effects should be interpreted as a correction to the electron density. Therefore the original MSW potential calculated in terms of the full electron density is correct.

Neutron-star radii provide useful information on the equation of state of neutron rich matter. Particularly interesting is the density dependence of the equation of state (EOS). For example, the softening of the EOS at high density, where the pressure rises slower than anticipated, could signal a transition to an exotic phase. However, extracting the density dependence of the EOS requires measuring the radii of neutron stars for a broad range of masses. A ``normal'' 1.4 solar mass neutron star has a central density of a few times nuclear-matter saturation density. In contrast, low mass (of the order of 0.5 solar masses) neutron stars have central densities near nuclear-matter saturation density so its radius provides information on the EOS at low density. Unfortunately, low-mass stars are rare because they may be hard to form. Instead, a precision measurement of nuclear radii on atomic nuclei may contain similar information. Indeed, we find a strong correlation between the neutron radius of 208Pb and the radius of a 0.5 solar-mass neutron star. Thus, the radius of such a neutron star can be inferred from a measurement of the the neutron radius of 208Pb. Comparing this value to the measured radius of a 1.4 solar-mass neutron star should provide the strongest constraint to date on the density dependence of the equation of state.

Recent observations by the Chandra observatory suggest that some neutron stars may cool rapidly, perhaps by the direct URCA process which requires a high proton fraction. The proton fraction is determined by the nuclear symmetry energy whose density dependence may be constrained by measuring the neutron radius of a heavy nucleus, such as 208Pb. Such a measurement is necessary for a reliable extrapolation of the proton fraction to the higher densities present in the neutron star. A large neutron radius in 208Pb implies a stiff symmetry energy that grows rapidly with density, thereby favoring a high proton fraction and allowing direct URCA cooling. Predictions for the neutron radius in 208Pb are correlated to the proton fraction in dense matter by using a variety of relativistic effective field-theory models. Models that predict a neutron (Rn) minus proton (Rp) root-mean-square radius in 208Pb to be Rn-Rp<0.20 fm have proton fractions too small to allow the direct URCA cooling of 1.4 solar-mass neutron stars. Conversely, if Rn-Rp>0.25 fm, the direct URCA process is allowed (by all models) to cool down a 1.4 solar-mass neutron star. The Parity Radius Experiment at Jefferson Laboratory aims to measure the neutron radius in 208Pb accurately and model independently via parity-violating electron scattering. Such a measurement would greatly enhance our ability to either confirm or dismiss the direct URCA cooling of neutron stars.

The neutrino driven wind during a core collapse supernova is an attractive site for r-process nucleosynthesis. The electron fraction $Y_e$ in the wind depends on observable neutrino energies and luminosities. The mean antineutrino energy is limited by supernova SN1987A data while lepton number conservation constrains the ratio of antineutrino to neutrino luminosities. If $Y_e$, in the wind, is to be suitable for rapid neutron capture nucleosynthesis, then the mean electron neutrino energy may be significantly lower then that predicted in present supernova simulations, or there may be new neutrino physics such as oscillations to sterile neutrinos.

Parity violating electron scattering can measure the neutron density of a heavy nucleus accurately and model independently. This is because the weak charge of the neutron is much larger then that of the proton. The Parity Radius Experiment (PREX) at Jefferson Laboratory aims to measure the root mean square neutron radius of $^{208}$Pb with an absolute accuracy of 1% ($\pm 0.05$ Fm). This is more accurate then past measurements with hadronic probes, which all suffer from controversial strong interaction uncertainties. PREX should clearly resolve the neutron-rich skin. Furthermore, this benchmark value for $^{208}$Pb will provide a calibration for hadronic probes, such as proton scattering, which can then be used to measure neutron densities of many exotic nuclei. The PREX result will also have many implications for neutron stars. The neutron radius of Pb depends on the pressure of neutron-rich matter: the greater the pressure, the larger the radius as neutrons are pushed out against surface tension. The same pressure supports a neutron star against gravity. The Pb radius is sensitive to the equation of state at normal densities while the radius of a 1.4 solar mass neutron star also depends on the equation of state at higher densities. Measurements of the radii of a number of isolated neutron stars such as Geminga and RX J185635-3754 should soon improve significantly. By comparing the equation of state information from the radii of both Pb and neutron stars one can search for a softening of the high density equation of state from a phase transition to an exotic state. Possibilities include kaon condensates, strange quark matter or color superconductors.

Weak magnetism increases antineutrino mean free paths in core collapse supernovae. The parity violating interference between axial and vector currents makes antineutrino-nucleon cross sections smaller then those for neutrinos. We calculate simple, exact correction factors to include recoil and weak magnetism in supernova simulations. Weak magnetism may significantly increase the neutrino energy flux. We calculate, in a diffusion approximation, an increase of order 15% in the total energy flux for temperatures near 10 MeV. This should raise the neutrino luminosity. Weak magnetism also changes the emitted spectrum of $ \bar\nu_x$ (with $x=\mu$ or $\tau$) and $\bar\nu_e$. We estimate that $\bar\nu_x$ will be emitted about 7% hotter than $\nu_x$ because $\bar\nu_x$ have longer mean free paths. Likewise weak magnetism may increase the $\bar\nu_e$ temperature by of order 10%. This increase in temperature coupled with the increase in neutrino luminosity should increase the heating in the low density region outside of the neutrino sphere. This, in turn, could be important for the success of an explosion. It is important to check our results with a full simulation that includes Boltzmann neutrino transport and weak magnetism corrections.

A new relation between the neutron skin of a heavy nucleus and the radius of a neutron star is proposed: the larger the neutron skin of the nucleus the larger the radius of the star. Relativistic models that reproduce a variety of ground-state observables can not determine uniquely the neutron skin of a heavy nucleus. Thus, a large range of neutron skins is generated by supplementing the models with nonlinear couplings between isoscalar and isovector mesons. We illustrate how the correlation between the neutron skin and the radius of the star can be used to place important constraints on the equation of state and how it may help elucidate the existence of a phase transition in the interior of the neutron star.

We study relationships between the neutron-rich skin of a heavy nucleus and the properties of neutron-star crusts. Relativistic effective field theories with a thicker neutron skin in $^{208}$Pb have a larger electron fraction and a lower liquid-to-solid transition density for neutron-rich matter. These properties are determined by the density dependence of the symmetry energy which we vary by adding nonlinear couplings between isoscalar and isovector mesons. An accurate measurement of the neutron radius in $^{208}$Pb---via parity violating electron scattering---may have important implications for the structure of neutron stars.

Core collapse supernovae are dominated by energy transport from neutrinos. Therefore, some supernova properties could depend on symetries and features of the standard model weak interactions. The cross section for neutrino capture is larger than that for antineutrino capture by one term of order the neutrino energy over the nucleon mass. This reduces the ratio of neutrons to protons in the $\nu$-driven wind above a protoneutron star by approximately 20 % and may significantly hinder r-process nucleosynthesis.

The well known charge conjugation violating interactions in the Standard Model increase neutrino- and decrease anti-neutrino- nucleon cross sections. This impacts neutrino transport in core collapse supernovae through "recoil" corrections of order the neutrino energy $k$ over the nucleon mass $M$. All $k/M$ corrections to neutrino transport deep inside a protoneutron star are calculated from angular integrals of the Boltzmann equation. We find these corrections significantly modify neutrino currents at high temperatures. This produces a large mu and tau number for the protoneutron star and can change the ratio of neutrons to protons. In addition, the relative size of neutrino mean free paths changes. At high temperatures, the electron anti-neutrino mean free path becomes {\it longer} than that for mu or tau neutrinos.

High energy solar neutrinos can come from the hep reaction ^3He(p,\nue^+)^4He with a large end point energy of 18.8 MeV. Understanding the hep reaction may be important for interpreting solar neutrino spectra. We calculate the contribution of the axial charge transition $^3P_0\to ^1S_0$ to the hep thermonuclear S factor using a one-body reaction model involving a nucleon moving in optical potentials. Our result is comparable to or larger than previous calculations of the s-wave Gamow Teller contribution. This indicates that the hep reaction may have p-wave strength leading to an enhancement of the S factor.

Core collapse supernovae are dominated by energy transport from neutrinos. Therefore, some supernova properties could depend on symetries and features of the standard model weak interactions. The cross section for neutrino capture is larger than that for antineutrino capture by one term of order the neutrino energy over the nucleon mass. This reduces the ratio of neutrons to protons in the $\nu$-driven wind above a protoneutron star by approximately 20 % and may significantly hinder r-process nucleosynthesis.

The neutrino-nucleon cross section is slightly larger than that for $\bar\nu-N$. Therefore, $\bar\nu$ will escape more quickly from core collapse supernovae leaving the star $\nu$ rich. A diffusion formalism is used to calculate the time evolution of the mu and tau lepton number densities. These quickly reach steady state equilibria. We estimate that a protoneutron star with a maximum temperature near 50 MeV will contain over 50 % more $\nu_\mu$ and $\nu_\tau$ than $\bar\nu_\mu$ and $\bar\nu_\tau$. Supernovae may be the only known systems with large mu and or tau lepton numbers.

The potential enhancement of the cross section $\nu\bar\nu\to e^+e^-$ in the presence of a magnetic field is of critical interest for the study of supernovae and as a possible mechanism for gamma ray bursts. While parity violation(PV) in this reaction in free space is forbidden by CP, the presence of a CP non-invariant background gas of electrons creates an asymmetry in the cross section which may contribute to the asymmetry of a supernova and the natal velocities of neutron stars. We calculate the cross section in the presence of fields as high as $B = 10^{16}$G and find no significant enhancement of the cross section with magnetic field strength. By studying the systematics of our results as a function of environmental variables(field strength, density,temperature), we extrapolate the relative strength of the parity violating terms to the weak field limit, and find the parity violation is insufficient to provide the ``kick'' required to explain the observed velocities of neutron stars.

Supernovae provide a unique opportunity for large scale parity violation because they are dominated by neutrinos. We calculate the parity violating asymmetry A of neutrino emission in a strong magnetic field. We assume the neutrinos elastically scatter many times from slightly polarized neutrons. Because of multple interactions, A grows with the optical thickness of the proto-neutron star and may be much larger than previous estimates. As a result, the neutron star could recoil at a significant velocity.

Core collapse supernovae are dominated by weakly interacting neutrinos. This provides a unique opportunity for macroscopic parity violation. We speculate that parity violation in a strong magnetic field can lead to an asymmetry in the explosion and a recoil of the newly formed neutron star. We estimate the asymmetry from neutrino-polarized-neutron elastic scattering, polarized electron capture and neutrino-nucleus elastic scattering in a (partially) polarized electron gas.

Neutrino-nucleus elastic scattering is reduced in dense matter because of correlations between ions. The static structure factor for a plasma of electrons and ions is calculated from Monte Carlo simulations and parameterized with a least squares fit. Our results imply a large increase in the neutrino mean free path. This strongly limits the trapping of neutrinos in a supernova by coherent neutral current interactions.

Neutrino-nucleus cross sections for the detection of atmospheric neutrinos are calculated in relativistic impulse and random phase approximations. Pion production is estimated via inclusive Delta-hole nuclear excitations. The number of pion events can confuse the identification of the neutrino flavor. As a further source of pions we calculate coherent pion production (where the nucleus remains in the ground state). We examine how these nuclear structure effects influence atmospheric neutrino experiments and study possible improvements in the present detector simulations.

The cosmological neutrino background will mediate long range forces between objects. For a background of temperature T, the potential decreases as 1/r^5 for r >> 1/T and as 1/r for r << 1/T. These forces have large spin-dependent components. If the neutrino background is nonrelativistic, the long range forces are enhanced by a factor of the inverse neutrino velocity. These long range forces may provide a method for observing the cosmological neutrino background.