Here, we study the dichotomy of the escaping atmosphere of the newly discovered close-in exoplanet AU Mic b. On one hand, the high EUV stellar flux is expected to cause a strong atmospheric escape in AU Mic b. On the other hand, the wind of this young star is believed to be very strong, which could reduce or even inhibit the planet's atmospheric escape. AU Mic is thought to have a wind mass-loss rate that is up to $1000$ times larger than the solar wind mass-loss rate ($\dot{M}_\odot$). To investigate this dichotomy, we perform 3D hydrodynamics simulations of the stellar wind--planetary atmosphere interactions in the AU Mic system and predict the synthetic Ly-$\alpha$ transits of AU Mic b. We systematically vary the stellar wind mass-loss rate from a `no wind' scenario to up to a stellar wind with a mass-loss rate of $1000~\dot{M}_\odot$. We find that, as the stellar wind becomes stronger, the planetary evaporation rate decreases from $6.5\times 10^{10}$ g/s to half this value. With a stronger stellar wind, the atmosphere is forced to occupy a smaller volume, affecting transit signatures. Our predicted Ly-$\alpha$ absorption drops from $\sim 20\%$, in the case of `no wind' to barely any Ly-$\alpha$ absorption in the extreme stellar wind scenario. Future Ly-$\alpha$ transits could therefore place constraints not only on the evaporation rate of AU Mic b, but also on the mass-loss rate of its host star.

Stellar high-energy radiation (X-ray and extreme ultraviolet, XUV) drives atmospheric escape in close-in exoplanets. Given that stellar irradiation depends on the stellar magnetism and that stars have magnetic cycles, we investigate how cycles affect the evolution of exoplanetary atmospheric escape. Firstly, we consider a hypothetical HD209458b-like planet orbiting the Sun. For that, we implement the observed solar XUV radiation available over one and a half solar cycles in a 1D hydrodynamic escape model of HD209458b. We find that atmospheric escape rates show a cyclic variation (from 7.6 to 18.5 $\times$ 10$^{10}$ g s$^{-1}$), almost proportional to the incident stellar radiation. To compare this with observations, we compute spectroscopic transits in two hydrogen lines. We find non-detectable cyclic variations in Ly$\alpha$ transits. Given the temperature sensitiveness of the H$\alpha$ line, its equivalent width has an amplitude of 1.9 mA variation over the cycle, which could be detectable in exoplanets such as HD209458b. We demonstrate that the XUV flux is linearly proportional to the magnetic flux during the solar cycle. Secondly, we apply this relation to derive the cyclic evolution of the XUV flux of HD189733 using the available magnetic flux observations of the star from Zeeman Doppler Imaging over nearly a decade. The XUV fluxes are then used to model escape in HD189733b, which shows escape rate varying from 2.8 to 6.5 $\times$ 10$^{10}$ g s$^{-1}$. Like in the HD209458b case, this introduces variations in Ly$\alpha$ and H$\alpha$ transits, with H$\alpha$ variations more likely to be observable. Finally, we show that a strong stellar flare would enhance significantly Ly$\alpha$ and H$\alpha$ transit depths.

We report the discovery of a warm Neptune and a hot sub-Neptune transiting TOI-421 (BD-14 1137, TIC 94986319), a bright (V=9.9) G9 dwarf star in a visual binary system observed by the TESS space mission in Sectors 5 and 6. We performed ground-based follow-up observations -- comprised of LCOGT transit photometry, NIRC2 adaptive optics imaging, and FIES, CORALIE, HARPS, HIRES, and PFS high-precision Doppler measurements -- and confirmed the planetary nature of the 16-day transiting candidate announced by the TESS team. We discovered an additional radial velocity signal with a period of 5 days induced by the presence of a second planet in the system, which we also found to transit its host star. We found that the outer warm Neptune, TOI-421 b, has an orbital period of Pb=16.06815+0.00034-0.00035 days, a mass of Mb=16.23+1.14-1.08 MEarth, a radius of Rb=5.17+-0.13 REarth and a density of rho_b=0.645+0.069-0.062 g/cm3, whereas the inner mini-Neptune, TOI-421 c, has a period of Pc=5.19676+0.00049-0.00048 days, a mass of Mc=7.05+0.71-0.70 MEarth and a radius of Rc=2.72+0.19-0.18 REarth. With its characteristics the latter planet (rho_c=1.93+0.49-0.39 g/cm3) is placed in the intriguing class of the super-puffy mini-Neptunes. TOI-421 b and c are found to be well suitable for atmospheric characterization. Our atmospheric simulations predict significant Ly-alpha transit absorption, due to strong hydrogen escape in both planets, and the presence of detectable CH4 in the atmosphere of TOI-421 b if equilibrium chemistry is assumed.

The growth of spectroscopic observations of exoplanetary systems allows the possibility of testing theoretical models and studying the interaction that exoplanetary atmospheres have with the wind and the energetic photons from the star. In this work, we present a set of numerical 3D simulations of HD 209458b for which spectral lines observations of their evaporative atmosphere are available. The different simulations aim to reproduce different scenarios for the star-planet interaction. With our models, we reconstruct the Ly$\alpha$ line during transit and compare with observations. The results allow us to analyse the shape of the line profile under these different scenarios and the comparison with the observations suggests that HD209458b may have a magnetic field off less than 1 G. We also explore the behaviour of the magnesium lines for models with and without magnetic fields.

We homogeneously analyse $\sim 3.2\times 10^5$ photometric measurements for $\sim 1100$ transit lightcurves belonging to $17$ exoplanet hosts. The photometric data cover $16$ years 2004--2019 and include amateur and professional observations. Old archival lightcurves were reprocessed using up-to-date exoplanetary parameters and empirically debiased limb-darkening models. We also derive self-consistent transit and radial-velocity fits for $13$ targets. We confirm the nonlinear TTV trend in the WASP-12 data at a high significance, and with a consistent magnitude. However, Doppler data reveal hints of a radial acceleration about $(-7.5\pm 2.2)$~m/s/yr, indicating the presence of unseen distant companions, and suggesting that roughly $10$ per cent of the observed TTV was induced via the light-travel (or Roemer) effect. For WASP-4, a similar TTV trend suspected after the recent TESS observations appears controversial and model-dependent. It is not supported by our homogeneus TTV sample, including $10$ ground-based EXPANSION lightcurves obtained in 2018 simultaneously with TESS. Even if the TTV trend itself does exist in WASP-4, its magnitude and tidal nature are uncertain. Doppler data cannot entirely rule out the Roemer effect induced by possible distant companions.

Lyman $\alpha$ observations of the transiting exoplanet HD 209458b enable the study of exoplanets exospheres exposed to stellar EUV fluxes, as well as the interacting stellar wind properties. In this study we present 3D hydrodynamical models for the stellar-planetary wind interaction including radiation pressure and charge exchange, together with photoionization, recombination and collisional ionization processes. Our models explore the contribution of the radiation pressure and charge exchange on the Ly$\alpha$ absorption profile in a hydrodynamical framework, and for a single set of stellar wind parameters appropriate for HD 209458. We find that most of the absorption is produced by the material from the planet, with a secondary contribution of neutralized stellar ions by charge exchange. At the same time, the hydrodynamic shock heats up the planetary material, resulting in a broad thermal profile. Meanwhile, the radiation pressure yielded a small velocity shift of the absorbing material. While neither charge exchange nor radiation pressure provide enough neutrals at the velocity needed to explain the observations at $-100~\mathrm{km~s^{-1}}$ individually, we find that the two effects combined with the broad thermal profile are able to explain the observations.

Although the present-day Sun rotates too slowly to exhibit centrifugally-supported 'slingshot prominences', at some time during its past it may have formed these clouds of cool gas and ejected them into the interplanetary medium. We determine the time period for this behaviour by using a rotation evolution code to derive the properties of the formation and ejection of slingshot prominences during the lifetime of a star similar to our Sun. The mass, mass loss rate and rate of ejection of these prominences are calculated using the analytical expression derived in our previous work. We find that for stars with an initial rotation rate larger than $4.6\, \Omega_\odot$, about half of all solar mass stars, slingshot prominences will be present even after the star reaches the main sequence phase. In a fast rotator, this means that prominences can form until the star reaches $\sim 800$ Myr old. Our results also indicate that the mass and lifetime of this type of prominence have maximum values when the star reaches the ZAMS at an age of $\sim 40$ Myr for a solar mass star.

Lyman $\alpha$ observations during an exoplanet transit have proved to be very useful to study the interaction between the stellar wind and the planetary atmosphere. They have been extensively used to constrain planetary system parameters that are not directly observed, such as the planetary mass loss rate. In this way, Ly $\alpha$ observations can be a powerful tool to infer the existence of a planetary magnetic field, since it is expected that the latter will affect the escaping planetary material. To explore the effect that magnetic fields have on the Ly $\alpha$ absorption of HD 209458b, we run a set of 3D MHD simulations including dipolar magnetic fields for the planet and the star. We assume values for the surface magnetic field at the poles of the planet in the range of [0-5] G, and from 1 to 5 G at the poles of the star. Our models also include collisional and photo-ionisation, radiative recombination, and an approximation for the radiation pressure. Our results show that the magnetic field of the planet and the star change the shape of the Ly $\alpha$ absorption profile, since it controls the extent of the planetary magnetosphere and the amount of neutral material inside it. The model that best reproduces the absorption observed in HD 209458b (with canonical values for the stellar wind parameters) corresponds to a dipole planetary field of <~ 1 gauss at the poles.

The observational signatures of prominences have been detected in single and binary G and K type stars for many years now, but recently this has been extended to the M dwarf regime. Prominences carry away both mass and angular momentum when they are ejected and the impact of this mass on any orbiting planets may be important for the evolution of exoplanetary atmospheres. By means of the classification used in the massive star community, that involves knowledge of two parameters (the co-rotation and Alfv\'en radii, $r_K$ and $r_A$), we have determined which cool stars could support prominences. From a model of mechanical support, we have determined that the prominence mass $m_p/M_\star=(E_M/E_G)(r_\star/r_K)^2 F$ where $E_MB_\star^2r_\star^3$ and $E_G = GM_\star^2/r_\star$ are magnetic and gravitational energies and $F$ is a geometric factor. Our calculated masses and ejection frequencies (typically $10^{16}-10^{17}$g and 0.4 d, respectively) are consistent with observations and are sufficient to ensure that an exoplanet orbiting in the habitable zone of an M dwarf could suffer frequent impacts.

White dwarfs are often found in binary systems with orbital periods ranging from tens of minutes to hours in which they can accrete gas from their companion stars. In about 15% of these binaries, the magnetic field of the white dwarf is strong enough ($\geq 10^6$ Gauss) to channel the accreted matter along field lines onto the magnetic poles. The remaining systems are referred to as "non-magnetic", since to date there has been no evidence that they have a dynamically significant magnetic field. Here we report an analysis of archival optical observations of the "non-magnetic" accreting white dwarf in the binary system MV Lyrae (hereafter MV Lyr), whose lightcurve displayed quasi-periodic bursts of $\approx 30$ minutes duration every $\approx 2$ hours. The observations indicate the presence of an unstable magnetically-regulated accretion mode, revealing the existence of magnetically gated accretion, where disk material builds up around the magnetospheric boundary (at the co-rotation radius) and then accretes onto the white dwarf, producing bursts powered by the release of gravitational potential energy. We infer a surface magnetic field strength for the white dwarf in MV Lyr between $2 \times 10^4 \leq B \leq 10^5$ Gauss, too low to be detectable by other current methods. Our discovery provides a new way of studying the strength and evolution of magnetic fields in accreting white dwarfs and extends the connections between accretion onto white dwarfs, young stellar objects and neutron stars, for which similar magnetically gated accretion cysles have been identified.

We present coordinated Chandra X-ray Observatory and Karl G. Jansky Very Large Array observations of the transitional millisecond pulsar PSR J1023+0038 in its low-luminosity accreting state. The unprecedented five hours of strictly simultaneous X-ray and radio continuum coverage for the first time unambiguously show a highly reproducible, anti-correlated variability pattern. The characteristic switches from the X-ray high mode into a low mode are always accompanied by a radio brightening with duration that closely matches the X-ray low mode interval. This behavior cannot be explained by a canonical inflow/outflow accretion model where the radiated emission and the jet luminosity are powered by, and positively correlated with, the available accretion energy. We interpret this phenomenology as alternating episodes of low-level accretion onto the neutron star during the X-ray high mode that are interrupted by rapid ejections of plasma by the active rotation-powered pulsar, possibly initiated by a reconfiguration of the pulsar magnetosphere, that cause a transition to a less luminous X-ray mode. The observed anti-correlation between radio and X-ray luminosity has an additional consequence: transitional MSPs can make excursions into a region of the radio/X-ray luminosity plane previously thought to be occupied solely by black hole X-ray binary sources. This complicates the use of this luminosity relation to identify candidate black holes, suggesting the need for additional discriminants when attempting to establish the true nature of the accretor.

Strongly magnetized accreting stars are often hypothesized to be in `spin equilibrium' with their surrounding accretion flows, which requires that the accretion rate changes more slowly than it takes the star to reach spin equilibrium. This is not true for most magnetically accreting stars, which have strongly variable accretion outbursts on time-scales much shorter than the time it would take to reach spin equilibrium. This paper examines how accretion outbursts affect the time a star takes to reach spin equilibrium and its final equilibrium spin period. I consider several different models for angular momentum loss -- either carried away in an outflow, lost to a stellar wind, or transferred back to the accretion disc (the `trapped disc'). For transient sources, the outflow scenario leads to significantly longer times to reach spin equilibrium ($\sim$10x), and shorter equilibrium spin periods than would be expected from spin equilibrium arguments, while the `trapped disc' does not. The results suggest that disc trapping plays a significant role in the spin evolution of strongly magnetic stars, with some caveats for young stellar objects.

In outburst, neutron star X-ray binaries produce less powerful jets than black holes at a given X-ray luminosity. This has made them more difficult to study as they fade towards quiescence. To explore whether neutron stars power jets at low accretion rates ($L_{\rm X} \lesssim 10^{36}$ erg s$^{-1}$), we investigate the radio and X-ray properties of three accreting millisecond X-ray pulsars (IGR J17511-3057, SAX J1808.4-3658 and IGR J00291+5934) during their outbursts in 2015, and of the non-pulsing neutron star Cen X-4 in quiescence (2015) and in outburst (1979). We did not detect the radio counterpart of IGR J17511-3057 in outburst or of Cen X-4 in quiescence, but did detect IGR J00291+5934 and SAX J1808.4-3658, showing that at least some neutron stars launch jets at low accretion rates. While the radio and X-ray emission in IGR J00291+5934 seem to be tightly correlated, the relationship in SAX J1808.4-3658 is more complicated. We find that SAX J1808.4-3658 produces jets during the reflaring tail, and we explore a toy model to ascertain whether the radio emission could be attributed to the onset of a strong propeller. The lack of a universal radio/X-ray correlation, with different behaviours in different neutron star systems (with various radio/X-ray correlations; some being radio faint and others not), points at distinct disk-jet interactions in individual sources, while always being fainter in the radio band than black holes at the same X-ray luminosity.

We observed the Rapid Burster with Chandra when it was in the 'banana' state that usually precedes the type-II X-ray bursting 'island' state for which the source is particularly known. We employed the High-Energy Transmission Grating Spectrometer in combination with the ACIS-S detector in continuous clocking mode. The observation yielded 20 thermonuclear type-I X-ray bursts emitted from the neutron star surface with recurrence times between 0.9 and 1.2 hr, and an e-folding decay time scale of 1 min. We searched for narrow spectral features in the burst emission that could constrain the composition of the ashes of the nuclear burning and the compactness of the neutron star, but found none. The upper limit on the equivalent width of narrow absorption lines between 2 and 6 keV is between 5 and 20 eV (single trial 3 sigma confidence level) and on those of absorption edges between 150 and 400 eV. The latter numbers are comparable to the levels predicted by Weinberg, Bildsten & Schatz (2006) for Eddington-limited thermonuclear bursts.

The low mass X-ray binary 2A 1822-371 is an eclipsing system with an accretion disc corona and with an orbital period of 5.57 hr. The primary is an 0.59 s X-ray pulsar with a proposed strong magnetic field of 10^10-10^12 G. In this paper we study the spin evolution of the pulsar and constrain the geometry of the system. We find that, contrary to previous claims, a thick corona is not required, and that the system characteristics could be best explained by a thin accretion outflow due to a super-Eddington mass transfer rate and a geometrically thick inner accretion flow. The orbital, spectral and timing observations can be reconciled in this scenario under the assumption that the mass transfer proceeds on a thermal timescale which would make 2A 1822-371, a mildly super-Eddington source viewed at high inclination angles. The timing analysis on 13 years of RXTE data show a remarkably stable spin-up that implies that 2A 1822-371, might quickly turn into a millisecond pulsar in the next few thousand years.

Transitional millisecond pulsars (tMSPs) switch, on roughly multi-year timescales, between rotation-powered radio millisecond pulsar (RMSP) and accretion-powered low-mass X-ray binary (LMXB) states. The tMSPs have raised several questions related to the nature of accretion flow in their LMXB state and the mechanism that causes the state switch. The discovery of coherent X-ray pulsations from PSR J1023+0038 (while in the LMXB state) provides us with the first opportunity to perform timing observations and to compare the neutron star's spin variation during this state to the measured spin-down in the RMSP state. Whereas the X-ray pulsations in the LMXB state likely indicate that some material is accreting onto the neutron star's magnetic polar caps, radio continuum observations indicate the presence of an outflow. The fraction of the inflowing material being ejected is not clear, but it may be much larger than that reaching the neutron star's surface. Timing observations can measure the total torque on the neutron star. We have phase-connected nine XMM-Newton observations of PSR J1023+0038 over the last 2.5 years of the LMXB state to establish a precise measurement of spin evolution. We find that the average spin-down rate as an LMXB is 26.8+/-0.4% faster than the rate (-2.39x10^-15 Hz s-1) determined during the RMSP state. This shows that negative angular momentum contributions (dipolar magnetic braking and outflow) exceed positive ones (accreted material), and suggests that the pulsar wind continues to operate at a largely unmodified level. We discuss implications of this tight observational constraint in the context of possible accretion models.

Close-in hot Jupiters are exposed to a tremendous photon flux that ionizes the neutral escaping material from the planet leaving an observable imprint that makes them an interesting laboratory for testing theoretical models. In this work we present 3D hydrodynamic simulations with radiation transfer calculations of a close-in exoplanet in a blow-off state. We calculate the Ly-$\alpha$ absorption and compare it with observations of HD 209458b an previous simplified model results.Our results show that the hydrodynamic interaction together with a proper calculation of the photoionization proccess are able to reproduce the main features of the observed Ly-$\alpha$ absorption, in particular at the blue-shifted wings of the line. We found that the ionizing stellar flux produce an almost linear effect on the amount of absorption in the wake. Varying the planetary mass loss rate and the radiation flux, we were able to reproduce the $10\%$ absorption observed at $-100~\mathrm{km~s^{-1}}$.

Accreting millisecond X-ray pulsars are an important subset of low-mass X-ray binaries in which coherent X-ray pulsations can be observed during occasional, bright outbursts (X-ray luminosity $L_X\sim 10^{36}$ erg s$^{-1}$). These pulsations show that matter is being channeled onto the neutron star's magnetic poles. However, such sources spend most of their time in a low-luminosity, quiescent state ($L_X\lesssim 10^{34}$ erg s$^{-1}$), where the nature of the accretion flow onto the neutron star (if any) is not well understood. Here we report that the millisecond pulsar/low-mass X-ray binary transition object PSR J1023+0038 intermittently shows coherent X-ray pulsations at luminosities nearly 100 times fainter than observed in any other accreting millisecond X-ray pulsar. We conclude that in spite of its low luminosity PSR J1023+0038 experiences episodes of channeled accretion, a discovery that challenges existing models for accretion onto magnetized neutron stars.

A recent joint XMM-Newton/Nuclear Spectroscopic Telescope Array (NuSTAR) observation of the accreting neutron star Cen X-4 ($L_{\rm X}\sim10^{33}{\rm~erg~s}^{-1}$) revealed a hard power-law component ($\Gamma\sim1$-$1.5$) with a relatively low cut-off energy (~10 keV), suggesting bremsstrahlung emission. The physical requirements for bremsstrahlung combined with other observed properties of Cen X-4 suggest the emission comes from a boundary layer rather than the accretion flow. The accretion flow itself is thus undetected (with an upper limit of $L_{\rm flow}\lesssim0.3 L_{\rm X}$). A deep search for coherent pulsations (which would indicate a strong magnetic field) places a 6 per cent upper limit on the fractional amplitude of pulsations, suggesting the flow is not magnetically regulated. Considering the expected energy balance between the accretion flow and the boundary layer for different values of the neutron star parameters (size, magnetic field, and spin) we use the upper limit on $L_{\rm flow}$ to set an upper limit of $\varepsilon\lesssim0.3$ for the intrinsic radiative efficiency of the accretion flow for the most likely model of a fast-spinning, non-magnetic neutron star. The non-detection of the accretion flow provides the first direct evidence that this flow is indeed 'radiatively inefficient', i.e. most of the gravitational potential energy lost by the flow before it hits the star is not emitted as radiation.

The accreting millisecond X-ray pulsar SAX J1808.4--3658 shows peculiar low luminosity states known as "reflares" after the end of the main outburst. During this phase the X-ray luminosity of the source varies by up to three orders of magnitude in less than 1-2 days. The lowest X-ray luminosity observed reaches a value of ~1e32 erg/s, only a factor of a few brighter than its typical quiescent level. We investigate the 2008 and 2005 reflaring state of SAX J1808.4-3658 to determine whether there is any evidence for a change in the accretion flow with respect to the main outburst. We perform a multiwavelength photometric and spectral study of the 2005 and 2008 reflares with data collected during an observational campaign covering the near-infrared, optical, ultra-violet and X-ray band. We find that the NIR/optical/UV emission, expected to some from the outer accretion disk shows variations in luminosity which are 1--2 orders of magnitude shallower than in X-rays. The X-ray spectral state observed during the reflares does not change substantially with X-ray luminosity indicating a rather stable configuration of the accretion flow. We investigate the most likely configuration of the innermost regions of the accretion flow and we infer an accretion disk truncated at or near the co-rotation radius. We interpret these findings as due to either a strong outflow (due to a propeller effect) or a trapped disk (with limited/no outflow) in the inner regions of the accretion flow.

The PSR J1023+0038 binary system hosts a neutron star and a low-mass, main-sequence-like star. It switches on year timescales between states as an eclipsing radio millisecond pulsar and a low-mass X-ray binary. We present a multi-wavelength observational campaign of PSR J1023+0038 in its most recent low-mass X-ray binary state. Two long XMM-Newton observations reveal that the system spends ~70% of the time in a $\approx$$3\times10^{33}$ erg/s X-ray luminosity mode, which, as shown in Archibald et al. (2014), exhibits coherent X-ray pulsations. This emission is interspersed with frequent lower flux mode intervals with $\approx$$5\times 10^{32}$ erg/s and sporadic flares reaching up to $\approx$$10^{34}$ erg/s, with neither mode showing significant X-ray pulsations. The switches between the three flux modes occur on timescales of order 10 s. In the UV and optical, we observe occasional intense flares coincident with those observed in X-rays. Our radio timing observations reveal no pulsations at the pulsar period during any of the three X-ray modes, presumably due to complete quenching of the radio emission mechanism by the accretion flow. Radio imaging detects highly variable, flat-spectrum continuum emission from PSR J1023+0038, consistent with an origin in a weak jet-like outflow. Our concurrent X-ray and radio continuum data sets do not exhibit any correlated behavior. The observational evidence we present bears qualitative resemblance to the behavior predicted by some existing propeller and trapped disk accretion models although none can fully account for all aspects of the rich phenomenology of this system.

We perform an analysis of ~80000 photometric measurements for the following 10 stars hosting transiting planets: WASP-2, -4, -5, -52, Kelt-1, CoRoT-2, XO-2, TrES-1, HD 189733, GJ 436. Our analysis includes mainly transit lightcurves from the Exoplanet Transit Database, public photometry from the literature, and some proprietary photometry privately supplied by other authors. Half of these lightcurves were obtained by amateurs. From this photometry we derive 306 transit timing measurements, as well as improved planetary transit parameters. Additionally, for 6 of these 10 stars we present a set of radial velocity measurements obtained from the spectra stored in the HARPS, HARPS-N, and SOPHIE archives using the HARPS-TERRA pipeline. Our analysis of these TTV and RV data did not reveal significant hints of additional orbiting bodies in almost all of the cases. In the WASP-4 case, we found hints of marginally significant TTV signals having amplitude 10-20 sec, although their parameters are model-dependent and uncertain, while radial velocities did not reveal statistically significant Doppler signals.

Type II bursts are thought to arise from instabilities in the accretion flow onto a neutron star in an X-ray binary. Despite having been known for almost 40 years, no model can yet satisfactorily account for all their properties. To shed light on the nature of this phenomenon and provide a reference for future theoretical work, we study the entire sample of Rossi X-ray Timing Explorer data of type II bursts from the Rapid Burster (MXB 1730-335). We find that type II bursts are Eddington-limited in flux, that a larger amount of energy goes in the bursts than in the persistent emission, that type II bursts can be as short as 0.130 s, and that the distribution of recurrence times drops abruptly below 15-18 s. We highlight the complicated feedback between type II bursts and the NS surface thermonuclear explosions known as type I bursts, and between type II bursts and the persistent emission. We review a number of models for type II bursts. While no model can reproduce all the observed burst properties and explain the source uniqueness, models involving a gating role for the magnetic field come closest to matching the properties of our sample. The uniqueness of the source may be explained by a special combination of magnetic field strength, stellar spin period and alignment between the magnetic field and the spin axis.

The discovery of quasi-periodic oscillations (QPOs) in magnetar giant flares has opened up prospects for neutron star asteroseismology. The scarcity of giant flares makes a search for QPOs in the shorter, far more numerous bursts from Soft Gamma Repeaters (SGRs) desirable. In Huppenkothen et al (2013), we developed a Bayesian method for searching for QPOs in short magnetar bursts, taking into account the effects of the complicated burst structure, and have shown its feasibility on a small sample of bursts. Here, we apply the same method to a much larger sample from a burst storm of 286 bursts from SGR J1550-5418. We report a candidate signal at 260 Hz in a search of the individual bursts, which is fairly broad. We also find two QPOs at 93 Hz and one at 127 Hz, when averaging periodograms from a number of bursts in individual triggers, at frequencies close to QPOs previously observed in magnetar giant flares. Finally, for the first time, we explore the overall burst variability in the sample, and report a weak anti-correlation between the power-law index of the broadband model characterising aperiodic burst variability, and the burst duration: shorter bursts have steeper power law indices than longer bursts. This indicates that longer bursts vary over a broader range of time scales, and are not simply longer versions of the short bursts.

We present a 3D hydrodynamic study of the effects that different stellar wind conditions and planetary wind structures have on the calculated Ly-$\alpha$ absorptions produced during the transit of HD 209458b. Considering a range of stellar wind speeds $\sim$[350-800] km s$^{-1}$, coronal temperature $\sim$[3-7] $\times10^{6}$ K and two values of the polytropic index $\Gamma$ $\sim$[1.01-1.13], while keeping fixed the stellar mass loss rate, we found a that a $\dot M_p$ range between $\sim$[3-5] $\times 10^{10}$g s$^{-1}$ give account for the observational absorption in Ly-$\alpha$ measured for the planetary system. Also, several models with anisotropic evaporation profiles for the planetary escaping atmosphere were carried out, showing that both, the escape through polar regions and through the night side yields larger absorptions than an isotropic planetary wind.

We conducted the first long-term (60 days), multiwavelength (optical, ultraviolet, and X-ray) simultaneous monitoring of Cen X-4 with daily Swift observations, with the goal of understanding variability in the low mass X-ray binary Cen X-4 during quiescence. We found Cen X-4 to be highly variable in all energy bands on timescales from days to months, with the strongest quiescent variability a factor of 22 drop in the X-ray count rate in only 4 days. The X-ray, UV and optical (V band) emission are correlated on timescales down to less than 110 s. The shape of the correlation is a power law with index gamma about 0.2-0.6. The X-ray spectrum is well fitted by a hydrogen NS atmosphere (kT=59-80 eV) and a power law (with spectral index Gamma=1.4-2.0), with the spectral shape remaining constant as the flux varies. Both components vary in tandem, with each responsible for about 50% of the total X-ray flux, implying that they are physically linked. We conclude that the X-rays are likely generated by matter accreting down to the NS surface. Moreover, based on the short timescale of the correlation, we also unambiguously demonstrate that the UV emission can not be due to either thermal emission from the stream impact point, or a standard optically thick, geometrically thin disc. The spectral energy distribution shows a small UV emitting region, too hot to arise from the accretion disk, that we identified as a hot spot on the companion star. Therefore, the UV emission is most likely produced by reprocessing from the companion star, indeed the vertical size of the disc is small and can only reprocess a marginal fraction of the X-ray emission. We also found the accretion disc in quiescence to likely be UV faint, with a minimal contribution to the whole UV flux.

We conducted the first long-term (60 days), multiwavelength (optical, ultraviolet, and X-ray) simultaneous monitoring of Cen X-4 with daily Swift observations from June to August 2012, with the goal of understanding variability in the low mass X-ray binary Cen X-4 during quiescence. We found Cen X-4 to be highly variable in all energy bands on timescales from days to months, with the strongest quiescent variability a factor of 22 drop in the X-ray count rate in only 4 days. The X-ray, UV and optical (V band) emission are correlated on timescales down to less than 110 s. The shape of the correlation is a power law with index gamma about 0.2-0.6. The X-ray spectrum is well fitted by a hydrogen NS atmosphere (kT=59-80 eV) and a power law (with spectral index Gamma=1.4-2.0), with the spectral shape remaining constant as the flux varies. Both components vary in tandem, with each responsible for about 50% of the total X-ray flux, implying that they are physically linked. We conclude that the X-rays are likely generated by matter accreting down to the NS surface. Moreover, based on the short timescale of the correlation, we also unambiguously demonstrate that the UV emission can not be due to either thermal emission from the stream impact point, or a standard optically thick, geometrically thin disc. The spectral energy distribution shows a small UV emitting region, too hot to arise from the accretion disk, that we identified as a hot spot on the companion star. Therefore, the UV emission is most likely produced by reprocessing from the companion star, indeed the vertical size of the disc is small and can only reprocess a marginal fraction of the X-ray emission. We also found the accretion disc in quiescence to likely be UV faint, with a minimal contribution to the whole UV flux.

We investigate the possibility of Photospheric Radius Expansion (PRE) during magnetar bursts. Identification of PRE would enable a determination of the magnetic Eddington limit (which depends on field strength and neutron star mass and radius), and shed light on the burst mechanism. To do this we model hydrostatic atmospheres in a strong radial magnetic field, determining both their maximum extent and photospheric temperatures. We find that spatially-extended atmospheres cannot exist in such a field configuration: typical maximum extent for magnetar-strength fields is ~10 m (as compared to 200 km in the non-magnetic case). Achieving balance of gravitational and radiative forces over a large range of radii, which is critical to the existence of extended atmospheres, is rendered impossible in strong fields due to the dependence of opacities on temperature and field strength. We conclude that high luminosity bursts in magnetars do not lead to expansion and cooling of the photosphere, as in the non-magnetic case. We also find the maximum luminosity that can propagate through a hydrostatic magnetar atmosphere to be lower than previous estimates. The proximity and small extent of the photospheres associated with the two different polarization modes also calls into question the interpretation of two blackbody fits to magnetar burst spectra as being due to extended photospheres.

The dynamics of the plasma in the inner regions of an accretion disk around accreting millisecond X-ray pulsars is controlled by the magnetic field of the neutron star. The interaction between an accretion disk and a strong magnetic field is not well-understood, particularly at low accretion rates (the so-called ``propeller regime'). This is due in part to the lack of clear observational diagnostics to constrain the physics of the disk-field interaction. Here we associate the strong ~1 Hz modulation seen in the accreting millisecond X-ray pulsar NGC 6440 X-2 with an instability that arises when the inner edge of the accretion disk is close to the corotation radius (where the stellar rotation rate matches the Keplerian speed in the disk). A similar modulation has previously been observed in another accreting millisecond X-ray pulsar (SAX J1808.4-3658) and we suggest that the two phenomena are related and that this may be a common phenomenon among other magnetized systems. Detailed comparisons with theoretical models suggest that when the instability is observed, the interaction region between the disk and the field is very narrow -- of the order of 1 km. Modelling further suggests that there is a transition region (~1-10 km) around the corotation radius where the disk-field torque changes sign from spin up to spin down. This is the first time that a direct observational constraint has been placed on the width of the disk-magnetosphere interaction region, in the frame of the trapped-disk instability model.

I present recent results studying flare emission in magnetars. Strong quasi-periodic oscillations observed in the tail of giant magnetar flares are frequently interpreted as evidence for global seismic oscillations. I demonstrate that such a global oscillation is not directly observable in the lightcurve. New work suggests the amplitude for the strongest QPO stays nearly constant in the rotation phases where it is observed, which I argue suggests it is produced by an additional emission process from the star.

Strong quasi-periodic oscillations in the tails of the giant gamma-ray flares seen in SGR 1806-20 and SGR 1900+14 are thought to be produced by starquakes in the flaring magnetar. However, the large fractional amplitudes (up to ~20%) observed are difficult to reconcile with predicted amplitudes of starquakes. Here we demonstrate that the steeply pulsed emission profile in the tail of the giant flare can enhance the observed amplitude of the underlying oscillation, analogously to a beam of light oscillating in and out of the line of sight. This mechanism will also broaden the feature in the power spectrum and introduce power at harmonics of the oscillation. The observed strength of the oscillation depends on the amplitude of the underlying starquake, the orientation and location of the emission on the surface of the star, and the gradient of the light curve profile. While the amplification of the signal can be significant, we demonstrate that even with uncertainties in the emission geometry, this effect is not sufficient to produce the observed QPOs. This result excludes the direct observation of a starquake, and suggests that the observed variations come from modulations in the intensity of the emission.

We show that discs accreting onto the magnetosphere of a rotating star can end up in a 'trapped' state, in which the inner edge of the disc stays near the corotation radius, even at low and varying accretion rates. The accretion in these trapped states can be steady or cyclic; we explore these states over wide range of parameter space. We find two distinct regions of instability, one related to the buildup and release of mass in the disk outside corotation, the other to mass storage within the transition region near corotation. With a set of calculations over long time scales we show how trapped states evolve from both nonaccreting and fully accreting initial conditions, and also calculate the effects of cyclic accretion on the spin evolution of the star. Observations of cycles such as found here would provide important clues on the physics of magnetospheric accretion. Recent observations of cyclic and other unusual variability in T Tauri stars (EXors) and X-ray binaries are discussed in this context.

In this Letter we re-examine the idea that gravitational waves are required as a braking mechanism to explain the observed maximum spin-frequency of neutron stars. We show that for millisecond X-ray pulsars, the existence of spin equilibrium as set by the disk/magnetosphere interaction is sufficient to explain the observations. We show as well that no clear correlation exists between the neutron star magnetic field B and the X-ray outburst luminosity Lx when considering an enlarged sample size of millisecond X-ray pulsars.

We investigate the evolution of a thin viscous disc surrounding magnetic star, including the spindown of the star by the magnetic torques it exerts on the disc. The transition from an accreting to a non-accreting state, and the change of the magnetic torque across the corotation radius are included in a generic way, the widths of the transition taken in the range suggested by numerical simulations. In addition to the standard accreting state, two more are found. An accreting state can develop into a 'dead' disc state, with inner edge well outside corotation. More often, a 'trapped' state develops, in which the inner disc edge stays close to corotation even at very low accretion rates. The long-term evolution of these two states is different. In the dead state the star spins down incompletely, retaining much of its initial spin. In the trapped state the star asymptotically can spin down to arbitarily low rates, its angular momentum transferred to the disc. We identify these outcomes with respectively the rapidly rotating and the very slowly rotating classes of Ap stars and magnetic white dwarfs.

Some accreting neutron stars and young stars show unexplained episodic flares in the form of quasi-periodic oscillations or recurrent outbursts. In a series of two papers we present new work on an instability that can lead to episodic outbursts when the accretion disc is truncated by the star's strong magnetic field close to the corotation radius (where the Keplerian frequency matches the star's rotational frequency). In this paper we outline the physics of the instability and use a simple parameterization of the disc-field interaction to explore the instability numerically, which we show can lead to repeated bursts of accretion as well as steady-state solutions, as first suggested by Sunyaev and Shakura. The cycle time of these bursts increases with decreasing accretion rate. These solutions show that the usually assumed `propeller' state, in which mass is ejected from the system, does not need to occur even at very low accretion rates.

Recent observations of two black hole candidates (GX 339-4 and J1753.5-0127) in the low-hard state (L_X/L_Edd ~ 0.003-0.05) suggest the presence of a cool accretion disk very close to the innermost stable orbit of the black hole. This runs counter to models of the low-hard state in which the cool disk is truncated at a much larger radius. We study the interaction between a moderately truncated disk and a hot inner flow. Ion-bombardment heats the surface of the disk in the overlap region between a two-temperature advection-dominated accretion flow and standard accretion disk, producing a hot (kT_e ~70 keV) layer on the surface of the cool disk. The hard X-ray flux from this layer heats the inner parts of the underlying cool disk, producing a soft X-ray excess. Together with interstellar absorption these effects mimic the thermal spectrum from a disk extending to the last stable orbit. The results show that soft excesses in the low-hard state are a natural feature of truncated disk models.

Precision timing is the key ingredient of ongoing pulsar-based gravitational wave searches and tests of general relativity using binary pulsars. The conventional approach to timing explicitly assumes that the radio emitting region is located at the center of the pulsar, while polarimetric observations suggest that radio emission is in fact produced at altitudes ranging from tens to thousands of kilometers above the neutron star surface. Here we present a calculation of the effects of finite emission height on the timing of binary pulsars using a simple model for the emitting region geometry. Finite height of emission changes the propagation path of radio photons through the binary and gives rise to a large spin velocity of the emission region co-rotating with the neutron star. Under favorable conditions these two effects introduce corrections to the conventional time delays at the microsecond level (for a millisecond pulsar in a double neutron star binary with a period of several hours and assuming the emission height of 100 km). Exploiting the dependence of the emission height on frequency (radius-to-frequency mapping) and using multi-frequency observations one should be able to detect these timing corrections even though they are formally degenerate with conventional time delays. Although even in the most accurately timed systems the magnitude of the finite emission height effects is currently somewhat below timing precision, longer-term observations and future facilities like SKA will make measurement of these effects possible, providing an independent check of existing emission height estimates.

It is unclear how very close binary stars form, given that during the pre-main-sequence phase the component stars would have been inside each other. One hypothesis is that they formed further apart, but were brought in closer after formation by gravitational interaction with a third member of the system. If so, all close binaries should be members of triple (or higher-order) systems. As a test of this prediction, we present a search for the signature of third components in archival spectra of close binaries. In our sample of 75 objects, 23 show evidence for the presence of a third component, down to a detection limit of tertiary flux contributions of about 0.8% at 5200 \AA\ (considering only contact and semi-detached binaries, we find 20 out of 66). In a homogeneous subset of 58 contact binaries, we are fairly confident that the 15 tertiaries we have detected are all tertiaries present with mass ratios $0.28\lesssim M_3/M_{12}\lesssim0.75$ and implied outer periods $P\lesssim10^6{\rm d}$. We find that if the frequency of tertiaries was the same as that of binary companions to solar-type stars, we would have expected to detect about 12 tertiaries. In contrast, if all contact binaries were in triple systems, one would expect about 20. Thus, our results are not conclusive, but sufficiently suggestive to warrant further studies.

We consider two types of mechanical coupling between the crust and the core of a precessing neutron star. First, we find that a hydromagnetic (MHD) coupling between the crust and the core strongly modifies the star's precessional modes when $t_a\le\sim (T_s\times T_p)^{1/2}$; here $t_a$ is the Alfven crossing timescale, and $T_s$ and $T_p$ are the star's spin and precession periods, respectively. We argue that in a precessing pulsar PSR B1828-11 the restoring MHD stress prevents a free wobble of the crust relative to the non-precessing core. Instead, the crust and the proton-electron plasma in the core must precess in unison, and their combined ellipticity determines the period of precession. Link has recently shown that the neutron superfluid vortices in the core of PSR B1828-11 cannot be pinned to the plasma; he has also argued that this lack of pinning is expected if the proton Fermi liquid in the core is type-I superconductor. In this case, the neutron superfluid is dynamically decoupled from the precessing motion. The pulsar's precession decays due to the mutual friction between the neutron superfluid and the plasma in the core. The decay is expected to occur over tens to hundreds of precession periods and may be measurable over a human lifetime. Such a measurement would provide information about the strong n-p interaction in the neutron-star core. Second, we consider the effect of gravitomagnetic coupling between the neutron superfluid in the core and the rest of the star and show that this coupling changes the rate of precession by about 10%. The general formalism developed in this paper may be useful for other applications.