Search

This note provides additional interpretation for the counterfactual outcome distribution and corresponding unconditional quantile "effects" defined and estimated by Firpo, Fortin, and Lemieux (2009) and Chernozhukov, Fern\'andez-Val, and Melly (2013). With conditional independence of the policy variable of interest, these methods estimate the policy effect for certain types of policies, but not others. In particular, they estimate the effect of a policy change that itself satisfies conditional independence., Comment: The main result is essentially the same as Proposition 1 of Rothe (2010), https://doi.org/10.1016/j.jeconom.2009.09.001

A new technique to search for new scalar and tensor interactions at the sub-micrometer scale is presented. The technique relies on small shifts of nuclear gamma lines produced by the coupling between matter and the nuclei in the source or absorber of a Mossbauer spectrometer. Remarkably, such energy shifts are rather insensitive to electromagnetic interactions that represent the largest background in searches for new forces using atomic matter. This is because nuclei are intrinsically shielded by the electron clouds. Additionally, electromagnetic interactions cause energy shifts by coupling to nuclear moments that are suppressed by the size of the nuclei, while new scalar interactions can directly affect these shifts. Finally, averaging over unpolarized nuclei, further reduces electromagnetic interactions. We discuss several possible configurations, using the traditional Mossbauer effect as well as nuclear resonant absorption driven by synchrotron radiation. For this purpose, we examine the viability of well known Mossbauer nuclides along with more exotic ones that result in substantially narrower resonances. We find that the technique introduced here could substantially improve the sensitivity to a variety of new interactions and could also be used, in conjunction with mechanical force measurements, to corroborate a discovery or explore the new physics that may be behind a discovery., Comment: 11 pages, 7 figures

We report the discovery of ZTF J2243+5242, an eclipsing double white dwarf binary with an orbital period of just $8.8$ minutes, the second known eclipsing binary with an orbital period less than ten minutes. The system likely consists of two low-mass white dwarfs, and will merge in approximately 400,000 years to form either an isolated hot subdwarf or an R Coronae Borealis star. Like its $6.91\, \rm min$ counterpart, ZTF J1539+5027, ZTF J2243+5242 will be among the strongest gravitational wave sources detectable by the space-based gravitational-wave detector The Laser Space Interferometer Antenna (LISA) because its gravitational-wave frequency falls near the peak of LISA's sensitivity. Based on its estimated distance of $d=2120^{+131}_{-115}\,\rm pc$, LISA should detect the source within its first few months of operation, and should achieve a signal-to-noise ratio of $87\pm5$ after four years. We find component masses of $M_A= 0.349^{+0.093}_{-0.074}\,M_\odot$ and $M_B=0.384^{+0.114}_{-0.074}\,M_\odot$, radii of $R_A=0.0308^{+0.0026}_{-0.0025}\,R_\odot$ and $R_B = 0.0291^{+0.0032}_{-0.0024}\,R_\odot$, and effective temperatures of $T_A=22200^{+1800}_{-1600}\,\rm K$ and $T_B=16200^{+1200}_{-1000}\,\rm K$. We determined all of these properties, and the distance to this system, using only photometric measurements, demonstrating a feasible way to estimate parameters for the large population of optically faint ($r>21 \, m_{\rm AB}$) gravitational-wave sources which the Vera Rubin Observatory (VRO) and LISA should identify., Comment: 12 pages, 7 figures, 3 tables, submitted

The current generation of all-sky surveys is rapidly expanding our ability to study variable and transient sources. These surveys, with a variety of sensitivities, cadences, and fields of view, probe many ranges of timescale and magnitude. Data from the Zwicky Transient Facility (ZTF) yields an opportunity to find variables on timescales from minutes to months. In this paper, we present the codebase, ztfperiodic, and the computational metrics employed for the catalogue based on ZTF's Second Data Release. We describe the publicly available, graphical-process-unit optimized period-finding algorithms employed, and highlight the benefit of existing and future graphical-process-unit clusters. We show how generating metrics as input to catalogues of this scale is possible for future ZTF data releases. Further work will be needed for future data from the Vera C. Rubin Observatory's Legacy Survey of Space and Time.

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project has the primary goal of detecting and characterizing low-frequency gravitational waves through high-precision pulsar timing. The mitigation of interstellar effects is crucial to achieve the necessary precision for gravitational wave detection. Effects like dispersion and scattering are more influential at lower observing frequencies, with the variation of these quantities over week-month timescales requiring high-cadence multi-frequency observations for pulsar timing projects. In this work, we utilize the dual-frequency observing capability of the Giant Metrewave Radio Telescope (GMRT) and evaluate the potential decrease in dispersion measure (DM) uncertainties when combined with existing pulsar timing array data. We present the timing analysis for four millisecond pulsars observed with the GMRT simultaneously at 322 and 607 MHz, and compare the DM measurements with those obtained through NANOGrav observations with the Green Bank Telescope (GBT) and Arecibo Observatory at 1400 to 2300 MHz frequencies. Measured DM values with the GMRT and NANOGrav program show significant offsets for some pulsars, which could be caused by pulse profile evolution in the two frequency bands. In comparison to the predicted DM uncertainties when incorporating these low-frequency data into the NANOGrav dataset, we find that higher-precision GMRT data is necessary to provide improved DM measurements. Through the detection and analysis of pulse profile baseline ripple in data on test pulsar B1929+10, we find that, while not important for this data, it may be relevant for other timing datasets. We discuss the possible advantages and challenges of incorporating GMRT data into NANOGrav and International Pulsar Timing Array datasets., Comment: 14 pages, 8 figures. Accepted to The Astrophysical Journal

We search for an isotropic stochastic gravitational-wave background (GWB) in the $12.5$-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. Our analysis finds strong evidence of a stochastic process, modeled as a power-law, with common amplitude and spectral slope across pulsars. The Bayesian posterior of the amplitude for an $f^{-2/3}$ power-law spectrum, expressed as the characteristic GW strain, has median $1.92 \times 10^{-15}$ and $5\%$--$95\%$ quantiles of $1.37$--$2.67 \times 10^{-15}$ at a reference frequency of $f_\mathrm{yr} = 1 ~\mathrm{yr}^{-1}$. The Bayes factor in favor of the common-spectrum process versus independent red-noise processes in each pulsar exceeds $10,000$. However, we find no statistically significant evidence that this process has quadrupolar spatial correlations, which we would consider necessary to claim a GWB detection consistent with general relativity. We find that the process has neither monopolar nor dipolar correlations, which may arise from, for example, reference clock or solar system ephemeris systematics, respectively. The amplitude posterior has significant support above previously reported upper limits; we explain this in terms of the Bayesian priors assumed for intrinsic pulsar red noise. We examine potential implications for the supermassive black hole binary population under the hypothesis that the signal is indeed astrophysical in nature., Comment: 25 pages, 14 figures, 5 tables, 3 appendices. Published in The Astrophysical Journal Letters. Please send any comments/questions to Joseph Simon (joe.simon@nanograv.org). Jupyter notebook tutorials and some MCMC chain files are available at https://github.com/nanograv/12p5yr_stochastic_analysis

Using photometry collected with the Zwicky Transient Facility (ZTF), we are conducting an ongoing survey for binary systems with short orbital periods ($P_{\rm b}<1\rm \,hr)$ with the goal of identifying new gravitational-wave sources detectable by the upcoming Laser Interferometer Space Antenna (LISA). Here, we present a sample of fifteen binary systems discovered thus far, with orbital periods ranging from $6.91\rm\,min$ to $56.35\rm\,min$. Of the fifteen systems, seven are eclipsing systems which do not show signs of significant mass transfer. Additionally, we have discovered two AM Canum Venaticorum (AM CVn) systems and six systems exhibiting primarily ellipsoidal variations in their light curves. We present follow-up spectroscopy and high-speed photometry confirming the nature of these systems, estimates of their LISA signal-to-noise ratios (SNR), and a discussion of their physical characteristics., Comment: 32 pages, 12 figures, 7 tables, submitted to ApJ

The first binary neutron star merger, GW170817, was accompanied by a radioactivity-powered optical/infrared transient called a kilonova. To date, no compelling kilonova has been found during optical surveys of the sky, independent of gravitational-wave triggers. In this work, we searched the first 23 months of the Zwicky Transient Facility (ZTF) data stream for candidate kilonovae in the form of rapidly evolving transients. We combined ZTF alert queries with forced point-spread-function photometry and nightly flux stacking to increase our sensitivity to faint and fast transients. Automatic queries yielded $>11,200$ candidates, 24 of which passed quality checks and strict selection criteria based on a grid of kilonova models tailored for both binary neutron star and neutron star-black hole mergers. None of the candidates in our sample was deemed a possible kilonova after thorough vetting, catalog cross-matching, and study of their color evolution. The sources that passed our selection criteria are dominated by Galactic cataclysmic variables. In addition, we identified two fast transients at high Galactic latitude, one of which is the confirmed afterglow of long-duration GRB190106A, and the other is a possible cosmological afterglow. Using a survey simulation code, we constrained the kilonova rate for a range of models including top-hat and linearly decaying light curves and synthetic light curves obtained with radiative transfer simulations. For prototypical GW170817-like kilonovae, we constrain the rate to be $R < 1775$ Gpc$^{-3}$ yr$^{-1}$ at 95% confidence level by requiring at least 2 high-significance detections. By assuming a population of kilonovae with the same geometry and composition of GW170817 observed under a uniform viewing angle distribution, we obtained a constraint on the rate of $R < 4029$ Gpc$^{-3}$ yr$^{-1}$., Comment: Submitted for publication in ApJ

A solution to the black hole information problem requires propagation of information from the interior of the black hole to the exterior. Such propagation violates general relativity and could conceivably be accomplished through "firewall" models. Based on the existence of similar firewalls at the inner horizons of charged and rotating black holes, a model of a firewall was recently constructed where the exterior spacetime reduces to that of the Schwarzschild metric but with a dramatically different interior. We investigate the radial and nonradial polar stability of these objects. We first study the dynamics of the shell under spherically symmetric perturbations, and impose constraints on the firewall model parameters by requiring a subluminal speed of sound on the firewall. We show that the demands of stability and subluminality impose significant constraints on the internal parameters of the firewall, narrowing down the range of objects that could be used to create such a structure., Comment: 18 pages, 11 figures

Present gravitational wave detectors are based on the measurement of linear displacement in stable optical cavities. Here, we instead suggest the measurement of the twist of a chiral mechanical element induced by a gravitational wave. The induced twist rotates a flat optical mirror on top of this chiral element, leading to the deflection of an incident laser beam. This angle change is enhanced by multiple bounces of light between the rotating mirror and an originally parallel nearby fixed flat mirror. Based on detailed continuum-mechanics calculations, we present a feasible design for the chiral mechanical element including the rotating mirror. Our approach is most useful for signals in the frequency band 1 -- 100 kHz where we show that fundamental metrological limits would allow for smaller shot noise in this setup in comparison to the detection of linear displacement. We estimate a gravitational wave strain sensitivity between 10^{-21}/\sqrt{Hz} and 10^{-23}/\sqrt{Hz} at around 10 kHz frequency. When appropriately scaling the involved geometrical parameters, the strain sensitivity is proportional to frequency., Comment: 6 pages, 6 figures

We report the results of the rapid follow-up observations of gamma-ray bursts (GRBs) detected by the Fermi satellite to search for associated fast radio bursts. The observations were conducted with the Australian Square Kilometre Array Pathfinder at frequencies from 1.2-1.4 GHz. A set of 20 bursts, of which four were short GRBs, were followed up with a typical latency of about one minute, for a duration of up to 11 hours after the burst. The data was searched using 4096 dispersion measure trials up to a maximum dispersion measure of 3763 pc cm$^{-3}$, and for pulse widths $w$ over a range of duration from 1.256 to 40.48 ms. No associated pulsed radio emission was observed above $26 {\rm Jy ms} (w/1 {\rm ms})^{-1/2}$ for any of the 20 GRBs., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged), Comment: Submitted to ApJ

We present the full panchromatic afterglow light curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data, and reprocessing a subset of it, we have evaluated the impact of differences in data processing or flux determination methods used by different groups, and attempted to mitigate these differences to provide a more uniform dataset. Simple power-law fits to the uniform afterglow light curve indicate a $t^{0.86\pm0.04}$ rise, a $t^{-1.92\pm0.12}$ decline, and a peak occurring at $155\pm4$ days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index ($-0.584\pm0.002$) across all epochs. This gives a precise and updated estimate of the electron power-law index, $p=2.168\pm0.004$. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized ISM density, $<$0.01 cm$^{-3}$, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having magnetic field strength between 10$^{10.4}$ G and 10$^{16}$ G. Our fits to the afterglow data using an analytical model that includes VLBI proper motion from Mooley et al. (2018), and a structured jet model that ignores the proper motion, indicates that the proper motion measurement needs to be considered while seeking an accurate estimate of the viewing angle., Comment: Accepted by ApJ. Revised X-ray and up-to-date afterglow measurements are available in ASCII format on GitHub (https://github.com/kmooley/GW170817/) and mirrored on tauceti.caltech.edu (http://www.tauceti.caltech.edu/kunal/gw170817/). The online dataset will be continuously updated as new measurements become available

We show that proton storage ring experiments designed to search for proton electric dipole moments can also be used to look for the nearly dc spin precession induced by dark energy and ultra-light dark matter. These experiments are sensitive to both axion-like and vector fields. Current technology permits probes of these phenomena up to three orders of magnitude beyond astrophysical limits. The relativistic boost of the protons in these rings allows this scheme to have sensitivities comparable to atomic co-magnetometer experiments that can also probe similar phenomena. These complementary approaches can be used to extract the micro-physics of a signal, allowing us to distinguish between pseudo-scalar, magnetic and electric dipole moment interactions., Comment: 19 pages, 7 figures

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences., Comment: 14 pages, 6 figures. Accepted by ApJ

We present a new analysis of the profile data from the 47 millisecond pulsars comprising the 12.5-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), which is presented in a parallel paper (Alam et al. 2021a; NG12.5). Our reprocessing is performed using "wideband" timing methods, which use frequency-dependent template profiles, simultaneous time-of-arrival (TOA) and dispersion measure (DM) measurements from broadband observations, and novel analysis techniques. In particular, the wideband DM measurements are used to constrain the DM portion of the timing model. We compare the ensemble timing results to NG12.5 by examining the timing residuals, timing models, and noise model components. There is a remarkable level of agreement across all metrics considered. Our best-timed pulsars produce encouragingly similar results to those from NG12.5. In certain cases, such as high-DM pulsars with profile broadening, or sources that are weak and scintillating, wideband timing techniques prove to be beneficial, leading to more precise timing model parameters by 10-15%. The high-precision, multi-band measurements of several pulsars indicate frequency-dependent DMs. Compared to the narrowband analysis in NG12.5, the TOA volume is reduced by a factor of 33, which may ultimately facilitate computational speed-ups for complex pulsar timing array analyses. This first wideband pulsar timing data set is a stepping stone, and its consistent results with NG12.5 assure us that such data sets are appropriate for gravitational wave analyses., Comment: 62 pages, 55 figures, 5 tables, 3 appendices. Data available at http://nanograv.org/data/ and via DOI 10.5281/zenodo.4312887

We present time-of-arrival (TOA) measurements and timing models of 47 millisecond pulsars (MSPs) observed from 2004 to 2017 at the Arecibo Observatory and the Green Bank Telescope by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). The observing cadence was three to four weeks for most pulsars over most of this time span, with weekly observations of six sources. These data were collected for use in low-frequency gravitational wave searches and for other astrophysical purposes. We detail our observational methods and present a set of TOA measurements, based on "narrowband" analysis, in which many TOAs are calculated within narrow radio-frequency bands for data collected simultaneously across a wide bandwidth. A separate set of "wideband" TOAs will be presented in a companion paper. We detail a number of methodological changes compared to our previous work which yield a cleaner and more uniformly processed data set. Our timing models include several new astrometric and binary pulsar measurements, including previously unpublished values for the parallaxes of PSRs J1832-0836 and J2322+2057, the secular derivatives of the projected semi-major orbital axes of PSRs J0613-0200 and J2229+2643, and the first detection of the Shapiro delay in PSR J2145-0750. We report detectable levels of red noise in the time series for 14 pulsars. As a check on timing model reliability, we investigate the stability of astrometric parameters across data sets of different lengths. We report flux density measurements for all pulsars observed. Searches for stochastic and continuous gravitational waves using these data will be subjects of forthcoming publications., Comment: 54 pages, 52 figures, 7 tables, 1 appendix. Data are available at http://nanograv.org/data/ and via the DOI 10.5281/zenodo.4312297

Recently, Kundu et al (2019) reported that the globular cluster NGC 5024 (M53) possesses five extra-tidal RR Lyrae. In fact, four of them were instead known members of a nearby globular cluster NGC 5053. The status of the remaining extra-tidal RR Lyrae is controversial depending on the adopted tidal radius of NGC 5024. We have also searched for additional extra-tidal RR Lyrae within an area of $\sim8$~deg$^2$ covering both globular clusters. This includes other known RR Lyrae within the search area, as well as stars that fall within the expected range of magnitudes and colors for RR Lyrae (and yet outside the cutoff of 2/3 of the tidal radii of each globular clusters for something to be called "extra-tidal") if they were extra-tidal RR Lyrae candidates for NGC 5024 or NGC 5053. Based on the the proper motion information and their locations on the color-magnitude diagram, none of the known RR Lyrae belong to the extra-tidal RR Lyrae of either globular clusters. In the cases where the stars satisfied the magnitude and color ranges of RR Lyrae, analysis of time series data taken from the Zwicky Transient Facility do not reveal periodicities, suggesting that none of these stars are RR Lyrae. We conclude that there are no extra-tidal RR Lyrae associated with either NGC 5024 or NGC 5053 located within our search area., Comment: 17 pages, 2 Tables, 12 Figures, AJ accepted

The discovery rate of fast radio bursts (FRBs) is increasing dramatically thanks to new radio facilities. Meanwhile, wide-field instruments such as the 47 deg$^2$ Zwicky Transient Facility (ZTF) survey the optical sky to study transient and variable sources. We present serendipitous ZTF observations of the CHIME repeating source FRB 180916.J0158+65, that was localized to a spiral galaxy 149 Mpc away and is the first FRB suggesting periodic modulation in its activity. While 147 ZTF exposures corresponded to expected high-activity periods of this FRB, no single ZTF exposure was at the same time as a CHIME detection. No $>3\sigma$ optical source was found at the FRB location in 683 ZTF exposures, totalling 5.69 hours of integration time. We combined ZTF upper limits and expected repetitions from FRB 180916.J0158+65 in a statistical framework using a Weibull distribution, agnostic of periodic modulation priors. The analysis yielded a constraint on the ratio between the optical and radio fluences of $\eta \lesssim 200$, corresponding to an optical energy $E_{\rm opt} \lesssim 3 \times 10^{46}$ erg for a fiducial 10 Jy ms FRB (90% confidence). A deeper (but less statistically robust) constraint of $\eta \lesssim 3$ can be placed assuming a rate of $r(>5$ Jy ms)= hr$^{-1}$ and $1.2\pm 1.1$ FRB occurring during exposures taken in high-activity windows. The constraint can be improved with shorter per-image exposures and longer integration time, or observing FRBs at higher Galactic latitudes. This work demonstrated how current surveys can statistically constrain multi-wavelength counterparts to FRBs even without deliberately scheduled simultaneous radio observation., Comment: Accepted for publication in ApJL, 9 pages, 4 figures, 1 table

PSR J1012+5307, a millisecond pulsar in orbit with a helium white dwarf (WD), has been timed with high precision for about 25 years. One of the main objectives of this long-term timing is to use the large asymmetry in gravitational binding energy between the neutron star and the WD to test gravitational theories. Such tests, however, will be eventually limited by the accuracy of the distance to the pulsar. Here, we present VLBI (very long baseline interferometry) astrometry results spanning approximately 2.5 years for PSR J1012+5307, obtained with the Very Long Baseline Array as part of the MSPSRPI project. These provide the first proper motion and absolute position for PSR J1012+5307 measured in a quasi-inertial reference frame. From the VLBI results, we measure a distance of $0.83^{+0.06}_{-0.02}$kpc (all the estimates presented in the abstract are at 68% confidence) for PSR J1012+5307, which is the most precise obtained to date. Using the new distance, we improve the uncertainty of measurements of the unmodeled contributions to orbital period decay, which, combined with three other pulsars, places new constraints on the coupling constant for dipole gravitational radiation $\kappa_D=(-1.7\pm1.7)\times 10^{-4}$ and the fractional time derivative of Newton's gravitational constant $\dot{G}/G = -1.8^{\,+5.6}_{\,-4.7}\times 10^{-13}\,{\rm yr^{-1}}$ in the local universe. As the uncertainties of the observed decays of orbital period for the four leading pulsar-WD systems become negligible in $\approx10$ years, the uncertainties for $\dot{G}/G$ and $\kappa_D$ will be improved to $\leq1.5\times10^{-13}\,{\rm yr^{-1}}$ and $\leq1.0\times10^{-4}$, respectively, predominantly limited by the distance uncertainties., Comment: published in ApJ (2020ApJ...896...85D)

Using selection criteria based on amplitude, time and color, we have identified 329 objects as known or candidate cataclysmic variable (CVs) during the first year of testing and operation of the Zwicky Transient Facility (ZTF). Of these, 90 are previously confirmed CVs, 218 are strong candidates based on the shape and color of their light curves obtained during 3-562 days of observations, and the remaining 21 are possible CVs but with too few data points to be listed as good candidates. Almost half the strong candidates are within 10 deg of the galactic plane, in contrast to most other large surveys which have avoided crowded fields. The available Gaia parallaxes are consistent with sampling the low mass transfer CVs, as predicted by population models. Our followup spectra have confirmed Balmer/helium emission lines in 27 objects, with four showing high excitation HeII emission, including candidates for an AM CVn, a polar and an intermediate polar. Our results demonstrate that a complete survey of the galactic plane is needed to accomplish an accurate determination of the number of CVs existing in the Milky Way., Comment: 27 pages, 9 figures, 3 tables

We report the discovery of the first short period binary in which a hot subdwarf star (sdOB) fills its Roche lobe and started mass transfer to its companion. The object was discovered as part of a dedicated high-cadence survey of the Galactic Plane named the Zwicky Transient Facility and exhibits a period of $P_{\rm orb}=39.3401(1)$ min, making it the most compact hot subdwarf binary currently known. Spectroscopic observations are consistent with an intermediate He-sdOB star with an effective temperature of $T_{\rm eff}=42,400\pm300$ K and a surface gravity of $\log(g)=5.77\pm0.05$. A high-signal-to noise GTC+HiPERCAM light curve is dominated by the ellipsoidal deformation of the sdOB star and an eclipse of the sdOB by an accretion disk. We infer a low-mass hot subdwarf donor with a mass $M_{\rm sdOB}=0.337\pm0.015$ M$_\odot$ and a white dwarf accretor with a mass $M_{\rm WD}=0.545\pm0.020$ M$_\odot$. Theoretical binary modeling indicates the hot subdwarf formed during a common envelope phase when a $2.5-2.8$ M$_\odot$ star lost its envelope when crossing the Hertzsprung Gap. To match its current $P_{\rm orb}$, $T_{\rm eff}$, $\log(g)$, and masses, we estimate a post-common envelope period of $P_{\rm orb}\approx150$ min, and find the sdOB star is currently undergoing hydrogen shell burning. We estimate that the hot subdwarf will become a white dwarf with a thick helium layer of $\approx0.1$ M$_\odot$ and will merge with its carbon/oxygen white dwarf companion after $\approx17$ Myr and presumably explode as a thermonuclear supernova or form an R CrB star., Comment: accepted for publication in ApJ, 21 pages, 12 figures and 4 tables

We construct a class of 2+1 dimensional relativistic quantum field theories which exhibit the Fractional Quantum Hall Effect in the infrared, both in the continuum and on the lattice. The UV completion consists of a perturbative $U(1)\times U(1)$ gauge theory with only integer-charged fields, while the low energy spectrum consists of nontrivial topological phases supporting fractional currents, fractionally charged chiral surface modes and bulk anyonic excitations. Exotic phenomena such as a Fractional Quantum Spin Hall Effect can arise in such models.

The jet opening angle and inclination of GW170817 -- the first detected binary neutron star merger -- were vital to understand its energetics, relation to short gamma-ray bursts, and refinement of the standard siren-based determination of the Hubble constant, $H_0$. These basic quantities were determined through a combination of the radio lightcurve and Very Long Baseline Interferometry (VLBI) measurements of proper motion. In this paper we discuss and quantify the prospects for the use of radio VLBI observations and observations of scintillation-induced variability to measure the source size and proper motion of merger afterglows, and thereby infer properties of the merger including inclination angle, opening angle and energetics. We show that these techniques are complementary as they probe different parts of the circum-merger density/inclination angle parameter space and different periods of the temporal evolution of the afterglow. We also find that while VLBI observations will be limited to the very closest events it will be possible to detect scintillation for a large fraction of events beyond the range of current gravitational wave detectors. Scintillation will also be detectable with next generation telescopes such as the Square Kilometre Array, 2000 antenna Deep Synoptic Array and the next generation Very Large Array, for a large fraction of events detected with third generation gravitational wave detectors. Finally, we discuss prospects for the measurement of the $H_0$ with VLBI observations of neutron star mergers and compare this technique to other standard siren methods.

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $\kappa < 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ \chi < 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv., Comment: Accepted to ApJ with no major changes

We present results from a search for a radio transient associated with the LIGO/Virgo source S190814bv, a likely neutron star-black hole (NSBH) merger, with the Australian Square Kilometre Array Pathfinder. We imaged a $30\,{\rm deg}^2$ field at $\Delta T$=2, 9 and 33 days post-merger at a frequency of 944\,MHz, comparing them to reference images from the Rapid ASKAP Continuum Survey observed 110 days prior to the event. Each epoch of our observations covers $89\%$ of the LIGO/Virgo localisation region. We conducted an untargeted search for radio transients in this field, resulting in 21 candidates. For one of these, \object[AT2019osy]{AT2019osy}, we performed multi-wavelength follow-up and ultimately ruled out the association with S190814bv. All other candidates are likely unrelated variables, but we cannot conclusively rule them out. We discuss our results in the context of model predictions for radio emission from neutron star-black hole mergers and place constrains on the circum-merger density and inclination angle of the merger. This survey is simultaneously the first large-scale radio follow-up of an NSBH merger, and the most sensitive widefield radio transients search to-date.

We report the discovery of a detached double white dwarf binary with an orbital period of $\approx20.6\,\rm minutes$, PTF J053332.05+020911.6. The visible object in this binary, PTF J0533+0209B, is a $\approx 0.17\,M_\odot$ mass white dwarf with a helium-dominated atmosphere containing traces of hydrogen (DBA). This object exhibits ellipsoidal variations due to tidal deformation, and is the visible component in a single-lined spectroscopic binary with a velocity semi-amplitude of $K_B=618.7\pm6.9 \, \rm km \, s^{-1}$. We have detected significant orbital decay due to the emission of gravitational radiation, and we expect that the Laser Interferometer Space Antenna (LISA) will detect this system with a signal to noise of $8.4^{+4.2}_{-3.0}$ after four years of operation. Because this system already has a well determined orbital period, radial velocity semi-amplitude, temperature, atmospheric composition, surface gravity, and orbital decay rate, a LISA signal will help fully constrain the properties of this system by providing a direct measurement of its inclination. Thus, this binary demonstrates the synergy between electromagnetic and gravitational radiation for constraining the physical properties of an astrophysical object., Comment: 12 pages, 9 figures, 3 tables. Published online by the Astrophysical Journal Letters November 15, 2019

Slow-roll inflation is a successful paradigm. However we find that even a small coupling of the inflaton to other light fields can dramatically alter the dynamics and predictions of inflation. As an example, the inflaton can generically have an axion-like coupling to gauge bosons. Even relatively small couplings will automatically induce a thermal bath during inflation. The thermal friction from this bath can easily be stronger than Hubble friction, significantly altering the usual predictions of any particular inflaton potential. Thermal effects suppress the tensor-to-scalar ratio $r$ significantly, and predict unique non-gaussianities. This axion-like coupling provides a minimal model of warm inflation which avoids the usual problem of thermal backreaction on the inflaton potential. As a specific example, we find that hybrid inflation with this axion-like coupling can easily fit the current cosmological data., Comment: 18 pages, 1 figure, v2: We added additional references and clarifying comments in the introduction. We added an estimate on thermalization in section III, and an additional comment on cosine-like potentials in section IV, and a footnote commenting on equation 12. v2 matches published version. v3 corrects several minor mistakes in equations in Sec.4

In this paper we build upon recent work in the radio astronomy community to experimentally demonstrate the viability of passive radar for Space Situational Awareness. Furthermore, we show that the six state parameters of objects in orbit may be measured and used to perform orbit characterisation/estimation., Comment: Published in: 2017 IEEE Radar Conference (RadarConf) URL: http://ieeexplore.ieee.org/stamp/stamp.jsp?tp=&arnumber=7944483&isnumber=7944108

The Green Bank North Celestial Cap (GBNCC) pulsar survey will cover the entire northern sky ($\delta > -40^\circ$) at 350 MHz, and is one of the most uniform and sensitive all-sky pulsar surveys to date. We have created a pipeline to re-analyze GBNCC survey data to take a 350MHz census of all pulsars detected by the survey, regardless of their discovery survey. Of the 1413 pulsars in the survey region, we were able to recover 661. For these we present measured signal-to-noise ratios (S/N), flux densities, pulse widths, profiles, and where appropriate, refined dispersion measurements (647 out of 661) and new or improved spectral indices (276 out of 661 total, 15 new, 261 improved). Detection scans for several hundred sources were reanalyzed in order to inspect pulsars' single pulse behavior and 223 were confirmed to exhibit evidence of nulling. With a detailed analysis of measured and expected S/N values and the evolving radio frequency interference environment at 350MHz, we assess the GBNCC survey's sensitivity as a function of spin period, dispersion measure, and sky position. We find the sky-averaged limiting flux density of the survey to be 0.74mJy. Combining this analysis with PsrPopPy pulsar population simulations, we predict 60/5 non-recycled/millisecond pulsar discoveries in the survey's remaining 21,000 pointings, and we begin to place constraints on population model parameters.

We demonstrate that rotational superradiance can be efficient in millisecond pulsars. Measurements from the two fastest known pulsars PSR J1748-2446ad and PSR B1937+21 can place bounds on bosons with masses below 10^{-11} eV. The bounds are maximally good at masses corresponding to the rotation rate of the star, where scalar interactions that mediate forces ~ 10^6 times weaker than gravity are ruled out, exceeding existing fifth force constraints by 3 orders of magnitude. For certain neutron star equations of state, these measurements would also constrain the QCD axion with masses between 5 10^{-13} and 3 10^{-12} eV. Despite the ability of most neutron star equations of state to support frequencies as high as ~ 1500 Hz, the observed absence of pulsars above ~ 700 Hz could be due to the existence of a new particle of mass ~ 10^{-11} eV with a Yukawa coupling to nucleons., Comment: 31 pages, 4 figures

We identified a highly-polarized, steep-spectrum radio source in a deep image with the Australian Square Kilometre Array Pathfinder (ASKAP) telescope at 888 MHz. After considering and rejecting a stellar origin for this source, we discovered a new millisecond pulsar (MSP) using observations from the Parkes radio telescope. This pulsar has period 2.77 ms and dispersion measure 228.27 pc/cm**3. Although this pulsar does not yet appear to be particularly remarkable, the short spin period, wide profile and high dispersion measure do make it relatively hard to discover through traditional blind periodicity searches. Over the course of several weeks we see changes in the barycentric period of this pulsar that are consistent with orbital motion in a binary system, but the properties of any binary need to be confirmed by further observations. While even a deep ASKAP survey may not identify large numbers of new MSPs compared to the existing population, it would be competitive with existing all-sky surveys and could discover interesting new MSPs at high Galactic latitude without the need for computationally-expensive all-sky periodicity searches., Comment: ApJ, in press

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

Low mass, low metallicity galaxies at low to moderate ($z\lesssim3$) redshifts offer the best opportunity for detailed examination of the interplay between massive stars, ionizing radiation and gas in sources similar to those that likely reionized the universe. We present new narrowband Hubble Space Telescope observations of Ly$\alpha$ emission and the adjacent ultraviolet (UV) continuum in the low mass ($M_{\star} = 2 \times 10^8$ M$_{\odot}$), low metallicity ($Z\sim1/20$ Z$_{\odot}$) and highly ionized gravitationally lensed galaxy SL2S J02176$-$0513 at $z=1.844$. The galaxy has strong Ly$\alpha$ emission with photometric equivalent width $W^{\rm phot}_{\rm{Ly}\alpha} = 218 \pm 12$ \AA, at odds with the Ly$\alpha$ escape fraction of 10%. However, the spectroscopic Ly$\alpha$ profile suggests the presence of broad absorption underlying the emission, and the total equivalent width is consistent with the escape fraction once this underlying absorption is included. The Ly$\alpha$ emission is more spatially extended than the UV continuum, and the 0.14" spatial resolution of HST coupled with the magnification of gravitational lensing enables us to examine the distribution of Ly$\alpha$ and the UV continuum on sub-kiloparsec scales. We find that the peaks of the Ly$\alpha$ emission and the UV continuum are offset by 650 pc, and there is no Ly$\alpha$ emission arising from the region with the strongest UV light. Our combined spectroscopic and imaging data imply a significant range in neutral hydrogen column density across the object. These observations offer indirect support for a model in which ionizing radiation escapes from galaxies through channels with low column density of neutral gas., Comment: 13 pages, 7 figures. Revised version resubmitted to ApJ

General relativity predicts that short orbital period binaries emit significant gravitational radiation, and the upcoming Laser Interferometer Space Antenna (LISA) is expected to detect tens of thousands of such systems; however, few have been identified, and only one is eclipsing--the double white dwarf binary SDSS J065133.338+284423.37, which has an orbital period of 12.75 minutes. Here, we report the discovery of an eclipsing double white dwarf binary system with an orbital period of only 6.91 minutes, ZTF J153932.16+502738.8. This system has an orbital period close to half that of SDSS J065133.338+284423.37 and an orbit so compact that the entire binary could fit within the diameter of the planet Saturn. The system exhibits a deep eclipse, and a double-lined spectroscopic nature. We observe rapid orbital decay, consistent with that expected from general relativity. ZTF J153932.16+502738.8 is a significant source of gravitational radiation close to the peak of LISA's sensitivity, and should be detected within the first week of LISA observations., Comment: 44 pages, 10 figures, 2 tables. Published online by Nature on July 24, 2019