Recent observations of neutron stars with gravitational waves and X-ray timing provide unprecedented access to the equation of state (EoS) of cold dense matter at densities difficult to realize in terrestrial experiments. At the same time, predictions for the EoS with reliable uncertainty estimates from chiral effective field theory ($\chi$EFT) bound our theoretical ignorance. In this work, we analyze astrophysical data using a nonparametric representation of the neutron star EoS conditioned on $\chi$EFT to directly constrain the underlying physical properties of the compact objects. We discuss how the data alone constrain the EoS at high densities when we condition on $\chi$EFT at low densities. We also demonstrate how to exploit astrophysical data to directly test the predictions of $\chi$EFT for the EoS up to twice nuclear saturation density, and estimate the density at which these predictions might break down. We find that the existence of massive pulsars, gravitational waves from GW170817, and NICER observations of PSR J0030+0451 favor $\chi$EFT predictions for the EoS up to nuclear saturation density over a more agnostic analysis by a factor of 7 for the quantum Monte Carlo (QMC) calculations used in this work. While $\chi$EFT predictions using QMC are fully consistent with gravitational-wave data up to twice nuclear saturation, NICER observations suggest that the EoS stiffens relative to these predictions at nuclear saturation density. Additionally, we marginalize over the uncertainty in the density at which $\chi$EFT begins to break down, constraining the radius of a $1.4\,M_\odot$ neutron star to $R_{1.4} = 11.40^{+1.38}_{-1.04}$ ($12.54^{+0.71}_{-0.63}$) km and the pressure at twice nuclear saturation density to $p(2n_\mathrm{sat}) = 14.2^{+18.1}_{-8.4}$ ($28.7^{+15.3}_{-15.0}$) $\mathrm{MeV}/\mathrm{fm}^3$ with massive pulsar and gravitational-wave (and NICER) data.

We study the role of pions in hot dense matter encountered in astrophysics. We find that strong interactions enhance the number density of negatively charged pions and that this enhancement can be calculated reliably for a relevant range of density and temperature using the virial expansion. We assess the influence of pions and muons on the equation of state (EOS) and weak interaction rates in hot dense matter. We find that thermal pions increase the proton fraction and soften the EOS. We also find that charged current weak reactions involving pions and muons $\bar{\nu}_\mu+\mu^- \rightarrow \pi$ and $\nu_\mu+\pi^- \rightarrow \mu^-$ make an important contribution to the opacity of muon neutrinos. This could influence the dynamics of core-collapse supernovae and neutron star mergers. Finally, we note that pion-nucleon reactions can alter the evolution of the proton fraction when weak interactions are not in equilibrium.

The properties of neutron stars are determined by the nature of the matter that they contain. These properties can be constrained by measurements of the star's size. We obtain stringent constraints on neutron-star radii by combining multimessenger observations of the binary neutron-star merger GW170817 with nuclear theory that best accounts for density-dependent uncertainties in the equation of state. We construct equations of state constrained by chiral effective field theory and marginalize over these using the gravitational-wave observations. Combining this with the electromagnetic observations of the merger remnant that imply the presence of a short-lived hyper-massive neutron star, we find that the radius of a $1.4\,\rm{M}_\odot$ neutron star is $R_{1.4\,\mathrm{M}_\odot} = 11.0^{+0.9}_{-0.6}~{\rm km}$ (90% credible interval). Using this constraint, we show that neutron stars are unlikely to be disrupted in neutron-star black-hole mergers; subsequently, such events will not produce observable electromagnetic emission.

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.

Future gravitational-wave observations will enable unprecedented and unique science in extreme gravity and fundamental physics answering questions about the nature of dynamical spacetimes, the nature of dark matter and the nature of compact objects.

We developed a tip-tilt system to compensate the turbulence induced image motion for 1.3 m telescope at Vainu Bappu Observatory, Kavalur. The instrument is designed to operate at visible wavelength band (480-700 nm) with a field of view $1^{\prime}\times1^{\prime}$. The tilt corrected images have shown up to $\approx$ 57% improvement in image resolution and a corresponding peak intensity increase by a factor of $\approx$ 2.8. A closed-loop correction bandwidth of $\approx$ 26 Hz has been achieved with on-sky tests and the root mean square motion of the star image has been reduced by a factor of $\sim$ 14. These results are consistent with theoretical and numerical predictions of wave-front aberrations caused by atmospheric turbulence and image quality improvement expected from a real-time control system. In this paper, we present the details of the instrument design, laboratory calibration studies and quantify its performance on the telescope.

Multi-messenger observations of neutron star (NS) mergers have the potential to revolutionize nuclear astrophysics. They will improve our understanding of nucleosynthesis, provide insights about the equation of state (EOS) of strongly-interacting matter at high densities, and enable tests of the theory of gravity and of dark matter. Here, we focus on the EOS, where both gravitational waves (GWs) from neutron-star mergers and X-ray observations from space-based detectors such as NICER will provide more stringent constraints on the structure of neutron stars. Furthermore, recent advances in nuclear theory have enabled reliable calculations of the EOS at low densities using effective field theory based Hamiltonians and advanced techniques to solve the quantum many-body problem. In this paper, we address how the first observation of GWs from GW170817 can be combined with modern calculations of the EOS to extract useful insights about the EOS of matter encountered inside neutron stars. We analyze the impact of various uncertainties, the role of phase transitions in the NS core, and discuss how future observations will improve our understanding of dense matter.

Gravitational waves from neutron-star mergers are expected to provide stringent constraints on the structure of neutron stars. At the same time, recent advances in nuclear theory have enabled reliable calculations of the low density equation of state using effective field theory based Hamiltonians and advanced techniques to solve the quantum many-body problem. In this paper, we address how the first observation of gravitational waves from GW170817 can be combined with modern calculations of the equation of state to extract useful insights about the equation of state of matter encountered inside neutron stars. We analyze the impact of various uncertainties and we show that the tidal deformability extracted from GW170817 is compatible, while less constraining, than modern nuclear physics knowledge.

We show how observations of gravitational waves from binary neutron star (BNS) mergers over the next few years can be combined with insights from nuclear physics to obtain useful constraints on the equation of state (EoS) of dense matter, in particular, constraining the neutron-matter EoS to within 20% between one and two times the nuclear saturation density $n_0\approx 0.16\ {\text{fm}^{-3}}$. Using Fisher information methods, we combine observational constraints from simulated BNS merger events drawn from various population models with independent measurements of the neutron star radii expected from x-ray astronomy (the Neutron Star Interior Composition Explorer (NICER) observations in particular) to directly constrain nuclear physics parameters. To parameterize the nuclear EoS, we use a different approach, expanding from pure nuclear matter rather than from symmetric nuclear matter to make use of recent quantum Monte Carlo (QMC) calculations. This method eschews the need to invoke the so-called parabolic approximation to extrapolate from symmetric nuclear matter, allowing us to directly constrain the neutron-matter EoS. Using a principal component analysis, we identify the combination of parameters most tightly constrained by observational data. We discuss sensitivity to various effects such as different component masses through population-model sensitivity, phase transitions in the core EoS, and large deviations from the central parameter values.

Observations with next-generation ground-based detectors further enhanced with multi-messenger (electromagnetic and neutrino) detections will allow us to probe new extreme astrophysics. Target sources included: core-collapse supernovae, continuous emission from isolated or accreting neutron stars, and bursts from magnetars and other pulsars.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

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.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

We study the effect of a dibaryon, S, in the mass range 1860 MeV < m_S < 2054 MeV, which is heavy enough not to disturb the stability of nuclei and light enough to possibly be cosmologically metastable. Such a deeply bound state can act as a baryon sink in regions of high baryon density and temperature. We find that the ambient conditions encountered inside a newly born neutron star are likely to sustain a sufficient population of hyperons to ensure that a population of S dibaryons can equilibrate in less than a few seconds. This would be catastrophic for the stability of neutron stars and the observation of neutrino emission from the proto-neutron star of Supernova 1987A over ~ O(10)s. A deeply bound dibaryon is therefore incompatible with the observed supernova explosion, unless the cross section for S production is severely suppressed.

We report the atmospheric turbulence parameters namely, atmospheric seeing, the tilt-anisoplanatic angle ($\theta_0$) and the coherence time ($\tau_0$), measured under various sky conditions, at Vainu Bappu Observatory in Kavalur. Bursts of short exposure images of selected stars were recorded with a high-speed, frame-transfer CCD mounted on the Cassegrain focus of a newly commissioned 1.3 m telescope. The estimated median seeing is $ \approx 1.85^{\prime\prime} $ at wavelength of $\sim 600$ nm, the image motion correlation between different pairs of stars is $\sim 44\%$ for $\theta_0 \approx 36^{\prime\prime}$ and mean $\tau_0$ is $\approx 2.4$ ms. This work was motivated by the design considerations and expected the performance of an adaptive optics system that is currently being planned for the telescope.

The correlation of the tidal polarizabilities $\Lambda_1$-$\Lambda_2$ for GW170817 is predicted by combining dense-matter equations of state (EOSs) that satisfy nuclear physics constraints with the chirp mass and mass asymmetry for this event. Our models are constrained by calculations of the neutron-matter EOS using chiral effective field theory Hamiltonians with reliable error estimates up to once or twice the nuclear saturation density. In the latter case, we find that GW170817 does not improve our understanding of the EOS. We contrast two distinct extrapolations to higher density: a minimal model (MM) which assumes that the EOS is a smooth function of density described by a Taylor expansion and a more general model parametrized by the speed of sound that admits phase transitions. This allows us to identify regions in the $\Lambda_1$-$\Lambda_2$ plots that could favor the existence of new phases of matter in neutron stars. We predict the combined tidal polarizability of the two neutron stars in GW170817 to be $80\le \tilde{\Lambda}\le 580$ ($280\le \tilde{\Lambda}\le 480$ for the MM), which is smaller than the range suggested by the LIGO-Virgo data analysis. Our analysis also shows that GW170817 requires a NS with $M=1.4M_\odot$ to have a radius $9.0<R_{1.4}<13.6$~km ($ 11.3 <R_{1.4}< 13.6$~km for the MM).

If the thermal evolution of the hot young neutron star in the supernova remnant HESS J1731-347 is driven by neutrino emission, it provides a stringent constraint on the coupling of light (mass $\ll 10$ keV) axion-like particles to neutrons. Using Markov-Chain Monte Carlo we find that for the values of axion-neutron coupling $g_{ann}^2 > 7.7 \times 10^{-20}$ (90% c.l.) the axion cooling from the bremsstrahlung reaction $n+n\rightarrow n+n +a$ is too rapid to account for the high observed surface temperature. This implies that the Pecci-Quinn scale or axion decay constant $f_a > 6.7 \times 10^7$ GeV for KSVZ axions and $f_a > 1.7 \times 10^9$ GeV for DFSZ axions. The high temperature of this neutron star also allows us to tighten constraints on the size of the nucleon pairing gaps.

We present a comprehensive and pedagogic discussion of the properties of photons in cold and dense nuclear matter based on the resummed one-loop photon self energy. Correlations between electrons, muons, protons and neutrons in beta equilibrium that arise due to electromagnetic and strong interactions are consistently taken into account within the random phase approximation. Screening effects and damping as well as collective excitations are systematically studied in a fully relativistic setup. Our study is relevant to linear response theory of dense nuclear matter, calculations of transport properties of cold dense matter and to investigations of the production and propagation of hypothetical vector bosons such as the dark photons.

Exotic new particles carrying baryon number and with mass of order the nucleon mass have been proposed for various reasons including baryogenesis, dark matter, mirror worlds, and the neutron lifetime puzzle. We show that the existence of neutron stars with mass greater than 0.7 $M_\odot$ places severe constraints on such particles, requiring them to be heavier than 1.2 GeV or to have strongly repulsive self-interactions.

We calculate the differential scattering rate for thermal neutrinos in a hot and dilute gas of interacting neutrons using linear response theory. The dynamical structure factors for density and spin fluctuations of the strongly interacting neutron matter, expected in the neutrino decoupling regions of supernovae and neutron star mergers, are calculated in the virial expansion for the first time. Correlations due to nucleon-nucleon interactions are taken into account using a pseudo-potential that reproduces measured nucleon-nucleon phase shifts, and we find that attractive s-wave interactions enhance the density response and suppress the spin response of neutron matter. The net effect of neutron correlations is to strongly suppress backscattering. Moreover, we find nearly exact scaling laws for the response functions, valid for the range $T = 5 - 10$ MeV and q < 30 MeV, allowing us to obtain analytic results for the dynamic structure factors at second-order in the fugacity of the neutron gas. We find that the modification of scattering rates depends on the energy and momentum exchanged, implying that dynamical structure factors are essential to describe neutrino decoupling in supernovae and neutron star mergers.

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.

We find that the abundance of large clusters of nucleons in neutron-rich matter at sub-nuclear density is greatly reduced by finite temperature effects when matter is close to beta-equilibrium. Large nuclei and exotic non-spherical nuclear configurations called pasta, favored in the vicinity of the transition to uniform matter at $T=0$, dissolve at relatively low temperature. For matter close to beta-equilibrium we find that the pasta melting temperature is $T_m^\beta\simeq 4\pm 1$~MeV for realistic equations of state. The mechanism for pasta dissolution is discussed, and in general $T_m^\beta$ is shown to be sensitive to the proton fraction. We find that coherent neutrino scattering from nuclei and pasta makes a modest contribution to the opacity under the conditions encountered in supernovae and neutron star mergers. Implications for neutrino signals from galactic supernovae are briefly discussed.

We analyze the chiral magnetic effect in a homogeneous neutral plasma from the point of view of energy conservation, and construct an effective potential for the growth of maximally helical perturbations of the electromagnetic field. We show that a negative curvature at the origin of the potential, indicating instability of the plasma, is induced by a chiral asymmetry in electron Fermi energy, as opposed to number density, while the potential grows at large field value. It follows that the ground state for a plasma has zero magnetic helicity; a nonzero electron mass will allow an excited state of a plasma with nonzero helicity to relax to that ground state quickly. We conclude that a chiral plasma instability triggered by weak interactions is not a viable mechanism for explaining magnetic fields in stars except possibly when dynamics drives the system far from equilibrium.

We derive the charged current absorption rate of electron and anti-electron neutrinos in dense matter using a fully relativistic approach valid at arbitrary matter degeneracy. We include mean field energy shifts due to nuclear interactions and the corrections due to weak magnetism. The rates are derived both from the familiar Fermi's Golden Rule, and from the techniques of finite temperature field theory, and shown to be equivalent. In various limits, these results can also be used to calculate neutral current opacities. We find that some pieces of the response have been left out in previous derivations and their contribution at high density can be significant. Useful formulae and detailed derivations are presented and we provide a new open source implementation of these opacities for use in radiation hydrodynamic simulations of core-collapse supernovae and neutron star mergers.

After a successful core collapse supernova (CCSN) explosion, a hot dense proto-neutron star (PNS) is left as a remnant. Over a time of twenty or so seconds, this PNS emits the majority of the neutrinos that come from the CCSN, contracts, and loses most of its lepton number. This is the process by which all neutron stars in our galaxy are likely born. The emitted neutrinos were detected from SN 1987A and they will be detected in much greater numbers from any future galactic CCSN. These detections can provide a direct window into the properties of the dense matter encountered inside neutron stars, and they can affect nucleosynthesis in the mate- rial ejected during the CCSN. In this chapter, we review the basic physics of PNS cooling, including the basic equations of PNS structure and neutrino diffusion in dense matter. We then discuss how the nuclear equation of state, neutrino opacities in dense matter, and convection can shape the temporal behavior of the neutrino signal. We also discuss what was learned from the late time SN 1987A neutrinos, the prospects for detection of these neutrinos from future galactic CCSNe, and the effects these neutrinos can have on nucleosynthesis.

An accretion outburst onto a neutron star transient heats the neutron star's crust out of thermal equilibrium with the core. After the outburst the crust thermally relaxes toward equilibrium with the neutron star core and the surface thermal emission powers the quiescent X-ray light curve. Crust cooling models predict that thermal equilibrium of the crust will be established $\approx 1000 \, \mathrm{d}$ into quiescence. Recent observations of the cooling neutron star transient MXB 1659-29, however, suggest that the crust did not reach thermal equilibrium with the core on the predicted timescale and continued to cool after $\approx 2500 \, \mathrm{d}$ into quiescence. Because the quiescent light curve reveals successively deeper layers of the crust, the observed late time cooling of MXB 1659-29 depends on the thermal transport in the inner crust. In particular, the observed late time cooling is consistent with a low thermal conductivity layer near the depth predicted for nuclear pasta that maintains a temperature gradient between the neutron star's inner crust and core for thousands of days into quiescence. As a result, the temperature near the crust-core boundary remains above the critical temperature for neutron superfluidity and a layer of normal neutrons forms in the inner crust. We find that the late time cooling of MXB 1659-29 is consistent with heat release from a normal neutron layer near the crust-core boundary with a long thermal time.

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.

We calculate the thermal conductivity of electrons for the strongly correlated multi-component ion plasma expected in the outer layers of neutron star's crust employing a Path Integral Monte Carlo (PIMC) approach. This allows us to isolate the low energy response of the ions and use it to calculate the electron scattering rate and the electron thermal conductivity. We find that the scattering rate is enhanced by a factor 2-4 compared to earlier calculations based on the simpler electron-impurity scattering formalism. This findings directly impacts the interpretation of thermal relaxation observed in transiently accreting neutron stars and has implications for the composition and nuclear reactions in the crust that occur during accretion.

This white paper informs the nuclear astrophysics community and funding agencies about the scientific directions and priorities of the field and provides input from this community for the 2015 Nuclear Science Long Range Plan. It summarizes the outcome of the nuclear astrophysics town meeting that was held on August 21-23, 2014 in College Station at the campus of Texas A&M University in preparation of the NSAC Nuclear Science Long Range Plan. It also reflects the outcome of an earlier town meeting of the nuclear astrophysics community organized by the Joint Institute for Nuclear Astrophysics (JINA) on October 9- 10, 2012 Detroit, Michigan, with the purpose of developing a vision for nuclear astrophysics in light of the recent NRC decadal surveys in nuclear physics (NP2010) and astronomy (ASTRO2010). The white paper is furthermore informed by the town meeting of the Association of Research at University Nuclear Accelerators (ARUNA) that took place at the University of Notre Dame on June 12-13, 2014. In summary we find that nuclear astrophysics is a modern and vibrant field addressing fundamental science questions at the intersection of nuclear physics and astrophysics. These questions relate to the origin of the elements, the nuclear engines that drive life and death of stars, and the properties of dense matter. A broad range of nuclear accelerator facilities, astronomical observatories, theory efforts, and computational capabilities are needed. With the developments outlined in this white paper, answers to long standing key questions are well within reach in the coming decade.

We calculate the rate of production of hypothetical light vector bosons (LVBs) from nucleon-nucleon bremsstrahlung reactions in the soft radiation limit directly in terms of the measured nucleon-nucleon elastic cross sections. We use these results and the observation of neutrinos from supernova SN1987a to deduce constraints on the couplings of vector bosons with masses $\lesssim 200$ MeV to either electric charge (dark photons) or to baryon number. We establish for the first time strong constraints on LVB that couple only to baryon number, and revise earlier constraints on the dark photon. For the latter, we find that the excluded region of parameter space is diminished by about a factor of 10.

We use the quantum Monte Carlo (QMC) techniques to calculate the static structure function $S(q)$ of a one-component ion lattice and use it to calculate the thermal conductivity $\kappa$ of high-density solid matter expected in the neutron star crust. By making detailed comparisons with the results for the thermal conductivity obtained using standard techniques based on the one-phonon approximation (OPA) valid at low temperature, and the multi-phonon harmonic approximation expected to be valid over a wide range of temperatures, we asses the temperature regime where $S(q)$ from QMC can be used directly to calculate $\kappa$. We also compare the QMC results to those obtained using classical Monte Carlo to quantitatively asses the magnitude of the quantum corrections. We find that quantum effects became relevant for the calculation of $\kappa$ at temperature $T \lesssim 0.3 ~\Omega_\mathrm{P}$, where $\Omega_\mathrm{P}$ is the ion plasma frequency. At $T \simeq 0.1 ~\Omega_\mathrm{P}$ the quantum effects suppress $\kappa$ by about $30\%$. The comparison with the results of the OPA indicates that dynamical information beyond the static structure is needed when $T \lesssim 0.1~ \Omega_\mathrm{P}$. These quantitative comparisons help to establish QMC as a viable technique to calculate $\kappa$ at moderate temperatures in the range $T=0.1-1~\Omega_\mathrm{P}$ of relevance to the study of accreting neutron stars. This finding is especially important because QMC is the only viable technique so far for calculating $\kappa$ in multi-component systems at low-temperatures.

We calculate neutrino absorption rates due to charged-current reactions $\nu_e+n \rightarrow e^- + p $ and $\bar{\nu}_e+p \rightarrow e^+ + n $ in the outer regions of a newly born neutron star called the neutrino-sphere. To improve on recent work which has shown that nuclear mean fields enhance the $\nu_e$ cross-section and suppress the $\bar{\nu}_e$ cross-section, we employ realistic nucleon-nucleon interactions that fit measured scattering phase shifts. Using these interactions we calculate the momentum-, density-, and temperature-dependent nucleon self-energies in the Hartree-Fock approximation. A potential derived from chiral effective field theory and a pseudo-potential constructed to reproduce nucleon-nucleon phase shifts at the mean-field level are used to study the equilibrium proton fraction and the charged-current rates are studied in detail. We compare our results to earlier calculations obtained using phenomenological mean-field models and to those obtained in the virial expansion valid at low density. We find that for typical ambient conditions in the neutrino-sphere, $T=5-10$ MeV and $\rho =10^{11}-10^{13}$ g/cm$^3$, the difference between the $\nu_{e}$ and $\bar{\nu}_e$ absorption rates is much larger than predicted earlier. Our results have important implications for heavy-element nucleosynthesis in supernovae and for supernova neutrino detection.

We show that electrons can couple to the neutron excitations in neutron stars and find that this can limit their contribution to the transport properties of dense matter, especially the shear viscosity. The coupling between electrons and neutrons is induced by protons in the core, and by ions in the crust. We calculate the effective electron-neutron interaction for the kinematics of relevance to the scattering of degenerate electrons at high density. We use this interaction to calculate the electron thermal conductivity, electrical conductivity, and shear viscosity in the neutron star inner crust, and in the core where we consider both normal and superfluid phases of neutron-rich matter. In some cases, particularly when protons are superconducting and neutrons are in their normal phase, we find that electron-neutron scattering can be more important than the other scattering mechanisms considered previously.

We show that the nonzero electron mass plays a critical role in determining the magnetic properties of neutron stars, making it impossible to generate the chiral charge density needed to trigger a strong chiral magnetic instability during the core collapse of supernovae. This instability has been proposed as a plausible mechanism for generating extremely large helical magnetic fields in neutron stars at their birth; the mechanism relies on the generation of a large non-equilibrium chiral charge density via electron capture reactions that selectively deplete left-handed electrons during core-collapse and the early evolution of the protoneutron star. Our calculation shows that the electron chirality violation rate induced by Rutherford scattering, despite being suppressed by the smallness of the electron mass relative to the electron chemical potential, is still fast compared to the weak interaction electron capture rate. The resulting asymmetry between right and left-handed electron densities is therefore never able to attain an astrophysically relevant magnitude.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

The Facility for Rare Isotope Beams (FRIB) will be a world-leading laboratory for the study of nuclear structure, reactions and astrophysics. Experiments with intense beams of rare isotopes produced at FRIB will guide us toward a comprehensive description of nuclei, elucidate the origin of the elements in the cosmos, help provide an understanding of matter in neutron stars, and establish the scientific foundation for innovative applications of nuclear science to society. FRIB will be essential for gaining access to key regions of the nuclear chart, where the measured nuclear properties will challenge established concepts, and highlight shortcomings and needed modifications to current theory. Conversely, nuclear theory will play a critical role in providing the intellectual framework for the science at FRIB, and will provide invaluable guidance to FRIB's experimental programs. This article overviews the broad scope of the FRIB theory effort, which reaches beyond the traditional fields of nuclear structure and reactions, and nuclear astrophysics, to explore exciting interdisciplinary boundaries with other areas. \keywords{Nuclear Structure and Reactions. Nuclear Astrophysics. Fundamental Interactions. High Performance Computing. Rare Isotopes. Radioactive Beams.

We study how many-body effects alter the dark matter (DM) thermalization time inside neutron stars. We find that Pauli blocking, kinematic constraints, and superfluidity and superconductivity in the neutron star significantly affect the DM thermalization time, in general lengthening it. This could change the final DM mass and DM-nucleon cross section constraints by considering black hole formation in neutron stars due to DM accretion. We consider the class of models in which DM is an asymmetric, complex scalar particle with a mass between 1 keV and 5 GeV which couples to regular matter via some heavy vector boson. Interestingly, we find that the discovery of asymmetric, bosonic DM could motivate the existence of exotic neutron star cores. We apply our results to the case of mixed sneutrino DM.

We explore the spectrum of low-energy collective excitations in the crust of a neutron star, especially in the inner region where neutron-proton clusters are immersed in a sea of superfluid neutrons. The speeds of the different modes are calculated systematically from the nuclear energy density functional theory using a Skyrme functional fitted to essentially all experimental atomic mass data.

We formulate a theory of Goldstone bosons and their interactions in the superfluid and superconducting phase of dense nucleonic matter at densities of relevance to the neutron star core. For typical neutron star temperatures in the range T = 10^6 to 10^9 K, the Goldstone mode associated with rotational symmetry, called angulons, couple weakly to each other and to electrons. Consequently, these modes have anomalously large mean free paths and can contribute to both diffusive and ballistic transport of heat and momentum. In contrast, the two Goldstone bosons associated with density oscillations of the neutron and electron + proton fluids, called superfluid phonons, mix and couple strongly to electrons. They have shorter mean free paths, and their contribution to transport is negligible. Long-wavelength superfluid phonons and angulons can play a role in neutron star seismology, and lead to interesting phenomenology as angulons couple to magnetic fields and have anisotropic dispersion relations.

We review the calculation of the equation of state of pure neutron matter using quantum Monte Carlo (QMC) methods. QMC algorithms permit the study of many-body nuclear systems using realistic two- and three-body forces in a nonperturbative framework. We present the results for the equation of state of neutron matter, and focus on the role of three-neutron forces at supranuclear density. We discuss the correlation between the symmetry energy, the neutron star radius and the symmetry energy. We also combine QMC and theoretical models of the three-nucleon interactions, and recent neutron star observations to constrain the value of the symmetry energy and its density dependence.

It is now possible to model thermal relaxation of neutron stars after bouts of accretion during which the star is heated out of equilibrium by nuclear reactions in its crust. Major uncertainties in these models can be encapsulated in modest variations of a handful of fudge parameters that change the crustal thermal conductivity, specific heat, and heating rates. Observations of thermal relaxation constrain these fudge parameters and allow us to predict longer term variability in terms of the neutron star core temperature. We demonstrate this explicitly by modeling ongoing thermal relaxation in the neutron star XTE J1701-462. Its future cooling, over the next 5 to 30 years, is strongly constrained and depends mostly on its core temperature, uncertainties in crust physics having essentially been pinned down by fitting to the first three years of observations.

We study the low-energy collective excitations of the neutron star inner crust, where a neutron superfluid coexists with a Coulomb lattice of nuclei. The dispersion relation of the modes is calculated systematically from a microscopic theory including neutron band structure effects. These effects are shown to lead to a strong mixing between the Bogoliubov-Anderson bosons of the neutron superfluid and the longitudinal crystal lattice phonons. In addition, the speed of the transverse shear mode is greatly reduced as a large fraction of superfluid neutrons are entrained by nuclei. Not only does the much smaller velocity of the transverse mode increase the specific heat of the inner crust, but it also decreases its electron thermal conductivity. These results may impact our interpretation of the thermal relaxation in accreting neutron stars. Due to strong mixing, the mean free path of the superfluid mode is found to be greatly reduced. Our results for the collective mode dispersion relations and their damping may also have implications for neutron star seismology.

We study the spin response of cold dense neutron matter in the limit of zero momentum transfer, and show that the frequency dependence of the long-wavelength spin response is well constrained by sum-rules and the asymptotic behavior of the two-particle response at high frequency. The sum-rules are calculated using Auxiliary Field Diffusion Monte Carlo technique and the high frequency two-particle response is calculated for several nucleon-nucleon potentials. At nuclear saturation density, the sum-rules suggest that the strength of the spin response peaks at $\omega \simeq$ 40--60 MeV, decays rapidly for $\omega \geq $100 MeV, and has a sizable strength below 40 MeV. This strength at relatively low energy may lead to enhanced neutrino production rates in dense neutron-rich matter at temperatures of relevance to core-collapse supernova.

Previous work on neutrino emission from proto-neutron stars which employed full solutions of the Boltzmann equation showed that the average energies of emitted electron neutrinos and antineutrinos are closer to one another than predicted by older, more approximate work. This in turn implied that the neutrino driven wind is proton rich during its entire life, precluding $r$-process nucleosynthesis and the synthesis of Sr, Y, and Zr. This work relied on charged current neutrino interaction rates that are appropriate for a free nucleon gas. Here, it is shown in detail that the inclusion of the nucleon potential energies and collisional broadening of the response significantly alters this conclusion. Iso-vector interactions, which give rise to the nuclear symmetry energy, produce a difference between the neutron and proton single-particle energies $\Delta U=U_n-U_p$ and alter the kinematics of the charged current reactions. In neutron-rich matter, and for a given neutrino/antineutrino energy, the rate for $\nu_e+n\rightarrow e^-+p$ is enhanced while $ \bar{\nu}_e+p\rightarrow n+e^+$ is suppressed because the $Q$ value for these reactions is altered by $\pm\Delta U$, respectively. In the neutrino decoupling region, collisional broadening acts to enhance both $\nu_e$ and $\bar{\nu}_e$ cross-sections and RPA corrections decrease the $\nu_e$ cross-section and increase the $\bar \nu_e$ cross-section, but mean field shifts have a larger effect. Therefore, electron neutrinos decouple at lower temperature than when the nucleons are assumed to be free and have lower average energies. The change is large enough to allow for a reasonable period of time when the neutrino driven wind is predicted to be neutron rich. It is also shown that the electron fraction in the wind is influenced by the nuclear symmetry energy.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

If the surface tension of quark matter is low enough, quark matter is not self bound. At sufficiently low pressure and temperature, it will take the form of a crystal of positively charged strangelets in a neutralizing background of electrons. In this case there will exist, in addition to the usual family of strange stars, a family of low-mass large-radius objects analogous to white dwarfs, which we call "strangelet dwarfs". Using a generic parametrization of the equation of state of quark matter, we calculate the mass-radius relationship of these objects.

We calculate the equation of state of neutron matter with realistic two- and three-nucleon interactions using quantum Monte Carlo techniques, and illustrate that the short-range three-neutron interaction determines the correlation between neutron matter energy at nuclear saturation density and higher densities relevant to neutron stars. Our model also makes an experimentally testable prediction for the correlation between the nuclear symmetry energy and its density dependence -- determined solely by the strength of the short-range terms in the three neutron force. The same force provides a significant constraint on the maximum mass and radius of neutron stars.

We review the nuclear and condensed matter physics underlying the thermal and transport properties of the neutron star inner crust. These properties play a key role in interpreting transient phenomena such as thermal relaxation in accreting neutron stars, superbursts, and magnetar flares. We emphasize simplifications that occur at low temperature where the inner crust can be described in terms of electrons and collective excitations. The heat conductivity and heat capacity of the solid and superfluid phase of matter is discussed in detail and we emphasize its role in interpreting observations of neutron stars in soft X-ray transients. We highlight recent theoretical and observational results, and identify future work needed to better understand a host of transient phenomena in neutron stars.

We model neutrino emission from a newly born neutron star subsequent to a supernova explosion to study its sensitivity to the equation of state, neutrino opacities, and convective instabilities at high baryon density. We find the time period and spatial extent over which convection operates is sensitive to the behavior of the nuclear symmetry energy at and above nuclear density. When convection ends within the proto-neutron star, there is a break in the predicted neutrino emission that may be clearly observable.

In cold compact stars, Cooper pairing between fermions in dense matter leads to the formation of a gap in their excitation spectrum and typically exponentially suppresses transport properties. However, we show here that weak Urca reactions become strongly enhanced and approach their ungapped level when the star undergoes density oscillations of sufficiently large amplitude. We study both the neutrino emissivity and the bulk viscosity due to direct Urca processes in hadronic, hyperonic and quark matter and discuss different superfluid and superconducting pairing patterns.

We formulate a low energy effective theory describing phases of matter that are both solid and superfluid. These systems simultaneously break translational symmetry and the phase symmetry associated with particle number. The symmetries restrict the combinations of terms that can appear in the effective action and the lowest order terms featuring equal number of derivatives and Goldstone fields are completely specified by the thermodynamic free energy, or equivalently by the long-wavelength limit of static correlation functions in the ground state. We show that the underlying interaction between particles that constitute the lattice and the superfluid gives rise to entrainment, and mixing between the Goldstone modes. As a concrete example we discuss the low energy theory for the inner crust of a neutron star, where a lattice of ionized nuclei coexists with a neutron superfluid.

A simple description of superfluid hydrodynamics in the inner crust of a neutron star is given. Particular attention is paid to the effect of the lattice of nuclei on the properties of the superfluid neutrons, and the effects of entrainment, the fact that some fraction of the neutrons are locked to the motion of the protons in nuclei.

We calculate the screening of the ion-ion potential due to electrons in the presence of a large background magnetic field, at densities of relevance to neutron star crusts. Using the standard approach to incorporate electron screening through the one-loop polarization function, we show that the magnetic field produces important corrections both at short and long distances. In extreme fields, realized in highly magnetized neutron stars called magnetars, electrons occupy only the lowest Landau levels in the relatively low density region of the crust. Here our results show that the screening length for Coulomb interactions between ions can be smaller than the inter-ion spacing. More interestingly, we find that the screening is anisotropic and the screened potential between two static charges exhibits long range Friedel oscillations parallel to the magnetic field. This long-range oscillatory behavior is likely to affect the lattice structure of ions, and can possibly create rod-like structures in the magnetar crusts. We also calculate the imaginary part of the electron polarization function which determines the spectrum of electron-hole excitations and plays a role in damping lattice phonon excitations. We demonstrate that even for modest magnetic fields this damping is highly anisotropic and will likely lead to anisotropic phonon heat transport in the outer neutron star crust.

We report on a new mechanism for heat conduction in the neutron star crust. We find that collective modes of superfluid neutron matter, called superfluid phonons (sPhs), can influence heat conduction in magnetized neutron stars. They can dominate the heat conduction transverse to magnetic field when the magnetic field $B \gsim 10^{13}$ G. At density $\rho \simeq 10^{12}-10^{14} $ g/cm$^3$ the conductivity due to sPhs is significantly larger than that due to lattice phonons and is comparable to electron conductivity when temperature $\simeq 10^8$ K. This new mode of heat conduction can limit the surface anisotropy in highly magnetized neutron stars. Cooling curves of magnetized neutron stars with and without superfluid heat conduction could show observationally discernible differences.

{We calculate the neutrino emissivity of superfluid neutron matter in the inner crust of neutron stars. We find that neutrino emission due to fluctuations resulting from the formation of Cooper pairs at finite temperature is highly suppressed in non-relativistic systems. This suppression of the pair breaking emissivity in a simplified model of neutron matter with interactions that conserve spin is of the order of $v_F^4$ for density fluctuations and $v_F^2$ for spin fluctuations, where $v_F$ is the Fermi velocity of neutrons. The larger suppression of density fluctuations arises because the dipole moment of the density distribution of a single component system does not vary in time. For this reason, we find that the axial current response (spin fluctuations) dominates. In more realistic models of neutron matter which include tensor interactions where the neutron spin is not conserved, neutrino radiation from bremsstrahlung reactions occurs at order $v_F^0$. Consequently, even with the suppression factors due to superfluidity, this rate dominates near $T_C$. Present calculations of the pair-breaking emissivity are incomplete because they neglect the tensor component of the nucleon-nucleon interaction.

We calculate the bulk viscosity of color-superconducting quark matter in the color-flavor-locked (CFL) phase. We assume that the lightest bosons are the superfluid mode H and the kaons K0 and K+, and that there is no kaon condensate. We calculate the rate of strangeness-equilibrating processes that convert kaons into superfluid modes, and the resultant bulk viscosity. We find that for oscillations with a timescale of milliseconds, at temperatures T << 1 MeV, the CFL bulk viscosity is much less than that of unpaired quark matter, but at higher temperatures the bulk viscosity of CFL matter can become larger.

Compact relativistic stars allow us to study the nature of matter under extreme conditions, probing regions of parameter space that are otherwise inaccessible. Nuclear theory in this regime is not well constrained: one key issue is whether neutron stars are in fact composed primarily of strange quark matter. Distinguishing the two possibilities, however, has been difficult. The recent detection of seismic vibrations in the aftermath of giant flares from two magnetars (highly magnetized compact stars) is a major breakthrough. The oscillations excited seem likely to involve the stellar crust, the properties of which differ dramatically for strange stars. We show that the resulting mode frequencies cannot be reconciled with the observations for reasonable magnetar parameters. Ruling out strange star models would place a strong constraint on models of dense quark matter.

The existence of deconfined quark matter in the superdense interior of neutron stars is a key question that has drawn considerable attention over the past few decades. Quark matter can comprise an arbitrary fraction of the star, from 0 for a pure neutron star to 1 for a pure quark star, depending on the equation of state of matter at high density. From an astrophysical viewpoint, these two extreme cases are generally expected to manifest different observational signatures. An intermediate fraction implies a hybrid star, where the interior consists of mixed or homogeneous phases of quark and nuclear matter, depending on surface and Coulomb energy costs, as well as other finite size and screening effects. In this brief review article, we discuss what we can deduce about quark matter in neutron stars in light of recent exciting developments in neutron star observations. We state the theoretical ideas underlying the equation of state of dense quark matter, including color superconducting quark matter. We also highlight recent advances stemming from re-examination of an old paradigm for the surface structure of quark stars and discuss possible evolutionary scenarios from neutron stars to quark stars, with emphasis on astrophysical observations.

We review theoretical developments in studies of dense matter and its phase structure of relevance to compact stars. Observational data on compact stars, which can constrain the properties of dense matter, are presented critically and interpreted.

We construct strangelets, taking into account electrostatic effects, including Debye screening, and arbitrary surface tension sigma of the interface between vacuum and quark matter. We find that there is a critical surface tension sigma_crit below which large strangelets are unstable to fragmentation and below which quark star surfaces will fragment into a crystalline crust made of charged strangelets immersed in an electron gas. We derive a model-independent relationship between sigma_crit and two parameters that characterize any quark matter equation of state. For reasonable model equations of state, we find sigma_crit typically of order a few MeV/fm^2. If sigma <= sigma_crit, the size-distribution of strangelets in cosmic rays could feature a peak corresponding to the stable strangelets that we construct.

We reexamine the surface composition of strange stars. Strange quark stars are hypothetical compact stars which could exist if strange quark matter was absolutely stable. It is widely accepted that they are characterized by an enormous density gradient ($~10^{26}$ g/cm$^4$) and large electric fields at surface. By investigating the possibility of realizing a heterogeneous crust, comprised of nuggets of strange quark matter embedded in an uniform electron background, we find that the strange star surface has a much reduced density gradient and negligible electric field. We comment on how our findings will impact various proposed observable signatures for strange stars.

We show that a hybrid (nuclear + quark matter) star can have a mass-radius relationship very similar to that predicted for a star made of purely nucleonic matter. We show this for a generic parameterization of the quark matter equation of state, and also for an MIT bag model, each including a phenomenological correction based on gluonic corrections to the equation of state. We obtain hybrid stars as heavy as 2 M_solar for reasonable values of the bag model parameters. For nuclear matter, we use the equation of state calculated by Akmal, Pandharipande, and Ravenhall using many-body techniques. Both mixed and homogeneous phases of nuclear and quark matter are considered.

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

We calculate the correlation functions needed to describe the linear response of superfluid matter, and go on to calculate the differential cross section for neutral-current neutrino scattering in superfluid neutron matter and in color-flavor locked quark matter (CFL). We report the first calculation of scattering rates that includes neutrino interactions with both pair-breaking excitations and low-lying collective excitations (Goldstone modes). Our results apply both above and below the critical temperature, allowing use in simulations of neutrino transport in supernovae and neutron stars.

Neutrino-matter cross sections and interaction rates are central to the core-collapse supernova phenomenon and, very likely, to the viability of the explosion mechanism itself. In this paper, we describe the major neutrino scattering, absorption, and production processes that together influence the outcome of core collapse and the cooling of protoneutron stars. One focus is on energy redistribution and many-body physics, but our major goal is to provide a useful resource for those interested in supernova neutrino microphysics.

Calculating the structure of white dwarf and neutron stars would be a suitable topic for an undergraduate thesis or an advanced special topics or independent study course. The subject is rich in many different areas of physics accessible to a junior or senior physics major, ranging from thermodynamics to quantum statistics to nuclear physics to special and general relativity. The computations for solving the coupled structure differential equations (both Newtonian and general relativistic) can be done using a symbolic computational package, such as Mathematica. In doing so, the student will develop computational skills and learn how to deal with dimensions. Along the way he or she will also have learned some of the physics of equations of state and of degenerate stars.

We study weak interactions involving Goldstone bosons in the neutral kaon condensed phase of color flavor locked quark matter. We calculate the rates for the dominant processes that contribute to the neutrino mean free p ath and to neutrino production. A light $K^+$ state, with a mass $\tilde{m}_{K^+} \propto (\Delta/\mu) (\Delta/m_s)(m_d-m_u)$, where $\mu$ and $\Delta$ are the quark chemical potential and superconducting gap respectively, is shown to play an important role. We identify unique characteristics of weak interaction rates in this novel phase and discuss how they might influence neutrino emission in core collapse supernova and neutron stars.

Quark matter at astrophysical densities may contain stable vortices due to the spontaneous breaking of hypercharge symmetry by kaon condensation. We argue that these vortices could be both charged and electrically superconducting. Current carrying loops (vortons) could be long lived and play a role in the magnetic and transport properties of this matter. We provide a scenario for vorton formation in protoneutron stars.

We study compact stars that contain quark matter. We look at the effect of color superconductivity in the quark matter on the nuclear-quark matter transition density, mass-radius relationship, and the density discontinuity at the boundary between nuclear and quark matter. We find that color superconducting quark matter will occur in compact stars at values of the bag constant where ordinary quark matter would not be allowed. We are able to construct ``hybrid'' stars with a color superconducting quark matter interior and nuclear matter surface, with masses in the range 1.3-1.6 M_solar and radii 8-11 km. Our results are consistent with recent mass-radius limits based on absorption lines from EXO0748-676.

We analyze the phases of dense quark matter at nonzero quark masses. We map out the phase diagrams for charged quark matter as well as charge neutral matter containing leptons. We find that in Color Flavor Locked (CFL) quark matter, the symmetric phase with equal numbers of $u$, $d$ and $s$ quarks is unlikely to occur in nature. Various less symmetric phases with additional meson condensates may play a role in neutron star cores, however, in agreement with the analysis of Bedaque and Sch\"afer. We discuss the implications of these novel color superconducting quark matter phases and their transitions for neutron stars.

We investigate the consequences of enforcing local color neutrality on the color superconducting phases of quark matter by utilizing the Nambu-Jona-Lasinio model supplemented by diquark and the t'Hooft six-fermion interactions. In neutrino free matter at zero temperature, color neutrality guarantees that the number densities of u, d, and s quarks in the Color-Flavor-Locked (CFL) phase will be equal even with physical current quark masses. Electric charge neutrality follows as a consequence and without the presence of electrons. In contrast, electric charge neutrality in the less symmetric 2-flavor superconducting (2SC) phase with ud pairing requires more electrons than the normal quark phase. The free energy density cost of enforcing color and electric charge neutrality in the CFL phase is lower than that in the 2SC phase, which favors the formation of the CFL phase. With increasing temperature and neutrino content, an unlocking transition occurs from the CFL phase to the 2SC phase with the order of the transition depending on the temperature, the quark and lepton number chemical potentials. The astrophysical implications of this rich structure in the phase diagram, including estimates of the effects from Goldstone bosons in the CFL phase, are discussed.

We study the weak interaction rates involving Goldstone bosons in the Color Flavor Locked (CFL) quark matter. Neutrino mean free path and the rate of energy loss due to neutrino emission in a thermal plasma of CFL pions and kaons is calculated. We find that in addition to neutrino scattering off thermal mesons, novel Cherenkov like processes wherein mesons are either emitted or absorbed contribute to the neutrino opacity. Lack of Lorentz invariance in the medium and loss of rotational invariance for processes involving mesons moving relative to the medium allow for novel processes such as $\pi^0 \to \nu \bar{\nu}$ and $e^- \pi^+ \to \nu_e$. We explore and comment on various astrophysical implications of our finding.

At nuclear matter density, electrically neutral strongly interacting matter in weak equilibrium is made of neutrons, protons and electrons. At sufficiently high density, such matter is made of up, down and strange quarks in the color-flavor locked phase, with no electrons. As a function of increasing density (or, perhaps, increasing depth in a compact star) other phases may intervene between these two phases which are guaranteed to be present. The simplest possibility, however, is a single first order phase transition between CFL and nuclear matter. Such a transition, in space, could take place either through a mixed phase region or at a single sharp interface with electron-free CFL and electron-rich nuclear matter in stable contact. Here we construct a model for such an interface. It is characterized by a region of separated charge, similar to an inversion layer at a metal-insulator boundary. On the CFL side, the charged boundary layer is dominated by a condensate of negative kaons. We then consider the energetics of the mixed phase alternative. We find that the mixed phase will occur only if the nuclear-CFL surface tension is significantly smaller than dimensional analysis would indicate.

The existence of compact dimensions which are accessible only to gravity represents an intriguing possible solution to the hierarchy problem. At present the strongest constraint on the existence of such compact Gravity-Only Dimensions (GODs) comes from SN1987a. Here we report on the first self-consistent simulations of the early, neutrino-emitting phase of a proto-neutron star which include energy losses due to the coupling of the Kaluza-Klein modes of the graviton which arise in a world with GODs. We compare the neutrino signals from these simulations to that from SN1987a and use a rigorous probabilistic analysis to derive improved bounds for the radii of such GODs. We find that the possibility that there are two compact extra dimensions with radii larger than 0.66 $\mu$m is excluded at the 95% confidence level---as is the possibility that there are three compact extra dimensions larger than 0.8 nm.

The main stages in the evolution of a neutron star, from its birth as a proto-neutron star, to its old age as a cold, catalyzed configuration, are described. A proto-neutron star is formed in the aftermath of a successful supernova explosion and its evolution is dominated by neutrino diffusion. Its neutrino signal is a valuable diagnostic of its internal structure and composition. During its transformation from a hot, lepton-rich to a cold, catalyzed remnant, the possibility exists that it can collapse into a black hole, which abruptly terminates neutrino emissions. The essential microphysics, reviewed herein, that controls its evolution are the equation of state of dense matter and its associated neutrino opacities. Several simulations of the proto-neutron star evolution, involving different assumptions about the composition of dense matter, are described. After its evolution into a nearly isothermal neutron star a hundred or so years after its birth, it may be observable through its thermal emission in X-rays during its life in the next million years. Its surface temperature will depend upon the rapidity of neutrino emission processes in its core, which depends on the composition of dense matter and whether or not its constituents exhibit superfluidity and superconductivity. Observations of thermal emission offer the best hope of a determination of the radius of a neutron star. The implications for the underlying dense matter equation of state of an accurate radius determination are explored.

We study the role of Coulomb and surface effects on the phase transition from dense nuclear matter to a mixed phase of nuclear and kaon-condensed matter. We calculate corrections to the bulk calculation of the equation of state (EOS) and the critical density for the transition by solving explicitly for spherical, cylindrical, and planar structures. The importance of Debye screening in the determination of the charged particle profiles is studied in some detail. We find that the surface and Coulomb contributions to the energy density are small, but that they play an important role in the determination of the critical pressure for the transition, as well as affecting the size and geometry of favored structures. This changes the EOS over a wide range of pressure and consequently increases the maximum mass by about 0.1 solar masses. Implications for transport properties of the mixed phase are also discussed.

We study the equation of state (EOS) of kaon-condensed matter including the effects of temperature and trapped neutrinos. It is found that the order of the phase transition to a kaon-condensed phase, and whether or not Gibbs' rules for phase equilibrium can be satisfied in the case of a first order transition, depend sensitively on the choice of the kaon-nucleon interaction. The main effect of finite temperature, for any value of the lepton fraction, is to mute the effects of a first order transition, so that the thermodynamics becomes similar to that of a second order transition. Above a critical temperature, found to be at least 30--60 MeV depending upon the interaction, the first order transition disappears. The phase boundaries in baryon density versus lepton number and baryon density versus temperature planes are delineated. We find that the thermal effects on the maximum gravitational mass of neutron stars are as important as the effects of trapped neutrinos, in contrast to previously studied cases in which the matter contained only nucleons or in which hyperons and/or quark matter were considered. Kaon-condensed EOSs permit the existence of metastable neutron stars, because the maximum mass of an initially hot, lepton-rich protoneutron star is greater than that of a cold, deleptonized neutron star. The large thermal effects imply that a metastable protoneutron star's collapse to a black hole could occur much later than in previously studied cases that allow metastable configurations.

One of the strongest constraints on the existence of large, compact, "gravity-only" dimensions comes from SN1987a. If the rate of energy loss into these putative extra dimensions is too high, then the neutrino pulse from the supernova will differ from that actually seen. The dominant mechanism for the production of Kaluza-Klein gravitons and dilatons in the supernova is via gravistrahlung and dilastrahlung from the nucleon-nucleon system. In this paper we compute the rates for these processes in a model-independent way using low-energy theorems which relate the emissivities to the measured nucleon-nucleon cross section. This is possible because for soft gravitons and dilatons the leading contribution to the energy-loss rate is from graphs in which the gravitational radiation is produced from external nucleon legs. Previous calculations neglected these mechanisms. We re-evaluate the bounds on toroidally-compactified "gravity-only" dimensions (GODs), and find that consistency with the observed SN1987a neutrino signal requires that if there are two such dimensions then their radius must be less than 1 micron.

Neutrino and axion production in neutron stars occurs mainly as bremsstrahlung from nucleon-nucleon (NN) scattering. The energy radiated via neutrinos or axions is typically very small compared to other scales in the two-nucleon system. The rate of emission of such "soft" radiation is directly related to the on-shell NN amplitude, and thereby to the NN experimental data. This facilitates the model-independent calculation of the neutrino and axion radiation rates which is presented here. We find that the resultant rates are a factor of four below earlier estimates based on a one-pion-exchange NN amplitude.

We calculate the neutrino mean free path in color superconducting quark matter, and employ it to study the cooling of matter via neutrino diffusion in the superconducting phase as compared to a free quark phase. The cooling process slows when quark matter undergoes a second order phase transition to a superconducting phase at the critical temperature $T_c$. Cooling subsequently accelerates as the temperature decreases below $T_c$. This will directly impact the early evolution of a newly born neutron star should its core contain quark matter. Consequently, there may be observable changes in the early neutrino emission which would provide evidence for superconductivity in hot and dense matter.

A first order phase transition at high baryon density implies that a mixed phase can occupy a significant region of the interior of a neutron star. In this article we investigate the effect of a droplet phase on neutrino transport inside the core. Two specific scenarios of the phase transition are examined, one having a kaon condensate and the other having quark matter in the high density phase. The coherent scattering of neutrinos off the droplets greatly increases the neutrino opacity of the mixed phase. We comment on how the existence of such a phase will affect a supernova neutrino signal.

An extensive study of the effects of correlations on both charged and neutral current weak interaction rates in dense matter is performed. Both strong and electromagnetic correlations are considered.The propagation of particle-hole interactions in the medium plays an important role in determining the neutrino mean free paths. The effects due to Pauli-Blocking and density, spin, and isospin correlations in the medium significantly reduce the neutrino cross sections. Due to the lack of experimental information at high density, these correlations are necessarily model dependent. For example, spin correlations in nonrelativistic models are found to lead to larger suppressions of neutrino cross sections compared to those of relativistic models. This is due to the tendency of the nonrelativistic models to develop spin instabilities. Notwithstanding the above caveats, and the differences between nonrelativistic and relativistic approaches such as the spin- and isospin-dependent interactions and the nucleon effective masses, suppressions of order 2--3, relative to the case in which correlations are ignored, are obtained. Neutrino interactions in dense matter are especially important for supernova and early neutron star evolution calculations. The effects of correlations for protoneutron star evolution are calculated. Large effects on the internal thermodynamic properties of protoneutron stars, such as the temperature, are found. These translate into significant early enhancements in the emitted neutrino energies and fluxes, especially after a few seconds. At late times, beyond about 10 seconds, the emitted neutrino fluxes decrease more rapidly compared to simulations without the effects of correlations, due to the more rapid onset of neutrino transparency in the protoneutron star

We study the thermal and chemical evolution during the Kelvin-Helmholtz phase of the birth of a neutron star, employing neutrino opacities that are consistently calculated with the underlying equation of state (EOS). Expressions for the diffusion coefficients appropriate for general relativistic neutrino transport in the equilibrium diffusion approximation are derived. The diffusion coefficients are evaluated using a field-theoretical finite temperature EOS that includes the possible presence of hyperons. The variation of the diffusion coefficients is studied as a function of EOS and compositional parameters. We present results from numerical simulations of protoneutron star cooling for internal stellar properties as well as emitted neutrino energies and luminosities. We discuss the influence of the initial stellar model, the total mass, the underlying EOS, and the addition of hyperons on the evolution of the protoneutron star and upon the expected signal in terrestrial detectors.

We study the diffusive transport of neutrinos in a newly born neutron star to explore its sensitivity to dense matter properties. Energy and lepton number which are trapped during the catastrophic implosion diffuse out on the time scale of a few tens of seconds. Results for different dense matter models are presented.

The early evolution of a protoneutron star depends both on the equation of state and neutrino interactions at high density. We identify the important sources of neutrino opacity and through model calculations show that in-medium effects play an important role in determining the neutrino mean free paths. The effects due to Pauli-blocking and many-body correlations due to strong interactions reduce the neutrino cross sections by large factors compared to the case in which these effects are ignored. We discuss these results in the context of neutrino transport in a protoneutron star.

We study the charged and neutral current weak interaction rates relevant for the determination of neutrino opacities in dense matter found in supernovae and neutron stars. We establish an efficient formalism for calculating differential cross sections and mean free paths for interacting, asymmetric nuclear matter at arbitrary degeneracy. The formalism is valid for both charged and neutral current reactions. Strong interaction corrections are incorporated through the in-medium single particle energies at the relevant density and temperature. The effects of strong interactions on the weak interaction rates are investigated using both potential and effective field-theoretical models of matter. We investigate the relative importance of charged and neutral currents for different astrophysical situations, and also examine the influence of strangeness-bearing hyperons. Our findings show that the mean free paths are significantly altered by the effects of strong interactions and the multi-component nature of dense matter. The opacities are then discussed in the context of the evolution of the core of a protoneutron star.

We calculate neutrino cross sections from neutral current reactions in the dense matter encountered in the evolution of a newly born neutron star. Effects of composition and of strong interactions in the deleptonization and cooling phases of the evolution are studied. The influence of the possible presence of strangeness-rich hyperons on the neutrino scattering cross sections is explored. Due to the large vector couplings of the Sigma-minus and Cascade-minus, |C_V|~2, these particles, if present in protoneutron star matter, give significant contributions to neutrino scattering. In the deleptonization phase, the presence of strangeness leads to large neutrino energies, which results in large enhancements in the cross sections compared to those in matter with nucleons only. In the cooling phase, in which matter is nearly neutrino free, the response of the Sigma-minus hyperons to thermal neutrinos is the most significant. Neutrinos couple relatively weakly to the Lambda hyperons and, hence, their contributions are significant only at high density.

A protoneutron star is formed immediately after the gravitational collapse of the core of a massive star. At birth, the hot and high density matter in such a star contains a large number of neutrinos trapped during collapse. Trapped neutrinos generally inhibit the presence of exotic matter -- hyperons, a kaon condensate, or quarks. However, as the neutrinos diffuse out in about 10-15 s, the threshold for the appearance of strangeness is reduced; hence, the composition and the structure of the star can change significantly. The effect of exotic, negatively-charged, strangeness-bearing components is always to soften the equation of state, and the possibility exists that the star collapses to a black hole at this time. This could explain why no neutron star has yet been seen in the remnant of supernova SN1987A, even though one certainly existed when neutrinos were detected on Feb. 23, 1987. With new generation neutrino detectors it is feasible to test different theoretical scenarios observationally.

We calculate neutrino cross sections from neutral current reactions in dense matter containing hyperons. We show that Sigma- hyperons give significant contributions. To lowest order, the contributions from the neutral Lambda and Sigma0, which have zero hyper-charge, vanish. However, their presence furnishes baryon number which decreases the relative concentrations of nucleons. This leads to significant reductions in the cross sections. Due to the uncertainty in strong interactions at high density, the neutrino opacity may vary by a factor of about 2 depending on the behavior of the effective masses.