The Robotic Observations of Microlensing Events/Reactive Event Assessment Survey was a Key Project at Las Cumbres Observatory (hereafter LCO) which continuously monitored 20 selected fields (3.76 sq.deg) in the Galactic Bulge throughout their seasonal visibility window over a three-year period, between 2017 March and 2020 March. Observations were made in three optical passbands (SDSS-g ', -r ', -i '), and LCO's multi-site telescope network enabled the survey to achieve a typical cadence of similar to 10 hr in i ' and similar to 15 hr in g ' and r ' . In addition, intervals of higher cadence (<1 hr) data were obtained during monitoring of key microlensing events within the fields. This paper describes the Difference Image Analysis data reduction pipeline developed to process these data, and the process for combining the photometry from LCO's three observing sites in the Southern Hemisphere. The full timeseries photometry for all similar to 8 million stars, down to a limiting magnitude of i similar to 18 mag is provided in the data release accompanying this paper, and samples of the data are presented for exemplar microlensing events, illustrating how the tri-band data are used to derive constraints on the microlensing source star parameters, a necessary step in determining the physical properties of the lensing object. The timeseries data also enables a wealth of additional science, for example in characterizing long-timescale stellar variability, and a few examples of the data for known variables are presented.

Aims. We present spectroscopy of the accreting X-ray binary and millisecond pulsar SAX J1808.4-3658. These observations are the first to be obtained during a reflaring phase. We collected spectroscopic data during the beginning of the reflaring of the 2019 outburst and compared them to previous datasets taken at different epochs, both of the same outburst and across the years. To this end, we also present spectra of the source taken during quiescence in 2007, one year before the next outburst. Methods. We made use of data taken by the Very Large Telescope (VLT) X-shooter spectrograph on August 31, 2019, three weeks after the outburst peak. For flux calibration, we used photometric data taken during the same night by the 1m telescopes from the Las Cumbres Observatory network that are located in Chile. We compare our spectra to the quiescent data taken by the VLT-FORS1 spectrograph in September, 2007. We inspected the spectral energy distribution by fitting our data with a multicolored accretion-disk model and sampled the posterior probability density function for the model parameters with a Markov chain Monte Carlo (MCMC) algorithm. Results. We find the optical spectra of the 2019 outburst to be unusually featureless, with no emission lines present despite the high resolution of the instrument. Fitting the UV-optical spectral energy distribution with a disk plus irradiated star model results in a very large value for the inner disk radius of similar to 5130 +/- 240 km, which could suggest that the disk was emptied of material during the outburst, possibly accounting for the emission-less spectra. Alternatively, the absence of emission lines could be due to a significant contribution of the jet emission at optical wavelengths.

We present long-term (2012-2022) optical monitoring of the candidate black hole X-ray binary Swift J1910.2-0546 with the Faulkes Telescopes and Las Cumbres Observatory (LCO) network. Following its initial bright 2012 outburst, we find that the source displayed a series of at least 7 quasi-periodic, high amplitude (~3 mags) optical reflares in 2013, with a recurrence time increasing from ~42 days to ~49 days. In 2014, the source experienced a mini-outburst with two peaks in the optical. We also study the recent 2022 outburst of the source at optical wavelengths, and perform a comparative analysis with the earlier rebrightenings. A single X-ray detection and only two radio detections were obtained during the 2013 reflaring period, and only optical detections were acquired in 2014. During the reflaring in both 2013 and 2014, the source showed bluer-when-brighter behavior, having optical colors consistent with a blackbody heating and cooling between 4500 and 9500 K, i.e. the temperature range in which hydrogen starts to ionize. Finally, we compare the flaring behavior of the source to re-brightening events in other X-ray binaries. We show that the repeated reflarings of Swift J1910.2-0546 are highly unusual, and propose that they arise from a sequence of repetitive heating and cooling front reflections travelling through the accretion disk.

Transitional millisecond pulsars are an emerging class of sources that link low-mass X-ray binaries to millisecond radio pulsars in binary systems. These pulsars alternate between a radio pulsar state and an active low-luminosity X-ray disc state. During the active state, these sources exhibit two distinct emission modes (high and low) that alternate unpredictably, abruptly, and incessantly. X-ray to optical pulsations are observed only during the high mode. The root cause of this puzzling behaviour remains elusive. This paper presents the results of the most extensive multi-wavelength campaign ever conducted on the transitional pulsar prototype, PSR J1023+0038, covering from the radio to X-rays. The campaign was carried out over two nights in June 2021 and involved 12 different telescopes and instruments, including XMM-Newton, HST, VLT/FORS2 (in polarimetric mode), ALMA, VLA, and FAST. By modelling the broadband spectral energy distributions in both emission modes, we show that the mode switches are caused by changes in the innermost region of the accretion disc. These changes trigger the emission of discrete mass ejections, which occur on top of a compact jet, as testified by the detection of at least one short-duration millimetre flare with ALMA at the high-to-low mode switch. The pulsar is subsequently re-enshrouded, completing our picture of the mode switches.

We present a detailed study of optical data from the 2012 outburst of the candidate black hole X-ray binary Swift J1910.2-0546 using the Faulkes Telescope and Las Cumbres Observatory (LCO). We analyse the peculiar spectral state changes of Swift J1910.2-0546 in different energy bands, and characterise how the optical and UV emission correlates with the unusual spectral state evolution. Using various diagnostic tools like the optical/X-ray correlation and spectral energy distributions, we disentangle the different emission processes contributing towards the optical flux of the system. When Swift J1910.2-0546 transitions to the pure hard state, we find significant optical brightening of the source along with a dramatic change in the optical colour due to the onset of a jet during the spectral state transition. For the rest of the spectral states, the optical/UV emission is mostly dominated by an X-ray irradiated disk. From our high cadence optical study, we have discovered a putative modulation. Assuming that this modulation arises from a superhump, we suggest Swift J1910.2-0546 to have an orbital period of 2.25-2.47 hr, which would make it the shortest orbital period black hole X-ray binary known to date. Finally, from the state transition luminosity of the source, we find that the distance to the source is likely to be ~4.5-20.8 kpc, which is also supported by the comparative position of the source in the global optical/X-ray correlation of a large sample of black hole and neutron star X-ray binaries.

Understanding the inter-relationships between traffic flow, density, and speed through the study of the fundamental diagram of road traffic is critical for traffic modelling and management. Consequently, over the last 85 years, a wealth of models have been developed for its functional form. However, there has been no clear answer as to which model is the most appropriate for observed (i.e. empirical) fundamental diagrams and under which conditions. A lack of data has been partly to blame. Motivated by shortcomings in previous reviews, we first present a comprehensive literature review on modelling the functional form of empirical fundamental diagrams. We then perform fits of 50 previously proposed models to a high quality sample of 10,150 empirical fundamental diagrams pertaining to 25 cities. Comparing the fits using information criteria, we find that the non-parametric Sun model greatly outperforms all of the other models. The Sun model maintains its winning position regardless of road type and congestion level. Our study, the first of its kind when considering the number of models tested and the amount of data used, finally provides a definitive answer to the question ``Which model for the functional form of an empirical fundamental diagram is currently the best?''. The word ``currently'' in this question is key, because previously proposed models adopt an inappropriate Gaussian noise model with constant variance. We advocate that future research should shift focus to exploring more sophisticated noise models. This will lead to an improved understanding of empirical fundamental diagrams and their underlying functional forms.

We report on the observations, analysis and interpretation of the microlensing event MOA-2019- BLG-008. The observed anomaly in the photometric light curve is best described through a binary lens model. In this model, the source did not cross caustics and no finite source effects were observed. Therefore the angular Einstein ring radius cannot be measured from the light curve alone. However, the large event duration, t E about 80 days, allows a precise measurement of the microlensing parallax. In addition to the constraints on the angular radius and the apparent brightness I s of the source, we employ the Besancon and GalMod galactic models to estimate the physical properties of the lens. We find excellent agreement between the predictions of the two Galactic models: the companion is likely a resident of the brown dwarf desert with a mass Mp about 30 MJup and the host is a main sequence dwarf star. The lens lies along the line of sight to the Galactic Bulge, at a distance of less then4 kpc. We estimate that in about 10 years, the lens and source will be separated by 55 mas, and it will be possible to confirm the exact nature of the lensing system by using high-resolution imaging from ground or space-based observatories.

ABSTRACT In traditional lucky imaging (TLI), many consecutive images of the same scene are taken with a high frame-rate camera, and all but the sharpest images are discarded before constructing the final shift-and-add image. Here, we present an alternative image analysis pipeline – The Thresher – for these kinds of data, based on online multi-frame blind deconvolution. It makes use of all available data to obtain the best estimate of the astronomical scene in the context of reasonable computational limits; it does not require prior estimates of the point-spread functions in the images, or knowledge of point sources in the scene that could provide such estimates. Most importantly, the scene it aims to return is the optimum of a justified scalar objective based on the likelihood function. Because it uses the full set of images in the stack, The Thresher outperforms TLI in signal-to-noise ratio; as it accounts for the individual-frame PSFs, it does this without loss of angular resolution. We demonstrate the effectiveness of our algorithm on both simulated data and real Electron-Multiplying CCD images obtained at the Danish 1.54-m telescope (hosted by ESO, La Silla). We also explore the current limitations of the algorithm, and find that for the choice of image model presented here, non-linearities in flux are introduced into the returned scene. Ongoing development of the software can be viewed at https://github.com/jah1994/TheThresher.

We present a detailed study of the evolution of the Galactic black hole transient GRS 1716-249 during its 2016-2017 outburst at optical (Las Cumbres Observatory), mid-infrared (Very Large Telescope), near-infrared (Rapid Eye Mount telescope), and ultraviolet (the Neil Gehrels Swift Observatory Ultraviolet/Optical Telescope) wavelengths, along with archival radio and X-ray data. We show that the optical/near-infrared and UV emission of the source mainly originates from a multi-temperature accretion disk, while the mid-infrared and radio emission are dominated by synchrotron emission from a compact jet. The optical/UV flux density is correlated with the X-ray emission when the source is in the hard state, consistent with an X-ray irradiated accretion disk with an additional contribution from the viscous disk during the outburst fade. We also report the long-term optical light curve of the source and find that the quiescent i-band magnitude is 21.39$\pm$0.15 mag. Furthermore, we discuss how previous estimates of the system parameters of the source are based on various incorrect assumptions, and so are likely to be inaccurate. By comparing our GRS 1716-249 dataset to those of other outbursting black hole X-ray binaries, we find that while GRS 1716-249 shows similar X-ray behaviour, it is noticeably optically fainter, if the literature distance of 2.4 kpc is adopted. Using several lines of reasoning, we argue that the source distance is further than previously assumed in the literature, likely within 4-17 kpc, with a most likely range of $\sim$4-8 kpc.

We report the first unambiguous detection and mass measurement of an isolated stellar-mass black hole (BH). We used the Hubble Space Telescope (HST) to carry out precise astrometry of the source star of the long-duration ( t E ≃ 270 days), high-magnification microlensing event MOA-2011-BLG-191/OGLE-2011-BLG-0462 (hereafter designated as MOA-11-191/OGLE-11-462), in the direction of the Galactic bulge. HST imaging, conducted at eight epochs over an interval of 6 yr, reveals a clear relativistic astrometric deflection of the background star’s apparent position. Ground-based photometry of MOA-11-191/OGLE-11-462 shows a parallactic signature of the effect of Earth’s motion on the microlensing light curve. Combining the HST astrometry with the ground-based light curve and the derived parallax, we obtain a lens mass of 7.1 ± 1.3 M ⊙ and a distance of 1.58 ± 0.18 kpc. We show that the lens emits no detectable light, which, along with having a mass higher than is possible for a white dwarf or neutron star, confirms its BH nature. Our analysis also provides an absolute proper motion for the BH. The proper motion is offset from the mean motion of Galactic disk stars at similar distances by an amount corresponding to a transverse space velocity of ∼45 km s −1 , suggesting that the BH received a “natal kick” from its supernova explosion. Previous mass determinations for stellar-mass BHs have come from radial velocity measurements of Galactic X-ray binaries and from gravitational radiation emitted by merging BHs in binary systems in external galaxies. Our mass measurement is the first for an isolated stellar-mass BH using any technique.

We report on a long-term optical monitoring of the neutron star X-ray binary Centaurus X-4 performed during the last 13.5 years. This source has been in quiescence since its outburst in 1979. Our monitoring reveals the overall evolution of the accretion disc; we detect short-duration flares, likely originating also in the disc, superimposed with a small-amplitude (< 0.1 mag) ellipsoidal modulation from the companion star due to geometrical effects. A long-term (~2300 days) downward trend, followed by a shorter (~1000 days) upward one, is observed in the disc light curve. Such a rise in the optical has been observed for other X-ray binaries preceding outbursts, as predicted by the disc instability model. For Cen X-4, the rise of the optical flux proceeded for ~3 years, and culminated in a flux increase at all wavelengths (optical-UV-X-rays) at the end of 2020. This increase faded after ~2 weeks, without giving rise to a full outburst. We suggest that the propagation of an inside-out heating front was ignited due to a partial ionization of hydrogen in the inner disc. The propagation might have stalled soon after the ignition due to the increasing surface density in the disc that the front encountered while propagating outwards. The stall was likely eased by the low level irradiation of the outer regions of the large accretion disc, as shown by the slope of the optical/X-ray correlation, suggesting that irradiation does not play a strong role in the optical, compared to other sources of emission.

We report an analysis of the planetary microlensing event OGLE-2018-BLG-1185, which was observed by a large number of ground-based telescopes and by the Spitzer Space Telescope. The ground-based light curve indicates a low planet–host star mass ratio of q = (6.9 ± 0.2) × 10−5, which is near the peak of the wide-orbit exoplanet mass-ratio distribution. We estimate the host star and planet masses with a Bayesian analysis using the measured angular Einstein radius under the assumption that stars of all masses have an equal probability of hosting the planet. The flux variation observed by Spitzer is marginal, but still places a constraint on the microlens parallax. Imposing a conservative constraint that this flux variation should be Δf Spz < 4 instrumental flux units yields a host mass of Mhost=0.37−0.21+0.35M⊙ and a planet mass of mp=8.4−4.7+7.9M⊕ . A Bayesian analysis including the full parallax constraint from Spitzer suggests smaller host star and planet masses of Mhost=0.091−0.018+0.064M⊙ and mp=2.1−0.4+1.5M⊕ , respectively. Future high-resolution imaging observations with the Hubble Space Telescope or Extremely Large Telescope could distinguish between these two scenarios and help reveal the planetary system properties in more detail.

Millisecond spinning, low-magnetic-field neutron stars are believed to attain their fast rotation in a 0.1–1-Gyr-long phase during which they accrete matter endowed with angular momentum from a low-mass companion star1. Despite extensive searches, coherent periodicities originating from accreting neutron star magnetospheres have been detected only at X-ray energies2 and in ~10% of the currently known systems3. Here we report the detection of optical and ultraviolet coherent pulsations at the X-ray period of the transient low-mass X-ray binary system SAX J1808.4−3658, during an accretion outburst that occurred in August 20194. At the time of the observations, the pulsar was surrounded by an accretion disk, displayed X-ray pulsations and its luminosity was consistent with magnetically funnelled accretion onto the neutron star. Current accretion models fail to account for the luminosity of both optical and ultraviolet pulsations; these are instead more likely to be driven by synchro-curvature radiation5,6 in the pulsar magnetosphere or just outside of it. This interpretation would imply that particle acceleration can take place even when mass accretion is going on, and opens up new perspectives in the study of coherent optical/ultraviolet pulsations from fast-spinning accreting neutron stars in low-mass X-ray binary systems. Outbursts from accreting pulsars encode much information on mass accretion in X-ray binary systems. Measuring optical as well as X-ray pulsations can constrain models and, indeed, point to particle acceleration taking place during accretion.

ABSTRACT IGR J17285−2922 is a known X-ray binary with a low peak 2–10 keV X-ray luminosity of ∼ 1036 erg s−1 during outburst. IGR J17285−2922 exhibited two outbursts in 2003 and 2010 and went into outburst again in 2019. We have monitored this ∼ 4-month long 2019 outburst with Swift in X-ray and the Very Large Array in radio. We have also obtained four optical spectra with the Gran Telescopio Canarias and Southern Astrophysical Research Telescope, three optical photometry measurements with the Las Cumbres Observatory, and one near-infrared spectrum with the Gemini South telescope. The ratio between its X-ray and radio luminosity is consistent with both samples of neutron star and black hole (BH) X-ray binaries, while the ratio between the X-ray and optical luminosity is consistent with BH X-ray binaries. Studying the evolution of its X-ray power-law index throughout the outburst, we find additional evidence for a BH as compact object. The four optical spectra show no H α emission and the nIR spectrum shows no Brγ emission, suggesting that the donor star could be hydrogen-poor and hence that IGR J17285−2922 might have an ultracompact binary orbit. The shape of the X-ray light curve is well described by an exponential, followed by a linear decay, from which we obtain a relation between the orbital period Porb and the binary mass ratio. We discuss how this relation is consistent with theoretical predictions and known ultracompact X-ray binaries. Lastly, we discuss how the observed properties are reminiscent of short-Porb BH X-ray binaries.