Search

X-ray binaries (XRBs) exhibit spectral hysteresis for luminosities in the range $10^{-2}\lesssim L/L_\mathrm{Edd}\lesssim 0.3$, with a hard X-ray spectral state that persists from quiescent luminosities up to $\gtrsim 0.3L_\mathrm{Edd}$, transitioning to a soft spectral state that survives with decreasing luminosities down to $\sim 10^{-2}L_\mathrm{Edd}$. We present a possible approach to explain this behavior based on the thermal properties of a magnetically arrested disk simulation. By post-processing the simulation to include radiative effects, we solve for all the thermal equilibrium solutions as the accretion rate, $\dot{M}$, varies along the XRB outburst. For an assumed scaling of the disk scale height and accretion speed with temperature, we find that there exists two solutions in the range of $ 10^{-3}\lesssim\dot{M}/\dot{M}_{\rm Eddington} \lesssim 0.1$ at $r=8\:r_g$ ($ 4\times10^{-2}\lesssim\dot{M}/\dot{M}_{\rm Eddington} \lesssim 0.5$ at $r=3\:r_g$) : a cold, optically thick one and a hot, optically thin one. This opens the possibility of a natural thermal hysteresis in the right range of luminosities for XRBs. We stress that our scenario for the hysteresis does not require to invoke the strong-ADAF principle nor does it require for the magnetization of the disk to change along the XRB outburst. In fact, our scenario requires a highly magnetized disk in the cold, soft state to reproduce the soft-to-hard state transition at the right luminosities. Hence, a prediction of our scenario is that there should be a jet, although possibly very weakly dissipative, in the soft state of XRBs. We also predict that if active galactic nuclei (AGN) have similar hysteresis cycles and are strongly magnetized, they should undergo a soft-to-hard state transition at much lower $L/L_\mathrm{Edd}$ than XRBs., Comment: 15 pages, 10 Figures. Accepted for publication in A&A

We present the largest 3D Particle-in-Cell shearing-box simulations of turbulence driven by the magnetorotational instability, for the first time employing the realistic proton-to-electron mass ratio. We investigate the energy partition between relativistically hot electrons and subrelativistic ions in turbulent accretion flows, a regime relevant to collisionless, radiatively inefficient accretion flows of supermassive black holes, such as those targeted by the Event Horizon Telescope. We provide a simple empirical formula to describe the measured heating ratio between ions and electrons, which can be used for more accurate global modeling of accretion flows with standard fluid approaches such as general-relativistic magnetohydrodynamics.

We present multi-frequency (5-345 GHz) and multi-resolution radio observations of 1ES 1927+654, widely considered one of the most unusual and extreme changing-look active galactic nuclei (CL-AGN). The source was first designated a CL-AGN after an optical outburst in late 2017 and has since displayed considerable changes in X-ray emission, including the destruction and rebuilding of the X-ray corona in 2019-2020. Radio observations prior to 2023 show a faint and compact radio source typical of radio-quiet AGN. Starting in February 2023, 1ES 1927+654 began exhibiting a radio flare with a steep exponential rise, reaching a peak 60 times previous flux levels, and has maintained this higher level of radio emission for over a year to date. The 5-23 GHz spectrum is broadly similar to gigahertz-peaked radio sources, which are understood to be young radio jets less than ~1000 years old. Recent high-resolution VLBA observations at 23.5 GHz now show resolved extensions on either side of the core, with a separation of ~0.15 pc, consistent with a new and mildly relativistic bipolar outflow. A steady increase in the soft X-ray band (0.3-2 keV) concurrent with the radio may be consistent with jet-driven shocked gas, though further observations are needed to test alternate scenarios. This source joins a growing number of CL-AGN and tidal disruption events which show late-time radio activity, years after the initial outburst., Comment: Submitted to ApJ Letters 24 June 2024; Accepted 14 October 2024

JWST observations demonstrate that supermassive black holes (SMBHs) exist by redshifts $z \gtrsim 10$, providing further evidence for "direct collapse" black hole (BH) formation, whereby massive ($\sim 10^{3-5} M_{\odot}$) SMBH seeds are generated within a few Myr as a byproduct of the rapid inflow of gas into the centers of protogalaxies. Here we analyze the intermediate "quasi-star" phase that accompanies some direct collapse models, during which a natal BH accretes mass from and energetically sustains (through accretion) an overlying gaseous envelope. We argue that previous estimates of the maximum BH mass that can be reached during this stage, $\sim 1\%$ of the total quasi-star mass, are unphysical, and arise from underestimating the efficiency with which energy can be transported outward from regions close to the BH. We construct new quasi-star models that consist of an inner, "saturated-convection" region (which conforms to a convection-dominated accretion flow near the BH) matched to an outer, adiabatic envelope. These solutions exist up to a BH mass of $\sim 60\%$ the total quasi-star mass, at which point the adiabatic envelope contains only 2\% of the mass (with the remaining $\sim 38\%$ in the saturated-convection region), and this upper limit is reached within a time of $20-40$ Myr. We conclude that quasi-stars remain a viable route for producing SMBHs at large redshifts, consistent with recent JWST observations., Comment: 14 pages, 6 figures, ApJ Accepted (Updated to match accepted version)

We present an analytic, quasi-local model for accretion discs threaded by net, vertical magnetic flux. In a simple slab geometry and ignoring stochastic mean-field dynamo effects, we calculate the large-scale field resulting from the balance between kinematic field amplification and turbulent diffusion. The ability of the disc to accumulate magnetic flux is sensitive to a single parameter dependent on the ratio of the vertical diffusion time to the Alfv\'en crossing time, and we show how the saturation levels of magnetorotational and other instabilities can govern disc structure and evolution. Under wide-ranging conditions, inflow is governed by large-scale magnetic stresses rather than internal viscous stress. We present models of such "magnetically boosted" discs and show that they lack a radiation pressure-dominated zone. Our model can account for "magnetically elevated" discs as well as instances of midplane outflow and field reversals with height that have been seen in some global simulations. Using the time-dependent features of our model, we find that the incorporation of global transport effects into disc structure can lead to steady or episodic "magnetically arrested discs" (MADs) that maximize the concentration of magnetic flux in their central regions., Comment: 12 pages, no figures, published in Monthly Notices of the Royal Astronomical Society

Double lobe radio sources associated with active galactic nuclei represent one of the longest studied groups in radio astronomy. A particular sub-group of double radio sources comprises the compact symmetric objects (CSOs). CSOs are distinguished by their prominent double structure and sub-kpc total size. It has been argued that the vast majority of high-luminosity CSOs (CSO 2s) represent a distinct class of active galactic nuclei with its own morphological structure and life-cycle. In this work, we present theoretical considerations regarding CSO 2s. We develop a semi-analytic evolutionary model, inspired by the results of large-scale numerical simulations of relativistic jets, that reproduces the features of the radio source population. We show that CSO 2s may be generated by finite energy injections and propose stellar tidal disruption events as a possible cause. We find that tidal disruption events of giant branch stars with masses $\gtrsim1$ M$_\odot$ can fuel these sources and discuss possible approaches to confirming this hypothesis. We predict that if the tidal disruption scenario holds, CSO 2s with sizes less than 400 pc should outnumber larger sources by more than a factor of $10$. Our results motivate future numerical studies to determine whether the scenarios we consider for fueling and source evolution can explain the observed radio morphologies. Multiwavelength observational campaigns directed at these sources will also provide critical insight into the origins of these objects, their environments, and their lifespans., Comment: 9 pages, 4 figures; Accepted for publication in MNRAS

We present 3D fully kinetic shearing-box simulations of pair-plasma magnetorotational turbulence with unprecedented macro-to-microscopic scale separation. While retrieving the expected fluid behavior of the plasma at large scales, we observe a steepening of turbulent spectra at kinetic scales and substantial angular-momentum transport linked with kinetic processes. For the first time, we provide a definitive demonstration of nonthermal particle acceleration in kinetic magnetorotational turbulence agnostically of initial conditions, by means of a novel strategy exploiting synchrotron cooling., Comment: Accepted in PRL

1ES1927+654 is a nearby active galactic nucleus that has shown an enigmatic outburst in optical/UV followed by X-rays, exhibiting strange variability patterns at timescales of months-years. Here we report the unusual X-ray, UV, and radio variability of the source in its post-flare state (Jan 2022- May 2023). Firstly, we detect an increase in the soft X-ray (0.3-2 keV) flux from May 2022- May 2023 by almost a factor of five, which we call the bright-soft-state. The hard X-ray 2-10 keV flux increased by a factor of two, while the UV flux density did not show any significant changes ($\le 30\%$) in the same period. The integrated energy pumped into the soft and hard X-ray during this period of eleven months is $\sim 3.57\times 10^{50}$ erg and $5.9\times 10^{49}$ erg, respectively. From the energetics, it is evident that whatever is producing the soft excess (SE) is pumping out more energy than either the UV or hard X-ray source. Since the energy source presumably is ultimately the accretion of matter onto the SMBH, the SE emitting region must be receiving the majority of this energy. In addition, the source does not follow the typical disc-corona relation found in AGNs, neither in the initial flare (in 2017-2019) nor in the current bright soft state (2022-2023). We found that the core (<1 pc) radio emission at 5 GHz gradually increased till March 2022 but showed a dip in August 2022. The G\"udel Benz relation ($L_{\rm radio}/L_{\rm X-ray}\sim 10^{-5}$), however, is still within the expected range for radio-quiet AGN and further follow-up radio observations are currently being undertaken., Comment: Accepted for publication in ApJ

We propose that the polarization of the emission from the X-ray binary Cygnus X-1, measured using the Imaging X-ray Polarimetry Explorer (IXPE), is imprinted by bulk Comptonization of coronal emission in a mildly relativistic wind or jet with a hollow-cone geometry. Models based on scattering in a static corona overlying a thin accretion disc have difficulty reproducing the relatively high polarization degree (PD ~4%) concurrently with the low inclination (~30 deg) of the binary orbit. We show that bulk outflow with a Lorentz factor > 1.5 is adequate to reproduce the observed PD, with position angle parallel to the large-scale jet, provided that the scattering occurs in a conical sheath offset from the jet axis and our line of sight aligns roughly with the opening angle of the cone. Physically, this flow geometry could represent the entrainment of dense material near the base of an accelerating jet as it passes through the disc corona, or a slow (but still relativistic) sheath around a fast jet. If similar outflows are present in other X-ray binaries at higher inclination, we might expect to see still higher degrees of linear polarization ~10% with an orientation perpendicular to the jet direction., Comment: 4 pages, 3 figures, submitted to MNRAS Letters

Large-scale magnetic fields in the nuclear regions of protogalaxies can promote the formation and early growth of supermassive black holes (SMBHs) by direct collapse and magnetically boosted accretion. Turbulence associated with gravitational infall and star formation can drive the rms field strength toward equipartition with the mean gas kinetic energy; this field has a generic tendency to self-organize into large, coherent structures. If the poloidal component of the field (relative to the rotational axis of a star-forming disc) becomes organized on scales $\lesssim r$ and attains an energy of order a few percent of the turbulent energy in the disc, then dynamo effects are expected to generate magnetic torques capable of increasing the inflow speed and thickening the disc. The accretion flow can transport matter toward the center of mass at a rate adequate to create and grow a massive direct-collapse black hole (DCBH) seed and fuel the subsequent AGN at a high rate, without becoming gravitationally unstable. Fragmentation and star formation are thus suppressed and do not necessarily deplete the mass supply for the accretion flow, in contrast to prevailing models for growing and fueling SMBHs through disc accretion., Comment: 6 pages, no figures, published in Monthly Notices of the Royal Astronomical Society Letters. Moderate revisions include additional references and comparative discussion of angular momentum transport mechanisms

Accretion discs properties should deviate from standard theory when magnetic pressure exceeds the thermal pressure. To quantify these deviations, we present a systematic study of the dynamical properties of magnetically arrested discs (MADs), the most magnetized type of accretion disc. Using an artificial cooling function to regulate the gas temperature, we study MADs of three different thermal thicknesses, $h_\mathrm{th}/r=0.3, 0.1$ and $0.03$. We find that the radial structure of the disc is never mostly supported by the magnetic field. In fact, thin MADs are very near Keplerian. However, as discs gets colder, they become more magnetized and the largest deviations from standard theory appear in our thinnest disc with $h_\mathrm{th}/r=0.03$. In this case, the disc is much more extended vertically and much less dense than in standard theory because of vertical support from the turbulent magnetic pressure and wind-driven angular momentum transport that enhances the inflow speed. The thin disc also dissipates a lot of thermal energy outside of $z/r = \pm 0.03$ and a significant fraction of this dissipation happens in mildly relativistic winds. The enhanced dissipation in low-density regions could possibly feed coronae in X-ray binaries (XRBs) and active galactic nuclei (AGN). Wind-driven accretion will also impact the dynamical evolution of accretion discs and could provide a mechanism to explain the rapid evolution of changing-look AGN and the secular evolution of XRBs. Finally, our MAD winds have terminal velocities and mass loss rates in good agreement with the properties of ultra-fast outflows observed in AGN., Comment: Accepted version. 19 pages, 18 Figures + 4 Appendix Figures

We propose that the strength of angular momentum transport in accretion discs threaded by net vertical magnetic field is determined by a self-regulation mechanism: the magnetorotational instability (MRI) grows until its own turbulent resistivity damps the fastest growing mode on the scale of the disc thickness. Given weak assumptions as to the structure of MRI-derived turbulence, supported by prior simulation evidence, the proposed mechanism reproduces the known scaling of the viscous $\alpha$-parameter, $\alpha \propto \beta_z^{-1/2}$. Here, $\beta_z = 8\pi p_g/B_{z0}^2$ is the initial plasma $\beta$-parameter on the disc midplane, $B_{z0}$ is the net field, and $p_g $ is the midplane gas pressure. We generalize the argument to discs with strong suprathermal toroidal magnetic fields, where the MRI growth rate is modified from the weak-field limit. Additional sources of turbulence are required if such discs are to become magnetically elevated, with the increased scale heights near the midplane that are seen in simulations. We speculate that tearing modes, associated with current sheets broadened by the effective resistivity, are a possible source of enhanced turbulence in elevated discs., Comment: 9 pages, 2 figures, to be published in Monthly Notices of the Royal Astronomical Society. Replaces arXiv:2211.09261, eliminating section 2.2 and with moderate revisions and clarifications in sections 3.1, 3.3, 4.2 and 4.3

The dichotomy between radio-loud (RL) and radio-quiet (RQ) active galactic nuclei (AGN) is thought to be intrinsically related to radio jet production. This difference may be explained by the presence of a strong magnetic field (B-field) that enhances, or is the cause of, the accretion activity and the jet power. Here, we report the first evidence of an intrinsic difference in the dust polarized emission cores of four RL and five RQ obscured AGN using 89 $\mu$m polarization with HAWC+/SOFIA. We find that the thermal polarized emission increases with the nuclear radio-loudness, $R_{20} = L_{\rm 5GHz}/ L_{\rm 20\mu m}$. The dust emission cores of RL AGN are measured to be polarized, $\sim5-11$%, while RQ AGN are unpolarized, $<1$%. For RQ AGN, our results are consistent with the observed region being filled with an unmagnetized or highly turbulent, disk and/or expanding outflow at scales of $5-130$ pc from the AGN. For RL AGN, the measured $89$ $\mu$m polarization arises primarily from magnetically aligned dust grains associated with a $5-130$ pc-scale dusty obscuring structure with a toroidal B-field orientation highly offset, $65\pm22^{\circ}$, with respect to the jet axis. Our results indicate that the size and strength of the B-fields surrounding the AGN are intrinsically related to the strength of the jet power -- the stronger the jet power is, the larger and stronger the toroidal B-field is. The detection of a $\le130$ pc-scale ordered toroidal B-field suggests that a) the infalling gas that fuels RL AGN is magnetized, b) there is a magnetohydrodynamic wind that collimates the jet, and/or c) the jet is able to magnetize its surroundings., Comment: 23 pages, 6 figures, 7 tables. Accepted for publication; ApJ

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e. kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic Particle-in-Cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e. MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron-ion simulations, on which we will focus in future work., Comment: Submitted to ApJ

Relativistic magnetized jets, such as those from AGN, GRBs and XRBs, are susceptible to current- and pressure-driven MHD instabilities that can lead to particle acceleration and non-thermal radiation. Here we investigate the development of these instabilities through 3D kinetic simulations of cylindrically symmetric equilibria involving toroidal magnetic fields with electron-positron pair plasma. Generalizing recent treatments by Alves et al. (2018) and Davelaar et al. (2020), we consider a range of initial structures in which the force due to toroidal magnetic field is balanced by a combination of forces due to axial magnetic field and gas pressure. We argue that the particle energy limit identified by Alves et al. (2018) is due to the finite duration of the fast magnetic dissipation phase. We find a rather minor role of electric fields parallel to the local magnetic fields in particle acceleration. In all investigated cases a kink mode arises in the central core region with a growth timescale consistent with the predictions of linearized MHD models. In the case of a gas-pressure-balanced (Z-pinch) profile, we identify a weak local pinch mode well outside the jet core. We argue that pressure-driven modes are important for relativistic jets, in regions where sufficient gas pressure is produced by other dissipation mechanisms., Comment: 22 pages, 20 figures, accepted for publication in ApJ

The nearby type-II AGN 1ES1927+654 went through a violent changing-look (CL) event beginning December 2017 during which the optical and UV fluxes increased by four magnitudes over a few months, and broad emission lines newly appeared in the optical/UV. By July 2018 the X-ray coronal emission had completely vanished, only to reappear a few months later. In this work we report the evolution of the radio, optical, UV and X-rays from the pre-flare state through mid-2021 with new and archival data from the Very Long Baseline Array (VLBA), the European VLBI Network, the Very Large Array (VLA), the Telescopio Nazionale Galileo (TNG), Gran Telescopio Canarias (GTC), The Neil Gehrels Swift observatory and XMM-Newton. The main results from our work are: (1) The source has returned to its pre-CL state in optical, UV, and X-ray; the disk-corona relation has been re-established as has been in the pre-CL state, with an $\alpha_{\rm OX}\sim 1.02$. The optical spectra are dominated by narrow emission lines. (2) The UV light curve follows a shallower slope of $\propto t^{-0.91\pm 0.04}$ compared to that predicted by a tidal disruption event. We conjecture that a magnetic flux inversion event is the possible cause for this enigmatic event. (3) The compact radio emission which we tracked in the pre-CL (2014), during CL (2018) and post-CL(2021) at spatial scales $<1$ pc was at its lowest level during the changing look event in 2018, nearly contemporaneous with a low $2-10$ keV emission. The radio to X-ray ratio of the compact source $L_{\rm Radio}/L_{\rm X-ray}\sim 10^{-5.5}$, follows the Gudel-Benz relation, typically found in coronally active stars, and several AGN. (4) We do not detect any presence of nascent jets at the spatial scales of $\sim 5-10$ pc., Comment: Resubmitted to ApJ following minor comments from the referee

Astrophysical jets, launched from the immediate vicinity of accreting black holes, carry away large amounts of power in a form of bulk kinetic energy of jet particles and electromagnetic flux. Here we consider a simple analytical model for relativistic jets at larger distances from their launching sites, assuming a cylindrical axisymmetric geometry with a radial velocity shear, and purely toroidal magnetic field. We argue that, as long as the jet plasma is in magnetohydrostatic equilibrium, such outflows tend to be particle dominated, i.e. the ratio of the electromagnetic to particle energy flux, integrated over the jet cross-sectional area, is typically below unity, $\sigma < 1$. At the same time, for particular magnetic and radial velocity profiles, magnetic pressure may still dominate over particle pressure for certain ranges of the jet radius, i.e. the local jet plasma parameter $\beta_{pl} < 1$, and this may be relevant in the context of particle acceleration and production of high-energy emission in such systems. The jet magnetization parameter can be elevated up to the modest values $\sigma \lesssim \mathcal{O}(10)$ only in the case of extreme gradients or discontinuities in the gaseous pressure, and a significantly suppressed velocity shear. Such configurations, which consist of a narrow, unmagnetized jet spine surrounded by an extended, force-free layer, may require an additional poloidal field component to stabilize them against current-driven oscillations, but even this will not elevate substantially their $\sigma$ parameter., Comment: The Astrophysical Journal, in press

Magnetically arrested accretion discs (MADs) around black holes (BH) have the potential to stimulate the production of powerful jets and account for recent ultra-high-resolution observations of BH environments. Their main properties are usually attributed to the accumulation of dynamically significant net magnetic (vertical) flux throughout the arrested region, which is then regulated by interchange instabilities. Here we propose instead that it is mainly a dynamically important {\it toroidal} field -- the result of dynamo action triggered by the significant but still relatively weak vertical field -- that defines and regulates the properties of MADs. We suggest that rapid convection-like instabilities, involving interchange of toroidal flux tubes and operating concurrently with the magnetorotatonal instability (MRI), can regulate the structure of the disc and the escape of net flux. We generalize the convective stability criteria and disc structure equations to include the effects of a strong toroidal field and show that convective flows could be driven towards two distinct marginally stable states, one of which we associate with MADs. We confirm the plausibility of our theoretical model by comparing its quantitative predictions to simulations of both MAD and SANE (strongly magnetized but not "arrested") discs, and suggest a set of criteria that could help to distinguish MADs from other accretion states. Contrary to previous claims in the literature, we argue that MRI is not suppressed in MADs and is probably responsible for the existence of the strong toroidal field., Comment: 12 pages, 11 Figures, Accepted for publication in MNRAS

We present global radiation GRMHD simulations of strongly magnetized accretion onto a spinning, stellar mass black hole at sub-Eddington rates. Using a frequency-dependent Monte Carlo procedure for Compton scattering, we self-consistently evolve a two-temperature description of the ion-electron fluid and its radiation field. For an Eddington ratio $L/L_{\rm Edd} \gtrsim 10^{-3}$, the emergent spectrum forms an apparent power law shape from thermal Comptonization up to a cutoff at $\simeq 100$ keV, characteristic of that seen in the hard spectral states of black hole X-ray binary systems. At these luminosities, the radiative efficiency is high ($\approx 24\%$) and results in a denser midplane region where magnetic fields are dynamically important. For $L/L_{\rm Edd} \sim 10^{-2}$, our hot accretion flow appears to undergo thermal runaway and collapse. Our simulations demonstrate that hot accretion flows can be radiatively efficient and provide an estimate of their maximum luminosity., Comment: 9 pages, 4 figures, ApJL in press

Sgr A* exhibits flares in the near-infrared and X-ray bands, with the luminosity in these bands increasing by factors of $10-100$ for $\approx 60$ minutes. One of the models proposed to explain these flares is synchrotron emission of non-thermal particles accelerated by magnetic reconnection events in the accretion flow. We use the results from PIC simulations of magnetic reconnection to post-process 3D two-temperature GRMHD simulations of a magnetically arrested disc (MAD). We identify current sheets, retrieve their properties, estimate their potential to accelerate non-thermal particles and compute the expected non-thermal synchrotron emission. We find that the flux eruptions of MADs can provide suitable conditions for accelerating non-thermal particles to energies $\gamma_e \lesssim 10^6$ and producing simultaneous X-ray and near-infrared flares. For a suitable choice of current-sheet parameters and a simpified synchrotron cooling prescription, the model can simultaneously reproduce the quiescent and flaring X-ray luminosities as well as the X-ray spectral shape. While the near-infrared flares are mainly due to an increase in the temperature near the black hole during the MAD flux eruptions, the X-ray emission comes from narrow current sheets bordering highly magnetized, low-density regions near the black hole and equatorial current sheets where the flux on the black hole reconnects. As a result, not all infrared flares are accompanied by X-ray ones. The non-thermal flaring emission can extend to very hard ($\lesssim 100$ keV) X-ray energies., Comment: 12 pages, 8 Figures, Accepted for publication in MNRAS

We argue that the changing-look event in the active galactic nucleus 1ES 1927+654, followed by a dip of 3 orders of magnitude in the X-ray luminosity, is controlled by a change in the accretion rate and an inversion of magnetic flux in a magnetically arrested disk (MAD). Before the changing-look event, strong magnetic flux on the black hole powers X-ray emission via the Blandford-Znajek process while the UV emission is produced by a radiatively inefficient magnetized disk. An advection event, bringing flux of the opposite polarity, propagates inward leading, first, to a rise in the UV/optical luminosity and, then, to a dip in the X-ray luminosity when it reaches the black hole. We estimate the timescale for magnetic flux advection and find that the observed timescale between the beginning of the changing-look event and the minimum in the X-ray luminosity, $\approx200$ days, is in agreement with the time needed to cancel the magnetic flux in a MAD extending to $\approx180\:r_g$. Although flux inversion events might be rare due to the large ratio of flux-to-mass that is needed, we argue that AGN showing an unusually high ratio of X-ray to UV luminosity are prime candidates for such events. We also suggest that similar events may lead to jet interruptions in radio-loud objects., Comment: Submitted to MNRAS Letters. 5 pages, 1 Figure

Dwarf nov\ae\ (DNe) and low mass X-ray binaries (LMXBs) are compact binaries showing variability on time scales from years to less than seconds. Here, we focus on explaining part of the rapid fluctuations in DNe, following the framework of recent studies on the monthly eruptions of DNe that use a hybrid disk composed of an outer standard disk and an inner magnetized disk. We show that the ionization instability, that is responsible for the monthly eruptions of DNe, is also able to operate in the inner magnetized disk. Given the low density and the fast accretion time scale of the inner magnetized disk, the ionization instability generates small, rapid heating and cooling fronts propagating back and forth in the inner disk. This leads to quasi-periodic oscillations (QPOs) with a period of the order of $1000$ s. A strong prediction of our model is that these QPOs can only develop in quiescence or at the beginning/end of an outburst. We propose that these rapid fluctuations might explain a subclass of already observed QPOs in DNe as well as a, still to observe, subclass of QPOs in LMXBs. We also extrapolate to the possibility that the radiation pressure instability might be related to Type B QPOs in LMXBs., Comment: Accepted for publication in MNRAS. 12 pages, 6 figures, 2 Tables

Relativistic jets, or highly collimated and fast-moving outflows, are endemic to many astrophysical phenomena. The jets produced by gamma-ray bursts and tidal disruption events are accompanied by the accretion of material onto a black hole or neutron star, with the accretion rate exceeding the Eddington limit of the compact object by orders of magnitude. In such systems, radiation dominates the energy-momentum budget of the outflow, and the dynamical evolution of the jet is governed by the equations of radiation hydrodynamics. Here we show that there are analytic solutions to the equations of radiation hydrodynamics in the viscous (i.e., diffusive) regime that describe structured, relativistic jets, which consist of a fast-moving, highly relativistic core surrounded by a slower-moving, less relativistic sheath. In these solutions, the slower-moving, outer sheath contains most of the mass, and the jet structure is mediated by local anisotropies in the radiation field. We show that, depending on the pressure and density profile of the ambient medium, the angular profile of the jet Lorentz factor is Gaussian or falls off even more steeply with angle. These solutions have implications for the nature of jet production and evolution in hyperaccreting systems, and demonstrate that such jets -- and the corresponding jet structure -- can be sustained entirely by radiative processes. We discuss the implications of these findings in the context of jetted tidal disruption events and short and long gamma-ray bursts., Comment: MNRAS Accepted

High-energy astrophysical systems frequently contain collisionless relativistic plasmas that are heated by turbulent cascades and cooled by emission of radiation. Understanding the nature of this radiative turbulence is a frontier of extreme plasma astrophysics. In this paper, we use particle-in-cell simulations to study the effects of external inverse Compton radiation on turbulence driven in an optically thin, relativistic pair plasma. We focus on the statistical steady state (where injected energy is balanced by radiated energy) and perform a parameter scan spanning from low magnetization to high magnetization ($0.04 \lesssim \sigma \lesssim 11$). We demonstrate that the global particle energy distributions are quasi-thermal in all simulations, with only a modest population of nonthermal energetic particles (extending the tail by a factor of $\sim 2$). This indicates that nonthermal particle acceleration (observed in similar non-radiative simulations) is quenched by strong radiative cooling. The quasi-thermal energy distributions are well fit by analytic models in which stochastic particle acceleration (due to, e.g., second-order Fermi mechanism or gyroresonant interactions) is balanced by the radiation reaction force. Despite the efficient thermalization of the plasma, nonthermal energetic particles do make a conspicuous appearance in the anisotropy of the global momentum distribution as highly variable, intermittent beams (for high magnetization cases). The beamed high-energy particles are spatially coincident with intermittent current sheets, suggesting that localized magnetic reconnection may be a mechanism for kinetic beaming. This beaming phenomenon may explain rapid flares observed in various astrophysical systems (such as blazar jets, the Crab nebula, and Sagittarius A*)., Comment: 26 pages, 30 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

We use global magnetohydrodynamic simulations to study the influence of net vertical magnetic fields on the structure of geometrically thin ($H/r \approx 0.05$) accretion disks in the Newtonian limit. We consider initial mid-plane gas to magnetic pressure ratios $\beta_0 = 1000,\, 300$ and $100$, spanning the transition between weakly and strongly magnetized accretion regimes. We find that magnetic pressure is important for the disks' vertical structure in all three cases, with accretion occurring at $z/R\approx 0.2$ in the two most strongly magnetized models. The disk midplane shows outflow rather than accretion. Accretion through the surface layers is driven mainly by stress due to coherent large scale magnetic field rather than by turbulent stress. Equivalent viscosity parameters measured from our simulations show similar dependencies on initial $\beta_0$ to those seen in shearing box simulations, though the disk midplane is not magnetic pressure dominated even for the strongest magnetic field case. Winds are present but are not the dominant driver of disk evolution. Over the (limited) duration of our simulations, we find evidence that the net flux attains a quasi-steady state at levels that can stably maintain a strongly magnetized disk. We suggest that geometrically thin accretion disks in observed systems may commonly exist in a magnetically "elevated" state, characterized by non-zero but modest vertical magnetic fluxes, with potentially important implications for disk phenomenology in X-ray binaries (XRBs) and active galactic nuclei (AGN).

Nonthermal relativistic plasmas are ubiquitous in astrophysical systems like pulsar wind nebulae and active galactic nuclei, as inferred from their emission spectra. The underlying nonthermal particle acceleration (NTPA) processes have traditionally been modeled with a Fokker-Planck (FP) diffusion-advection equation in momentum space. In this paper, we directly test the FP framework in ab-initio kinetic simulations of driven magnetized turbulence in relativistic pair plasma. By statistically analyzing the motion of tracked particles, we demonstrate the diffusive nature of NTPA and measure the FP energy diffusion ($D$) and advection ($A$) coefficients as functions of particle energy $\gamma m_e c^2$. We find that $D(\gamma)$ scales as $\gamma^2$ in the high-energy nonthermal tail, in line with 2nd-order Fermi acceleration theory, but has a much weaker scaling at lower energies. We also find that $A$ is not negligible and reduces NTPA by tending to pull particles towards the peak of the particle energy distribution. This study provides strong support for the FP picture of turbulent NTPA, thereby enhancing our understanding of space and astrophysical plasmas., Comment: 7 pages, 5 figures, Accepted to Astrophysical Journal Letters

Electron and ion energization (i.e., heating and nonthermal acceleration) is a fundamental, but poorly understood, outcome of plasma turbulence. In this work, we present new results on this topic from particle-in-cell simulations of driven turbulence in collisionless, relativistic electron-ion plasma. We focus on temperatures such that ions (protons) are sub-relativistic and electrons are ultra-relativistic, a regime relevant for high-energy astrophysical systems such as hot accretion flows onto black holes. We find that ions tend to be preferentially heated, gaining up to an order of magnitude more energy than electrons, and propose a simple empirical formula to describe the electron-ion energy partition as a function of the ratio of electron-to-ion gyroradii (which in turn is a function of initial temperatures and plasma beta). We also find that while efficient nonthermal particle acceleration occurs for both species in the ultra-relativistic regime, nonthermal electron populations are diminished with decreasing temperature whereas nonthermal ion populations are essentially unchanged. These results have implications for modeling and interpreting observations of hot accretion flows., Comment: 6 pages, 4 figures, accepted for publication in Physical Review Letters

We carry out an extensive linear stability analysis of magnetized cylindrical jets in a global framework. Foregoing the commonly invoked force-free limit, we focus on the small-scale, internal instabilities triggered in regions of the jet dominated by a toroidal magnetic field, with a weak vertical field and finite thermal pressure gradient. Such regions are likely to occur far from the jet source and boundaries, and are potential sites of magnetic energy dissipation that is essential to explain the particle acceleration and radiation observed from astrophysical jets. We validate the local stability analysis of Begelman by verifying that the eigenfunctions of the most unstable modes are radially localized. This finding allows us to propose a generic stability criterion in the presence of a weak vertical field. A stronger vertical field with a radial gradient complicates the stability criterion, due to the competition between the destabilizing thermal pressure gradient and stabilizing magnetic pressure gradients. Nevertheless, we argue that the jet interiors generically should be subject to rapidly growing, small-scale instabilities, capable of producing current sheets that lead to dissipation. We identify some new instabilities, not predicted by the local analysis, which are sensitive to the background radial profiles but have smaller growth rates than the local instabilities, and discuss the relevance of our work to the findings of recent numerical jet simulations., Comment: 31 pages, 15 figures and 6 tables. Revised version, accepted for publication in MNRAS

Rapid, large amplitude variability at optical to X-ray wavelengths is now seen in an increasing number of Seyfert galaxies and luminous quasars. The variations imply a global change in accretion power, but are too rapid to be communicated by inflow through a standard thin accretion disc. Such discs are long known to have difficulty explaining the observed optical/UV emission from active galactic nuclei. Here we show that alternative models developed to explain these observations have larger scale heights and shorter inflow times. Accretion discs supported by magnetic pressure in particular are geometrically thick at all luminosities, with inflow times as short as the observed few year timescales in extreme variability events to date. Future time-resolved, multi-wavelength observations can distinguish between inflow through a geometrically thick disc as proposed here, and alternative scenarios of extreme reprocessing of a central source or instability-driven limit cycles., Comment: 5 pages, 2 figures, submitted to MNRAS letters

We apply collisionless particle-in-cell simulations of relativistic pair plasmas to explore whether driven turbulence is a viable high-energy astrophysical particle accelerator. We characterize nonthermal particle distributions for varying system sizes up to $L/2\pi\rho_{e0} = 163$, where $L/2\pi$ is the driving scale and $\rho_{e0}$ is the initial characteristic Larmor radius. We show that turbulent particle acceleration produces power-law energy distributions that, when compared at a fixed number of large-scale dynamical times, slowly steepen with increasing system size. We demonstrate, however, that convergence is obtained by comparing the distributions at different times that increase with system size (approximately logarithmically). We suggest that the system-size dependence arises from the time required for particles to reach the highest accessible energies via Fermi acceleration. The converged power-law index of the energy distribution, $\alpha \approx 3.0$ for magnetization $\sigma = 3/8$, makes turbulence a possible explanation for nonthermal spectra observed in systems such as the Crab nebula., Comment: 7 pages, 4 figures, submitted to Astrophysical Journal Letters

We have modeled direct collapse of a primordial gas within dark matter halos in the presence of radiative transfer, in high-resolution zoom-in simulations in a cosmological framework, down to the formation of the photosphere and the central object. Radiative transfer has been implemented in the flux-limited diffusion (FLD) approximation. Adiabatic models were run for comparison. We find that (a) the FLD flow forms an irregular central structure and does not exhibit fragmentation, contrary to adiabatic flow which forms a thick disk, driving a pair of spiral shocks, subject to Kelvin-Helmholtz shear instability forming fragments; (b) the growing central core in the FLD flow quickly reaches ~10 Mo and a highly variable luminosity of 10^{38}-10^{39} erg/s, comparable to the Eddington luminosity. It experiences massive recurrent outflows driven by radiation force and thermal pressure gradients, which mix with the accretion flow and transfer the angular momentum outwards; and (c) the interplay between these processes and a massive accretion, results in photosphere at ~10 AU. We conclude that in the FLD model (1) the central object exhibits dynamically insignificant rotation and slower than adiabatic temperature rise with density; (2) does not experience fragmentation leading to star formation, thus promoting the fast track formation of a supermassive black hole (SMBH) seed; (3) inclusion of radiation force leads to outflows, resulting in the mass accumulation within the central 10^{-3} pc, which is ~100 times larger than characteristic scale of star formation. The inclusion of radiative transfer reveals complex early stages of formation and growth of the central structure in the direct collapse scenario of SMBH seed formation., Comment: 19 pages, 16 figures, MNRAS, accepted for publication

We examine the properties of strongly magnetized accretion discs in a global framework, with particular focus on the evolution of magnetohydrodynamic instabilities such as the magnetorotational instability (MRI). Work by Pessah and Psaltis showed that MRI is stabilized beyond a critical toroidal field in compressible, differentially rotating flows and, also, reported the appearance of two new instabilities beyond this field. Their results stemmed from considering geometric curvature effects due to the suprathermal background toroidal field, which had been previously ignored in weak-field studies. However, their calculations were performed under the local approximation, which poses the danger of introducing spurious behavior due to the introduction of global geometric terms in an otherwise local framework. In order to avoid this, we perform a global eigenvalue analysis of the linearized MHD equations in cylindrical geometry. We confirm that MRI indeed tends to be highly suppressed when the background toroidal field attains the Pessah-Psaltis limit. We also observe the appearance of two new instabilities that emerge in the presence of highly suprathermal toroidal fields. These results were additionally verified using numerical simulations in PLUTO. There are, however, certain differences between the the local and global results, especially in the vertical wavenumber occupancies of the various instabilities, which we discuss in detail. We also study the global eigenfunctions of the most unstable modes in the suprathermal regime, which are inaccessible in the local analysis. Overall, our findings emphasize the necessity of a global treatment for accurately modeling strongly magnetized accretion discs., Comment: 23 pages, 7 Figures and 1 Table. Accepted for publication in MNRAS

We describe results from particle-in-cell simulations of driven turbulence in collisionless, magnetized, relativistic pair plasma. This physical regime provides a simple setting for investigating the basic properties of kinetic turbulence and is relevant for high-energy astrophysical systems such as pulsar wind nebulae and astrophysical jets. In this paper, we investigate the statistics of turbulent fluctuations in simulations on lattices of up to $1024^3$ cells and containing up to $2 \times 10^{11}$ particles. Due to the absence of a cooling mechanism in our simulations, turbulent energy dissipation reduces the magnetization parameter to order unity within a few dynamical times, causing turbulent motions to become sub-relativistic. In the developed stage, our results agree with predictions from magnetohydrodynamic turbulence phenomenology at inertial-range scales, including a power-law magnetic energy spectrum with index near $-5/3$, scale-dependent anisotropy of fluctuations described by critical balance, log-normal distributions for particle density and internal energy density (related by a $4/3$ adiabatic index, as predicted for an ultra-relativistic ideal gas), and the presence of intermittency. We also present possible signatures of a kinetic cascade by measuring power-law spectra for the magnetic, electric, and density fluctuations at sub-Larmor scales., Comment: 24 pages, 26 figures, submitted for publication

We propose that the accretion disks fueling active galactic nuclei are supported vertically against gravity by a strong toroidal ($\phi-$direction) magnetic field that develops naturally as the result of an accretion disk dynamo. The magnetic pressure elevates most of the gas carrying the accretion flow at $R$ to large heights $z > 0.1 R$ and low densities, while leaving a thin dense layer containing most of the mass --- but contributing very little accretion --- around the equator. We show that such a disk model leads naturally to the formation of a broad emission line region through thermal instability. Extrapolating to larger radii, we demonstrate that local gravitational instability and associated star formation are strongly suppressed compared to standard disk models for AGN, although star formation in the equatorial zone is predicted for sufficiently high mass supply rates. This new class of accretion disk models thus appears capable of resolving two longstanding puzzles in the theory of AGN fueling: the formation of broad emission line regions and the suppression of fragmentation thought to inhibit accretion at the required rates. We show that the disk of stars that formed in the Galactic Center a few million years ago could have resulted from an episode of magnetically elevated accretion at $> 0.1$ of the Eddington limit., Comment: 8 pages, no figures, to appear in Monthly Notices of the Royal Astronomical Society

The rate at which matter flows into a galactic nucleus during early phases of galaxy evolution can sometimes exceed the Eddington limit of the growing central black hole by several orders of magnitude. We discuss the necessary conditions for the black hole to actually accrete this matter at such a high rate, and consider the observational appearance and detectability of a hyperaccreting black hole. In order to be accreted at a hyper-Eddington rate, the infalling gas must have a sufficiently low angular momentum. Although most of the gas is accreted, a significant fraction accumulates in an optically thick envelope with luminosity $\sim L_{\rm Edd}$, probably pierced by jets of much higher power. If $\dot M > 10^3 \dot M_{\rm Edd}$, the envelope spectrum resembles a blackbody with a temperature of a few thousand K, but for lower (but still hyper-Eddington) accretion rates the spectrum becomes a very dilute and hard Wien spectrum. We consider the likelihood of various regimes of hyperaccretion, and discuss its possible observational signatures., Comment: 6 pages, 1 figure, to appear in Monthly Notices of the Royal Astronomical Society

We present results from particle-in-cell simulations of driven turbulence in magnetized, collisionless, and relativistic pair plasma. We find that fluctuations are consistent with the classical $k_\perp^{-5/3}$ magnetic energy spectrum at fluid scales and a steeper $k_\perp^{-4}$ spectrum at sub-Larmor scales, where $k_\perp$ is the wavevector perpendicular to the mean field. We demonstrate the development of a non-thermal, power-law particle energy distribution, $f(E) \sim E^{-\alpha}$, with index $\alpha$ that decreases with increasing magnetization and increases with increasing system size (relative to the characteristic Larmor radius). Our simulations indicate that turbulence can be a viable source of energetic particles in high-energy astrophysical systems, such as pulsar wind nebulae, if scalings asymptotically become insensitive to the system size., Comment: 6 pages, 4 figures, to appear in PRL

We investigate the pre-disruption gravitational dynamics and post-disruption hydrodynamics of the tidal disruption of stars by supermassive black hole (SMBH) binaries. We focus on binaries with relatively low mass primaries ($10^6M_{\odot}$), moderate mass ratios, and separations with reasonably long gravitational wave inspiral times (tens of Myr). First, we generate a large ensemble (between 1 and 10 million) of restricted three-body integrations to quantify the statistical properties of tidal disruptions by circular SMBH binaries of initially-unbound stars. Compared to the reference case of a disruption by a single SMBH, the binary potential induces significant variance into the specific energy and angular momentum of the star at the point of disruption. Second, we use Newtonian numerical hydrodynamics to study the detailed evolution of the fallback debris from 120 disruptions randomly selected from the three-body ensemble (excluding only the most deeply penetrating encounters). We find that the overall morphology of the debris is greatly altered by the presence of the second black hole, and the accretion rate histories display a wide range of behaviors, including order of magnitude dips and excesses relative to control simulations that include only one black hole. Complex evolution persists, in some cases, for many orbital periods of the binary. We find evidence for power in the accretion curves on timescales related to the binary orbital period, though there is no exact periodicity. We discuss our results in the context of future wide-field surveys, and comment on the prospects of identifying and characterizing the subset of events occurring in nuclei with binary SMBHs., Comment: 31 pages, 25 figures, MNRAS Accepted; updated to incorporate referee comments. Movies of simulations can be found here: http://w.astro.berkeley.edu/~eric_coughlin/movies.html

We present a method for analyzing the interaction between radiation and matter in regions of intense, relativistic shear that can arise in many astrophysical situations. We show that there is a simple velocity profile that should be manifested in regions of large shear that have "lost memory" of their boundary conditions, and we use this self-similar velocity profile to construct the surface of last scattering, or $\tau \simeq 1$ surface, as viewed from any comoving point within the flow. We demonstrate that a simple treatment of scattering from this $\tau \simeq 1$ surface exactly conserves photon number, and derive the rate at which the radiation field is heated due to the shear present in the flow. The components of the comoving radiation energy-momentum tensor are calculated, and we show that they have relatively simple, approximate forms that interpolate between the viscous (small shear) and streaming (large shear) limits. We put our expression for the energy-momentum tensor in a covariant form that does not depend on the explicit velocity profile within the fluid and, therefore, represents a natural means for analyzing general, radiation-dominated, relativistic shear flows., Comment: 15 pages, 5 figures. Accepted for publication in ApJ

The power mechanism and accretion geometry for low-power FR1 radio galaxies is poorly understood in comparison to Seyfert galaxies and QSOs. In this paper we use the diagnostic power of the Lya recombination line observed using the Cosmic Origins Spectrograph (COS) aboard the Hubble Space Telescope (HST) to investigate the accretion flows in three well-known, nearby FR1s: M87, NGC4696, and HydraA. The Ly$\alpha$ emission line's luminosity, velocity structure and the limited knowledge of its spatial extent provided by COS are used to assess conditions within a few parsecs of the super-massive black hole (SMBH) in these radio-mode AGN. We observe strong Ly$\alpha$ emission in all three objects with similar total luminosity to that seen in BL Lacertae objects. M87 shows a complicated emission line profile in Lya which varies spatially across the COS aperture and possibly temporally over several epochs of observation. In both NGC 4696 and M87, the Ly$\alpha$ luminosities $\sim 10^{40}$ ergs/s are closely consistent with the observed strength of the ionizing continuum in Case B recombination theory and with the assumption of near unity covering factor. It is possible that the Ly$\alpha$ emitting clouds are ionized largely by beamed radiation associated with the jets. Long-slit UV spectroscopy can be used to test this hypothesis. Hydra A and the several BL Lac objects studied in this and previous papers have Lya luminosities larger than M87 but their extrapolated, non-thermal continua are so luminous that they over-predict the observed strength of Ly$\alpha$, a clear indicator of relativistic beaming in our direction. Given their substantial space density ($\sim 4\times10^{-3} Mpc^{-3}$) the unbeamed Lyman continuum radiation of FR1s may make a substantial minority contribution (~10%) to the local UV background if all FR1s are similar to M87 in ionizing flux level., Comment: Eight figures, 15 pages in ApJ format. ApJ accepted

A tidal disruption event (TDE) -- when a star is destroyed by the immense gravitational field of a supermassive black hole -- transforms a star into a stream of tidally-shredded debris. The properties of this debris ultimately determine the observable signatures of TDEs. Here we derive a simple, self-similar solution for the velocity profile of the debris streams produced from TDEs, and show that this solution agrees extremely well with numerical results. Using this self-similar solution, we calculate an analytic, approximate expression for the radial density profile of the stream. We show that there is a critical adiabatic index that varies as a function of position along the stream above (below) which the stream is unstable (stable) to gravitational fragmentation. We also calculate the impact of heating and cooling on this stability criterion., Comment: 17 pages, 10 figures, Resubmitted to MNRAS after first referee report

Motivated by indirect observational evidence for strongly magnetized accretion discs around black holes, and the novel theoretical properties of such solutions, we investigate how a strong magnetization state can develop and persist. To this end, we perform local simulations of accretion discs with an initially purely toroidal magnetic field of equipartition strength. We demonstrate that discs with zero net vertical magnetic flux and realistic boundary conditions cannot sustain a strong toroidal field. However, a magnetic pressure-dominated disc can form from an initial configuration with a sufficient amount of net vertical flux and realistic boundary conditions. Our results suggest that poloidal flux is a necessary prerequisite for the sustainability of strongly magnetized accretion discs., Comment: 7 pages, 3 figures, 1 table; Accepted to MNRAS

We have developed a quasi-analytical model for the production of radiation in strong-line blazars, assuming a spine-sheath jet structure. The model allows us to study how the spine and sheath spectral components depend on parameters describing the geometrical and physical structure of "the blazar zone". We show that typical broad-band spectra of strong-line blazars can be reproduced by assuming the magnetization parameter to be of order unity and reconnection to be the dominant dissipation mechanism. Furthermore, we demonstrate that the spine-sheath model can explain why gamma-ray variations are often observed to have much larger amplitudes than the corresponding optical variations. The model is also less demanding of jet power than one-zone models, and can reproduce the basic features of extreme gamma-ray events., Comment: 7 pages, 9 figures. Minor changes are made to match the version accepted for publication in MNRAS

Strongly magnetized accretion discs around black holes have attractive features that may explain enigmatic aspects of X-ray binary behaviour. The structure and evolution of these discs are governed by a dynamo-like mechanism, which channels part of the accretion power liberated by the magnetorotational instability (MRI) into an ordered toroidal magnetic field. To study dynamo activity, we performed three-dimensional, stratified, isothermal, ideal magnetohydrodynamic shearing box simulations. The strength of the self-sustained toroidal magnetic field depends on the net vertical magnetic flux, which we vary across almost the entire range over which the MRI is linearly unstable. We quantify disc structure and dynamo properties as a function of the initial ratio of mid-plane gas pressure to vertical magnetic field pressure, $\beta_0^{\rm mid} = p_{\rm gas} / p_B$. For $10^5 \geq \beta_0^{\rm mid} \geq 10$ the effective $\alpha$-viscosity parameter scales as a power-law. Dynamo activity persists up to and including $\beta_0^{\rm mid} = 10^2$, at which point the entire vertical column of the disc is magnetic pressure-dominated. Still stronger fields result in a highly inhomogeneous disc structure, with large density fluctuations. We show that the turbulent steady state $\beta^{\rm mid}$ in our simulations is well-matched by the analytic model of Begelman et al. (2015) describing the creation and buoyant escape of toroidal field, while the vertical structure of the disc can be broadly reproduced using this model. Finally, we discuss the implications of our results for observed properties of X-ray binaries., Comment: 19 pages, 13 figures, 3 tables; Accepted to MNRAS

A tidal disruption event, which occurs when a star is destroyed by the gravitational field of a supermassive black hole, produces a stream of debris, the evolution of which ultimately determines the observational properties of the event. Here we show that a post-periapsis caustic -- a location where the locus of gas parcels comprising the stream would collapse into a two-dimensional surface if they evolved solely in the gravitational field of the hole -- occurs when the pericenter distance of the star is on the order of the tidal radius of the hole. It is demonstrated that this "pancake" induces significant density perturbations in the debris stream, and, for stiffer equations of state (adiabatic index $\gamma \gtrsim 5/3$), these fluctuations are sufficient to gravitationally destabilize the stream, resulting in its fragmentation into bound clumps. The results of our findings are discussed in the context of the observational properties of tidal disruption events., Comment: 18 pages, 14 figures. Accepted for publication in MNRAS

We use cosmological adaptive mesh refinement (AMR) code Enzo zoom-in simulations to study the long term evolution of the collapsing gas within dark matter (DM) halos at high redshifts. This direct collapse process is a leading candidate for rapid formation of supermassive black hole (SMBH) seeds at high z. To circumvent the Courant condition at small radii, we have used the sink particle method, and focus on the evolution on scales ~0.01-10 pc. The collapse proceeds in two stages, with the secondary runaway happening within the central 10 pc, and with no detected fragmentation. The sink particles form when the collapsing gas requires additional refinement of the grid size at the highest refinement level. Their mass never exceeds ~10^3 Mo, with the sole exception of the central seed which grows dramatically to ~ 2 x 10^6 Mo in ~2 Myr, confirming the feasibility of this path to the SMBH. The time variability of angular momentum axis in the accreted gas results in the formation of two misaligned disks --- a small inner disk, and a more massive, outer disk which is inclined by ~45^o to the inner disk. The self-gravity of these disks is heavily diluted --- both disks lie within the Roche limit of the central seed. While the inner disk is geometrically thin and weakly asymmetric, the outer disk flares due to turbulent motions as a result of the massive inflow along a pair of penetrating filaments. The geometry of inflow via filaments determines the dominant and secondary Fourier modes in this disk --- these modes have a non-self-gravitational origin. We do not confirm that m=1 is a principal mode that drives the inflow in the presence of a central massive object. While the positions of the disks depend on the scale chosen to break the self-similar collapse, the overall configuration appears to be generic, and is expected to form when the central seed becomes sufficiently massive., Comment: 14 pages, 11 figures, MNRAS, in press, typos corrected

We investigate the distribution of particle acceleration sites, independently of the actual acceleration mechanism, during plasmoid-dominated, relativistic collisionless magnetic reconnection by analyzing the results of a particle-in-cell numerical simulation. The simulation is initiated with Harris-type current layers in pair plasma with no guide magnetic field, negligible radiative losses, no initial perturbation, and using periodic boundary conditions. We find that the plasmoids develop a robust internal structure, with colder dense cores and hotter outer shells, that is recovered after each plasmoid merger on a dynamical time scale. We use spacetime diagrams of the reconnection layers to probe the evolution of plasmoids, and in this context we investigate the individual particle histories for a representative sample of energetic electrons. We distinguish three classes of particle acceleration sites associated with (1) magnetic X-points, (2) regions between merging plasmoids, and (3) the trailing edges of accelerating plasmoids. We evaluate the contribution of each class of acceleration sites to the final energy distribution of energetic electrons -- magnetic X-points dominate at moderate energies, and the regions between merging plasmoids dominate at higher energies. We also identify the dominant acceleration scenarios, in order of decreasing importance -- (1) single acceleration between merging plasmoids, (2) single acceleration at a magnetic X-point, and (3) acceleration at a magnetic X-point followed by acceleration in a plasmoid. Particle acceleration is absent only in the vicinity of stationary plasmoids. The effect of magnetic mirrors due to plasmoid contraction does not appear to be significant in relativistic reconnection., Comment: 14 pages, 6 figures, 2 tables, accepted in ApJ

Magnetohydrodynamic accretion disk simulations suggest that much of the energy liberated by the magnetorotational instability (MRI) can be channeled into large-scale toroidal magnetic fields through dynamo action. Under certain conditions, this field can dominate over gas and radiation pressure in providing vertical support against gravity, even close to the midplane. Using a simple model for the creation of this field, its buoyant rise, and its coupling to the gas, we show how disks could be driven into this magnetically dominated state and deduce the resulting vertical pressure and density profiles. Applying an established criterion for MRI to operate in the presence of a toroidal field, we show that magnetically supported disks can have two distinct MRI-active regions, separated by a "dead zone" where local MRI is suppressed, but where magnetic energy continues to flow upward from the dynamo region below. We suggest that the relative strengths of the MRI zones, and the local poloidal flux, determine the spectral states of X-ray binaries. Specifically, "intermediate" and "hard" accretion states occur when MRI is triggered in the hot, upper zone of the corona, while disks in "soft" states do not develop the upper MRI zone. We discuss the conditions under which various transitions should take place and speculate on the relationship of dynamo activity to the various types of quasi-periodic oscillations that sometimes appear in the hard spectral components. The model also explains why luminous accretion disks in the "soft" state show no signs of the thermal/viscous instability predicted by standard alpha models., Comment: ApJ, in press. Minor corrections to Figure 1 and section 4.4 to match proof version

We analyze the interaction of a radiation-dominated jet and its surroundings using the equations of radiation hydrodynamics in the viscous limit. In a previous paper we considered the two-stream scenario, which treats the jet and its surroundings as distinct media interacting through radiation viscous forces. Here we present an alternative boundary layer model, known as the free-streaming jet model -- where a narrow stream of fluid is injected into a static medium -- and present solutions where the flow is ultrarelativistic and the boundary layer is dominated by radiation. It is shown that these jets entrain material from their surroundings and that their cores have a lower density of scatterers and a harder spectrum of photons, leading to observational consequences for lines of sight that look "down the barrel of the jet." These jetted outflow models may be applicable to the jets produced during long gamma-ray bursts and super-Eddington phases of tidal disruption events., Comment: 10 pages, 10 figures. Accepted for publication in ApJ

Using the relativistic equations of radiation hydrodynamics in the viscous limit, we analyze the boundary layers that develop between radiation-dominated jets and their environments. In this paper we present the solution for the self-similar, 2-D, plane-parallel two-stream problem, wherein the jet and the ambient medium are considered to be separate, interacting fluids, and we compare our results to those of previous authors. (In a companion paper we investigate an alternative scenario, known as the free-streaming jet model.) Consistent with past findings, we show that the boundary layer that develops between the jet and its surroundings creates a region of low-density material. These models may be applicable to sources such as super-Eddington tidal disruption events and long gamma-ray bursts., Comment: 9 pages, 6 figures, paper II to appear tomorrow. Accepted for publication in ApJ

We study the early stage of the formation of seed supermassive black holes via direct collapse in dark matter (DM) halos, in the cosmological context. We perform high-resolution zoom-in simulations of such collapse at high-$z$. Using the adaptive mesh refinement code ENZO, we resolve the formation and growth of a DM halo, until its virial temperature reaches $\sim 10^4$K, atomic cooling turns on, and collapse ensues. We demonstrate that direct collapse proceeds in two stages, although they are not well separated. The first stage is triggered by the onset of atomic cooling, and leads to rapidly increasing accretion rate with radius, from $\dot M\sim 0.1\,M_\odot {\rm yr^{-1}}$ at the halo virial radius to few $M_\odot \,{\rm yr^{-1}}$, around the scale radius $R_{\rm s}\sim 30$pc of the NFW DM density profile. The second stage of the collapse commences when the gas density takes precedence over the DM density. This is associated with the gas decoupling from the DM gravitational potential. The ensuing collapse approximates that of an isothermal sphere with $\dot M ( r )\sim $const. We confirm that the gas loses its angular momentum through non-axisymmetric perturbations and gravitational torques, to overcome the centrifugal barrier. During the course of the collapse, this angular momentum transfer process happens on nearly all spatial scales, and the angular momentum vector of the gas varies with position and time. Collapsing gas also exhibits supersonic turbulent motions which suppress gas fragmentation, and are characterized by density PDF consisting of a lognormal part and a high-density power law tail., Comment: 14 pages and 11 figures accepted for publication in MNRAS