New views of the distant stellar halo

Currently, only a small number of Milky Way (MW) stars are known to exist beyond 100 kpc from the Galactic Centre. Though the distribution of these stars in the outer halo is believed to be sparse, they can provide evidence of more recent accretion events than in the inner halo and help map out the MW's dark matter halo to its virial radius. We have re-examined the outermost regions of 11 existing stellar halo models with two synthetic surveys: one mimicking present-day searches for distant M giants and another mimicking RR Lyra (RRL) projections for the Large Synoptic Survey Telescope (LSST). Our models suggest that colour and proper motion cuts currently used to select M giant candidates for follow-up successfully remove nearly all self-contamination from foreground halo dwarf stars and are useful for focusing observations on distant M giants, of which there are thousands to tens of thousands beyond 100 kpc in our models. We likewise expect that LSST will identify comparable numbers of RRLe at these distances. We demonstrate that several observable properties of both tracers, such as proximity of neighbouring stars, proper motions and distances (for RRLe), could help us separate different accreted dwarf galaxies from one another in the distant MW halo. We also discuss prospects for using ratios of M giants to RRLe as a proxy for accretion time, which in the future could provide new constraints on the recent accretion history of our Galaxy

Simulating X-ray Reverberation in the UV-Emitting Regions of Active Galactic Nuclei Accretion Disks with 3D Multi-Frequency Magnetohydrodynamic Simulations

Active galactic nuclei (AGN) light curves observed with different wavebands show that the variability in longer wavelength bands lags the variability in shorter wavelength bands. Measuring these lags, or reverberation mapping, is used to measure the radial temperature profile and extent of AGN disks, typically with a reprocessing model that assumes X-rays are the main driver of the variability in other wavelength bands. To demonstrate how this reprocessing works with realistic accretion disk structures, we use 3D local shearing box multi-frequency radiation magnetohydrodynamic (MHD) simulations to model the UV-emitting region of an AGN disk, which is unstable to the magnetorotational instability (MRI) and convection. At the same time, we inject hard X-rays ($>1$~keV) into the simulation box to study the effects of X-ray irradiation on the local properties of the turbulence and the resulting variability of the emitted UV light curve. We find that disk turbulence is sufficient to drive intrinsic variability in emitted UV light curves and that a damped random walk (DRW) model is a good fit to this UV light curve for timescales $>5$~days. Meanwhile, the injected X-rays have almost no impact on the power spectrum of the emitted UV light curve. In addition, the injected X-ray and emitted UV light curves are only correlated if there is X-ray variability on timescales $>1$~day, in which case we find a correlation coefficient $r=0.52$. These results suggest that hard X-rays with scattering dominated opacity are likely not the main driver of the reverberation signals.Comment: 9 pages, 3 figures, submitted to ApJ

Magnetorotational Instability in a Swirling Partially Ionized Gas

The magnetorotational instability (MRI) has been proposed as the method of angular momentum transport that enables accretion in astrophysical discs. However, for weakly-ionized discs, such as protoplanetary discs, it remains unclear whether the combined non-ideal magnetohydrodynamic (MHD) effects of Ohmic resistivity, ambipolar diffusion, and the Hall effect make these discs MRI-stable. While much effort has been made to simulate non-ideal MHD MRI, these simulations make simplifying assumptions and are not always in agreement with each other. Furthermore, it is difficult to directly observe the MRI astrophysically because it occurs on small scales. Here, we propose the concept of a swirling gas experiment of weakly-ionized argon gas between two concentric cylinders threaded with an axial magnetic field that can be used to study non-ideal MHD MRI. For our proposed experiment, we derive the hydrodynamic equilibrium flow and a dispersion relation for MRI that includes the three non-ideal effects. We solve this dispersion relation numerically for the parameters of our proposed experiment. We find it should be possible to produce non-ideal MRI in such an experiment because of the Hall effect, which increases the MRI growth rate when the vertical magnetic field is anti-aligned with the rotation axis. As a proof of concept, we also present experimental results for a hydrodynamic flow in an unmagnetized prototype. We find that our prototype has a small, but non-negligible, $\alpha$-parameter that could serve as a baseline for comparison to our proposed magnetized experiment, which could be subject to additional turbulence from the MRI.Comment: 14 pages, 13 figures, submitted to MNRA

Negative Lags on the Viscous Timescale in Quasar Photometry and Prospects for Detecting More with LSST

The variability of quasar light curves can be used to study the structure of quasar accretion disks. For example, continuum reverberation mapping uses delays between variability in short and long wavelength bands ("short" lags) to measure the radial extent and temperature profile of the disk. Recently, a potential reverse lag, where variations in shorter wavelength bands lag the longer wavelength bands at the much longer viscous timescale, was detected for Fairall 9. Inspired by this detection, we derive a timescale for these "long" negative lags from fluctuation propagation models and recent simulations. We use this timescale to forecast our ability to detect long lags using the Vera Rubin Legacy Survey of Space and Time (LSST). After exploring several methods, including the interpolated cross-correlation function, a Von-Neumann estimator, javelin, and a maximum-likelihood Fourier method, we find that our two main methods, javelin and the maximum-likelihood method, can together detect long lags of up to several hundred days in mock LSST light curves. Our methods work best on proposed LSST cadences with long season lengths, but can also work for the current baseline LSST cadence, especially if we add observations from other optical telescopes during seasonal gaps. We find that LSST has the potential to detect dozens to hundreds of additional long lags. Detecting these long lags can teach us about the vertical structure of quasar disks and how it scales with different quasar properties.Comment: 40 pages, 36 figures, submitted to Ap

A Negative Long Lag from the Optical to the UV Continuum in Fairall 9

We report the detection of a long-timescale negative lag, where the blue bands lag the red bands, in the nearby Seyfert 1 galaxy Fairall 9. Active Galactic Nuclei (AGN) light curves show variability over a wide range of timescales. By measuring time lags between different wavelengths, the otherwise inaccessible structure and kinematics of the accretion disk can be studied. One common approach, reverberation mapping, quantifies the continuum and line lags moving outwards through the disk at the light-travel time, revealing the size and temperature profile of the disk. Inspired by numerical simulations, we expect longer lags to exist in AGN light curves that travel inward on longer timescales, tracing the accretion process itself. By analyzing AGN light curves in both temporal and frequency space, we report the detection of long-timescale lags ($\sim -70$ days) in Fairall 9 which propagate in the opposite direction to the reverberation lag. The short continuum lag ($<10$ days) is also detected and is consistent with reverberation lags reported in the literature. When fitting the longer lag as a function of frequency with a model motivated by the thin disk model, we find that the disk scale height likely increases outward in the disk. This detection raises the exciting prospect of mapping accretion disk structures across a wide range of AGN parameters.Comment: 24 pages, 14 figures, submitted to Ap

Aligning Retrograde Nuclear Cluster Orbits with an Active Galactic Nucleus Accretion Disc

Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination $i<90^{\circ}$) decrease inclination, as well as semi-major axis ($a$) and eccentricity ($e$) until orbital alignment with the gas disc ("disc capture"). Captured stellar-origin black holes (sBH) add to the embedded AGN population which drives sBH-sBH mergers detectable in gravitational waves using LIGO-Virgo-KAGRA (LVK) or sBH-SMBH mergers detectable with LISA (Laser Interferometer Space Antenna). Captured stars can be tidally disrupted by sBH or the SMBH or rapidly grow into massive 'immortal' stars. Here, we investigate the behaviour of polar and retrograde orbiters ($i \geq 90^{\circ}$) interacting with the disc. We show that retrograde stars are captured faster than prograde stars, flip to prograde ($i<90^{\circ}$) during capture and decrease $a$ dramatically towards the SMBH. For sBH, we find a critical angle $i_{\rm ret} \sim 113^{\circ}$, below which retrograde sBH decay towards embedded prograde orbits ($i \rightarrow 0^{\circ}$), while for $i_{\rm o}>i_{\rm ret}$ sBH decay towards embedded retrograde orbits ($i \rightarrow 180^{\circ}$). sBH near polar orbits ($i \sim 90^{\circ}$) and stars on nearly embedded retrograde orbits ($i \sim 180^{\circ}$) show the greatest decreases in $a$. Whether a star is captured by the disc within an AGN lifetime depends primarily on disc density, and secondarily on stellar type and initial $a$. For sBH, disc capture time is longest for polar orbits, low mass sBH and lower density discs. Larger mass sBH should typically spend more time in AGN discs, with implications for the embedded sBH spin distribution.Comment: 10 pages, 7 figures, 1 table; submitted to MNRA