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.

Context:Spider pulsars are binary systems containing an energetic millisecond pulsar that intensely irradiates a closely orbiting low-mass companion. Modelling their companion's optical light curves is essential to the study of the orbital properties of the binary, including the determination of the pulsar mass, characterising the pulsar wind and the star itself. We aim to generalise the traditional direct heating model of irradiation, whereby energy deposited by the pulsar wind into the stellar envelope is locally re-emitted, by introducing heat redistribution via diffusion and convection within the outer stellar envelope. We approximate the irradiated stellar envelope as a two-dimensional shell. This allows us to propose an effective equation of energy conservation that can be solved at a reduced computational cost. We then implement this model in the \texttt{Icarus} software and use evidence sampling to determine the most likely convection and diffusion laws for the light curve of the redback companion of PSR J2215+5135. Redistribution effects concentrate near the terminator line of pulsar irradiation, and can create apparent hot and cold spots. Among the models tested for PSR J2215+5135, we find that all models with heat redistribution are more likely than symmetric direct heating. The best-fitting redistribution model involves diffusion together with a uniformly rotating envelope. However, we caution that all models still present serious systematic effects, and that prior knowledge from pulsar timing, spectroscopy and distance are key to determine with certainty the most accurate redistribution law. We propose an extension of the direct heating framework that allows for exploring a variety of heat redistribution effects. Future work is necessary to determine the relevant laws from first principles and empirically using complementary observations.

Between May 2016 and September 2018, the intermediate polar (IP) FO Aquarii exhibited two distinct low states and one failed low state. We present optical spectroscopy of FO Aquarii throughout this period, making this the first detailed study of an accretion disc during a low state in any IP. Analysis of these data confirm that the low states are the result of a drop in the mass transfer rate between the secondary star and the magnetic white dwarf primary, and are characterised by a decrease in the system's brightness coupled with a change of the system's accretion structures from an accretion disc-fed geometry to a combination of disc-fed and ballistic stream-fed accretion, and that effects from accretion onto both magnetic poles become detectable. The failed low state only displays a decrease in brightness, with the accretion geometry remaining primarily disc-fed. We also find that the WD appears to be exclusively accretion disc-fed during the high state. There is evidence for an outflow close to the impact region between the ballistic stream and the disc which is detectable in all of the states. Finally, there is marginal evidence for narrow high velocity features in the H$\alpha$ emission line during the low states which may arise due to an outflow from the WD. These features may be evidence of a collimated jet, a long predicted yet elusive feature of cataclysmic variables.

The intermediate polar FO Aquarii (FO Aqr) experienced its first-reported low-accretion states in 2016, 2017, and 2018, and using newly available photographic plates, we identify pre-discovery low states in 1965, 1966, and 1974. The primary focus of our analysis, however, is an extensive set of time-series photometry obtained between 2002 and 2018, with particularly intensive coverage of the 2016-2018 low states. After computing an updated spin ephemeris for the white dwarf (WD), we show that its spin period began to increase in 2014 after having spent 27 years decreasing; no other intermediate polar has experienced a sign change of its period derivative, but FO Aqr has now done so twice. Our central finding is that the recent low states all occurred shortly after the WD began to spin down, even though no low states were reported in the preceding quarter-century, when it was spinning up. Additionally, the system's mode of accretion is extremely sensitive to the mass-transfer rate, with accretion being almost exclusively disk-fed when FO Aqr is brighter than V~14 and substantially stream-fed when it is not. Even in the low states, a grazing eclipse remains detectable, confirming the presence of a disk-like structure (but not necessarily a Keplerian accretion disk). We relate these various observations to theoretical predictions that during the low state, the system's accretion disk dissipates into a non-Keplerian ring of diamagnetic blobs. Finally, a new XMM-Newton observation from a high state in 2017 reveals an anomalously soft X-ray spectrum and diminished X-ray luminosity compared to pre-2016 observations.

The Kepler spacecraft observed a total of only four AM Herculis cataclysmic variable stars during its lifetime. We analyze the short-cadence K2 light curve of one of those systems, Tau 4 (RX J0502.8+1624), which underwent a serendipitous jump from a low-accretion state into a high state during the final days of the observation. Apart from one brief flare, there was no evidence of accretion during the 70 d of observations of the low state. As Tau 4 transitioned into a high state, the resumption of accretion was very gradual, taking approximately six days (~90 binary orbits). We supplement Tau 4's K2 light curve with time-resolved spectroscopy obtained in both high and low states of accretion. High-excitation lines, such as He II 468.6 nm, were extraordinarily weak, even when the system was actively accreting. This strongly suggests the absence of an accretion shock, placing Tau 4 in the bombardment regime predicted for AM Herculis systems with low accretion rates. In both the high-state and low-state spectra, Zeeman absorption features from the white dwarf's photosphere are present and reveal a surface-averaged field strength of $15\pm2$ MG. Remarkably, the high-state spectra also show Zeeman-split emission lines produced in a region with a field strength of $12\pm1$ MG. Zeeman emission has not been previously reported in an AM Herculis system, and we propose that the phenomenon is caused by a temperature inversion in the WD's atmosphere near the accretion region.

We present the first measurement of the gravitational quadrupole moment of the companion star of a spider pulsar, namely the black widow PSR J2051$-$0827. To this end we have re-analysed radio timing data using a new model which is able to account for periastron precession caused by tidal and centrifugal deformations of the star as well as by general relativity. The model allows for a time-varying component of the quadrupole moment, thus self-consistently accounting for the ill-understood orbital period variations observed in these systems. Our analysis results in the first detection of orbital precession in a spider system at $\dot{\omega} = -68.6_{-0.5}^{+0.9}$ deg/yr and the most accurate determination of orbital eccentricity for PSR J2051$-$0827 with $e = (4.2 \pm 0.1) \times 10^{-5}$. We show that the variable quadrupole component is about 100 times smaller than the average quadrupole moment $\bar{Q} = -2.2_{-1}^{+0.6} \times 10^{41}$ kg.m$^2$. We discuss how accurate modelling of high precision optical light curves of the companion star will allow its apsidal motion constant to be derived from our results.

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.

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 present results from a study of TESS observations of the eclipsing dwarf nova system Z Cha, covering both an outburst and a superoutburst. We discover that Z Cha undergoes hysteretic loops in eclipse depth - out-of-eclipse flux space in both the outburst and the superoutburst. The direction that these loops are executed in indicates that the disk size increases during an outburst before the mass transfer rate through the disk increases, placing constraints on the physics behind the triggering of outbursts and superoutbursts. By fitting the signature of the superhump period in a flux-phase diagram, we find the rate at which this period decreases in this system during a superoutburst for the first time. We find that the superhumps in this source skip evolutionary stage "A" seen during most dwarf nova superoutbursts, even though this evolutionary stage has been seen during previous superoutbursts of the same object. Finally, O-C values of eclipses in our sample are used to calculate new ephemerides for the system, strengthening the case for a third body in Z Cha and placing new constraints on its orbit.

The TESS spacecraft observed the asynchronous polar CD Ind at a two-minute cadence almost continuously for 28 days in 2018, covering parts of 5 consecutive cycles of the system's 7.3-day beat period. These observations provide the first uninterrupted photometry of a full spin-orbit beat cycle of an asynchronous polar. Twice per beat cycle, the accretion flow switched between magnetic poles on the white dwarf, causing the spin pulse of the white dwarf (WD) to alternate between two waveforms after each pole-switch. An analysis of the waveforms suggests that one accretion region is continuously visible when it is active, while the other region experiences lengthy self-eclipses by the white dwarf. We argue that the previously accepted periods for both the binary orbit and the WD spin have been misidentified, and while the cause of this misidentification is a subtle and easily overlooked effect, it has profound consequences for the interpretation of the system's accretion geometry and doubles the estimated time to resynchronization. Moreover, our timings of the photometric maxima do not agree with the quadratic ephemeris from Myers et al. (2017), and it is possible that the optical spin pulse might be an unreliable indicator of the white dwarf's rotation. Finally, we use Doppler tomography of archival time-resolved spectra from 2006 to study the accretion flow. While the accretion flow showed a wider azimuthal extent than is typical for synchronous polars, it was significantly less extended than in the three other asynchronous polars for which Doppler tomography has been reported.

We obtained high temporal resolution spectroscopy of the unusual binary system AR Sco covering nearly an orbit. The H$\alpha$ emission shows a complex line structure similar to that seen in some polars during quiescence. Such emission is thought to be due to long-lived prominences originating on the red dwarf. A difference between AR Sco and these other systems is that the white dwarf in AR Sco is rapidly spinning relative to the orbital period. "Slingshot" prominences stable at 3 to 5 stellar radii require surface magnetic fields between 100 and 500 G. This is comparable to the estimated WD magnetic field strength near the surface of the secondary. Our time-resolved spectra also show emission fluxes, line equivalent widths, and continuum color varying over the orbit and the beat/spin periods of the system. During much of the orbit, the optical spectral variations are consistent with synchrotron emission with the highest energy electrons cooling between pulses. On the time-scale of the beat/spin period we detect red and blue-shifted H$\alpha$ emission flashes that reach velocities of 700 km/s. Red-shifted Balmer emission flashes are correlated with the bright phases of the continuum beat pulses while blue-shifted flashes appear to prefer the time of minimum in the beat light curve. We propose that much of the energy generated in AR Sco comes from fast magnetic reconnection events occurring near the inward face of the secondary and we show that the energy generated by magnetic reconnection can account for the observed excess luminosity from the system.

We used the light curve code XRBinary to model the quiescent K2 light curves of three low-inclination cataclysmic variables (CVs): 1RXS\,J0632+2536 (J0632+2536), RZ\,Leo, TW\,Vir and the pre-CV WD\,1144+011. Optimized light curve models were obtained using a nonlinear fitting code NMfit and visualized by Phoebe 2.0. The disk model of J0632+2536 shows that one hotspot at the edge of the disk is enough to describe its light curve, while the other two dwarf nova (DN): RZ\,Leo and TW\,Vir require two hotspots. A typical pre-CV model with a weak irradiation effect for WD\,1144+011 can explain its single-hump modulation, and the newly observed spectrum confirms its previous classification. The synthetic analyses for the DN clearly indicate that phase zero of the double-hump modulations occurs around the secondary minimum and the primary hump is mainly caused by the hotspot at the edge of the disk. The quiescent disk has a flat temperature distribution with a power index of $\sim0.11$. The disk model of RZ\,Leo implies a truncated disk, supporting its previously speculated classification as an intermediate polar (IP). Except for the IP model of RZ\,Leo, which lacks a component related to the inferred accretion curtain, the models of J0632+2536, TW\,Vir and WD\,1144+011 are consistent with results from the Gaia mission. The derived masses and radii of the secondaries of the three DN are consistent with the semi-empirical relations for CV donor stars, while their effective temperatures are higher than the predictions. Irradiation of the donor stars is investigated to explain this discrepancy.

We present a long-term light curve of the precataclysmic variable (CV) V1082 Sgr obtained by the K2-mission over the course of 81 days. We analyze the entire complex light curve as well as explore several sections in detail with a sliding periodogram. The long dataset allows the first detection of the orbital period in the light curve, as well as the confirmation of cyclical variability on a longer timescale of about a month. A portion of the light curve in deep minimum reveals a clean, near-sinusoidal variability attributed to the rotation of the spotted surface of the donor star. We model that portion of the light curve assuming that the donor star grossly under-fills its Roche lobe, has cool spots similar to a chromospherically active, slightly evolved early K-star, and might be irradiated by the X-ray beam from the magnetically accreting white dwarf. The fast variability of the object in the active phases resembles the light curves of magnetic CVs (polars).

CRTS J035905.9+175034 is the first eclipsing SU UMa system for which a superoutburst has been observed by Kepler in the short-cadence mode. The light curve contains one superoutburst, eight normal outbursts (including a precursor to the superoutburst), and several minioutbursts that are present before -- but not after -- the superoutburst. The superoutburst began with a precursor normal outburst, and shortly after the peak of the precursor, the system developed large-amplitude superhumps that achieved their maximum amplitude after just three superhump cycles. The period excess of the initial superhump period relative to the orbital period implies a mass ratio of 0.281 +/- 0.015, placing it marginally above most theoretical predictions of the highest-possible mass ratio for superhump formation. In addition, our analysis of the variations in eclipse width and depth, as well as the hot spot amplitudes, generally provides substantiation of the thermal-tidal instability model. The K2 data, in conjunction with our ground-based time-resolved spectroscopy and photometry from 2014-2016, allows us to determine many of the fundamental parameters of this system.

For 80 days in 2017, the Kepler Space Telescope continuously observed the transitional millisecond pulsar system PSR J1023+0038 in its accreting state. We present analyses of the 59 second cadence data, focusing on investigations of the orbital light curve of the irradiated companion star, and of flaring activity in the neutron star's accretion disc. The underlying orbital modulation from the companion star retains a similar amplitude and asymmetric heating profile as seen in previous photometric observations of the system in its radio pulsar state, suggesting that the heating mechanism has not been affected by the state change. We also find tentative evidence that this asymmetry may vary with time. The light curve also exhibits "flickering" activity, evident as short time-scale flux correlations throughout the observations, and periods of rapid mode-switching activity on time scales shorter than the observation cadence. Finally, the system spent 20% of the observations in a flaring state, with the length of these flares varying from 2 minutes up to several hours. The flaring behaviour is consistent with a self-organised criticality mechanism, most likely related to the build up and release of mass at the inner edge of the accretion disc.

We present time-resolved photometry and spectroscopy of the recently classified polar MASTER OT J132104.04+560957.8. The spectrum shows a smooth, non-thermal continuum at the time of maximum light, without any individually discernible cyclotron harmonics. Using homogeneous cyclotron modeling, we interpret this as cyclotron radiation whose individual harmonics have blended together, and on this basis, we loosely constrain the magnetic field strength to be less than ~30 MG. In addition, for about one-tenth of the orbital period, the Balmer and He I emission lines transition into absorption features, with He II developing an absorption core. We use our observations of this phenomenon to test theoretical models of the accretion curtain and conclude that the H and He I lines are produced throughout the curtain, in contravention of theoretical predictions of separate H and He I line-forming regions. Moreover, a significant amount of He II emission originates within the accretion curtain, implying that the curtain is significantly hotter than expected from theory. Finally, we comment on the object's long-term photometry, including evidence that it recently transitioned into a prolonged, exceptionally stable high state following a potentially decades-long low state.

The quiescent periodic photometric modulations of two low-inclination cataclysmic variables observed in Kepler K2 Campaigns 0 and 1, KZ Gem and TW Vir, are investigated. A phase-correcting method was successfully used to detect the orbital modulations of KZ Gem and TW Vir and improve their orbital periods. The light curve morphologies of both CVs were further analyzed by defining flux ratios and creating colormaps. KZ Gem shows ellipsoidal modulations with an orbital period of 0.22242(1) day, twice the period listed in the updated RK catalogue (Edition 7.24). With this newly determined period, KZ Gem is no longer a CV in the period gap, but a long-period CV. A part of the quiescent light curve of TW Vir that had the highest stability was used to deduce its improved orbital period of 0.182682(3) day. The flat patterns shown in the colormaps of the flux ratios for KZ Gem demonstrate the stability of their orbital modulations, while TW Vir show variable orbital modulations during the K2 datasets. In TW Vir, the single versus double-peaked nature of the quiescent orbital variations before and after superoutburst may be related to the effect of the superoutburst on the accretion disk.

We analyze long-cadence Kepler K2 observations of AR Sco from 2014, along with survey photometry obtained between 2005 and 2016 by the Catalina Real-Time Sky Survey and the All-Sky Automated Survey for Supernovae. The K2 data show the orbital modulation to have been fairly stable during the 78 days of observations, but we detect aperiodic deviations from the average waveform with an amplitude of ~2% on a timescale of a few days. A comparison of the K2 data with the survey photometry reveals that the orbital waveform gradually changed between 2005 and 2010, with the orbital maximum shifting to earlier phases. We compare these photometric variations with proposed models of this unusual system.

We present the first ever X-ray data taken of an intermediate polar, FO Aqr, when in a low accretion state and during the subsequent recovery. The Swift and Chandra X-ray data taken during the low accretion state in July 2016 both show a softer spectrum when compared to archival data taken when FO Aqr was in a high state. The X-ray spectrum in the low state showed a significant increase in the ratio of the soft X-ray flux to the hard X-ray flux due to a change in the partial covering fraction of the white dwarf from $>85\%$ to $70^{+5}_{-8}\%$ and a change in the hydrogen column density within the disc from 19$^{+1.2}_{-0.9}\times 10^{22}$ cm$^{-2}$ to 1.3$^{+0.6}_{-0.3}\times 10^{22}$ cm$^{-2}$. XMM-Newton observations of FO Aqr during the subsequent recovery suggest that the system had not yet returned to its typical high state by November 2016, with the hydrogen column density within the disc found to be 15$^{+3.0}_{-2.0}$ cm$^{-2}$. The partial covering fraction varied in the recovery state between $85\%$ and $95\%$. The spin period of the white dwarf in 2014 and 2015 has also been refined to 1254.3342(8) s. Finally, we find an apparent phase difference between the high state X-ray pulse and recovery X-ray pulse of 0.17, which may be related to a restructuring of the X-ray emitting regions within the system.

We present the first X-ray observations of the eclipsing cataclysmic variables Lanning 386 and MASTER OTJ192328.22+612413.5, possible SW Sextantis systems. The X-ray light curve of Lanning 386 shows deep eclipses, similar to the eclipses seen in the optical light curve, confirming the high inclination of the system. There is evidence of a periodicity between 17-22 min in the X-ray and optical light curves of Lanning 386, which is associated with quasi-periodic oscillations. This system also displays a hard X-ray spectrum which is well fit by a partially covered, absorbed 2 temperature plasma. The cool plasma temperature (0.24$^{+0.17}_{-0.08}$ keV) and hot plasma temperature (9$^{+4}_{-2}$ keV) are not atypical plasma temperatures of known intermediate polar systems. Based on this model, we suggest that Lanning 386 is an intermediate polar with a high accretion rate. The hot plasma temperature limits the white dwarf mass to $>$0.5 M$_{\odot}$. From the optical spectrum obtained using the Large Binocular Telescope, we find that the secondary in the system is consistent with an M5V star, and refine the distance to Lanning 386 to be 160$\pm$50 pc. Finally, we use the high time resolution of the optical spectra to crudely constrain the magnetic moment of the white dwarf in Lanning 386. J1923 was not detected during the observations, but the upper limit on the flux is inline with J1923 and Lanning 386 being related.

In 2016 May, the intermediate polar FO~Aqr was detected in a low state for the first time in its observational history. We report time-resolved photometry of the system during its initial recovery from this faint state. Our data, which includes high-speed photometry with cadences of just 2 sec, shows the existence of very strong periodicities at 22.5 min and 11.26 min, equivalent to the spin-orbit beat frequency and twice its value, respectively. A pulse at the spin frequency is also present but at a much lower amplitude than is normally observed in the bright state. By comparing our power spectra with theoretical models, we infer that a substantial amount of accretion was stream-fed during our observations, in contrast to the disk-fed accretion that dominates the bright state. In addition, we find that FO~Aqr's rate of recovery has been unusually slow in comparison to rates of recovery seen in other magnetic cataclysmic variables, with an $e$-folding time of 115$\pm7$ days. The recovery also shows irregular variations in the median brightness of as much as 0.2~mag over a 10-day span. Finally, we show that the arrival times of the spin pulses are dependent upon the system's overall brightness.

We present optical photometry and spectroscopy of the new eclipsing Cataclysmic Variable MASTER OTJ192328.22+612413.5, discovered by the MASTER team. We find the orbital period to be P=0.16764612(5) day /4.023507(1) hour. The depth of the eclipse (2.9$\pm$0.1 mag) suggests that the system is nearly edge on, and modeling of the system confirms the inclination to be between 81.3-83.6 degree. The brightness outside of eclipse varies between observations, with a change of 1.6$\pm$0.1 mag. Spectroscopy reveals double-peaked Balmer emission lines. By using spectral features matching a late M-type companion, we bound the distance to be 750$\pm$250 pc, depending on the companion spectral type. The source displays 2 mag brightness changes on timescales of days. The amplitude of these changes, along with the spectrum at the faint state, suggest the system is possibly a dwarf nova. The lack of any high excitation HeII lines suggests this system is not magnetically dominated. The light curve in both quiescence and outburst resembles that of Lanning 386, implying MASTER OTJ192328.22+612413.5 is a possible cross between a dwarf nova and a SW Sextantis star.

We present photometry of the intermediate polar FO Aquarii obtained as part of the K2 mission using the Kepler space telescope. The amplitude spectrum of the data confirms the orbital period of 4.8508(4) h, and the shape of the light curve is consistent with the outer edge of the accretion disk being eclipsed when folded on this period. The average flux of FO Aquarii changed during the observations, suggesting a change in the mass accretion rate. There is no evidence in the amplitude spectrum of a longer period that would suggest disk precession. The amplitude spectrum also shows the white dwarf spin period of 1254.3401(4) s, the beat period of 1351.329(2) s, and 31 other spin and orbital harmonics. The detected period is longer than the last reported period of 1254.284(16) s, suggesting that FO Aqr is now spinning down, and has a positive $\dot{P}$. There is no detectable variation in the spin period over the course of the K2 observations, but the phase of the spin cycle is correlated with the system brightness. We also find the amplitude of the beat signal is correlated with the system brightness.

There are 15 cataclysmic variables (CVs) observed in the first two campaigns of the K2 mission. In this paper, the eight CVs showing distinct features are analyzed in detail. Among these eight, modulations during quiescence are evident at the known orbital periods in the SU UMa stars QZ Vir and RZ Leo, and at our newly determined orbital periods of 1RXS J0632+2536 and WD 1144+011. The periodogram analysis for the quiescent light curve of QZ Vir reveals multi-period modulations and the coexistence of orbital and superhump periods. The phased orbital light curves for the other 3 CVs in quiescence display wide (about half cycle) and shallow (< 0.5 mag) eclipse features. Besides these modulations, their quiescent light curves reveal several transient events: a sudden decrease of system light in 1RXS J0632+2536, a low level flare-like event in QZ Vir, a short brightening event in RZ Leo and a temporary disappearance of the orbital modulation in WD 1144+011. The two known dwarf novae UV Gem and TW Vir and the CVs USNO-B1.01144-00115322 and CSS130516:111236.7+002807 show outbursts, including 1 complete and 3 incomplete normal outbursts and 2 complete superoutbursts. An incomplete but typical normal outburst confirms the dwarf nova identification of the USNO-B1.01144-00115322. The one complete normal outburst in UV Gem possibly provides the orbital period, since its modulations are shorter than the previously observed superhump period. The superoutbursts of TW Vir and CSS130516:111236.7+002807, along with their corresponding superhump periods, indicate that both objects are SU UMa stars. The derived superhump period of CSS130516:111236.7+002807 is 1.44 hr, implying that this new SU UMa star is close to the period minimum.

The "evolved main-sequence" channel is thought to contribute significantly to the population of AM CVn type systems in the Galaxy, and also to the number of cataclysmic variables detected below the period minimum for hydrogen rich systems. CSS120422:J111127+571239 was discovered by the Catalina Sky Survey in April 2012. Its period was found to be 56 minutes, well below the minimum, and the optical spectrum is clearly depleted in hydrogen relative to helium, but still has two orders of magnitude more hydrogen than AM CVn stars. Doppler tomography of the H$\alpha$ line hinted at a spiral structure existing in the disk. Here we present spectroscopy of CSS120422:J111127+571239 using the COS FUV instrument on the Hubble Space Telescope and using the MODS spectrograph on the Large Binocular Telescope. The UV spectrum shows SiIV, NV and HeII, but no detectable CIV. The anomalous nitrogen/carbon ratio is seen in a small number of other CVs and confirms a unique binary evolution. We also present and compare the optical spectrum of V418 Ser and advocate that it is also an evolved main-sequence system.

The asymmetries observed in the line profiles of solar flares can provide important diagnostics of the properties and dynamics of the flaring atmosphere. In this paper the evolution of the Halpha and Ca II 8542 {\AA} lines are studied using high spatial, temporal and spectral resolution ground-based observations of an M1.1 flare obtained with the Swedish 1-m Solar Telescope. The temporal evolution of the Halpha line profiles from the flare kernel shows excess emission in the red wing (red asymmetry) before flare maximum, and excess in the blue wing (blue asymmetry) after maximum. However, the Ca II 8542 {\AA} line does not follow the same pattern, showing only a weak red asymmetry during the flare. RADYN simulations are used to synthesise spectral line profiles for the flaring atmosphere, and good agreement is found with the observations. We show that the red asymmetry observed in Halpha is not necessarily associated with plasma downflows, and the blue asymmetry may not be related to plasma upflows. Indeed, we conclude that the steep velocity gradients in the flaring chromosphere modifies the wavelength of the central reversal in the Halpha line profile. The shift in the wavelength of maximum opacity to shorter and longer wavelengths generates the red and blue asymmetries, respectively.

PSR J102347.6+003841 is a radio pulsar system with a spin period of 1.69 ms and an orbital period of 4.75 hours. Uniquely, it undergoes periods of transient accretion from its companion star: it occupies an important position in the evolutionary track from X-ray binary to isolated millisecond radio pulsar. Here we present a spectroscopic study of this system showing late-type absorption features which match those of a G2V star. We find a semiamplitude of $286 \pm 3$ kms$^{-1}$ and a best fit orbital period of 0.1980966(1) days. We combine these measurements with optical photometry which suggests the secondary star may be underfilling its Roche lobe by between 15\% and 20\%. We weakly constrain the mass of the neutron star to be $\leq$ 2.2 M$_\odot$ at the 2$\sigma$ level. We also discuss the possible origins of the H$\alpha$ emission line in our template subtracted, averaged spectrum. Finally we present and discuss new optical photometry of J1023 taken during the recent outburst of the system.

The structure and properties of the diffuse interstellar medium (ISM) on small scales, sub-au to 1 pc, are poorly understood. We compare interstellar absorption-lines, observed towards a selection of O- and B-type stars at two or more epochs, to search for variations over time caused by the transverse motion of each star combined with changes in the structure in the foreground ISM. Two sets of data were used: 83 VLT- UVES spectra with approximately 6 yr between epochs and 21 McDonald observatory 2.7m telescope echelle spectra with 6 - 20 yr between epochs, over a range of scales from 0 - 360 au. The interstellar absorption-lines observed at the two epochs were subtracted and searched for any residuals due to changes in the foreground ISM. Of the 104 sightlines investigated with typically five or more components in Na I D, possible temporal variation was identified in five UVES spectra (six components), in Ca II, Ca I and/or Na I absorption-lines. The variations detected range from 7\% to a factor of 3.6 in column density. No variation was found in any other interstellar species. Most sightlines show no variation, with 3{\sigma} upper limits to changes of the order 0.1 - 0.3 dex in Ca II and Na I. These variations observed imply that fine-scale structure is present in the ISM, but at the resolution available in this study, is not very common at visible wavelengths. A determination of the electron densities and lower limits to the total number density of a sample of the sightlines implies that there is no striking difference between these parameters in sightlines with, and sightlines without, varying components.

We investigated the response of the solar atmosphere to non-thermal electron beam heating using the radiative transfer and hydrodynamics modelling code RADYN. The temporal evolution of the parameters that describe the non-thermal electron energy distribution were derived from hard X-ray observations of a particular flare, and we compared the modelled and observed parameters. The evolution of the non-thermal electron beam parameters during the X1.5 solar flare on 2011 March 9 were obtained from analysis of RHESSI X-ray spectra. The RADYN flare model was allowed to evolve for 110 seconds, after which the electron beam heating was ended, and was then allowed to continue evolving for a further 300s. The modelled flare parameters were compared to the observed parameters determined from extreme-ultraviolet spectroscopy. The model produced a hotter and denser flare loop than that observed and also cooled more rapidly, suggesting that additional energy input in the decay phase of the flare is required. In the explosive evaporation phase a region of high-density cool material propagated upward through the corona. This material underwent a rapid increase in temperature as it was unable to radiate away all of the energy deposited across it by the non-thermal electron beam and via thermal conduction. A narrow and high-density ($n_{e} \le 10^{15}$ cm$^{-3}$) region at the base of the flare transition region was the source of optical line emission in the model atmosphere. The collision-stopping depth of electrons was calculated throughout the evolution of the flare, and it was found that the compression of the lower atmosphere may permit electrons to penetrate farther into a flaring atmosphere compared to a quiet Sun atmosphere.

We present spectral classifications for 438 B-type stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 30 Doradus region of the Large Magellanic Cloud. Radial velocities are provided for 307 apparently single stars, and for 99 targets with radial-velocity variations which are consistent with them being spectroscopic binaries. We investigate the spatial distribution of the radial velocities across the 30 Dor region, and use the results to identify candidate runaway stars. Excluding potential runaways and members of two older clusters in the survey region (SL 639 and Hodge 301), we determine a systemic velocity for 30 Dor of 271.6 +/- 12.2 km/s from 273 presumed single stars. Employing a 3-sigma criterion we identify nine candidate runaway stars (2.9% of the single stars with radial-velocity estimates). The projected rotational velocities of the candidate runaways appear to be significantly different to those of the full B-type sample, with a strong preference for either large (>345 km/s) or small (<65 km/s) rotational velocities. Of the candidate runaways, VFTS 358 (classified B0.5: V) has the largest differential radial velocity (-106.9 +/- 16.2 km/s), and a preliminary atmospheric analysis finds a significantly enriched nitrogen abundance of 12+log(N/H) > ~8.5. Combined with a large rotational velocity (vsini = 345 +/- 22 km/s), this is suggestive of past binary interaction for this star.

Differential emission measures (DEMs) during the impulsive phase of solar flares were constructed using observations from the EUV Variability Experiment (EVE) and the Markov-Chain Monte Carlo method. Emission lines from ions formed over the temperature range log T = 5.8 - 7.2 allow the evolution of the DEM to be studied over a wide temperature range at 10s cadence. The technique was applied to several M- and X-class flares, where impulsive phase EUV emission is observable in the disk-integrated EVE spectra from emission lines formed up to 3 - 4 MK, and we use spatially-unresolved EVE observations to infer the thermal structure of the emitting region. For the nine events studied the DEMs exhibited a two component distribution during the impulsive phase, a low temperature component with peak temperature of 1 - 2 MK, and a broad high temperature one from 7 - 30 MK. A bimodal high temperature component is also found for several events, with peaks at 8 and 25 MK during the impulsive phase. The origin of the emission was verified using AIA images to be the flare ribbons and footpoints, indicating that the constructed DEMs represent the spatially-average thermal structure of the chromospheric flare emission during the impulsive phase.

Aims: Projected rotational velocities (\vsini) have been estimated for 334 targets in the VLT-FLAMES Tarantula survey that do not manifest significant radial velocity variations and are not supergiants. They have spectral types from approximately O9.5 to B3. The estimates have been analysed to infer the underlying rotational velocity distribution, which is critical for understanding the evolution of massive stars. Methods: Projected rotational velocities were deduced from the Fourier transforms of spectral lines, with upper limits also being obtained from profile fitting. For the narrower lined stars, metal and non-diffuse helium lines were adopted, and for the broader lined stars, both non-diffuse and diffuse helium lines; the estimates obtained using the different sets of lines are in good agreement. The uncertainty in the mean estimates is typically 4% for most targets. The iterative deconvolution procedure of Lucy has been used to deduce the probability density distribution of the rotational velocities. Results: Projected rotational velocities range up to approximately 450 \kms and show a bi-modal structure. This is also present in the inferred rotational velocity distribution with 25% of the sample having $0\leq$\ve$\leq$100\,\kms and the high velocity component having \ve$\sim 250$\,\kms. There is no evidence from the spatial and radial velocity distributions of the two components that they represent either field and cluster populations or different episodes of star formation. Be-type stars have also been identified. Conclusions: The bi-modal rotational velocity distribution in our sample resembles that found for late-B and early-A type stars. While magnetic braking appears to be a possible mechanism for producing the low-velocity component, we can not rule out alternative explanations.

Temporally-resolved electron density measurements of solar flare plasmas are presented using data from the EUV Variability Experiment (EVE) onboard the Solar Dynamics Observatory (SDO). The EVE spectral range contains emission lines formed between 10^4-10^7 K, including transitions from highly ionized iron (>10 MK). Using three density-sensitive Fe XXI ratios, peak electron densities of 10^(11.2)-10^(12.1) cm^(-3) were found during four X-class flares. While previous measurements of densities at such high temperatures were made at only one point during a flaring event, EVE now allows the temporal evolution of these high-temperature densities to be determined at 10 s cadence. A comparison with GOES data revealed that the peak of the density time profiles for each line ratio correlated well with that of the emission measure time profile for each of the events studied.

Archival observations from the Very Large Array (VLA) at frequencies between 1.4 GHz and 43 GHz of the 6.6-day O6.5-7+O5.5-6 binary Cyg OB2 #5 over 20 years are re-examined. The aim is to determine the location and character of its known variable radio emission. The radio emission consists of a primary component associated with the binary, and a non-thermal source (NE), $0.8" to the NE. This work reveals that NE shows no evidence of variation demonstrating that the variable emission arises in the primary component. With NE constant, the radio flux from the primary can now be well determined for the first time, especially in observations that do not resolve both the primary and NE components. The variable radio emission from the primary has a period of 6.7+/-0.3 years which is described by a simple model of a non-thermal source orbiting within the stellar wind envelope of the binary. Such a model implies the presence of a third, unresolved stellar companion (Star C) orbiting the 6.6-day binary with a period of 6.7 years. The non-thermal emission arises from either a WCR between Star C and the binary system, or possibly from Star C directly. Examination of radial velocity observations suggests reflex motion of the binary due to Star C, for which a mass of 23^+22_-14 solar mass is deduced. Together with the star associated with NE, this implies that Cyg OB2 #5 is a quadruple system.

Fifty archival radio observations of the supergiant binary Cyg OB2 #5 using the Very Large Array over 20 years are re-examined to determine the location and character of the previously detected variable radio emission. The radio emission from the system consists of a primary component that is associated with the binary, and a non-thermal source (NE) that has been ascribed to a wind-collision region (WCR) between the stellar winds of the binary and that of a B-type star (Star D) to the NE. NE shows no evidence of variation in 23 epochs where it is resolved separately from the primary radio component, demonstrating that the variable emission arises in the primary. Since NE is non-variable, the radio flux from the primary can now be well determined for the first time, most especially in observations that do not resolve both the primary and NE components. The variable radio emission from the primary has a period of 6.7+/-0.3 years which is described by a simple model of a non-thermal source orbiting within the stellar wind envelope of the binary. Such a model implies the presence of a third, unresolved stellar companion (Star C) orbiting the 6.6-day binary with a period of 6.7 years. The variable non-thermal emission arises from either a WCR between Star C and the binary system, or possibly from Star C directly. The model gives a mass-loss rate of 3.4 x 10^{-5} solar mass/yr for Cyg OB2 #5, unusually high for an Of supergiant and comparable to that of WR stars, and consistent with an unusually strong He I 1.083-micron emission line, also redolent of WR stars. An examination of radial velocity observations suggests reflex motion of the binary due to Star C. The natures of NE and Star D are also examined. (abridged)