We report results from searches for new physics with low-energy electronic recoil data recorded with the XENON1T detector, with an exposure of 0.65 t-y and an unprecedentedly low background rate of $76 \pm 2_{stat}$ events/(t y keV) between 1-30 keV. An excess over known backgrounds is observed below 7 keV, rising towards lower energies and prominent between 2-3 keV. The solar axion model has a 3.5$\sigma$ significance, and a three-dimensional 90% confidence surface is reported for axion couplings to electrons, photons, and nucleons. This surface is inscribed in the cuboid defined by $g_{ae} < 3.7 \times 10^{-12}$, $g_{ae}g_{an}^{eff} < 4.6 \times 10^{-18}$, and $g_{ae}g_{a\gamma} < 7.6\times10^{-22}~{GeV}^{-1}$, and excludes either $g_{ae}=0$ or $g_{ae}g_{a\gamma}=g_{ae}g_{an}^{eff}=0$. The neutrino magnetic moment signal is similarly favored over background at 3.2$\sigma$ and a confidence interval of $\mu_{\nu} \in (1.4,~2.9)\times10^{-11}\mu_B$ (90% C.L.) is reported. Both results are in strong tension with stellar constraints. The excess can also be explained by $\beta$ decays of tritium at 3.2$\sigma$ significance with a corresponding tritium concentration in xenon of $(6.2 \pm 2.0) \times 10^{-25}$ mol/mol. Such a trace amount can be neither confirmed nor excluded with current knowledge of production and reduction mechanisms. The significances of the solar axion and neutrino magnetic moment hypotheses are decreased to 2.1$\sigma$ and 0.9$\sigma$, respectively, if an unconstrained tritium component is included in the fitting. With respect to bosonic dark matter, the excess favors a monoenergetic peak at ($2.3 \pm 0.2$) keV (68% C.L.) with a 3.0$\sigma$ global (4.0$\sigma$ local) significance over background. This analysis sets the most restrictive direct constraints to date on pseudoscalar and vector bosonic dark matter for most masses between 1 and 210 keV/c$^2$.

We detail the sensitivity of the liquid xenon (LXe) DARWIN observatory to solar neutrinos via elastic electron scattering. We find that DARWIN will have the potential to measure the fluxes of five solar neutrino components: $pp$, $^7$Be, $^{13}$N, $^{15}$O and $pep$. The precision of the $^{13}$N, $^{15}$O and $pep$ components is hindered by the double-beta decay of $^{136}$Xe and, thus, would benefit from a depleted target. A high-statistics observation of $pp$ neutrinos would allow us to infer the values of the weak mixing angle, $\sin^2\theta_w$, and the electron-type neutrino survival probability, $P_e$, in the electron recoil energy region from a few keV up to 200 keV for the first time, with relative precision of 5% and 4%, respectively, at an exposure of 300 ty. An observation of $pp$ and $^7$Be neutrinos would constrain the neutrino-inferred solar luminosity down to 0.2%. A combination of all flux measurements would distinguish between the high (GS98) and low metallicity (AGS09) solar models with 2.1-2.5$\sigma$ significance, independent of external measurements from other experiments or a measurement of $^8$B neutrinos through coherent elastic neutrino-nucleus scattering in DARWIN. Finally, we demonstrate that with a depleted target DARWIN may be sensitive to the neutrino capture process of $^{131}$Xe.

Xenon dual-phase time projection chambers designed to search for Weakly Interacting Massive Particles have so far shown a relative energy resolution which degrades with energy above $\sim$200 keV. This has limited their sensitivity in the search for rare events like the neutrinoless double-beta decay of $^{136}$Xe at its $Q$-value, $Q_{\beta\beta}\simeq$ 2.46 MeV. For the XENON1T dual-phase time projection chamber, we demonstrate that the relative energy resolution at 1 $\sigma/\mu$ is as low as (0.79$\pm$0.02) % in its one-ton fiducial mass, and for single-site interactions at $Q_{\beta\beta}$. We achieve this by using a signal correction method to rectify the saturation effects of the signal readout system. The very good energy resolution from keV to MeV energies demonstrated in XENON1T opens up new windows for the xenon dual-phase dark matter detectors to simultaneously search for other rare events.

Direct dark matter detection experiments based on a liquid xenon target are leading the search for dark matter particles with masses above $\sim$ 5 GeV/c$^2$, but have limited sensitivity to lighter masses because of the small momentum transfer in dark matter-nucleus elastic scattering. However, there is an irreducible contribution from inelastic processes accompanying the elastic scattering, which leads to the excitation and ionization of the recoiling atom (the Migdal effect) or the emission of a Bremsstrahlung photon. In this letter, we report on a probe of low-mass dark matter with masses down to about 85 MeV/c$^2$ by looking for electronic recoils induced by the Migdal effect and Bremsstrahlung, using data from the XENON1T experiment. Besides the approach of detecting both scintillation and ionization signals, we exploit an approach that uses ionization signals only, which allows for a lower detection threshold. This analysis significantly enhances the sensitivity of XENON1T to light dark matter previously beyond its reach.

We report constraints on light dark matter (DM) models using ionization signals in the XENON1T experiment. We mitigate backgrounds with strong event selections, rather than requiring a scintillation signal, leaving an effective exposure of $(22 \pm 3)$ tonne-days. Above $\sim\!0.4\,\mathrm{keV}_\mathrm{ee}$, we observe $<1 \, \text{event}/(\text{tonne} \times \text{day} \times \text{keV}_\text{ee})$, which is more than one thousand times lower than in similar searches with other detectors. Despite observing a higher rate at lower energies, no DM or CEvNS detection may be claimed because we cannot model all of our backgrounds. We thus exclude new regions in the parameter spaces for DM-nucleus scattering for DM masses $m_\chi$ within $3-6\,\mathrm{GeV}/\mathrm{c}^2$, DM-electron scattering for $m_\chi > 30\,\mathrm{MeV}/\mathrm{c}^2$, and absorption of dark photons and axion-like particles for $m_\chi$ within $0.186 - 1 \, \mathrm{keV}/\mathrm{c}^2$.

In this paper we report on the characterization of the Hamamatsu VUV4 (S/N: S13370-6152) Vacuum Ultra-Violet (VUV) sensitive Silicon Photo-Multipliers (SiPMs) as part of the development of a solution for the detection of liquid xenon scintillation light for the nEXO experiment. Various SiPM features, such as: dark noise, gain, correlated avalanches, direct crosstalk and Photon Detection Efficiency (PDE) were measured in a dedicated setup at TRIUMF. SiPMs were characterized in the range $163 \text{ } \text{K} \leq \text{T}\leq 233 \text{ } \text{K}$. At an over voltage of $3.1\pm0.2$ V and at $\text{T}=163 \text{ }\text{K}$ we report a number of Correlated Avalanches (CAs) per pulse in the $1 \upmu\text{s}$ interval following the trigger pulse of $0.161\pm0.005$. At the same settings the Dark-Noise (DN) rate is $0.137\pm0.002 \text{ Hz/mm}^{2}$. Both the number of CAs and the DN rate are within nEXO specifications. The PDE of the Hamamatsu VUV4 was measured for two different devices at $\text{T}=233 \text{ }\text{K}$ for a mean wavelength of $189\pm7\text{ nm}$. At $3.6\pm0.2$ V and $3.5\pm0.2$ V of over voltage we report a PDE of $13.4\pm2.6\text{ }\%$ and $11\pm2\%$, corresponding to a saturation PDE of $14.8\pm2.8\text{ }\%$ and $12.2\pm2.3\%$, respectively. Both values are well below the $24\text{ }\%$ saturation PDE advertised by Hamamatsu. More generally, the second device tested at $3.5\pm0.2$ V of over voltage is below the nEXO PDE requirement. The first one instead yields a PDE that is marginally close to meeting the nEXO specifications. This suggests that with modest improvements the Hamamatsu VUV4 MPPCs could be considered as an alternative to the FBK-LF SiPMs for the final design of the nEXO detector.

Neutron stars cooling after sustained accretion outbursts provide unique information about the neutron star crust and underlying dense matter. Comparisons between astronomical observations of these cooling transients and model calculations of neutron star crust cooling have frequently been used to constrain neutron star properties such as the mass, radius, crust composition, and presence of nuclear pasta. These calculations often use a fixed pressure at which the crust-core transition happens, though this quantity depends on the dense matter equation of state. We demonstrate that varying the crust-core transition pressure in a manner consistent with adopting various equations of state results in modest changes to the crust cooling light curve. This validates the approach adopted in most crust cooling studies to date, where the neutron star mass and radius are varied while leaving the crust-core transition pressure constant.

We report the first experimental results on spin-dependent elastic weakly interacting massive particle (WIMP) nucleon scattering from the XENON1T dark matter search experiment. The analysis uses the full ton year exposure of XENON1T to constrain the spin-dependent proton-only and neutron-only cases. No significant signal excess is observed, and a profile likelihood ratio analysis is used to set exclusion limits on the WIMP-nucleon interactions. This includes the most stringent constraint to date on the WIMP-neutron cross section, with a minimum of $6.3\times10^{-42}$ cm$^2$ at 30 GeV/c${}^2$ and 90% confidence level. The results are compared with those from collider searches and used to exclude new parameter space in an isoscalar theory with an axial-vector mediator.

We present first results on the scalar WIMP-pion coupling from 1 t$\times$yr of exposure with the XENON1T experiment. This interaction is generated when the WIMP couples to a virtual pion exchanged between the nucleons in a nucleus. In contrast to most non-relativistic operators, these pion-exchange currents can be coherently enhanced by the total number of nucleons, and therefore may dominate in scenarios where spin-independent WIMP-nucleon interactions are suppressed. Moreover, for natural values of the couplings, they dominate over the spin-dependent channel due to their coherence in the nucleus. Using the signal model of this new WIMP-pion channel, no significant excess is found, leading to an upper limit cross section of $6.4\times10^{-46}$ cm$^2$ (90 % confidence level) at 30 GeV/c$^2$ WIMP mass.

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.

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

Monitoring the cooling of neutron-star crusts heated during accretion outbursts allows us to infer the physics of the dense matter present in the crust. We examine the crust cooling evolution of the low-mass X-ray binary MXB 1659-29 up to ~505 days after the end of its 2015 outburst (hereafter outburst II) and compare it with what we observed after its previous 1999 outburst (hereafter outburst I) using data obtained from the Swift, XMM-Newton, and Chandra observatories. The observed effective surface temperature of the neutron star in MXB 1659-29 dropped from ~92 eV to ~56 eV from ~12 days to ~505 days after the end of outburst II. The most recently performed observation after outburst II suggests that the crust is close to returning to thermal equilibrium with the core. We model the crust heating and cooling for both its outbursts collectively to understand the effect of parameters that may change for every outburst (e.g., the average accretion rate, the length of outburst, the envelope composition of the neutron star at the end of the outburst) and those which can be assumed to remain the same during these two outbursts (e.g., the neutron star mass, its radius). Our modelling indicates that all parameters were consistent between the two outbursts with no need for any significant changes. In particular, the strength and the depth of the shallow heating mechanism at work (in the crust) were inferred to be the same during both outbursts, contrary to what has been found when modelling the cooling curves after multiple outburst of another source, MAXI J0556-332. This difference in source behaviour is not understood. We discuss our results in the context of our current understanding of cooling of accretion-heated neutron-star crusts, and in particular with respect to the unexplained shallow heating mechanism.

We report on a search for Weakly Interacting Massive Particles (WIMPs) using 278.8 days of data collected with the XENON1T experiment at LNGS. XENON1T utilizes a liquid xenon time projection chamber with a fiducial mass of $(1.30 \pm 0.01)$ t, resulting in a 1.0 t$\times$yr exposure. The energy region of interest, [1.4, 10.6] $\mathrm{keV_{ee}}$ ([4.9, 40.9] $\mathrm{keV_{nr}}$), exhibits an ultra-low electron recoil background rate of $(82\substack{+5 \\ -3}\textrm{ (sys)}\pm3\textrm{ (stat)})$ events/$(\mathrm{t}\times\mathrm{yr}\times\mathrm{keV_{ee}})$. No significant excess over background is found and a profile likelihood analysis parameterized in spatial and energy dimensions excludes new parameter space for the WIMP-nucleon spin-independent elastic scatter cross-section for WIMP masses above 6 GeV/c${}^2$, with a minimum of $4.1\times10^{-47}$ cm$^2$ at 30 GeV/c${}^2$ and 90% confidence level.

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.

X-ray observations of transiently accreting neutron stars during quiescence provide information about the structure of neutron star crusts and the properties of dense matter. Interpretation of the observational data requires an understanding of the nuclear reactions that heat and cool the crust during accretion, and define its nonequilibrium composition. We identify here in detail the typical nuclear reaction sequences down to a depth in the inner crust where the mass density is 2E12 g/cm^3 using a full nuclear reaction network for a range of initial compositions. The reaction sequences differ substantially from previous work. We find a robust reduction of crust impurity at the transition to the inner crust regardless of initial composition, though shell effects can delay the formation of a pure crust somewhat to densities beyond 2E12 g/cm^3. This naturally explains the small inner crust impurity inferred from observations of a broad range of systems. The exception are initial compositions with A >= 102 nuclei, where the inner crust remains impure with an impurity parameter of Qimp~20 due to the N = 82 shell closure. In agreement with previous work we find that nuclear heating is relatively robust and independent of initial composition, while cooling via nuclear Urca cycles in the outer crust depends strongly on initial composition. This work forms a basis for future studies of the sensitivity of crust models to nuclear physics and provides profiles of composition for realistic crust models.

In this paper, we describe the XENON100 data analyses used to assess the target-intrinsic background sources radon ($^{222}$Rn), thoron ($^{220}$Rn) and krypton ($^{85}$Kr). We detail the event selections of high-energy alpha particles and decay-specific delayed coincidences. We derive distributions of the individual radionuclides inside the detector and quantify their abundances during the main three science runs of the experiment over a period of $\sim$ 4 years, from January 2010 to January 2014. We compare our results to external measurements of radon emanation and krypton concentrations where we find good agreement. We report an observed reduction in concentrations of radon daughters that we attribute to the plating-out of charged ions on the negatively biased cathode.

We report on the response of liquid xenon to low energy electronic recoils below 15 keV from beta decays of tritium at drift fields of 92 V/cm, 154 V/cm and 366 V/cm using the XENON100 detector. A data-to-simulation fitting method based on Markov Chain Monte Carlo is used to extract the photon yields and recombination fluctuations from the experimental data. The photon yields measured at the two lower fields are in agreement with those from literature; additional measurements at a higher field of 366 V/cm are presented. The electronic and nuclear recoil discrimination as well as its dependence on the drift field and photon detection efficiency are investigated at these low energies. The results provide new measurements in the energy region of interest for dark matter searches using liquid xenon.

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.

The transient neutron star (NS) low-mass X-ray binary MAXI J0556$-$332 provides a rare opportunity to study NS crust heating and subsequent cooling for multiple outbursts of the same source. We examine {\it MAXI}, {\it Swift}, {\it Chandra}, and {\it XMM-Newton} data of MAXI J0556$-$332 obtained during and after three accretion outbursts of different durations and brightness. We report on new data obtained after outburst III. The source has been tracked up to $\sim$1800 d after the end of outburst I. Outburst I heated the crust strongly, but no significant reheating was observed during outburst II. Cooling from $\sim$333 eV to $\sim$146 eV was observed during the first $\sim$1200 d. Outburst III reheated the crust up to $\sim$167 eV, after which the crust cooled again to $\sim$131 eV in $\sim$350 d. We model the thermal evolution of the crust and find that this source required a different strength and depth of shallow heating during each of the three outbursts. The shallow heating released during outburst I was $\sim$17 MeV nucleon$^{-1}$ and outburst III required $\sim$0.3 MeV nucleon$^{-1}$. These cooling observations could not be explained without shallow heating. The shallow heating for outburst II was not well constrained and could vary from $\sim$0--2.2 MeV nucleon$^{-1}$, i.e., this outburst could in principle be explained without invoking shallow heating. We discuss the nature of the shallow heating and why it may occur at different strengths and depths during different outbursts.

We report the first dark matter search results from XENON1T, a $\sim$2000-kg-target-mass dual-phase (liquid-gas) xenon time projection chamber in operation at the Laboratori Nazionali del Gran Sasso in Italy and the first ton-scale detector of this kind. The blinded search used 34.2 live days of data acquired between November 2016 and January 2017. Inside the (1042$\pm$12) kg fiducial mass and in the [5, 40] $\mathrm{keV}_{\mathrm{nr}}$ energy range of interest for WIMP dark matter searches, the electronic recoil background was $(1.93 \pm 0.25) \times 10^{-4}$ events/(kg $\times$ day $\times \mathrm{keV}_{\mathrm{ee}}$), the lowest ever achieved in a dark matter detector. A profile likelihood analysis shows that the data is consistent with the background-only hypothesis. We derive the most stringent exclusion limits on the spin-independent WIMP-nucleon interaction cross section for WIMP masses above 10 GeV/c${}^2$, with a minimum of 7.7 $\times 10^{-47}$ cm${}^2$ for 35-GeV/c${}^2$ WIMPs at 90% confidence level.

We present the first search for dark matter-induced delayed coincidence signals in a dual-phase xenon time projection chamber, using the 224.6 live days of the XENON100 science run II. This very distinct signature is predicted in the framework of magnetic inelastic dark matter which has been proposed to reconcile the modulation signal reported by the DAMA/LIBRA collaboration with the null results from other direct detection experiments. No candidate event has been found in the region of interest and upper limits on the WIMP's magnetic dipole moment are derived. The scenarios proposed to explain the DAMA/LIBRA modulation signal by magnetic inelastic dark matter interactions of WIMPs with masses of 58.0 GeV/c$^2$ and 122.7 GeV/c$^2$ are excluded at 3.3 $\sigma$ and 9.3 $\sigma$, respectively.

We present results of searches for vector and pseudo-scalar bosonic super-WIMPs, which are dark matter candidates with masses at the keV-scale, with the XENON100 experiment. XENON100 is a dual-phase xenon time projection chamber operated at the Laboratori Nazionali del Gran Sasso. A profile likelihood analysis of data with an exposure of 224.6 live days $\times$ 34\,kg showed no evidence for a signal above the expected background. We thus obtain new and stringent upper limits in the $(8-125)$\,keV/c$^2$ mass range, excluding couplings to electrons with coupling constants of $g_{ae} > 3\times10^{-13}$ for pseudo-scalar and $\alpha'/\alpha > 2\times10^{-28}$ for vector super-WIMPs, respectively. These limits are derived under the assumption that super-WIMPs constitute all of the dark matter in our galaxy.

The XENON1T experiment at the Laboratori Nazionali del Gran Sasso (LNGS) is the first WIMP dark matter detector operating with a liquid xenon target mass above the ton-scale. Out of its 3.2t liquid xenon inventory, 2.0t constitute the active target of the dual-phase time projection chamber. The scintillation and ionization signals from particle interactions are detected with low-background photomultipliers. This article describes the XENON1T instrument and its subsystems as well as strategies to achieve an unprecedented low background level. First results on the detector response and the performance of the subsystems are also presented.

We report on WIMP search results in the XENON100 detector using a non-relativistic effective field theory approach. The data from science run II (34 kg $\times$ 224.6 live days) was re-analyzed, with an increased recoil energy interval compared to previous analyses, ranging from $(6.6 - 240)~\mathrm{keV_\mathrm{nr}}$. The data is found to be compatible with the background-only hypothesis. We present 90% confidence level exclusion limits on the coupling constants of WIMP-nucleon effective operators using a binned profile likelihood method. We also consider the case of inelastic WIMP scattering, where incident WIMPs may up-scatter to a higher mass state, and set exclusion limits on this model as well.

We present the first constraints on the spin-dependent, inelastic scattering cross section of Weakly Interacting Massive Particles (WIMPs) on nucleons from XENON100 data with an exposure of 7.64$\times$10$^3$\,kg\,day. XENON100 is a dual-phase xenon time projection chamber with 62\,kg of active mass, operated at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy and designed to search for nuclear recoils from WIMP-nucleus interactions. Here we explore inelastic scattering, where a transition to a low-lying excited nuclear state of $^{129}$Xe is induced. The experimental signature is a nuclear recoil observed together with the prompt de-excitation photon. We see no evidence for such inelastic WIMP-$^{129}$Xe interactions. A profile likelihood analysis allows us to set a 90\% C.L. upper limit on the inelastic, spin-dependent WIMP-nucleon cross section of $3.3 \times 10^{-38}$\,cm$^{2}$ at 100\,GeV/c$^2$. This is the most constraining result to date, and sets the pathway for an analysis of this interaction channel in upcoming, larger dual-phase xenon detectors.

We report on WIMP search results of the XENON100 experiment, combining three runs summing up to 477 live days from January 2010 to January 2014. Data from the first two runs were already published. A blind analysis was applied to the last run recorded between April 2013 and January 2014 prior to combining the results. The ultra-low electromagnetic background of the experiment, ~$5 \times 10^{-3}$ events/(keV$_{\mathrm{ee}}\times$kg$\times$day) before electronic recoil rejection, together with the increased exposure of 48 kg $\times$ yr improves the sensitivity. A profile likelihood analysis using an energy range of (6.6 - 43.3) keV$_{\mathrm{nr}}$ sets a limit on the elastic, spin-independent WIMP-nucleon scattering cross section for WIMP masses above 8 GeV/$c^2$, with a minimum of 1.1 $\times 10^{-45}$ cm$^2$ at 50 GeV/$c^2$ and 90% confidence level. We also report updated constraints on the elastic, spin-dependent WIMP-nucleon cross sections obtained with the same data. We set upper limits on the WIMP-neutron (proton) cross section with a minimum of 2.0 $\times 10^{-40}$ cm$^2$ (52$\times 10^{-40}$ cm$^2$) at a WIMP mass of 50 GeV/$c^2$, at 90% confidence level.

We report on a search for electronic recoil event rate modulation signatures in the XENON100 data accumulated over a period of 4 years, from January 2010 to January 2014. A profile likelihood method, which incorporates the stability of the XENON100 detector and the known electronic recoil background model, is used to quantify the significance of periodicity in the time distribution of events. There is a weak modulation signature at a period of $431^{+16}_{-14}$ days in the low energy region of $(2.0-5.8)$ keV in the single scatter event sample, with a global significance of $1.9\,\sigma$, however no other more significant modulation is observed. The expected annual modulation of a dark matter signal is not compatible with this result. Single scatter events in the low energy region are thus used to exclude the DAMA/LIBRA annual modulation as being due to dark matter electron interactions via axial vector coupling at $5.7\,\sigma$.

We present the first results of searches for axions and axion-like-particles with the XENON100 experiment. The axion-electron coupling constant, $g_{Ae}$, has been probed by exploiting the axio-electric effect in liquid xenon. A profile likelihood analysis of 224.6 live days $\times$ 34 kg exposure has shown no evidence for a signal. By rejecting $g_{Ae}$, larger than $7.7 \times 10^{-12}$ (90\% CL) in the solar axion search, we set the best limit to date on this coupling. In the frame of the DFSZ and KSVZ models, we exclude QCD axions heavier than 0.3 eV/c$^2$ and 80 eV/c$^2$, respectively. For axion-like-particles, under the assumption that they constitute the whole abundance of dark matter in our galaxy, we constrain $g_{Ae}$, to be lower than $1 \times 10^{-12}$ (90\% CL) for mass range from 1 to 40 keV/c$^2$, and set the best limit to date as well.

In AM CVn binaries, a white dwarf primary accretes material from a helium-rich white dwarf or stellar companion. The unstable ignition of nuclear burning via the $3\alpha$ reaction in an accumulated helium layer powers a thermonuclear runaway near accretion rates $\dot{M} \lesssim 10^{-6} \, \mathrm{M_{\odot} \ yr^{-1}}$ that may be observed as helium nova or .Ia supernova. Helium burning in the primary's envelope at temperatures $T \lesssim 10^{8} \, \mathrm{K}$ may proceed via the direct three-body fusion of $\alpha$-particles. Here we show that the direct three-body rate by Nguyen et al. (2012) -- which is reduced relative to the extrapolated resonant rate at temperatures $T \gtrsim 5 \times 10^{7} \, \mathrm{K}$ -- results in novae with longer recurrence times and larger ignition masses. By contrast, we find that the enhancement in the direct three-body rate at temperatures below $T \lesssim 5 \times 10^{7} \, \mathrm{K}$ does not result in significant differences in nova outburst properties. The most massive envelopes in our models are near the density threshold for detonation of the helium layer, where an increase in the density at ignition due to the $3\alpha$ rate may be important.

An accretion outburst onto a neutron star deposits hydrogen-rich and/or helium-rich material into the neutron star's envelope. Thermonuclear burning of accreted material robustly produces Urca pairs\textrm{---}pairs of nuclei that undergo cycles of $e^-$-capture and $\beta ^-$-decay. The strong $T^5$ dependence of the Urca cooling neutrino luminosity means that Urca pairs in the neutron star interior potentially remove heat from accretion-driven nuclear reactions. In this study, we identify Urca pairs in the neutron star's ocean \textrm{---} a plasma of ions and electrons overlaying the neutron star crust \textrm{---} and demonstrate that Urca cooling occurs at all depths in the ocean. We find that Urca pairs in the ocean and crust lower the ocean's steady-state temperature during an accretion outburst and that unstable carbon ignition, which is thought to trigger superbursts, occurs deeper than it would otherwise. Cooling superburst light curves, however, are only marginally impacted by cooling from Urca pairs because the superburst peak radiative luminosity $L_{\rm peak}$ is always much greater than the Urca pair neutrino luminosity $L_{\nu}$ in the hot post-superburst ocean.

Crustal cooling of accretion-heated neutron stars provides insight into the stellar interior of neutron stars. The neutron star X-ray transient, KS~1731$-$260, was in outburst for 12.5 years before returning to quiescence in 2001. We have monitored the cooling of this source since then through {\it Chandra} and {\it XMM-Newton} observations. Here, we present a 150 ks {\it Chandra} observation of KS~1731$-$260 taken in August 2015, about 14.5 years into quiescence, and 6 years after the previous observation. We find that the neutron star surface temperature is consistent with the previous observation, suggesting that crustal cooling has likely stopped and the crust has reached thermal equilibrium with the core. Using a theoretical crust thermal evolution code, we fit the observed cooling curves and constrain the core temperature (T$_c = 9.35\pm0.25\times10^7$ K), composition (Q$_{imp} = 4.4^{+2.2}_{-0.5}$) and level of extra shallow heating required (Q$_{sh} = 1.36\pm0.18$ MeV/nucleon). We find that the presence of a low thermal conductivity layer, as expected from nuclear pasta, is not required to fit the cooling curve well, but cannot be excluded either.

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.

DARk matter WImp search with liquid xenoN (DARWIN) will be an experiment for the direct detection of dark matter using a multi-ton liquid xenon time projection chamber at its core. Its primary goal will be to explore the experimentally accessible parameter space for Weakly Interacting Massive Particles (WIMPs) in a wide mass-range, until neutrino interactions with the target become an irreducible background. The prompt scintillation light and the charge signals induced by particle interactions in the xenon will be observed by VUV sensitive, ultra-low background photosensors. Besides its excellent sensitivity to WIMPs above a mass of 5 GeV/c2, such a detector with its large mass, low-energy threshold and ultra-low background level will also be sensitive to other rare interactions. It will search for solar axions, galactic axion-like particles and the neutrinoless double-beta decay of 136-Xe, as well as measure the low-energy solar neutrino flux with <1% precision, observe coherent neutrino-nucleus interactions, and detect galactic supernovae. We present the concept of the DARWIN detector and discuss its physics reach, the main sources of backgrounds and the ongoing detector design and R&D efforts.

We characterize two 40 kBq sources of electrodeposited Th-228 for use in low-background experiments. The sources efficiently emanate Rn-220, a noble gas that can diffuse in a detector volume. Rn-220 and its daughter isotopes produce alpha, beta, and gamma-radiation, which may used to calibrate a variety of detector responses and features, before decaying completely in only a few days. We perform various tests to place limits on the release of other long-lived isotopes. In particular, we find an emanation of <0.008 atoms/min/kBq (90% CL) for Th-228 and 1.53 atoms/min/kBq for Ra-224. The sources lend themselves in particular to the calibration of detectors employing liquid noble elements such as argon and xenon. With the source mounted in a noble gas system, we demonstrate that filters are highly efficient in reducing the activity of these longer-lived isotopes further. We thus confirm the suitability of these sources even for use in next-generation experiments, such as XENON1T/XENONnT, LZ, and nEXO.

We refine our previously introduced parameterized model for explosive carbon-oxygen fusion during thermonuclear supernovae (SN Ia) by adding corrections to post-processing of recorded Lagrangian fluid element histories to obtain more accurate isotopic yields. Deflagration and detonation products are verified for propagation in a uniform density medium. A new method is introduced for reconstructing the temperature-density history within the artificially thick model deflagration front. We obtain better than 5\% consistency between the electron capture computed by the burning model and yields from post-processing. For detonations, we compare to a benchmark calculation of the structure of driven steady-state planar detonations performed with a large nuclear reaction network and error-controlled integration. We verify that, for steady-state planar detonations down to a density of 5x10^6 g/cc, our post processing matches the major abundances in the benchmark solution typically to better than 10% for times greater than 0.01 s after the shock front passage. As a test case to demonstrate the method, presented here with post-processing for the first time, we perform a two dimensional simulation of a SN Ia in the Chandrasekhar-mass deflagration-detonation transition (DDT) scenario. We find that reconstruction of deflagration tracks leads to slightly more complete silicon burning than without reconstruction. The resulting abundance structure of the ejecta is consistent with inferences from spectroscopic studies of observed SNe Ia. We confirm the absence of a central region of stable Fe-group material for the multi-dimensional DDT scenario. Detailed isotopic yields are tabulated and only change modestly when using deflagration reconstruction.

In Type Ia Supernovae (\sneia), the relative abundances of chemical elements are affected by the neutron excess in the composition of the progenitor white dwarf. Since these products leave signatures in the spectra near maximum light, spectral features may be used to constrain the composition of the progenitor. We calculate the nucleosynthetic yields for three \snia simulations, assuming single degenerate, Chandrasekhar mass progenitors, for a wide range of progenitor metallicities, and calculate synthetic light curves and spectra to explore correlations between progenitor metallicity and the strength of spectral features. We use two 2D simulations of the deflagration-detonation-transition scenario with different $^{56}$Ni yields and the W7 simulation to control for differences between explosion models and total yields. While the overall yields of intermediate mass elements (16 $<$ A $\leq$ 40) differ between the three cases, trends in the yields are similar. With increasing metallicity, $^{28}$Si yields remain nearly constant, $^{40}$Ca yields decline, and Ti and $^{54}$Fe yields increase. In the synthetic spectra, we identify two features at 30 days post explosion that appear to deepen with progenitor metallicity: a Ti feature around 4200\,\AA\ and a Fe feature around 5200\,\AA\@. In all three simulations, their pseudo equivalent widths show a systematic trend with progenitor metallicity. This suggests that these two features may allow differentiation among progenitor metallicities of observed \sneia and potentially help reduce the intrinsic Hubble scatter.

The XENON1T experiment is currently in the commissioning phase at the Laboratori Nazionali del Gran Sasso, Italy. In this article we study the experiment's expected sensitivity to the spin-independent WIMP-nucleon interaction cross section, based on Monte Carlo predictions of the electronic and nuclear recoil backgrounds. The total electronic recoil background in $1$ tonne fiducial volume and ($1$, $12$) keV electronic recoil equivalent energy region, before applying any selection to discriminate between electronic and nuclear recoils, is $(1.80 \pm 0.15) \cdot 10^{-4}$ ($\rm{kg} \cdot day \cdot keV)^{-1}$, mainly due to the decay of $^{222}\rm{Rn}$ daughters inside the xenon target. The nuclear recoil background in the corresponding nuclear recoil equivalent energy region ($4$, $50$) keV, is composed of $(0.6 \pm 0.1)$ ($\rm{t} \cdot y)^{-1}$ from radiogenic neutrons, $(1.8 \pm 0.3) \cdot 10^{-2}$ ($\rm{t} \cdot y)^{-1}$ from coherent scattering of neutrinos, and less than $0.01$ ($\rm{t} \cdot y)^{-1}$ from muon-induced neutrons. The sensitivity of XENON1T is calculated with the Profile Likelihood Ratio method, after converting the deposited energy of electronic and nuclear recoils into the scintillation and ionization signals seen in the detector. We take into account the systematic uncertainties on the photon and electron emission model, and on the estimation of the backgrounds, treated as nuisance parameters. The main contribution comes from the relative scintillation efficiency $\mathcal{L}_\mathrm{eff}$, which affects both the signal from WIMPs and the nuclear recoil backgrounds. After a $2$ y measurement in $1$ t fiducial volume, the sensitivity reaches a minimum cross section of $1.6 \cdot 10^{-47}$ cm$^2$ at m$_\chi$=$50$ GeV/$c^2$.

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.

Aql X-1 is a prolific transient neutron star low-mass X-ray binary that exhibits an accretion outburst approximately once every year. Whether the thermal X-rays detected in intervening quiescent episodes are the result of cooling of the neutron star or due to continued low-level accretion remains unclear. In this work we use Swift data obtained after the long and bright 2011 and 2013 outbursts, as well as the short and faint 2015 outburst, to investigate the hypothesis that cooling of the accretion-heated neutron star crust dominates the quiescent thermal emission in Aql X-1. We demonstrate that the X-ray light curves and measured neutron star surface temperatures are consistent with the expectations of the crust cooling paradigm. By using a thermal evolution code, we find that ~1.2-3.2 MeV/nucleon of shallow heat release describes the observational data well, depending on the assumed mass-accretion rate and temperature of the stellar core. We find no evidence for varying strengths of this shallow heating after different outbursts, but this could be due to limitations of the data. We argue that monitoring Aql X-1 for up to ~1 year after future outbursts can be a powerful tool to break model degeneracies and solve open questions about the magnitude, depth and origin of shallow heating in neutron star crusts.

We investigate constraints on neutron star structure arising from the assumptions that neutron stars have crusts, that recent calculations of pure neutron matter limit the equation of state of neutron star matter near the nuclear saturation density, that the high-density equation of state is limited by causality and the largest high-accuracy neutron star mass measurement, and that general relativity is the correct theory of gravity. We explore the role of prior assumptions by considering two classes of equation of state models. In a first, the intermediate- and high-density behavior of the equation of state is parameterized by piecewise polytropes. In the second class, the high-density behavior of the equation of state is parameterized by piecewise continuous line segments. The smallest density at which high-density matter appears is varied in order to allow for strong phase transitions above the nuclear saturation density. We critically examine correlations among the pressure of matter, radii, maximum masses, the binding energy, the moment of inertia, and the tidal deformability, paying special attention to the sensitivity of these correlations to prior assumptions about the equation of state. It is possible to constrain the radii of $1.4~\mathrm{M}_{\odot}$ neutron stars to a be larger than 10 km, even without consideration of additional astrophysical observations, for example, those from photospheric radius expansion bursts or quiescent low-mass X-ray binaries. We are able to improve the accuracy of known correlations between the moment of inertia and compactness as well as the binding energy and compactness. We also demonstrate the existence of a correlation between the neutron star binding energy and the moment of inertia.

We have searched for periodic variations of the electronic recoil event rate in the (2-6) keV energy range recorded between February 2011 and March 2012 with the XENON100 detector, adding up to 224.6 live days in total. Following a detailed study to establish the stability of the detector and its background contributions during this run, we performed an un-binned profile likelihood analysis to identify any periodicity up to 500 days. We find a global significance of less than 1 sigma for all periods suggesting no statistically significant modulation in the data. While the local significance for an annual modulation is 2.8 sigma, the analysis of a multiple-scatter control sample and the phase of the modulation disfavor a dark matter interpretation. The DAMA/LIBRA annual modulation interpreted as a dark matter signature with axial-vector coupling of WIMPs to electrons is excluded at 4.8 sigma.

Laboratory experiments searching for galactic dark matter particles scattering off nuclei have so far not been able to establish a discovery. We use data from the XENON100 experiment to search for dark matter interacting with electrons. With no evidence for a signal above the low background of our experiment, we exclude a variety of representative dark matter models that would induce electronic recoils. For axial-vector couplings to electrons, we exclude cross-sections above 6x10^(-35) cm^2 for particle masses of m_chi = 2 GeV/c^2. Independent of the dark matter halo, we exclude leptophilic models as explanation for the long-standing DAMA/LIBRA signal, such as couplings to electrons through axial-vector interactions at a 4.4 sigma confidence level, mirror dark matter at 3.6 sigma, and luminous dark matter at 4.6 sigma.

An accretion outburst in an X-ray transient deposits material onto the neutron star primary; this accumulation of matter induces reactions in the neutron star's crust. During the accretion outburst these reactions heat the crust out of thermal equilibrium with the core. When accretion halts, the crust cools to its long-term equilibrium temperature on observable timescales. Here we examine the accreting neutron star transient MAXI J0556-332, which is the hottest transient, at the start of quiescence, observed to date. Models of the quiescent light curve require a large deposition of heat in the shallow outer crust from an unknown source. The additional heat injected is $\approx 4\textrm{-}10\,\mathrm{MeV}$ per accreted nucleon; when the observed decline in accretion rate at the end of the outburst is accounted for, the required heating increases to $\approx 6\textrm{-}16\,\mathrm{MeV}$. This shallow heating is still required to fit the lightcurve even after taking into account a second accretion episode, uncertainties in distance, and different surface gravities. The amount of shallow heating is larger than that inferred for other neutron star transients and is larger than can be supplied by nuclear reactions or compositionally driven convection; but it is consistent with stored mechanical energy in the accretion disk. The high crust temperature ($T_b \gtrsim 10^{9} \, {\rm K}$) makes its cooling behavior in quiescence largely independent of the crust composition and envelope properties, so that future observations will probe the gravity of the source. Fits to the lightcurve disfavor the presence of Urca cooling pairs in the crust.

When neutron stars reside in transient X-ray binaries, their crustal layers become heated during accretion outbursts and subsequently cool in quiescence. Observing and modeling this thermal response has yielded valuable insight into the physics of neutron star crusts. However, one unresolved problem is the evidence in several neutron stars for an extra energy source, located at shallow depth in the crust, that is not accounted for by standard heating models. Its origin remains puzzling, and it is currently unclear whether this additional heating occurs in all neutron stars, and if the magnitude is always the same. Here, we report on Chandra observations that cover two years after the 2012 outburst of the transient neutron star X-ray binary Swift J174805.3-244637 in the globular cluster Terzan 5. The temperature of the neutron star was elevated during the first two months following its ~8 week accretion episode, but had decayed to the pre-outburst level within ~100 days. Interpreting this as rapid cooling of the accretion-heated crust, we model the observed temperature curve with a thermal evolution code. We find that there is no need to invoke shallow heating for this neutron star, although an extra energy release up to ~1.4 MeV/nucleon is allowed by the current data (2-sigma confidence). We also present two new data points on the crust cooling curve of the 11-Hz X-ray pulsar IGR J17480-2446 in Terzan 5, which was active in 2010. The temperature of this neutron star remained significantly above its pre-outburst level, but we detect no temperature change since the previous measurements of 2013 February. This is consistent with the slower cooling expected several years post-outburst.

The low-background, VUV-sensitive 3-inch diameter photomultiplier tube R11410 has been developed by Hamamatsu for dark matter direct detection experiments using liquid xenon as the target material. We present the results from the joint effort between the XENON collaboration and the Hamamatsu company to produce a highly radio-pure photosensor (version R11410-21) for the XENON1T dark matter experiment. After introducing the photosensor and its components, we show the methods and results of the radioactive contamination measurements of the individual materials employed in the photomultiplier production. We then discuss the adopted strategies to reduce the radioactivity of the various PMT versions. Finally, we detail the results from screening 216 tubes with ultra-low background germanium detectors, as well as their implications for the expected electronic and nuclear recoil background of the XENON1T experiment.

XENON is a dark matter direct detection project, consisting of a time projection chamber (TPC) filled with liquid xenon as detection medium. The construction of the next generation detector, XENON1T, is presently taking place at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy. It aims at a sensitivity to spin-independent cross sections of $2 \cdot 10^{-47} ~ \mathrm{cm}^{\mathrm{2}}$ for WIMP masses around 50 GeV/c$^{2}$, which requires a background reduction by two orders of magnitude compared to XENON100, the current generation detector. An active system that is able to tag muons and muon-induced backgrounds is critical for this goal. A water Cherenkov detector of $\sim$10 m height and diameter has been therefore developed, equipped with 8 inch photomultipliers and cladded by a reflective foil. We present the design and optimization study for this detector, which has been carried out with a series of Monte Carlo simulations. The muon veto will reach very high detection efficiencies for muons ($>99.5%$) and showers of secondary particles from muon interactions in the rock ($>70%$). Similar efficiencies will be obtained for XENONnT, the upgrade of XENON1T, which will later improve the WIMP sensitivity by another order of magnitude. With the Cherenkov water shield studied here, the background from muon-induced neutrons in XENON1T is negligible.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of thermonuclear X-ray bursts on accreting neutron stars. For a summary, we refer to the paper.

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.

If the observed quasi-periodic oscillations in magnetar flares are partially confined to the crust, then the oscillation frequencies are unique probes of the nuclear physics of the neutron star crust. We study crustal oscillations in magnetars including corrections for a finite Alfven velocity. Our crust model uses a new nuclear mass formula that predicts nuclear masses with an accuracy very close to that of the finite range droplet model. This mass model for equilibrium nuclei also includes shell corrections and an updated neutron-drip line. We perturb our crust model to predict axial crust modes and assign them to observed giant flare quasi-periodic oscillation frequencies from the soft gamma-ray repeater SGR 1806-20. We find magnetar crusts that match observations for various magnetic field strengths, entrainment of the free neutron gas in the inner crust, and crust-core transition densities. We find that observations can be reconciled with smaller values of the symmetry energy slope parameter, L, if there is a significant amount of entrainment of the neutrons by the superfluid or if the crust-core transition density is large. We also find neutron star masses and radii which are in agreement with expectations from what is known about low-density matter from nuclear experiment. Matching observations with a field-free model we obtain the approximate values of M =1.35 Msun and R = 11.9 km. Matching observations using a model with the surface dipole field of SGR 1806-20 (B=2.4x10^15 G) we obtain the approximate values of M = 1.25 Msun and R = 12.4 km.

The radioactive isomer $^{83\mathrm{m}}$Kr has many properties that make it very useful for various applications. Its low energy decay products, like conversion, shake-off and Auger electrons as well as X- and $\gamma$-rays are used for calibration purposes in neutrino mass experiments and direct dark matter detection experiments. Thanks to the short half-life of 1.83 h and the decay to the ground state $^{83}$Kr, one does not risk contamination of any low-background experiment with long- lived radionuclides. In this paper, we present two new applications of $^{83\mathrm{m}}$Kr. It can be used as a radioactive tracer in noble gases to characterize the particle flow inside of gas routing systems. A method of doping $^{83\mathrm{m}}$Kr into xenon gas and its detection, using special custom-made detectors, based on a photomultiplier tube, is described. This technique has been used to determine the circulation speed of gas particles inside of a gas purification system for xenon. Furthermore, 83m Kr can be used to rapidly estimate separation performance of a distillation system.

X-ray observations of quiescent X-ray binaries have the potential to provide insight into the structure and the composition of neutron stars. EXO 0748-676 had been actively accreting for over 24 yr before its outburst ceased in late 2008. Subsequent X-ray monitoring revealed a gradual decay of the quiescent thermal emission that can be attributed to cooling of the accretion-heated neutron star crust. In this work, we report on new Chandra and Swift observations that extend the quiescent monitoring to ~5 yr post-outburst. We find that the neutron star temperature remained at ~117 eV between 2009 and 2011, but had decreased to ~110 eV in 2013. This suggests that the crust has not fully cooled yet, which is supported by the lower temperature of ~95 eV that was measured ~4 yr prior to the accretion phase in 1980. Comparing the data to thermal evolution simulations reveals that the apparent lack of cooling between 2009 and 2011 could possibly be a signature of convection driven by phase separation of light and heavy nuclei in the outer layers of the neutron star.

Superbursts from accreting neutron stars probe nuclear reactions at extreme densities ($\rho \approx 10^{9}~g\,cm^{-3}$) and temperatures ($T>10^9~K$). These bursts ($\sim$1000 times more energetic than type I X-ray bursts) are most likely triggered by unstable ignition of carbon in a sea of heavy nuclei made during the rp-process of regular type I X-ray bursts (where the accumulated hydrogen and helium are burned). An open question is the origin of sufficient amounts of carbon, which is largely destroyed during the rp-process in X-ray bursts. We explore carbon production in steady-state burning via the rp-process, which might occur together with unstable burning in systems showing superbursts. We find that for a wide range of accretion rates and accreted helium mass fractions large amounts of carbon are produced, even for systems that accrete solar composition. This makes stable hydrogen and helium burning a viable source of carbon to trigger superbursts. We also investigate the sensitivity of the results to nuclear reactions. We find that the $^{14}$O($\alpha$,p)$^{17}$F reaction rate introduces by far the largest uncertainties in the $^{12}$C yield.

There is compelling evidence that the peak brightness of a Type Ia supernova is affected by the electron fraction Ye at the time of the explosion. The electron fraction is set by the aboriginal composition of the white dwarf and the reactions that occur during the pre explosive convective burning. To date, determining the makeup of the white dwarf progenitor has relied on indirect proxies, such as the average metallicity of the host stellar population. In this paper, we present analytical calculations supporting the idea that the electron fraction of the progenitor systematically influences the nucleosynthesis of silicon group ejecta in Type Ia supernovae. In particular, we suggest the abundances generated in quasi nuclear statistical equilibrium are preserved during the subsequent freezeout. This allows one to potential recovery of Ye at explosion from the abundances recovered from an observed spectra. We show that measurement of 28Si, 32S, 40Ca, and 54Fe abundances can be used to construct Ye in the silicon rich regions of the supernovae. If these four abundances are determined exactly, they are sufficient to recover Ye to 6 percent. This is because these isotopes dominate the composition of silicon-rich material and iron rich material in quasi nuclear statistical equilibrium. Analytical analysis shows that the 28Si abundance is insensitive to Ye, the 32S abundance has a nearly linear trend with Ye, and the 40Ca abundance has a nearly quadratic trend with Ye. We verify these trends with post-processing of 1D models and show that these trends are reflected in model synthetic spectra.

The XENON100 dark matter experiment uses liquid xenon in a time projection chamber (TPC) to measure xenon nuclear recoils resulting from the scattering of dark matter Weakly Interacting Massive Particles (WIMPs). In this paper, we report the observation of single-electron charge signals which are not related to WIMP interactions. These signals, which show the excellent sensitivity of the detector to small charge signals, are explained as being due to the photoionization of impurities in the liquid xenon and of the metal components inside the TPC. They are used as a unique calibration source to characterize the detector. We explain how we can infer crucial parameters for the XENON100 experiment: the secondary-scintillation gain, the extraction yield from the liquid to the gas phase and the electron drift velocity.

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 temperature in the crust of an accreting neutron star, which comprises its outermost kilometer, is set by heating from nuclear reactions at large densities, neutrino cooling, and heat transport from the interior. The heated crust has been thought to affect observable phenomena at shallower depths, such as thermonuclear bursts in the accreted envelope. Here we report that cycles of electron capture and its inverse, $\beta^-$ decay, involving neutron-rich nuclei at a typical depth of about 150 m, cool the outer neutron star crust by emitting neutrinos while also thermally decoupling the surface layers from the deeper crust. This Urca mechanism has been studied in the context of white dwarfs and Type Ia supernovae, but hitherto was not considered in neutron stars, because previous models computed the crust reactions using a zero-temperature approximation and assumed that only a single nuclear species was present at any given depth. This thermal decoupling means that X-ray bursts and other surface phenomena are largely independent of the strength of deep crustal heating. The unexpectedly short recurrence times, of the order of years, observed for very energetic thermonuclear superbursts are therefore not an indicator of a hot crust, but may point instead to an unknown local heating mechanism near the neutron star surface.

We conducted the first long-term (60 days), multiwavelength (optical, ultraviolet, and X-ray) simultaneous monitoring of Cen X-4 with daily Swift observations, with the goal of understanding variability in the low mass X-ray binary Cen X-4 during quiescence. We found Cen X-4 to be highly variable in all energy bands on timescales from days to months, with the strongest quiescent variability a factor of 22 drop in the X-ray count rate in only 4 days. The X-ray, UV and optical (V band) emission are correlated on timescales down to less than 110 s. The shape of the correlation is a power law with index gamma about 0.2-0.6. The X-ray spectrum is well fitted by a hydrogen NS atmosphere (kT=59-80 eV) and a power law (with spectral index Gamma=1.4-2.0), with the spectral shape remaining constant as the flux varies. Both components vary in tandem, with each responsible for about 50% of the total X-ray flux, implying that they are physically linked. We conclude that the X-rays are likely generated by matter accreting down to the NS surface. Moreover, based on the short timescale of the correlation, we also unambiguously demonstrate that the UV emission can not be due to either thermal emission from the stream impact point, or a standard optically thick, geometrically thin disc. The spectral energy distribution shows a small UV emitting region, too hot to arise from the accretion disk, that we identified as a hot spot on the companion star. Therefore, the UV emission is most likely produced by reprocessing from the companion star, indeed the vertical size of the disc is small and can only reprocess a marginal fraction of the X-ray emission. We also found the accretion disc in quiescence to likely be UV faint, with a minimal contribution to the whole UV flux.

The XENON100 experiment, situated in the Laboratori Nazionali del Gran Sasso, aims at the direct detection of dark matter in the form of weakly interacting massive particles (WIMPs), based on their interactions with xenon nuclei in an ultra low background dual-phase time projection chamber. This paper describes the general methods developed for the analysis of the XENON100 data. These methods have been used in the 100.9 and 224.6 live days science runs from which results on spin-independent elastic, spin-dependent elastic and inelastic WIMP-nucleon cross-sections have already been reported.

We conducted the first long-term (60 days), multiwavelength (optical, ultraviolet, and X-ray) simultaneous monitoring of Cen X-4 with daily Swift observations from June to August 2012, with the goal of understanding variability in the low mass X-ray binary Cen X-4 during quiescence. We found Cen X-4 to be highly variable in all energy bands on timescales from days to months, with the strongest quiescent variability a factor of 22 drop in the X-ray count rate in only 4 days. The X-ray, UV and optical (V band) emission are correlated on timescales down to less than 110 s. The shape of the correlation is a power law with index gamma about 0.2-0.6. The X-ray spectrum is well fitted by a hydrogen NS atmosphere (kT=59-80 eV) and a power law (with spectral index Gamma=1.4-2.0), with the spectral shape remaining constant as the flux varies. Both components vary in tandem, with each responsible for about 50% of the total X-ray flux, implying that they are physically linked. We conclude that the X-rays are likely generated by matter accreting down to the NS surface. Moreover, based on the short timescale of the correlation, we also unambiguously demonstrate that the UV emission can not be due to either thermal emission from the stream impact point, or a standard optically thick, geometrically thin disc. The spectral energy distribution shows a small UV emitting region, too hot to arise from the accretion disk, that we identified as a hot spot on the companion star. Therefore, the UV emission is most likely produced by reprocessing from the companion star, indeed the vertical size of the disc is small and can only reprocess a marginal fraction of the X-ray emission. We also found the accretion disc in quiescence to likely be UV faint, with a minimal contribution to the whole UV flux.

The XENON100 experiment, installed underground at the Laboratori Nazionali del Gran Sasso (LNGS), aims to directly detect dark matter in the form of Weakly Interacting Massive Particles (WIMPs) via their elastic scattering off xenon nuclei. This paper presents a study on the nuclear recoil background of the experiment, taking into account neutron backgrounds from ($\alpha$,n) and spontaneous fission reactions due to natural radioactivity in the detector and shield materials, as well as muon-induced neutrons. Based on Monte Carlo simulations and using measured radioactive contaminations of all detector components, we predict the nuclear recoil backgrounds for the WIMP search results published by the XENON100 experiment in 2011 and 2012, 0.11$^{+0.08}_{-0.04}$ events and 0.17$^{+0.12}_{-0.07}$ events, respectively, and conclude that they do not limit the sensitivity of the experiment.

The transient neutron star low-mass X-ray binary and 11 Hz X-ray pulsar IGR J17480-2446 in the globular cluster Terzan 5 exhibited an 11-week accretion outburst in 2010. Chandra observations performed within five months after the end of the outburst revealed evidence that the crust of the neutron star became substantially heated during the accretion episode and was subsequently cooling in quiescence. This provides the rare opportunity to probe the structure and composition of the crust. Here, we report on new Chandra observations of Terzan 5 that extend the monitoring to ~2.2 yr into quiescence. We find that the thermal flux and neutron star temperature have continued to decrease, but remain significantly above the values that were measured before the 2010 accretion phase. This suggests that the crust has not thermally relaxed yet, and may continue to cool. Such behavior is difficult to explain within our current understanding of heating and cooling of transiently accreting neutron stars. Alternatively, the quiescent emission may have settled at a higher observed equilibrium level (for the same interior temperature), in which case the neutron star crust may have fully cooled.

Results from the nuclear recoil calibration of the XENON100 dark matter detector installed underground at the Laboratori Nazionali del Gran Sasso (LNGS), Italy are presented. Data from measurements with an external 241AmBe neutron source are compared with a detailed Monte Carlo simulation which is used to extract the energy dependent charge-yield Qy and relative scintillation efficiency Leff. A very good level of absolute spectral matching is achieved in both observable signal channels - scintillation S1 and ionization S2 - along with agreement in the 2-dimensional particle discrimination space. The results confirm the validity of the derived signal acceptance in earlier reported dark matter searches of the XENON100 experiment.

The quasi-persistent neutron star low-mass X-ray binary MXB 1659-29 went into quiescence in 2001, and we have followed its quiescent X-ray evolution since. Observations over the first 4 years showed a rapid drop in flux and temperature of the neutron star atmosphere, interpreted as cooling of the neutron star crust which had been heated during the 2.5 year outburst. However, observations taken approximately 1400 and 2400 days into quiescence were consistent with each other, suggesting the crust had reached thermal equilibrium with the core. Here we present a new Chandra observation of MXB 1659-29 taken 11 years into quiescence and 4 years since the last Chandra observation. This new observation shows an unexpected factor of ~3 drop in count rate and change in spectral shape since the last observation, which cannot be explained simply by continued cooling. Two possible scenarios are that either the neutron star temperature has remained unchanged and there has been an increase in the column density, or, alternatively the neutron star temperature has dropped precipitously and the spectrum is now dominated by a power-law component. The first scenario may be possible given that MXB 1659-29 is a near edge-on system, and an increase in column density could be due to build-up of material in, and a thickening of, a truncated accretion disk during quiescence. But, a large change in disk height may not be plausible if standard accretion disk theory holds during quiescence. Alternatively, the disk may be precessing, leading to a higher column density during this latest observation.

We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA Star. Improvements in MESA Star's ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA Star. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESA Star calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 Msun stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESA Star solves the fully coupled stellar structure and composition equations, and we show how this has improved MESA's performance scaling on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit (SDK) that packages all the required components needed to form a unified and maintained build environment for MESA. [Abridged]

Quiescent emission from the neutron star low-mass X-ray binary Cen X-4 is seen to be variable on timescales from hundreds of seconds to years, suggesting that at least in this object, low-level accretion is important during quiescence. Here we present results from recent XMM-Newton and Swift observations of Cen X-4, where the X-ray flux (0.5 - 10 keV) varies by a factor of 6.5 between the brightest and faintest states. We find a positive correlation between the X-ray flux and the simultaneous near-UV flux, where as there is no significant correlation between the X-ray and simultaneous optical (V, B) fluxes. This suggests that while the X-ray and UV emitting regions are somehow linked, the optical region originates elsewhere. Comparing the luminosities, it is plausible that the UV emission originates due to reprocessing of the X-ray flux by the accretion disk, with the hot inner region of the disk being a possible location for the UV emitting region. The optical emission, however, could be dominated by the donor star. The X-ray/UV correlation does not favour the accretion stream-impact point as the source of the UV emission.

We report on a search for particle dark matter with the XENON100 experiment, operated at the Laboratori Nazionali del Gran Sasso (LNGS) for 13 months during 2011 and 2012. XENON100 features an ultra-low electromagnetic background of (5.3 \pm 0.6) \times 10^-3 events (kg day keVee)^-1 in the energy region of interest. A blind analysis of 224.6 live days \times 34 kg exposure has yielded no evidence for dark matter interactions. The two candidate events observed in the pre-defined nuclear recoil energy range of 6.6-30.5 keVnr are consistent with the background expectation of (1.0 \pm 0.2) events. A Profile Likelihood analysis using a 6.6-43.3 keVnr energy range sets the most stringent limit on the spin-independent elastic WIMP-nucleon scattering cross section for WIMP masses above 8 GeV/c^2, with a minimum of 2 \times 10^-45 cm^2 at 55 GeV/c^2 and 90% confidence level.

The XENON100 experiment, located at the Laboratori Nazionali del Gran Sasso (LNGS), aims to directly detect dark matter in the form of Weakly Interacting Massive Particles (WIMPs) via their elastic scattering off xenon nuclei. We present a comprehensive study of the predicted electronic recoil background coming from radioactive decays inside the detector and shield materials, and intrinsic contamination. Based on GEANT4 Monte Carlo simulations using a detailed geometry together with the measured radioactivity of all detector components, we predict an electronic recoil background in the WIMP-search energy range (0-100 keV) in the 30 kg fiducial mass of less than 10e-2 events/(kg-day-keV), consistent with the experiment's design goal. The predicted background spectrum is in very good agreement with the data taken during the commissioning of the detector, in Fall 2009.

We present new experimental constraints on the elastic, spin-dependent WIMP-nucleon cross section using recent data from the XENON100 experiment, operated in the Laboratori Nazionali del Gran Sasso in Italy. An analysis of 224.6 live days x 34 kg of exposure acquired during 2011 and 2012 revealed no excess signal due to axial-vector WIMP interactions with 129-Xe and 131-Xe nuclei. This leads to the most stringent upper limits on WIMP-neutron cross sections for WIMP masses above 6 GeV, with a minimum cross section of 3.5 x 10^{-40} cm^2 at a WIMP mass of 45 GeV, at 90% confidence level.

The XENON100 experiment, in operation at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy, was designed to search for evidence of dark matter interactions inside a volume of liquid xenon using a dual-phase time projection chamber. This paper describes the Slow Control System (SCS) of the experiment with emphasis on the distributed architecture as well as on its modular and expandable nature. The system software was designed according to the rules of Object-Oriented Programming and coded in Java, thus promoting code reusability and maximum flexibility during commissioning of the experiment. The SCS has been continuously monitoring the XENON100 detector since mid 2008, remotely recording hundreds of parameters on a few dozen instruments in real time, and setting emergency alarms for the most important variables.

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.

In a recent manuscript (arXiv:1208.5046) Peter Sorensen claims that XENON100's upper limits on spin-independent WIMP-nucleon cross sections for WIMP masses below 10 GeV "may be understated by one order of magnitude or more". Having performed a similar, though more detailed analysis prior to the submission of our new result (arXiv:1207.5988), we do not confirm these findings. We point out the rationale for not considering the described effect in our final analysis and list several potential problems with his study.

We present a study exploring a systematic effect on the brightness of type Ia supernovae using numerical models that assume the single-degenerate paradigm. Our investigation varied the central density of the progenitor white dwarf at flame ignition, and considered its impact on the explosion yield, particularly the production and distribution of radioactive Ni-56, which powers the light curve. We performed a suite of two-dimensional simulations with randomized initial conditions, allowing us to characterize the statistical trends that we present. The simulations indicate that production of Fe-group material is statistically independent of progenitor central density, but the mass of stable Fe-group isotopes is tightly correlated with central density, with a decrease in the production of Ni-56 at higher central densities. These results imply progenitors with higher central densities produce dimmer events. We provide details of the post-explosion distribution of Ni-56 in the models, including the lack of a consistent centrally-located deficit of Ni-56, which may be compared to observed remnants. By performing a self-consistent extrapolation of our model yields and considering the main-sequence lifetime of the progenitor star and the elapsed time between the formation of the white dwarf and the onset of accretion, we develop a brightness-age relation that improves our prediction of the expected trend for single degenerates and we compare this relation with observations.

A new full three-body method is introduced to compute the rate of the triple-alpha capture reaction which is the primary source of $^{12}$C in stars. In this work, we combine the Faddeev hyperspherical harmonics and the R-matrix method to obtain a full solution to the three-body $\alpha+\alpha+\alpha$ continuum. Particular attention is paid to the long range effects caused by the pairwise Coulomb interactions. The new rate agrees with the NACRE rate for temperatures greater than 0.07 GK, but a large enhancement at lower temperature is found ($\approx 10^{14}$ at 0.02 GK). Our results are compared to previous calculations where additional approximations were made. We show that the new rate does not significantly change the evolution of stars around one solar mass. In particular, such stars still undergo a red-giant phase consistent with observations, and no significant differences are found in the final white dwarfs.

The equation of state (EOS) of dense matter has been a long-sought goal of nuclear physics. Equations of state generate unique mass versus radius (M-R) relations for neutron stars, the ultra-dense remnants of stellar evolution. In this work, we determine the neutron star mass-radius relation and, based on recent observations of both transiently accreting and bursting sources, we show that the radius of a 1.4 solar mass neutron star lies between 10.4 and 12.9 km, independent of assumptions about the composition of the core. We show, for the first time, that these constraints remain valid upon removal from our sample of the most extreme transient sources or of the entire set of bursting sources; our constraints also apply even if deconfined quark matter exists in the neutron star core. Our results significantly constrain the dense matter EOS and are, furthermore, consistent with constraints from both heavy-ion collisions and theoretical studies of neutron matter. We predict a relatively weak dependence of the symmetry energy on the density and a value for the neutron skin thickness of lead which is less than 0.20 fm, results that are testable in forthcoming experiments.

Type Ia supernovae are bright stellar explosions distinguished by standardizable light curves that allow for their use as distance indicators for cosmological studies. Despite their highly successful use in this capacity, the progenitors of these events are incompletely understood. We describe simulating type Ia supernovae in the paradigm of a thermonuclear runaway occurring in a massive white dwarf star. We describe the multi-scale physical processes that realistic models must incorporate and the numerical models for these that we employ. In particular, we describe a flame-capturing scheme that addresses the problem of turbulent thermonuclear combustion on unresolved scales. We present the results of our study of the systematics of type Ia supernovae including trends in brightness following from properties of the host galaxy that agree with observations. We also present performance results from simulations on leadership-class architectures.

The XENON100 dark matter experiment uses liquid xenon (LXe) in a time projection chamber (TPC) to search for Xe nuclear recoils resulting from the scattering of dark matter Weakly Interacting Massive Particles (WIMPs). In this paper we present a detailed description of the detector design and present performance results, as established during the commissioning phase and during the first science runs. The active target of XENON100 contains 62 kg of LXe, surrounded by an LXe veto of 99 kg, both instrumented with photomultiplier tubes (PMTs) operating inside the liquid or in Xe gas. The LXe target and veto are contained in a low-radioactivity stainless steel vessel, embedded in a passive radiation shield. The experiment is installed underground at the Laboratori Nazionali del Gran Sasso (LNGS), Italy and has recently published results from a 100 live-days dark matter search. The ultimate design goal of XENON100 is to achieve a spin-independent WIMP-nucleon scattering cross section sensitivity of \sigma = 2x10^-45 cm^2 for a 100 GeV/c^2 WIMP.

The temporal heating and subsequent cooling of the crusts of transiently accreting neutron stars carries unique information about their structure and a variety of nuclear reaction processes. We report on a new Chandra Director's Discretionary Time observation of the globular cluster Terzan 5, aimed to monitor the transiently accreting 11-Hz X-ray pulsar IGR J17480-2446 after the cessation of its recent 10-week long accretion outburst. During the observation, which was performed ~125 days into quiescence, the source displays a thermal spectrum that fits to a neutron star atmosphere model with a temperature for an observer at infinity of kT~92 eV. This is ~10% lower than found ~75 days earlier, yet ~20% higher than the quiescent base level measured prior to the recent outburst. This can be interpreted as cooling of the accretion-heated neutron star crust, and implies that crust cooling is observable after short accretion episodes. Comparison with neutron star thermal evolution simulations indicates that substantial heat must be released at shallow depth inside the neutron star, which is not accounted for in current nuclear heating models.

The location of electron capture heat sources in the crust of accreting neutron stars depends on the masses of extremely neutron-rich nuclei. We present first results from a new implementation of the time-of-flight technique to measure nuclear masses of rare isotopes at the National Superconducting Cyclotron Laboratory. The masses of 16 neutron-rich nuclei in the scandium -- nickel range were determined simultaneously, improving the accuracy compared to previous data in 12 cases. The masses of $^{61}${V}, $^{63}${Cr}, $^{66}${Mn}, and $^{74}${Ni} were measured for the first time with mass excesses of $-30.510(890)$ MeV, $-35.280(650)$ MeV, $-36.900(790)$ MeV, and $-49.210(990)$ MeV, respectively. With the measurement of the $^{66}$Mn mass, the locations of the two dominant electron capture heat sources in the outer crust of accreting neutron stars that exhibit superbursts are now experimentally constrained. We find that the location of the $^{66}$Fe$\rightarrow^{66}$Mn electron capture transition occurs significantly closer to the surface than previously assumed because our new experimental Q-value is 2.1 MeV (2.6$\sigma$) smaller than predicted by the FRDM mass model.

The XENON100 experiment has recently completed a dark matter run with 100.9 live-days of data, taken from January to June 2010. Events in a 48kg fiducial volume in the energy range between 8.4 and 44.6 keVnr have been analyzed. A total of three events have been found in the predefined signal region, compatible with the background prediction of (1.8 \pm 0.6) events. Based on this analysis we present limits on the WIMP-nucleon cross section for inelastic dark matter. With the present data we are able to rule out the explanation for the observed DAMA/LIBRA modulation as being due to inelastic dark matter scattering off iodine at a 90% confidence level.

We present results from the direct search for dark matter with the XENON100 detector, installed underground at the Laboratori Nazionali del Gran Sasso of INFN, Italy. XENON100 is a two-phase time projection chamber with a 62 kg liquid xenon target. Interaction vertex reconstruction in three dimensions with millimeter precision allows to select only the innermost 48 kg as ultra-low background fiducial target. In 100.9 live days of data, acquired between January and June 2010, no evidence for dark matter is found. Three candidate events were observed in a pre-defined signal region with an expected background of 1.8 +/- 0.6 events. This leads to the most stringent limit on dark matter interactions today, excluding spin-independent elastic WIMP-nucleon scattering cross-sections above 7.0x10^-45 cm^2 for a WIMP mass of 50 GeV/c^2 at 90% confidence level.

Many experiments that aim at the direct detection of Dark Matter are able to distinguish a dominant background from the expected feeble signals, based on some measured discrimination parameter. We develop a statistical model for such experiments using the Profile Likelihood ratio as a test statistic in a frequentist approach. We take data from calibrations as control measurements for signal and background, and the method allows the inclusion of data from Monte Carlo simulations. Systematic detector uncertainties, such as uncertainties in the energy scale, as well as astrophysical uncertainties, are included in the model. The statistical model can be used to either set an exclusion limit or to make a discovery claim, and the results are derived with a proper treatment of statistical and systematic uncertainties. We apply the model to the first data release of the XENON100 experiment, which allows to extract additional information from the data, and place stronger limits on the spin-independent elastic WIMP-nucleon scattering cross-section. In particular, we derive a single limit, including all relevant systematic uncertainties, with a minimum of 2.4x10^-44 cm^2 for WIMPs with a mass of 50 GeV/c^2.

We present the results of continued monitoring of the quiescent neutron star low-mass X-ray binary XTE J1701-462 with Chandra and Swift. A new Chandra observation from 2010 October extends our tracking of the neutron star surface temperature from ~800 days to ~1160 days since the end of an exceptionally luminous 19 month outburst. This observation indicates that the neutron star crust may still be slowly cooling toward thermal equilibrium with the core; another observation further into quiescence is needed to verify this. The shape of the overall cooling curve is consistent with that of a broken power law, although an exponential decay to a constant level cannot be excluded with the present data. To investigate possible low-level activity, we conducted a monitoring campaign of XTE J1701-462 with Swift during 2010 April-October. Short-term flares - presumably arising from episodic low-level accretion - were observed up to a luminosity of ~1e35 erg/s, ~20 times higher than the normal quiescent level. We conclude that flares of this magnitude are not likely to have significantly affected the equilibrium temperature of the neutron star and are probably not able to have a measurable impact on the cooling curve. However, it is possible that brighter and longer periods of low-level activity have had an appreciable effect on the equilibrium temperature.

This paper reports the search for quiescent low-mass X-ray binaries (qLMXBs) in the globular cluster (GC) NGC 6553 using an XMM-Newton observation designed specifically for that purpose. We spectrally identify one candidate qLMXB in the core of the cluster, based on the consistency of the spectrum with a neutron star H-atmosphere model at the distance of NGC 6553. Specifically, the best-fit radius found using the three XMM European Photon Imaging Camera spectra is R_NS=6.3(+2.3)(-0.8) km (for M_NS=1.4 Msun) and the best-fit temperature is kTeff=136 (+21)(-34) eV. Both physical parameters are in accordance with typical values of previously identified qLMXBs in GC and in the field, i.e., R_NS~5-20 km and kTeff~50-150 eV. A power-law (PL) component with a photon index Gamma=2.1(+0.5)(-0.8) is also required for the spectral fit and contributes to ~33% of the total flux of the X-ray source. A detailed analysis supports the hypothesis that the PL component originates from nearby sources in the core, unresolved with XMM. The analysis of an archived Chandra observation provides marginal additional support to the stated hypothesis. Finally, a catalog of all the sources detected within the XMM field of view is presented here.

Dark Matter and Double Beta Decay experiments require extremely low radioactivity within the detector materials. For this purpose, the University of California, Los Angeles and Hamamatsu Photonics have developed the QUartz Photon Intensifying Detector (QUPID), an ultra-low background photodetector based on the Hybrid Avalanche Photo Diode (HAPD) and entirely made of ultraclean synthetic fused silica. In this work we present the basic concept of the QUPID and the testing measurements on QUPIDs from the first production line. Screening of radioactivity at the Gator facility in the Laboratori Nazionali del Gran Sasso has shown that the QUPIDs safely fulfill the low radioactive contamination requirements for the next generation zero background experiments set by Monte Carlo simulations. The quantum efficiency of the QUPID at room temperature is > 30% at the xenon scintillation wavelength. At low temperatures, the QUPID shows a leakage current less than 1 nA and a global gain of 10^5. In these conditions, the photocathode and the anode show > 95% linearity up to 1 uA for the cathode and 3 mA for the anode. The photocathode and collection efficiency are uniform to 80% over the entire surface. In parallel with single photon counting capabilities, the QUPIDs have a good timing response: 1.8 +/- 0.1 ns rise time, 2.5 +/- 0.2 ns fall time, 4.20 +/- 0.05 ns pulse width, and 160 +/- 30 ps transit time spread. The QUPIDs have also been tested in a liquid xenon environment, and scintillation light from 57Co and 210Po radioactive sources were observed.

This paper reports the spectral and timing analyses of the quiescent low-mass X-ray binary U24 observed during five archived Chandra-ACIS exposures of the nearby globular cluster NGC 6397, for a total of 350 ksec. We find that the X-ray flux and the parameters of the hydrogen atmosphere spectral model are consistent with those previously published for this source. On short timescales, we find no evidence of aperiodic intensity variability, with 90% confidence upper limits during five observations ranging between <8.6% rms and <19% rms, in the 0.0001-0.1 Hz frequency range (0.5-8.0 keV); and no evidence of periodic variability, with maximum observed powers in this frequency range having a chance probability of occurrence from a Poisson-deviated light curve in excess of 10%. We also report the improved neutron star physical radius measurements, with statistical accuracy of the order of ~10%: R_ns = 8.9(+0.9)(-0.6) km for M_ns = 1.4 Msun. Alternatively, we provide the confidence regions in mass-radius space as well as the best-fit projected radius R_infinity= 11.9(+1.0)(-0.8)km, as seen by an observer at infinity. The best-fit effective temperature, kTeff = 80(+4)(-5) eV, is used to estimate the neutron star core temperature which falls in the range T_core = (3.0 - 9.8) x10 7 K, depending on the atmosphere model considered. This makes U24 the third most precisely measured neutron star radius among qLMXBs, after those in OmCen and in M13.

The XENON100 experiment, in operation at the Laboratori Nazionali del Gran Sasso in Italy, is designed to search for dark matter WIMPs scattering off 62 kg of liquid xenon in an ultra-low background dual-phase time projection chamber. In this letter, we present first dark matter results from the analysis of 11.17 live days of non-blind data, acquired in October and November 2009. In the selected fiducial target of 40 kg, and within the pre-defined signal region, we observe no events and hence exclude spin-independent WIMP-nucleon elastic scattering cross-sections above 3.4 x 10^-44 cm^2 for 55 GeV/c^2 WIMPs at 90% confidence level. Below 20 GeV/c^2, this result constrains the interpretation of the CoGeNT and DAMA signals as being due to spin-independent, elastic, light mass WIMP interactions.

Type Ia supernovae are bright stellar explosions thought to occur when a thermonuclear runaway consumes roughly a solar mass of degenerate stellar material. These events produce and disseminate iron-peak elements, and properties of their light curves allow for standardization and subsequent use as cosmological distance indicators. The explosion mechanism of these events remains, however, only partially understood. Many models posit the explosion beginning with a deflagration born near the center of a white dwarf that has gained mass from a stellar companion. In order to match observations, models of this single-degenerate scenario typically invoke a subsequent transition of the (subsonic) deflagration to a (supersonic) detonation that rapidly consumes the star. We present an investigation into the systematics of thermonuclear supernovae assuming this paradigm. We utilize a statistical framework for a controlled study of two-dimensional simulations of these events from randomized initial conditions. We investigate the effect of the composition and thermal history of the progenitor on the radioactive yield, and thus brightness, of an event. Our results offer an explanation for some observed trends of mean brightness with properties of the host galaxy.

In late 2008, the quasi-persistent neutron star X-ray transient and eclipsing binary EXO 0748-676 started a transition from outburst to quiescence, after it had been actively accreting for more than 24 years. In a previous work, we discussed Chandra and Swift observations obtained during the first five months after this transition. Here, we report on further X-ray observations of EXO 0748-676, extending the quiescent monitoring to 1.6 years. Chandra and XMM-Newton data reveal quiescent X-ray spectra composed of a soft, thermal component that is well-fitted by a neutron star atmosphere model. An additional hard powerlaw tail is detected that changes non-monotonically over time, contributing between 4 and 20 percent to the total unabsorbed 0.5-10 keV flux. The combined set of Chandra, XMM-Newton and Swift data reveals that the thermal bolometric luminosity fades from ~1E34 to 6E33 (D/7.4 kpc)^2 erg/s, whereas the inferred neutron star effective temperature decreases from ~124 to 109 eV. We interpret the observed decay as cooling of the neutron star crust and show that the fractional quiescent temperature change of EXO 0748-676 is markedly smaller than observed for three other neutron star X-ray binaries that underwent prolonged accretion outbursts.

Some neutron star low-mass X-ray binaries have very long outbursts (lasting several years) which can generate a significant amount of heat in the neutron star crust. After the system has returned to quiescence, the crust then thermally relaxes. This provides a rare opportunity to study the thermal properties of neutron star crusts, putting constraints on the thermal conductivity and hence the structure and composition of the crust. KS 1731-260 is one of only four systems where this crustal cooling has been observed. Here, we present a new Chandra observation of this source approximately 8 years after the end of the last outburst, and 4 years since the last observation. We find that the source has continued to cool, with the cooling curve displaying a simple power-law decay. This suggests that the crust has not fully thermally relaxed yet, and may continue to cool further. A simple power law decay is in contrast to theoretical cooling models of the crust, which predict that the crust should now have cooled to the same temperature as the neutron star core.

We determine an empirical dense matter equation of state from a heterogeneous dataset of six neutron stars: three type I X-ray bursters with photospheric radius expansion, studied by Ozel et al., and three transient low-mass X-ray binaries. We critically assess the mass and radius determinations from the X-ray burst sources and show explicitly how systematic uncertainties, such as the photospheric radius at touchdown, affect the most probable masses and radii. We introduce a parameterized equation of state and use a Markov Chain Monte Carlo algorithm within a Bayesian framework to determine nuclear parameters such as the incompressibility and the density dependence of the bulk symmetry energy. Using this framework we show, for the first time, that these parameters, predicted solely on the basis of astrophysical observations, all lie in ranges expected from nuclear systematics and laboratory experiments. We find significant constraints on the mass-radius relation for neutron stars, and hence on the pressure-density relation of dense matter. The predicted symmetry energy and the equation of state near the saturation density are soft, resulting in relatively small neutron star radii around 11-12 km for M=1.4 Msun. The predicted equation of state stiffens at higher densities, however, and our preferred model for X-ray bursts suggests that the neutron star maximum mass is relatively large, 1.9-2.2 Msun. Our results imply that several commonly used equations of state are inconsistent with observations.

The nearby neutron star low-mass X-ray binary, Cen X-4, has been in a quiescent state since its last outburst in 1979. Typically, quiescent emission from these objects consists of thermal emission (presumably from the neutron star surface) with an additional hard power-law tail of unknown nature. Variability has been observed during quiescence in Cen X-4 on both timescales as short as hundreds of seconds and as long as years. However, the nature of this variability is still unknown. Early observations seemed to show it was all due to a variable hard X-ray tail. Here, we present new and archival observations that contradict this. The most recent Suzaku observation of Cen X-4 finds it in a historically low state, a factor of 4.4 fainter than the brightest quiescent observation. As the spectrum during the brightest observation was comprised of approximately 60% from the thermal component and 40% from the power-law component, such a large change cannot be explained by just power-law variability. Spectral fits with a variable thermal component fit the data well, while spectral fits allowing both the column density and the power-law to vary do not, leading to the conclusion that the thermal component must be variable. Interestingly, we also find that the thermal fraction remains consistent between all epochs, implying that the thermal and power-law fluxes vary by approximately the same amount. If the emitting area remains unchanged between observations, then the effective surface temperature must change. Alternatively, if the temperature remains constant, then the emitting area must change. The nature of this thermal variability is unclear, but may be explained by variable low-level accretion.

We explore the effects of the deflagration to detonation transition (DDT) density on the production of Ni-56 in thermonuclear supernova explosions (type Ia supernovae). Within the DDT paradigm, the transition density sets the amount of expansion during the deflagration phase of the explosion and therefore the amount of nuclear statistical equilibrium (NSE) material produced. We employ a theoretical framework for a well-controlled statistical study of two-dimensional simulations of thermonuclear supernovae with randomized initial conditions that can, with a particular choice of transition density, produce a similar average and range of Ni-56 masses to those inferred from observations. Within this framework, we utilize a more realistic "simmered" white dwarf progenitor model with a flame model and energetics scheme to calculate the amount of Ni-56 and NSE material synthesized for a suite of simulated explosions in which the transition density is varied in the range 1-3x10^7 g/cc. We find a quadratic dependence of the NSE yield on the log of the transition density, which is determined by the competition between plume rise and stellar expansion. By considering the effect of metallicity on the transition density, we find the NSE yield decreases by 0.055 +/- 0.004 solar masses for a 1 solar metallicity increase evaluated about solar metallicity. For the same change in metallicity, this result translates to a 0.067 +/- 0.004 solar mass decrease in the Ni-56 yield, slightly stronger than that due to the variation in electron fraction from the initial composition. Observations testing the dependence of the yield on metallicity remain somewhat ambiguous, but the dependence we find is comparable to that inferred from some studies.

Recent observational studies of type Ia supernovae (SNeIa) suggest correlations between the peak brightness of an event and the age of the progenitor stellar population. This trend likely follows from properties of the progenitor white dwarf (WD), such as central density, that follow from properties of the host stellar population. We present a statistically well-controlled, systematic study utilizing a suite of multi-dimensional SNeIa simulations investigating the influence of central density of the progenitor WD on the production of Fe-group material, particularly radioactive Ni-56, which powers the light curve. We find that on average, as the progenitor's central density increases, production of Fe-group material does not change but production of Ni-56 decreases. We attribute this result to a higher rate of neutronization at higher density. The central density of the progenitor is determined by the mass of the WD and the cooling time prior to the onset of mass transfer from the companion, as well as the subsequent accretion heating and neutrino losses. The dependence of this density on cooling time, combined with the result of our central density study, offers an explanation for the observed age-luminosity correlation: a longer cooling time raises the central density at ignition thereby producing less Ni-56 and thus a dimmer event. While our ensemble of results demonstrates a significant trend, we find considerable variation between realizations, indicating the necessity for averaging over an ensemble of simulations to demonstrate a statistically significant result.

We study the stability properties of the Cauchy horizon for two different self-dual black hole solutions obtained in a model inspired by Loop Quantum Gravity. The self-dual spacetimes depend on a free dimensionless parameter called a polymeric parameter P. For the first metric the Cauchy horizon is stable for supermassive black holes only if this parameter is sufficiently small. For small black holes, however the stability is easily implemented. The second metric analyzed is not only self-dual but also "form-invariant" under the transformation r -> r*^2/r and r* = 2 m P. We find that this symmetry protects the Cauchy horizon for any value of the polymeric parameter.

We present Rossi X-Ray Timing Explorer and Swift observations made during the final three weeks of the 2006-2007 outburst of the super-Eddington neutron star (NS) transient XTE J1701-462, as well as Chandra and XMM-Newton observations covering the first ~800 days of the subsequent quiescent phase. The source transitioned quickly from active accretion to quiescence, with the luminosity dropping by over 3 orders of magnitude in ~13 days. The spectra obtained during quiescence exhibit both a thermal component, presumed to originate in emission from the NS surface, and a non-thermal component of uncertain origin, which has shown large and irregular variability. We interpret the observed decay of the inferred effective surface temperature of the NS in quiescence as the cooling of the NS crust after having been heated and brought out of thermal equilibrium with the core during the outburst. The interpretation of the data is complicated by an apparent temporary increase in temperature ~220 days into quiescence, possibly due to an additional spurt of accretion. We derive an exponential decay timescale of ~120 (+30/-20) days for the inferred temperature (excluding observations affected by the temporary increase). This short timescale indicates a highly conductive NS crust. Further observations are needed to confirm whether the crust is still slowly cooling or has already reached thermal equilibrium with the core at a surface temperature of ~125 eV. The latter would imply a high equilibrium bolometric thermal luminosity of ~5x10^{33} erg/s for an assumed distance of 8.8 kpc.

Observationally inferred superburst ignition depths are shallower than models predict. We address this discrepancy by reexamining the superburst trigger mechanism. We first explore the hypothesis of Kuulkers et al. that exothermic electron captures trigger superbursts. We find that all electron capture reactions are thermally stable in accreting neutron star oceans and thus are not a viable trigger mechanism. Fusion reactions other than 12C + 12C are infeasible as well since the possible reactants either deplete at much shallower depths or have prohibitively large Coulomb barriers. Thus we confirm the proposal of Cumming & Bildsten and Strohmayer & Brown that 12C + 12C triggers superbursts. We then examine the 12C + 12C fusion rate. The reaction cross-section is experimentally unknown at astrophysically relevant energies, but resonances exist in the 12C + 12C system throughout the entire measured energy range. Thus it is likely, and in fact has been predicted, that a resonance exists near the Gamow peak energy ~ 1.5 MeV. For such a hypothetical 1.5 MeV resonance, we derive both a fiducial value and upper limit to the resonance strength and find that such a resonance could decrease the theoretically predicted superburst ignition depth by up to a factor of 4; in this case, observationally inferred superburst ignition depths would accord with model predictions for a range of plausible neutron star parameters. Said differently, such a resonance would decrease the temperature required for unstable 12C ignition at a column depth 10^12 g/cm^2 from 6 x 10^8 K to 5 x 10^8 K. Determining the existence of a strong resonance in the Gamow window requires measurements of the 12C + 12C cross-section down to a center-of-mass energy near 1.5 MeV, which is within reach of the proposed DUSEL facility.

We present a theoretical framework for formal study of systematic effects in Supernovae Type Ia (SN Ia) that utilizes 2-d simulations to implement a form of the deflagration-detonation transition (DDT) explosion scenario. The framework is developed from a randomized initial condition that leads to a sample of simulated SN Ia whose Ni56 masses have a similar average and range to those observed, and have many other modestly realistic features such as the velocity extent of intermediate mass elements. The intended purpose is to enable statistically well-defined studies of both physical and theoretical parameters of the SN Ia explosion simulation. We present here a thorough description of the outcome of the SN Ia explosions produced by our current simulations. A first application of this framework is utilized to study the dependence of the SN Ia on the Ne22 content, which is known to be directly influenced by the progenitor stellar population's metallicity. Our study is very specifically tailored to measure how the Ne22 content influences the competition between the rise of plumes of burned material and the expansion of the star before these plumes reach DDT conditions. This competition controls the amount of material in nuclear statistical equilibrium (NSE) and therefore Ni56 produced by setting the density at which nucleosynthesis takes place during the detonation phase of the explosion. Although the outcome following from any particular ignition condition can change dramatically with Ne22 content, with a sample of 20 ignition conditions we find that the systematic change in the expansion of the star prior to detonation is not large enough to compete with the dependence on initial neutron excess discussed by Timmes, Brown & Truran (2003). (Abridged)

This paper presents the analysis of candidate quiescent low mass xray binarie (qLMXBs) observed during a short Chandra/ACIS observation of the globular cluster (GC) NGC 6304. Two out of the three candidate qLMXBs of this cluster, XMMU 171433-292747 and XMMU 171421-292917, lie within the field of view. This permits comparison with the discovery observation of these sources. The one in the GC core -- XMMU 171433-292747 -- is spatially resolved into two separate X-ray sources, one of which is consistent with a pure H-atmosphere qLMXB, and the other is an X-ray power-law spectrum source. These two spectral components separately account for those observed from XMMU 171433-292747 in its discovery observation. We find that the observed flux and spectral parameters of the H-atmosphere spectral components are consistent with the previous observation, as expected from a qLMXB powered by deep crustal heating. XMMU 171421-292917 also has neutron star atmosphere spectral parameters consistent with those in the XMM-Newton observation and the observed flux has decreased by a factor 0.54^{+0.30}_{-0.24}.

The monitoring of quiescent emission from neutron star transients with accretion outbursts long enough to significantly heat the neutron star crust has opened a new vista onto the physics of dense matter. In this paper we construct models of the thermal relaxation of the neutron star crust following the end of a protracted accretion outburst. We confirm the finding of Shternin et al., that the thermal conductivity of the neutron star crust is high, consistent with a low impurity parameter. We describe the basic physics that sets the broken power-law form of the cooling lightcurve. The initial power law decay gives a direct measure of the temperature profile, and hence the thermal flux during outburst, in the outer crust. The time of the break, at hundreds of days post-outburst, corresponds to the thermal time where the solid transitions from a classical to quantum crystal, close to neutron drip. We calculate in detail the constraints on the crust parameters of both KS 1731-260 and MXB 1659-29 from fitting their cooling lightcurves. Our fits to the lightcurves require that the neutrons do not contribute significantly to the heat capacity in the inner crust, and provide evidence in favor of the existence of a neutron superfluid throughout the inner crust. Our fits to both sources indicate an impurity parameter of order unity in the inner crust.

The quasi-persistent neutron star X-ray transient and eclipsing binary EXO 0748-676 recently started the transition to quiescence following an accretion outburst that lasted more than 24 years. We report on two Chandra and twelve Swift observations performed within five months after the end of the outburst. The Chandra spectrum is composed of a soft, thermal component that fits to a neutron star atmosphere model with kT^inf~0.12 keV, joined by a hard powerlaw tail that contributes ~20% of the total 0.5-10 keV unabsorbed flux. The combined Chandra/Swift data set reveals a relatively hot and luminous quiescent system with a temperature of kT^inf~0.11-0.13 keV and a bolometric thermal luminosity of ~8.1E33-1.6E34 (d/7.4 kpc)^2 erg/s. We discuss our results in the context of cooling neutron star models.

In the next decade Type Ia supernovae (SNe Ia) will be used to test theories predicting changes in the Dark Energy equation of state with time. Ultimately this requires a dedicated space mission like JDEM. SNe Ia are mature cosmological probes --- their limitations are well characterized, and a path to improvement is clear. Dominant systematic errors include photometric calibration, selection effects, reddening, and population-dependent differences. Building on past lessons, well-controlled new surveys are poised to make strides in these areas: the Palomar Transient Factory, Skymapper, La Silla QUEST, Pan-STARRS, the Dark Energy Survey, LSST, and JDEM. They will obviate historical calibrations and selection biases, and allow comparisons via large subsamples. Some systematics follow from our ignorance of SN Ia progenitors, which there is hope of determining with SN Ia rate studies from 0<z<4. Aside from cosmology, SNe Ia regulate galactic and cluster chemical evolution, inform stellar evolution, and are laboratories for extreme physics. Essential probes of SNe Ia in these contexts include spectroscopy from the UV to the IR, X-ray cluster and SN remnant observations, spectropolarimetry, and advanced theoretical studies. While there are an abundance of discovery facilities planned, there is a deficit of follow-up resources. Living in the systematics era demands deep understanding rather than larger statistics. NOAO ReSTAR initiative to build 2-4m telescopes would provide necessary follow-up capability. Finally, to fully exploit LSST, well-matched wide-field spectroscopic capabilities are desirable.

We highlight recent theoretical and observational progress in several areas of neutron star astrophysics, and discuss the prospect for advances in the next decade.

This white paper, directed to the Stars and Stellar Evolution panel, has three objectives: 1) to provide the Astro2010 Decadal Survey with a vista into the goals of the nuclear physics and nuclear astrophysics community; 2) to alert the astronomical community of joint opportunities for discoveries at the interface between nuclear physics and astronomy; and 3) to delineate efforts in nuclear physics and describe the observational and theoretical advances in astrophysics necessary to make progress towards answering the following questions in the Nuclear Science 2007 Long Range Plan: 1) What is the origin and distribution of the elements? 2) What are the nuclear reactions that power stars and stellar explosions? 3) What is the nature of dense matter? The scope of this white paper concerns the specific area of "low energy" nuclear astrophysics. We define this as the area of overlap between astrophysics and the study of nuclear structure and reactions. Of the questions listed above, two -- What is the origin of the elements? and What is the nature of dense matter? -- were specifically listed in the National Academies Study, "Connecting Quarks with the Cosmos".

A multi-target detection system XAX, comprising concentric 10 ton targets of 136Xe and 129/131Xe, together with a geometrically similar or larger target of liquid Ar, is described. Each is configured as a two-phase scintillation/ionization TPC detector, enhanced by a full 4pi array of ultra-low radioactivity Quartz Photon Intensifying Detectors (QUPIDs) replacing the conventional photomultipliers for detection of scintillation light. It is shown that background levels in XAX can be reduced to the level required for dark matter particle (WIMP) mass measurement at a 10^-10 pb WIMP-nucleon cross section, with single-event sensitivity below 10^-11 pb. The use of multiple target elements allows for confirmation of the A^2 dependence of a coherent cross section, and the different Xe isotopes provide information on the spin-dependence of the dark matter interaction. The event rates observed by Xe and Ar would modulate annually with opposite phases from each other for WIMP mass >~100 GeV/c^2. The large target mass of 136Xe and high degree of background reduction allow neutrinoless double beta decay to be observed with lifetimes of 10^27-10^28 years, corresponding to the Majorana neutrino mass range 0.01-0.1 eV, the most likely range from observed neutrino mass differences. The use of a 136Xe-depleted 129/131Xe target will also allow measurement of the pp solar neutrino spectrum to a precision of 1-2%.

We report the search for low-mass X-ray binaries in quiescence (qLMXBs) in the globular cluster NGC 6304 using XMM observations. We present the spectral analysis leading to the identification of three candidate qLMXBs within the field of this globular cluster (GC), each consistent with the X-ray spectral properties of previously identified qLMXBs in the field and in other globular clusters -- specifically, with a hydrogen atmosphere neutron star with radius between 5--20\km. One (source 4, with R=11.7^{+8.3}_{-0.4} (D/5.97 kpc) km and kT_eff=117^{+59}_{-44} eV) is located within one core radius (r_c) of the centre of NGC 6304. This candidate also presents a spectral power-law component contributing 49 per cent of the 0.5-10 keV flux. A second one (source 9 with R=15.3^{+11.2}_{-6.5} (D/5.97 kpc) km and kT_eff=100^{+24}_{-19} eV) is found well outside the optical core (at 32 r_c) but still within the tidal radius. From spatial coincidence, we identify a bright 2MASS infrared counterpart which, at the distance of NGC 6304, seems to be a post-asymptotic giant branch star. The third qLMXB (source 5 with R=23^{+38}_{-14} (D/5.97 kpc) km and kT_eff=70^{+28}_{-20} eV) is a low signal-to-noise candidate for which we also identify from spatial coincidence a bright 2MASS infrared counterpart, with 99.916 per cent confidence. Three qLMXBs from this GC is marginally consistent with that expected from the encounter rate of NGC 6304. We also report a low signal-to-noise source with an unusually hard photon index (\alpha=-2.0^{+1.2}_{-2.2}). Finally, we present an updated catalogue of the X-ray sources lying in the field of NGC 6304, and compare this with the previous catalogue compiled from ROSAT observations.

Timmes, Brown & Truran found that metallicity variations could theoretically account for a 25% variation in the mass of 56Ni synthesized in Type Ia supernovae (SNe Ia), and thus account for a large fraction of the scatter in observed SN Ia luminosities. Higher-metallicity progenitors are more neutron-rich, producing more stable burning products relative to radioactive 56Ni. We develop a new method for estimating bolometric luminosity and 56Ni yield in SNe Ia and use it to test the theory with data from the Supernova Legacy Survey. We find that the average 56Ni yield does drop in SNe Ia from high metallicity environments, but the theory can only account for 7%--10% of the dispersion in SN Ia 56Ni mass, and thus luminosity. This is because the effect is dominant at metallicities significantly above solar, whereas we find that SN hosts have predominantly subsolar or only moderately above-solar metallicities. We also show that allowing for changes in O/Fe with the metallicity [Fe/H] does not have a major effect on the theoretical prediction of Timmes, Brown & Truran, so long as one is using the O/H as the independent variable. Age may have a greater effect than metallicity -- we find that the luminosity weighted age of the host galaxy is correlated with 56Ni yield, and thus more massive progenitors give rise to more luminous explosions. This is hard to understand if most SNe Ia explode when the primaries reach the Chandrasekhar mass. Finally, we test the findings of Gallagher et al., that the residuals of SNe Ia from the Hubble diagram are correlated with host galaxy metallicity, and we find no such correlation.

In quasi-persistent neutron star transients, long outbursts cause the neutron star crust to be heated out of thermal equilibrium with the rest of the star. During quiescence, the crust then cools back down. Such crustal cooling has been observed in two quasi-persistent sources: KS 1731-260 and MXB 1659-29. Here we present an additional Chandra observation of MXB 1659-29 in quiescence, which extends the baseline of monitoring to 6.6 yr after the end of the outburst. This new observation strongly suggests that the crust has thermally relaxed, with the temperature remaining consistent over 1000 days. Fitting the temperature cooling curve with an exponential plus constant model we determine an e-folding timescale of 465 +/- 25 days, with the crust cooling to a constant surface temperature of kT = 54 +/- 2 eV (assuming D=10 kpc). From this, we infer a core temperature in the range 3.5E7-8.3E7 K (assuming D=10 kpc), with the uncertainty due to the surface composition. Importantly, we tested two neutron star atmosphere models as well as a blackbody model, and found that the thermal relaxation time of the crust is independent of the chosen model and the assumed distance.

Proton-rich material in a state of nuclear statistical equilibrium (NSE) is one of the least studied regimes of nucleosynthesis. One reason for this is that after hydrogen burning, stellar evolution proceeds at conditions of equal number of neutrons and protons or at a slight degree of neutron-richness. Proton-rich nucleosynthesis in stars tends to occur only when hydrogen-rich material that accretes onto a white dwarf or neutron star explodes, or when neutrino interactions in the winds from a nascent proto-neutron star or collapsar-disk drive the matter proton-rich prior to or during the nucleosynthesis. In this paper we solve the NSE equations for a range of proton-rich thermodynamic conditions. We show that cold proton-rich NSE is qualitatively different from neutron-rich NSE. Instead of being dominated by the Fe-peak nuclei with the largest binding energy per nucleon that have a proton to nucleon ratio close to the prescribed electron fraction, NSE for proton-rich material near freeze-out temperature is mainly composed of Ni56 and free protons. Previous results of nuclear reaction network calculations rely on this non-intuitive high proton abundance, which this paper will explain. We show how the differences and especially the large fraction of free protons arises from the minimization of the free energy as a result of a delicate competition between the entropy and the nuclear binding energy.

Prior to the explosion of a carbon-oxygen white dwarf in a Type Ia supernova there is a long "simmering," during which the 12C + 12C reaction gradually heats the white dwarf on a long (~ 1000 yr) timescale. Piro & Bildsten showed that weak reactions during this simmering set a maximum electron abundance Ye at the time of the explosion. We investigate the nuclear reactions during this simmering with a series of self-heating, at constant pressure, reaction network calculations. Unlike in AGB stars, proton captures onto 22Ne and heavier trace nuclei do not play a significant role. The 12C abundance is sufficiently high that the neutrons preferentially capture onto 12C, rather than iron group nuclei. As an aid to hydrodynamical simulations of the simmering phase, we present fits to the rates of heating, electron capture, change in mean atomic mass, and consumption of 12C in terms of the screened thermally averaged cross section for 12C + 12C. Our evaluation of the net heating rate includes contributions from electron captures into the 3.68 MeV excited state of 13C. This results in a slightly larger energy release, per 12C consumed, than that found by Piro & Bildsten, but less than that released for a burn to only 20Ne and 23Na. We compare our one-zone results to more accurate integrations over the white dwarf structure to estimate the amount of 12C that must be consumed to raise the white dwarf temperature, and hence to determine the net reduction of Ye during simmering.

We report the analysis of the first superburst from a transiently accreting neutron star system with the All-Sky Monitor (ASM) on the Rossi X-ray Timing Explorer. The superburst occurred 55 days after the onset of an accretion outburst in 4U 1608-522. During that time interval, the accretion rate was at least 7% of the Eddington limit. The peak flux of the superburst is 22 to 45% of the Eddington limit, and its radiation energy output is between 4e41 and 9e41 erg for a distance of 3.2 kpc. Fits of cooling models to the superburst light curve indicate an ignition column depth between 1.5e12 and 4.1e12 g/cm2. Extrapolating the accretion history observed by the ASM, we derive that this column was accreted over a period of 26 to 72 years. The superburst characteristics are consistent with those seen in other superbursting low-mass X-ray binaries. However, the transient nature of the hosting binary presents significant challenges for superburst theory, requiring additional ingredients for the models. The carbon that fuels the superburst is thought to be produced mostly during the accretion outbursts and destroyed in the frequent type-I X-ray bursts. Mixing and sedimentation of the elements in the neutron star envelope may significantly influence the balance between the creation and destruction of carbon. Furthermore, predictions for the temperature of the neutron star crust fail to reach the values required for the ignition of carbon at the inferred column depth.

We perform a self-consistent calculation of the thermal structure in the crust of a superbursting neutron star. In particular, we follow the nucleosynthetic evolution of an accreted fluid element from its deposition into the atmosphere down to a depth where the electron Fermi energy is 20 MeV. We include temperature-dependent continuum electron capture rates and realistic sources of heat loss by thermal neutrino emission from the crust and core. We show that, in contrast to previous calculations, electron captures to excited states and subsequent gamma-emission significantly reduce the local heat loss due to weak-interaction neutrinos. Depending on the initial composition these reactions release up to a factor of 10 times more heat at densities < 10^{11} g/cc than obtained previously. This heating reduces the ignition depth of superbursts. In particular, it reduces the discrepancy noted by Cumming et al. between the temperatures needed for unstable 12C ignition on timescales consistent with observations and the reduction in crust temperature from Cooper pair neutrino emission.

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.

Carbon-oxygen white dwarfs contain neon-22 formed from alpha-captures onto nitrogen during core He burning in the progenitor star. In a white dwarf (type Ia) supernova, the neon-22 abundance determines, in part, the neutron-to-proton ratio and hence the abundance of radioactive nickel-56 that powers the lightcurve. The neon-22 abundance also changes the burning rate and hence the laminar flame speed. We tabulate the flame speedup for different initial carbon and neon-22 abundances and for a range of densities. This increase in the laminar flame speed--about 30% for a neon-22 mass fraction of 6%--affects the deflagration just after ignition near the center of the white dwarf, where the laminar speed of the flame dominates over the buoyant rise, and in regions of lower density ~ 10^7 g/cm3 where a transition to distributed burning is conjectured to occur. The increase in flame speed will decrease the density of any transition to distributed burning.

Neutron stars, with their strong surface gravity, have interestingly short timescales for the sedimentation of heavy elements. Recent observations of unstable thermonuclear burning (observed as X-ray bursts) on the surfaces of slowly accreting neutron stars ($< 0.01$ of the Eddington rate) motivate us to examine how sedimentation of CNO isotopes affects the ignition of these bursts. We further estimate the burst development using a simple one-zone model with a full reaction network. We report a region of mass accretion rates for weak H flashes. Such flashes can lead to a large reservoir of He, the unstable burning of which may explain some observed long bursts (duration $\sim 1000$ s).

We develop and calibrate a realistic model flame for hydrodynamical simulations of deflagrations in white dwarf (Type Ia) supernovae. Our flame model builds on the advection-diffusion-reaction model of Khokhlov and includes electron screening and Coulomb corrections to the equation of state in a self-consistent way. We calibrate this model flame--its energetics and timescales for energy release and neutronization--with self-heating reaction network calculations that include both these Coulomb effects and up-to-date weak interactions. The burned material evolves post-flame due to both weak interactions and hydrodynamic changes in density and temperature. We develop a scheme to follow the evolution, including neutronization, of the NSE state subsequent to the passage of the flame front. As a result, our model flame is suitable for deflagration simulations over a wide range of initial central densities and can track the temperature and electron fraction of the burned material through the explosion and into the expansion of the ejecta.

Neutron stars, with their strong surface gravity, have interestingly short timescales for the sedimentation of heavy elements. Motivated by observations of Type I X-ray bursts from sources with extremely low persistent accretion luminosities, $L_X < 10^{36}\usp\ergspersecond (\simeq 0.01\ensuremath{L_{\mathrm{Edd}}}$), we study how sedimentation affects the distribution of isotopes and the ignition of H and He in the envelope of an accreting neutron star. For local mass accretion rates $\mdot \lesssim 10^{-2}\medd$ (for which the ignition of H is unstable), where $\medd = 8.8\times 10^{4}\nsp\gpscps$, the helium and CNO elements sediment out of the accreted fuel before reaching a temperature where H would ignite. Using one-zone calculations of the thermonuclear burning, we find a range of accretion rates for which the unstable H ignition does not trigger unstable He burning. This range depends on the emergent flux from reactions in the deep neutron star crust; for $F = 0.1\nsp\MeV(\dot{m}/\mb)$, the range is $3\times 10^{-3}\medd\lesssim\mdot\lesssim 10^{-2}\medd$. We speculate that sources accreting in this range will build up a massive He layer that later produces an energetic and long X-ray burst. At mass accretion rates lower than this range, we find that the H flash leads to a strong mixed H/He flash. Surprisingly, even at accretion rates $\mdot \gtrsim 0.1\medd$, although the H and He do not completely segregate, the H abundance at the base of the accumulated layer is still reduced. While following the evolution of the X-ray burst is beyond the scope of this introductory paper, we note that the reduced proton-to-seed ratio favors the production of \iso{12}{C}--an important ingredient for subsequent superbursts.

We calculate the thermal and dynamical evolution of the surface layers of an accreting neutron star during the rise of a superburst. For the first few hours following unstable 12C ignition, the nuclear energy release is transported by convection. However, as the base temperature rises, the heating time becomes shorter than the eddy turnover time and convection becomes inefficient. This results in a hydrodynamic nuclear runaway, in which the heating time becomes shorter than the local dynamical time. Such hydrodynamic burning can drive shock waves into the surrounding layers and may be the trigger for the normal X-ray burst found to immediately precede the onset of the superburst in both cases where the Rossi X-Ray Timing Explorer was observing.

In this Paper we report on our analysis of a ~25 ksec. Chandra X-ray observation of the neutron star soft X-ray transient (SXT) 1H1905+000 in quiescence. Furthermore, we discuss our findings of the analysis of optical photometric observations which we obtained using the Magellan telescope and photometric and spectroscopic observations which we obtained using the Very Large Telescope at Paranal. The X-ray counterpart of 1H1905+000 was not detected in our Chandra data, with a 95 per cent confidence limit to the source count rate of 1.2x10^-4 counts s^-1. For different spectral models this yields an upper limit on the luminosity of 1.8x10^31 erg s^-1 (for an upper limit on the distance of 10 kpc.) This luminosity limit makes 1H1905+000 the faintest neutron star SXT in quiescence observed to date. The neutron star luminosity is so low that it is similar to the lowest luminosities derived for black hole SXTs in quiescence. This low luminosity for a neutron star SXT challanges the hypothesis presented in the literature that black hole SXTs in quiescence have lower luminosities than neutron star SXTs as a result of the presence of a black hole event horizon. Furthermore, the limit on the neutron star luminosity obtained less than 20 years after the outburst has ceased, constrains the thermal conductivity of the neutron star crust. Finally, the neutron star core must be so cold that unless the time averaged mass accretion rate is lower than 2x10^-12 M_sun yr^-1, core cooling has to proceed via enhanced neutrino emission processes. We derive a limit on the absolute I-band magnitude of the quiescent counterpart of M_I>7.8 assuming the source is at 10 kpc. This is in line with 1H1905+000 being an ultra-compact X-ray binary, as has been proposed based on the low outburst V-band absolute magnitude.

We present our first nucleosynthesis results from a numerical simulation of the thermonuclear disruption of a static cold Chandrasekhar-mass C/O white dwarf. The two-dimensional simulation was performed with an adaptive-mesh Eulerian hydrodynamics code, FLASH, that uses as a flame capturing scheme the evolution of a passive scaler. To compute the isotopic yields and their velocity distribution, 10,000 massless tracer particles are embedded in the star. The particles are advected along streamlines and provide a Lagrangian description of the explosion. We briefly describe our verification tests and preliminary results from post-processing the particle trajectories with a modest (214 isotopes) reaction network.

Superbursts are thought to be powered by the unstable ignition of a carbon-enriched layer formed from the burning of accreted hydrogen and helium. As shown by Cumming & Bildsten, the short recurrence time hinges on the crust being sufficiently hot at densities > 1.0e9 g/cc. In this Letter, we self-consistently solve for the flux coming from the deep crust and core. The temperature where the carbon unstably ignites is only weakly sensitive to the composition of the ashes of H/He burning, but does depend on the thermal conductivity of the inner crust and the neutrino emissivity of the core. The observed superburst recurrence times and energetics suggest that the crust thermal conductivity is low, as if the crust were amorphous instead of crystalline. If the conductivity is higher, such as from a lattice with impurities, then matching the superburst properties require that the neutrino emissivity be not stronger than modified Urca. Observations of superbursts-energetics, recurrence times, and cooling times-therefore complement observations of isolated cooling neutron stars and soft X-ray transients in constraining properties of dense matter. Perhaps the most interesting object in this regard is KS 1731-260, which produced a superburst during its protracted accretion outburst but had a rapidly declining quiescent luminosity.

This paper reports on the detailed rp-process reaction flow on an accreting neutron star and the resulting ashes of a type I X-ray burst. It is obtained by coupling a 298 isotope reaction network to a self-consistent one-dimensional model calculation with a constant accretion rate of dM/dt=1.0e17g/s (0.09 Eddington).

The light curves of type I X-ray bursts (XRBs) result from energy released from the atmosphere of a neutron star when accreted hydrogen and helium ignite and burn explosively via the rp-process. Since charged particle reaction rates are both density and very temperature dependent, a simulation model must provide accurate values of these variables to predict the reaction flow. This paper uses a self-consistent one-dimensional model calculation with a constant accretion rate of dM/dt=5e16g/s (0.045 Eddington) and reports on the detailed rp-process reaction flow of a given burst.

Many classical nova ejecta are enriched in CNO and Ne. Rosner et al. recently suggested that the enrichment might originate in the resonant interaction between large-scale shear flows in the accreted H/He envelope and gravity waves at the interface between the envelope and the underlying C/O white dwarf. The shear flow amplifies the waves, which eventually form cusps and break. This wave breaking injects a spray of C/O into the superincumbent H/He. Using two-dimensional simulations, we formulate a quantitative expression for the amount of C/O per unit area that can be entrained into the H/He at saturation. The fraction of the envelope that is enriched depends on the horizontal distribution of shear velocity and the density contrast between the C/O white dwarf and the H/He layer but is roughly independent of the vertical shape of the shear profile. Using this parameterization for the mixed mass, we then perform several one-dimensional Lagrangian calculations of an accreting white dwarf envelope and consider two scenarios: that the wave breaking and mixing is driven by the convective flows; and that the mixing occurs prior to the onset of convection. In the absence of enrichment prior to ignition, the base of the convective zone, as calculated from mixing-length theory with the Ledoux instability criterion, does not reach the C/O interface. As a result, there is no additional mixing, and the runaway is slow. In contrast, the formation of a mixed layer during the accretion of H/He, prior to ignition, causes a more violent runaway. The envelope can be enriched by <25% of C/O by mass (consistent with that observed in some ejecta) for shear velocities, over the surface, with Mach numbers <0.4.

We explore the idea that the observed variations in the peak luminosities of Type Ia supernovae originate in part from a scatter in metallicity of the main-sequence stars that become white dwarfs. Previous, numerical, studies have not self-consistently explored metallicities greater than solar. One-dimensional Chandrasekhar mass models of SNe Ia produce most of their 56Ni in a burn to nuclear statistical equilibrium between the mass shells 0.2 and 0.8 solar masses, for which the electron to nucleon ratio is constant during the burn. We show analytically that, under these conditions, charge and mass conservation constrain the mass of 56Ni produced to depend linearly on the original metallicity of the white dwarf progenitor. Detailed post-processing of W7-like models confirms this linear dependence. The effect that we identify is most evident at metallicities larger than solar, and is in agreement with previous self-consistent calculations over the metallicity range common to both calculations. The observed scatter in the metallicity (1/3--3 times solar) of the solar neighborhood is enough to induce a 25% variation in the mass of 56Ni ejected by Type Ia supernovae. This is sufficient to vary the peak V-band brightness by approximately 0.2. This scatter in metallicity is present out to the limiting redshifts of current observations (z < 1). Sedimentation of 22Ne can possibly amplify the variation in 56Ni mass up to 50%. Further numerical studies can determine if other metallicity-induced effects, such as a change in the mass of the 56Ni-producing region, offset or enhance this variation.

We obtained four Chandra/ACIS-S observations beginning two weeks after the end of the November 2000 outburst of the neutron star (NS) transient Aql X-1. Over the five month span in quiescence, the X-ray spectra are consistent with thermal emission from a NS with a pure hydrogen photosphere and R_{\infty}=15.9+{0.8}-{2.9} (d/5 kpc) km at the optically implied X-ray column density. We also detect a hard power-law tail during two of the four observations. The intensity of Aql X-1 first decreased by 50+/-4% over three months, then increased by 35+/-5% in one month, and then remained constant (<6% change) over the last month. These variations in the first two observations cannot be explained by a change in the power-law spectral component, nor in the X-ray column density. Presuming that R_{\infty} is not variable and a pure hydrogen atmosphere, the long-term changes can only be explained by variations in the NS effective temperature, from kT_{eff, \infty}=130+3-5 eV, down to 113+3-4 eV, finally increasing to 118+9-4 eV for the final two observations. During one of these observations, we observe two phenomena which were previously suggested as indicators of quiescent accretion onto the NS: short-timescale (<1e4 sec) variability (at 32+8-6% rms), and a possible absorption feature near 0.5 keV. The possible absorption feature can potentially be explained as due to a time-variable response in the ACIS detector. Even so, such a feature has not been detected previously from a NS, and if confirmed and identified, can be exploited for simultaneous measurements of the photospheric redshift and NS radius.

The composition of the outer 100 m of a neutron star sets the heat flux that flows outwards from the core. For an accreting neutron star in an X-ray transient, the thermal quiescent flux depends sensitively on the amount of hydrogen and helium remaining on the surface after an accretion outburst and on the composition of the underlying ashes of previous H/He burning. Because H/He has a higher thermal conductivity, a larger mass of H/He implies a shallower thermal gradient through the low density envelope and hence a higher effective temperature for a given core temperature. The mass of residual H and He varies from outburst to outburst, so the thermal quiescent flux is variable even though the core temperature is constant for timescales < 10 000 yr. Heavy elements settle from a H/He envelope in a few hours; we therefore model the quiescent envelope as two distinct layers, H/He over heavier elements, and treat the mass of H/He as a free parameter. We find that the emergent thermal quiescent flux can vary by a factor of 2 to 3 between different quiescent epochs. The variation is more pronounced at lower interior temperatures, making systems with low quiescent luminosities and frequent outbursts, such as SAX J1808.4-3658, ideal candidates from which to observe this effect. We compute, for different ash compositions, the interior temperatures of Cen X-4, Aql X-1, and SAX J1808.4-3658. In the case of Aql X-1, the inferred high interior temperature suggests that neutrino cooling contributes to the neutron star's thermal balance.

We present a detailed observational and theoretical study of a ~3 hr long X-ray burst (the ``super burst'') observed by the Rossi X-ray Timing Explorer (RXTE) from the low mass X-ray binary (LMXB) 4U 1820-30. This is the longest X-ray burst ever observed from this source, and perhaps one of the longest ever observed in great detail from any source. We show that the super burst is thermonuclear in origin. The level of the accretion driven flux as well as the total energy release of ~1.5 x 10^{42} ergs indicate that helium could not be the energy source for the super burst. We outline the physics relevant to carbon production and burning on helium accreting neutron stars and present calculations of the thermal evolution and stability of a carbon layer and show that this process is the most likely explanation for the super burst. We show that for large columns of accreted carbon fuel, a substantial fraction of the energy released in the carbon burning layer is radiated away as neutrinos, and the heat that is conducted from the burning layer in large part flows inward, only to be released on timescales longer than the observed burst. Thus the energy released possibly exceeds that observed in X-rays by more than a factor of ten. Spectral analysis during the super burst reveals the presence of a broad emission line between 5.8 - 6.4 keV and an edge at 8 - 9 keV likely due to reflection of the burst flux from the inner accretion disk in 4U 1820-30. We believe this is the first time such a signature has been unambiguously detected in the spectrum of an X-ray burst.

(Abridged). The type-I X-ray bursting low mass X-ray binary KS 1731-260 was recently detected for the first time in quiescence by Wijnands et al., following an approximately 13 yr outburst which ended in Feb 2001. Unlike all other known transient neutron stars, the duration of this recent outburst is as long as the thermal diffusion time of the crust. The large amount of heat deposited by reactions in the crust will have heated the crust to temperatures much higher than the equilibrium core temperature. As a result, the thermal luminosity currently observed from the neutron star is dominated not by the core, but by the crust. Moreover, the level and the time evolution of quiescent luminosity is determined mostly by the amount of heat deposited in the crust during the most recent outburst. Using estimates of the outburst mass accretion rate, our calculations of the quiescent flux immediately following the end of the outburst agree with the observed quiescent flux to within a factor of a few. We present simulations of the evolution of the quiescent lightcurve for different scenarios of the crust microphysics, and demonstrate that monitoring observations (with currently flying instruments) spanning from 1--30 yr can measure the crust cooling timescale and the total amount of heat stored in the crust. These quantities have not been directly measured for any neutron star.

Using the Chandra/ACIS-I detector, we have identified an X-ray source (CXOU 132619.7-472910.8) in the globular cluster NGC 5139 with a thermal spectrum identical to that observed from transiently accreting neutron stars in quiescence. The absence of intensity variability on timescales as short as 4 seconds (< 25% rms variability) and as long as 5 years (<50% variability) supports the identification of this source as a neutron star, most likely maintained at a high effective temperature (approximately 1e6 K) by transient accretion from a binary companion. The ability to spectrally identify quiescent neutron stars in globular clusters (where the distance and interstellar column densities are known) opens up new opportunities for precision neutron star radius measurements.

We report on the quiescent spectrum measured with Chandra/ACIS-S of the transient, type-I X-ray bursting neutron star Aql X-1, immediately following an accretion outburst. The neutron star radius, assuming a pure hydrogen atmosphere and hard power-law spectrum, is $R_\infty$=13.4{+5}{-4} (d/5 \kpc) km. Based on the historical outburst record of RXTE/ASM, the quiescent luminosity is consistent with that predicted by Brown, Bildsten and Rutledge from deep crustal heating, lending support to this theory for providing a minimum quiescent luminosity of transient neutron stars. While not required by the data, the hard power-law component can account for 18+/-8% of the 0.5-10 keV thermal flux. Short-timescale intensity variability during this observation is less than 15% rms (3 sigma; 0.0001-1 Hz, 0.2-8 keV). Comparison between the Chandra spectrum and three X-ray spectral observations made between Oct 1992 and Oct 1996 find all spectra consistent with a pure H atmosphere, but with temperatures ranging from 145--168 eV, spanning a factor of 1.87+/-0.21 in observed flux. The source of variability in the quiescent luminosity on long timescales (greater than years) remains a puzzle. If from accretion, then it remains to be explained why the quiescent accretion rate provides a luminosity so nearly equal to that from deep crustal heating.

We assess the stability of Kruskal-Schwarzschild (magnetic Rayleigh-Taylor) type modes for accreted matter on the surface of a neutron star confined by a strong (>= 1.E12 G) magnetic field. Employing the energy principle to analyze the stability of short-wavelength ballooning modes, we find that line-tying to the neutron star crust stabilizes these modes until the overpressure at the top of the neutron star crust exceeds the magnetic pressure by a factor ~ 8(a/h), where a and h are respectively the lateral extent of the accretion region and the density scale height. The most unstable modes are localized within a density scale height above the crust. We calculate the amount of mass that can be accumulated at the polar cap before the onset of instability.

We report on spectral and intensity variability analysis from a Chandra/ACIS-S observation of the transient, type-I X-ray bursting low-mass X-ray binary Cen X-4. The quiescent X-ray spectrum during this observation is statistically identical to one observed previously with Beppo/SAX, and close, but not identical, to one observed previously with ASCA. The X-ray spectrum is best described as a pure Hydrogen atmosphere thermal spectrum plus a power-law component that dominates the spectrum above 2 keV. The best-fit radius of the neutron star is r=12.9+/-2.6 (d/1.2 kpc) km if the interstellar absorption is fixed at the value implied by the optical reddening. Allowing the interstellar absorption to be a free parameter yields r=19+45-10 (d/1.2 kpc) km (90% confidence). The thermal spectrum from the neutron star surface is inconsistent with a solar metallicity. We find a 3sigma upper-limit of root-mean-square variability <18% (0.2-2.0 keV; 0.0001-1 Hz) during the observation. On the other hand, the 0.5-10.0 keV luminosity decreased by 40+/-8% in the 4.9 years between the Asca and Chandra observations. This variability can be attributed to the power-law component. Moreover, we limit the variation in thermal temperature to <10% over these 4.9 years. The stability of the thermal temperature and emission area radius supports the interpretation that the quiescent thermal emission is due to the hot neutron star core.

Recent observations suggest that neutron stars in low-mass X-ray binaries rotate within a narrow range of spin frequencies clustered around 300 Hz. A proposed explanation for this remarkable fact is that gravitational radiation from a steady-state r-mode oscillation in the neutron star's core halts the spin-up due to accretion. For the neutron star transients, balancing the time-averaged accretion torque with gravitational wave emission from steady-state, constant amplitude r-mode pulsations implies a quiescent luminosity too bright to be consistent with observations (in particular of Aql X-1). The viscous dissipation (roughly 10 MeV per accreted nucleon for a spin of 300 Hz) from such an r-mode makes the core sufficiently hot to power a thermal luminosity of order 1e34 erg/s when accretion halts. This is the minimum quiescent luminosity that the neutron star must emit when viscous heating in the core is balanced by radiative cooling from the surface, as is the case when the core of the star is superfluid. We therefore conclude that either the accretion torque is much less than (dM/dt)(GMR)^{1/2}, or that a steady-state r-mode does not limit the spin rate of the neutron star transients. Future observations with AXAF and XMM promise to further constrain the amount of viscous dissipation in the neutron star core.

A neutron star in a long-lived, low-mass binary can easily accrete enough matter to replace its entire crust. Previous authors noted that an accreted crust, being formed from the burning of accreted hydrogen and helium, allows a series of non-equilibrium reactions, at densities >6e11 g/cc, which release a substantial amount of heat (1 MeV per accreted nucleon). Recent calculations by Schatz et al. showed that the crystalline lattice of an accreted crust is also likely to be quite impure. This paper discusses the thermal structure of such a neutron star and surveys how the crust reactions and impurities affect the crust temperature. During accretion rapid enough to make the accreted hydrogen and helium burn stably (near the Eddington accretion rate; typical of the brightest low-mass neutron star binaries), most of the heat released in the crust is conducted into the core, where neutrino emission regulates the temperature. As a result there is an inversion of the thermal gradient: the temperature decreases with depth in the inner crust. The thermal structure in the crust at these high accretion rates is insensitive to the temperature in the hydrogen/helium burning shell. When the crust is very impure, the temperature can reach approximately 8e8 K at densities > 6e11 g/cc. This peak temperature depends mostly on the amount of heat released and the thermal conductivity and in particular is roughly independent of the core temperature. The high crust temperatures are sufficient to melt the crystalline lattice in thin layers where electron captures have substantially reduced the nuclear charge.

We fit hydrogen atmosphere models to the X-ray data for four neutron stars (three from a previous paper, plus 4U 2129+47) and six black hole candidates (A0620-00, GS 2000+25, GS 1124-68, GS 2023+33, GRO J1655-40, and GRO J0422+32). While the neutron stars are similar in their intrinsic X-ray spectra (similar effective temperatures and emission area radii ~10 km), the spectra of two black hole candidates are significantly different, and the spectra of the remaining four are consistent with a very large parameter space that includes the neutron stars. The spectral differences between the neutron stars and black hole candidates favors the interpretation that the quiescent neutron star emission is predominantly thermal emission from the neutron star surface. Our work suggests that an X-ray spectral comparison in quiescence provides an additional means for distinguishing between neutron stars and black holes. The faint X-ray sources in globular clusters are also a class of objects which can be investigated in this manner.

We re-analyze the available X-ray spectral data of the type I bursting neutron star transients Aql X-1, Cen X-4, and 4U 1608-522 using realistic hydrogen atmosphere models. Previous spectral fits assumed a blackbody spectrum; because the free-free dominated photospheric opacity decreases with increasing frequency, blackbody spectral fits overestimate the effective temperature and underestimate, by as much as two orders of magnitude, the emitting area. Hydrogen atmosphere spectral models, when fit to the available observational data, imply systematically larger emission area radii, consistent with the canonical 10 km radius of a neutron star. This suggests that a substantial fraction of the quiescent luminosity is thermal emission from the surface of the neutron star. The magnitude of the equivalent hydrogen column density toward these systems, however, presents a considerable systematic uncertainty, which can only be eliminated by high signal-to-noise X-ray spectral measurements (e.g., with AXAF or XMM) which would permit simultaneous determination of the equivalent hydrogen column density, emission area, and thermal temperature.

Nuclear reactions occurring deep in the crust of a transiently accreting neutron star efficiently maintain the core at a temperature >5e7 K. When accretion halts, the envelope relaxes to a thermal equilibrium set by the flux from the hot core, as if the neutron star were newly born. For the time-averaged accretion rates typical of low-mass X-ray transients, standard neutrino cooling is unimportant and the core thermally re-radiates the deposited heat. The resulting luminosity has the same magnitude as that observed from several transient neutron stars in quiescence. Confirmation of this mechanism would strongly constrain rapid neutrino cooling mechanisms for neutron stars. Thermal emission had previously been dismissed as a predominant source of quiescent emission since blackbody spectral fits implied an emitting area much smaller than a neutron star's surface. However, as with thermal emission from radio pulsars, fits with realistic emergent spectra will imply a substantially larger emitting area. Other emission mechanisms, such as accretion or a pulsar shock, can also operate in quiescence and generate intensity and spectral variations over short timescales. Indeed, quiescent accretion may produce gravitationally redshifted metal photoionization edges in the quiescent spectra (detectable with AXAF and XMM). We discuss past observations of Aql~X-1 and note that the low luminosity X-ray sources in globular clusters and the Be star/X-ray transients are excellent candidates for future study.

Abundance measurements of the light elements lithium, beryllium, and boron are playing an increasingly important role in the study of stellar physics. Because these elements are easily destroyed in stars at temperatures 2--4 million K, the abundances in the surface convective zone are diagnostics of the star's internal workings. Standard stellar models cannot explain depletion patterns observed in low mass stars, and so are not accounting for all the relevant physical processes. These processes have important implications for stellar evolution and primordial lithium production in big bang nucleosynthesis. Because beryllium is destroyed at slightly higher temperatures than lithium, observations of both light elements can differentiate between the various proposed depletion mechanisms. Unfortunately, the reaction rate for the main destruction channel, 9Be(p,alpha)6Li, is uncertain. A level in the compound nucleus 10B is only 25.7 keV below the reaction's energetic threshold. The angular momentum and parity of this level are not well known; current estimates indicate that the resonance entrance channel is either s- or d-wave. We show that an s-wave resonance can easily increase the reaction rate by an order of magnitude at temperatures of approximately 4 million K. Observations of sub-solar mass stars can constrain the strength of the resonance, as can experimental measurements at lab energies lower than 30 keV.

We investigate the atmosphere, ocean, and crust of neutron stars accreting at rates sufficiently high (typically in excess of the local Eddington limit) to stabilize the burning of accreted hydrogen and helium. For hydrogen-rich accretion at global rates in excess of 10^-8 solar masses per year (typical of a few neutron stars), we discuss the thermal state of the deep ocean and crust and their coupling to the neutron star core, which is heated by conduction (from the crust) and cooled by neutrino emission. We estimate the Ohmic diffusion time in the hot, deep crust and find that it is noticeably shortened (to less than 10^8 yr) from the values characteristic of the colder crusts in slowly accreting neutron stars. We speculate on the implications of these calculations for magnetic field evolution in the bright accreting X-ray sources. We also explore the consequences of rapid compression at local accretion rates exceeding ten times the Eddington rate. This rapid accretion heats the atmosphere/ocean to temperatures of order 10^9 K at relatively low densities; for stars accreting pure helium, this causes unstable ignition of the ashes (mostly carbon) resulting from stable helium burning. This unstable burning can re-occur on timescales as short as hours to days, and might be the cause of some flares on helium accreting pulsars, in particular 4U~1626--67. Such rapid local accretion rates are common on accreting X-ray pulsars, where the magnetic field focuses the accretion flow onto a small fraction of the stellar area.

We present an analytic calculation of the thermonuclear depletion of the light elements lithium, beryllium, and boron in fully convective, low-mass stars. Under the presumption that the pre--main-sequence star is always fully mixed during contraction, we find that the burning of these rare light elements can be computed analytically, even when the star is degenerate. Using the effective temperature as a free parameter, we constrain the properties of low-mass stars from observational data, independently of the uncertainties associated with modeling their atmospheres and convection. Our analytic solution explains the dependence of the age at a given level of elemental depletion on the stellar effective temperature, nuclear cross sections, and chemical composition. Most importantly, our results allow observers to translate lithium non-detections in young cluster members into a model-independent minimum age for that cluster. Using this procedure, we have found lower limits to the ages of the Pleiades (100 Myr) and Alpha Persei (60 Myr) clusters. Recent experimental work on the low energy resonance in the ^10B(p,\alpha)^7Be reaction has greatly enhanced estimates of the destruction rate of ^10B, making it possible for stars with M>0.1 M_sun to deplete both ^10B and ^11B before reaching the main sequence. Moreover, there is an interesting range of masses, 0.085 M_sun < M < 0.13 M_sun, where boron depletion occurs on the main sequence in less than a Hubble time, providing a potential ``clock'' for dating low-mass stars.

We present an analytic calculation of the thermonuclear depletion of lithium in contracting, fully convective, pre-main sequence stars of mass M < 0.5 M_sun. Previous numerical work relies on still-uncertain physics (atmospheric opacities and convection, in particular) to calculate the effective temperature as a unique function of stellar mass. We assume that the star's effective temperature, T_eff, is fixed during Hayashi contraction and allow its actual value to be a free parameter constrained by observation. Using this approximation, we compute lithium burning analytically and explore the dependence of lithium depletion on T_eff, M, and composition. Our calculations yield the radius, age, and luminosity of a pre-main sequence star as a function of lithium depletion. This allows for more direct comparisons to observations of lithium depleted stars. Our results agree with those numerical calculations that explicitly determine stellar structure during Hayashi contraction. In agreement with Basri, Marcy, and Graham (1996), we show that the absence of lithium in the Pleiades star HHJ 3 implies that it is older than 100 Myr. We also suggest a generalized method for dating galactic clusters younger than 100 Myr (i.e., those with contracting stars of M > 0.08 M_sun) and for constraining the masses of lithium depleted stars.

We investigate the thermal stability of nuclear burning on the accreting X-ray pulsar GRO J1744-28. The neutron star's dipolar magnetic field is <3\times 10^{11} G if persistent spin-up implies that the magnetospheric radius is less than the co-rotation radius. After inferring the properties of the neutron star, we study the thermal stability of hydrogen/helium burning and show that thermonuclear instabilities are unlikely causes of the hourly bursts seen at very high accretion rates. We then discuss how the stability of the thermonuclear burning depends on both the global accretion rate and the neutron star's magnetic field strength. We emphasize that the appearance of the instability (i.e., whether it looks like a Type I X-ray burst or a flare lasting a few minutes) will yield crucial information on the neutron star's surface magnetic field and the role of magnetic fields in convection. We suggest that a thermal instability in the accretion disk is the origin of the long (~300 days) outburst and that the recurrence time of these outbursts is >50 years. We also discuss the nature of the binary and point out that a velocity measurement of the stellar companion (most likely a Roche-lobe filling giant with m_K>17) will constrain the neutron star mass.