Search

The determination of the full population of Resident Space Objects (RSOs) in Low Earth Orbit (LEO) is a key issue in the field of space situational awareness that will only increase in importance in the coming years. We endeavour to describe a novel method of inferring the population of RSOs as a function of orbital height and inclination for a range of magnitudes. The method described uses observations of an orbit of known height and inclination to detect RSOs on neighbouring orbits. These neighbouring orbit targets move slowly relative to our tracked orbit, and are thus detectable down to faint magnitudes. We conduct simulations to show that, by observing multiple passes of a known orbit, we can infer the population of RSOs within a defined region of orbital parameter space. Observing a range of orbits from different orbital sites will allow for the inference of a population of LEO RSOs as a function of their orbital parameters and object magnitude., Comment: 15 pages, 12 figures. Accepted for publication in RAS Techniques and Instruments (RASTI)

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

We present the methodology and results of a simulation to determine the recoverability of LEO objects using a blind stacking technique. The method utilises sCMOS and GPU technology to inject and recover LEO objects in real observed data. We explore the target recovery fraction and pipeline run-time as a function of three optimisation parameters; number of frames per data-set, exposure time, and binning factor. Results are presented as a function of magnitude and velocity. We find that target recovery using blind stacking is significantly more complete, and can reach fainter magnitudes, than using individual frames alone. We present results showing that, depending on the combination of optimisation parameters, recovery fraction is up to 90% of detectable targets for magnitudes up to 13.5, and then falls off steadily up to a magnitude limit around 14.5. Run-time is shown to be a few multiples of the observing time for the best combinations of optimisation parameters, approaching real-time processing., Comment: 14 pages, 14 figures. Accepted for publication in Advances in Space Research (ASR)

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.

Large constellations of orbiting communication satellites will become an important source of noise for present and future astronomical observatories. Mitigation measures rely on high quality predictive models of the position and expected brightness of these objects. Optical linear imaging polarimetry holds promise as a quantitative tool to improve our understanding of the physics of reflection of sunlight off satellite components and through which models of expected brightness can be improved. We present the first simultaneous short-timescale linear polarimetry and optical photometry observations of a geostationary satellite, using the new MOPTOP imaging polarimeter on the 2m Liverpool Telescope. Our target, telecommunication satellite Thor-6, shows prominent short timescale glint-like features in the lightcurve, some as short as seconds. Our polarimetric observations overlap with several of these micro-glints, and have the cadence required to resolve them. We find that the polarisation lightcurve is remarkably smooth, the short time scale glints are not seen to produce strong polarimetric features in our observation. We show how short timescale polarimetry can further constrain the properties of the components responsible for these micro-glints., Comment: 6 figures, 30 pages; accepted for publication in Advances in Space Research

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

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

Light curves of the accreting white dwarf pulsator GW Librae spanning a 7.5 month period in 2017 were obtained as part of the Next Generation Transit Survey. This data set comprises 787 hours of photometry from 148 clear nights, allowing the behaviour of the long (hours) and short period (20min) modulation signals to be tracked from night to night over a much longer observing baseline than has been previously achieved. The long period modulations intermittently detected in previous observations of GW Lib are found to be a persistent feature, evolving between states with periods ~83min and 2-4h on time-scales of several days. The 20min signal is found to have a broadly stable amplitude and frequency for the duration of the campaign, but the previously noted phase instability is confirmed. Ultraviolet observations obtained with the Cosmic Origin Spectrograph onboard the Hubble Space Telescope constrain the ultraviolet-to-optical flux ratio to ~5 for the 4h modulation, and <=1 for the 20min period, with caveats introduced by non-simultaneous observations. These results add further observational evidence that these enigmatic signals must originate from the white dwarf, highlighting our continued gap in theoretical understanding of the mechanisms that drive them., Comment: Accepted for publication in MNRAS, 8 pages, 7 figures

The Gravitational-wave Optical Transient Observer (GOTO) is a wide-field telescope project focused on detecting optical counterparts to gravitational-wave sources. The GOTO Telescope Control System (G-TeCS) is a custom robotic control system which autonomously manages the GOTO telescope hardware and nightly operations. Since the commissioning the GOTO prototype on La Palma in 2017, development of the control system has focused on the alert handling and scheduling systems. These allow GOTO to receive and process transient alerts and then schedule and carry out observations, all without the need for human involvement. GOTO is ultimately anticipated to include multiple telescope arrays on independent mounts, both on La Palma and at a southern site in Australia. When complete these mounts will be linked to form a single multi-site observatory, requiring more advanced scheduling systems to best optimise survey and follow-up observations., Comment: 12 pages, 6 figures, 2 tables, submitted to SPIE Astronomical Telescopes + Instrumentation 2020

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

The pulsating hydrogen atmosphere white dwarf star G 117-B15A has been observed since 1974. Its main pulsation period at 215.19738823(63) s, observed in optical light curves, varies by only (5.12+/-0.82)x10^{-15} s/s and shows no glitches, as pulsars do. The observed rate of period change corresponds to a change of the pulsation period by 1 s in 6.2 million years. We demonstrate that this exceptional optical clock can continue to put stringent limits on fundamental physics, such as constraints on interaction from hypothetical dark matter particles, as well as to search for the presence of external substellar companions., Comment: 14 pages, 2 figures, 2 tables. arXiv admin note: substantial text overlap with arXiv:astro-ph/0507487; replaced because Fig. 2 was overlaying text

We present time-resolved optical and ultraviolet spectroscopy and photometry of V1460~Her, an eclipsing cataclysmic variable with a 4.99\,h orbital period and an overluminous K5-type donor star. The optical spectra show emission lines from an accretion disc along with absorption lines from the donor. We use these to measure radial velocities, which, together with constraints upon the orbital inclination from photometry, imply masses of $M_1=0.869\pm0.006\,\mathrm{M}_\odot$ and $M_2=0.295\pm0.004\,\mathrm{M}_\odot$ for the white dwarf and the donor. The radius of the donor, $R_2=0.43\pm0.002\,\mathrm{R}_\odot$, is $\approx 50$ per cent larger than expected given its mass, while its spectral type is much earlier than the M3.5 type that would be expected from a main sequence star with a similar mass. HST spectra show strong $\mathrm{N{\small V}}$ 1240 A emission but no $\mathrm{C{\small IV}}$ 1550 A emission, evidence for CNO-processed material. The donor is therefore a bloated, over-luminous remnant of a thermal-timescale stage of high mass transfer and has yet to re-establish thermal equilibrium. Remarkably, the HST ultraviolet data also show a strong 30 per cent peak-to-peak, $38.9\,$s pulsation that we explain as being due to the spin of the white dwarf, potentially putting V1460 Her in a similar category to the propeller system AE Aqr in terms of its spin frequency and evolutionary path. AE Aqr also features a post-thermal timescale mass donor, and V1460 Her may therefore be its weak magnetic field analogue since the accretion disc is still present, with the white dwarf spin-up a result of a recent high accretion rate.

Recent anomalies exhibited by satellites and rocket bodies have highlighted that a population of faint debris exists at geosynchronous (GEO) altitudes, where there are no natural removal mechanisms. Despite previous optical surveys probing to around 10-20 cm in size, regular monitoring of faint sources at GEO is challenging, thus our knowledge remains sparse. It is essential that we continue to explore the faint debris population using large telescopes to better understand the risk posed to active GEO satellites. To this end, we present photometric results from a survey of the GEO region carried out with the 2.54 m Isaac Newton Telescope in La Palma, Canary Islands. We probe to 21st visual magnitude (around 10 cm, assuming Lambertian spheres with an albedo of 0.1), uncovering 129 orbital tracks with GEO-like motion across the eight nights of dark-grey time comprising the survey. The faint end of our brightness distribution continues to rise until the sensitivity limit of the sensor is reached, suggesting that the modal brightness could be even fainter. We uncover a number of faint, uncatalogued objects that show photometric signatures of rapid tumbling, many of which straddle the limiting magnitude of our survey over the course of a single exposure, posing a complex issue when estimating object size. This work presents the first instalment of DebrisWatch, an ongoing collaboration between the University of Warwick and the Defence Science and Technology Laboratory (UK) investigating the faint population of GEO debris., Comment: 24 pages, 8 figures. Accepted for publication in Adv. Space Res

The amount of observational data produced by time-domain astronomy is exponentially in-creasing. Human inspection alone is not an effective way to identify genuine transients fromthe data. An automatic real-bogus classifier is needed and machine learning techniques are commonly used to achieve this goal. Building a training set with a sufficiently large number of verified transients is challenging, due to the requirement of human verification. We presentan approach for creating a training set by using all detections in the science images to be thesample of real detections and all detections in the difference images, which are generated by the process of difference imaging to detect transients, to be the samples of bogus detections. This strategy effectively minimizes the labour involved in the data labelling for supervised machine learning methods. We demonstrate the utility of the training set by using it to train several classifiers utilizing as the feature representation the normalized pixel values in 21-by-21pixel stamps centered at the detection position, observed with the Gravitational-wave Optical Transient Observer (GOTO) prototype. The real-bogus classifier trained with this strategy can provide up to 95% prediction accuracy on the real detections at a false alarm rate of 1%., Comment: 9 pages, 8 figures

We report the discovery of two transiting exoplanets from the WASP survey, WASP-150b and WASP-176b. WASP-150b is an eccentric ($e$ = 0.38) hot Jupiter on a 5.6 day orbit around a $V$ = 12.03, F8 main-sequence host. The host star has a mass and radius of 1.4 $\rm M_{\odot}$ and 1.7 $\rm R_{\odot}$ respectively. WASP-150b has a mass and radius of 8.5 $\rm M_J$ and 1.1 $\rm R_J$, leading to a large planetary bulk density of 6.4 $\rm \rho_J$. WASP-150b is found to be $\sim3$ Gyr old, well below its circularisation timescale, supporting the eccentric nature of the planet. WASP-176b is a hot Jupiter planet on a 3.9 day orbit around a $V$ = 12.01, F9 sub-giant host. The host star has a mass and radius of 1.3 $\rm M_{\odot}$ and 1.9 $\rm R_{\odot}$. WASP-176b has a mass and radius of 0.86 $\rm M_J$ and 1.5 $\rm R_J$ respectively, leading to a planetary bulk density of 0.23 $\rm \rho_J$., Comment: 15 pages, 14 figures. Accepted for publication in The Astronomical Journal

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

White dwarfs are dense, cooling stellar embers consisting mostly of carbon and oxygen, or oxygen and neon (with a few percent carbon) at higher initial stellar masses. These stellar cores are enveloped by a shell of helium which in turn is usually surrounded by a layer of hydrogen, generally prohibiting direct observation of the interior composition. However, carbon is observed at the surface of a sizeable fraction of white dwarfs, sometimes with traces of oxygen, and it is thought to be dredged-up from the core by a deep helium convection zone. In these objects only traces of hydrogen are found as large masses of hydrogen are predicted to inhibit hydrogen/helium convective mixing within the envelope. We report the identification of WDJ055134.612+413531.09, an ultra-massive (1.14 $M_\odot$) white dwarf with a unique hydrogen/carbon mixed atmosphere (C/H=0.15 in number ratio). Our analysis of the envelope and interior indicates that the total hydrogen and helium mass fractions must be several orders of magnitude lower than predictions of single star evolution: less than $10^{-9.5}$ and $10^{-7.0}$, respectively. Due to the fast kinematics ($129\pm5$ km/s relative to the local standard of rest), large mass, and peculiar envelope composition, we argue that WDJ0551+4135 is consistent with formation from the merger of two white dwarfs in a tight binary system., Comment: Published in Nature Astronomy Letters on March 2nd 2020, https://www.nature.com/articles/s41550-020-1028-0

We report the discovery of a new ultra-short period transiting hot Jupiter from the Next Generation Transit Survey (NGTS). NGTS-10b has a mass and radius of $2.162\,^{+0.092}_{-0.107}$ M$_{\rm J}$ and $1.205\,^{+0.117}_{-0.083}$ R$_{\rm J}$ and orbits its host star with a period of $0.7668944\pm0.0000003$ days, making it the shortest period hot Jupiter yet discovered. The host is a $10.4\pm2.5$ Gyr old K5V star ($T_\mathrm{eff}$=$4400\pm100$\,K) of Solar metallicity ([Fe/H] = $-0.02\pm0.12$\,dex) showing moderate signs of stellar activity. NGTS-10b joins a short list of ultra-short period Jupiters that are prime candidates for the study of star-planet tidal interactions. NGTS-10b orbits its host at just $1.46\pm0.18$ Roche radii, and we calculate a median remaining inspiral time of $38$\,Myr and a potentially measurable transit time shift of $7$\,seconds over the coming decade, assuming a stellar tidal quality factor $Q'_{\rm s}=2\times10^{7}$., Comment: 16 pages, 19 figures and 5 tables. Submitted 27 Sept 2019. Accepted 10 Jan 2020. Published 20 Feb 2020

PDS 110 is a young disk-hosting star in the Orion OB1A association. Two dimming events of similar depth and duration were seen in 2008 (WASP) and 2011 (KELT), consistent with an object in a closed periodic orbit. In this paper we present data from a ground-based observing campaign designed to measure the star both photometrically and spectroscopically during the time of predicted eclipse in September 2017. Despite high-quality photometry, the predicted eclipse did not occur, although coherent structure is present suggesting variable amounts of stellar flux or dust obscuration. We also searched for RV oscillations caused by any hypothetical companion and can rule out close binaries to 0.1 $M_\odot$. A search of Sonneberg plate archive data also enabled us to extend the photometric baseline of this star back more than 50 years, and similarly does not re-detect any deep eclipses. Taken together, they suggest that the eclipses seen in WASP and KELT photometry were due to aperiodic events. It would seem that PDS 110 undergoes stochastic dimmings that are shallower and shorter-duration than those of UX Ori variables, but may have a similar mechanism., Comment: Accepted to MNRAS; 12 pages, 7 figures; Supplementary photometric data in zipped latex source as all_photometry.csv

We present new Hubble Space Telescope (HST) ultraviolet and ground-based optical observations of the hot, metal-rich white dwarf GD 394. Extreme-ultraviolet (EUV) observations in 1992-1996 revealed a 1.15d periodicity with a 25 percent amplitude, hypothesised to be due to metals in a surface accretion spot. We obtained phase-resolved HST/Space Telescope Imaging Spectrograph (STIS) high-resolution far-ultraviolet (FUV) spectra of GD 394 that sample the entire period, along with a large body of supplementary data. We find no evidence for an accretion spot, with the flux, accretion rate and radial velocity of GD 394 constant over the observed timescales at ultraviolet and optical wavelengths. We speculate that the spot may have no longer been present when our observations were obtained, or that the EUV variability is being caused by an otherwise undetected evaporating planet. The atmospheric parameters obtained from separate fits to optical and ultraviolet spectra are inconsistent, as is found for multiple hot white dwarfs. We also detect non-photospheric, high-excitation absorption lines of multiple volatile elements, which could be evidence for a hot plasma cocoon surrounding the white dwarf., Comment: Accepted for publication in MNRAS

We report the discovery of NGTS-4b, a sub-Neptune-sized planet transiting a 13th magnitude K-dwarf in a 1.34d orbit. NGTS-4b has a mass M=$20.6\pm3.0$M_E and radius R=$3.18\pm0.26$R_E, which places it well within the so-called "Neptunian Desert". The mean density of the planet ($3.45\pm0.95$g/cm^3) is consistent with a composition of 100% H$_2$O or a rocky core with a volatile envelope. NGTS-4b is likely to suffer significant mass loss due to relatively strong EUV/X-ray irradiation. Its survival in the Neptunian desert may be due to an unusually high core mass, or it may have avoided the most intense X-ray irradiation by migrating after the initial activity of its host star had subsided. With a transit depth of $0.13\pm0.02$%, NGTS-4b represents the shallowest transiting system ever discovered from the ground, and is the smallest planet discovered in a wide-field ground-based photometric survey., Comment: Submitted to MNRAS

WD 1145+017 is currently the only white dwarf known to exhibit periodic transits of planetary debris as well as absorption lines from circumstellar gas. We present the first simultaneous fast optical spectrophotometry and broad-band photometry of the system, obtained with the Gran Telescopio Canarias (GTC) and the Liverpool Telescope (LT), respectively. The observations spanned $5.5$ h, somewhat longer than the $4.5$-h orbital period of the debris. Dividing the GTC spectrophotometry into five wavelength bands reveals no significant colour differences, confirming grey transits in the optical. We argue that absorption by an optically thick structure is a plausible alternative explanation for the achromatic nature of the transits that can allow the presence of small-sized ($\sim\mu$m) particles. The longest ($87$ min) and deepest ($50$ per cent attenuation) transit recorded in our data exhibits a complex structure around minimum light that can be well modelled by multiple overlapping dust clouds. The strongest circumstellar absorption line, Fe II $\lambda$5169, significantly weakens during this transit, with its equivalent width reducing from a mean out-of-transit value of $2$ \AA\ to $1$ \AA\ in-transit, supporting spatial correlation between the circumstellar gas and dust. Finally, we made use of the Gaia Data Release 2 and archival photometry to determine the white dwarf parameters. Adopting a helium-dominated atmosphere containing traces of hydrogen and metals, and a reddening $E(B-V)=0.01$ we find $T_\mathrm{eff}=15\,020 \pm 520$ K, $\log g=8.07\pm0.07$, corresponding to $M_\mathrm{WD}=0.63\pm0.05\ \mbox{$\mathrm{M}_{\odot}$}$ and a cooling age of $224\pm30$ Myr., Comment: 13 pages, 9 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society (2018 Aug 22)

The Gravitational-wave Optical Transient Observer (GOTO) is a wide-field telescope project aimed at detecting optical counterparts to gravitational wave sources. The prototype instrument was inaugurated in July 2017 on La Palma in the Canary Islands. We describe the GOTO Telescope Control System (G-TeCS), a custom robotic control system written in Python which autonomously manages the telescope hardware and nightly operations. The system comprises of multiple independent control daemons, which are supervised by a master control program known as the "pilot". Observations are decided by a "just-in-time" scheduler, which instructs the pilot what to observe in real time and provides quick follow-up of transient events., Comment: 14 pages, 6 figures, 1 table, SPIE Astronomical Telescopes + Instrumentation, Austin, Texas, 10-15 June 2018

AR Scorpii is unique amongst known white dwarf binaries in showing powerful pulsations extending to radio frequencies. Here we aim to investigate the multi-frequency radio emission of AR Sco in detail, in order to constrain its origin and emission mechanisms. We present interferometric radio frequency imaging of AR Sco at 1.5, 5 and 9 GHz, analysing the total flux and polarization behaviour of this source at high time resolution (10, 3 and 3 s), across a full 3.6 hr orbital period in each band. We find strong modulation of the radio flux on the orbital period and the orbital sideband of the white dwarf's spin period (also known as the "beat" period). This indicates that, like the optical flux, the radio flux arises predominantly from on or near the inner surface of the M-dwarf companion star. The beat-phase pulsations of AR Sco decrease in strength with decreasing frequency. They are strongest at 9 GHz and at an orbital phase ~0.5. Unlike the optical emission from this source, radio emission from AR Sco shows weak linear polarization but very strong circular polarization, reaching ~30% at an orbital phase ~0.8. We infer the probable existence of a non-relativistic cyclotron emission component, which dominates at low radio frequencies. Given the required magnetic fields, this also likely arises from on or near the M-dwarf., Comment: 14 pages. A&A accepted

We describe the Next Generation Transit Survey (NGTS), which is a ground-based project searching for transiting exoplanets orbiting bright stars. NGTS builds on the legacy of previous surveys, most notably WASP, and is designed to achieve higher photometric precision and hence find smaller planets than have previously been detected from the ground. It also operates in red light, maximising sensitivity to late K and early M dwarf stars. The survey specifications call for photometric precision of 0.1 per cent in red light over an instantaneous field of view of 100 square degrees, enabling the detection of Neptune-sized exoplanets around Sun-like stars and super-Earths around M dwarfs. The survey is carried out with a purpose-built facility at Cerro Paranal, Chile, which is the premier site of the European Southern Observatory (ESO). An array of twelve 20cm f/2.8 telescopes fitted with back-illuminated deep-depletion CCD cameras are used to survey fields intensively at intermediate Galactic latitudes. The instrument is also ideally suited to ground-based photometric follow-up of exoplanet candidates from space telescopes such as TESS, Gaia and PLATO. We present observations that combine precise autoguiding and the superb observing conditions at Paranal to provide routine photometric precision of 0.1 per cent in 1 hour for stars with I-band magnitudes brighter than 13. We describe the instrument and data analysis methods as well as the status of the survey, which achieved first light in 2015 and began full survey operations in 2016. NGTS data will be made publicly available through the ESO archive., Comment: 20 pages, 17 Figures, Accepted for publication in MNRAS

We present photometry and spectroscopy for 27 pulsating hydrogen-atmosphere white dwarfs (DAVs, a.k.a. ZZ Ceti stars) observed by the Kepler space telescope up to K2 Campaign 8, an extensive compilation of observations with unprecedented duration (>75 days) and duty cycle (>90%). The space-based photometry reveals pulsation properties previously inaccessible to ground-based observations. We observe a sharp dichotomy in oscillation mode linewidths at roughly 800 s, such that white dwarf pulsations with periods exceeding 800 s have substantially broader mode linewidths, more reminiscent of a damped harmonic oscillator than a heat-driven pulsator. Extended Kepler coverage also permits extensive mode identification: We identify the spherical degree of 61 out of 154 unique radial orders, providing direct constraints of the rotation period for 20 of these 27 DAVs, more than doubling the number of white dwarfs with rotation periods determined via asteroseismology. We also obtain spectroscopy from 4m-class telescopes for all DAVs with Kepler photometry. Using these homogeneously analyzed spectra we estimate the overall mass of all 27 DAVs, which allows us to measure white dwarf rotation as a function of mass, constraining the endpoints of angular momentum in low- and intermediate-mass stars. We find that 0.51-to-0.73-solar-mass white dwarfs, which evolved from 1.7-to-3.0-solar-mass ZAMS progenitors, have a mean rotation period of 35 hr with a standard deviation of 28 hr, with notable exceptions for higher-mass white dwarfs. Finally, we announce an online repository for our Kepler data and follow-up spectroscopy, which we collect at http://www.k2wd.org., Comment: 33 pages, 31 figures, 5 tables; accepted for publication in ApJS. All raw and reduced data are collected at http://www.k2wd.org

The upcoming NASA Transiting Exoplanet Survey Satellite (TESS) will obtain space-based uninterrupted light curves for a large sample of bright white dwarfs distributed across the entire sky, providing a very rich resource for asteroseismological studies and the search for transits from planetary debris. We have compiled an all-sky catalogue of ultraviolet, optical, and infrared photometry as well as proper motions, which we propose as an essential tool for the preliminary identification and characterisation of potential targets. We present data for 1864 known white dwarfs and 305 high-probability white dwarf candidates brighter than 17 mag. We describe the spectroscopic follow-up of 135 stars, of which 82 are white dwarfs and 25 are hot subdwarfs. The new confirmed stars include six pulsating white dwarf candidates (ZZ Cetis), and nine white dwarf binaries with a cool main-sequence companion. We identify one star with a spectroscopic distance of only 25 pc from the Sun. Around the time TESS is launched, we foresee that all white dwarfs in this sample will have trigonometric parallaxes measured by the ESA Gaia mission next year., Comment: Accepted for publication on MNRAS. 21 pages, 8 tables. Machine readable tables, which will be made also available through VizieR, are also included

The unprecedented extent of coverage provided by Kepler observations recently revealed outbursts in two hydrogen-atmosphere pulsating white dwarfs (DAVs) that cause hours-long increases in the overall mean flux of up to 14%. We have identified two new outbursting pulsating white dwarfs in K2, bringing the total number of known outbursting white dwarfs to four. EPIC 211629697, with T_eff = 10,780 +/- 140 K and log(g) = 7.94 +/- 0.08, shows outbursts recurring on average every 5.0 d, increasing the overall flux by up to 15%. EPIC 229227292, with T_eff = 11,190 +/- 170 K and log(g) = 8.02 +/- 0.05, has outbursts that recur roughly every 2.4 d with amplitudes up to 9%. We establish that only the coolest pulsating white dwarfs within a small temperature range near the cool, red edge of the DAV instability strip exhibit these outbursts., Comment: 11 pages, 9 figures, accepted for publication in ApJ

We report continued Hubble Space Telescope (HST) ultraviolet spectra and ground-based optical photometry and spectroscopy of GW Librae eight years after its largest known dwarf nova outburst in 2007. This represents the longest cooling timescale measured for any dwarf nova. The spectra reveal that the white dwarf still remains about 3000 K hotter than its quiescent value. Both ultraviolet and optical light curves show a short period of 364-373 s, similar to one of the non-radial pulsation periods present for years prior to the outburst, and with a similar large UV/optical amplitude ratio. A large modulation at a period of 2 h (also similar to that observed prior to outburst) is present in the optical data preceding and during the HST observations, but the satellite observation intervals did not cover the peaks of the optical modulation so it is not possible to determine its corresponding UV amplitude. The similarity of the short and long periods to quiescent values implies the pulsating, fast spinning white dwarf in GW Lib may finally be nearing its quiescent configuration., Comment: 6 figures, accepted in AJ

We employ an asteroseismic model of L19-2, a relatively massive ($M_{\star} \sim 0.75 M_{\sun}$) and hot ($T_{\rm eff} \sim 12\,100$ K) pulsating DA (H-rich atmosphere) white dwarf star (DAV or ZZ Ceti variable), and use the observed values of the temporal rates of period change of its dominant pulsation modes ($\Pi \sim 113$ s and $\Pi \sim 192$ s), to derive a new constraint on the mass of the axion, the hypothetical non-barionic particle considered as a possible component of the dark matter of the Universe. If the asteroseismic model employed is an accurate representation of L19-2, then our results indicate hints of extra cooling in this star, compatible with emission of axions of mass $m_{\rm a} \cos^2 \beta \lesssim 25$ meV or an axion-electron coupling constant of $g_{\rm ae} \lesssim 7 \times 10^{-13}$., Comment: 20 pages, 6 figures, 3 tables. Revised version. Accepted for publication in the Journal of Cosmology and Astroparticle Physics. arXiv admin note: text overlap with arXiv:1205.6180

We present new time-series photometry of the accreting pulsating white dwarf system GW Librae obtained in 2012 and 2013 at the University of Canterbury Mt John Observatory in New Zealand. Our 2012 data show the return of a $\sim$19 minute periodicity that was previously detected in 2008. This pulsation mode was a dominant feature of our quality May 2012 data set, which consisted of six contiguous nights; a detailed analysis indicated a degree of frequency variability. We show by comparison with the previously identified pulsation modes that this periodicity is best explained as a new mode, and that the quasi-stability of the periods appears to be a general feature of the pulsations in these systems. We also find a previously unreported 3-hour modulation period, which we believe to be related to the known two and four hour periods of so far unknown origin., Comment: 10 pages, 8 figures, 3 tables

We obtained high-speed photometry of the disintegrating planetesimals orbiting the white dwarf WD1145+017, spanning a period of four weeks. The light curves show a dramatic evolution of the system since the first observations obtained about seven months ago. Multiple transit events are detected in every light curve, which have varying durations(~3-12min) and depths (~10-60%). The time-averaged extinction is ~11%, much higher than at the time of the Kepler observations. The shortest-duration transits require that the occulting cloud of debris has a few times the size of the white dwarf, longer events are often resolved into the superposition of several individual transits. The transits evolve on time scales of days, both in shape and in depth, with most of them gradually appearing and disappearing over the course of the observing campaign. Several transits can be tracked across multiple nights, all of them recur on periods of ~4.49h, indicating multiple planetary debris fragments on nearly identical orbits. Identifying the specific origin of these bodies within this planetary system, and the evolution leading to their current orbits remains a challenging problem., Comment: Astrophysical Journal Letters in press

We present observations of a new phenomenon in pulsating white dwarf stars: large-amplitude outbursts at timescales much longer than the pulsation periods. The cool (Teff = 11,010 K), hydrogen-atmosphere pulsating white dwarf PG 1149+057 was observed nearly continuously for more than 78.8 d by the extended Kepler mission in K2 Campaign 1. The target showed 10 outburst events, recurring roughly every 8 d and lasting roughly 15 hr, with maximum flux excursions up to 45% in the Kepler bandpass. We demonstrate that the outbursts affect the pulsations and therefore must come from the white dwarf. Additionally, we argue that these events are not magnetic reconnection flares, and are most likely connected to the stellar pulsations and the relatively deep surface convection zone. PG 1149+057 is now the second cool pulsating white dwarf to show this outburst phenomenon, after the first variable white dwarf observed in the Kepler mission, KIC 4552982. Both stars have the same effective temperature, within the uncertainties, and are among the coolest known pulsating white dwarfs of typical mass. These outbursts provide fresh observational insight into the red edge of the DAV instability strip and the eventual cessation of pulsations in cool white dwarfs., Comment: 6 pages, 5 figure, accepted for publication in the Astrophysical Journal Letters

We present an analysis of the binary and physical parameters of a unique pulsating white dwarf with a main-sequence companion, SDSS J1136+0409, observed for more than 77 d during the first pointing of the extended Kepler mission: K2 Campaign 1. Using new ground-based spectroscopy, we show that this post-common-envelope binary has an orbital period of 6.89760103(60) hr, which is also seen in the photometry as a result of Doppler beaming and ellipsoidal variations of the secondary. We spectroscopically refine the temperature of the white dwarf to 12330(260) K and its mass to 0.601(36) Msun. We detect seven independent pulsation modes in the K2 light curve. A preliminary asteroseismic solution is in reasonable agreement with the spectroscopic atmospheric parameters. Three of the pulsation modes are clearly rotationally split multiplets, which we use to demonstrate that the white dwarf is not synchronously rotating with the orbital period but has a rotation period of 2.49(53) hr. This is faster than any known isolated white dwarf, but slower than almost all white dwarfs measured in non-magnetic cataclysmic variables, the likely future state of this binary., Comment: 13 pages, 7 figures, accepted for publication in MNRAS

We have designed a realistic simulation of astronomical observing using a relatively low-cost commercial CCD camera and a microcontroller-based circuit that drives LEDs inside a light-tight box with time-varying intensities. As part of a laboratory experiment, students can acquire sequences of images using the camera, and then perform data analysis using a language such as MATLAB or Python to: (a) extract the intensity of the imaged LEDs, (b) perform basic calibrations on the time-series data, and (c) convert their data into the frequency domain where they can then identify the frequency structure. The primary focus is on studying light curves produced by the pulsating white dwarf stars. The exercise provides an introduction to CCD observing, a framework for teaching concepts in numerical data analysis and Fourier techniques, and connections with the physics of white dwarf stars., Comment: 10 pages, 9 figures, submitted to American Journal of Physics

We describe a portable CCD-based instrumentation system designed to efficiently undertake high time precision fast photometry. The key components of the system are (1) an externally triggered commercial frame-transfer CCD, (2) a custom GPS-derived time source, and (3) flexible software for both instrument control and online analysis/display. Two working instruments that implement this design are described. The New Zealand based instrument employs a Princeton Instruments (PI) 1k x 1k CCD and has been used with the 1 m telescope at Mt John University Observatory, while the other uses a newer 1k x 1k electron-multiplying CCD supplied by PI and is based at the University of Texas at Austin. We include some recent observations that illustrate the capabilities of the instruments., Comment: 8 pages, 9 figures, 2 tables

We present the analysis of 71 hours of high quality time-series CCD photometry of the helium atmosphere pulsating white dwarf (DBV) EC 04207-4748 obtained using the facilities at Mt John University Observatory in New Zealand. The photometric data set consists of four week-long observing sessions covering the period March to November 2011. A Fourier analysis of the lightcurves yielded clear evidence of four independent eigenmodes in the star with the dominant mode having a period of 447 s. The lightcurve variations exhibit distinct nonsinusoidal shapes, which results in significant harmonics of the dominant frequency appearing in the Fourier transforms. These observed variations are interpreted in terms of nonlinear contributions from the energy flux transmission through the subsurface convection zone in the star. Our modelling of this mechanism, using the methods first introduced by Montgomery (2005), yields a time-averaged convective response time of tau_0 ~ 150 s for the star, and this is shown to be broadly consistent with a MLT/alpha parameter value between 0.8 and 1.2. It is argued that for the DBV pulsators the measured value of tau_0 is a better estimate of the relative stellar surface temperatures than those obtained via spectroscopic techniques., Comment: 9 pages, 3 figures, 6 tables

We carry out high-speed photometry on 20 of the shortest-period, detached white dwarf binaries known and discover systems with eclipses, ellipsoidal variations (due to tidal deformations of the visible white dwarf), and Doppler beaming. All of the binaries contain low-mass white dwarfs with orbital periods less than 4 hr. Our observations identify the first eight tidally distorted white dwarfs, four of which are reported for the first time here, which we use to put empirical constraints on the mass-radius relationship for extremely low-mass (<0.30 Msun) white dwarfs. We also detect Doppler beaming in several of these binaries, which confirms the high-amplitude radial-velocity variability. All of these systems are strong sources of gravitational radiation, and long-term monitoring of those that display ellipsoidal variations can be used to detect spin-up of the tidal bulge due to orbital decay., Comment: 14 pages, 5 figures, accepted for publication in The Astrophysical Journal

Characterizing a microlensing planet is done from modeling an observed lensing light curve. In this process, it is often confronted that solutions of different lensing parameters result in similar light curves, causing difficulties in uniquely interpreting the lens system, and thus understanding the causes of different types of degeneracy is important. In this work, we show that incomplete coverage of a planetary perturbation can result in degenerate solutions even for events where the planetary signal is detected with a high level of statistical significance. We demonstrate the degeneracy for an actually observed event OGLE-2012-BLG-0455/MOA-2012-BLG-206. The peak of this high-magnification event $(A_{\rm max}\sim400)$ exhibits very strong deviation from a point-lens model with $\Delta\chi^{2}\gtrsim4000$ for data sets with a total number of measurement 6963. From detailed modeling of the light curve, we find that the deviation can be explained by four distinct solutions, i.e., two very different sets of solutions, each with a two-fold degeneracy. While the two-fold (so-called "close/wide") degeneracy is well-understood, the degeneracy between the radically different solutions is not previously known. The model light curves of this degeneracy differ substantially in the parts that were not covered by observation, indicating that the degeneracy is caused by the incomplete coverage of the perturbation. It is expected that the frequency of the degeneracy introduced in this work will be greatly reduced with the improvement of the current lensing survey and follow-up experiments and the advent of new surveys., Comment: 5 pages, 3 figures, ApJ accepted

We present optical and X-ray observations of two tidally distorted, extremely low-mass white dwarfs (WDs) with massive companions. There is no evidence of neutron stars in our Chandra and XMM observations of these objects. SDSS J075141.18$-$014120.9 (J0751) is an eclipsing double WD binary containing a 0.19 Msol WD with a 0.97 Msol companion in a 1.9 h orbit. J0751 becomes the fifth eclipsing double WD system currently known. SDSS J174140.49+652638.7 (J1741) is another binary containing a 0.17 Msol WD with an unseen M > 1.11 Msol WD companion in a 1.5 h orbit. With a mass ratio of ~0.1, J1741 will have stable mass transfer through an accretion disk and turn into an interacting AM Canum Venaticorum (AM CVn) system in the next ~160 Myr. With a mass ratio of 0.2, J0751 is likely to follow a similar evolutionary path. These are the first known AM CVn progenitor binary systems and they provide important constraints on the initial conditions for AM CVn. Theoretical studies suggest that both J0751 and J1741 may create thermonuclear supernovae in ~10^8 yr, either .Ia or Ia. Such explosions can account for ~1% of the Type Ia supernova rate., Comment: MNRAS Letters, in press

We analyze the planetary microlensing event MOA-2010-BLG-328. The best fit yields host and planetary masses of Mh = 0.11+/-0.01 M_{sun} and Mp = 9.2+/-2.2M_Earth, corresponding to a very late M dwarf and sub-Neptune-mass planet, respectively. The system lies at DL = 0.81 +/- 0.10 kpc with projected separation r = 0.92 +/- 0.16 AU. Because of the host's a-priori-unlikely close distance, as well as the unusual nature of the system, we consider the possibility that the microlens parallax signal, which determines the host mass and distance, is actually due to xallarap (source orbital motion) that is being misinterpreted as parallax. We show a result that favors the parallax solution, even given its close host distance. We show that future high-resolution astrometric measurements could decisively resolve the remaining ambiguity of these solutions., Comment: 30 pages, 6 figures. accepted for publication in ApJ. Figure 1 and 2 are updated

A planetary microlensing signal is generally characterized by a short-term perturbation to the standard single lensing light curve. A subset of binary-source events can produce perturbations that mimic planetary signals, thereby introducing an ambiguity between the planetary and binary-source interpretations. In this paper, we present analysis of the microlensing event MOA-2012-BLG-486, for which the light curve exhibits a short-lived perturbation. Routine modeling not considering data taken in different passbands yields a best-fit planetary model that is slightly preferred over the best-fit binary-source model. However, when allowed for a change in the color during the perturbation, we find that the binary-source model yields a significantly better fit and thus the degeneracy is clearly resolved. This event not only signifies the importance of considering various interpretations of short-term anomalies, but also demonstrates the importance of multi-band data for checking the possibility of false-positive planetary signals., Comment: 5 pages, 4 figures and ApJ submitted

Observations of accretion disks around young brown dwarfs have led to the speculation that they may form planetary systems similar to normal stars. While there have been several detections of planetary-mass objects around brown dwarfs (2MASS 1207-3932 and 2MASS 0441-2301), these companions have relatively large mass ratios and projected separations, suggesting that they formed in a manner analogous to stellar binaries. We present the discovery of a planetary-mass object orbiting a field brown dwarf via gravitational microlensing, OGLE-2012-BLG-0358Lb. The system is a low secondary/primary mass ratio (0.080 +- 0.001), relatively tightly-separated (~0.87 AU) binary composed of a planetary-mass object with 1.9 +- 0.2 Jupiter masses orbiting a brown dwarf with a mass 0.022 M_Sun. The relatively small mass ratio and separation suggest that the companion may have formed in a protoplanetary disk around the brown dwarf host, in a manner analogous to planets., Comment: 6 pages, 3 figures, 3 tables, Submitted in ApJ