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.

We present time series photometry in the field of Stock 8 and identified 130 short period variable stars. Twenty eight main-sequence and 23 pre-main-sequence variables are found to be part of cluster Stock 8. The main-sequence variables are classified as slow pulsator of the B type, $\beta$ Cep and $\delta$ Scuti stars. Fourteen main-sequence stars could be new class variables as discussed by Mowlavi et al. (2013) and Lata et al. (2014). The age and mass of pre-main-sequence variables are found to be $\lesssim$ 5 Myr and in the mass range of 0.5 to 2.8 M$_{\odot}$, respectively. These pre-main-sequence stars could be T-Tauri variables. We have found 5 and 2 of 23 PMS variables as classical T-Tauri stars and Herbig Ae/Be stars, respectively, whereas 16 PMS stars are classified as weak-line T Tauri stars.

We present measurements of the component masses in 15 Cataclysmic Variables (CVs) - 6 new estimates and 9 improved estimates. We provide new calibrations of the relationship between superhump period excess and mass ratio, and use this relation to estimate donor star masses for 225 superhumping CVs. With an increased sample of donor masses we revisit the implications for CV evolution. We confirm the high mass of white dwarfs in CVs, but find no trend in white dwarf mass with orbital period. We argue for a revision in the location of the orbital period minimum of CVs to $79.6 \pm 0.2$ min, significantly shorter than previous estimates. We find that CV donors below the gap have an intrinsic scatter of only 0.005 R$_{\odot}$ around a common evolutionary track, implying a correspondingly small variation in angular momentum loss rates. In contrast to prior studies, we find that standard CV evolutionary tracks - without additional angular momentum loss - are a reasonable fit to the donor masses just below the period gap, but that they do not reproduce the observed period minimum, or fit the donor radii below 0.1 M$_{\odot}$.

M dwarfs are prime targets in the hunt for habitable worlds around other stars. This is due to their abundance as well as their small radii and low masses and temperatures, which facilitate the detection of temperate, rocky planets in orbit around them. However, the fundamental properties of M dwarfs are difficult to constrain, often limiting our ability to characterise the planets they host. Here we test several theoretical relationships for M dwarfs by measuring 23 high precision, model-independent masses and radii for M dwarfs in binaries with white dwarfs. We find a large scatter in the radii of these low-mass stars, with 25 per cent having radii consistent with theoretical models while the rest are up to 12 per cent over-inflated. This scatter is seen in both partially- and fully-convective M dwarfs. No clear trend is seen between the over-inflation and age or metallicity, but there are indications that the radii of slowly rotating M dwarfs are more consistent with predictions, albeit with a similar amount of scatter in the measurements compared to more rapidly rotating M dwarfs. The sample of M dwarfs in close binaries with white dwarfs appears indistinguishable from other M dwarf samples, implying that common envelope evolution has a negligible impact on their structure. We conclude that theoretical and empirical mass-radius relationships lack the precision and accuracy required to measure the fundamental parameters of M dwarfs well enough to determine the internal structure and bulk composition of the planets they host.

Two narrow and dense rings (called C1R and C2R) were discovered around the Centaur object (10199) Chariklo during a stellar occultation observed on 2013 June 3. Following this discovery, we planned observations of several occultations by Chariklo's system in order to better characterize the physical properties of the ring and main body. Here, we use 12 successful occulations by Chariklo observed between 2014 and 2016. They provide ring profiles (physical width, opacity, edge structure) and constraints on the radii and pole position. Our new observations are currently consistent with the circular ring solution and pole position, to within the $\pm 3.3$ km formal uncertainty for the ring radii derived by Braga-Ribas et al. The six resolved C1R profiles reveal significant width variations from $\sim 5$ to 7.5 km. The width of the fainter ring C2R is less constrained, and may vary between 0.1 and 1 km. The inner and outer edges of C1R are consistent with infinitely sharp boundaries, with typical upper limits of one kilometer for the transition zone between the ring and empty space. No constraint on the sharpness of C2R's edges is available. A 1$\sigma$ upper limit of $\sim 20$ m is derived for the equivalent width of narrow (physical width <4 km) rings up to distances of 12,000 km, counted in the ring plane.

We present the second paper of a series of publications aiming at obtaining a better understanding regarding the nature of type Ia supernovae (SNIa) progenitors by studying a large sample of detached F, G and K main sequence stars in close orbits with white dwarf companions (i.e. WD+FGK binaries). We employ the LAMOST (Large Sky Area Multi-Object Fibre Spectroscopic Telescope) data release 4 spectroscopic data base together with GALEX (Galaxy Evolution Explorer) ultraviolet fluxes to identify 1,549 WD+FGK binary candidates (1,057 of which are new), thus doubling the number of known sources. We measure the radial velocities of 1,453 of these binaries from the available LAMOST spectra and/or from spectra obtained by us at a wide variety of different telescopes around the globe. The analysis of the radial velocity data allows us to identify 24 systems displaying more than 3sigma radial velocity variation that we classify as close binaries. We also discuss the fraction of close binaries among WD+FGK systems, which we find to be ~10 per cent, and demonstrate that high-resolution spectroscopy is required to efficiently identify double-degenerate SNIa progenitor candidates.

Using data from the extended Kepler mission in K2 Campaign 10 we identify two eclipsing binaries containing white dwarfs with cool companions that have extremely short orbital periods of only 71.2 min (SDSS J1205-0242, a.k.a. EPIC 201283111) and 72.5 min (SDSS J1231+0041, a.k.a. EPIC 248368963). Despite their short periods, both systems are detached with small, low-mass companions, in one case a brown dwarf, and the other case either a brown dwarf or a low-mass star. We present follow-up photometry and spectroscopy of both binaries, as well as phase-resolved spectroscopy of the brighter system, and use these data to place preliminary estimates on the physical and binary parameters. SDSS J1205-0242 is composed of a $0.39\pm0.02$M$_\odot$ helium-core white dwarf which is totally eclipsed by a $0.049\pm0.006$M$_\odot$ ($51\pm6$M$_J$) brown dwarf companion, while SDSS J1231+0041 is composed of a $0.56\pm0.07$M$_\odot$ white dwarf which is partially eclipsed by a companion of mass $\lesssim 0.095$M$_\odot$. In the case of SDSS J1205-0242 we look at the combined constraints from common-envelope evolution and brown dwarf models; the system is compatible with similar constraints from other post common-envelope binaries given the current parameter uncertainties, but has potential for future refinement.

We present high precision, model independent, mass and radius measurements for 16 white dwarfs in detached eclipsing binaries and combine these with previously published data to test the theoretical white dwarf mass-radius relationship. We reach a mean precision of 2.4 per cent in mass and 2.7 per cent in radius, with our best measurements reaching a precision of 0.3 per cent in mass and 0.5 per cent in radius. We find excellent agreement between the measured and predicted radii across a wide range of masses and temperatures. We also find the radii of all white dwarfs with masses less than 0.48M$_\odot$ to be fully consistent with helium core models, but they are on average 9 per cent larger than those of carbon-oxygen core models. In contrast, white dwarfs with masses larger than 0.52M$_\odot$ all have radii consistent with carbon-oxygen core models. Moreover, we find that all but one of the white dwarfs in our sample have radii consistent with possessing thick surface hydrogen envelopes ($10^{-5} \ge M_\mathrm{H}/M_\mathrm{WD} \ge 10^{-4}$), implying that the surface hydrogen layers of these white dwarfs are not obviously affected by common envelope evolution.

We present results derived from four stellar occultations by the plutino object (208996) 2003~AZ$_{84}$, detected at January 8, 2011 (single-chord event), February 3, 2012 (multi-chord), December 2, 2013 (single-chord) and November 15, 2014 (multi-chord). Our observations rule out an oblate spheroid solution for 2003~AZ$_{84}$'s shape. Instead, assuming hydrostatic equilibrium, we find that a Jacobi triaxial solution with semi axes $(470 \pm 20) \times (383 \pm 10) \times (245 \pm 8)$~km % axis ratios $b/a= 0.82 \pm 0.05$ and $c/a= 0.52 \pm 0.02$, can better account for all our occultation observations. Combining these dimensions with the rotation period of the body (6.75~h) and the amplitude of its rotation light curve, we derive a density $\rho=0.87 \pm 0.01$~g~cm$^{-3}$ a geometric albedo $p_V= 0.097 \pm 0.009$. A grazing chord observed during the 2014 occultation reveals a topographic feature along 2003~AZ$_{84}$'s limb, that can be interpreted as an abrupt chasm of width $\sim 23$~km and depth $> 8$~km or a smooth depression of width $\sim 80$~km and depth $\sim 13$~km (or an intermediate feature between those two extremes).

We present the discovery of the hot subdwarf B star (sdB) binary PTF1 J082340.04+081936.5. The system has an orbital period P$_{\rm orb}=87.49668(1)$ min (0.060761584(10) days), making it the second-most compact sdB binary known. The lightcurve shows ellipsoidal variations. Under the assumption that the sdB primary is synchronized with the orbit, we find a mass $M_{\rm sdB}=0.45^{+0.09}_{-0.07}$ M$_\odot$, a companion white dwarf mass $M_{\rm WD}=0.46^{+0.12}_{-0.09}$ M$_\odot$ and a mass ratio $q = \frac{M_{\rm WD}}{M_{\rm sdB}}=1.03^{+0.10}_{-0.08}$. The future evolution was calculated using the MESA stellar evolution code. Adopting a canonical sdB mass of $M_{\rm sdB}=0.47$ M$_\odot$, we find that the sdB still burns helium at the time it will fill its Roche lobe if the orbital period was less than 106 min at the exit from the last common envelope phase. For longer common envelope exit periods the sdB will have stopped burning helium and turned into a C/O white dwarf at the time of contact. Comparing the spectroscopically derived log(g) and $T_{\rm eff}$ with our MESA models, we find that an sdB model with a hydrogen envelope mass of $5\times10^{-4} M_\odot$ matches the measurements at a post-common envelope age of 94 Myr, corresponding to a post-common envelope orbital period of 109 min which is close to the limit to start accretion while the sdB is still burning helium.

We present new time-resolved photometry of 74 cataclysmic variables (CVs), 47 of which are eclipsing. 13 of these eclipsing systems are newly discovered. For all 47 eclipsing systems we show high cadence (1-20 seconds) light curves obtained with the high-speed cameras ultracam and ultraspec. We provide new or refined ephemerides, and supply mid-eclipse times for all observed eclipses. We assess the potential for light curve modelling of all 47 eclipsing systems to determine their system parameters, finding 20 systems which appear to be suitable for future study.

We present a multiwavelength investigation of the young stellar population and star formation activities around the \hii region Sharpless 311. Using our deep near-infrared observations and archival {\it Spitzer}-IRAC observations, we have detected a total of 125 young stellar objects (YSOs) in an area of $\sim$86 arcmin$^2$. The YSOs sample include 8 Class I and 117 Class II candidate YSOs. The mass completeness of the identified YSOs sample is estimated to be 1.0 \msun. The ages and masses of the majority of the candidate YSOs are estimated to be in the range of $\sim$0.1$-$5 Myr and $\sim$0.3$-$6 \msun, respectively. The 8 \mum image of S311 displays an approximately spherical cavity around the ionizing source which is possibly created due to the expansion of the \hii region. The spatial distribution of the candidate YSOs reveals that a significant number of them are distributed systematically along the 8 $\mu$m emission with a majority clustered around the eastern border of the \hii region. Four clumps/compact \hii regions are detected in the radio continuum observations at 1280 MHz, which might have been formed during the expansion of the \hii region. The estimated dynamical age of the region, main-sequence lifetime of the ionizing source, the spatial distribution and ages of the candidate YSOs indicate triggered star formation in the complex.

We present a long-term programme for timing the eclipses of white dwarfs in close binaries to measure apparent and/or real variations in their orbital periods. Our programme includes 67 close binaries, both detached and semi-detached and with M-dwarfs, K-dwarfs, brown dwarfs or white dwarfs secondaries. In total, we have observed more than 650 white dwarf eclipses. We use this sample to search for orbital period variations and aim to identify the underlying cause of these variations. We find that the probability of observing orbital period variations increases significantly with the observational baseline. In particular, all binaries with baselines exceeding 10 yrs, with secondaries of spectral type K2 -- M5.5, show variations in the eclipse arrival times that in most cases amount to several minutes. In addition, among those with baselines shorter than 10 yrs, binaries with late spectral type (>M6), brown dwarf or white dwarf secondaries appear to show no orbital period variations. This is in agreement with the so-called Applegate mechanism, which proposes that magnetic cycles in the secondary stars can drive variability in the binary orbits. We also present new eclipse times of NN Ser, which are still compatible with the previously published circumbinary planetary system model, although only with the addition of a quadratic term to the ephemeris. Finally, we conclude that we are limited by the relatively short observational baseline for many of the binaries in the eclipse timing programme, and therefore cannot yet draw robust conclusions about the cause of orbital period variations in evolved, white dwarf binaries.

We report further results from the program of lunar occultation (LO) observations started at the 2.4-m Thai National Telescope (TNT) in 2014. We have recorded LO events of 18 stellar sources, leading to the detection of four angular diameters and two binary stars. With two exceptions, these are first-time determinations. We could resolve angular diameters as small as 2 milliarcseconds (mas) and projected separations as small as 4mas. We discuss the individual results, in the context of previous observations when available. The first-time angular diameters for o Psc, HR 6196 and 75 Leo are in good agreement with expected values, while that of Pi Leo agrees with the average of previous determinations but has a higher accuracy. We find a new secondary in o Psc, as previously suspected from Hipparcos data. We also obtain an accuratemeasurement of the companion in 31 Ari, revealing inconsistencies in the currently available orbital parameters. The TNT, equipped with the fast ULTRASPEC imager, is the leading facility in Southeast Asia for high time resolution observations. The LO technique at this telescope achieves a sensitivity of i up to 10mag, with a potential to detect several hundreds of LO events per year.

SDSS J1021+1744 is a detached, eclipsing white dwarf / M dwarf binary discovered in the Sloan Digital Sky Survey. Outside the primary eclipse, the light curves of such systems are usually smooth and characterised by low-level variations caused by tidal distortion and heating of the M star component. Early data on SDSS J1021+1744 obtained in June 2012 was unusual in showing a dip in flux of uncertain origin shortly after the white dwarf's eclipse. Here we present high-time resolution, multi-wavelength observations of 35 more eclipses over 1.3 years, showing that the dip has a lifetime extending over many orbits. Moreover the "dip" is in fact a series of dips that vary in depth, number and position, although they are always placed in the phase interval 1.06 to 1.26 after the white dwarf's eclipse, near the L5 point in this system. Since SDSS J1021+1744 is a detached binary, it follows that the dips are caused by the transit of the white dwarf by material around the Lagrangian L5 point. A possible interpretation is that they are the signatures of prominences, a phenomenon already known from H-alpha observations of rapidly rotating single stars as well as binaries. What makes SDSS J1021+1744 peculiar is that the material is dense enough to block continuum light. The dips appear to have finally faded out around 2015 May after the first detection by Parsons et al. in 2012, suggesting a lifetime of years.

We study the mid-egress eclipse timing data gathered for the cataclysmic binary HU Aquarii during the years 1993-2014. The (O-C) residuals were previously attributed to a single ~7 Jupiter mass companion in ~5 au orbit or to a stable 2-planet system with an unconstrained outermost orbit. We present 22 new observations gathered between June, 2011 and July, 2014 with four instruments around the world. They reveal a systematic deviation of ~60 - 120 seconds from the older ephemeris. We re-analyse the whole set of the timing data available. Our results provide an erratum to the previous HU Aqr planetary models, indicating that the hypothesis for a third and fourth body in this system is uncertain. The dynamical stability criterion and a particular geometry of orbits rule out coplanar 2-planet configurations. A putative HU Aqr planetary system may be more complex, e.g., highly non-coplanar. Indeed, we found examples of 3-planet configurations with the middle planet in a retrograde orbit, which are stable for at least 1Gyr, and consistent with the observations. The (O-C) may be also driven by oscillations of the gravitational quadrupole moment of the secondary, as predicted by the Lanza et al. modification of the Applegate mechanism. Further systematic, long-term monitoring of HU Aqr is required to interpret the (O-C) residuals.

We present high time resolution SDSS-$g'$ and SDSS-$z'$ light curves of the primary eclipse in SDSS J141126.20+200911.1, together with time-resolved X-Shooter spectroscopy and near-infrared $JHK_{s}$ photometry. Our observations confirm the substellar nature of the companion, making SDSS J141126.20+200911.1 the first eclipsing white dwarf/brown dwarf binary known. We measure a (white dwarf model dependent) mass and radius for the brown dwarf companion of $M_{2} = 0.050 \pm 0.002$ $M_{\odot}$ and $R_{2} = 0.072 \pm 0.004$ $M_{\odot}$, respectively. The lack of a robust detection of the companion light in the $z'$-band eclipse constrains the spectral type of the companion to be later than L5. Comparing the NIR photometry to the expected white dwarf flux reveals a clear $K_s$-band excess, suggesting a spectral type in the range L7-T1. The radius measurement is consistent with the predictions of evolutionary models, and suggests a system age in excess of three Gyr. The low companion mass is inconsistent with the inferred spectral type of L7-T1, instead predicting a spectral type nearer T5. This indicates that irradiation of the companion in SDSS J1411 could be causing a significant temperature increase, at least on one hemisphere.

ULTRASPEC is a high-speed imaging photometer mounted permanently at one of the Nasmyth focii of the 2.4-m Thai National Telescope (TNT) on Doi Inthanon, Thailand's highest mountain. ULTRASPEC employs a 1024x1024 pixel frame-transfer, electron-multiplying CCD (EMCCD) in conjunction with re-imaging optics to image a field of 7.7'x7.7' at (windowed) frame rates of up to ~200 Hz. The EMCCD has two outputs - a normal output that provides a readout noise of 2.3 e- and an avalanche output that can provide essentially zero readout noise. A six-position filter wheel enables narrow-band and broad-band imaging over the wavelength range 330-1000 nm. The instrument saw first light on the TNT in November 2013 and will be used to study rapid variability in the Universe. In this paper we describe the scientific motivation behind ULTRASPEC, present an outline of its design and report on its measured performance on the TNT.

The recently inaugurated 2.4\,m Thai National Telescope (TNT) is equipped, among other instruments, with the ULTRASPEC low-noise, frame-transfer EMCCD camera. At the end of its first official observing season, we report on the use of this facility to record high time resolution imaging using small detector subarrays with sampling as fast as several $10^2$\,Hz. In particular, we have recorded lunar occultations of several stars which represent the first contribution to this area of research made from South-East Asia with a telescope of this class. Among the results, we discuss an accurate measurement of $\alpha $~Cnc, which has been reported previously as a suspected close binary. Attempts to resolve this star by several authors have so far met with a lack of unambiguous confirmation. With our observation we are able to place stringent limits on the projected angular separation ($<0\farcs003$) and brightness ($\Delta{\rm m} > 5$) of a putative companion. We also present a measurement of the binary {HR~7072}, which extends considerably the time coverage available for its yet undetermined orbit. We discuss our precise determination of the flux ratio and projected separation in the context of other available data. We conclude by providing an estimate of the performance of ULTRASPEC at TNT for lunar occultation work. This facility can help to extend the lunar occultation technique in a geographical area where no comparable resources were available until now.