Large scale simulations of the jet-IGM interaction

In a parameter study extending to jet densities of $10^{-5}$ times the ambient one, I have recently shown that light large scale jets start their lives in a spherical bow shock phase. This allows an easy description of the sideways bow shock propagation in that phase. Here, I present new, bipolar, simulations of very light jets in 2.5D and 3D, reaching the observationally relevant scale of $>200$ jet radii. Deviations from the early bow shock propagation law are expected because of various effects. The net effect is, however, shown to remain small. I calculate the X-ray appearance of the shocked cluster gas and compare it to Cygnus A and 3C 317. Rings, bright spots and enhancements inside the radio cocoon may be explained.Comment: 8 pages, 5 figures, ApSS accepted, proceedings of the virtual jets 2003 conference in Dogliani/Italy, v3: funny and unimportant bug corrected, one reference adde

Probing gaseous halos of galaxies with radio jets

Reproduced with permission from Astronomy & Astrophysics. © 2019 ESOContext. Gaseous halos play a key role in understanding inflow, feedback, and the overall baryon budget in galaxies. Literature models predict transitions of the state of the gaseous halo between cold and hot accretion, winds, fountains, and hydrostatic halos at certain galaxy masses. Since luminosities of radio AGN are sensitive to halo densities, any significant transition would be expected to show up in the radio luminosities of large samples of galaxies. The LOw Frequency ARray (LOFAR) Two-Metre Sky Survey (LoTSS) has identified a galaxy stellar mass scale, 10 11 M ⊙, above which the radio luminosities increase disproportionately. Aims. We investigate if radio luminosities of galaxies, especially the marked rise at galaxy masses around 10 11 M ⊙, can be explained with standard assumptions regarding jet powers, scaling between black hole mass and galaxy mass, and gaseous halos. Methods. Based on observational data and theoretical constraints, we developed models for the radio luminosity of radio AGN in halos under infall, galactic wind, and hydrostatic conditions. We compared these models to LoTSS data for a large sample of galaxies in the mass range between 10 8.5 M ⊙ and 10 12 M ⊙. Results. Under the assumption that the same characteristic upper limit to jet powers known from high galaxy masses holds at all masses, we find the maximum radio luminosities for the hydrostatic gas halos to lie close to the upper envelope of the distribution of the LOFAR data. The marked rise in radio luminosity at 10 11 M ⊙ is matched in our model and is related to a significant change in halo gas density around this galaxy mass, which is a consequence of lower cooling rates at a higher virial temperature. Wind and infall models overpredict the radio luminosities for small galaxy masses and have no particular steepening of the run of the radio luminosities predicted at any galaxy mass. Conclusions. Radio AGN could have the same characteristic Eddington-scaled upper limit to jet powers in galaxies of all masses in the sample if the galaxies have hydrostatic gas halos in phases when radio AGN are active. We find no evidence of a change of the type of galaxy halo with the galaxy mass. Galactic winds and quasi-spherical cosmological inflow phases cannot frequently occur at the same time as powerful jet episodes unless the jet properties in these phases are significantly different from what we assumed in our model.Peer reviewedFinal Accepted Versio

Jets and multi-phase turbulence

Jets are observed to stir up multi-phase turbulence in the inter-stellar medium as well as far beyond the host galaxy. Here we present detailed simulations of this process. We evolve the hydrodynamics equations with optically thin cooling for a 3D Kelvin Helmholtz setup with one initial cold cloud. The cloud is quickly disrupted, but the fragments remain cold and are spread throughout our simulation box. A scale free isotropic Kolmogorov power spectrum is built up first on the large scales, and reaches almost down to the grid scale after the simulation time of ten million years. We find a pronounced peak in the temperature distribution at 14,000K. The luminosity of the gas in this peak is correlated with the energy. We interpret this as a realisation of the shock ionisation scenario. The interplay between shock heating and radiative cooling establishes the equilibrium temperature. This is close to the observed emission in some Narrow Line Regions. We also confirm the shift of the phase equilibrium, i.e. a lower (higher) level of turbulence produces a higher (lower) abundance of cold gas. The effect could plausibly lead to a high level of cold gas condensation in the cocoons of extragalactic jets, explaining the so called Alignment Effect.Comment: 4 pages, 3 figures, conference: "The central kpc", Crete, June 2008, Volume 79 of the Memorie della Societa Astronomica Italian

A Markov Chain Monte Carlo approach for measurement of jet precession in radio-loud active galactic nuclei

© 2020 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved.Jet precession can reveal the presence of binary systems of supermassive black holes. The ability to accurately measure the parameters of jet precession from radio-loud AGN is important for constraining the binary supermassive black hole population, which are expected as a result of hierarchical galaxy evolution. The age, morphology, and orientation along the line of sight of a given source often result in uncertainties regarding jet path. This paper presents a new approach for efficient determination of precession parameters using a 2D MCMC curve-fitting algorithm which provides us a full posterior probability distribution on the fitted parameters. Applying the method to Cygnus A, we find evidence for previous suggestions that the source is precessing. Interpreted in the context of binary black holes leads to a constraint of parsec scale and likely sub-parsec orbital separation for the putative supermassive binary.Peer reviewe

Preparation, Structure, and Reactivity of Nonstabilized Organoiron Compounds. Implications for Iron-Catalyzed Cross Coupling Reactions

A series of unprecedented organoiron complexes of the formal oxidation states −2, 0, +1, +2, and +3 is presented, which are largely devoid of stabilizing ligands and, in part, also electronically unsaturated (14-, 16-, 17- and 18-electron counts). Specifically, it is shown that nucleophiles unable to undergo β-hydride elimination, such as MeLi, PhLi, or PhMgBr, rapidly reduce Fe(3+) to Fe(2+) and then exhaustively alkylate the metal center. The resulting homoleptic organoferrate complexes [(Me4Fe)(MeLi)][Li(OEt2)]2 (3) and [Ph4Fe][Li(Et2O)2][Li(1,4-dioxane)] (5) could be characterized by X-ray crystal structure analysis. However, these exceptionally sensitive compounds turned out to be only moderately nucleophilic, transferring their organic ligands to activated electrophiles only, while being unable to alkylate (hetero)aryl halides unless they are very electron deficient. In striking contrast, Grignard reagents bearing alkyl residues amenable to β-hydride elimination reduce FeXn (n = 2, 3) to clusters of the formal composition [Fe(MgX)2]n. The behavior of these intermetallic species can be emulated by structurally well-defined lithium ferrate complexes of the type [Fe(C2H4)4][Li(tmeda)]2 (8), [Fe(cod)2][Li(dme)]2 (9), [CpFe(C2H4)2][Li(tmeda)] (7), [CpFe(cod)][Li(dme)] (11), or [Cp*Fe(C2H4)2][Li(tmeda)] (14). Such electron-rich complexes, which are distinguished by short intermetallic Fe−Li bonds, were shown to react with aryl chlorides and allyl halides; the structures and reactivity patterns of the resulting organoiron compounds provide first insights into the elementary steps of low valent iron-catalyzed cross coupling reactions of aryl, alkyl, allyl, benzyl, and propargyl halides with organomagnesium reagents. However, the acquired data suggest that such C−C bond formations can occur, a priori, along different catalytic cycles shuttling between metal centers of the formal oxidation states Fe(+1)/Fe(+3), Fe(0)/Fe(+2), and Fe(−2)/Fe(0). Since these different manifolds are likely interconnected, an unambiguous decision as to which redox cycle dominates in solution remains difficult, even though iron complexes of the lowest accessible formal oxidation states promote the reactions most effectively

Numerical modelling of the lobes of radio galaxies – Paper V: universal pressure profile cluster atmospheres

© 2023 The Author(s). Published by Oxford University Press on behalf of Royal Astronomical Society. This is the final manuscript version of an article which has been published in final form at https://doi.org/10.1093/mnras/stad2933We present relativistic magnetohydrodynamic modelling of jets running into hydrostatic, spherically symmetric cluster atmospheres. For the first time in a numerical simulation, we present model cluster atmospheres based upon the universal pressure profile (UPP), incorporating a temperature profile for a ‘typical’ self-similar atmosphere described by only one parameter – M 500. We explore a comprehensive range of realistic atmospheres and jet powers and derive dynamic, energetic, and polarimetric data which provide insight into what we should expect of future high-resolution studies of AGN outflows. From the simulated synchrotron emission maps which include Doppler beaming we find sidedness distributions that agree well with observations. We replicated a number of findings from our previous work, such as higher power jets inflating larger aspect-ratio lobes, and the cluster environment impacting the distribution of energy between the lobe and shocked regions. Comparing UPP and β-profiles we find that the cluster model chosen results in a different morphology for the resultant lobes with the UPP more able to clear lobe material from the core; and that these different atmospheres influence the ratio between the various forms of energy in the fully developed lobes. This work also highlights the key role played by Kelvin–Helmholtz instabilities in the formation of realistic lobe aspect ratios. Our simulations point to the need for additional lobe-widening mechanisms at high jet powers, for example jet precession. Given that the UPP is our most representative general cluster atmosphere, these numerical simulations represent the most realistic models yet for spherically symmetric atmospheres.Peer reviewe

New mechanisms for forming multiple hotspots in radio jets

Hotspots of radio galaxies are regions of shock-driven particle acceleration. Multiple hotspots have long been identified as potential indicators of jet movement or precession. Two frequent explanations describe a secondary hotspot as either the location of a prior jet termination point, or a deflected backflow-driven shock: the so-called Dentist's Drill and Splatter Spot models. We created high-resolution simulations of precessing jets with a range of parameters. In addition to the existing mechanisms, our results show three additional mechanisms for multiple hotspot formation: (1) the splitting of a large terminal hotspots into passive and active components; (2) jet stream splitting resulting in two active hotspots; (3) dynamic multiple hotspot complexes that form as a result of jet termination in a turbulent cocoon, linked here to rapid precession. We show that these distinct types of multiple hotspots are difficult to differentiate in synthetic radio maps, particularly hotspot complexes which can easily be mistaken for the jet itself. We discuss the implication for hypothesised binary supermassive black hole systems where jet precession is a key component of the morphology, and show a selection of potential precession candidates found using the LOFAR Two-Metre Sky Survey Data Release 2 (LoTSS DR2).Comment: 14 pages, 10 figure

Could kilomasers pinpoint supermassive stars?

A strong nuclear kilomaser, W1, has been found in the nearby galaxy NGC 253, associated with a forming super star cluster. Kilomasers could arise from the accretion disc around supermassive stars (>10^3 Msun), hypothetical objects that have been proposed as polluters responsible for the chemical peculiarities in globular clusters. The supermassive stars would form via runaway collisions, simultaneously with the cluster. Their discs are perturbed by stellar flybys, inspiralling and colliding stars. This raises the question if an accretion disc would at all be able to survive in such a dynamic environment and mase water lines. We investigated what the predicted maser spectrum of such a disc would look like using 2D hydrodynamic simulations and compared this to the W1 kilomaser. We derived model maser spectra from the simulations by using a general maser model for appropriate disc temperatures. All our model discs survived. The model maser spectra for the most destructive case for the simulations of M = 1000 Msun are a reasonable match with the W1 kilomaser spectrum in terms of scaling, flux values and some of the signal trends. Details in the spectrum suggest that a star of a few 1000 Msun might fit even better, with 10,000 Msun clearly giving too large velocities. Our investigations thus support the hypothesis that kilomasers could pinpoint supermassive stars.Comment: 14 pages, 7 figure

RAiSE III : 3C radio AGN energetics and composition

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society. © 2017 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved.Kinetic jet power estimates based exclusively on observed monochromatic radio luminosities are highly uncertain due to confounding variables and a lack of knowledge about some aspects of the physics of active galactic nuclei (AGNs). We propose a new methodology to calculate the jet powers of the largest, most powerful radio sources based on combinations of their size, lobe luminosity, and shape of their radio spectrum; this approach avoids the uncertainties encountered by previous relationships. The outputs of our model are calibrated using hydrodynamical simulations and tested against independent X-ray inverse-Compton measurements. The jet powers and lobe magnetic field strengths of radio sources are found to be recovered using solely the lobe luminosity and spectral curvature, enabling the intrinsic properties of unresolved high-redshift sources to be inferred. By contrast, the radio source ages cannot be estimated without knowledge of the lobe volumes. The monochromatic lobe luminosity alone is incapable of accurately estimating the jet power or source age without knowledge of the lobe magnetic field strength and size, respectively. We find that, on average, the lobes of the Third Cambridge Catalogue of Radio Sources (3C) have magnetic field strengths approximately a factor three lower than the equipartition value, inconsistent with equal energy in the particles and the fields at the 5σ level. The particle content of 3C radio lobes is discussed in the context of complementary observations; we do not find evidence favouring an energetically dominant proton population.Peer reviewe

Constraints on positron annihilation kinematics in the inner Galaxy

Context. The annihilation of cosmic positrons with electrons in the interstellar medium results in the strongest persistent γ-ray line signal in the sky. For the past 50 yr, this 511 keV emission - predominantly from the galactic bulge region and from a low surface-brightness disk - has puzzled observers and theoreticians. A key issue for understanding positron astrophysics is found in cosmic-ray propagation, especially at low kinetic energies (≲ 10 MeV). Aims. We want to shed light on how positrons propagate and the resulting morphology of the annihilation emission. We approach this "positron puzzle" by inferring kinematic information of the 511 keV line in the inner radian of the Galaxy. This constrains propagation scenarios and positron source populations in the Milky Way. Methods. By dissecting the positron annihilation emission as measured with INTEGRAL/SPI, we derived spectra for individual and independent regions in the sky. The centroid energies of these spectra around the 511 keV line are converted into Doppler shifts, representing the line-of-sight velocity along different galactic longitudes. This results in a longitude-velocity diagram of positron annihilation. From high-resolution spectra, we also determined Doppler-broadening from γ-ray line shape parameters to study annihilation conditions as they vary with galactic longitude. Results. We found line-of-sight velocities in the 511 keV line that are consistent with zero, as well as with galactic rotation from CO measurements (2-3 km s -1 deg -1), and measurements of radioactive 26Al (7.5-9.5 km s -1 deg -1). The velocity gradient in the inner ±30° is determined to be 4 ± 6 km s -1 deg -1. The width of the 511 keV line is constant as a function of longitude at 2.43 ± 0.14 keV, with possibly different values towards the disk. The positronium fraction is found to be 1.0 along the galactic plane. Conclusions. The weak signals in the disk leave the question open of whether positron annihilation is associated with the high velocities seen in 26Al or rather with ordinarily rotating components of the Milky Way's interstellar medium. We confirm previous results that positrons are slowed down to the 10 eV energy scale before annihilation and constrain bulk Doppler-broadening contributions to ≲ 1.25 keV in the inner radian. Consequently, the true annihilation conditions remain unclear.Peer reviewedFinal Accepted Versio