Accelerating AGN jets to parsec scales using general relativistic MHD simulations

Accreting black holes produce collimated outflows, or jets, that traverse many orders of magnitude in distance, accelerate to relativistic velocities, and collimate into tight opening angles. Of these, perhaps the least understood is jet collimation due to the interaction with the ambient medium. In order to investigate this interaction, we carried out axisymmetric general relativistic magnetohydrodynamic simulations of jets produced by a large accretion disc, spanning over 5 orders of magnitude in time and distance, at an unprecedented resolution. Supported by such a disc, the jet attains a parabolic shape, similar to the M87 galaxy jet, and the product of the Lorentz factor and the jet half-opening angle, $\gamma\theta\ll 1$, similar to values found from very long baseline interferometry (VLBI) observations of active galactic nuclei (AGN) jets; this suggests extended discs in AGN. We find that the interaction between the jet and the ambient medium leads to the development of pinch instabilities, which produce significant radial and lateral variability across the jet by converting magnetic and kinetic energy into heat. Thus pinched regions in the jet can be detectable as radiating hotspots and may provide an ideal site for particle acceleration. Pinching also causes gas from the ambient medium to become squeezed between magnetic field lines in the jet, leading to enhanced mass-loading of the jet and potentially contributing to the spine-sheath structure observed in AGN outflows.Comment: 18 pages, 24 figures, submitted to MNRAS, revised version. See our Youtube channel for accompanying animations: https://www.youtube.com/playlist?list=PLjldVlE2vDFzHMGn75tgc2Lod0kcTWZd

Do the Herschel cold clouds in the Galactic halo embody its dark matter?

Recent Herschel/SPIRE maps of the Small and Large Magellanic Clouds (SMC, LMC) exhibit in each thousands of clouds. Observed at 250 microns, they must be cold, T ~ 15 K, hence the name "Herschel cold clouds" (HCCs). From the observed rotational velocity profile and the assumption of spherical symmetry, the Galactic mass density is modeled in a form close to that of an isothermal sphere. If the HCCs constitute a certain fraction of it, their angular size distribution has a specified shape. A fit to the data deduced from the SMC/LMC maps supports this and yields for their radius 2.5 pc, with a small change when allowing for a spread in HCC radii. There are so many HCCs that they will make up all the missing Halo mass density if there is spherical symmetry and their average mass is of order 15,000 Mo. This compares well with the Jeans mass of circa 40,000 Mo and puts forward that the HCCs are in fact Jeans clusters, constituting all the Galactic dark matter and much of its missing baryons, a conclusion deduced before from a different field of the sky (Nieuwenhuizen, Schild and Gibson 2011). A preliminary analysis of the intensities yields that the Jeans clusters themselves may consist of some billion MACHOs of a few dozen Earth masses. With a size of dozens of solar radii, they would mostly obscure stars in the LMC, SMC and towards the Galactic center, and may thus have been overlooked in microlensing.Comment: Revised and corrected version, matches published version. Conclusions unchange

Formation of Precessing Jets by Tilted Black-hole Discs in 3D General Relativistic MHD Simulations

Gas falling into a black hole (BH) from large distances is unaware of BH spin direction, and misalignment between the accretion disc and BH spin is expected to be common. However, the physics of tilted discs (e.g., angular momentum transport and jet formation) is poorly understood. Using our new GPU-accelerated code H-AMR, we performed 3D general relativistic magnetohydrodynamic simulations of tilted thick accretion discs around rapidly spinning BHs, at the highest resolution to date. We explored the limit where disc thermal pressure dominates magnetic pressure, and showed for the first time that, for different magnetic field strengths on the BH, these flows launch magnetized relativistic jets propagating along the rotation axis of the tilted disc (rather than of the BH). If strong large-scale magnetic flux reaches the BH, it bends the inner few gravitational radii of the disc and jets into partial alignment with the BH spin. On longer time scales, the simulated disc-jet system as a whole undergoes Lense-Thirring precession and approaches alignment, demonstrating for the first time that jets can be used as probes of disc precession. When the disc turbulence is well-resolved, our isolated discs spread out, causing both the alignment and precession to slow down.Comment: MNRAS Letters, accepted. Animations available at https://www.youtube.com/playlist?list=PL39mDr1uU6a5RYZdXLAjKE1C_GAJkQJN

Simulation of the hydrogen ground state in Stochastic Electrodynamics

Stochastic electrodynamics is a classical theory which assumes that the physical vacuum consists of classical stochastic fields with average energy $\frac{1}{2}\hbar \omega$ in each mode, i.e., the zero-point Planck spectrum. While this classical theory explains many quantum phenomena related to harmonic oscillator problems, hard results on nonlinear systems are still lacking. In this work the hydrogen ground state is studied by numerically solving the Abraham -- Lorentz equation in the dipole approximation. First the stochastic Gaussian field is represented by a sum over Gaussian frequency components, next the dynamics is solved numerically using OpenCL. The approach improves on work by Cole and Zou 2003 by treating the full $3d$ problem and reaching longer simulation times. The results are compared with a conjecture for the ground state phase space density. Though short time results suggest a trend towards confirmation, in all attempted modelings the atom ionises at longer times.Comment: 20 pages, 9 figures. Published version, minor change

Spectral and Imaging properties of Sgr A* from High-Resolution 3D GRMHD Simulations with Radiative Cooling

The candidate supermassive black hole in the Galactic Centre, Sagittarius A* (Sgr A*), is known to be fed by a radiatively inefficient accretion flow (RIAF), inferred by its low accretion rate. Consequently, radiative cooling has in general been overlooked in the study of Sgr A*. However, the radiative properties of the plasma in RIAFs are poorly understood. In this work, using full 3D general-relativistic magneto-hydrodynamical simulations, we study the impact of radiative cooling on the dynamical evolution of the accreting plasma, presenting spectral energy distributions and synthetic sub-millimeter images generated from the accretion flow around Sgr A*. These simulations solve the approximated equations for radiative cooling processes self-consistently, including synchrotron, bremsstrahlung, and inverse Compton processes. We find that radiative cooling plays an increasingly important role in the dynamics of the accretion flow as the accretion rate increases: the mid-plane density grows and the infalling gas is less turbulent as cooling becomes stronger. The changes in the dynamical evolution become important when the accretion rate is larger than $10^{-8}\,M_{\odot}~{\rm yr}^{-1}$ ($\gtrsim 10^{-7} \dot{M}_{\rm Edd}$, where $\dot{M}_{\rm Edd}$ is the Eddington accretion rate). The resulting spectra in the cooled models also differ from those in the non-cooled models: the overall flux, including the peak values at the sub-mm and the far-UV, is slightly lower as a consequence of a decrease in the electron temperature. Our results suggest that radiative cooling should be carefully taken into account in modelling Sgr A* and other low-luminosity active galactic nuclei that have a mass accretion rate of $\dot{M} > 10^{-7}\,\dot{M}_{\rm Edd}$.Comment: 16 pages, 16 figures, accepted for publication in MNRA

Winds and Disk Turbulence Exert Equal Torques on Thick Magnetically Arrested Disks

The conventional accretion disk lore is that magnetized turbulence is the principal angular momentum transport process that drives accretion. However, when dynamically important magnetic fields thread an accretion disk, they can produce mass and angular momentum outflows that also drive accretion. Yet, the relative importance of turbulent and wind-driven angular momentum transport is still poorly understood. To probe this question, we analyze a long-duration ($1.2 \times 10^5 r_{\rm g}/c$) simulation of a rapidly rotating ($a=0.9$) black hole (BH) feeding from a thick ($H/r\sim0.3$), adiabatic, magnetically arrested disk (MAD), whose dynamically-important magnetic field regulates mass inflow and drives both uncollimated and collimated outflows (e.g., "winds" and "jets", respectively). By carefully disentangling the various angular momentum transport processes occurring within the system, we demonstrate the novel result that both disk winds and disk turbulence extract roughly equal amounts of angular momentum from the disk. We find cumulative angular momentum and mass accretion outflow rates of $\dot{L}\propto r^{0.9}$ and $\dot{M}\propto r^{0.4}$, respectively. This result suggests that understanding both turbulent and laminar stresses is key to understanding the evolution of systems where geometrically thick MADs can occur, such as the hard state of X-ray binaries, low-luminosity active galactic nuclei, some tidal disruption events, and possibly gamma ray bursts.Comment: 15 pages, 6 figures. Submitted to ApJ. Comments welcom

Millimeter Observational Signatures of Flares in Magnetically Arrested Black Hole Accretion Models

In general relativistic magneto-hydrodynamic (GRMHD) simulations, accreted magnetic flux on the black hole horizon episodically decays, during which magnetic reconnection heats up the plasma near the horizon, potentially powering high-energy flares like those observed in M87* and Sgr A*. We study the mm observational counterparts of such flaring episodes. The change in 230 GHz flux during the expected high energy flares depends primarily on the efficiency of accelerating $\gamma \gtrsim 100$ ($T_e \gtrsim 10^{11}$ K) electrons. For models in which the electrons are heated to $T_e \sim 10^{11}$ K during flares, the hot plasma produced by reconnection significantly enhances 230 GHz emission and increases the size of the 230 GHz image. By contrast, for models in which the electrons are heated to higher temperatures (which we argue are better motivated), the reconnection-heated plasma is too hot to produce significant 230 GHz synchrotron emission, and the 230 GHz flux decreases during high energy flares. We do not find a significant change in the mm polarization during flares as long as the emission is Faraday thin. We also present expectations for the ring-shaped image as observed by the Event Horizon Telescope during flares, as well as multi-wavelength synchrotron spectra. Our results highlight several limitations of standard post-processing prescriptions for the electron temperature in GRMHD simulations. We also discuss the implications of our results for current and future observations of flares in Sgr A*, M87*, and related systems. Appendices contain detailed convergence studies with respect to resolution and plasma magnetization.Comment: 11+7 pages, 9+7 figures, 1 table, accepted by MNRA

Jets with a Twist: Emergence of FR0 Jets in 3D GRMHD Simulation of Zero Angular Momentum Black Hole Accretion

Spinning supermassive black holes (BHs) in active galactic nuclei (AGN) magnetically launch relativistic collimated outflows, or jets. Without angular momentum supply, such jets are thought to perish within $3$ orders of magnitude in distance from the BH, well before reaching kpc-scales. We study the survival of such jets at the largest scale separation to date, via 3D general relativistic magnetohydrodynamic simulations of rapidly spinning BHs immersed into uniform zero-angular-momentum gas threaded by weak vertical magnetic field. We place the gas outside the BH sphere of influence, or the Bondi radius, chosen much larger than the BH gravitational radius, $R_\text{B}=10^3R_\text{g}$. The BH develops dynamically-important large-scale magnetic fields, forms a magnetically-arrested disk (MAD), and launches relativistic jets that propagate well outside $R_\text{B}$ and suppress BH accretion to $1.5\%$ of the Bondi rate, $\dot{M}_\text{B}$. Thus, low-angular-momentum accretion in the MAD state can form large-scale jets in Fanaroff-Riley (FR) type I and II galaxies. Subsequently, the disk shrinks and exits the MAD state: barely a disk (BAD), it rapidly precesses, whips the jets around, globally destroys them, and lets $5-10\%$ of $\dot{M}_\text{B}$ reach the BH. Thereafter, the disk starts rocking back and forth by angles $90-180^\circ$: the rocking accretion disk (RAD) launches weak intermittent jets that spread their energy over a large area and suppress BH accretion to $\lesssim 2 \% ~ \dot{M}_\text{B}$. Because BAD and RAD states tangle up the jets and destroy them well inside $R_\text{B}$, they are promising candidates for the more abundant, but less luminous, class of FR0 galaxies

Validation Study Report: Performance assessment of the AR-CALUX® in vitro method: to support the development of an international test guideline for Androgen Receptor Transactivation Assays (ARTA) for the detection of compounds with (anti)androgenic potential

The JRC’s EU Reference Laboratory for alternatives to animal testing (EURL ECVAM) conducted a validation study of the AR-CALUX in vitro method. The method is applied to the detection of compounds with endocrine disrupting potential and more specifically (anti)androgen activity. The objectives of the study included assessing the reproducibility (within and between laboratories) and the relevance of the in vitro method. The participating laboratories included three test facilities from the European Union Network of Validation Laboratories for alternative methods (EU-NETVAL), being RISE, Covance, and Charles River, and, the test method submitter Biodetection Systems. The validation study report presents the results of the method performed on 46 test chemicals, evaluated for reproducibility within and between laboratories, variability within the measurements, and classification. A comparison of the classifications was made with publicly available ARTA classifications. It was concluded that this test method is reliable, has high reproducibility, low variability and merits proposal to OECD for the development of a test guideline. The statistical report and the final SOP are part of this report and can be found as separate files.JRC.F.3-Chemicals Safety and Alternative Method