Search

We present an adaptation of the exoplanet transit model code batman, in order to permit the generation of X-ray transits. Our underlying extended coronal model assumes an isothermal plasma that is radially symmetric. While this ignores the effect of bright, active regions, observations of transits in X-rays will require averaging across multiple epochs of data for the foreseeable future, significantly reducing the importance of more complex modelling. Our publicly available code successfully generates the predicted W-shaped transit profile in X-rays due to the optically thin nature of the emission, which concentrates the expected observational emission around the limb of the photospheric stellar disc. We provide some examples based on the best known X-ray transit target, HD 189733b, and examine the effect of varying the planet size, coronal temperature, and impact parameter on the resulting transit profile. We also derived scaling relationships for how the overall transit detectability is affected by changing these parameters. Over most of the parameter space, we find that the detectability scales linearly with the cross-sectional area of the planet in X-rays. The relationship with increasing coronal temperature is less fixed, but averages out to a power law with slope $-1/4$ except when the impact parameter is high. Indeed, varying impact parameter has little effect on detectability at all until it approaches unity., Comment: 16 pages, 8 figures, 2 tables. Accepted for publication in the Astronomical Journal. Code available at https://github.com/georgewking/batmanX-rays

One of the leading mechanisms invoked to explain the existence of the radius valley is atmospheric mass loss driven by X-ray and extreme-ultraviolet irradiation, with this process stripping the primordial envelopes of young, small planets to produce the observed bimodal distribution. We present an investigation into the TOI-431 and $\nu^2$ Lupi planetary systems, both of which host planets either side of the radius valley, to determine if their architectures are consistent with evolution by the XUV mechanism. With $\textit{XMM-Newton}$, we measure the current X-ray flux of each star, and see evidence for a stellar flare in the TOI-431 observations. We then simulate the evolution of all of the transiting planets across the two systems in response to the high-energy irradiation over their lifetimes. We use the measured X-ray fluxes as an anchor point for the XUV time evolution in our simulations, and employ several different models of estimating mass loss rates. While the simulations for TOI-431b encountered a problem with the initial calculated radii, we estimate a likely short ($\sim$ Myr) timespan for primordial envelope removal using reasonable assumptions for the initial planet. $\nu^2$ Lupi b is likely harder to strip, but is achieved in a moderate fraction of our simulations. None of our simulations stripped any of the lower density planets of their envelope, in line with prediction. We conclude that both systems are consistent with expectations for generation of the radius valley through XUV photoevaporation., Comment: 14 pages, 12 figures, 5 tables; accepted for publication in MNRAS

We present an analysis of XMM-Newton observations of four stars in the young (670 Myr) open cluster Praesepe. The planets hosted by these stars all lie close in radius-period space to the radius-period valley and/or the Neptunian desert, two features that photoevaporation by X-ray and extreme ultraviolet (EUV) photons could be driving. Although the stars are no longer in the saturated regime, strong X-ray and extreme ultraviolet irradiation is still ongoing. Based on EUV time evolution slopes we derived in a previous paper, in all four cases, two-thirds of their EUV irradiation is still to come. We compare the XMM-Newton light curves to those simultaneously measured with K2 at optical wavelengths, allowing us to search for correlated variability between the X-ray and optical light curves. We find that the X-ray flux decreases and flattens off while the optical flux rises throughout for K2-100, something that could result from active regions disappearing from view as the star spins. Finally, we also investigate possible futures for the four planets in our sample with simulations of their atmosphere evolution still to come, finding that complete photoevaporative stripping of the envelope of three of the four planets is possible, depending on the current planet masses., Comment: 15 pages, 12 figures, 7 tables; accepted for publication in MNRAS

We analyze the two brightest Chandra X-ray flares detected from Sagittarius A*, with peak luminosities more than 600 x and 245 x greater than the quiescent X-ray emission. The brightest flare has a distinctive double-peaked morphology --- it lasts 5.7 ksec ($\sim 2$ hours), with a rapid rise time of 1500 sec and a decay time of 2500 sec. The second flare lasts 3.4 ksec, with rise and decay times of 1700 sec and 1400 sec. These luminous flares are significantly harder than quiescence: the first has a power law spectral index $\Gamma = 2.06\pm 0.14$ and the second has $\Gamma = 2.03\pm 0.27$, compared to $\Gamma = 3.0\pm0.2$ for the quiescent accretion flow. These spectral indices (as well as the flare hardness ratios) are consistent with previously-detected Sgr A* flares, suggesting that bright and faint flares arise from similar physical processes. Leveraging the brightest flare's long duration and high signal-to-noise, we search for intraflare variability and detect excess X-ray power at a frequency of $\nu \approx 3$ mHz, but show that it is an instrumental artifact and not of astrophysical origin. We find no other evidence (at the 95% confidence level) for periodic or quasi-periodic variability in either flares' time series. We also search for non-periodic excess power but do not find compelling evidence in the power spectrum. Bright flares like these remain our most promising avenue for identifying Sgr A*'s short timescale variability in the X-ray, which may probe the characteristic size scale for the X-ray emission region., Comment: Updated to match published version; 19 pages, 7 figures, 3 tables

Finding electromagnetic (EM) counterparts of future gravitational wave (GW) sources would bring rich scientific benefits. A promising possibility, in the case of the coalescence of a super-massive black hole binary (SMBHB), is that prompt emission from merger-induced disturbances in a supersonic circumbinary disk may be detectable. We follow the post-merger evolution of a thin, zero-viscosity circumbinary gas disk with two-dimensional simulations, using the hydrodynamic code FLASH. We analyze perturbations arising from the 530 km/s recoil of a 10^6 M_sun binary, oriented in the plane of the disk, assuming either an adiabatic or a pseudo-isothermal equation of state for the gas. We find that a single-armed spiral shock wave forms and propagates outward, sweeping up about 20% of the mass of the disk. The morphology and evolution of the perturbations agrees well with those of caustics predicted to occur in a collisionless disk. Assuming that the disk radiates nearly instantaneously to maintain a constant temperature, we estimate the amount of dissipation and corresponding post-merger light-curve. The luminosity rises steadily on the time-scale of months, and reaches few times 10^{43} erg/s, corresponding to about 10% of the Eddington luminosity of the central SMBHB. We also analyze the case in which gravitational wave emission results in a 5% mass loss in the merger remnant. The mass-loss reduces the shock overdensities and the overall luminosity of the disk by 15-20%, without any other major effects on the spiral shock pattern., Comment: 16 pages with 14 figures, submitted to MNRAS