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.

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.

YZ Phe is a very short-period contact binary (Sp.= $K2\,V$) with an orbital period of 0.2347 days near the short period limit (0.22 d). We present the complete light curves in $VRI$ bands, which photometric data were obtained with the 0.61-m telescope of PROMPT-8 at CTIO in Chile during June to October 2016 and August 2017. The photometric solutions were determined by using the W-D method and the results reveal that YZ Phe is a W-subtype shallow contact binary ($f\sim$ 10%, $q$ = 2.635 or $1/q$ = 0.379 for W subtype) with rotational motion of a large hot spot on the more massive component, showing a strong O'Connell effect with variation of maxima in photometric time series at period of 4.20 yr and stellar cycle at period of 1.28 yr. By compiling all available eclipse times, the result shows a long-term period decrease at a rate of $\mathrm{d}P/\mathrm{d}t = -2.64(\pm 0.02)\times 10^{-8}$ d $yr^{-1}$, superimposed on a cyclic variation ($A_3$ = 0.0081 days and $P_3$ = 40.76 years). This variation cannot be explained by Applegate mechanism. Thus, the cyclic change may be interpreted as light-travel time effect via the presence of a cool third body. Based on photometric solutions, the third light was detected with 2% contribution of total light in $V$ and $I$ bands. Those support the existence of a third body. For the long-term period decrease, it can be explained by mass transfer from the more massive component ($M_2 \sim 0.74 M_{\odot}$) to the less massive one ($M_1 \sim 0.28 M_{\odot}$) or plus AML via magnetic braking. With $1/q$ $<$ 0.4 and long-term period decrease, all suggest that YZ Phe is on the AML-controlled state and its fill-out factor will increase, as well as the system will evolve into a deeper normal contact binary.

TY Pup is a well-known bright eclipsing binary in southern hemisphere with an orbital period of 0.8192 days. New light curves in $B, V, (RI)_C$ bands were obtained with the 0.61-m reflector robotic telescope (PROMPT-8) at CTIO in Chile from January to February 2015 and from March to April 2017. By analyzing those photometric data with the W-D method, it is found that TY Pup is a low-mass-ratio ($q \sim$ 0.184) and deep contact binary with a high fill-out factor ($84.3\,\%$). An investigation of all available times of minimum light including three new ones obtained with the 60-cm and the 1.0-m telescopes at Yunnan Observatories in China indicates that the period change of TY Pup is complex. An upward parabolic variation in the $O-C$ diagram is detected to be superimposed on a cyclic oscillation. The upward parabolic change reveals a long-term continuous increase in the orbital period at a rate of $\mathrm{d}P/\mathrm{d}t = 5.57(\pm 0.08)\times10^{-8}$ d $\textrm{yr}^{-1}$. The period increase can be explained by mass transfer from the less massive component ($M_2 \sim 0.3 M_{\odot}$) to the more massive one ($M_1 \sim 1.65 M_{\odot}$) and the mass ratio of the binary system will become more smaller. In this way, the binary will be merging when it meets the criterion that the orbital angular momentum is less than 3 times the total spin angular momentum, i.e., $J_{orb} < 3J_{rot}$. This suggests that the system will finally merge into a rapid-rotating single star and may produce a luminous red nova. The cyclic oscillation in the $O-C$ diagram can be interpreted by the light-travel time effect (LITE) via the presence of an additional companion.

New CCD photometric light curves of short period (P=0.285d) eclipsing binary RW Dor are presented. The observations performed with the PROMPT-8 robotic telescope at CTIO in Chile from March 2015 to March 2017. The other eclipse timings were obtained from the 2.15-m JS telescope at CASLEO, San Juan, Argentina in December 2011. By light-curve analysis, it is found that RW Dor is a W-type shallow contact binary with a fill-out factor $f \sim 11\%$ and high mass ratio $q \sim 1.587$ (1/q = 0.63), where the hotter component is the less massive one ($M_1 \sim 0.52M_{\odot}$ and $M_2 \sim 0.82M_{\odot}$). For orbital period investigation, the new fifteen eclipse times and those in previous published were compiled. Based on $O-C$ analysis with very weak evidence suggests that a long-term period decrease with a rate of $\mathrm{d}P/\mathrm{d}t = -9.61\times10^{-9}$ d $\textrm{yr}^{-1}$ is superimposed on a cyclic variation ($A_3$ = 0.0054 days and $P_3$ = 49.9 yrs). The long-term period decrease can be interpreted as mass transfer from the more massive component to the less massive one or combine with the angular momentum loss (AML) via magnetic braking. In addition, with the marginal contact phase, high mass ratio (1/q $>$ 0.4) and the long-term period decrease, all suggest that RW Dor is a newly formed contact binary via a Case A mass transfer and it will evolve into a deeper normal contact binary. If the cyclic change is correct, the light-travel time effect via the presence of a cool third body will be more plausible to explain for this.

GU Ori was observed with the 1m telescope at Yunnan Observatories in 2005. To determine its physical properties, the Wilson-Devinney program is used. The results reveal that GU Ori is a W-subtype shallow contact binary with a more massive but cooler star 2. The mass of its two component stars are estimated to be $M_1 = 0.45M_\odot$, $M_2 = 1.05M_\odot$. The O'Connell effect was reported to be negative on the light curves observed in 2005. However, it changed to a positive one on the light curves observed from 2011 to 2012. The mean surface temperatures of star 2 ($T_2$) determined by the two sets of light curves were different, which may result from stellar activities. The O - C diagram shows that the period of GU Ori is decreasing at a rate of $dP/dt=-6.24\times{10^{-8}}day\cdot year^{-1}$, which may be caused by mass transfer from star 2 to star 1 with a rate of $\frac{dM_{2}}{dt}= - 2.98\times{10^{-8}}M_\odot/year$. GU Ori is a contact binary with quite high metallicity.

We present a timing analysis of the eclipsing post-common envelope binary (PCEB) DE CVn. Based on new CCD photometric observations and the published data, we found that the orbital period in DE CVn has a cyclic period oscillation with an amplitude of $28.08$ s and a period of $11.22$ years plus a rapid period decrease at a rate of $\dot{P}=-3.35\times10^{-11}ss^{-1}$. According to the evolutionary theory, secular period decreases in PCEBs arise from angular momentum losses (AMLs) driven by gravitational radiation (GR) and magnetic braking (MB). However, the observed orbital decay is too fast to be produced by AMLs via GR and MB, indicating that there could be other AML mechanism. We suggest that a circumbinary disk around DE CVn may be responsible for the additional AML. The disk mass was derived as a few$\times$$10^{-4}$-$10^{-3}$$M_{\odot}$ , which is in agreement with that inferred from previous studies in the order of magnitude. The cyclic change is most likely result of the gravitational perturbation by a circumbinary object due to the Applegate's mechanism fails to explain such a large period oscillation. The mass of the potential third body is calculated as $M_{3}\sin{i'}=0.011(\pm0.003)M_{\odot}$. Supposing the circumbinary companion and the eclipsing binary is coplanar, its mass would correspond to a giant planet. This hypothetical giant planet is moving in a circular orbit of radius $\sim5.75(\pm2.02)$ AU around its host star.

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.

GJ3470b is a hot Neptune exoplanet orbiting an M dwarf and the first sub-Jovian planet to exhibit Rayleigh scattering. We present transit timing variation (TTV) and transmission spectroscopy analyses of multi-wavelength optical photometry from 2.4-m and 0.5-m telescopes at the Thai National Observatory, and the 0.6-m PROMPT-8 telescope in Chile. Our TTV analysis allows us to place an upper mass limit for a second planet in the system. The presence of a hot Jupiter with a period of less than 10 days or a planet with an orbital period between 2.5 and 4.0 days are excluded. Combined optical and near-infrared transmission spectroscopy favour an H/He dominated haze (mean molecular weight 1.08$\pm$0.20) with high particle abundance at high altitude. We also argue that previous near-infrared data favour the presence of methane in the atmosphere of GJ3470b.

We present time-resolved optical photometry of the binary millisecond `redback' pulsar PSR J1023+0038 (=AY Sex) during its low-mass X-ray binary phase. The light curves taken between 2014 January and April show an underlying sinusoidal modulation due to the irradiated secondary star and accretion disc. We also observe superimposed rapid flaring on time-scales as short as ~20 s with amplitudes of ~0.1-0.5 mag and additional large flare events on time-scales of ~5-60 min with amplitudes ~0.5-1.0 mag. The power density spectrum of the optical flare light curves is dominated by a red-noise component, typical of aperiodic activity in X-ray binaries. Simultaneous X-ray and UV observations by the Swift satellite reveal strong correlations that are consistent with X-ray reprocessing of the UV light, most likely in the outer regions of the accretion disc. On some nights we also observe sharp-edged, rectangular, flat-bottomed dips randomly distributed in orbital phase, with a median duration of ~250 s and a median ingress/egress time of ~20 s. These rectangular dips are similar to the mode-switching behaviour between disc `active' and `passive' luminosity states, observed in the X-ray light curves of other redback millisecond pulsars. This is the first time that the optical analogue of the X-ray mode-switching has been observed. The properties of the passive and active state light curves can be explained in terms of clumpy accretion from a trapped inner accretion disc near the corotation radius, resulting in rectangular, flat-bottomed optical and X-ray light curves.

The final assembly of terrestrial planets occurs via massive collisions, which can launch copious clouds of dust that are warmed by the star and glow in the infrared. We report the real-time detection of a debris-producing impact in the terrestrial planet zone around a 35-million year-old solar analog star. We observed a substantial brightening of the debris disk at 3-5 {\mu}m, followed by a decay over a year, with quasi-periodic modulations of the disk flux. The behavior is consistent with the occurrence of a violent impact that produced vapor out of which a thick cloud of silicate spherules condensed that were ground into dust by collisions. These results demonstrate how the time domain can become a new dimension for the study of terrestrial planet formation.

Luminous debris disks of warm dust in the terrestrial planet zones around solar-like stars are recently found to vary, indicative of ongoing large-scale collisions of rocky objects. We use Spitzer 3.6 and 4.5 {\mu}m time-series observations in 2012 and 2013 (extended to 2014 in one case) to monitor 5 more debris disks with unusually high fractional luminosities ("extreme debris disk"), including P1121 in the open cluster M47 (80 Myr), HD 15407A in the AB Dor moving group (80 Myr), HD 23514 in the Pleiades (120 Myr), HD 145263 in the Upper Sco Association (10 Myr), and the field star BD+20 307 (>1 Gyr). Together with the published results for ID8 in NGC 2547 (35 Myr), this makes the first systematic time-domain investigation of planetary impacts outside the solar system. Significant variations with timescales shorter than a year are detected in five out of the six extreme debris disks we have monitored. However, different systems show diverse sets of characteristics in the time domain, including long-term decay or growth, disk temperature variations, and possible periodicity.

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.