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.

We describe the design and assembly of the LUX-ZEPLIN experiment, a direct detection search for cosmic WIMP dark matter particles. The centerpiece of the experiment is a large liquid xenon time projection chamber sensitive to low energy nuclear recoils. Rejection of backgrounds is enhanced by a Xe skin veto detector and by a liquid scintillator Outer Detector loaded with gadolinium for efficient neutron capture and tagging. LZ is located in the Davis Cavern at the 4850' level of the Sanford Underground Research Facility in Lead, South Dakota, USA. We describe the major subsystems of the experiment and its key design features and requirements.

We present calculated rate coefficients for ro-vibrational transitions of CO in collisions with H atoms for a gas temperature range of 10 K $\leq T \leq$ 3000 K, based on the recent three-dimensional ab initio H-CO interaction potential of Song et al(2013). Rate coefficients for ro-vibrational $v=1,j=0-30 \rightarrow v'=0, j'$ transitions were obtained from scattering cross sections previously computed with the close-coupling method by Song et al(2015). Combining these with the rate coefficients for vibrational $v=1-5 \rightarrow v' < v$ quenching obtained with the infinite-order sudden approximation, we propose a new extrapolation scheme that yields the rate coefficients for ro-vibrational $v=2-5,j=0-30 \rightarrow v',j'$ de-excitation. Cross sections and rate coefficients for ro-vibrational $v=2, j=0-30 \rightarrow v'=1,j'$ transitions calculated with the close-coupling method confirm the effectiveness of this extrapolation scheme. Our calculated and extrapolated rates are very different from those that have been adopted in the modeling of many astrophysical environments. The current work provides the most comprehensive and accurate set of ro-vibrational de-excitation rate coefficients for the astrophysical modeling of the H-CO collision system. Application of the previously available and new data sets in astrophysical slab models shows that the line fluxes typically change by 20-70% in high temperature environments (800 K) with an H/H$_2$ ratio of 1; larger changes occur for lower temperatures.

Carbon monoxide is a simple molecule present in many astrophysical environments, and collisional excitation rate coefficients due to the dominant collision partners are necessary to accurately predict spectral line intensities and extract astrophysical parameters. We report new quantum scattering calculations for rotational deexcitation transitions of CO induced by H using the three-dimensional potential energy surface~(PES) of Song et al. (2015). State-to-state cross sections for collision energies from 10$^{-5}$ to 15,000~cm$^{-1}$ and rate coefficients for temperatures ranging from 1 to 3000~K are obtained for CO($v=0$, $j$) deexcitation from $j=1-45$ to all lower $j'$ levels, where $j$ is the rotational quantum number. Close-coupling and coupled-states calculations were performed in full-dimension for $j$=1-5, 10, 15, 20, 25, 30, 35, 40, and 45 while scaling approaches were used to estimate rate coefficients for all other intermediate rotational states. The current rate coefficients are compared with previous scattering results using earlier PESs. Astrophysical applications of the current results are briefly discussed.

Collisional excitation rate coefficients play an important role in the dynamics of energy transfer in the interstellar medium. In particular, accurate rotational excitation rates are needed to interpret microwave and infrared observations of the interstellar gas for nonlocal thermodynamic equilibrium line formation. Theoretical cross sections and rate coefficients for collisional deexcitation of rotationally excited HF in the vibrational ground state are reported. The quantum-mechanical close-coupling approach implemented in the nonreactive scattering code MOLSCAT was applied in the cross section and rate coefficient calculations on an accurate 2D HF-He potential energy surface. Estimates of rate coefficients for H and H$_2$ colliders were obtained from the HF-He collisional data with a reduced-potential scaling approach. The calculation of state-to-state rotational quenching cross sections for HF due to He with initial rotational levels up to $j=20$ were performed for kinetic energies from 10$^{-5}$ to 15000 cm$^{-1}$. State-to-state rate coefficients for temperatures between 0.1 and 3000 K are also presented. The comparison of the present results with previous work for lowly-excited rotational levels reveals significant differences. In estimating HF-H$_2$ rate coefficients, the reduced-potential method is found to be more reliable than the standard reduced-mass approach. The current state-to-state rate coefficient calculations are the most comprehensive to date for HF-He collisions. We attribute the differences between previously reported and our results to differences in the adopted interaction potential energy surfaces. The new He rate coefficients can be used in a variety of applications. The estimated H$_2$ and H collision rates can also augment the smaller datasets previously developed for H$_2$ and electrons.

Rate coefficients for collisional processes such as rotational and vibrational excitation are essential inputs in many astrophysical models. When rate coefficients are unknown, they are often estimated using known values from other systems. The most common example is to use He-collider rate coefficients to estimate values for other colliders, typically H$_2$, using scaling arguments based on the reduced mass of the collision system. This procedure is often justified by the assumption that the inelastic cross section is independent of the collider. Here we explore the validity of this approach focusing on rotational inelastic transitions for collisions of H, para-H$_2$, $^3$He, and $^4$He with CO in its vibrational ground state. We compare rate coefficients obtained via explicit calculations to those deduced by standard reduced-mass scaling. Not surprisingly, inelastic cross sections and rate coefficients are found to depend sensitively on both the reduced mass and the interaction potential energy surface. We demonstrate that standard reduced-mass scaling is not valid on physical and mathematical grounds, and as a consequence, the common approach of multiplying a rate coefficient for a molecule-He collision system by the constant factor of ~1.4 to estimate the rate coefficient for para-H$_2$ collisions is deemed unreliable. Furthermore, we test an alternative analytic scaling approach based on the strength of the interaction potential and the reduced mass of the collision systems. Any scaling approach, however, may be problematic when low-energy resonances are present; explicit calculations or measurements of rate coefficients are to be preferred.

We present observations of interstellar rubidium toward o Per, zeta Per, AE Aur, HD 147889, chi Oph, zeta Oph, and 20 Aql. Theory suggests that stable 85Rb and long-lived 87Rb are produced predominantly by high-mass stars, through a combination of the weak s- and r-processes. The 85Rb/87Rb ratio was determined from measurements of the Rb I line at 7800 angstroms and was compared to the solar system meteoritic ratio of 2.59. Within 1-sigma uncertainties all directions except HD 147889 have Rb isotope ratios consistent with the solar system value. The ratio toward HD 147889 is much lower than the meteoritic value and similar to that toward rho Oph A (Federman et al. 2004); both lines of sight probe the Rho Ophiuchus Molecular Cloud. The earlier result was attributed to a deficit of r-processed 85Rb. Our larger sample suggests instead that 87Rb is enhanced in these two lines of sight. When the total elemental abundance of Rb is compared to the K elemental abundance, the interstellar Rb/K ratio is significantly lower than the meteoritic ratio for all the sight lines in this study. Available interstellar samples for other s- and r- process elements are used to help interpret these results.

The results of new spectroscopic analyses of 30 stars with giant planet and/or brown dwarf companions are presented. Values for Teff and [Fe/H] are used in conjunction with Hipparcos data and Padova isochrones to derive masses, ages, and theoretical surface gravities. These new data are combined with spectroscopic and photometric metallicity estimates of other stars harboring planets and published samples of F, G, and K dwarfs to compare several subsets of planet bearing stars with similarly well-constrained control groups. The distribution of [Fe/H] values continues the trend uncovered in previous studies in that stars hosting planetary companions have a higher mean value than otherwise similar nearby stars. We also investigate the relationship between stellar mass and the presence of giant planets and find statistically marginal but suggestive evidence of a decrease in the incidence of radial velocity companions orbiting relatively less massive stars. If confirmed with larger samples, this would represent a critical constraint to both planetary formation models as well as to estimates of the distribution of planetary systems in our galaxy.

The polarized mm/sub-mm radiation from Sgr A* is apparently produced by a Keplerian structure whose peak emission occurs within several Schwarzschild radii (r_S=2GM/c^2) of the black hole. The Chandra X-ray counterpart, if confirmed, is presumably the self-Comptonized component from this region. In this paper, we suggest that sub-mm timing observations could yield a signal corresponding to the period P_0 of the marginally stable orbit, and therefore point directly to the black hole's spin a. Sgr A*'s mass is now known to be (2.6\pm 0.2)\times 10^6 M_\odot (an unusually accurate value for supermassive black hole candidates), for which 2.7 min<P_0<36 min, depending on the value of a and whether the Keplerian flow is prograde or retrograde. A Schwarzschild black hole (a=0) should have P_0 ~ 20 min. The identification of the orbital frequency with the innermost stable circular orbit is made feasible by the transition from optically thick to thin emission at sub-mm wavelengths. With stratification in the emitter, the peak of the sub-mm bump in Sgr A*'s spectrum is thus produced at the smallest radius. We caution, however, that theoretical uncertainties in the structure of the emission region may still produce some ambiguity in the timing signal. Given that Sgr A*'s flux at $\nu\sim 1$ mm is several Jy, these periods should lie within the temporal-resolving capability of sub-mm telescopes using bolometric detectors. A determination of P_0 should provide not only a value of a, but it should also define the angular momentum vector of the orbiting gas in relation to the black hole's spin axis. In addition, since the X-ray flux detected by Chandra appears to be the self-Comptonized mm to sub-mm component, these temporal fluctuations may also be evident in the X-ray signal.

We have performed searches for isolated flares and for steady flickering in Gamma Ray Burst light curves on the microsecond to millisecond time scales. Two bursts out of our sample of 20 revealed four isolated flares with time scales from $256 \mu s$ to $2048 \mu s$. A wavelet analysis for our sample showed low level flickering for all bursts on time scales from $256 \mu s$ to 33 ms, with the majority of bursts containing rise times faster than four milliseconds and 30% having rise times faster than one millisecond. These results show that millisecond variability is common in classical bursts and not some exceptional activity by a possibly separate class of bursts. These fast rise times can be used to place severe limits on burst models: (1) The characteristic thickness of the energy generation region must be less than 1200 km along the line of sight. (2) The angular size of the gamma ray emission region as subtended from the central source must be less than 42 arc-seconds. (3) The expanding ejecta must have a range of Lorentz factors along a radius line with a dispersion of less than roughly 2%. (4) Within the external shock scenario, the characteristic dimension of the impacted cloud must be smaller than 16 Astronomical Units on average. (5) Within the collimated jet scenario, the collimation angle must be smaller than 42 arc-seconds.

The Gamma Ray Burst GRB920229 has four extreme and unprecedented properties; a rise in brightness with an e-folding time scale of $220 \pm 30 \mu s$, a fall in brightness with an e-folding time scale of $400 \pm 100 \mu s$, a large change in spectral shape over a time of $768 \mu s$, and a sharp spectral cutoff to high energies with $\Delta E/E = 18 %$. The rapid changes occur during a spike in the light curve which was seen 0.164 s after the start of the burst. The spectrum has a peak $\nu F_{\nu}$ at 200 keV with no significant flux above 239 keV, although the cutoff energy shifts to less than 100 keV during the decay of the spike. These numbers can be used to place severe limits on fireball models of bursts. The thickness of the energy production region must be smaller than $\sim 66 km$, ejected shells must have a dispersion of the Lorentz factor of less than roughly 1% along a particular radius, and the angular size of the radiation emission region is of order 1 arc-minute as viewed from the burst center. The physical mechanism that caused the sharp spectral cutoff has not been determined.