Using the NOvA neutrino detectors, a broad search has been performed for any signal coincident with 28 gravitational wave events detected by the LIGO/Virgo Collaboration between September 2015 and July 2019. For all of these events, NOvA is sensitive to possible arrival of neutrinos and cosmic rays of GeV and higher energies. For five (seven) events in the NOvA Far (Near) Detector, timely public alerts from the LIGO/Virgo Collaboration allowed recording of MeV-scale events. No signal candidates were found.

We present the first detection of the 13C17O J=3-2 transition toward the HL Tau protoplanetary disc. We find significantly more gas mass (at least a factor of ten higher) than has been previously reported using C18O emission. This brings the observed total disc mass to 0.2 M, which we consider to be a conservative lower limit. Our analysis of the Toomre Q profile suggests that this brings the disc into the regime of gravitational instability. The radial region of instability (50-110 au) coincides with the location of a proposed planet-carved gap in the dust disc and a spiral in the gas. We, therefore, propose that if the origin of the gap is confirmed to be due to a forming giant planet, then it is likely to have formed via the gravitational fragmentation of the protoplanetary disc.

Measurements of the gas mass are necessary to determine the planet formation potential of protoplanetary disks. Observations of rare CO isotopologues are typically used to determine disk gas masses; however, if the line emission is optically thick this will result in an underestimated disk mass. With ALMA we have detected the rarest stable CO isotopologue, 13C17O, in a protoplanetary disk for the first time. We compare our observations with the existing detections of 12CO, 13CO, C18O and C17O in the HD163296 disk. Radiative transfer modelling using a previously benchmarked model, and assuming interstellar isotopic abundances, significantly underestimates the integrated intensity of the 13C17O J=3-2 line. Reconciliation between the observations and the model requires a global increase in CO gas mass by a factor of 3.5. This is a factor of 2-6 larger than previous gas mass estimates using C18O. We find that C18O emission is optically thick within the CO snow line, while the 13C17O emission is optically thin and is thus a robust tracer of the bulk disk CO gas mass.

Observations of different molecular lines in protoplanetary disks provide valuable information on the gas kinematics, as well as constraints on the radial density and temperature structure of the gas. With ALMA we have detected H13CO+ (J=4-3) and HC15N (J=4-3) in the HD97048 protoplanetary disk for the first time. We compare these new detections to the ringed continuum mm-dust emission and the spatially resolved CO (J=3-2) and HCO+ (J=4-3) emission. The radial distributions of the H13CO+ and HC15N emission show hints of ringed sub-structure whereas, the optically thick tracers, CO and HCO+, do not. We calculate the HCO+/H13CO+ intensity ratio across the disk and find that it is radially constant (within our uncertainties). We use a physio-chemical parametric disk structure of the HD97048 disk with an analytical prescription for the HCO+ abundance distribution to generate synthetic observations of the HCO+ and H13CO+ disk emission assuming LTE. The best by-eye fit models require radial variations in the HCO+/H13CO+ abundance ratio and an overall enhancement in H13CO+ relative to HCO+. This highlights the need to consider isotope selective chemistry and in particular low temperature carbon isotope exchange reactions. This also points to the presence of a reservoir of cold molecular gas in the outer disk (T < 10K, R > 200au). Chemical models are required to confirm that isotope-selective chemistry alone can explain the observations presented here. With these data, we cannot rule out that the known dust substructure in the HD97048 disk is responsible for the observed trends in molecular line emission, and higher spatial resolution observations are required to fully explore the potential of optically thin tracers to probe planet-carved dust gaps. We also report non-detections of H13CO+ and HC15N in the HD100546 protoplanetary disk.

Using two years of data from the NOvA Near Detector at Fermilab, we report a seasonal variation of cosmic ray induced multiple-muon event rates which has an opposite phase to the seasonal variation in the atmospheric temperature. The strength of the seasonal multipl$ increase as a function of the muon multiplicity. However, no significant dependence of the strength of the seasonal variation of the multiple-muon variation is seen as a function of the muon zenith angle, or the spatial or angular separation between the correlated muons.

In this paper, we analyse the upper limit fluxes of sub-millimeter ortho-H$_{2}$$^{16}$O 321 GHz, para-H$_{2}$$^{18}$O 322 GHz, and HDO 335 GHz lines from the protoplanetary disk around the Herbig Ae star HD 163296, using the Atacama Large Millimeter/Submillimeter Array (ALMA). These water lines are considered to be the best candidate sub-millimeter lines to locate the position of the H$_{2}$O snowline, on the basis of our previous model calculations. We compare the upper limit fluxes with the values calculated by our models with dust emission included, and we constrain the line emitting region and the dust opacity from the observations. We conclude that, if the outer edge of the region with high water vapor abundance and if the position of the water snowline are beyond 8 au also, the mm dust opacity will have a value larger than 2.0 cm$^{2}$ g$^{-1}$. In addition, the position of the water snowline will be inside 20 au, if the mm dust opacity is 2.0 cm$^{2}$ g$^{-1}$. Future observations of the dust continuum emission at higher angular resolution and sub-millimeter water lines with longer observation time are required to clarify the detailed structures and the position of the H$_{2}$O snowline in the disk midplane.

Sulphur-bearing volatiles are observed to be significantly depleted in interstellar and circumstellar regions. This missing sulphur is postulated to be mostly locked up in refractory form. With ALMA we have detected sulphur monoxide (SO), a known shock tracer, in the HD 100546 protoplanetary disk. Two rotational transitions: $J=7_{7}-6_{6}$ (301.286 GHz) and $J=7_8-6_7$ (304.078 GHz) are detected in their respective integrated intensity maps. The stacking of these transitions results in a clear 5$\sigma$ detection in the stacked line profile. The emission is compact but is spectrally resolved and the line profile has two components. One component peaks at the source velocity and the other is blue-shifted by 5 km s$^{-1}$. The kinematics and spatial distribution of the SO emission are not consistent with that expected from a purely Keplerian disk. We detect additional blue-shifted emission that we attribute to a disk wind. The disk component was simulated using LIME and a physical disk structure. The disk emission is asymmetric and best fit by a wedge of emission in the north east region of the disk coincident with a `hot-spot' observed in the CO $J=3-2$ line. The favoured hypothesis is that a possible inner disk warp (seen in CO emission) directly exposes the north-east side of the disk to heating by the central star, creating locally the conditions to launch a disk wind. Chemical models of a disk wind will help to elucidate why the wind is particularly highlighted in SO emission and whether a refractory source of sulphur is needed. An alternative explanation is that the SO is tracing an accretion shock from a circumplanetary disk associated with the proposed protoplanet embedded in the disk at 50 au. We also report a non-detection of SO in the protoplanetary disk around HD 97048.

We set out to study the atmosphere of WASP-80b, a warm inflated gas giant with an equilibrium temperature of $\sim$800~K, using ground-based transmission spectroscopy covering the spectral range from 520~to~910~nm. The observations allow us to probe the existence and abundance of K and Na in WASP-80b's atmosphere, existence of high-altitude clouds, and Rayleigh-scattering in the blue end of the spectrum. We observed two spectroscopic time series of WASP-80b transits with the OSIRIS spectrograph installed in the Gran Telescopio CANARIAS, and use the observations to estimate the planet's transmission spectrum between 520~nm and 910~nm in 20~nm-wide passbands, and around the K~I and Na~I resonance doublets in 6~nm-wide passbands. We model three previously published broadband datasets consisting of 27 light curves jointly prior to the transmission spectroscopy analysis in order to obtain improved prior estimates for the planet's orbital parameters, average radius ratio, and stellar density. We recover a flat transmission spectrum with no evidence of Rayleigh scattering or K~I or Na~I absorption, and obtain an improved system characterisation as a by-product of the broadband- and GTC-dataset modelling. The transmission spectra estimated separately from the two observing runs are consistent with each other, as are the transmission spectra estimated using either a parametric or nonparametric systematics models. The flat transmission spectrum favours an atmosphere model with high-altitude clouds over cloud-free models with stellar or sub-stellar metallicities.

The Keck Interferometer Nuller (KIN) was used to survey 25 nearby main sequence stars in the mid-infrared, in order to assess the prevalence of warm circumstellar (exozodiacal) dust around nearby solar-type stars. The KIN measures circumstellar emission by spatially blocking the star but transmitting the circumstellar flux in a region typically 0.1 - 4 AU from the star. We find one significant detection (eta Crv), two marginal detections (gamma Oph and alpha Aql), and 22 clear non-detections. Using a model of our own Solar System's zodiacal cloud, scaled to the luminosity of each target star, we estimate the equivalent number of target zodis needed to match our observations. Our three zodi detections are eta Crv (1250 +/- 260), gamma Oph (200 +/- 80) and alpha Aql (600 +/- 200), where the uncertainties are 1-sigma. The 22 non-detected targets have an ensemble weighted average consistent with zero, with an average individual uncertainty of 160 zodis (1-sigma). These measurements represent the best limits to date on exozodi levels for a sample of nearby main sequence stars. A statistical analysis of the population of 23 stars not previously known to contain circumstellar dust (excluding eta Crv and gamma Oph) suggests that, if the measurement errors are uncorrelated (for which we provide evidence) and if these 23 stars are representative of a single class with respect to the level of exozodi brightness, the mean exozodi level for the class is <150 zodis (3-sigma upper-limit, corresponding to 99% confidence under the additional assumption that the measurement errors are Gaussian). We also demonstrate that this conclusion is largely independent of the shape and mean level of the (unknown) true underlying exozodi distribution.

We report observations of the nova RS Ophiuchi (RS Oph) using the Keck Interferometer Nuller (KIN), approximately 3.8 days following the most recent outburst that occurred on 2006 February 12. These observations represent the first scientific results from the KIN, which operates in N-band from 8 to 12.5 microns in a nulling mode. By fitting the unique KIN data, we have obtained an angular size of the mid-infrared continuum of 6.2, 4.0, or 5.4 mas for a disk profile, gaussian profile (FWHM), and shell profile respectively. The data show evidence of enhanced neutral atomic hydrogen emission and atomic metals including silicon located in the inner spatial regime near the white dwarf (WD) relative to the outer regime. There are also nebular emission lines and evidence of hot silicate dust in the outer spatial region, centered at ! 17 AU from the WD, that are not found in the inner regime. Our evidence suggests that these features have been excited by the nova flash in the outer spatial regime before the blast wave reached these regions. These identifications support a model in which the dust appears to be present between outbursts and is not created during the outburst event. We further discuss the present results in terms of a unifying model of the system that includes an increase in density in the plane of the orbit of the two stars created by a spiral shock wave caused by the motion of the stars through the cool wind of the red giant star. These data show the power and potential of the nulling technique which has been developed for the detection of Earth-like planets around nearby stars for the Terrestrial Planet Finder Mission and Darwin missions.

New angular diameter determinations for the bright southern F8 supergiant Delta CMa enable the bolometric emergent flux and effective temperature of the star to be determined with improved accuracy. The spectral flux distribution and bolometric flux have been determined from published photometry and spectrophotometry and combined with the angular diameter to derive the bolometric emergent flux F = (6.50 plus/minus 0.24) x 10^7 W/m^2 and the effective temperature Teff = 5818 plus/minus 53 K. The new value for the effective temperature is compared with previous interferometric and infrared flux method determinations. The accuracy of the effective temperature is now limited by the uncertainty in the bolometric flux rather than by the uncertainty in the angular diameter.

Using a new grism at the Keck Interferometer, we obtained spectrally dispersed (R ~ 230) interferometric measurements of the Mira star R Vir. These data show that the measured radius of the emission varies substantially from 2.0-2.4 microns. Simple models can reproduce these wavelength-dependent variations using extended molecular layers, which absorb stellar radiation and re-emit it at longer wavelengths. Because we observe spectral regions with and without substantial molecular opacity, we determine the stellar photospheric radius, uncontaminated by molecular emission. We infer that most of the molecular opacity arises at approximately twice the radius of the stellar photosphere.

We present new K-band long baseline interferometer observations of three young stellar objects of the FU Orionis class, V1057 Cyg, V1515 Cyg and Z CMa-SE, obtained at the Keck Interferometer during its commissioning science period. The interferometer clearly resolves the source of near-infrared emission in all three objects. Using simple geometrical models we derive size scales (0.5-4.5 AU) for this emission. All three objects appear significantly more resolved than expected from simple models of accretion disks tuned to fit the broadband optical and infrared spectro-photometry. We explore variations in the key parameters that are able to lower the predicted visibility amplitudes to the measured levels, and conclude that accretion disks alone do not reproduce the spectral energy distributions and K-band visibilities simultaneously. We conclude that either disk models are inadequate to describe the near-infrared emission, or additional source components are needed. We hypothesize that large scale emission (10s of AU) in the interferometer field of view is responsible for the surprisingly low visibilities. This emission may arise in scattering by large envelopes believed to surround these objects.

We report the results of a sensitive K-band survey of Herbig Ae/Be disk sizes using the 85-m baseline Keck Interferometer. Targets were chosen to span the maximum range of stellar properties to probe the disk size dependence on luminosity and effective temperature. For most targets, the measured near-infrared sizes (ranging from 0.2 to 4 AU) support a simple disk model possessing a central optically-thin (dust-free) cavity, ringed by hot dust emitting at the expected sublimation temperatures (T_sub~1000-1500K). Furthermore, we find a tight correlation of disk size with source luminosity R propto L^(1/2) for Ae and late Be systems (valid over more than 2 decades in luminosity), confirming earlier suggestions based on lower-quality data. Interestingly, the inferred dust-free inner cavities of the highest luminosity sources (Herbig B0-B3 stars) are under-sized compared to predictions of the optically-thin cavity model, likely due to optically-thick gas within the inner AU.

The bright southern binary star beta Centauri (HR 5267) has been observed with the Sydney University Stellar Interferometer (SUSI) and spectroscopically with the ESO CAT and Swiss Euler telescopes at La Silla. The interferometric observations have confirmed the binary nature of the primary component and have enabled the determination of the orbital parameters of the system. At the observing wavelength of 442 nm the two components of the binary system have a magnitude difference of 0.15. The combination of interferometric and spectroscopic data gives the following results: orbital period 357 days, semi-major axis 25.30 mas, inclination 67.4 degrees, eccentricity 0.821, distance 102.3 pc, primary and secondary masses M1 = M2 = 9.1 solar masses and absolute visual magnitudes of the primary and secondary M1V = -3.85 and M2V = -3.70. The high accuracy of the results offers a fruitful starting point for future asteroseismic modelling of the pulsating binary components.

We report novel, high-angular resolution interferometric measurements that imply the near-infrared nuclear emission in NGC 4151 is unexpectedly compact. We have observed the nucleus of NGC 4151 at 2.2 microns using the two 10-meter Keck telescopes as an interferometer and find a marginally resolved source ~0.1 pc in diameter. Our measurements rule out models in which a majority of the K band nuclear emission is produced on scales larger than this size. The interpretation of our measurement most consistent with other observations is that the emission mainly originates directly in the central accretion disk. This implies that AGN unification models invoking hot, optically thick dust may not be applicable to NGC 4151.