Binaries harbouring millisecond pulsars enable a unique path to determine neutron star masses: radio pulsations reveal the motion of the neutron star, while that of the companion can be characterised through studies in the optical range. PSR J1012+5307 is a millisecond pulsar in a 14.5-h orbit with a helium-core white dwarf companion. In this work we present the analysis of an optical spectroscopic campaign, where the companion star absorption features reveal one of the lightest known white dwarfs. We determine a white dwarf radial velocity semi-amplitude of K_2 = 218.9 +- 2.2 km/s, which combined with that of the pulsar derived from the precise radio timing, yields a mass ratio of q=10.44+- 0.11. We also attempt to infer the white dwarf mass from observational constraints using new binary evolution models for extremely low-mass white dwarfs, but find that they cannot reproduce all observed parameters simultaneously. In particular, we cannot reconcile the radius predicted from binary evolution with the measurement from the photometric analysis (R_WD=0.047+-0.003 Rsun). Our limited understanding of extremely low-mass white dwarf evolution, which results from binary interaction, therefore comes as the main factor limiting the precision with which we can measure the mass of the white dwarf in this system. Our conservative white dwarf mass estimate of M_WD = 0.165 +- 0.015 Msun, along with the mass ratio enables us to infer a pulsar mass of M_NS = 1.72 +- 0.16 Msun. This value is clearly above the canonical 1.4 Msun, therefore adding PSR J1012+5307 to the growing list of massive millisecond pulsars.

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 an optical, X-ray, and $\gamma$-ray study of 1SXPS J042749.2-670434, an eclipsing X-ray binary which has an associated $\gamma$-ray counterpart, 4FGL J0427.8-6704. This association has led to the source being classified as a transitional millisecond pulsar (tMSP) in an accreting state. We analyse 10.5 years of Fermi LAT data, and detect a $\gamma$-ray eclipse at the same phase as optical and X-ray eclipses at the >5$\sigma$ level, a significant improvement on the 2.8$\sigma$level of the previous detection. The confirmation of this eclipse solidifies the association between the X-ray source and the $\gamma$-ray source, strengthening the tMSP classification. However, analysis of several optical data sets and an X-ray observation do not reveal a change in the source's median brightness over long timescales or a bi-modality on short timescales. Instead, the light curve is dominated by flickering which has a correlation time of 2.6 min alongside a potential quasi-periodic oscillation at $\sim$21 min. The mass of the primary and secondary star are constrained to be $M_1=1.43^{+0.33}_{-0.19}$ M$_{\odot}$ and $M_2=0.3^{+0.17}_{-0.12}$ M$_{\odot}$ through modelling of the optical light curve. While this is still consistent with a white dwarf primary, we favour the transitional millisecond pulsar in a low accretion state classification due to the significance of the $\gamma$-ray eclipse detection.

We present the discovery of an optical accretion disk wind in the X-ray transient Swift J1858.6-0814. Our 90-spectrum data set, taken with the 10.4m GTC telescope over 8 different epochs and across five months, reveals the presence of conspicuous P-Cyg profiles in He I at 5876 Angs and Halpha. These features are detected throughout the entire campaign, albeit their intensity and main observational properties are observed to vary on time-scales as short as five minutes. In particular, we observe significant variations in the wind velocity, between a few hundreds and ~ 2400 km/s. In agreement with previous reports, our observations are characterised by the presence of frequent flares, although the relation between the continuum flux variability and the presence/absence of wind features is not evident. The reported high activity of the system at radio waves indicates that the optical wind of Swift J1858.6-0814 is contemporaneous with the radio-jet, as is the case for the handful of X-ray binary transients that have shown so far optical P-Cyg profiles. Finally, we compare our results with those of other sources showing optical accretion disk winds, with emphasis on V404 Cyg and V4641 Sgr, since they also display strong and variable optical wind features as well as similar flaring behaviour.

The Low-Frequency Array radio telescope discovered the $707$ Hz binary millisecond pulsar (MSP) J0952$-$0607 in a targeted radio pulsation search of an unidentified $\textit{Fermi}$ gamma-ray source. This source shows a weak energy flux of $F_\gamma = 2.6 \times 10^{-12}\,\text{erg}\,\text{cm}^{-2}\,\text{s}^{-1}$ in the energy range between $100\,\text{MeV}$ and $100\,\text{GeV}$. Here we report the detection of pulsed gamma-ray emission from PSR$\,$J0952$-$0607 in a very sensitive gamma-ray pulsation search. The pulsar's rotational, binary, and astrometric properties are measured over seven years of $\textit{Fermi}$-Large Area Telescope data. For this we take into account the uncertainty on the shape of the gamma-ray pulse profile. We present an updated radio-timing solution now spanning more than two years and show results from optical modeling of the black-widow-type companion based on new multi-band photometric data taken with HiPERCAM on the Gran Telescopio Canarias on La Palma and ULTRACAM on the New Technology Telescope at ESO La Silla. PSR$\,$J0952$-$0607 is now the fastest-spinning pulsar for which the intrinsic spin-down rate has been reliably constrained ($\dot{P}_\text{int} \lesssim 4.6 \times 10^{-21}\,\text{s}\,\text{s}^{-1}$). The inferred surface magnetic field strength of $B_\text{surf} \lesssim 8.2 \times 10^{7}\,\text{G}$ is among the ten lowest of all known pulsars. This discovery is another example of an extremely fast spinning black-widow pulsar hiding within an unidentified $\textit{Fermi} gamma-ray source. In the future such systems might help to pin down the maximum spin frequency and the minimum surface magnetic field strength of MSPs.

We have imaged the Sunyaev-Zel'dovich (SZ) effect signals at 140 and 270 GHz towards ten galaxy clusters with Bolocam and AzTEC/ASTE. We also used Planck data to constrain the signal at large angular scales, Herschel-SPIRE images to subtract the brightest galaxies that comprise the cosmic infrared background (CIB), Chandra imaging to map the electron temperature $T_e$ of the intra-cluster medium (ICM), and HST imaging to derive models of each galaxy cluster's mass density. The galaxy clusters gravitationally lens the background CIB, which produced an on-average reduction in brightness towards the galaxy clusters' centers after the brightest galaxies were subtracted. We corrected for this deficit, which was between 5-25% of the 270 GHz SZ effect signal within $R_{2500}$. Using the SZ effect measurements, along with the X-ray constraint on $T_e$, we measured each galaxy cluster's average line of sight (LOS) velocity $v_z$ within $R_{2500}$, with a median per-cluster uncertainty of +-700 km/s. We found an ensemble-mean <$v_z$> of 430+-210 km/s, and an intrinsic cluster-to-cluster scatter $\sigma_{int}$ of 470+-340 km/s. We also obtained maps of $v_z$ over each galaxy cluster's face with an angular resolution of 70". All four galaxy clusters previously identified as having a merger oriented along the LOS showed an excess variance in these maps at a significance of 2-4$\sigma$, indicating an internal $v_z$ rms of $\gtrsim$1000 km/s. None of the six galaxy clusters previously identified as relaxed or plane of sky mergers showed any such excess variance.

We report on a detailed optical spectroscopic follow-up of the black hole transient MAXI J1820+070 (ASASSN-18ey). The observations cover the main part of the X-ray binary outburst, when the source alternated between hard and soft states following the classical pattern widely seen in other systems. We focus the analysis on the He I emission lines at 5876 and 6678 Angs, as well as on Halpha. We detect clear accretion disk wind features (P-Cyg profiles and broad emission line wings) in the hard state, both during outburst rise and decay. These are not witnessed during the several months long soft state. However, our data suggest that the visibility of the outflow might be significantly affected by the ionisation state of the accretion disk. The terminal velocity of the wind is above ~ 1200 km/s, which is similar to outflow velocities derived from (hard-state) optical winds and (soft-state) X-ray winds in other systems. The wind signatures, in particular the P-Cyg profiles, are very shallow, and their detection has only been possible thanks to a combination of source brightness and intense monitoring at very high signal-to-noise. This study indicates that cold, optical winds are most likely a common feature of black hole accretion, and therefore, that wind-like outflows are a general mechanism of mass and angular momentum removal operating throughout the entire X-ray binary outburst.

MAXI J1807+132 is an X-ray transient discovered during the decay of an outburst in 2017. We present optical and X-ray monitoring of the source over more than 125 days, from outburst to quiescence. The outburst decay is characterized by the presence of several re-flares with a quasi-periodic recurrence time of $\sim 6.5$ days. We detect broad H and He emission lines during outburst, characteristic of transient low mass X-ray binaries. These emission lines show strong variability from epoch to epoch and, in particular, during the early stages are found embedded into deep and very broad absorption features. The quiescent spectrum shows H$\alpha$ in emission and no obvious signatures of the donor star. XMM-Newton and Swift spectra can be fitted with standard X-ray models for accreting black-holes and neutron stars, although the obtained spectral parameters favour the latter scenario. Conversely, other observables such as the optical/X-ray flux ratio, the likely systemic velocity ($\gamma \sim -150$ km s$^{-1}$) and the re-flares recurrence time suggest a black hole nature. We discuss all the above possibilities with emphasis on the strong similarities of MAXI J1807+132 with short orbital period systems.

The black hole transient V404 Cygni exhibited a bright outburst in June 2015 that was intensively followed over a wide range of wavelengths. Our team obtained high time resolution optical spectroscopy (~90 s), which included a detailed coverage of the most active phase of the event. We present a database consisting of 651 optical spectra obtained during this event, that we combine with 58 spectra gathered during the fainter December 2015 sequel outburst, as well as with 57 spectra from the 1989 event. We previously reported the discovery of wind-related features (P-Cygni and broad-wing line profiles) during both 2015 outbursts. Here, we build diagnostic diagrams that enable us to study the evolution of typical emission line parameters, such as line fluxes and equivalent widths, and develop a technique to systematically detect outflow signatures. We find that these are present throughout the outburst, even at very low optical fluxes, and that both types of outflow features are observed simultaneously in some spectra, confirming the idea of a common origin. We also show that the nebular phases depict loop patterns in many diagnostic diagrams, while P-Cygni profiles are highly variable on time-scales of minutes. The comparison between the three outbursts reveals that the spectra obtained during June and December 2015 share many similarities, while those from 1989 exhibit narrower emission lines and lower wind terminal velocities. The diagnostic diagrams presented in this work have been produced using standard measurement techniques and thus may be applied to other active low-mass X-ray binaries.

We present 5 years of optical and infrared data of the black hole candidate MAXI J1659-152 covering its 2010 outburst, decay and quiescence. Combining optical data taken during the outburst decay, we obtain an orbital period of 2.414 $\pm$ 0.005 h, in perfect agreement with the value previously measured from X-ray dips. In addition, we detect a clear H$\alpha$ excess in MAXI J1659-152 with data taken during the outburst decay. We also detect a single hump modulation most likely produced by irradiation. Assuming that the maximum occurs at orbital phase 0.5, we constrain the phase of the X-ray dips to be ~ 0.65. We also detect the quiescent optical counterpart at r' = 24.20 $\pm$ 0.08, I = 23.32 $\pm$ 0.02 and H = 20.7 $\pm$ 0.1. These magnitudes provide colour indices implying an M2-M5 donor star assuming 60% contribution from a disc component in the r'-band.

We present a detailed spectroscopic study of the optical counterpart of the neutron star X-ray transient Aquila X-1 during its 2011, 2013 and 2016 outbursts. We use 65 intermediate resolution GTC-10.4m spectra with the aim of detecting irradiation-induced Bowen blend emission from the donor star. While Gaussian fitting does not yield conclusive results, our full phase coverage allows us to exploit Doppler mapping techniques to independently constrain the donor star radial velocity. By using the component N III 4640.64/4641.84 A we measure Kem = 102 +-6 km/s. This highly significant detection ( >13 sigma) is fully compatible with the true companion star radial velocity obtained from near-infrared spectroscopy during quiescence. Combining these two velocities we determine, for the first time, the accretion disc opening angle and its associated error from direct spectroscopic measurements and detailed modelling, obtaining alpha = 15.5 +2.5 -5 deg. This value is consistent with theoretical work if significant X-ray irradiation is taken into account and is important in the light of recent observations of GX339-4, where discrepant results were obtained between the donor's intrinsic radial velocity and the Bowen-inferred value. We also discuss the limitations of the Bowen technique when complete phase coverage is not available.

The new recurrent X-ray transient MAXI J1957+032 has had four X-ray outbursts within 16 months, all very briefly detected (lasting <5 days). During the most recent event (Sep/Oct 2016), we obtained with SALT the first optical spectrum of the transient counterpart, showing the classic blue continuum of an X-ray irradiated disc in an LMXB and no other features. At high Galactic latitude below the plane (-13 deg) reddening is low but there is no quiescent counterpart visible on any of the existing sky surveys, nor any other known X-ray source in the region. Swift monitoring of 3 of the 4 events is presented, showing rapidly fading X-ray outbursts together with significant UVOT detections in the UV (W1,M2,W2), U and B bands. The optical properties are most like those of the short-period LMXBs, which, combined with the softening witnessed during the decay to quiescence would place the system at d <13 kpc. The short duration and short recurrence time of the outbursts are reminiscent of the AMXPs, which exhibit peak luminosities of 1 % L_Edd. Assuming this peak luminosity would place MAXI J1957+032 at a distance of d= 5-6 kpc.

After more than 26 years in quiescence, the black hole transient V404 Cyg went into a luminous outburst in June 2015, and additional activity was detected in late December of the same year. Here, we present an optical spectroscopic follow-up of the December mini-outburst, together with X-ray, optical and radio monitoring that spanned more than a month. Strong flares with gradually increasing intensity are detected in the three spectral ranges during the ~10 days following the Swift trigger. Our optical spectra reveal the presence of a fast outflowing wind, as implied by the detection of a P-Cyg profile (He I - 5876 A) with a terminal velocity of ~2500 km/s. Nebular-like spectra -- with an H_alpha equivalent width of ~500 A -- are also observed. All these features are similar to those seen during the main June 2015 outburst. Thus, the fast optical wind simultaneous with the radio jet is most likely present in every V404 Cyg outburst. Finally, we report on the detection of a strong radio flare in late January 2016, when X-ray and optical monitoring had stopped due to Sun constraints.

The low mass X-ray binary Aquila X-1 is one of the most active neutron star X-ray transients. Despite it has a relatively bright quiescent optical counterpart, the detection of its companion has been hampered by the presence of a nearby interloper star. Using the infrared integral field spectrograph SINFONI on the VLT-8.2m telescope, we unambiguously single out Aquila X-1 from the interloper. Phase-resolved near infrared spectroscopy reveals absorption features from a K4 +- 2 companion star moving at a projected velocity of K_2= 136 +- 4 km/s. We here present the first dynamical solution and associated fundamental parameters of Aquila X-1, imposing new constraints to the orbital inclination (36 deg < i < 47 deg) and the distance (d = 6 +- 2 kpc) to this prototypical neutron star transient.

Accretion of matter onto black holes is universally associated with strong radiative feedback and powerful outflows. In particular, black hole transients show outflows whose properties are strongly coupled to those of the accretion flow. This includes X-ray winds of ionized material, expelled from the accretion disc encircling the black hole, and collimated radio jets. Very recently, a distinct optical variability pattern has been reported in the transient black hole transient V404 Cyg, and interpreted as disrupted mass flow into the inner regions of its large accretion disc. Here, we report on the discovery of a sustained outer accretion disc wind in V404 Cyg, which is unlike any seen previously. We find that the outflowing wind is neutral, has a large covering factor, expands at 1% of the speed of light and triggers a nebular phase once accretion sharply drops and the ejecta become optically thin. The large expelled mass (> 10^-8 Msun) indicates that the outburst was prematurely ended when a sizeable fraction of the outer disc was depleted by the wind, detaching the inner regions from the rest of the disc. The luminous, but brief, accretion phases shown by transients with large accretion discs imply that this outflow is most likely a new fundamental ingredient regulating mass accretion onto black holes.

The radio emission from Sgr A$^\ast$ is thought to be powered by accretion onto a supermassive black hole of $\sim\! 4\times10^6~ \rm{M}_\odot$ at the Galactic Center. At millimeter wavelengths, Very Long Baseline Interferometry (VLBI) observations can directly resolve the bright innermost accretion region of Sgr A$^\ast$. Motivated by the addition of many sensitive, long baselines in the north-south direction, we developed a full VLBI capability at the Large Millimeter Telescope Alfonso Serrano (LMT). We successfully detected Sgr A$^\ast$ at 3.5~mm with an array consisting of 6 Very Long Baseline Array telescopes and the LMT. We model the source as an elliptical Gaussian brightness distribution and estimate the scattered size and orientation of the source from closure amplitude and self-calibration analysis, obtaining consistent results between methods and epochs. We then use the known scattering kernel to determine the intrinsic two dimensional source size at 3.5 mm: $(147\pm7~\mu\rm{as}) \times (120\pm12~\mu\rm{as})$, at position angle $88^\circ\pm7^\circ$ east of north. Finally, we detect non-zero closure phases on some baseline triangles, but we show that these are consistent with being introduced by refractive scattering in the interstellar medium and do not require intrinsic source asymmetry to explain.

Swift J1357.2-0933 is one of the shortest orbital period black hole X-ray transients (BHTs). It exhibited deep optical dips together with an extremely broad H$\alpha$ line during outburst. We present 10.4-m GTC time-resolved spectroscopy during quiescence searching for donor star absorption features. The large contribution of the accretion flow to the total luminosity prevents the direct detection of the companion. Nevertheless, we constrain the non-stellar contribution to be larger than $\sim 80\%$ of the total optical light, which sets new lower limits to the distance ($d > 2.29\, \rm{kpc}$) and the height over the Galactic plane ($z>1.75\, \rm{kpc}$). This places the system in the galactic thick disc. We measure a modulation in the centroid of the H$\alpha$ line with a period of $P=0.11\pm0.04\, \rm{d}$ which, combined with the recently presented FWHM-$K_2$ correlation, results in a massive black hole ($M_1>9.3 \, \rm{M_\odot}$) and a $\sim$ M2V companion star ($M_2\sim 0.4\, \rm{M_\odot}$). We also present further evidence supporting a very high orbital inclination ($i\gtrsim 80^\circ$).

More than 50 years after the dawn of X-ray astronomy, the dynamical parameters of the prototypical X-ray binary Sco X-1 are still unknown. We combine a Monte Carlo analysis, which includes all the previously known orbital parameters of the system, along with the K-correction to set dynamical constraints to the masses of the compact object (M1<1.73 Msun) and the companion star (0.28 Msun<M2<0.70 Msun). For the case of a canonical neutron star mass of M1= 1.4 Msun, the orbital inclination is found to be lower than 40 degree. We also present the best near-infrared spectrum of the source to date. There is no evidence of donor star features on it, but we are able to constrain the veiling factor as a function of the spectral type of the secondary star. The combination of both techniques restricts the spectral type of the donor to be later than K4 and luminosity class IV. It also constrains the contribution of the companion light to the infrared emission of Sco X-1 to be lower than 33 percent. This implies that the accretion related luminosity of the system in the K band is larger than 4 x 10E35 erg/s.

It is well-known that stars with giant planets are on average more metal-rich than stars without giant planets, whereas stars with detected low-mass planets do not need to be metal-rich. With the aim of studying the weak boundary that separates giant planets and brown dwarfs (BDs) and their formation mechanism, we analyze the spectra of a sample of stars with already confirmed BD companions both by radial velocity and astrometry. We employ standard and automatic tools to perform an EW-based analysis and to derive chemical abundances from CORALIE spectra of stars with BD companions. We compare these abundances with those of stars without detected planets and with low-mass and giant-mass planets. We find that stars with BDs do not have metallicities and chemical abundances similar to those of giant-planet hosts but they resemble the composition of stars with low-mass planets. The distribution of mean abundances of $\alpha$-elements and iron peak elements of stars with BDs exhibit a peak at about solar abundance whereas for stars with low-mass and high-mass planets the [X$_\alpha$/H] and [X$_{\rm Fe}$/H] peak abundances remain at $\sim -0.1$~dex and $\sim +0.15$~dex, respectively. We display these element abundances for stars with low-mass and high-mass planets, and BDs versus the minimum mass, $m_C \sin i$, of the most-massive substellar companion in each system, and we find a maximum in $\alpha$-element as well as Fe-peak abundances at $m_C \sin i \sim 1.35\pm 0.20$ jupiter masses. We discuss the implication of these results in the context of the formation scenario of BDs in comparison with that of giant planets.