Search

VY Scl binaries are a sub-class of cataclysmic variable (CV) which show extended low states, but do not show outbursts which are seen in other classes of CV. To better determine how often these systems spend in low states and to resolve the state transitions we have analysed ZTF data on eight systems and TESS data on six systems. Half of the sample spent most of the time in a high state; three show a broad range and one spends roughly half the time transitioning between high and low states. Using the ZTF data we explore the colour variation as a function of brightness. In KR Aur, we identify a series of repeating outburst events whose brightness appears to increase over time. Using TESS data we searched for periods other than the orbital. In LN UMa we find evidence for a peak whose period varies between 3-6 d. We outline the current models which aim to explain the observed properties of VY Scl systems which includes disc irradiation and a white dwarf having a significant magnetic field., Comment: 16 pages, 8 Figures, Accepted to MNRAS

SN 2023tsz is a Type Ibn supernova (SNe Ibn) discovered in an extremely low-mass host. SNe Ibn are an uncommon subtype of stripped-envelope core-collapse SNe. They are characterised by narrow helium emission lines in their spectra and are believed to originate from the collapse of massive Wolf-Rayet (WR) stars, though their progenitor systems still remain poorly understood. In terms of energetics and spectrophotometric evolution, SN 2023tsz is largely a typical example of the class, although line profile asymmetries in the nebular phase are seen, which may indicate the presence of dust formation or unshocked circumstellar material. Intriguingly, SN 2023tsz is located in an extraordinarily low-mass host galaxy that is in the 2nd percentile for SESN host masses and star formation rates (SFR). The host has a radius of 1.0 kpc, a $g$-band absolute magnitude of $-12.73$, and an estimated metallicity of $\log(Z_{*}/Z_{\odot}$) = $-1.56$. The SFR and metallicity of the host galaxy raise questions about the progenitor of SN 2023tsz. The low SFR suggests that a star with sufficient mass to evolve into a WR would be uncommon in this galaxy. Further, the very low-metallicity is a challenge for single stellar evolution to enable H and He stripping of the progenitor and produce a SN Ibn explosion. The host galaxy of SN 2023tsz adds another piece to the ongoing puzzle of SNe Ibn progenitors, and demonstrates that they can occur in hosts too faint to be observed in contemporary sky surveys at a more typical SN Ibn redshift., Comment: 14 pages, 10 figures, submitted to MNRAS

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

The main science aim of the BlackGEM array is to detect optical counterparts to gravitational wave mergers. Additionally, the array will perform a set of synoptic surveys to detect Local Universe transients and short time-scale variability in stars and binaries, as well as a six-filter all-sky survey down to ~22nd mag. The BlackGEM Phase-I array consists of three optical wide-field unit telescopes. Each unit uses an f/5.5 modified Dall-Kirkham (Harmer-Wynne) design with a triplet corrector lens, and a 65cm primary mirror, coupled with a 110Mpix CCD detector, that provides an instantaneous field-of-view of 2.7~square degrees, sampled at 0.564\arcsec/pixel. The total field-of-view for the array is 8.2 square degrees. Each telescope is equipped with a six-slot filter wheel containing an optimised Sloan set (BG-u, BG-g, BG-r, BG-i, BG-z) and a wider-band 440-720 nm (BG-q) filter. Each unit telescope is independent from the others. Cloud-based data processing is done in real time, and includes a transient-detection routine as well as a full-source optimal-photometry module. BlackGEM has been installed at the ESO La Silla observatory as of October 2019. After a prolonged COVID-19 hiatus, science operations started on April 1, 2023 and will run for five years. Aside from its core scientific program, BlackGEM will give rise to a multitude of additional science cases in multi-colour time-domain astronomy, to the benefit of a variety of topics in astrophysics, such as infant supernovae, luminous red novae, asteroseismology of post-main-sequence objects, (ultracompact) binary stars, and the relation between gravitational wave counterparts and other classes of transients, Comment: Submitted to PASP

Shock waves are excited by coronal mass ejections (CMEs) and large-scale extreme-ultraviolet (EUV) wave fronts and can result in low-frequency radio emission under certain coronal conditions. In this work, we investigate a moving source of low-frequency radio emission as a CME and an associated EUV wave front move along a channel of a lower density, magnetic field, and Alfv\'en speed in the solar corona. Observations from the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory, the Nan\c{c}ay Radio Heliograph (NRH), and the Irish Low Frequency Array(I-LOFAR) were analysed. Differential emission measure maps were generated to determine densities and Alfv\'en maps, and the kinematics of the EUV wave front was tracked using CorPITA. The radio sources' positions and velocity were calculated from NRH images and I-LOFAR dynamic spectra. The EUV wave expanded radially with a uniform velocity of $\sim$ 500 km s$^{-1}$. However, the radio source was observed to be deflected and appeared to move along a channel of a lower Alfv\'en speed, abruptly slowing from 1700 km s$^{-1}$ to 250 km s$^{-1}$ as it entered a quiet-Sun region. A shock wave with an apparent radial velocity of > 420 km s$^{-1}$ was determined from the drift rate of the associated Type II radio burst. The apparent motion of the radio source may have resulted from a wave front moving along a coronal wave guide or by different points along the wave front emitting at locations with favourable conditions for shock formation., Comment: 6 pages, 3 figures

GRB 201015A is a peculiarly low luminosity, spectrally soft gamma-ray burst (GRB), with $T_{\rm 90} = 9.8 \pm 3.5$ s (time interval of detection of 90\% of photons from the GRB), and an associated supernova (likely to be type Ic or Ic-BL). GRB 201015A has an isotropic energy $E_{\gamma,\rm iso} = 1.75 ^{+0.60} _{-0.53} \times 10^{50}$ erg, and photon index $\Gamma = 3.00 ^{+0.50} _{-0.42}$ (15-150 keV). It follows the Amati relation, a correlation between $E_{\gamma,\rm iso}$ and spectral peak energy $E_{\rm p}$ followed by long GRBs. It appears exceptionally soft based on $\Gamma$, the hardness ratio of HR = $0.47 \pm 0.24$, and low-$E_{\rm p}$, so we have compared it to other GRBs sharing these properties. These events can be explained by shock breakout, poorly collimated jets, and off-axis viewing. Follow-up observations of the afterglow taken in the X-ray, optical, and radio, reveal a surprisingly late flattening in the X-ray from $t = (2.61 \pm 1.27)\times 10^4$ s to $t = 1.67 ^{+1.14} _{-0.65} \times 10^6$ s. We fit the data to closure relations describing the synchrotron emission, finding the electron spectral index to be $p = 2.42 ^{+0.44} _{-0.30}$, and evidence of late-time energy injection with coefficient $q = 0.24 ^{+0.24} _{-0.18}$. The jet half opening angle lower limit ($\theta_{j} \ge 16^{\circ}$) is inferred from the non-detection of a jet break. The launch of SVOM and Einstein Probe in 2023, should enable detection of more low luminosity events like this, providing a fuller picture of the variety of GRBs., Comment: 15 pages, 4 figures

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.

Polars are highly magnetic cataclysmic variables which have been long observed to have both high and low brightness states. The duration of these states has been previously seen to vary from a number of days up to years. Despite this; these states and their physical origin has not been explained in a consistent manner. We present observations of the shortest duration states of a number of Polars observed by ZTF and TESS. This has allowed us to determine that short duration states are a relatively common feature across the population of Polars. Furthermore we have been able to generalise the model of star spot migration to explain both short lived high and low states in Polars by incorporating the interaction between the magnetic field of the white dwarf and that of the star spots., Comment: 16 pages, 9 Figures. Final accepted version to MNRAS

Blue Large Amplitude Pulsators (BLAPs) are a relatively new class of blue variable stars showing periodic variations in their light curves with periods shorter than a few tens of mins and amplitudes of more than ten percent. We report nine blue variable stars identified in the OmegaWhite survey conducted using ESO's VST, which show a periodic modulation in the range 7-37 min and an amplitude in the range 0.11-0.28 mag. We have obtained a series of followup photometric and spectroscopic observations made primarily using SALT and telescopes at SAAO. We find four stars which we identify as BLAPs, one of which was previously known. One star, OW J0820--3301, appears to be a member of the V361 Hya class of pulsating stars and is spatially close to an extended nebula. One further star, OW J1819--2729, has characteristics similar to the sdAV pulsators. In contrast, OW J0815--3421 is a binary star containing an sdB and a white dwarf with an orbital period of 73.7 min, making it only one of six white dwarf-sdB binaries with an orbital period shorter than 80 min. Finally, high cadence photometry of four of the candidate BLAPs show features which we compare with notch-like features seen in the much longer period Cepheid pulsators., Comment: Accepted by MNRAS

Using light curves with $\sim$3 min cadence and a duration of 2 hrs made using the OmegaWhite survey, we present the results of a search for short-period variable stars in the field of 20 open clusters. We identified 92 variable stars in these fields. Using a range of cluster member catalogues and Gaia EDR3 data, we have determined that 10 are cluster members and 2 more are probable members. Based on their position on the Gaia HRD and their photometric periods, we find that most of these are $\delta$ Sct stars. We obtained low-resolution optical spectroscopy of some of these cluster members and field stars. We discuss the cluster variable stars in the context of $\delta$ Sct stars in other open clusters., Comment: Accepted for publication in MNRAS. 20 pages, 14 figures

DW Cnc is an intermediate polar which has previously been observed in both high and low states. Observations of the high state of DW Cnc have previously revealed a spin period at ~ 38.6 min, however observations from the 2018/19 low state showed no evidence of the spin period. We present results from our analysis of 12 s cadence photometric data collected by NGTS of DW Cnc during the high state which began in 2019. Following the previously reported suppression of the spin period signal we identify the return of this signal during the high state, consistent with previous observations of it. We identify this as the restarting of accretion during the high state. We further identified three short outbursts lasting ~ 1 d in DW Cnc with a mean recurrence time of ~ 60 d and an amplitude of ~ 1 mag. These are the first outbursts identified in DW Cnc since 2008. Due to the short nature of these events we identify them not as a result of accretion instabilities but instead either from instabilities originating from the interaction of the magnetorotational instability in the accretion disc and the magnetic field generated by the white dwarf or the result of magnetic gating., Comment: Accepted to MNRAS; 8 pages, 4 figues

PLATO (PLAnetary Transits and Oscillations of stars) is an ESA M-class satellite planned for launch by end 2026 and dedicated to the wide-field search of transiting planets around bright and nearby stars, with a strong focus on discovering habitable rocky planets hosted by solar-like stars. The choice of the fields to be pointed at is a crucial task since it has a direct impact on the scientific return of the mission. In this paper we describe and discuss the formal requirements and the key scientific prioritization criteria that have to be taken into account in the Long-duration Observation Phase (LOP) field selection, and apply a quantitative metric to guide us in this complex optimization process. We identify two provisional LOP fields, one for each hemisphere (LOPS1, LOPN1), and discuss their properties and stellar content. While additional fine-tuning shall be applied to LOP selection before the definitive choice (to be made two years before launch), we expect their position will not move by more than a few degrees with respect to what is proposed in this paper., Comment: 17 pages, 10 figures, 3 tables, A&A in press (accepted: October 26, 2021). Some typos corrected. Author list corrected

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

Context. The ESA PLAnetary Transits and Oscillations of stars (PLATO) mission will search for terrestrial planets in the habitable zone of solar-type stars. Because of telemetry limitations, PLATO targets need to be pre-selected. Aims. In this paper, we present an all sky catalogue that will be fundamental to selecting the best PLATO fields and the most promising target stars, deriving their basic parameters, analysing the instrumental performances, and then planing and optimising follow-up observations. This catalogue also represents a valuable resource for the general definition of stellar samples optimised for the search of transiting planets. Methods. We used Gaia Data Release 2 (DR2) astrometry and photometry and 3D maps of the local interstellar medium to isolate FGK (V$\leq$13) and M (V$\leq$16) dwarfs and subgiant stars. Results. We present the first public release of the all-sky PLATO Input Catalogue (asPIC1.1) containing a total of 2 675 539 stars including 2 378 177 FGK dwarfs and subgiants and 297 362 M dwarfs. The median distance in our sample is 428 pc for FGK stars and 146 pc for M dwarfs, respectively. We derived the reddening of our targets and developed an algorithm to estimate stellar fundamental parameters (Teff, radius, mass) from astrometric and photometric measurements. Conclusions. We show that the overall (internal+external) uncertainties on the stellar parameter determined in the present study are $\sim$230 K (4%) for the effective temperatures, $\sim$0.1 R$_{\odot}$ (9%) for the stellar radii, and $\sim$0.1 M$_{\odot}$ (11%) for the stellar mass. We release a special target list containing all known planet hosts cross-matched with our catalogue., Comment: Accepted for publication in Astronomy & Astrophysics

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.

The advent of wide-field sky surveys has led to the growth of transient and variable source discoveries. The data deluge produced by these surveys has necessitated the use of machine learning (ML) and deep learning (DL) algorithms to sift through the vast incoming data stream. A problem that arises in real-world applications of learning algorithms for classification is imbalanced data, where a class of objects within the data is underrepresented, leading to a bias for over-represented classes in the ML and DL classifiers. We present a recurrent neural network (RNN) classifier that takes in photometric time-series data and additional contextual information (such as distance to nearby galaxies and on-sky position) to produce real-time classification of objects observed by the Gravitational-wave Optical Transient Observer (GOTO), and use an algorithm-level approach for handling imbalance with a focal loss function. The classifier is able to achieve an Area Under the Curve (AUC) score of 0.972 when using all available photometric observations to classify variable stars, supernovae, and active galactic nuclei. The RNN architecture allows us to classify incomplete light curves, and measure how performance improves as more observations are included. We also investigate the role that contextual information plays in producing reliable object classification., Comment: 16 pages, 12 figures, to be published in Monthly Notices of the Royal Astronomical Society

We have adapted the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) Science Pipelines to process data from the Gravitational-Wave Optical Transient Observer (GOTO) prototype. In this paper, we describe how we used the Rubin Observatory LSST Science Pipelines to conduct forced photometry measurements on nightly GOTO data. By comparing the photometry measurements of sources taken on multiple nights, we find that the precision of our photometry is typically better than 20~mmag for sources brighter than 16 mag. We also compare our photometry measurements against colour-corrected PanSTARRS photometry, and find that the two agree to within 10~mmag (1$\sigma$) for bright (i.e., $\sim14^{\rm th}$~mag) sources to 200~mmag for faint (i.e., $\sim18^{\rm th}$~mag) sources. Additionally, we compare our results to those obtained by GOTO's own in-house pipeline, {\sc gotophoto}, and obtain similar results. Based on repeatability measurements, we measure a $5\sigma$ L-band survey depth of between 19 and 20 magnitudes, depending on observing conditions. We assess, using repeated observations of non-varying standard SDSS stars, the accuracy of our uncertainties, which we find are typically overestimated by roughly a factor of two for bright sources (i.e., $<15^{\rm th}$~mag), but slightly underestimated (by roughly a factor of 1.25) for fainter sources ($>17^{\rm th}$~mag). Finally, we present lightcurves for a selection of variable sources, and compare them to those obtained with the Zwicky Transient Factory and GAIA. Despite the Rubin Observatory LSST Science Pipelines still undergoing active development, our results show that they are already delivering robust forced photometry measurements from GOTO data., Comment: Accepted for publication in PASA

Large-scale sky surveys have played a transformative role in our understanding of astrophysical transients, only made possible by increasingly powerful machine learning-based filtering to accurately sift through the vast quantities of incoming data generated. In this paper, we present a new real-bogus classifier based on a Bayesian convolutional neural network that provides nuanced, uncertainty-aware classification of transient candidates in difference imaging, and demonstrate its application to the datastream from the GOTO wide-field optical survey. Not only are candidates assigned a well-calibrated probability of being real, but also an associated confidence that can be used to prioritise human vetting efforts and inform future model optimisation via active learning. To fully realise the potential of this architecture, we present a fully-automated training set generation method which requires no human labelling, incorporating a novel data-driven augmentation method to significantly improve the recovery of faint and nuclear transient sources. We achieve competitive classification accuracy (FPR and FNR both below 1%) compared against classifiers trained with fully human-labelled datasets, whilst being significantly quicker and less labour-intensive to build. This data-driven approach is uniquely scalable to the upcoming challenges and data needs of next-generation transient surveys. We make our data generation and model training codes available to the community., Comment: 17 pages, 12 figures, resubmitted to MNRAS following reviewer comments

We present results of our analysis of up to 15 years of photometric data from eight AM CVn systems with orbital periods between 22.5 and 26.8 min. Our data has been collected from the GOTO, ZTF, Pan-STARRS, ASAS-SN and Catalina all-sky surveys and amateur observations collated by the AAVSO. We find evidence that these interacting ultra-compact binaries show a similar diversity of long term optical properties as the hydrogen accreting dwarf novae. We found that AM CVn systems in the previously identified accretion disc instability region are not a homogenous group. Various members of the analysed sample exhibit behaviour reminiscent of Z Cam systems with long super outbursts and standstills, SU UMa systems with regular, shorter super outbursts, and nova-like systems which appear only in a high state. The addition of TESS full frame images of one of these systems, KL Dra, reveals the first evidence for normal outbursts appearing as a precursor to super outbursts in an AM CVn system. Our results will inform theoretical modelling of the outbursts of hydrogen deficient systems., Comment: 11 Pages, 7 Figures, 2 Tables. Accepted for publication in MNRAS. Author's final submitted version

The past few decades have seen the burgeoning of wide field, high cadence surveys, the most formidable of which will be the Legacy Survey of Space and Time (LSST) to be conducted by the Vera C. Rubin Observatory. So new is the field of systematic time-domain survey astronomy, however, that major scientific insights will continue to be obtained using smaller, more flexible systems than the LSST. One such example is the Gravitational-wave Optical Transient Observer (GOTO), whose primary science objective is the optical follow-up of Gravitational Wave events. The amount and rate of data production by GOTO and other wide-area, high-cadence surveys presents a significant challenge to data processing pipelines which need to operate in near real-time to fully exploit the time-domain. In this study, we adapt the Rubin Observatory LSST Science Pipelines to process GOTO data, thereby exploring the feasibility of using this "off-the-shelf" pipeline to process data from other wide-area, high-cadence surveys. In this paper, we describe how we use the LSST Science Pipelines to process raw GOTO frames to ultimately produce calibrated coadded images and photometric source catalogues. After comparing the measured astrometry and photometry to those of matched sources from PanSTARRS DR1, we find that measured source positions are typically accurate to sub-pixel levels, and that measured L-band photometries are accurate to $\sim50$ mmag at $m_L\sim16$ and $\sim200$ mmag at $m_L\sim18$. These values compare favourably to those obtained using GOTO's primary, in-house pipeline, GOTOPHOTO, in spite of both pipelines having undergone further development and improvement beyond the implementations used in this study. Finally, we release a generic "obs package" that others can build-upon should they wish to use the LSST Science Pipelines to process data from other facilities., Comment: 17 pages, 13 figures. Accepted for publication in PASA

In our previous study of low mass stars using TESS, we found a handful which show a periodic modulation on a period <1 d but also displayed no flaring activity. Here we present the results of a systematic search for Ultra Fast Rotators (UFRs) in the southern ecliptic hemisphere which were observed in 2 min cadence with TESS. Using data from Gaia DR2, we obtain a sample of over 13,000 stars close to the lower main sequence. Of these, we identify 609 stars which lie on the lower main sequence and have a periodic modulation <1 d. The fraction of stars which show flares appears to drop significantly at periods <0.2 d. If the periods are a signature of the rotation rate, this would be a surprise, since faster rotators would be expected to have a stronger magnetic field and, therefore, produce more flares. We explore possible reasons for our finding: the flare inactive stars are members of binaries, in which case the stars rotation rate could have increased as the binary orbital separation reduced due to angular momentum loss over time, or that enhanced emission occurs at blue wavelengths beyond the pass band of TESS. Follow-up spectroscopy and flare monitoring at blue/ultraviolet wavelengths of these flare inactive stars are required to resolve this question., Comment: Accepted for publication in MNRAS

We report the results of optical follow-up observations of 29 gravitational-wave triggers during the first half of the LIGO-Virgo Collaboration (LVC) O3 run with the Gravitational-wave Optical Transient Observer (GOTO) in its prototype 4-telescope configuration (GOTO-4). While no viable electromagnetic counterpart candidate was identified, we estimate our 3D (volumetric) coverage using test light curves of on- and off-axis gamma-ray bursts and kilonovae. In cases where the source region was observable immediately, GOTO-4 was able to respond to a GW alert in less than a minute. The average time of first observation was $8.79$ hours after receiving an alert ($9.90$ hours after trigger). A mean of $732.3$ square degrees were tiled per event, representing on average $45.3$ per cent of the LVC probability map, or $70.3$ per cent of the observable probability. This coverage will further improve as the facility scales up alongside the localisation performance of the evolving gravitational-wave detector network. Even in its 4-telescope prototype configuration, GOTO is capable of detecting AT2017gfo-like kilonovae beyond 200~Mpc in favourable observing conditions. We cannot currently place meaningful electromagnetic limits on the population of distant ($\hat{D}_L = 1.3$~Gpc) binary black hole mergers because our test models are too faint to recover at this distance. However, as GOTO is upgraded towards its full 32-telescope, 2 node (La Palma \& Australia) configuration, it is expected to be sufficiently sensitive to cover the predicted O4 binary neutron star merger volume, and will be able to respond to both northern and southern triggers., Comment: 15 pages, 7 figures, 3 tables. Accepted for publication in MNRAS. Author's final submitted version

On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event., Comment: 52 pages, revised version now accepted for publication in A&A. Abstract abridged to meet arXiv requirements

Detailed studies of the Sun have shown that sunspots and solar flares are closely correlated. Photometric data from Kepler/K2 has allowed similar studies to be carried out on other stars. Here, we utilise TESS photometric 2-min cadence of 167 low mass stars from Sectors 1 - 3 to investigate the relationship between starspots and stellar flares. From our sample, 90 percent show clear rotational modulation likely due to the presence of a large, dominant starspot and we use this to determine a rotational period for each star. Additionally, each low mass star shows one or more flares in its lightcurve and using Gaia DR2 parallaxes and SkyMapper magnitudes we can estimate the energy of the flares in the TESS band-pass. Overall, we have 1834 flares from the 167 low mass stars with energies from $6.0\times 10^{29}$ - $2.4\times 10^{35}$~erg. We find none of the stars in our sample show any preference for rotational phase suggesting the lack of a correlation between the large, dominant star spot and flare number. We discuss this finding in greater detail and present further scenarios to account for the origin of flares on these low mass stars., Comment: 10 pages, 9 figures, 2 tables, accepted for publication in MNRAS

AM CVn binaries are hydrogen deficient compact binaries with an orbital period in the 5-65 min range and are predicted to be strong sources of persistent gravitational wave radiation. Using Gaia Data Release 2, we present the parallaxes and proper motions of 41 out of the 56 known systems. Compared to the parallax determined using the HST Fine Guidance Sensor we find that the archetype star, AM CVn, is significantly closer than previously thought. This resolves the high luminosity and mass accretion rate which models had difficulty in explaining. Using Pan-STARRS1 data we determine the absolute magnitude of the AM CVn stars. There is some evidence that donor stars have a higher mass and radius than expected for white dwarfs or that the donors are not white dwarfs. Using the distances to the known AM CVn stars we find strong evidence that a large population of AM CVn stars have still to be discovered. As this value sets the background to the gravitational wave signal of LISA, this is of wide interest. We determine the mass transfer rate for 15 AM CVn stars and find that the majority have a rate significantly greater than expected from standard models. This is further evidence that the donor star has a greater size than expected., Comment: Accepted by A&A in main journal

Asteroid (3200) Phaethon is a Near-Earth Apollo asteroid with an unusual orbit that brings it closer to the Sun than any other known asteroid. Its last close approach to the Earth was in mid-December 2017 and the next one will be on October 2026. Previous rotationally time-resolved spectroscopy of Phaethon showed that its spectral slope is slightly bluish, in agreement with its B/F taxonomic classification, but at some rotational phases, it changes to slightly reddish. Motivated by this result we performed time-resolved imaging polarimetry of Phaethon during its recent close approach to the Earth. Phaethon has a spin period of 3.604 hours and we found a variation of the linear polarisation with rotation. This seems to be a rare case in which such variation is unambiguously found, also a consequence of its fairly large amplitude. Combining this new information with the brightness and colour variation, as well as previously reported results from Arecibo radar observations, we conclude that there is no variation of the mineralogy across the surface of Phaeton. However, the observed change in the linear polarisation may be related to differences in the thickness of the surface regolith in different areas or local topographic features., Comment: 5 pages, 2 figures, MNRAS Letters

We present an analysis of K2 short cadence data of 34 M dwarfs which have spectral types in the range M0 - L1. Of these stars, 31 showed flares with a duration between $\sim$10-90 min. Using distances obtained from Gaia DR2 parallaxes, we determined the energy of the flares to be in the range $\sim1.2\times10^{29}-6\times10^{34}$ erg. In agreement with previous studies we find rapidly rotating stars tend to show more flares, with evidence for a decline in activity in stars with rotation periods longer than $\sim$10 days. The rotational modulation seen in M dwarf stars is widely considered to result from a starspot which rotates in and out of view. Flux minimum is therefore the rotation phase where we view the main starspot close to the stellar disk center. Surprisingly, having determined the rotational phase of each flare in our study we find none show any preference for rotational phase. We outline three scenarios which could account for this unexpected finding. The relationship between rotation phase and flare rate will be explored further using data from wide surveys such as NGTS and TESS., Comment: Accepted main Journal MNRAS

We present high speed optical photometric observations made using the NTT and ULTRACAM of the optical counterpart of 3XMMJ051034.6--670353, which was recently identified as an X-ray source showing a modulation on a period of 23.6 min. Although the optical counterpart is faint (g=21.4), we find that the u'g'r' light curves show a periodic modulation on a period which is consistent with the X-ray period. We also obtained three low resolution spectra of 3XMMJ051034.6--670353 using the Gemini South Telescope and GMOS. There is no evidence for strong emission lines in the optical spectrum of 3XMMJ051034.6--670353. We compare and contrast the optical and X-ray observations of 3XMMJ051034.6--670353 with the ultra compact binaries HM Cnc and V407 Vul. We find we can identify a distribution of binary masses in which stable direct impact accretion can occur., Comment: Accepted A&A

Ultracompact binaries with orbital periods less than a few hours will dominate the gravitational wave signal in the mHz regime. Until recently, 10 systems were expected have a predicted gravitational wave signal strong enough to be detectable by the Laser Interferometer Space Antenna (LISA), the so-called `verification binaries'. System parameters, including distances, are needed to provide an accurate prediction of the expected gravitational wave strength to be measured by LISA. Using parallaxes from {\sl Gaia} Data Release 2 we calculate signal-to-noise ratios (SNR) for $\approx$50 verification binary candidates. We find that 11 binaries reach a SNR$\geq$20, two further binaries reaching a SNR$\geq$5 and three more systems are expected to have a SNR$\approx$5 after four years integration with LISA. For these 16 systems we present predictions of the gravitational wave amplitude ($\mathcal{A}$) and parameter uncertainties from Fisher information matrix on the amplitude ($\mathcal{A}$) and inclination ($\iota$)., Comment: accepted for publication in MNRAS, 9 pages, 4 figures

We report the discovery of the ultracompact hot subdwarf (sdOB) binary OW J074106.0-294811.0 with an orbital period of P$_{\rm orb}=44.66279\pm1.16\times10^{-4}$ min, making it the most compact hot subdwarf binary known. Spectroscopic observations using the VLT, Gemini and Keck telescopes revealed a He-sdOB primary with an intermediate helium abundance, T$_{\rm eff}=39 400\pm500$ K and log(g)=$5.74\pm0.09$. High signal-to-noise ratio lightcurves show strong ellipsoidal modulation resulting in a derived sdOB mass $M_{\rm sdOB}=0.23\pm0.12$ M$_\odot$ with a WD companion ($M_{\rm WD}=0.72\pm0.17$ M$_\odot$). The mass ratio was found to be $q = M_{\rm sdOB}/M_{\rm WD}=0.32\pm0.10$. The derived mass for the He-sdOB is inconsistent with the canonical mass for hot sbudwarfs of $\approx0.47$ M$_\odot$. To put constraints on the structure and evolutionary history of the sdOB star we compared the derived T$_{\rm eff}$, log(g) and sdOB mass to evolutionary tracks of helium stars and helium white dwarfs calculated with Modules for Experiments in Stellar Astrophysics (MESA). We find that the best fitting model is a helium white dwarf with a mass of $0.320$ M$_\odot$, which left the common envelope ${\approx}1.1$ Myr ago, is consistent with the observations. As a helium white dwarf with a massive white dwarf companion the object will reach contact in 17.6 Myr at an orbital period of 5 min. Depending on the spin-orbit synchronization timescale the object will either merge to form an R CrB star or end up as a stably accreting AM CVn-type system with a helium white dwarf donor., Comment: Accepted for publication in ApJ, 13 pages, 7 figures

We present the discovery and follow-up observations of the second known variable warm DQ white dwarf OW J175358.85-310728.9 (OW J1753-3107). OW J1753-3107 is the brightest of any of the currently known warm or hot DQ and was discovered in the OmegaWhite Survey as exhibiting optical variations on a period of 35.5452 (2) mins, with no evidence for other periods in its light curves. This period has remained constant over the last two years and a single-period sinusoidal model provides a good fit for all follow-up light curves. The spectrum consists of a very blue continuum with strong absorption lines of neutral and ionised carbon, a broad He I 4471 A line, and possibly weaker hydrogen lines. The C I lines are Zeeman split, and indicate the presence of a strong magnetic field. Using spectral Paschen-Back model descriptions, we determine that OW J1753-3107 exhibits the following physical parameters: T_eff = 15430 K, log(g) = 9.0, log(N(C)/N(He)) = -1.2, and the mean magnetic field strength is B_z =2.1 MG. This relatively low temperature and carbon abundance (compared to the expected properties of hot DQs) is similar to that seen in the other warm DQ SDSS J1036+6522. Although OW J1753-3107 appears to be a twin of SDSS J1036+6522, it exhibits a modulation on a period slightly longer than the dominant period in SDSS J1036+6522 and has a higher carbon abundance. The source of variations is uncertain, but they are believed to originate from the rotation of the magnetic white dwarf., Comment: 11 pages, 8 figures, 7 tables. Accepted for publication by MNRAS

We present photometric and spectroscopic follow-up observations of short-period variables discovered in the OmegaWhite survey: a wide-field high-cadence g-band synoptic survey targeting the Galactic Plane. We have used fast photometry on the SAAO 1.0-m and 1.9-m telescopes to obtain light curves of 27 variables, and use these results to validate the period and amplitude estimates from the OmegaWhite processing pipeline. Furthermore, 57 sources (44 unique, 13 also with new light curves) were selected for spectroscopic follow-up using either the SAAO 1.9-m telescope or the Southern African Large Telescope. We find many of these variables have spectra which are consistent with being delta Scuti type pulsating stars. At higher amplitudes, we detect four possible pulsating white dwarf/subdwarf sources and an eclipsing cataclysmic variable. Due to their rarity, these targets are ideal candidates for detailed follow-up studies. From spectroscopy, we confirm the symbiotic binary star nature of two variables identified as such in the SIMBAD database. We also report what could possibly be the first detection of the `Bump Cepheid' phenomena in a delta Scuti star, with OW J175848.21-271653.7 showing a pronounced 22% amplitude dip lasting 3 minutes during each pulsational cycle peak. However, the precise nature of this target is still uncertain as it exhibits the spectral features of a B-type star., Comment: 28 pages, 12 figures, 6 tables. Accepted for publication in MNRAS

We present the results of a series of short radio observations of six ultracool dwarfs made using the upgraded VLA in S (2--4GHz) and C (4--7GHz) bands. LSR J1835+3259 exhibits a 100 percent right-hand circularly polarised burst which shows intense narrowband features with a fast negative frequency drift of about $-30$ MHz $\textrm{s}^{-1}$. They are superimposed on a fainter broadband emission feature with a total duration of about 20 minutes, bandwidth of about 1 GHz, centred at about 3.5 GHz, and a slow positive frequency drift of about 1 MHz $\textrm{s}^{-1}$. This makes it the first such event detected below 4 GHz and the first one exhibiting both positive and negative frequency drifts. Polarised radio emission is also seen in 2MASS J00361617+1821104 and NLTT 33370, while LP 349-25 and TVLM 513-46546 have unpolarised emission and BRI B0021-0214 was not detected. We can reproduce the main characteristics of the burst from LSR J1835+3259 using a model describing the magnetic field of the dwarf as a tilted dipole. We also analyse the origins of the quiescent radio emission and estimate the required parameters of the magnetic field and energetic electrons. Although our results are non-unique, we find a set of models which agree well with the observations., Comment: 16 pages, 9 figures, 5 tables

We report the discovery of 42 white dwarfs in the original Kepler mission field, including nine new confirmed pulsating hydrogen-atmosphere white dwarfs (ZZ Ceti stars). Guided by the Kepler-INT Survey (KIS), we selected white dwarf candidates on the basis of their U-g, g-r, and r-H_alpha photometric colours. We followed up these candidates with high-signal-to-noise optical spectroscopy from the 4.2-m William Herschel Telescope. Using ground-based, time-series photometry, we put our sample of new spectroscopically characterized white dwarfs in the context of the empirical ZZ Ceti instability strip. Prior to our search, only two pulsating white dwarfs had been observed by Kepler. Ultimately, four of our new ZZ Cetis were observed from space. These rich datasets are helping initiate a rapid advancement in the asteroseismic investigation of pulsating white dwarfs, which continues with the extended Kepler mission, K2., Comment: 9 pages, 6 figures, accepted for publication in MNRAS

We present the goals, strategy and first results of the OmegaWhite survey: a wide-field high-cadence $g$-band synoptic survey which aims to unveil the Galactic population of short-period variable stars (with periods $<$ 80 min), including ultracompact binary star systems and stellar pulsators. The ultimate goal of OmegaWhite is to cover 400 square degrees along the Galactic Plane reaching a depth of $g = $ 21.5 mag (10$\sigma$), using OmegaCam on the VLT Survey Telescope (VST). The fields are selected to overlap with surveys such as the Galactic Bulge Survey (GBS) and the VST Photometric H$\alpha$ Survey of the Southern Galactic Plane (VPHAS+) for multi-band colour information. Each field is observed using 38 exposures of 39 s each, with a median cadence of $\sim$2.7 min for a total duration of two hours. Within an initial 26 square degrees, we have extracted the light curves of 1.6 million stars, and have identified 613 variable candidates which satisfy our selection criteria. Furthermore, we present the light curves and statistical properties of 20 sources which have the highest-likelihood of being variable stars. One of these candidates exhibits the colours and light curve properties typically associated with ultracompact AM CVn binaries, although its spectrum exhibits weak Balmer absorption lines and is thus not likely to be such a binary system. We also present follow-up spectroscopy of five other variable candidates, which identifies them as likely low-amplitude $\delta$ Sct pulsating stars., Comment: 27 pages, 10 figures, 9 tables. Accepted for publication by MNRAS, 24 August 2015

IGR J04571+4527 and Swift J0525.6+2416 are two hard X-ray sources detected in the Swift/BAT and INTEGRAL/IBIS surveys. They were proposed to be magnetic cataclysmic variables of the Intermediate Polar (IP) type, based on optical spectroscopy. IGR J04571+4527 also showed a 1218 s optical periodicity, suggestive of the rotational period of a white dwarf, further pointing towards an IP classification. We here present detailed X-ray (0.3-10 keV) timing and spectral analysis performed with XMM-Newton, complemented with hard X-ray coverage (15-70 keV) from Swift/BAT. These are the first high signal to noise observations in the soft X-ray domain for both sources, allowing us to identify the white dwarf X-ray spin period of Swift J0525.6+2416 (226.28 s), and IGR J04571+4527 (1222.6 s). A model consisting of multi-temperature optically thin emission with complex absorption adequately fits the broad-band spectrum of both sources. We estimate a white dwarf mass of about 1.1 and 1.0 solar masses for IGR J04571+4527 and Swift J0525.6+2416, respectively. The above characteristics allow us to unambiguously classify both sources as IPs, confirming the high incidence of this subclass among hard X-ray emitting Cataclysmic Variables., Comment: 8 pages, 4 figures, 3 tables. Accepted for publication in MNRAS

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of accreting white dwarfs. For a summary, we refer to the paper., Comment: White Paper in Support of the Mission Concept of the Large Observatory for X-ray Timing