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

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

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

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

Comparing simulated 26^{26}Al maps to gamma-ray measurements

© ESO 2019.Context. The diffuse gamma-ray emission of $^{26}{\rm Al}$ at 1.8 MeV reflects ongoing nucleosynthesis in the Milky Way, and traces massive-star feedback in the interstellar medium due to its 1 Myr radioactive lifetime. Interstellar-medium morphology and dynamics are investigated in astrophysics through 3D hydrodynamic simulations in fine detail, as only few suitable astronomical probes are available. Aims. We compare a galactic-scale hydrodynamic simulation of the Galaxy's interstellar medium, including feedback and nucleosynthesis, with gamma-ray data on $^{26}{\rm Al}$ emission in the Milky Way extracting constraints that are only weakly dependent on the particular realisation of the simulation or Galaxy structure. Methods. Due to constraints and biases in both the simulations and the gamma-ray observations, such comparisons are not straightforward. For a direct comparison, we perform maximum likelihood fits of simulated sky maps as well as observation-based maximum entropy maps to measurements with INTEGRAL/SPI. To study general morphological properties, we compare the scale heights of $^{26}{\rm Al}$ emission produced by the simulation to INTEGRAL/SPI measurements.} Results. The direct comparison shows that the simulation describes the observed inner Galaxy well, but differs significantly from the observed full-sky emission morphology. Comparing the scale height distribution, we see similarities for small scale height features and a mismatch at larger scale heights. We attribute this to the prominent foreground emission sites that are not captured by the simulation.Peer reviewedFinal Accepted Versio

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

How frequent are close supermassive binary black holes in powerful jet sources?

24 pages, 36 figures. © 2018 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/journals/pages/open_access/funder_policies/chorus/standard_publication_model)Supermassive black hole binariesmay be detectable by an upcoming suite of gravitationalwave experiments. Their binary nature can also be revealed by radio jets via a short-period precession driven by the orbital motion as well as the geodetic precession at typically longer periods. We have investigated Karl G. Jansky Very Large Array and Multi-Element Radio Linked Interferometer Network (MERLIN) radio maps of powerful jet sources for morphological evidence of geodetic precession. For perhaps the best-studied source, Cygnus A, we find strong evidence for geodetic precession. Projection effects can enhance precession features, for which we find indications in strongly projected sources. For a complete sample of 33 3CR radio sources, we find strong evidence for jet precession in 24 cases (73 per cent). The morphology of the radio maps suggests that the precession periods are of the order of 10 6- 10 7 yr. We consider different explanations for the morphological features and conclude that geodetic precession is the best explanation. The frequently observed gradual jet angle changes in samples of powerful blazars can be explained by orbital motion. Both observations can be explained simultaneously by postulating that a high fraction of powerful radio sources have subparsec supermassive black hole binaries.We consider complementary evidence and discuss if any jetted supermassive black hole with some indication of precession could be detected as individual gravitational wave source in the near future. This appears unlikely, with the possible exception of M87.Peer reviewedFinal Published versio

26Al kinematics: superbubbles following the spiral arms? : Constraints from the statistics of star clusters and HI supershells

High energy resolution spectroscopy of the 1.8 MeV radioactive decay line of 26Al with the SPI instrument on board the INTEGRAL satellite has recently revealed that diffuse 26Al has large velocities in comparison to other components of the interstellar medium in the Milky Way. 26Al shows Galactic rotation in the same sense as the stars and other gas tracers, but reaches excess velocities up to 300 km s−1Peer reviewe

The non-thermal superbubble in IC 10 : the generation of cosmic ray electrons caught in the act

Superbubbles are crucial for stellar feedback, with supposedly high (of the order of 10 per cent) thermalization rates. We combined multiband radio continuum observations from the Very Large Array (VLA) with Effelsberg data to study the non-thermal superbubble (NSB) in IC 10, a starburst dwarf irregular galaxy in the Local Group. Thermal emission was subtracted using a combination of Balmer Hα and VLA 32 GHz continuum maps. The bubble’s nonthermal spectrum between 1.5 and 8.8 GHz displays curvature and can be well fitted with a standard model of an ageing cosmic ray electron population. With a derived equipartition magnetic field strength of 44 ±8 μG, and measuring the radiation energy density from Spitzer MIPS maps as 5±1×10−11 erg cm−3, we determine, based on the spectral curvature, a spectral age of the bubble of 1.0 ± 0.3 Myr. Analysis of the LITTLE THINGS HI data cube shows an expanding HI hole with 100 pc diameter and a dynamical age of 3.8 ± 0.3 Myr, centred to within 16 pc on IC 10 X-1, a massive stellar mass black hole (M > 23 M⊙). The results are consistent with the expected evolution for a superbubble with a few massive stars, where a very energetic event like a Type Ic supernova/hypernova has taken place about 1 Myr ago. We discuss alternatives to this interpretationPeer reviewe