Search

JWST, despite not being designed to observe astrophysical phenomena that vary on rapid time scales, can be an unparalleled tool for such studies. If timing systematics can be controlled, JWST will be able to open up the sub-second infrared timescale regime. Rapid time-domain studies, such as lag measurements in accreting compact objects and Solar System stellar occultations, require both precise inter-frame timing and knowing when a time series begins to an absolute accuracy significantly below 1s. In this work we present two long-duration observations of the deeply eclipsing double white dwarf system ZTF J153932.16+502738.8, which we use as a natural timing calibrator to measure the absolute timing accuracy of JWST's clock. From our two epochs, we measure an average clock accuracy of $0.12\pm0.06$s, implying that JWST can be used for sub-second time-resolution studies down to the $\sim100$ms level, a factor $\sim5$ improvement upon the pre-launch clock accuracy requirement. We also find an asymmetric eclipse profile in the F322W2 band, which we suggest has a physical origin., Comment: 16 pages, 13 figures, accepted for publication in AJ

Detections of fast X-ray transients (FXTs) have been accrued over the last few decades. However, their origin has remained mysterious. There is now rapid progress thanks to timely discoveries and localisations with the Einstein Probe mission. Early results indicate that FXTs may frequently, but not always, be associated with gamma-ray bursts (GRBs). Here, we report on the multi-wavelength counterpart of FXT EP240414a, which has no reported gamma-ray counterpart. The transient is located 25.7~kpc in projection from a massive galaxy at $z=0.40$. We perform comprehensive photometric and spectroscopic follow-up. The optical light curve shows at least three distinct emission episodes with timescales of $\sim 1, 4$ and 15 days and peak absolute magnitudes of $M_R \sim -20$, $-21$, and $-19.5$, respectively. The optical spectrum at early times is extremely blue, inconsistent with afterglow emission. It may arise from the interaction of both jet and supernova shock waves with the stellar envelope and a dense circumstellar medium, as has been suggested for some Fast Blue Optical Transients (LFBOTs). At late times, the spectrum evolves to a broad-lined~Type~Ic supernova, similar to those seen in collapsar long-GRBs. This implies that the progenitor of EP240414a is a massive star creating a jet-forming supernova inside a dense envelope, resulting in an X-ray outburst with a luminosity of $\sim 10^{48}$ erg s$^{-1}$, and the complex observed optical/IR light curves. If correct, this argues for a causal link between the progenitors of long-GRBs, FXTs and LFBOTs., Comment: 36 pages, 13 figures, submitted to ApJ

Due to their short orbital periods and relatively high flux ratios, irradiated brown dwarfs in binaries with white dwarfs offer better opportunities to study irradiated atmospheres than hot Jupiters, which have lower planet-to-star flux ratios. WD1032+011 is an eclipsing, tidally locked white dwarf-brown dwarf binary with a 9950 K white dwarf orbited by a 69.7 M$_{Jup}$ brown dwarf in a 0.09 day orbit. We present time-resolved Hubble Space Telescope Wide Field Camera 3 spectrophotometric data of WD1032+011. We isolate the phase-dependent spectra of WD1032+011B, finding a 210 K difference in brightness temperature between the dayside and nightside. The spectral type of the brown dwarf is identified as L1 peculiar, with atmospheric retrievals and comparison to field brown dwarfs showing evidence for a cloud-free atmosphere. The retrieved temperature of the dayside is $1748^{+66}_{-67}$ K, with a nightside temperature of $1555^{+76}_{-62}$ K, showing an irradiation-driven temperature contrast coupled with inefficient heat redistribution from the dayside to the nightside. The brown dwarf radius is inflated, likely due to the constant irradiation from the white dwarf, making it the only known inflated brown dwarf in an eclipsing white dwarf-brown dwarf binary., Comment: 20 pages, 23 figures. Accepted for Publication in MNRAS

UW Coronae Borealis (UW CrB) is a low mass X-ray binary that shows both Type 1 X-ray and optical bursts, which typically last for 20 s. The system has a binary period of close to 2 hours and is thought to have a relatively high inclination due to the presence of an eclipse in the optical light curve. There is also evidence that an asymmetric disc is present in the system, which precesses every 5.5 days based on changes in the depth of the eclipse. In this paper, we present optical photometry and spectroscopy of UW CrB taken over 2 years. We update the orbital ephemeris using observed optical eclipses and refine the orbital period to 110.97680(1) min. A total of 17 new optical bursts are presented, with 10 of these bursts being resolved temporally. The average $e$-folding time of $19\pm3$s for the bursts is consistent with the previously found value. Optical bursts are observed during a previously identified gap in orbital phase centred on $\phi=0.967$, meaning the reprocessing site is not eclipsed as previously thought. Finally, we find that the apparent P-Cygni profiles present in some of the atomic lines in the optical spectra are due to transient absorption., Comment: 14 pages, 17 figures, submitted to OJAp

Time-domain astrophysics continues to grow rapidly, with the inception of new surveys drastically increasing data volumes. Democratised, distributed approaches to training sets for machine learning classifiers are crucial to make the most of this torrent of discovery -- with citizen science approaches proving effective at meeting these requirements. In this paper, we describe the creation of and the initial results from the $\textit{Kilonova Seekers}$ citizen science project, built to find transient phenomena from the GOTO telescopes in near real-time. $\textit{Kilonova Seekers}$ launched in July 2023 and received over 600,000 classifications from approximately 2,000 volunteers over the course of the LIGO-Virgo-KAGRA O4a observing run. During this time, the project has yielded 20 discoveries, generated a `gold-standard' training set of 17,682 detections for augmenting deep-learned classifiers, and measured the performance and biases of Zooniverse volunteers on real-bogus classification. This project will continue throughout the lifetime of GOTO, pushing candidates at ever-greater cadence, and directly facilitate the next-generation classification algorithms currently in development., Comment: 20 pages, 15 figures. Accepted in MNRAS

Long gamma-ray bursts (GRBs) are believed to originate from core collapse of massive stars. High-redshift GRBs can probe the star formation and reionization history of the early universe, but their detection remains rare. Here we report the detection of a GRB triggered in the 0.5--4 keV band by the Wide-field X-ray Telescope (WXT) on board the Einstein Probe (EP) mission, designated as EP240315a, whose bright peak was also detected by the Swift Burst Alert Telescope and Konus-Wind through off-line analyses. At a redshift of $z=4.859$, EP240315a showed a much longer and more complicated light curve in the soft X-ray band than in gamma-rays. Benefiting from a large field-of-view ($\sim$3600 deg$^2$) and a high sensitivity, EP-WXT captured the earlier engine activation and extended late engine activity through a continuous detection. With a peak X-ray flux at the faint end of previously known high-$z$ GRBs, the detection of EP240315a demonstrates the great potential for EP to study the early universe via GRBs., Comment: 41 pages, 8 figures, 7 tables

The nature of the minute-to-hour long Fast X-ray Transients (FXTs) localised by telescopes such as Chandra, Swift, and XMM-Newton remains mysterious, with numerous models suggested for the events. Here, we report multi-wavelength observations of EP240315a, a 1600 s long transient detected by the Einstein Probe, showing it to have a redshift of z=4.859. We measure a low column density of neutral hydrogen, indicating that the event is embedded in a low-density environment, further supported by direct detection of leaking ionising Lyman-continuum. The observed properties are consistent with EP240315a being a long-duration gamma-ray burst, and these observations support an interpretation in which a significant fraction of the FXT population are lower-luminosity examples of similar events. Such transients are detectable at high redshifts by the Einstein Probe and, in the (near) future, out to even larger distances by SVOM, THESEUS, and Athena, providing samples of events into the epoch of reionisation., Comment: 41 pages, 7 figures, submitted

A key feature of active galactic nuclei (AGN) is their variability across all wavelengths. Typically, AGN vary by a few tenths of a magnitude or more over periods lasting from hours to years. By contrast, extreme variability of AGN -- large luminosity changes that are a significant departure from the baseline variability -- are known as AGN flares. These events are rare and their timescales poorly constrained, with most of the literature focusing on individual events. It has been suggested that extreme AGN variability including flares can provide insights into the accretion processes in the disk. With surveys such as the Legacy Survey of Space and Time (LSST) promising millions of transient detections per night in the coming decade, there is a need for fast and efficient classification of AGN flares. The problem with the systematic detection of AGN flares is the requirement to detect them against a stochastically variable baseline; the ability to define a signal as a significant departure from the ever-present variability is a statistical challenge. Recently, Gaussian Processes (GPs) have revolutionised the analysis of time-series data in many areas of astronomical research. They have, however, seen limited uptake within the field of transient detection and classification. Here we investigate the efficacy of Gaussian Processes to detect AGN flares in both simulated and real optical light curves. We show that GP analysis can successfully detect AGN flares with a false-positive rate of less than seven per cent, and we present examples of AGN light curves that show extreme variability., Comment: 14 pages, 25 figures, accepted for publication in MNRAS

Radio emission has been detected from tens of white dwarfs, in particular in accreting systems. Additionally, radio emission has been predicted as a possible outcome of a planetary system around a white dwarf. We searched for 3 GHz radio continuum emission in 846,000 candidate white dwarfs previously identified in Gaia using the Very Large Array Sky Survey (VLASS) Epoch 1 Quick Look Catalogue. We identified 13 candidate white dwarfs with a counterpart in VLASS within 2". Five of those were found not to be white dwarfs in follow-up or archival spectroscopy, whereas seven others were found to be chance alignments with a background source in higher-resolution optical or radio images. The remaining source, WDJ204259.71+152108.06, is found to be a white dwarf and M-dwarf binary with an orbital period of 4.1 days and long-term stochastic optical variability, as well as luminous radio and X-ray emission. For this binary, we find no direct evidence of a background contaminant, and a chance alignment probability of only ~2 per cent. However, other evidence points to the possibility of an unfortunate chance alignment with a background radio and X-ray emitting quasar, including an unusually poor Gaia DR3 astrometric solution for this source. With at most one possible radio emitting white dwarf found, we conclude that strong (> 1-3 mJy) radio emission from white dwarfs in the 3 GHz band is virtually nonexistent outside of interacting binaries., Comment: 18 pages, 17 figures. Updated to match version accepted for publication in MNRAS

Approximately 150 low-mass white dwarfs, with masses below 0.4Msun, have been discovered. The majority of these low-mass WDs are observed in binary systems as they cannot be formed through single-star evolution within the Hubble time. In this study, we present a comprehensive analysis of the double low-mass WD eclipsing binary system J2102-4145. Our investigation involved an extensive observational campaign, resulting in the acquisition of approximately 28 hours of high-speed photometric data across multiple nights using NTT/ULTRACAM, SOAR/Goodman, and SMARTS-1m telescopes. These observations have provided critical insights into the orbital characteristics of this system, including parameters such as inclination and orbital period. To disentangle the binary components of J2102-4145, we employed the XT GRID spectral fitting method with GMOS/Gemini-South and X-Shooter data. Additionally, we used the PHOEBE package for light curve analysis on NTT/ULTRACAM high-speed time-series photometry data to constrain the binary star properties. Our analysis reveals remarkable similarities between the two components of this binary system. For the primary star, we determined Teff1 = 13688 +- 65 K, log g1 = 7.36 +- 0.01, R1 = 0.0211 +- 0.0002 Rsun, and M1 = 0.375 +- 0.003 Msun, while the secondary star is characterized by Teff2 = 12952 +- 53 K, log g2 = 7.32 +- 0.01, R2 = 0.0203 +- 0.0002 Rsun, and M2 = 0.31 +- 0.003 Msun. Furthermore, we observe a notable discrepancy between Teff and R of the less massive WD compared to evolutionary sequences for WDs from the literature, which has significant implications for our understanding of WD evolution. We discuss a potential formation scenario for this system that might explain this discrepancy and explore its future evolution. We predict that this system will merge in about 800 Myr, evolving into a helium-rich hot subdwarf star and later into a hybrid He/CO WD.

A large number of white dwarf and donor masses in cataclysmic variables have been found via modelling the primary eclipse, a method that relies on untested assumptions. Recent measurements of the mass of the white dwarf in the cataclysmic variable GY Cnc, obtained via modelling its ultraviolet spectrum, conflict with the mass obtained via modelling the eclipse light curve. Here we measure the radial velocity of the absorption lines from the donor star in GY Cnc to be $K_{\rm abs} = 280 \pm 2$ kms$^{-1}$, in excellent agreement with the prediction based on the masses derived from modelling the eclipse light curve. It is possible that the white dwarf mass derived from the ultraviolet spectrum of GY Cnc is affected by the difficulty of disentangling the white dwarf spectrum from the accretion disc spectrum., Comment: 7 pages, 6 figures, accepted for Publication in Monthly Notices of the Royal Astronomical Society

AM CVn-type systems are ultracompact, helium-accreting binary systems which are evolutionarily linked to the progenitors of thermonuclear supernovae and are expected to be strong Galactic sources of gravitational waves detectable to upcoming space-based interferometers. AM CVn binaries with orbital periods $\lesssim$ 20--23 min exist in a constant high state with a permanently ionised accretion disc. We present the discovery of TIC 378898110, a bright ($G=14.3$ mag), nearby ($309.3 \pm 1.8$ pc), high-state AM CVn binary discovered in TESS two-minute-cadence photometry. At optical wavelengths this is the third-brightest AM CVn binary known. The photometry of the system shows a 23.07172(6) min periodicity, which is likely to be the `superhump' period and implies an orbital period in the range 22--23 min. There is no detectable spectroscopic variability. The system underwent an unusual, year-long brightening event during which the dominant photometric period changed to a shorter period (constrained to $20.5 \pm 2.0$ min), which we suggest may be evidence for the onset of disc-edge eclipses. The estimated mass transfer rate, $\log (\dot{M} / \mathrm{M_\odot} \mathrm{yr}^{-1}) = -6.8 \pm 1.0$, is unusually high and may suggest a high-mass or thermally inflated donor. The binary is detected as an X-ray source, with a flux of $9.2 ^{+4.2}_{-1.8} \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 0.3--10 keV range. TIC 378898110 is the shortest-period binary system discovered with TESS, and its large predicted gravitational-wave amplitude makes it a compelling verification binary for future space-based gravitational wave detectors., Comment: 15 pages, 14 figures. Accepted to MNRAS

We present the results of our analysis of Gaia19dke, an extraordinary microlensing event in the Cygnus constellation that was first spotted by the {\gaia} satellite. This event featured a strong microlensing parallax effect, which resulted in multiple peaks in the light curve. We conducted extensive photometric, spectroscopic, and high-resolution imaging follow-up observations to determine the mass and the nature of the invisible lensing object. Using the Milky Way priors on density and velocity of lenses, we found that the dark lens is likely to be located at a distance of $D_L =(3.05^{+4.10}_{-2.42})$kpc, and has a mass of $M_L =(0.51^{+3.07}_{-0.40}) M_\odot$. Based on its low luminosity and mass, we propose that the lens in Gaia19dke event is an isolated white dwarf., Comment: 16 pages, submitted to A&A

White dwarfs, the extremely dense remnants left behind by most stars after their death, are characterised by a mass comparable to that of the Sun compressed into the size of an Earth-like planet. In the resulting strong gravity, heavy elements sink toward the centre and the upper layer of the atmosphere contains only the lightest element present, usually hydrogen or helium. Several mechanisms compete with gravitational settling to change a white dwarf's surface composition as it cools, and the fraction of white dwarfs with helium atmospheres is known to increase by a factor ~2.5 below a temperature of about 30,000 K; therefore, some white dwarfs that appear to have hydrogen-dominated atmospheres above 30,000 K are bound to transition to be helium-dominated as they cool below it. Here we report observations of ZTF J203349.8+322901.1, a transitioning white dwarf with two faces: one side of its atmosphere is dominated by hydrogen and the other one by helium. This peculiar nature is likely caused by the presence of a small magnetic field, which creates an inhomogeneity in temperature, pressure or mixing strength over the surface. ZTF J203349.8+322901.1 might be the most extreme member of a class of magnetic, transitioning white dwarfs -- together with GD 323, a white dwarf that shows similar but much more subtle variations. This new class could help shed light on the physical mechanisms behind white dwarf spectral evolution., Comment: 45 pages, 11 figures

We present the discovery of the eclipsing double white dwarf (WD) binary WDJ 022558.21-692025.38 that has an orbital period of 47.19 min. Following identification with the Transiting Exoplanet Survey Satellite, we obtained time-series ground based spectroscopy and high-speed multi-band ULTRACAM photometry which indicate a primary DA WD of mass 0.40 +- 0.04 Msol and a 0.28 +- 0.02 Msol mass secondary WD, which is likely of type DA as well. The system becomes the third-closest eclipsing double WD binary discovered with a distance of approximately 400 pc and will be a detectable source for upcoming gravitational wave detectors in the mHz frequency range. Its orbital decay will be measurable photometrically within 10 yrs to a precision of better than 1%. The fate of the binary is to merge in approximately 41 Myr, likely forming a single, more massive WD., Comment: Accepted for publication in MNRAS, 8 pages + 2 appendix pages, 6 figures

This paper reports the ULTRACAM discovery of dipolar surface spots in two cool magnetic white dwarfs with Balmer emission lines, while a third system exhibits a single spot, similar to the prototype GD 356. The light curves are modeled with simple, circular, isothermal dark spots, yielding relatively large regions with minimum angular radii of 20 deg. For those stars with two light curve minima, the dual spots are likely observed at high inclination (or colatitude), however, identical and antipodal spots cannot simultaneously reproduce both the distinct minima depths and the phases of the light curve maxima. The amplitudes of the multi-band photometric variability reported here are all several times larger than that observed in the prototype GD 356; nevertheless, all DAHe stars with available data appear to have light curve amplitudes that increase toward the blue in correlated ratios. This behavior is consistent with cool spots that produce higher contrasts at shorter wavelengths, with remarkably similar spectral properties given the diversity of magnetic field strengths and rotation rates. These findings support the interpretation that some magnetic white dwarfs generate intrinsic chromospheres as they cool, and that no external source is responsible for the observed temperature inversion. Spectroscopic time-series data for DAHe stars is paramount for further characterization, where it is important to obtain well-sampled data, and consider wavelength shifts, equivalent widths, and spectropolarimetry., Comment: 9 pages, 5 figures, 3 tables, accepted to MNRAS

The mergers of binary compact objects such as neutron stars and black holes are of central interest to several areas of astrophysics, including as the progenitors of gamma-ray bursts (GRBs), sources of high-frequency gravitational waves and likely production sites for heavy element nucleosynthesis via rapid neutron capture (the r-process). These heavy elements include some of great geophysical, biological and cultural importance, such as thorium, iodine and gold. Here we present observations of the exceptionally bright gamma-ray burst GRB 230307A. We show that GRB 230307A belongs to the class of long-duration gamma-ray bursts associated with compact object mergers, and contains a kilonova similar to AT2017gfo, associated with the gravitational-wave merger GW170817. We obtained James Webb Space Telescope mid-infrared (mid-IR) imaging and spectroscopy 29 and 61 days after the burst. The spectroscopy shows an emission line at 2.15 microns which we interpret as tellurium (atomic mass A=130), and a very red source, emitting most of its light in the mid-IR due to the production of lanthanides. These observations demonstrate that nucleosynthesis in GRBs can create r-process elements across a broad atomic mass range and play a central role in heavy element nucleosynthesis across the Universe., Comment: Submitted. Comments welcome! Nature (2023)

White dwarf stars are the most common stellar fossils. When in binaries, they make up the dominant form of compact object binary within the Galaxy and can offer insight into different aspects of binary formation and evolution. One of the most remarkable white dwarf binary systems identified to date is AR Scorpii (henceforth AR Sco). AR Sco is composed of an M-dwarf star and a rapidly-spinning white dwarf in a 3.56-hour orbit. It shows pulsed emission with a period of 1.97 minutes over a broad range of wavelengths, which led to it being known as a white dwarf pulsar. Both the pulse mechanism and the evolutionary origin of AR Sco provide challenges to theoretical models. Here we report the discovery of the first sibling of AR Sco, J191213.72-441045.1 (henceforth J1912-4410), which harbours a white dwarf in a 4.03-hour orbit with an M-dwarf and exhibits pulsed emission with a period of 5.30 minutes. This discovery establishes binary white dwarf pulsars as a class and provides support for proposed formation models for white dwarf pulsars., Comment: Authors' version of article published in Nature Astronomy (DOI 10.1038/s41550-023-01995-x)

We present a catalogue of 507 cataclysmic variables (CVs) observed in SDSS I to IV including 70 new classifications collated from multiple archival data sets. This represents the largest sample of CVs with high-quality and homogeneous optical spectroscopy. We have used this sample to derive unbiased space densities and period distributions for the major sub-types of CVs. We also report on some peculiar CVs, period bouncers and also CVs exhibiting large changes in accretion rates. We report 70 new CVs, 59 new periods, 178 unpublished spectra and 262 new or updated classifications. From the SDSS spectroscopy, we also identified 18 systems incorrectly identified as CVs in the literature. We discuss the observed properties of 13 peculiar CVS, and we identify a small set of eight CVs that defy the standard classification scheme. We use this sample to investigate the distribution of different CV sub-types, and we estimate their individual space densities, as well as that of the entire CV population. The SDSS I to IV sample includes 14 period bounce CVs or candidates. We discuss the variability of CVs across the Hertzsprung-Russell diagram, highlighting selection biases of variability-based CV detection. Finally, we searched for, and found eight tertiary companions to the SDSS CVs. We anticipate that this catalogue and the extensive material included in the Supplementary Data will be useful for a range of observational population studies of CVs., Comment: Published by MNRAS (including 2 corrections published on 7 December 2023). Includes supplementary material including machine readable version

Fast X-ray Transients (FXTs) are X-ray flares with a duration ranging from a few hundred seconds to a few hours. Possible origins include the tidal disruption of a white dwarf by an intermediate-mass black hole, a supernova shock breakout, and a binary neutron star merger. We present the X-ray light curve and spectrum, and deep optical imaging of the FXT XRT 210423, which has been suggested to be powered by a magnetar produced in a binary neutron star merger. Our Very Large Telescope and Gran Telescopio Canarias (GTC) observations began on May 6, 2021, thirteen days after the onset of the flare. No transient optical counterpart is found in the 1" (3$\sigma$) X-ray uncertainty region of the source to a depth $g_{s}$=27.0 AB mag. A candidate host lies within the 1" X-ray uncertainty region with a magnitude of 25.9 $\pm$ 0.1 in the GTC/HiPERCAM $g_s$-filter. Due to its faintness, it was not detected in other bands, precluding a photometric redshift determination. We detect two additional candidate host galaxies; one with $z_{\rm spec}=1.5082 \pm 0.0001$ and an offset of 4.2$\pm$1" (37$\pm$9 kpc) from the FXT and another one with $z_{\rm phot}=1.04^{+0.22}_{-0.14}$, at an offset of 3.6$\pm$1" (30$\pm$8 kpc). Based on the properties of all the prospective hosts we favour a binary neutron star merger, as previously suggested in the literature, as explanation for XRT 210423., Comment: 14 pages, 6 figures. Accepted for publication in ApJ

Wide-field time-domain photometric sky surveys are now finding hundreds of eclipsing white dwarf plus M dwarf binaries, a population encompassing a wealth of information and potential insight into white dwarf and close binary astrophysics. Precise follow-up observations are essential in order to fully constrain these systems and capitalise on the power of this sample. We present the first results from our program of high-speed, multi-band photometric follow-up. We develop a method to measure temperatures, (model-dependent) masses, and radii for both components from the eclipse photometry alone and characterize 34 white dwarf binaries, finding general agreement with independent estimates using an alternative approach while achieving around a factor of two increase in parameter precision. In addition to these parameter estimates, we discover a number of interesting systems -- finding four with sub-stellar secondaries, doubling the number of eclipsing examples, and at least six where we find the white dwarf to be strongly magnetic, making these the first eclipsing examples of such systems and key to investigating the mechanism of magnetic field generation in white dwarfs. We also discover the first two pulsating white dwarfs in detached and eclipsing post-common-envelope binaries -- one with a low-mass, likely helium core, and one with a relatively high mass, towards the upper end of the known sample of ZZ Cetis. Our results demonstrate the power of eclipse photometry, not only as a method of characterising the population, but as a way of discovering important systems that would have otherwise been missed by spectroscopic follow-up., Comment: 12 pages with a 5 page appendix and 14 figures. Accepted for publication in MNRAS

Black widows are extreme millisecond pulsar binaries where the pulsar wind ablates their low-mass companion stars. Their optical light curves vary periodically due to the high irradiation and tidal distortion of the companion, which allows us to infer the binary parameters. We present simultaneous multi-band observations obtained with the HIPERCAM instrument at the 10.4-m GTC telescope for six of these systems. The combination of this five-band fast photometer with the world's largest optical telescope enables us to inspect the light curve range near minima. We present the first light curve for PSR J1641+8049, as well as attain a significant increase in signal-to-noise and cadence compared with previous publications for the remaining 5 targets: PSR J0023+0923, PSR J0251+2606, PSR J0636+5129, PSR J0952-0607 and PSR J1544+4937. We report on the results of the light curve modelling with the Icarus code for all six systems, which reveals some of the hottest and densest companion stars known. We compare the parameters derived with the limited but steadily growing black widow population for which optical modelling is available. We find some expected correlations, such as that between the companion star mean density and the orbital period of the system, but also a puzzling positive correlation between the orbital inclination and the irradiation temperature of the companion. We propose such a correlation would arise if pulsars with magnetic axis orthogonal to their spin axis are capable of irradiating their companions to a higher degree., Comment: 18 pages (+12 pages for appendix), 12 figures (+13 in the appendix), 3 tables (1 in the appendix). Accepted for publication in MNRAS on 2023 January 17th

We present a rapid timing analysis of optical (HiPERCAM and ULTRACAM) and X-ray (NICER) observations of the X-ray transient Swift J1858.6-0814 during 2018 and 2019. The optical light curves show relatively slow, large amplitude (~1 mags in g$_s$) `blue' flares (i.e. stronger at shorter wavelengths) on time-scales of ~minutes as well as fast, small amplitude (~0.1 mag in g$_s$) `red' flares (i.e. stronger at longer wavelengths) on time-scales of ~seconds. The `blue' and `red' flares are consistent with X-ray reprocessing and optically thin synchrotron emission, respectively, similar to what is observed in other X-ray binaries. The simultaneous optical versus soft- and hard-band X-ray light curves show time- and energy dependent correlations. The 2019 March 4 and parts of the June data show a nearly symmetric positive cross correlations (CCFs) at positive lags consistent with simple X-ray disc reprocessing. The soft- and hard-band CCFs are similar and can be reproduced if disc reprocessing dominates in the optical and one component (disc or synchrotron Comptonization) dominates both the soft and hard X-rays. A part of the 2019 June data shows a very different CCFs. The observed positive correlation at negative lag in the soft-band can be reproduced if the optical synchrotron emission is correlated with the hot flow X-ray emission. The observed timing properties are in qualitative agreement with the hybrid inner hot accretion flow model, where the relative role of the different X-ray and optical components that vary during the course of the outburst, as well as on shorter time-scales, govern the shape of the optical/X-ray CCFs., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 18 pages, 14 figures

We present fast (~200s sampling) ugriz photometry of the low mass AGN NGC 4395 with the Liverpool Telescope, followed by very fast (3s sampling) us, gs, rs, is and zs simultaneous monitoring with HiPERCAM on the 10.4m GTC. These observations provide the fastest ever AGN multiband photometry and very precise lag measurements. Unlike in all other AGN, gs lags us by a large amount, consistent with disc reprocessing but not with reprocessing in the Broad Line Region (BLR). There is very little increase in lag with wavelength at long wavelengths, indicating an outer edge (Rout) to the reprocessor. We have compared truncated disc reprocessing models to the combined HiPERCAM and previous X-ray/UV lags. For the normally accepted mass of 3.6E5 solar, we obtain reasonable agreement with zero spin, Rout ~1700 Rg, and the DONE physically-motivated temperature-dependent disc colour correction factor (fcol). A smaller mass of 4E4 solar can only be accomodated if fcol=2.4, which is probably unrealistically high. Disc self gravity is probably unimportant in this low mass AGN but an obscuring wind may provide an edge. For the small mass the dust sublimation radius is similar to Rout, so the wind could be dusty. However for the more likely large mass the sublimation radius is further out so the optically-thick base of a line-driven gaseous wind is more likely. The inner edge of the BLR is close to Rout in both cases. These observations provide the first good evidence for a truncated AGN disc and caution that truncation should be included in reverberation lag modelling., Comment: 19 pages, 27 figures

The shortest-period binary star system known to date, RX J0806.3+1527 (HM Cancri), has now been observed in the optical for more than two decades. Although it is thought to be a double degenerate binary undergoing mass transfer, an early surprise was that its orbital frequency, $f_0$, is currently increasing as the result of gravitational wave radiation. This is unusual since it was expected that the mass donor was degenerate and would expand on mass loss, leading to a decreasing $f_0$. We exploit two decades of high-speed photometry to precisely quantify the trajectory of HM Cancri, allowing us to find that $\ddot f_0$ is negative, where $\ddot f_0~=~(-5.38\pm2.10)\times10^{-27}$ Hz s$^{-2}$. Coupled with our positive frequency derivative, we show that mass transfer is counteracting gravitational-wave dominated orbital decay and that HM Cancri will turn around within $2100\pm800\,$yrs from now. We present Hubble Space Telescope ultra-violet spectra which display Lyman-$\alpha$ absorption, indicative of the presence of hydrogen accreted from the donor star. We use these pieces of information to explore a grid of permitted donor and accretor masses with the Modules for Experiments in Stellar Astrophysics suite, finding models in good accordance with many of the observed properties for a cool and initially hydrogen-rich extremely-low-mass white dwarf ($\approx0.17\,$M$_\odot$) coupled with a high accretor mass white dwarf ($\approx 1.0\,$M$_\odot$). Our measurements and models affirm that HM~Cancri is still one of the brightest verification binaries for the Laser Interferometer Space Antenna spacecraft., Comment: 12 pages (+5 pages appendix), 9 figures, 6 tables. Accepted for publication in MNRAS

In over a thousand known cataclysmic variables (CVs), where a white dwarf is accreting from a hydrogen-rich star, only a dozen have orbital periods below 75 minutes. One way to achieve these short periods requires the donor star to have undergone substantial nuclear evolution prior to interacting with the white dwarf, and it is expected that these objects will transition to helium accretion. These transitional CVs have been proposed as progenitors of helium CVs. However, no known transitional CV is expected to reach an orbital period short enough to account for most of the helium CV population, leaving the role of this evolutionary pathway unclear. Here we report observations of ZTF J1813+4251, a 51-minute orbital period, fully eclipsing binary system consisting of a star with a temperature comparable to that of the Sun but a density 100 times greater due to its helium-rich composition, accreting onto a white dwarf. Phase-resolved spectra, multi-band light curves and the broadband spectral energy distribution allow us to obtain precise and robust constraints on the masses, radii and temperatures of both components. Evolutionary modeling shows that ZTF J1813+4251 is destined to become a helium CV binary, reaching an orbital period under 20 minutes, rendering ZTF J1813+4251 a previously missing link between helium CV binaries and hydrogen-rich CVs., Comment: 48 Pages, 12 figures, 2 tables, Published online by Nature on Oct 5, 2022

Developing an effective automatic classifier to separate genuine sources from artifacts is essential for transient follow-ups in wide-field optical surveys. The identification of transient detections from the subtraction artifacts after the image differencing process is a key step in such classifiers, known as real-bogus classification problem. We apply a self-supervised machine learning model, the deep-embedded self-organizing map (DESOM) to this "real-bogus" classification problem. DESOM combines an autoencoder and a self-organizing map to perform clustering in order to distinguish between real and bogus detections, based on their dimensionality-reduced representations. We use 32x32 normalized detection thumbnails as the input of DESOM. We demonstrate different model training approaches, and find that our best DESOM classifier shows a missed detection rate of 6.6% with a false positive rate of 1.5%. DESOM offers a more nuanced way to fine-tune the decision boundary identifying likely real detections when used in combination with other types of classifiers, for example built on neural networks or decision trees. We also discuss other potential usages of DESOM and its limitations.

We present simultaneous, multi-colour optical light curves of the companion star to the black-widow pulsar PSR J2051-0827, obtained approximately 10 years apart using ULTRACAM and HiPERCAM, respectively. The ULTRACAM light curves confirm the previously reported asymmetry in which the leading hemisphere of the companion star appears to be brighter than the trailing hemisphere. The HiPERCAM light curves, however, do not show this asymmetry, demonstrating that whatever mechanism is responsible for it varies on timescales of a decade or less. We fit the symmetrical HiPERCAM light curves with a direct-heating model to derive the system parameters, finding an orbital inclination of $55.9^{+4.8}_{-4.1}$ degrees, in good agreement with radio-eclipse constraints. We find that approximately half of the pulsar's spin-down energy is converted to optical luminosity, resulting in temperatures ranging from approximately $5150^{+190}_{-190}$ K on the day side to $2750^{+130}_{-150}$ K on the night side of the companion star. The companion star is close to filling its Roche lobe ($f_{\rm RL} =0.88^{+0.02}_{-0.02}$) and has a mass of $0.039^{+0.010}_{-0.011}$ M$_{\odot}$, giving a mean density of $20.24^{+0.59}_{-0.44}$ g cm$^{-3}$ and an apsidal motion constant in the range $0.0036 < k_2 < 0.0047$. The companion mass and mean density values are consistent with those of brown dwarfs, but the apsidal motion constant implies a significantly more centrally-condensed internal structure than is typical for such objects., Comment: 10 pages, 4 figures, 3 tables, accepted for publication in Monthly Notices of the Royal Astronomical Society

AR Scorpii (AR Sco) is the only radio-pulsing white dwarf known to date. It shows a broad-band spectrum extending from radio to X-rays whose luminosity cannot be explained by thermal emission from the system components alone, and is instead explained through synchrotron emission powered by the spin-down of the white dwarf. We analysed NTT/ULTRACAM, TNT/ULTRASPEC, and GTC/HiPERCAM high-speed photometric data for AR Sco spanning almost seven years and obtained a precise estimate of the spin frequency derivative, now confirmed with 50-sigma significance. Using archival photometry, we show that the spin down rate of P/Pdot = 5.6e6 years has remained constant since 2005. As well as employing the method of pulse-arrival time fitting used for previous estimates, we also found a consistent value via traditional Fourier analysis for the first time. In addition, we obtained optical time-resolved spectra with WHT/ISIS and VLT/X-shooter. We performed modulated Doppler tomography for the first time for the system, finding evidence of emission modulated on the orbital period. We have also estimated the projected rotational velocity of the M-dwarf as a function of orbital period and found that it must be close to Roche lobe filling. Our findings provide further constraints for modelling this unique system., Comment: 15 pages, 18 figures. Accepted for publication in MNRAS

We present a detailed time-resolved photometric study of the ultra-compact X-ray binary candidate 4U 1812-12. The multicolor light curves obtained with HiPERCAM on the 10.4-m Gran Telescopio Canarias show an aprox 114 min modulation similar to a superhump. Under this interpretation, this period should lie very close to the orbital period of the system. Contrary to what its other observational properties suggest (namely, persistent dim luminosity, low optical-to-X-ray flux ratio and lack of hydrogen features in the optical spectrum), this implies that 4U1812-12 is most likely not an ultra-compact X-ray binary, which are usually defined as systems with orbital periods lower than 80 min. We discuss the nature of the system, showing that a scenario in which 4U 1812-12 is the progenitor of an ultra-compact X-ray binary may reconcile all the observables., Comment: Accepted for publication in ApJ Letters

Over a dozen millisecond pulsars are ablating low-mass companions in close binary systems. In the original "black widow", the 8-hour orbital period eclipsing pulsar PSR J1959+2048 (PSR B1957+20), high energy emission originating from the pulsar is irradiating and may eventually destroy a low-mass companion. These systems are not only physical laboratories that reveal the dramatic result of exposing a close companion star to the relativistic energy output of a pulsar, but are also believed to harbour some of the most massive neutron stars, allowing for robust tests of the neutron star equation of state. Here, we report observations of ZTF J1406+1222, a wide hierarchical triple hosting a 62-minute orbital period black widow candidate whose optical flux varies by a factor of more than 10. ZTF J1406+1222 pushes the boundaries of evolutionary models, falling below the 80 minute minimum orbital period of hydrogen-rich systems. The wide tertiary companion is a rare low metallicity cool subdwarf star, and the system has a Galactic halo orbit consistent with passing near the Galactic center, making it a probe of formation channels, neutron star kick physics, and binary evolution., Comment: 63 pages, 15 figures, 2 tables, published in Nature on May 4, 2022

Over the past few years, $\sim$30 extragalactic fast X-ray transients (FXRTs) have been discovered, mainly in Chandra and XMM-Newton data. Their nature remains unclear, with proposed origins including a double neutron star merger, a tidal disruption event involving an intermediate-mass black hole and a white dwarf, or a supernova shock breakout. A decisive differentiation between these three promising mechanisms for their origin requires an understanding of the FXRT energetics, environments, and/or host properties. We present optical observations obtained with the Very Large Telescope for the FXRTs XRT 000519 and XRT 110103 and Gran Telescopio Canarias observations for XRT 000519 designed to search for host galaxies of these FXRTs. In the $g_s$, $r_s$ and $R$-band images, we detect an extended source on the North-West side of the $\sim$ $1^{\prime\prime}$ (68% confidence) error circle of the X-ray position of XRT 000519 with a Kron magnitude of $g_s=$26.29$\pm$0.09 (AB magnitude). We discuss the XRT 000519 association with the probable host candidate for various possible distances, and we conclude that if XRT 000519 is associated with the host candidate a supernova shock breakout scenario is likely excluded. No host galaxy is found near XRT 110103 down to a limiting magnitude of $R>25.8$., Comment: 11 pages, 11 figures, submitted to MNRAS

With the prevalence of wide-field, time-domain photometric sky surveys, the number of eclipsing white dwarf systems being discovered is increasing dramatically. An efficient method to follow these up will be key to determining any population trends and finding any particularly interesting examples. We demonstrate that multi-band eclipse photometry of binaries containing a white dwarf and an M~dwarf can be used to determine the masses and temperatures of the white dwarfs to better than 5 per cent. For the M~dwarfs we measure their parameters to a precision of better than 6 per cent with the uncertainty dominated by the intrinsic scatter of the M~dwarf mass-radius relationship. This precision is better than what can typically be achieved with low-resolution spectroscopy. The nature of this method means that it will be applicable to LSST data in the future, enabling direct characterisation without follow-up spectroscopy. Additionally, we characterise three new post-common-envelope binaries from their eclipse photometry, finding two systems containing hot helium-core white dwarfs with low-mass companions (one near the brown dwarf transition regime) and a possible detached cataclysmic variable at the lower edge of the period gap., Comment: 10 pages (and appendix), 5 figures, accepted to MNRAS

We present a detailed study of the stellar and orbital parameters of the post-common envelope binary central star of the planetary nebula Ou~5. Low-resolution spectra obtained during the primary eclipse -- to our knowledge the first isolated spectra of the companion to a post-common-envelope planetary nebula central star -- were compared to catalogue spectra, indicating that the companion star is a late K- or early M-type dwarf. Simultaneous modelling of multi-band photometry and time-resolved radial velocity measurements was then used to independently determine the parameters of both stars as well as the orbital period and inclination. The modelling indicates that the companion star is low mass ($\sim$0.25~M$_\odot$) and has a radius significantly larger than would be expected for its mass. Furthermore, the effective temperature and surface gravity of nebular progenitor, as derived by the modelling, do not lie on single-star post-AGB evolutionary tracks, instead being more consistent with a post-RGB evolution. However, an accurate determination of the component masses is challenging. This is principally due to the uncertainty on the locus of the spectral lines generated by the irradiation of the companion's atmosphere by the hot primary (used to derive companion star's radial velocities), as well as the lack of radial velocities of the primary., Comment: 9 pages, 6 figures, MNRAS

A promising progenitor scenario for Type Ia supernovae (SNeIa) is the thermonuclear detonation of a white dwarf in a close binary system with another white dwarf. After the primary star explodes, the surviving donor can be spontaneously released as a hypervelocity runaway. One such runaway donor candidate is LP 398-9, whose orbital trajectory traces back $\approx 10^5$ years to a known supernova remnant. Here we report the discovery of carbon-rich circumstellar material around LP 398-9, revealed by a strong infrared excess and analyzed with follow-up spectroscopy. The circumstellar material is most plausibly composed of inflated layers from the star itself, mechanically and radioactively heated by the past companion's supernova. We also detect a 15.4 hr periodic signal in the UV and optical light curves of LP 398-9, which we interpret as surface rotation. The rotation rate is consistent with theoretical predictions from this supernova mechanism, and the brightness variations could originate from surface inhomogeneity deposited by the supernova itself. Our observations strengthen the case for this double-degenerate SNIa progenitor channel, and motivate the search for more runaway SNIa donors., Comment: Accepted to MNRAS

The Gravitational-wave Optical Transient Observer (GOTO) is an array of wide-field optical telescopes, designed to exploit new discoveries from the next generation of gravitational wave detectors (LIGO, Virgo, KAGRA), study rapidly evolving transients, and exploit multi-messenger opportunities arising from neutrino and very high energy gamma-ray triggers. In addition to a rapid response mode, the array will also perform a sensitive, all-sky transient survey with few day cadence. The facility features a novel, modular design with multiple 40-cm wide-field reflectors on a single mount. In June 2017 the GOTO collaboration deployed the initial project prototype, with 4 telescope units, at the Roque de los Muchachos Observatory (ORM), La Palma, Canary Islands. Here we describe the deployment, commissioning, and performance of the prototype hardware, and discuss the impact of these findings on the final GOTO design. We also offer an initial assessment of the science prospects for the full GOTO facility that employs 32 telescope units across two sites., Comment: 19 pages, 16 Figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Accreting magnetic white dwarfs offer an opportunity to understand the interplay between spin-up and spin-down torques in binary systems. Monitoring of the white dwarf spin may reveal whether the white dwarf spin is currently in a state of near-equilibrium, or of uni-directional evolution towards longer or shorter periods, reflecting the recent history of the system and providing constraints for evolutionary models. This makes the monitoring of the spin history of magnetic white dwarfs of high interest. In this paper we report the results of a campaign of follow-up optical photometry to detect and track the 39 sec white dwarf spin pulses recently discovered in Hubble Space Telescope data of the cataclysmic variable V1460 Her. We find the spin pulsations to be present in g-band photometry at a typical amplitude of 0.4%. Under favourable observing conditions, the spin signal is detectable using 2-meter class telescopes. We measured pulse-arrival times for all our observations, which allowed us to derive a precise ephemeris for the white dwarf spin. We have also derived an orbital modulation correction that can be applied to the measurements. With our limited baseline of just over four years, we detect no evidence yet for spin-up or spin-down of the white dwarf, obtaining a lower limit of |P/Pdot|> 4e7 years, which is already 4 to 8 times longer than the timescales measured in two other cataclysmic variable systems containing rapidly rotating white dwarfs, AE Aqr and AR Sco., Comment: 8 pages, 6 figures. Accepted for publication in MNRAS

The typical detection rate of $\sim1$ gamma-ray burst (GRB) per day by the \emph{Fermi} Gamma-ray Burst Monitor (GBM) provides a valuable opportunity to further our understanding of GRB physics. However, the large uncertainty of the \emph{Fermi} localization typically prevents rapid identification of multi-wavelength counterparts. We report the follow-up of 93 \emph{Fermi} GRBs with the Gravitational-wave Optical Transient Observer (GOTO) prototype on La Palma. We selected 53 events (based on favourable observing conditions) for detailed analysis, and to demonstrate our strategy of searching for optical counterparts. We apply a filtering process consisting of both automated and manual steps to 60\,085 candidates initially, rejecting all but 29, arising from 15 events. With $\approx3$ GRB afterglows expected to be detectable with GOTO from our sample, most of the candidates are unlikely to be related to the GRBs. Since we did not have multiple observations for those candidates, we cannot confidently confirm the association between the transients and the GRBs. Our results show that GOTO can effectively search for GRB optical counterparts thanks to its large field of view of $\approx40$ square degrees and its depth of $\approx20$ mag. We also detail several methods to improve our overall performance for future follow-up programs of \emph{Fermi} GRBs.

We present optical photometry of the cataclysmic variable LAMOST J024048.51+195226.9 taken with the high-speed, five-band CCD camera HiPERCAM on the 10.4 m Gran Telescopio Canarias (GTC). We detect pulsations originating from the spin of its white dwarf, finding a spin period of 24.9328(38)s. The pulse amplitude is of the order of 0.2% in the g-band, below the detection limits of previous searches. This detection establishes LAMOST J024048.51+195226.9 as only the second white dwarf magnetic propeller system, a twin of its long-known predecessor, AE Aquarii. At 24.93s, the white dwarf in LAMOST J024048.51+195226.9 has the shortest known spin period of any cataclysmic variable star. The white dwarf must have a mass of at least 0.7MSun to sustain so short a period. The observed faintest u-band magnitude sets an upper limit on the white dwarf's temperature of ~25000K. The pulsation amplitudes measured in the five HiPERCAM filters are consistent with an accretion spot of ~30000K covering ~2% of the white dwarf's visible area, although spots that are hot and smaller, or cooler and larger cannot be ruled out., Comment: 5 pages, 5 figures. Updated to match the version published in MNRAS Letters

HiPERCAM is a portable, quintuple-beam optical imager that saw first light on the 10.4-m Gran Telescopio Canarias (GTC) in 2018. The instrument uses re-imaging optics and 4 dichroic beamsplitters to record $u_s g_s r_s i_s z_s$ ($320-1060$ nm) images simultaneously on its five CCD cameras, each of 3.1 arcmin (diagonal) field of view. The detectors in HiPERCAM are frame-transfer devices cooled thermo-electrically to 183 K, thereby allowing both long-exposure, deep imaging of faint targets, as well as high-speed (over 1000 windowed frames per second) imaging of rapidly varying targets. A comparison-star pick-off system in the telescope focal plane increases the effective field of view to 6.7 arcmin for differential photometry. Combining HiPERCAM with the world's largest optical telescope enables the detection of astronomical sources to $g_s \sim 23$ in 1 s and $g_s \sim 28$ in 1 h. In this paper we describe the scientific motivation behind HiPERCAM, present its design, report on its measured performance, and outline some planned enhancements., Comment: 18 pages, 12 figures, 1 table, accepted for publication in Monthly Notices of the Royal Astronomical Society

AM CVn systems are ultra-compact, helium-rich, accreting binaries with degenerate or semi-degenerate donors. We report the discovery of five new eclipsing AM CVn systems with orbital periods of 61.5, 55.5, 53.3, 37.4, and 35.4 minutes. These systems were discovered by searching for deep eclipses in the Zwicky Transient Facility (ZTF) lightcurves of white dwarfs selected using Gaia parallaxes. We obtained phase-resolved spectroscopy to confirm that all systems are AM CVn binaries, and we obtained high-speed photometry to confirm the eclipse and characterize the systems. The spectra of two long-period systems (61.5 and 53.3 minutes) show many emission and absorption lines, indicating the presence of N, O, Na, Mg, Si, and Ca, and also the K and Zn, elements which have never been detected in AM CVn systems before. By modelling the high-speed photometry, we measured the mass and radius of the donor star, potentially constraining the evolutionary channel that formed these AM CVn systems. We determined that the average mass of the accreting white dwarf is $\approx0.8$$\mathrm{M_{\odot}}$, and that the white dwarfs in long-period systems are hotter than predicted by recently updated theoretical models. The donors have a high entropy and are a factor of $\approx$ 2 more massive compared to zero-entropy donors at the same orbital period. The large donor radius is most consistent with He-star progenitors, although the observed spectral features seem to contradict this. The discovery of 5 new eclipsing AM~CVn systems is consistent with the known observed AM CVn space density and estimated ZTF recovery efficiency. Based on this estimate, we expect to find another 1--4 eclipsing AM CVn systems as ZTF continues to obtain data. This will further increase our understanding of the population, but will require high precision data to better characterize these 5 systems and any new discoveries., Comment: submitted to MNRAS

Using photometric ULTRACAM observations of three new short period cataclysmic variables, we model the primary eclipse lightcurves to extract the orbital separation, masses, and radii of their component stars. We find donor masses of 0.060 +/- 0.008 solar masses, 0.042 +/- 0.001 solar masses, and 0.042 +/- 0.004 solar masses, two being very low-mass sub-stellar donors, and one within 2 sigma of the hydrogen burning limit. All three of the new systems lie close to the modified, "optimal" model evolutionary sequence of Knigge et al. (2011). We briefly re-evaluate the long-standing discrepancy between observed donor mass and radius data, and theoretical CV evolutionary tracks. By looking at the difference in the observed period at each mass and the period predicted by the Knigge et al. (2011) evolutionary sequence, we qualitatively examine the form of excess angular momentum loss that is missing from the models below the period gap. We show indications that the excess angular momentum loss missing from CV models grows in importance relative to gravitational losses as the period decreases. Detailed CV evolutionary models are necessary to draw more quantitative conclusions in the future., Comment: 17 pages, 14 figures, 5 tables