Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.

Radioactive isotopes produced in core-collapse supernovae (CCSNe) provide useful insights into the underlying processes driving the collapse mechanism and the origins of elemental abundances. Their study generates a confluence of major physics research, including experimental measurements of nuclear reaction rates, astrophysical modeling, and $\gamma$-ray observations. Here we identify the key nuclear reaction rates to the nucleosynthesis of observable radioactive isotopes in explosive silicon-burning during CCSNe. Using the nuclear reaction network calculator SkyNet and current REACLIB reaction rates, we evolve temperature-density-time profiles of the innermost $0.45~M_\odot$ ejecta from the core collapse and explosion of a $12~M_\odot$ star. Individually varying 3403 reaction rates by factors of 100, we identify 141 reactions which cause significant differences in the isotopes of interest, namely, $^{43}$K, $^{47}$Ca, $^{44,47}$Sc, $^{44}$Ti, $^{48,51}$Cr, $^{48,49}$V, $^{52,53}$Mn, $^{55,59}$Fe, $^{56,57}$Co, and $^{56,57,59}$Ni. For each of these reactions, we present a novel method to extract the temperature range pertinent to the nucleosynthesis of the relevant isotope; the resulting temperatures lie within the range $T = 0.47$ to $6.15~$GK. Limiting the variations to within $1\sigma$ of STARLIB reaction rate uncertainties further reduces the identified reactions to 48 key rates, which can be used to guide future experimental research. Complete results are presented in tabular form.

We report mass measurements of neutron-rich Ga isotopes $^{80-85}$Ga with TRIUMF's Ion Trap for Atomic and Nuclear science (TITAN). The measurements determine the masses of $^{80-83}$Ga in good agreement with previous measurements. The masses of $^{84}$Ga and $^{85}$Ga were measured for the first time. Uncertainties between $25-48$ keV were reached. The new mass values reduce the nuclear uncertainties associated with the production of A $\approx$ 84 isotopes by the \emph{r}-process for astrophysical conditions that might be consistent with a binary neutron star (BNS) merger producing a blue kilonova. Our nucleosynthesis simulations confirm that BNS merger may contribute to the first abundance peak under moderate neutron-rich conditions with electron fractions $Y_e=0.35-0.38$.

Our understanding of the core-collapse supernova explosion mechanism is incomplete. While the favoured scenario is delayed revival of the stalled shock by neutrino heating, it is difficult to reliably compute explosion outcomes and energies, which depend sensitively on the complex radiation hydrodynamics of the post-shock region. The dynamics of the (non-)explosion depend sensitively on how energy is transported from inside and near the proto-neutron star (PNS) to material just behind the supernova shock. Although most of the PNS energy is lost in the form of neutrinos, hydrodynamic and hydromagnetic waves can also carry energy from the PNS to the shock. We show that gravity waves excited by core PNS convection can couple with outgoing acoustic waves that present an appreciable source of energy and pressure in the post-shock region. Using one-dimensional simulations, we estimate the gravity wave energy flux excited by PNS convection and the fraction of this energy transmitted upward to the post-shock region as acoustic waves. We find wave energy fluxes near $10^{51}\,\mathrm{erg}\,\mathrm{s}^{-1}$ are likely to persist for $\sim1\,\mathrm{s}$ post-bounce. The wave pressure on the shock may exceed $10\%$ of the thermal pressure, potentially contributing to shock revival and, subsequently, a successful and energetic explosion. We also discuss how future simulations can better capture the effects of waves, and more accurately quantify wave heating rates.

Uncertainties in our knowledge of the properties of dense matter near and above nuclear saturation density are among the main sources of variations in multi-messenger signatures predicted for core-collapse supernovae (CCSNe) and the properties of neutron stars (NSs). We construct 97 new finite-temperature equations of state (EOSs) of dense matter that obey current experimental, observational, and theoretical constraints and discuss how systematic variations in the EOS parameters affect the properties of cold nonrotating NSs and the core collapse of a $20\,M_\odot$ progenitor star. The core collapse of the $20\,M_\odot$ progenitor star is simulated in spherical symmetry using the general-relativistic radiation-hydrodynamics code GR1D where neutrino interactions are computed for each EOS using the NuLib library. We conclude that the effective mass of nucleons at densities above nuclear saturation density is the largest source of uncertainty in the CCSN neutrino signal and dynamics even though it plays a subdominant role in most properties of cold NS matter. Meanwhile, changes in other observables affect the properties of cold NSs, while having little effect in CCSNe. To strengthen our conclusions, we perform six octant three-dimensional CCSN simulations varying the effective mass of nucleons at nuclear saturation density. We conclude that neutrino heating and, thus, the likelihood of explosion is significantly increased for EOSs where the effective mass of nucleons at nuclear saturation density is large.

We present a systematic numerical relativity study of the mass ejection and the associated electromagnetic transients and nucleosynthesis from binary neutron star (NS) mergers. We find that a few $10^{-3}\, M_\odot$ of material are ejected dynamically during the mergers. The amount and the properties of these outflow depend on binary parameters and on the NS equation of state (EOS). A small fraction of these ejecta, typically ${\sim}10^{-6}\, M_\odot$, is accelerated by shocks formed shortly after merger to velocities larger than $0.6\, {\rm c}$ and produces bright radio flares on timescales of weeks, months, or years after merger. Their observation could constrain the strength with which the NSs bounce after merger and, consequently, the EOS of matter at extreme densities. The dynamical ejecta robustly produce second and third $r$-process peak nuclei with relative isotopic abundances close to solar. The production of light $r$-process elements is instead sensitive to the binary mass ratio and the neutrino radiation treatment. Accretion disks of up to ${\sim}0.2\, M_\odot$ are formed after merger, depending on the lifetime of the remnant. In most cases, neutrino- and viscously-driven winds from these disks dominate the overall outflow. Finally, we generate synthetic kilonova light curves and find that kilonovae depend on the merger outcome and could be used to constrain the NS EOS.

Over the past three decades instruments on the ground and in space have discovered thousands of planets outside the solar system. These observations have given rise to an astonishingly detailed picture of the demographics of short-period planets, but are incomplete at longer periods where both the sensitivity of transit surveys and radial velocity signals plummet. Even more glaring is that the spectra of planets discovered with these indirect methods are either inaccessible (radial velocity detections) or only available for a small subclass of transiting planets with thick, clear atmospheres. Direct detection can be used to discover and characterize the atmospheres of planets at intermediate and wide separations, including non-transiting exoplanets. Today, a small number of exoplanets have been directly imaged, but they represent only a rare class of young, self-luminous super-Jovian-mass objects orbiting tens to hundreds of AU from their host stars. Atmospheric characterization of planets in the <5 AU regime, where radial velocity (RV) surveys have revealed an abundance of other worlds, is technically feasible with 30-m class apertures in combination with an advanced AO system, coronagraph, and suite of spectrometers and imagers. There is a vast range of unexplored science accessible through astrometry, photometry, and spectroscopy of rocky planets, ice giants, and gas giants. In this whitepaper we will focus on one of the most ambitious science goals --- detecting for the first time habitable-zone rocky (<1.6 R_Earth) exoplanets in reflected light around nearby M-dwarfs

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.

General-relativistic simulations of binary neutron star mergers with viscosity reveal a new outflow mechanism operating in unequal mass binaries on dynamical timescales and enabled by turbulent viscosity. These "viscous-dynamical" ejecta are launched during merger due to the thermalization of mass exchange streams between the secondary and the primary neutron star. They are characterized by asymptotic velocities extending up to ${\sim} 0.8\, c$, and have masses that depend on the efficiency of the viscous mechanism. Depending on the unknown strength of the effective viscosity arising from magnetohydrodynamics instabilities operating during merger, the overall mass of the dynamical ejecta could be enhanced by a factor of a few and the mass of the fast tail of the ejecta having asymptotic velocities $\geq 0.6\, c$ by up to four orders of magnitude. The radioactive decay of the expanding viscous-dynamical ejecta could produce bright kilonova transients with signatures of free neutron decay in the first hour and enhanced near infrared flux on a timescale of a few days. The synchrotron remnant produced by the interaction between the ejecta and the interstellar medium could also be significantly enhanced by viscosity. Such remnant could be detected in the case of GW170817 as a rebrightening of the radio signal in the next months to years.

In July 2018 an FRIB Theory Alliance program was held on the implications of GW170817 and its associated kilonova for r-process nucleosynthesis. Topics of discussion included the astrophysical and nuclear physics uncertainties in the interpretation of the GW170817 kilonova, what we can learn about the astrophysical site or sites of the r process from this event, and the advances in nuclear experiment and theory most crucial to pursue in light of the new data. Here we compile a selection of scientific contributions to the workshop, broadly representative of progress in r-process studies since the GW170817 event.

Direct detection, also known as direct imaging, is a method for discovering and characterizing the atmospheres of planets at intermediate and wide separations. It is the only means of obtaining spectra of non-transiting exoplanets. Characterizing the atmospheres of planets in the <5 AU regime, where RV surveys have revealed an abundance of other worlds, requires a 30-m-class aperture in combination with an advanced adaptive optics system, coronagraph, and suite of spectrometers and imagers - this concept underlies planned instruments for both TMT (the Planetary Systems Imager, or PSI) and the GMT (GMagAO-X). These instruments could provide astrometry, photometry, and spectroscopy of an unprecedented sample of rocky planets, ice giants, and gas giants. For the first time habitable zone exoplanets will become accessible to direct imaging, and these instruments have the potential to detect and characterize the innermost regions of nearby M-dwarf planetary systems in reflected light. High-resolution spectroscopy will not only illuminate the physics and chemistry of exo-atmospheres, but may also probe rocky, temperate worlds for signs of life in the form of atmospheric biomarkers (combinations of water, oxygen and other molecular species). By completing the census of non-transiting worlds at a range of separations from their host stars, these instruments will provide the final pieces to the puzzle of planetary demographics. This whitepaper explores the science goals of direct imaging on 30-m telescopes and the technology development needed to achieve them.

We summarize the red channel (2-5 micron) of the Planetary Systems Imager (PSI), a proposed second-generation instrument for the TMT. Cold exoplanets emit the majority of their light in the thermal infrared, which means these exoplanets can be detected at a more modest contrast than at other wavelengths. PSI-Red will be able to detect and characterize a wide variety of exoplanets, including radial-velocity planets on wide orbits, accreting protoplanets in nearby star-forming regions, and reflected-light planets around the nearest stars. PSI-Red will feature an imager, a low-resolution lenslet integral field spectrograph, a medium-resolution lenslet+slicer integral field spectrograph, and a fiber-fed high-resolution spectrograph.

Magnetohydrodynamic (MHD) turbulence is of key importance in many high-energy astrophysical systems, where MHD instabilities can amplify local magnetic field over very short time scales. Specifically, the magnetorotational instability (MRI) and dynamo action have been suggested as a mechanism to grow magnetar-strength magnetic field ($\ge10^{15} G$) and magnetorotationally power the explosion of a rotating massive star. Such stars are progenitor candidates for type Ic-bl hypernova explosions and make up all supernovae connected to long gamma-ray bursts (GRBs). The MRI has been studied with local high-resolution shearing box simulations in 3D or with global 2D simulations, but it is an open question whether MRI-driven turbulence can result in the creation of a large-scale ordered and dynamically relevant field. Here we report results from global 3D general-relativistic magnetohydrodynamic (GRMHD) turbulence simulations and show that MRI-driven MHD turbulence in rapidly rotating protoneutron stars produces an inverse cascade of energy. We find a large-scale ordered toroidal field that is consistent with the formation of bipolar magnetorotationally driven outflows. Our results demonstrate that rapidly rotating massive stars are plausible progenitors for both type Ic-bl supernovae and long GRBs, present a viable formation scenario for magnetars, and may account for potentially magnetar-powered superluminous supernovae.

Multidimensional simulations show that non-radial, turbulent, fluid motion is a fundamental component of the core-collapse supernova (CCSN) explosion mechanism. Neutrino-driven convection, the standing accretion shock instability, and relic-perturbations from advanced stages of nuclear burning can all impact the outcome of core collapse in a qualitative and quantitative way. Here, we review the current understanding of these phenomena and their role in the explosion of massive stars. We also discuss the role of protoneutron star convection and of magnetic fields in the context of the delayed neutrino mechanism.

Almost all of the elements heavier than hydrogen that are present in our solar system were produced by nuclear burning processes either in the early universe or at some point in the life cycle of stars. In all of these environments, there are dozens to thousands of nuclear species that interact with each other to produce successively heavier elements. In this paper, we present SkyNet, a new general-purpose nuclear reaction network that evolves the abundances of nuclear species under the influence of nuclear reactions. SkyNet can be used to compute the nucleosynthesis evolution in all astrophysical scenarios where nucleosynthesis occurs. SkyNet is free and open-source and aims to be easy to use and flexible. Any list of isotopes can be evolved and SkyNet supports various different types of nuclear reactions. SkyNet is modular so that new or existing physics, like nuclear reactions or equations of state, can easily be added or modified. Here, we present in detail the physics implemented in SkyNet with a focus on a self-consistent transition to and from nuclear statistical equilibrium (NSE) to non-equilibrium nuclear burning, our implementation of electron screening, and coupling of the network to an equation of state. We also present comprehensive code tests and comparisons with existing nuclear reaction networks. We find that SkyNet agrees with published results and other codes to an accuracy of a few percent. Discrepancies, where they exist, can be traced to differences in the physics implementations.

We investigate r-process nucleosynthesis in three-dimensional (3D) general-relativistic magnetohydrodynamic simulations of rapidly rotating strongly magnetized core collapse. The simulations include a microphysical finite-temperature equation of state and a leakage scheme that captures the overall energetics and lepton number exchange due to postbounce neutrino emission and absorption. We track the composition of the ejected material using the nuclear reaction network SkyNet. Our results show that the 3D dynamics of magnetorotational core-collapse supernovae (CCSN) are important for their nucleosynthetic signature. We find that production of r-process material beyond the second peak is reduced by a factor of 100 when the magnetorotational jets produced by the rapidly rotating core undergo a kink instability. Our results indicate that 3D magnetorotationally powered CCSNe are a robust r-process source only if they are obtained by the collapse of cores with unrealistically large precollapse magnetic fields of order $10^{13}$G. Additionally, a comparison simulation that we restrict to axisymmetry, results in overly optimistic r-process production for lower magnetic field strengths.

We present a first study of the progenitor star dependence of the three-dimensional (3D) neutrino mechanism of core-collapse supernovae. We employ full 3D general-relativistic multi-group neutrino radiation-hydrodynamics and simulate the post-bounce evolutions of progenitors with zero-age main sequence masses of $12$, $15$, $20$, $27$, and $40\,M_\odot$. All progenitors, with the exception of the $12\,M_\odot$ star, experience shock runaway by the end of their simulations. In most cases, a strongly asymmetric explosion will result. We find three qualitatively distinct evolutions that suggest a complex dependence of explosion dynamics on progenitor density structure, neutrino heating, and 3D flow. (1) Progenitors with massive cores, shallow density profiles, and high post-core-bounce accretion rates experience very strong neutrino heating and neutrino-driven turbulent convection, leading to early shock runaway. Accretion continues at a high rate, likely leading to black hole formation. (2) Intermediate progenitors experience neutrino-driven, turbulence-aided explosions triggered by the arrival of density discontinuities at the shock. These occur typically at the silicon/silicon-oxygen shell boundary. (3) Progenitors with small cores and density profiles without strong discontinuities experience shock recession and develop the 3D standing-accretion shock instability (SASI). Shock runaway ensues late, once declining accretion rate, SASI, and neutrino-driven convection create favorable conditions. These differences in explosion times and dynamics result in a non-monotonic relationship between progenitor and compact remnant mass.

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

We investigate the nucleosynthesis of heavy elements in the winds ejected by accretion disks formed in neutron star mergers. We compute the element formation in disk outflows from hypermassive neutron star (HMNS) remnants of variable lifetime, including the effect of angular momentum transport in the disk evolution. We employ long-term axisymmetric hydrodynamic disk simulations to model the ejecta, and compute r-process nucleosynthesis with tracer particles using a nuclear reaction network containing $\sim 8000$ species. We find that the previously known strong correlation between HMNS lifetime, ejected mass, and average electron fraction in the outflow is directly related to the amount of neutrino irradiation on the disk, which dominates mass ejection at early times in the form of a neutrino-driven wind. Production of lanthanides and actinides saturates at short HMNS lifetimes ($\lesssim 10$ ms), with additional ejecta contributing to a blue optical kilonova component for longer-lived HMNSs. We find good agreement between the abundances from the disk outflow alone and the solar r-process distribution only for short HMNS lifetimes ($\lesssim 10$ ms). For longer lifetimes, the rare-earth and third r-process peaks are significantly under-produced compared to the solar pattern, requiring additional contributions from the dynamical ejecta. The nucleosynthesis signature from a spinning black hole (BH) can only overlap with that from a HMNS of moderate lifetime ($\lesssim 60$ ms). Finally, we show that angular momentum transport not only contributes with a late-time outflow component, but that it also enhances the neutrino-driven component by moving material to shallower regions of the gravitational potential, in addition to providing additional heating.

The equation of state (EOS) of dense matter is an essential ingredient for numerical simulations of core-collapse supernovae and neutron star mergers. The properties of matter near and above nuclear saturation density are uncertain, which translates into uncertainties in astrophysical simulations and their multi-messenger signatures. Therefore, a wide range of EOSs spanning the allowed range of nuclear interactions are necessary for determining the sensitivity of these astrophysical phenomena and their signatures to variations in input microphysics. We present a new set of finite temperature EOSs based on experimentally allowed Skyrme forces. We employ a liquid drop model of nuclei to capture the non-uniform phase of nuclear matter at sub-saturation density, which is blended into a nuclear statistical equilibrium EOS at lower densities. We also provide a new, open-source code for calculating EOSs for arbitrary Skyrme parametrizations. We then study the effects of different Skyrme parametrizations on thermodynamical properties of dense astrophysical matter, the neutron star mass-radius relationship, and the core collapse of 15 and 40 solar mass stars.

We investigate the ejecta from black hole - neutron star mergers by modeling the formation and interaction of mass ejected in a tidal tail and a disk wind. The outflows are neutron-rich, giving rise to optical/infrared emission powered by the radioactive decay of $r$-process elements (a kilonova). Here we perform an end-to-end study of this phenomenon, where we start from the output of a fully-relativistic merger simulation, calculate the post-merger hydrodynamical evolution of the ejecta and disk winds including neutrino physics, determine the final nucleosynthetic yields using post-processing nuclear reaction network calculations, and compute the kilonova emission with a radiative transfer code. We study the effects of the tail-to-disk mass ratio by scaling the tail density. A larger initial tail mass results in fallback matter becoming mixed into the disk and ejected in the subsequent disk wind. Relative to the case of a disk without dynamical ejecta, the combined outflow has lower mean electron fraction, faster speed, larger total mass, and larger absolute mass free of high-opacity Lanthanides or Actinides. In most cases, the nucleosynthetic yield is dominated by the heavy $r$-process contribution from the unbound part of the tidal tail. A Solar-like abundance distribution can however be obtained when the total mass of the dynamical ejecta is comparable to the mass of the disk outflows. The kilonova has a characteristic duration of 1 week and a luminosity of ~$10^{41}$ erg/s, with orientation effects leading to variations of a factor ~2 in brightness. At early times (< 1 day) the emission includes an optical component from the (hot) Lanthanide-rich material, but the spectrum evolves quickly to the infrared thereafter.

We present a proof-of-concept study, based on numerical-relativity simulations, of how gravitational waves (GWs) from neutron star merger remnants can probe the nature of matter at extreme densities. Phase transitions and extra degrees of freedom can emerge at densities beyond those reached during the inspiral, and typically result in a softening of the equation of state (EOS). We show that such physical effects change the qualitative dynamics of the remnant evolution, but they are not identifiable as a signature in the GW frequency, with the exception of possible black-hole formation effects. The EOS softening is, instead, encoded in the GW luminosity and phase and is in principle detectable up to distances of the order of several Mpcs with advanced detectors and up to hundreds of Mpcs with third generation detectors. Probing extreme-density matter will require going beyond the current paradigm and developing a more holistic strategy for modeling and analyzing postmerger GW signals.

We present the first near infra-red spectrum of the faint white dwarf companion HD 114174 B, obtained with Project 1640. Our spectrum, covering the Y, J and H bands, combined with previous TRENDS photometry measurements, allows us to place further constraints on this companion. We suggest two possible scenarios; either this object is an old, low mass, cool H atmosphere white dwarf with Teff ~ 3800 K or a high mass white dwarf with Teff > 6000 K, potentially with an associated cool (Teff ~ 700 K) brown dwarf or debris disk resulting in an infra-red excess in the L band. We also provide an additional astrometry point for 2014 June 12 and use the modelled companion mass combined with the RV and direct imaging data to place constraints on the orbital parameters for this companion.

We report on a set of long-term general-relativistic three-dimensional (3D) multi-group (energy-dependent) neutrino-radiation hydrodynamics simulations of core-collapse supernovae. We employ a full 3D two-moment scheme with the local M1 closure, three neutrino species, and 12 energy groups per species. With this, we follow the post-core-bounce evolution of the core of a nonrotating $27$-$M_\odot$ progenitor in full unconstrained 3D and in octant symmetry for $\gtrsim$$ 380\,\mathrm{ms}$. We find the development of an asymmetric runaway explosion in our unconstrained simulation. We test the resolution dependence of our results and, in agreement with previous work, find that low resolution artificially aids explosion and leads to an earlier runaway expansion of the shock. At low resolution, the octant and full 3D dynamics are qualitatively very similar, but at high resolution, only the full 3D simulation exhibits the onset of explosion.

During the merger of a black hole and a neutron star, baryonic mass can become unbound from the system. Because the ejected material is extremely neutron-rich, the r-process rapidly synthesizes heavy nuclides as the material expands and cools. In this work, we map general relativistic models of black hole-neutron star (BHNS) mergers into a Newtonian smoothed particle hydrodynamics (SPH) code and follow the evolution of the thermodynamics and morphology of the ejecta until the outflows become homologous. We investigate how the subsequent evolution depends on our mapping procedure and find that the results are robust. Using thermodynamic histories from the SPH particles, we then calculate the expected nucleosynthesis in these outflows while varying the level of neutrino irradiation coming from the postmerger accretion disk. We find that the ejected material robustly produces r-process nucleosynthesis even for unrealistically high neutrino luminosities, due to the rapid velocities of the outflow. Nonetheless, we find that neutrinos can have an impact on the detailed pat- tern of the r-process nucleosynthesis. Electron neutrinos are captured by neutrons to produce protons while neutron capture is occurring. The produced protons rapidly form low mass seed nuclei for the r-process. These low mass seeds are eventually incorporated into the first r-process peak at A~78, producing mainly Ge and Se. We consider the mechanism of this process in detail and discuss if it can impact galactic chemical evolution of the first peak r-process nuclei.

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

The nearby Sun-like star GJ 758 hosts a cold substellar companion, GJ 758 B, at a projected separation of $\lesssim$30 AU, previously detected in high-contrast multi-band photometric observations. In order to better constrain the companion's physical characteristics, we acquired the first low-resolution ($R \sim 50$) near-infrared spectrum of it using the high-contrast hyperspectral imaging instrument Project 1640 on Palomar Observatory's 5-m Hale telescope. We obtained simultaneous images in 32 wavelength channels covering the $Y$, $J$, and $H$ bands ($\sim$952-1770 nm), and used data processing techniques based on principal component analysis to efficiently subtract chromatic background speckle-noise. GJ 758 B was detected in four epochs during 2013 and 2014. Basic astrometric measurements confirm its apparent northwest trajectory relative to the primary star, with no clear signs of orbital curvature. Spectra of SpeX/IRTF observed T dwarfs were compared to the combined spectrum of GJ 758 B, with ${\chi}^2$ minimization suggesting a best fit for spectral type T7.0$\pm$1.0, but with a shallow minimum over T5-T8. Fitting of synthetic spectra from the BT-Settl13 model atmospheres gives an effective temperature $T_{\text{eff}}=741 \pm 25$ K and surface gravity $\log g = 4.3 \pm 0.5$ dex (cgs). Our derived best-fit spectral type and effective temperature from modeling of the low-resolution spectrum suggest a slightly earlier and hotter companion than previous findings from photometric data, but do not rule out current results, and confirm GJ 758 B as one of the coolest sub-stellar companions to a Sun-like star to date.

We explore heavy element nucleosynthesis in neutrino-driven winds from rapidly-rotating, strongly magnetized proto-neutron stars for which the magnetic dipole is aligned with the rotation axis, and the field is assumed to be a static force-free configuration. We process the proto-magnetar wind trajectories calculated by Vlasov et al 2014 through the r-process nuclear reaction network SkyNet using contemporary models for the evolution of the wind electron fraction during the proto-neutron star cooling phase. Although we do not find a successful second or third peak r-process for any rotation period P, we show that proto-magnetars with P around 1-5 ms produce heavy element abundance distributions that extend to higher nuclear mass number than from otherwise equivalent spherical winds (with the mass fractions of some elements enhanced by factors of 100-1000). The heaviest elements are synthesized by outflows emerging along flux tubes which graze the closed zone and pass near the equatorial plane outside the light cylinder. Due to dependence of the nucleosynthesis pattern on the magnetic field strength and rotation rate of the proto-neutron star, natural variations in these quantities between core collapse events could contribute to the observed diversity of the abundances of weak r-process nuclei in metal-poor stars. Further diversity, including possibly even a successful third-peak r-process, could be achieved for misaligned rotators with non-zero magnetic inclination with respect to the rotation axis. If proto-magnetars are central engines for GRBs, their relativistic jets should contain a high mass fraction of heavy nuclei of characteristic average mass number A of order 100, providing a possible source for ultra-high energy cosmic rays comprised of heavy nuclei with an energy spectrum that extends beyond the nominal GZK cut-off for protons or iron nuclei.

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

Binary neutron star mergers are promising sources of gravitational waves for ground-based detectors such as Advanced LIGO. Neutron-rich material ejected by these mergers may also be the main source of r-process elements in the Universe, while radioactive decays in the ejecta can power bright electromagnetic post-merger signals. Neutrino-matter interactions play a critical role in the evolution of the composition of the ejected material, which significantly impacts the outcome of nucleosynthesis and the properties of the associated electromagnetic signal. In this work, we present a simulation of a binary neutron star merger using an improved method for estimating the average neutrino energies in our energy-integrated neutrino transport scheme. These energy estimates are obtained by evolving the neutrino number density in addition to the neutrino energy and flux densities. We show that significant changes are observed in the composition of the polar ejecta when comparing our new results with earlier simulations in which the neutrino spectrum was assumed to be the same everywhere in optically thin regions. In particular, we find that material ejected in the polar regions is less neutron rich than previously estimated. Our new estimates of the composition of the polar ejecta make it more likely that the color and timescale of the electromagnetic signal depend on the orientation of the binary with respect to an observer's line-of-sight. These results also indicate that important observable properties of neutron star mergers are sensitive to the neutrino energy spectrum, and may need to be studied through simulations including a more accurate, energy-dependent neutrino transport scheme.

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

We present results from the first large parameter study of neutron star mergers using fully general relativistic simulations with finite-temperature microphysical equations of state and neutrino cooling. We consider equal and unequal-mass binaries drawn from the galactic population and simulate each binary with three different equations of state. Our focus is on the emission of energy and angular momentum in gravitational waves in the postmerger phase. We find that the emitted gravitational-wave energy in the first $\sim$$10\,\mathrm{ms}$ of the life of the resulting hypermassive neutron star (HMNS) is about twice the energy emitted over the entire inspiral history of the binary. The total radiated energy per binary mass is comparable to or larger than that of nonspinning black hole inspiral-mergers. About $0.8-2.5\%$ of the binary mass-energy is emitted at kHz frequencies in the early HMNS evolution. We find a clear dependence of the postmerger GW emission on binary configuration and equation of state and show that it can be encoded as a broad function of the binary tidal coupling constant $\kappa^T_2$. Our results also demonstrate that the dimensionless spin of black holes resulting from subsequent HMNS collapse are limited to $\lesssim0.7-0.8$. This may significantly impact the neutrino pair annihilation mechanism for powering short gamma-ray bursts (sGRB).

We present fully general-relativistic simulations of binary neutron star mergers with a temperature and composition dependent nuclear equation of state. We study the dynamical mass ejection from both quasi-circular and dynamical-capture eccentric mergers. We systematically vary the level of our treatment of the microphysics to isolate the effects of neutrino cooling and heating and we compute the nucleosynthetic yields of the ejecta. We find that eccentric binaries can eject significantly more material than quasi-circular binaries and generate bright infrared and radio emission. In all our simulations the outflow is composed of a combination of tidally- and shock-driven ejecta, mostly distributed over a broad $\sim 60^\circ$ angle from the orbital plane, and, to a lesser extent, by thermally driven winds at high latitudes. Ejecta from eccentric mergers are typically more neutron rich than those of quasi-circular mergers. We find neutrino cooling and heating to affect, quantitatively and qualitatively, composition, morphology, and total mass of the outflows. This is also reflected in the infrared and radio signatures of the binary. The final nucleosynthetic yields of the ejecta are robust and insensitive to input physics or merger type in the regions of the second and third r-process peaks. The yields for elements on the first peak vary between our simulations, but none of our models is able to explain the Solar abundances of first-peak elements without invoking additional first-peak contributions from either neutrino and viscously-driven winds operating on longer timescales after the mergers, or from core-collapse supernovae.

$\mu$ Her is a nearby quadruple system with a G-subgiant primary and several low mass companions arranged in a 2+2 architecture. While the BC components have been well characterized, the Ab component has been detected astrometrically and with direct imaging but there has been some confusion over its nature, in particular whether the companion is stellar or substellar. Using near-infrared spectroscopy we are able to estimate the spectral type of the companion as a M4$\pm$1V star. In addition, we have measured the astrometry of the system for over a decade. We combined the astrometry with archival radial velocity measurements to compute an orbit of the system. From the combined orbit, we are able to compute the mass sum of the system. Using the estimated mass of the primary, we estimate the mass of the secondary as 0.32 M_sun, which agrees with the estimated spectral type. Our computed orbit is preliminary due to the incomplete orbital phase coverage, but it should be sufficient to predict ephemerides over the next decade.

Neutron star mergers are among the most promising sources of gravitational waves for advanced ground-based detectors. These mergers are also expected to power bright electromagnetic signals, in the form of short gamma-ray bursts, infrared/optical transients, and radio emission. Simulations of these mergers with fully general relativistic codes are critical to understand the merger and post-merger gravitational wave signals and their neutrinos and electromagnetic counterparts. In this paper, we employ the SpEC code to simulate the merger of low-mass neutron star binaries (two $1.2M_\odot$ neutron stars) for a set of three nuclear-theory based, finite temperature equations of state. We show that the frequency peaks of the post-merger gravitational wave signal are in good agreement with predictions obtained from simulations using a simpler treatment of gravity. We find, however, that only the fundamental mode of the remnant is excited for long periods of time: emission at the secondary peaks is damped on a millisecond timescale in the simulated binaries. For such low-mass systems, the remnant is a massive neutron star which, depending on the equation of state, is either permanently stable or long-lived. We observe strong excitations of l=2, m=2 modes, both in the massive neutron star and in the form of hot, shocked tidal arms in the surrounding accretion torus. We estimate the neutrino emission of the remnant using a neutrino leakage scheme and, in one case, compare these results with a gray two-moment neutrino transport scheme. We confirm the complex geometry of the neutrino emission, also observed in previous simulations with neutrino leakage, and show explicitly the presence of important differences in the neutrino luminosity, disk composition, and outflow properties between the neutrino leakage and transport schemes.

r-Process nucleosynthesis in material ejected during neutron star mergers may lead to radioactively powered transients called kilonovae. The timescale and peak luminosity of these transients depend on the composition of the ejecta, which determines the local heating rate from nuclear decays and the opacity. Kasen et al. (2013, ApJ, 774, 25) and Tanaka & Hotokezaka (2013, ApJ, 775, 113) pointed out that lanthanides can drastically increase the opacity in these outflows. We use the new general-purpose nuclear reaction network SkyNet to carry out a parameter study of r-process nucleosynthesis for a range of initial electron fractions $Y_e$, initial specific entropies $s$, and expansion timescales $\tau$. We find that the ejecta is lanthanide-free for $Y_e \gtrsim 0.22 - 0.30$, depending on $s$ and $\tau$. The heating rate is insensitive to $s$ and $\tau$, but certain, larger values of $Y_e$ lead to reduced heating rates, due to individual nuclides dominating the heating. We calculate approximate light curves with a simplified gray radiative transport scheme. The light curves peak at about a day (week) in the lanthanide-free (-rich) cases. The heating rate does not change much as the ejecta becomes lanthanide-free with increasing $Y_e$, but the light curve peak becomes about an order of magnitude brighter because it peaks much earlier when the heating rate is larger. We also provide parametric fits for the heating rates between 0.1 and $100\,\text{days}$, and we provide a simple fit in $Y_e$, $s$, and $\tau$ to estimate whether the ejecta is lanthanide-rich or not.

We present the SuperNova Explosion Code (SNEC), an open-source Lagrangian code for the hydrodynamics and equilibrium-diffusion radiation transport in the expanding envelopes of supernovae. Given a model of a progenitor star, an explosion energy, and an amount and distribution of radioactive nickel, SNEC generates the bolometric light curve, as well as the light curves in different broad bands assuming black body emission. As a first application of SNEC, we consider the explosions of a grid of 15 Msun (at zero-age main sequence) stars whose hydrogen envelopes are stripped to different extents and at different points in their evolution. The resulting light curves exhibit plateaus with durations of ~20-100 days if >~1.5-2 Msun of hydrogen-rich material is left and no plateau if less hydrogen-rich material is left. If these shorter plateau lengths are not seen for Type IIP supernovae in nature, it suggests that, at least for zero-age main sequence masses <~ 20 Msun, hydrogen mass loss occurs as an all or nothing process. This perhaps points to the important role binary interactions play in generating the observed mass-stripped supernovae (i.e., Type Ib/c events). These light curves are also unlike what is typically seen for Type IIL supernovae, arguing that simply varying the amount of mass loss cannot explain these events. The most stripped models begin to show double-peaked light curves similar to what is often seen for Type IIb supernovae, confirming previous work that these supernovae can come from progenitors that have a small amount of hydrogen and a radius of ~500 Rsun.

HD 177830 is an evolved K0IV star with two known exoplanets. In addition to the planetary companions it has a late-type stellar companion discovered with adaptive optics imagery. We observed the binary star system with the PHARO near-IR camera and the Project 1640 coronagraph. Using the Project 1640 coronagraph and integral field spectrograph we extracted a spectrum of the stellar companion. This allowed us to determine that the spectral type of the stellar companion is a M4$\pm$1V. We used both instruments to measure the astrometry of the binary system. Combining these data with published data, we determined that the binary star has a likely period of approximately 800 years with a semi-major axis of 100-200 AU. This implies that the stellar companion has had little or no impact on the dynamics of the exoplanets. The astrometry of the system should continue to be monitored, but due to the slow nature of the system, observations can be made once every 5-10 years.

We conduct a series of numerical experiments into the nature of three-dimensional (3D) hydrodynamics in the postbounce stalled-shock phase of core-collapse supernovae using 3D general-relativistic hydrodynamic simulations of a $27$-$M_\odot$ progenitor star with a neutrino leakage/heating scheme. We vary the strength of neutrino heating and find three cases of 3D dynamics: (1) neutrino-driven convection, (2) initially neutrino-driven convection and subsequent development of the standing accretion shock instability (SASI), (3) SASI dominated evolution. This confirms previous 3D results of Hanke et al. 2013, ApJ 770, 66 and Couch & Connor 2014, ApJ 785, 123. We carry out simulations with resolutions differing by up to a factor of $\sim$4 and demonstrate that low resolution is artificially favorable for explosion in the 3D convection-dominated case, since it decreases the efficiency of energy transport to small scales. Low resolution results in higher radial convective fluxes of energy and enthalpy, more fully buoyant mass, and stronger neutrino heating. In the SASI-dominated case, lower resolution damps SASI oscillations. In the convection-dominated case, a quasi-stationary angular kinetic energy spectrum $E(\ell)$ develops in the heating layer. Like other 3D studies, we find $E(\ell) \propto \ell^{-1}$ in the "inertial range," while theory and local simulations argue for $E(\ell) \propto \ell^{-5/3}$. We argue that current 3D simulations do not resolve the inertial range of turbulence and are affected by numerical viscosity up to the energy containing scale, creating a "bottleneck" that prevents an efficient turbulent cascade.

We present a new and computationally efficient method for characterizing very low mass companions using low resolution ($R\sim$30) near-infrared ($YJH$) spectra from high contrast imaging campaigns with integral field spectrograph (IFS) units. We conduct a detailed quantitative comparison of the efficacy of this method through tests on simulated data comparable in spectral coverage and resolution to the currently operating direct imaging systems around the world. In particular, we simulate Project 1640 data as an example of the use, accuracy, and precision of this technique. We present results from comparing simulated spectra of M, L, and T dwarfs with a large and finely-sampled grid of synthetic spectra using Markov Chain Monte Carlo techniques. We determine the precision and accuracy of effective temperature and surface gravity inferred from fits to PHOENIX dusty and cond, which we find reproduce the low-resolution spectra of all objects within the adopted flux uncertainties. Uncertainties in effective temperature decrease from $\pm$100-500 K for M dwarfs to as small as $\pm$30 K for some L and T spectral types. Surface gravity is constrained to within 0.2-0.4 dex for mid-L through T dwarfs, but uncertainties are as large as 1.0 dex or more for M dwarfs. Results for effective temperature from low-resolution $YJH$ spectra generally match predictions from published spectral type-temperature relationships except for L-T transition objects and young objects. Single-band spectra (i.e., narrower wavelength coverage) result in larger uncertainties and often discrepant results, suggesting that high contrast IFS observing campaigns can compensate for low spectral resolution by expanding the wavelength coverage for reliable characterization of detected companions. [Abstract truncated]

Investigations of elemental abundances in the ancient and most metal deficient stars are extremely important because they serve as tests of variable nucleosynthesis pathways and can provide critical inferences of the type of stars that lived and died before them. The presence of r-process elements in a handful of carbon-enhanced metal-poor (CEMP-r) stars, which are assumed to be closely connected to the chemical yield from the first stars, is hard to reconcile with standard neutron star mergers. Here we show that the production rate of dynamically assembled compact binaries in high-z nuclear star clusters can attain a sufficient high value to be a potential viable source of heavy r-process material in CEMP-r stars. The predicted frequency of such events in the early Galaxy, much lower than the frequency of Type II supernovae but with significantly higher mass ejected per event, can naturally lead to a high level of scatter of Eu as observed in CEMP-r stars.

We present a first simulation of the post-merger evolution of a black hole-neutron star binary in full general relativity using an energy-integrated general relativistic truncated moment formalism for neutrino transport. We describe our implementation of the moment formalism and important tests of our code, before studying the formation phase of a disk after a black hole-neutron star merger. We use as initial data an existing general relativistic simulation of the merger of a neutron star of 1.4 solar mass with a black hole of 7 solar mass and dimensionless spin a/M=0.8. Comparing with a simpler leakage scheme for the treatment of the neutrinos, we find noticeable differences in the neutron to proton ratio in and around the disk, and in the neutrino luminosity. We find that the electron neutrino luminosity is much lower in the transport simulations, and that the remnant is less neutron-rich. The spatial distribution of the neutrinos is significantly affected by relativistic effects. Over the short timescale evolved, we do not observe purely neutrino-driven outflows. However, a small amount of material (3e-4Msun) is ejected in the polar region during the circularization of the disk. Most of that material is ejected early in the formation of the disk, and is fairly neutron rich. Through r-process nucleosynthesis, that material should produce high-opacity lanthanides in the polar region, and could thus affect the lightcurve of radioactively powered electromagnetic transients. We also show that by the end of the simulation, while the bulk of the disk is neutron-rich, its outer layers have a higher electron fraction. As that material would be the first to be unbound by disk outflows on longer timescales, the changes in Ye experienced during the formation of the disk could have an impact on the nucleosynthesis outputs from neutrino-driven and viscously-driven outflows. [Abridged]

We present an analysis of the orbital motion of the four sub-stellar objects orbiting HR8799. Our study relies on the published astrometric history of this system augmented with an epoch obtained with the Project 1640 coronagraph + Integral Field Spectrograph (IFS) installed at the Palomar Hale telescope. We first focus on the intricacies associated with astrometric estimation using the combination of an Extreme Adaptive Optics system (PALM-3000), a coronagraph and an IFS. We introduce two new algorithms. The first one retrieves the stellar focal plane position when the star is occulted by a coronagraphic stop. The second one yields precise astrometric and spectro-photometric estimates of faint point sources even when they are initially buried in the speckle noise. The second part of our paper is devoted to studying orbital motion in this system. In order to complement the orbital architectures discussed in the literature, we determine an ensemble of likely Keplerian orbits for HR8799bcde, using a Bayesian analysis with maximally vague priors regarding the overall configuration of the system. While the astrometric history is currently too scarce to formally rule out coplanarity, HR8799d appears to be misaligned with respect to the most likely planes of HR8799bce orbits. This misalignment is sufficient to question the strictly coplanar assumption made by various authors when identifying a Laplace resonance as a potential architecture. Finally, we establish a high likelihood that HR8799de have dynamical masses below 13 M_Jup using a loose dynamical survival argument based on geometric close encounters. We illustrate how future dynamical analyses will further constrain dynamical masses in the entire system.

The merger of two neutron stars leaves behind a rapidly spinning hypermassive object whose survival is believed to depend on the maximum mass supported by the nuclear equation of state, angular momentum redistribution by (magneto-)rotational instabilities, and spindown by gravitational waves. The high temperatures (~5-40 MeV) prevailing in the merger remnant may provide thermal pressure support that could increase its maximum mass and, thus, its life on a neutrino-cooling timescale. We investigate the role of thermal pressure support in hypermassive merger remnants by computing sequences of spherically-symmetric and axisymmetric uniformly and differentially rotating equilibrium solutions to the general-relativistic stellar structure equations. Using a set of finite-temperature nuclear equations of state, we find that hot maximum-mass critically spinning configurations generally do not support larger baryonic masses than their cold counterparts. However, subcritically spinning configurations with mean density of less than a few times nuclear saturation density yield a significantly thermally enhanced mass. Even without decreasing the maximum mass, cooling and other forms of energy loss can drive the remnant to an unstable state. We infer secular instability by identifying approximate energy turning points in equilibrium sequences of constant baryonic mass parametrized by maximum density. Energy loss carries the remnant along the direction of decreasing gravitational mass and higher density until instability triggers collapse. Since configurations with more thermal pressure support are less compact and thus begin their evolution at a lower maximum density, they remain stable for longer periods after merger.

The Debris Disk Explorer (DDX) is a proposed balloon-borne investigation of debris disks around nearby stars. Debris disks are analogs of the Asteroid Belt (mainly rocky) and Kuiper Belt (mainly icy) in our Solar System. DDX will measure the size, shape, brightness, and color of tens of disks. These measurements will enable us to place the Solar System in context. By imaging debris disks around nearby stars, DDX will reveal the presence of perturbing planets via their influence on disk structure, and explore the physics and history of debris disks by characterizing the size and composition of disk dust. The DDX instrument is a 0.75-m diameter off-axis telescope and a coronagraph carried by a stratospheric balloon. DDX will take high-resolution, multi-wavelength images of the debris disks around tens of nearby stars. Two flights are planned; an overnight test flight within the United States followed by a month-long science flight launched from New Zealand. The long flight will fully explore the set of known debris disks accessible only to DDX. It will achieve a raw contrast of 10^-7, with a processed contrast of 10^-8. A technology benefit of DDX is that operation in the near-space environment will raise the Technology Readiness Level of internal coronagraphs, deformable mirrors, and wavefront sensing and control, all potentially needed for a future space-based telescope for high-contrast exoplanet imaging.

Kappa Andromedae is a B9IVn star at 52 pc for which a faint substellar companion separated by 55 AU was recently announced. In this work, we present the first spectrum of the companion, "kappa And B," using the Project 1640 high-contrast imaging platform. Comparison of our low-resolution YJH-band spectra to empirical brown dwarf spectra suggests an early-L spectral type. Fitting synthetic spectra from PHOENIX model atmospheres to our observed spectrum allows us to constrain the effective temperature to ~2000K, as well as place constraints on the companion surface gravity. Further, we use previously reported log(g) and effective temperature measurements of the host star to argue that the kappa And system has an isochronal age of 220 +/- 100 Myr, older than the 30 Myr age reported previously. This interpretation of an older age is corroborated by the photometric properties of kappa And B, which appear to be marginally inconsistent with other 10-100 Myr low-gravity L-dwarfs for the spectral type range we derive. In addition, we use Keck aperture masking interferometry combined with published radial velocity measurements to rule out the existence of any tight stellar companions to kappa And A that might be responsible for the system's overluminosity. Further, we show that luminosity enhancements due to a nearly "pole-on" viewing angle coupled with extremely rapid rotation is unlikely. Kappa And A is thus consistent with its slightly evolved luminosity class (IV) and we propose here that kappa And, with a revised age of 220 +/- 100 Myr, is an interloper to the 30 Myr Columba association with which it was previously associated. The photometric and spectroscopic evidence for kappa And B combined with our re-assesment of the system age implies a substellar companion mass of 50^{+16}_{-13} Jupiter Masses, consistent with a brown dwarf rather than a planetary mass companion.

We present adaptive optics photometry and spectra in the JHKL-bands along with high spectral resolution K-band spectroscopy for each component of the Z Canis Majoris system. Our high angular resolution photometry of this very young (<1 Myr) binary, comprised of an FU Ori object and a Herbig Ae/Be star, were gathered shortly after the 2008 outburst while our high resolution spectroscopy was gathered during a quiescent phase. Our photometry conclusively determine that the outburst was due solely to the embedded Herbig Ae/Be member, supporting results from earlier works, and that the optically visible FU Ori component decreased slightly (~30%) in luminosity during the same period, consistent with previous works on the variability of FU Ori type systems. Further, our high-resolution K-band spectra definitively demonstrate that the 2.294 micron CO absorption feature seen in composite spectra of the system is due solely to the FU Ori component, while a prominent CO emission feature at the same wavelength, long suspected to be associated with the innermost regions of a circumstellar accretion disk, can be assigned to the Herbig Ae/Be member. These findings are in contrast to previous analyses (e.g. Malbet et al 2010, Benisty et al. 2010) of this complex system which assigned the CO emission to the FU Ori component.

We report low resolution near infrared spectroscopic observations of the eruptive star FU Orionis using the Integral Field Spectrograph Project 1640 installed at the Palomar Hale telescope. This work focuses on elucidating the nature of the faint source, located 0.5" south of FU Ori, and identified in 2003 as FU Ori S. We first use our observations in conjunction with published data to demonstrate that the two stars are indeed physically associated and form a true binary pair. We then proceed to extract J and H band spectro-photometry using the damped LOCI algorithm, a reduction method tailored for high contrast science with IFS. This is the first communication reporting the high accuracy of this technique, pioneered by the Project 1640 team, on a faint astronomical source. We use our low resolution near infrared spectrum in conjunction with 10.2 micron interferometric data to constrain the infrared excess of FU Ori S. We then focus on estimating the bulk physical properties of FU Ori S. Our models lead to estimates of an object heavily reddened, A_V =8-12, with an effective temperature of ~ 4000-6500 K . Finally we put these results in the context of the FU Ori N-S system and argue that our analysis provides evidence that FU Ori S might be the more massive component of this binary system

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

A new code for following the evolution and emissions of proto-neutron stars during the first minute of their lives is developed and tested. The code is one dimensional, fully implicit, and general relativistic. Multi-group, multi-flavor neutrino transport is incorporated that makes use of variable Eddington factors obtained from a formal solution of the static general relativistic Boltzmann equation with linearized scattering terms. The timescales of neutrino emission and spectral evolution obtained using the new code are broadly consistent with previous results. Unlike other recent calculations, however, the new code predicts that the neutrino-driven wind will be characterized, at least for part of its existence, by a neutron excess. This change, potentially consequential for nucleosynthesis in the wind, is due to an improved treatment of the charged-current interactions of electron flavored neutrinos and anti-neutrinos with nucleons. A comparison is also made between the results obtained using either variable Eddington factors or simple equilibrium flux-limited diffusion. The latter approximation, which has been frequently used in previous studies of proto-neutron star cooling, accurately describes the total neutrino luminosities (to within 10%) for most of the evolution, until the proto-neutron star becomes optically thin.

We used the Project 1640 near-infrared coronagraph and integral field spectrograph to observe 19 young solar type stars. Five of these stars are known binary stars and we detected the late-type secondaries and were able to measure their JH spectra with a resolution of R\sim30. The reduced, extracted, and calibrated spectra were compared to template spectra from the IRTF spectral library. With this comparison we test the accuracy and consistency of spectral type determination with the low-resolution near-infrared spectra from P1640. Additionally, we determine effective temperature and surface gravity of the companions by fitting synthetic spectra calculated with the PHOENIX model atmosphere code. We also present several new epochs of astrometry of each of the systems. Together these data increase our knowledge and understanding of the stellar make up of these systems. In addition to the astronomical results, the analysis presented helps validate the Project 1640 data reduction and spectral extraction processes and the utility of low-resolution, near-infrared spectra for characterizing late-type companions in multiple systems.

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

High-contrast imaging instruments are now being equipped with integral field spectrographs (IFS) to facilitate the detection and characterization of faint substellar companions. Algorithms currently envisioned to handle IFS data, such as the Locally Optimized Combination of Images (LOCI) algorithm, rely upon aggressive point-spread-function (PSF) subtraction, which is ideal for initially identifying companions but results in significantly biased photometry and spectroscopy due to unwanted mixing with residual starlight. This spectro-photometric issue is further complicated by the fact that algorithmic color response is a function of the companion's spectrum, making it difficult to calibrate the effects of the reduction without using iterations involving a series of injected synthetic companions. In this paper, we introduce a new PSF calibration method, which we call "damped LOCI", that seeks to alleviate these concerns. By modifying the cost function that determines the weighting coefficients used to construct PSF reference images, and also forcing those coefficients to be positive, it is possible to extract companion spectra with a precision that is set by calibration of the instrument response and transmission of the atmosphere, and not by post-processing. We demonstrate the utility of this approach using on-sky data obtained with the Project 1640 IFS at Palomar. Damped-LOCI does not require any iterations on the underlying spectral type of the companion, nor does it rely upon priors involving the chromatic and statistical properties of speckles. It is a general technique that can readily be applied to other current and planned instruments that employ IFS's.

A study of the host stars to exoplanets is important to understanding their environment. To that end, we report new speckle observations of a sample of exoplanet host primaries. The bright exoplanet host HD 8673 (= HIP 6702) is revealed to have a companion, although at this time we cannot definitively establish the companion as physical or optical. The observing lists for planet searches and for these observations have for the most part been pre-screened for known duplicity, so the detected binary fraction is lower than what would otherwise be expected. Therefore, a large number of double stars were observed contemporaneously for verification and quality control purposes, to ensure the lack of detection of companions for exoplanet hosts was valid. In these additional observations, ten pairs are resolved for the first time and sixty pairs are confirmed. These observations were obtained with the USNO speckle camera on the NOAO 4m telescopes at both KPNO and CTIO from 2001 to 2010.

The results of an adaptive optics survey of exoplanet host stars for stellar companions is presented. We used the AEOS Telescope and its adaptive optics system to collect deep images of the stars in $I$-band. Sixty-two exoplanet host stars were observed and fifteen multiple star systems were resolved. Of these eight are known multiples, while seven are new candidate binaries. For all binaries, we measured the relative astrometry of the pair and the differential magnitude in $I$-band. We improved the orbits of HD 19994 and $\tau$ Boo. These observations will provide improved statistics on the duplicity of exoplanet hosts stars and provide an increased understanding of the dynamics of known binary star exoplanet hosts.

The possibility that long tidal tails formed during compact object mergers may power optical transients through the decay of freshly synthesized r-process material is investigated. Precise modeling of the merger dynamics allows for a realistic determination of the thermodynamic conditions in the ejected debris. The results of hydrodynamic and full nuclear network calculations are combined to calculate the resultant r-process abundances and the heating of the material by their decays. The subsequent homologous structure is mapped into a radiative transfer code to synthesize emergent model light curves and determine how their properties (variability and color evolution) depend on the mass ratio and orientation of the merging binary. The radiation emanating from the ejected debris, though less spectacular than a typical supernova, should be observable in transient surveys and we estimate the associated detection rates. The case for (or against) compact object mergers as the progenitors of short gamma-ray bursts can be tested if such electromagnetic transients are detected (or not) in coincidence with some bursts, although they may be obscured by on-axis afterglows.

Project 1640 is a high-contrast imaging instrument recently commissioned at Palomar observatory. A combination of a coronagraph with an integral field spectrograph (IFS), Project 1640 is designed to detect and characterize extrasolar planets, brown dwarfs, and circumstellar material orbiting nearby stars. In this paper, we present our data processing techniques for improving upon instrument raw sensitivity via the removal of quasi-static speckles. Our approach utilizes the chromatic image diversity provided by the IFS in combination with the locally-optimized combination of images (LOCI) algorithm to suppress the intensity of residual contaminating light in close angular proximity to target stars. We describe the Project 1640 speckle suppression pipeline (PSSP) and demonstrate the ability to detect companions with brightness comparable to and below that of initial speckle intensities using on-sky commissioning data. Our preliminary results indicate that suppression factors of at least one order of magnitude are consistently possible, reaching $5\sigma$ contrast levels of $2.1\times10^{-5}$ at $1\arcsec$ in the H-band in 20 minutes of on-source integration time when non-common-path errors are reasonably well-calibrated. These results suggest that near-infrared contrast levels of order $\approx10^{-7}$ at subarcsecond separations will soon be possible for Project 1640 and similarly designed instruments that receive a diffraction-limited beam corrected by adaptive optics (AO) systems employing deformable mirrors with high actuator-density.

The electron-to-nucleon ratio or electron fraction is a key parameter in many astrophysical studies. Its value is determined by weak-interaction rates that are based on theoretical calculations subject to several nuclear physics uncertainties. Consequently, it is important to have a model independent way of constraining the electron fraction value in different astrophysical environments. Here we show that nuclear statistical equilibrium combined with beta equilibrium can provide such a constraint. We test the validity of this approximation in presupernova models and give lower limits for the electron fraction in type Ia supernova and accretion-induced collapse.

The nearby star Alpha Oph (Ras Alhague) is a rapidly rotating A5IV star spinning at ~89% of its breakup velocity. This system has been imaged extensively by interferometric techniques, giving a precise geometric model of the star's oblateness and the resulting temperature variation on the stellar surface. Fortuitously, Alpha Oph has a previously known stellar companion, and characterization of the orbit provides an independent, dynamically-based check of both the host star and the companion mass. Such measurements are crucial to constrain models of such rapidly rotating stars. In this study, we combine eight years of Adaptive Optics imaging data from the Palomar, AEOS, and CFHT telescopes to derive an improved, astrometric characterization of the companion orbit. We also use photometry from these observations to derive a model-based estimate of the companion mass. A fit was performed on the photocenter motion of this system to extract a component mass ratio. We find masses of 2.40^{0.23}_{0.37} solar masses and 0.85^{0.06}_{0.04} solar masses for Alpha Oph A and Alpha Oph B, respectively. Previous orbital studies of this system found a mass too high for this system, inconsistent with stellar evolutionary calculations. Our measurements of the host star mass are more consistent with these evolutionary calculations, but with slightly higher uncertainties. In addition to the dynamically-derived masses, we use IJHK photometry to derive a model-based mass for Alpha Oph B, of 0.77 +/- 0.05 solar masses marginally consistent with the dynamical masses derived from our orbit. Our model fits predict a periastron passage on 2012 April 19, with the two components having a ~50 milliarcsec separation from March to May 2012. A modest amount of interferometric and radial velocity data during this period could provide a mass determination of this star at the few percent level.

The Lyot project used an optimized Lyot coronagraph with Extreme Adaptive Optics at the 3.63m Advanced Electro-Optical System telescope (AEOS) to observe 86 stars from 2004 to 2007. In this paper we give an overview of the survey results and a statistical analysis of the observed nondetections around 58 of our targets to place constraints on the population of substellar companions to nearby stars. The observations did not detect any companion in the substellar regime. Since null results can be as important as detections, we analyzed each observation to determine the characteristics of the companions that can be ruled out. For this purpose we use a Monte Carlo approach to produce artificial companions, and determine their detectability by comparison with the sensitivity curve for each star. All the non-detection results are combined using a Bayesian approach and we provide upper limits on the population of giant exoplanets and brown dwarfs for this sample of stars. Our nondetections confirm the rarity of brown dwarfs around solar-like stars and we constrain the frequency of massive substellar companions (M>40Mjup) at orbital separation between and 10 and 50 AU to be <20%.

Although they are but a small fraction of the mass ejected in core-collapse supernovae, neutrino-driven winds (NDWs) from nascent proto-neutron stars (PNSs) have the potential to contribute significantly to supernova nucleosynthesis. In previous works, the NDW has been implicated as a possible source of r-process and light p-process isotopes. In this paper we present time-dependent hydrodynamic calculations of nucleosynthesis in the NDW which include accurate weak interaction physics coupled to a full nuclear reaction network. Using two published models of PNS neutrino luminosities, we predict the contribution of the NDW to the integrated nucleosynthetic yield of the entire supernova. For the neutrino luminosity histories considered, no true r-process occurs in the most basic scenario. The wind driven from an older $1.4 M_\odot$ model for a PNS is moderately neutron-rich at late times however, and produces $^{87}$Rb, $^{88}$Sr, $^{89}$Y, and $^{90}$Zr in near solar proportions relative to oxygen. The wind from a more recently studied $1.27 M_\odot$ PNS is proton-rich throughout its entire evolution and does not contribute significantly to the abundance of any element. It thus seems very unlikely that the simplest model of the NDW can produce the r-process. At most, it contributes to the production of the N = 50 closed shell elements and some light p-nuclei. In doing so, it may have left a distinctive signature on the abundances in metal poor stars, but the results are sensitive to both uncertain models for the explosion and the masses of the neutron stars involved.

Through the combination of high-order Adaptive Optics and coronagraphy, we report the discovery of a faint stellar companion to the A3V star zeta Virginis. This companion is ~7 magnitudes fainter than its host star in the H-band, and infrared imaging spanning 4.75 years over five epochs indicates this companion has common proper motion with its host star. Using evolutionary models, we estimate its mass to be 0.168+/-.016 solar masses, giving a mass ratio for this system q = 0.082. Assuming the two objects are coeval, this mass suggests a M4V-M7V spectral type for the companion, which is confirmed through integral field spectroscopic measurements. We see clear evidence for orbital motion from this companion and are able to constrain the semi-major axis to be greater than 24.9 AU, the period > 124$ yrs, and eccentricity > 0.16. Multiplicity studies of higher mass stars are relatively rare, and binary companions such as this one at the extreme low end of the mass ratio distribution are useful additions to surveys incomplete at such a low mass ratio. Moreover, the frequency of binary companions can help to discriminate between binary formation scenarios that predict an abundance of low-mass companions forming from the early fragmentation of a massive circumstellar disk. A system such as this may provide insight into the anomalous X-ray emission from A stars, hypothesized to be from unseen late-type stellar companions. Indeed, we calculate that the presence of this M-dwarf companion easily accounts for the X-ray emission from this star detected by ROSAT.

The A5V star Alcor has an M3-M4 dwarf companion, as evidenced by a novel astrometric technique. Imaging spectroscopy combined with adaptive optics coronagraphy allowed for the detection and spectrophotometric characterization of the point source at a contrast of ~6 J- and H-band magnitudes and separation of 1" from the primary star. The use of an astrometric pupil plane grid allowed us to determine the projected separations between the companion and the coronagraphically occulted primary star to <=3 milliarcsecond precision at two observation epochs. Our measurements demonstrate common parallactic and proper motion over the course of 103 days, significantly shorter than the period of time needed for most companion confirmations through proper motion measurements alone. This common parallax method is potentially more rigorous than common proper motion, ensuring that the neighboring bodies lie at the same distance, rather than relying on the statistical improbability that two objects in close proximity to each other on the sky move in the same direction. The discovery of a low-mass (~0.25M_sun) companion around a bright (V = 4.0), nearby (d = 25 pc) star highlights a region of binary star parameter space that to date has not been fully probed.

We demonstrate the versatility of a dual imaging polarimeter working in tandem with a Lyot coronagraph and Adaptive Optics to suppress the highly static speckle noise pattern--the greatest hindrance to ground-based direct imaging of planets and disks around nearby stars. Using a double difference technique with the polarimetric data, we quantify the level of speckle suppression, and hence improved sensitivity, by placing an ensemble of artificial faint companions into real data, with given total brightness and polarization. For highly polarized sources within 0.5 arcsec, we show that we achieve 3 to 4 magnitudes greater sensitivity through polarimetric speckle suppression than simply using a coronagraph coupled to a high-order Adaptive Optics system. Using such a polarimeter with a classical Lyot coronagraph at the 3.63m AEOS telescope, we have obtained a 6.5 sigma detection in the H-band of the 76 AU diameter circumstellar debris disk around the star HR 4796A. Our data represent the first definitive, ground-based, near-IR polarimetric image of the HR 4796A debris disk and clearly show the two outer ansae of the disk, evident in Hubble Space Telescope NICMOS/STIS imaging. We derive a lower limit to the fractional linear polarization of 29% caused by dust grains in the disk. In addition, we fit simple morphological models of optically thin disks to our data allowing us to constrain the dust disk scale height to 2.5{+5.0}_{-1.3} AU and scattering asymmetry parameter (g=0.20^{+.07}_{-.10}). These values are consistent with several lines of evidence suggesting that the HR 4796A disk is dominated by a micron-sized dust population, and are indeed typical of disks in transition between those surrounding the Herbig Ae stars to those associated with Vega-like stars.

We present the results of an adaptive optics survey for faint companions among Galactic O-type star systems (V < 8) using the Advanced Electro-Optical System (AEOS) 3.6-meter telescope on Haleakala. We surveyed these O star systems in I-band, typically being able to detect a companion with a magnitude difference of about 6 in the projected separation range of 0.5 to 1.0 arcseconds, and of about 9.5 in the range of 1.0 to 5.0 arcseconds. In the course of the survey, we discovered 40 new companions among 31 of the 116 objects examined and made astrometric and differential magnitude measurements of 24 additional known pairs, several of them being confirmation detections. We present new astrometric orbits for two binaries, BU 1032AB (WDS 05387-0236; sig Ori AB) and SEE 322 (WDS 17158-3344; HD 155889AB). We lack magnitude differences for other filter bands, so it is difficult to determine physical from line-of-sight companions, but we present empirical arguments for the limiting magnitude difference where field contamination is significant.

We have investigated the temporal variability and statistics of the "instantaneous" Strehl ratio. The observations were carried out with the 3.63-m AEOS telescope equipped with a high-order adaptive optics system. In this paper Strehl ratio is defined as the peak intensity of a single short exposure. We have also studied the behaviour of the phase variance computed on the reconstructed wavefronts. We tested the Marechal approximation and used it to explain the observed negative skewness of the Strehl ratio distribution. The estimate of the phase variance is shown to fit a three-parameter Gamma distribution model. We show that simple scaling of the reconstructed wavefronts has a large impact on the shape of the Strehl ratio distribution.

The major obstacle to the direct detection of companions to nearby stars is the overwhelming brightness of the host star. Current instruments employing the combination of adaptive optics (AO) and coronagraphy can typically detect objects within 2'' of the star that are 10^{4-5} times fainter. Correlated speckle noise is one of the biggest obstacles limiting such high-contrast imaging. We have obtained a series of 284 8 s, AO-corrected, coronagraphically occulted H-band images of the star Vega at the 3.63 m AEOS telescope located on Haleakala, Hawaii. This dataset is unique for studying the temporal behavior of speckle noise, and represents the first time such a study on highly corrected _coronagraphic_ AO images has been carried out in a quantitative way. We find the speckle pattern to be highly stable in both position and time in our data. This is due to the fact that the AO system corrects disturbances to the stellar wave front at the level where the instrumental wave front errors dominate. Because of this, we find that our detection limit is not significantly improved simply with increased exposure time alone. However, we are able to improve our dynamic range by 1.5-2 magnitudes through subtraction of static/quasi-static speckles in two rotating frames: the telescope pupil frame and the deformable mirror frame. Furthermore, from our data, we are able to constrain the mass of any purported companion to Vega to be less than 45 M_J at 8 AU and less than 30 M_J at 16 AU, radii not previously probed at these sensitivities.