Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

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 entertain the exotic possibility that dark matter (DM) decays or annihilations taking place in our galaxy may produce a flux of relativistic very weakly-coupled bosons, axions or dark photons. We show that there exist several upper bounds for this flux on Earth assuming generic minimal requirements for DM, such as a lifetime longer than the age of the Universe or an annihilation rate that leaves unaffected the background evolution during matter domination. These bounds do not depend on the identity or the couplings of the bosons. We then show that this new flux cannot be large enough to explain the recent XENON1T excess, while assuming that the bosons' couplings to the Standard Model are consistent with all current experimental and observational constraints. We also discuss a possible caveat to these bounds and a route to explain the excess.

The advent of the Gaia era has led to potentially revolutionary understanding of dark matter (DM) dynamics in our galaxy, which has important consequences for direct detection (DD) experiments. In this paper, we study how the empirical DM velocity distribution inferred from Gaia-Sausage, a dominant substructure in the solar neighborhood, affects the interpretation of DD data. We survey different classes of operators in the non-relativistic effective field theory that could arise from several relativistic benchmark models and emphasize that the Gaia velocity distribution could modify both the total number of events as well as the shape of the differential recoil spectra, the two primary observables in DD experiments. Employing the euclideanized signal method, we investigate the effects of the Gaia distribution on reconstructing DM model parameters and identifying the correct DM model given a positive signal at future DD experiments. We find that for light DM with mass ~10 GeV, the Gaia distribution poses an additional challenge for characterizing DM interactions with ordinary matter, which may be addressed by combining complementary DD experiments with different targets and lowering the detection threshold. Meanwhile, for heavy DM with mass above ~30 GeV, depending on the type of DM model, there could be a (moderate) improvement in the sensitivity at future DD experiments.

We present a publicly-available toolkit of flight-proven hardware and software to retrieve 5 TB of data or small physical samples from a stratospheric balloon platform. Before launch, a capsule is attached to the balloon, and rises with it. Upon remote command, the capsule is released and descends via parachute, continuously transmitting its location. Software to predict the trajectory can be used to select a safe but accessible landing site. We dropped two such capsules from the SuperBIT telescope, in September 2019. The capsules took ~37 minutes to descend from ~30 km altitude. They drifted 32 km and 19 km horizontally, but landed within 300 m and 600 m of their predicted landing sites. We found them easily, and successfully recovered the data. We welcome interest from other balloon teams for whom the technology would be useful.

At a fraction the total cost of an equivalent orbital mission, scientific balloon-borne platforms, operating above 99.7% of the Earth's atmosphere, offer attractive, competitive, and effective observational capabilities -- namely space-like resolution, transmission, and backgrounds -- that are well suited for modern astronomy and cosmology. SuperBIT is a diffraction-limited, wide-field, 0.5 m telescope capable of exploiting these observing conditions in order to provide exquisite imaging throughout the near-IR to near-UV. It utilizes a robust active stabilization system that has consistently demonstrated a 1 sigma sky-fixed pointing stability at 48 milliarcseconds over multiple 1 hour observations at float. This is achieved by actively tracking compound pendulations via a three-axis gimballed platform, which provides sky-fixed telescope stability at < 500 milliarcseconds and corrects for field rotation, while employing high-bandwidth tip/tilt optics to remove residual disturbances across the science imaging focal plane. SuperBIT's performance during the 2019 commissioning flight benefited from a customized high-fidelity science-capable telescope designed with exceptional thermo- and opto-mechanical stability as well as tightly constrained static and dynamic coupling between high-rate sensors and telescope optics. At the currently demonstrated level of flight performance, SuperBIT capabilities now surpass the science requirements for a wide variety of experiments in cosmology, astrophysics and stellar dynamics.

We use stellar kinematics from the latest Gaia data release (DR2) to measure the local dark matter density $\rho_{\rm DM}$ in a heliocentric cylinder of radius $R= 150 \ {\rm pc}$ and half-height $z= 200 \ {\rm pc}$. We also explore the prospect of using our analysis to estimate the DM density in local substructure by setting constraints on the surface density and scale height of a thin dark disk aligned with the baryonic disk and formed due to dark matter self-interaction. Performing the statistical analysis within a Bayesian framework for three types of tracers, we obtain ${\rho_{\rm DM}= 0.016 \pm 0.010}$ M$_\odot$/pc$^3$ for A stars; early G stars give a similar result, while F stars yield a significantly higher value. For a thin dark disk, A stars set the strongest constraint: excluding surface densities (5-12) M$_\odot$/pc$^2$ for scale heights below 100 pc with 95% confidence. Comparing our results with those derived using Tycho-Gaia Astrometric Solution (TGAS) data, we find that the uncertainty in our measurements of the local DM content is dominated by systematic errors that arise from assumptions of our kinematic analysis in the low $z$ region. Furthermore, there will only be a marginal reduction in these uncertainties with more data in the Gaia era. We comment on the robustness of our method and discuss potential improvements for future work.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

We use effective field theory (EFT) formalism to forecast the constraint on Horndeski class of dark energy models with future supernova and galaxy surveys. Previously (Gleyzes {\it et al.}) computed unmarginalized constraints (68\% CL error $\sim 10^{-3}$--$10^{-2}$) on EFT dark energy parameters by fixing all other parameters. We extend the previous work by allowing all cosmological parameters and nuisance parameters to vary and marginalizing over them. We find that (i) the constraints on EFT dark energy parameters are typically worsen by a factor of few after marginalization, and (ii) the constraint on the dark energy equation of state $w$ is not significantly affected by the inclusion of EFT dark energy parameters.