GRANDMA is a network of 25 telescopes of different sizes, including both photometric and spectroscopic facilities. The network aims to coordinate follow-up observations of gravitational-wave candidate alerts, especially those with large localisation uncertainties, to reduce the delay between the initial detection and the optical confirmation. In this paper, we detail GRANDMA's observational performance during Advanced LIGO/Advanced Virgo Observing Run 3 (O3), focusing on the second part of O3; this includes summary statistics pertaining to coverage and possible astrophysical origin of the candidates. To do so, we quantify our observation efficiency in terms of delay between gravitational-wave candidate trigger time, observations, and the total coverage. Using an optimised and robust coordination system, GRANDMA followed-up about 90 % of the gravitational-wave candidate alerts, i.e. 49 out of 56 candidates. This led to coverage of over 9000 deg2 during O3. The delay between the gravitational-wave candidate trigger and the first observation was below 1.5 hour for 50 % of the alerts. We did not detect any electromagnetic counterparts to the gravitational-wave candidates during O3, likely due to the very large localisation areas (on average thousands of degrees squares) and relatively large distance of the candidates (above 200 Mpc for 60 % of BNS candidates). We derive constraints on potential kilonova properties for two potential binary neutron star coalescences (GW190425 and S200213t), assuming that the events' locations were imaged.

Galactic microlensing datasets now comprise in excess of $10^4$ events, and with the advent of next generation microlensing surveys that may be undertaken with facilities such as the Vera Rubin Observatory (formerly LSST) and Nancy Grace Roman Space Telescope (formerly WFIRST), this number will increase significantly. So too will the fraction of events with measurable higher order information such as finite source effects and lens-source relative proper motion. Analysing such data requires a more sophisticated Galactic microlens modeling approach. We present a new second-generation Manchester-Besan\c{c}on Microlensing Simulator (MaB$\mu$lS-2), which uses a version of the Besan\c{c}on population synthesis Galactic model that provides good agreement with stellar kinematics observed by HST towards the bulge. MaB$\mu$lS-2 provides high-fidelity signal-to-noise limited maps of the microlensing optical depth, rate and average timescale towards a 400 sq. degree region of the Galactic bulge in several optical to near-infrared pass-bands. The maps take full account of the unresolved stellar background as well as limb-darkened source profiles. Comparing MaB$\mu$lS-2 to the efficiency-corrected OGLE-IV 8,000 event sample shows a much improved agreement over the previous version of MaB$\mu$lS, and succeeds in matching even small-scale structural features in the OGLE-IV event rate map. However, there remains evidence for a small under-prediction in the event rate per source and over-prediction in timescale. MaB$\mu$lS-2 is available online (\<www.mabuls.net>) to provide on-the-fly maps for user supplied cuts in survey magnitude, event timescale and relative proper motion.

We present an optical-to-infrared transmission spectrum of the inflated sub-Saturn-mass exoplanet KELT-11b measured with the Transiting Exoplanet Survey Satellite (TESS), the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) G141 spectroscopic grism, and the Spitzer Space Telescope (Spitzer) at 3.6 $\mu$m, in addition to a Spitzer 4.5 $\mu$m secondary eclipse. The precise HST transmission spectrum notably reveals a low-amplitude water feature with an unusual shape. We apply a suite of modeling tools to the transmission spectrum to investigate the planet's properties. Based on a retrieval analysis with varying molecular abundances, we find strong evidence for water absorption in the spectrum and tentative evidence for other absorbers (HCN, TiO, and AlO) depending on model assumptions. The retrieved water abundance is generally $\lesssim 0.1\times$ solar (0.001--0.7$\times$ solar over a range of model assumptions), which is several orders of magnitude lower than expected from planet formation models based on the solar system metallicity trend. We also consider chemical equilibrium and self-consistent 1D radiative-convective equilibrium model fits and find that they too prefer low metallicities ($[M/H] \lesssim -2$, consistent with the free retrieval results); however, the equilibrium models fit poorly and should be interpreted with caution. Finally, we find that the dayside flux measured from the Spitzer secondary eclipse is indicative of full heat redistribution from KELT-11b's dayside to nightside, assuming the dayside is not cloudy. Altogether, the unusual results for KELT-11b's composition are suggestive of new challenges on the horizon for atmosphere and formation models in the face of increasingly precise measurements of exoplanet spectra.

One of the principal bottlenecks to atmosphere characterisation in the era of all-sky surveys is the availability of fast, autonomous and robust atmospheric retrieval methods. We present a new approach using unsupervised machine learning to generate informed priors for retrieval of exoplanetary atmosphere parameters from transmission spectra. We use principal component analysis (PCA) to efficiently compress the information content of a library of transmission spectra forward models generated using the PLATON package. We then apply a $k$-means clustering algorithm in PCA space to segregate the library into discrete classes. We show that our classifier is almost always able to instantaneously place a previously unseen spectrum into the correct class, for low-to-moderate spectral resolutions, $R$, in the range $R~=~30-300$ and noise levels up to $10$~per~cent of the peak-to-trough spectrum amplitude. The distribution of physical parameters for all members of the class therefore provides an informed prior for standard retrieval methods such as nested sampling. We benchmark our informed-prior approach against a standard uniform-prior nested sampler, finding that our approach is up to a factor two faster, with negligible reduction in accuracy. We demonstrate the application of this method to existing and near-future observatories, and show that it is suitable for real-world application. Our general approach is not specific to transmission spectroscopy and should be more widely applicable to cases that involve repetitive fitting of trusted high-dimensional models to large data catalogues, including beyond exoplanetary science.

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 report the analysis of additional multiband photometry and spectroscopy and new adaptive optics (AO) imaging of the nearby planetary microlensing event TCP J05074264+2447555 (Kojima-1), which was discovered toward the Galactic anticenter in 2017 (Nucita et al.). We confirm the planetary nature of the light-curve anomaly around the peak while finding no additional planetary feature in this event. We also confirm the presence of apparent blending flux and the absence of significant parallax signal reported in the literature. The AO image reveals no contaminating sources, making it most likely that the blending flux comes from the lens star. The measured multiband lens flux, combined with a constraint from the microlensing model, allows us to narrow down the previously unresolved mass and distance of the lens system. We find that the primary lens is a dwarf on the K/M boundary (0.581 \pm 0.033 M_sun) located at 505 \pm 47 pc and the companion (Kojima-1Lb) is a Neptune-mass planet (20.0 \pm 2.0 M_earth) with a semi-major axis of 1.08 ^{+0.62}_{-0.18} au. This orbit is a few times smaller than those of typical microlensing planets and is comparable to the snow-line location at young ages. We calculate that the a priori detection probability of Kojima-1Lb is only \sim 35%, which may imply that Neptunes are common around the snow line, as recently suggested by the transit and radial velocity techniques. The host star is the brightest among the microlensing planetary systems (Ks = 13.7), offering a great opportunity to spectroscopically characterize this system, even with current facilities.

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.

We present the discovery of KELT-24 b, a massive hot Jupiter orbiting a bright (V=8.3 mag, K=7.2 mag) young F-star with a period of 5.6 days. The host star, KELT-24 (HD 93148), has a $T_{\rm eff}$ =$6509^{+50}_{-49}$ K, a mass of $M_{*}$ = $1.460^{+0.055}_{-0.059}$ $M_{\odot}$, radius of $R_{*}$ = $1.506\pm0.022$ $R_{\odot}$, and an age of $0.78^{+0.61}_{-0.42}$ Gyr. Its planetary companion (KELT-24 b) has a radius of $R_{\rm P}$ = $1.272\pm0.021$ $R_{\rm J}$, a mass of $M_{\rm P}$ = $5.18^{+0.21}_{-0.22}$ $M_{\rm J}$, and from Doppler tomographic observations, we find that the planet's orbit is well-aligned to its host star's projected spin axis ($\lambda$ = $2.6^{+5.1}_{-3.6}$). The young age estimated for KELT-24 suggests that it only recently started to evolve from the zero-age main sequence. KELT-24 is the brightest star known to host a transiting giant planet with a period between 5 and 10 days. Although the circularization timescale is much longer than the age of the system, we do not detect a large eccentricity or significant misalignment that is expected from dynamical migration. The brightness of its host star and its moderate surface gravity make KELT-24b an intriguing target for detailed atmospheric characterization through spectroscopic emission measurements since it would bridge the current literature results that have primarily focused on lower mass hot Jupiters and a few brown dwarfs.

The large number of new planets expected from wide-area transit surveys means that follow-up transmission spectroscopy studies of their atmospheres will be limited by the availability of telescope assets. We argue that telescopes covering a broad range of apertures will be required, with even 1m-class instruments providing a potentially important contribution. Survey strategies that employ automated target selection will enable robust population studies. As part of such a strategy, we propose a decision metric to pair the best target to the most suitable telescope, and demonstrate its effectiveness even when only primary transit observables are available. Transmission spectroscopy target selection need not therefore be impeded by the bottle-neck of requiring prior follow-up observations to determine the planet mass. The decision metric can be easily deployed within a distributed heterogeneous network of telescopes equipped to undertake either broadband photometry or spectroscopy. We show how the metric can be used either to optimise the observing strategy for a given telescope (e.g. choice of filter) or to enable the selection of the best telescope to optimise the overall sample size. Our decision metric can also provide the basis for a selection function to help evaluate the statistical completeness of follow-up transmission spectroscopy datasets. Finally, we validate our metric by comparing its ranked set of targets against lists of planets that have had their atmospheres successfully probed, and against some existing prioritised exoplanet lists.

GY Cnc is a deeply eclipsing cataclysmic variable star with an orbital period of 4.21 hours that has shown several dwarf nova outbursts. The variable was continuously observed by the K2/Kepler satellite with a short cadence in Campaign 5 (C05) for 75 days during 2015. The star was again observed in 2017/2018 for 80 consecutive days during Campaign 16 (C16). 419 well-observed eclipses were measured over C5 and 446 timings were determined in C16. A new ephemeris was calculated combining the K2 data, previously published timings, new light curves from the KELT Follow-Up Network, and additional observations in the AAVSO database. We have refined the orbital period of GY Cnc and improved its ephemeris which had accumulated an error of about 300s over 18 years. A quadratic term was not found to be significant indicating that there is currently no detectable orbital period derivative. We observed a correlation between the quiescent system brightness and the eclipse timing residuals that are as large as +/-15s in the K2 data.

The Kilodegree Extremely Little Telescope (KELT) project has been conducting a photometric survey for transiting planets orbiting bright stars for over ten years. The KELT images have a pixel scale of ~23"/pixel---very similar to that of NASA's Transiting Exoplanet Survey Satellite (TESS)---as well as a large point spread function, and the KELT reduction pipeline uses a weighted photometric aperture with radius 3'. At this angular scale, multiple stars are typically blended in the photometric apertures. In order to identify false positives and confirm transiting exoplanets, we have assembled a follow-up network (KELT-FUN) to conduct imaging with higher spatial resolution, cadence, and photometric precision than the KELT telescopes, as well as spectroscopic observations of the candidate host stars. The KELT-FUN team has followed-up over 1,600 planet candidates since 2011, resulting in more than 20 planet discoveries. Excluding ~450 false alarms of non-astrophysical origin (i.e., instrumental noise or systematics), we present an all-sky catalog of the 1,128 bright stars (6<V<10) that show transit-like features in the KELT light curves, but which were subsequently determined to be astrophysical false positives (FPs) after photometric and/or spectroscopic follow-up observations. The KELT-FUN team continues to pursue KELT and other planet candidates and will eventually follow up certain classes of TESS candidates. The KELT FP catalog will help minimize the duplication of follow-up observations by current and future transit surveys such as TESS.

We present 45 ground-based photometric observations of the K2-22 system collected between December 2016 and May 2017, which we use to investigate the evolution of the transit of the disintegrating planet K2-22b. Last observed in early 2015, in these new observations we recover the transit at multiple epochs and measure a typical depth of <1.5%. We find that the distribution of our measured transit depths is comparable to the range of depths measured in observations from 2014 and 2015. These new observations also support ongoing variability in the K2-22b transit shape and time, although the overall shallowness of the transit makes a detailed analysis of these transit parameters difficult. We find no strong evidence of wavelength-dependent transit depths for epochs where we have simultaneous coverage at multiple wavelengths, although our stacked Las Cumbres Observatory data collected over days-to-months timescales are suggestive of a deeper transit at blue wavelengths. We encourage continued high-precision photometric and spectroscopic monitoring of this system in order to further constrain the evolution timescale and to aid comparative studies with the other few known disintegrating planets.

Proxima Centauri is known as the closest star from the Sun. Recently, radial velocity observations revealed the existence of an Earth-mass planet around it. With an orbital period of ~11 days, the surface of Proxima Centauri b is temperate and might be habitable. We took a photometric monitoring campaign to search for its transit, using the Bright Star Survey Telescope at the Zhongshan Station in Antarctica. A transit-like signal appearing on 2016 September 8th, is identified tentatively. Its midtime, $T_{C}=2,457,640.1990\pm0.0017$ HJD, is consistent with the predicted ephemeris based on RV orbit in a 1$\sigma$ confidence interval. Time-correlated noise is pronounced in the light curve of Proxima Centauri, affecting detection of transits. We develop a technique, in a Gaussian process framework, to gauge the statistical significance of potential transit detection. The tentative transit signal reported here, has a confidence level of $2.5\sigma$. Further detection of its periodic signals is necessary to confirm the planetary transit of Proxima Centauri b. We plan to monitor Proxima Centauri in next Polar night at Dome A in Antarctica, taking the advantage of continuous darkness. \citet{Kipping17} reported two tentative transit-like signals of Proxima Centauri b, observed by the Microvariability and Oscillation of Stars space Telescope in 2014 and 2015, respectively. The midtransit time of our detection is 138 minutes later than that predicted by their transit ephemeris. If all the signals are real transits, the misalignment of the epochs plausibly suggests transit timing variations of Proxima Centauri b induced by an outer planet in this system.

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.

Galactic bulge microlensing surveys provide a probe of Galactic structure. We present the first field-by-field comparison between microlensing observations and the Besan\c{c}on population synthesis Galactic model. Using an updated version of the model we provide maps of optical depth, average event duration and event rate for resolved source populations and for difference imaging (DIA) events. We also compare the predicted event timescale distribution to that observed. The simulation follows the selection criteria of the MOA-II survey (Sumi et al. 2013). We modify the Besan\c{c}on model to include M dwarfs and brown dwarfs. Our best fit model requires a brown dwarf mass function slope of $-0.4$. The model provides good agreement with the observed average duration, and respectable consistency with the shape of the timescale distribution (reduced $\chi^2 \simeq 2.2$). The DIA and resolved source limiting yields bracket the observed number of events by MOA-II ($2.17\times$ and $0.83\times$ the number observed, respectively). We perform a 2-dimensional fit to the event spatial distribution to predict the optical depth and event rate across the Galactic bulge. The most serious difficulty for the model is that it provides only $\sim 50\%$ of the measured optical depth and event rate per star at low Galactic latitude around the inner bulge ($|b|<3{^\circ}$). This discrepancy most likely is associated with known under-estimated extinction and star counts in the innermost regions and therefore provides additional support for a missing inner stellar population.

In this paper, the detectability of habitable exomoons orbiting around giant planets in M-dwarf systems using Transit Timing Variations (TTVs) and Transit Timing Durations (TDVs) with Kepler-class photometry is investigated. Light curves of systems with various configurations were simulated around M-dwarf hosts of mass 0.5 Msun and radius 0.55 Rsun. Jupiter-like giant planets which offer the best potential for hosting habitable exomoons were considered with rocky super-Earth-mass moons. The detectability is measured by using the phase-correlation between TTV and TDV signals. Since the TDV signal is typically weaker than the TTV signal, confirmation of an exomoon detection will depend on being able to detect a TDV signal. We find that exomoons around planets orbiting within the habitable zone of an M-dwarf host star can produce both detectable TTV and TDV signatures with Kepler-class photometry. While aliasing between the planet period and moon period may hinder exomoon detection, we also find some strong correlation signatures in our simulation (eg. correlation: >0.7) which would provide convincing exomoon signatures. With the addition of red noise stellar variability, correlations generally weaken. However simulated examples with planet masses less than around 25 Mearth, moons of mass 8-10 Mearth and specific values of planet and moon periods still yield detectable correlation in 25-50% of cases. Our simulation indicates that Kepler provides one of the best available opportunities for exomoon detection.