Search

The precession of astrophysical jets produced by active-galactic nuclei is likely related to the dynamics of the accretion disks surrounding the central supermassive black holes (BHs) from which jets are launched. The two main mechanisms that can drive jet precession arise from Lense-Thirring precession and tidal torquing. These can explain direct and indirect observations of precessing jets; however, such explanations often utilize crude approximations of the disk evolution and observing jet precession can be challenging with electromagnetic facilities. Simultaneously, the Laser Interferometer Space Antenna (LISA) is expected to measure gravitational waves from the mergers of massive binary BHs with high accuracy and probe their progenitor evolution. In this paper, we connect the LISA detectability of binary BH mergers to the possible jet precession during their progenitor evolution. We make use of a semi-analytic model that self-consistently treats disk-driven BH alignment and binary inspiral and includes the possibility of disk breaking. We find that tidal torquing of the accretion disk provides a wide range of jet precession timescales depending on the binary separation and the spin direction of the BH from which the jet is launched. Efficient disk-driven BH alignment results in shorter timescales of $\sim 1$ yr which are correlated with higher LISA signal-to-noise ratios. Disk breaking results in the longest possible times of $\sim 10^7$ yrs, suggesting a deep interplay between the disk critical obliquity (i.e. where the disk breaks) and jet precession. Studies such as ours will help to reveal the cosmic population of precessing jets that are detectable with gravitational waves.

Odd Radio Circles (ORCs) are a class of low surface brightness, circular objects approximately one arcminute in diameter. ORCs were recently discovered in the Australian Square Kilometre Array Pathfinder (ASKAP) data, and subsequently confirmed with follow-up observations on other instruments, yet their origins remain uncertain. In this paper, we suggest that ORCs could be remnant lobes of powerful radio galaxies, re-energised by the passage of a shock. Using relativistic hydrodynamic simulations with synchrotron emission calculated in post-processing, we show that buoyant evolution of remnant radio lobes is alone too slow to produce the observed ORC morphology. However, the passage of a shock can produce both filled and edge-brightnened ORC-like morphologies for a wide variety of shock and observing orientations. Circular ORCs are predicted to have host galaxies near the geometric centre of the radio emission, consistent with observations of these objects. Significantly offset hosts are possible for elliptical ORCs, potentially causing challenges for accurate host galaxy identification. Observed ORC number counts are broadly consistent with a paradigm in which moderately powerful radio galaxies are their progenitors., Comment: 25 pages, 15 figures. Accepted for publication in PASA

During the first four all-sky surveys eRASS:4 carried out from December 2019 to 2021, the extended Roentgen Survey with an Imaging Telescope Array (eROSITA) on board Spektrum-Roentgen-Gamma (Spektr-RG, SRG) observed the Galactic HII region Carina nebula. We analysed the eRASS:4 data to study the distribution and the spectral properties of the hot interstellar plasma and the bright stellar sources in the Carina nebula. Spectral extraction regions of the diffuse emission were defined based on X-ray spectral morphology and multi-wavelength data. The spectra were fit with a combination of thermal and non-thermal emission models. X-ray bright point sources in the Carina nebula are the colliding wind binary $\eta$ Car, several O stars, and Wolf-Rayet (WR) stars. We extracted the spectrum of the brightest stellar sources, which can be well fit with a multi-component thermal plasma model. The spectra of the diffuse emission in the brighter parts of the Carina nebula is well reproduced by two thermal models, a lower-temperature component ($\sim$0.2 keV) and a higher-temperature component (0.6 - 0.8 keV). An additional non-thermal component dominates the emission above $\sim$1 keV in the central region around $\eta$ Car and the other massive stars. Significant orbital variation of the X-ray flux was measured for $\eta$ Car, WR22 and WR25. $\eta$ Car requires an additional time-variable thermal component in the spectral model, which is associated to the wind-wind-collision zone. Properties like temperature, pressure, and luminosity of the X-ray emitting plasma in the Carina nebula derived from the eROSITA data are consistent with theoretical calculations of emission from superbubbles. It confirms that the X-ray emission is caused by the hot plasma inside the Carina nebula which has been shocked-heated by the stellar winds of the massive stars, in particular, of $\eta$ Car., Comment: Accepted for publication in Astronomy & Astrophysics

Active galactic nuclei (AGN) play an integral role in galaxy formation and evolution by influencing galaxies and their environments through radio jet feedback. Historically, interpreting observations of radio galaxies and quantifying radio jet feedback has been challenging due to degeneracies between their physical parameters. In particular, it is well-established that different combinations of jet kinetic power and environment density can yield indistinguishable radio continuum properties, including apparent size and Stokes I luminosity. We present an approach to breaking this degeneracy by probing the line-of-sight environment with Faraday rotation. We study this effect in simulations of three-dimensional relativistic magnetohydrodynamic AGN jets in idealised environments with turbulent magnetic fields. We generate synthetic Stokes I emission and Faraday rotation measure (RM) maps, which enable us to distinguish between our simulated sources. We find enhanced RMs near the jet head and lobe edges and an RM reversal across the jet axis. We show that increasing the environment density and the average cluster magnetic field strength broadens the distribution of Faraday rotation measure values. We study the depolarisation properties of our sources, finding that the hotspot regions depolarise at lower frequencies than the lobes. We quantify the effect of depolarisation on the RM distribution, finding that the frequency at which the source is too depolarised to measure the RM distribution accurately is a probe of environmental properties. This technique offers a range of new opportunities for upcoming surveys, including probing radio galaxy environments and determining more accurate estimates of the AGN feedback budget., Comment: 17 pages, 11 figures. Submitted for publication in MNRAS

We 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 $\beta$-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 (KH) 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., Comment: 20 pages, 18 figures, 2 tables, accepted for publication in MNRAS, minor correction to references in section 4.2

The conditions required for massive star formation are debated, particularly whether massive stars must form in conjunction with massive clusters. Some authors have advanced the view that stars of any mass (below the total cluster mass) can form in clusters of any mass with some probability (random sampling). Others pointed out that the scatter in the determinations of the most massive star mass for a given cluster mass was consistent with the measurement error, such that the mass of the most massive star was determined by the total cluster mass (optimal sampling). Here we investigate the relation between cluster mass (M\textsubscript{ecl}) and the maximum stellar mass (M\textsubscript{max}) using a suite of SPH simulations. Varying cloud mass and turbulence random seed results in a range of cluster masses which we compare with their respective maximum star masses. We find that more massive clusters will have, on average, higher mass stars with this trend being steeper at lower cluster masses ($M\textsubscript{max} \propto M\textsubscript{ecl}^{0.31}$ for $M\textsubscript{ecl}<500M\,_{\odot}$) and flattening at higher cluster masses ($M\textsubscript{max} \propto M\textsubscript{ecl}^{0.11}$ for $M\textsubscript{ecl}>500M\,_{\odot}$). This rules out purely stochastic star formation in our simulations. Significant scatter in the maximum masses with identical initial conditions also rules out the possibility that the relation is purely deterministic (that is that a given cluster mass will result in a specific maximum stellar mass). In conclusion our simulations disagree with both random and optimal sampling of the initial mass function., Comment: Accepted in MNRAS

The Compton Spectrometer and Imager (COSI) is a NASA Small Explorer (SMEX) satellite mission in development with a planned launch in 2027. COSI is a wide-field gamma-ray telescope designed to survey the entire sky at 0.2-5 MeV. It provides imaging, spectroscopy, and polarimetry of astrophysical sources, and its germanium detectors provide excellent energy resolution for emission line measurements. Science goals for COSI include studies of 0.511 MeV emission from antimatter annihilation in the Galaxy, mapping radioactive elements from nucleosynthesis, determining emission mechanisms and source geometries with polarization measurements, and detecting and localizing multimessenger sources. The instantaneous field of view for the germanium detectors is >25% of the sky, and they are surrounded on the sides and bottom by active shields, providing background rejection as well as allowing for detection of gamma-ray bursts and other gamma-ray flares over most of the sky. In the following, we provide an overview of the COSI mission, including the science, the technical design, and the project status., Comment: 8 pages

The Compton Spectrometer and Imager (COSI) is a selected Small Explorer (SMEX) mission launching in 2027. It consists of a large field-of-view Compton telescope that will probe with increased sensitivity the under-explored MeV gamma-ray sky (0.2-5 MeV). We will present the current status of cosipy, a Python library that will perform spectral and polarization fits, image deconvolution, and all high-level analysis tasks required by COSI's broad science goals: uncovering the origin of the Galactic positrons, mapping the sites of Galactic nucleosynthesis, improving our models of the jet and emission mechanism of gamma-ray bursts (GRBs) and active galactic nuclei (AGNs), and detecting and localizing gravitational wave and neutrino sources. The cosipy library builds on the experience gained during the COSI balloon campaigns and will bring the analysis of data in the Compton regime to a modern open-source likelihood-based code, capable of performing coherent joint fits with other instruments using the Multi-Mission Maximum Likelihood framework (3ML). In this contribution, we will also discuss our plans to receive feedback from the community by having yearly software releases accompanied by publicly-available data challenges.

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 figures

Diffuse gamma-ray line emission traces freshly produced radioisotopes in the interstellar gas, providing a unique perspective on the entire Galactic cycle of matter from nucleosynthesis in massive stars to their ejection and mixing in the interstellar medium. We aim at constructing a model of nucleosynthesis ejecta on galactic scale which is specifically tailored to complement the physically most important and empirically accessible features of gamma-ray measurements in the MeV range, in particular for decay gamma-rays such as $^{26}$Al, $^{60}$Fe or $^{44}$Ti. Based on properties of massive star groups, we developed a Population Synthesis Code which can instantiate galaxy models quickly and based on many different parameter configurations, such as the star formation rate, density profiles, or stellar evolution models. As a result, we obtain model maps of nucleosynthesis ejecta in the Galaxy which incorporate the population synthesis calculations of individual massive star groups. Based on a variety of stellar evolution models, supernova explodabilities, and density distributions, we find that the measured $^{26}$Al distribution from INTEGRAL/SPI can be explained by a Galaxy-wide population synthesis model with a star formation rate of $4$-$8\,\mathrm{M_{\odot}\,yr^{-1}}$ and a spiral-arm dominated density profile with a scale height of at least 700 pc. Our model requires that most massive stars indeed undergo a supernova explosion. This corresponds to a supernova rate in the Milky Way of $1.8$-$2.8$ per century, with quasi-persistent $^{26}$Al and $^{60}$Fe masses of $1.2$-$2.4\,\mathrm{M_{\odot}}$ and $1$-$6\,\mathrm{M_{\odot}}$, respectively. Comparing the simulated morphologies to SPI data suggests that a frequent merging of superbubbles may take place in the Galaxy, and that an unknown but strong foreground emission at 1.8 MeV could be present., Comment: 16 pages, 12 figures, 2 page appendix, accepted for publication in A&A

Diffuse gamma-ray emission from the decay of radioactive $^{26}$Al is a messenger from the nucleosynthesis activity in our current-day galaxy. Because this material is attributed to ejections from massive stars and their supernovae, the gamma-ray signal includes information about nucleosynthesis in massive star interiors as it varies with evolutionary stages, and about their feedback on the surrounding interstellar medium. Our method of population synthesis of massive-star groups has been refined as a diagnostic tool for this purpose. It allows to build a bottom-up prediction of the diffuse gamma-ray sky when known massive star group distributions and theoretical models of stellar evolution and core-collapse supernova explosions are employed. We find general consistency of an origin in such massive-star groups, in particular we also find support for the clumpy distribution of such source regions across the Galaxy, and characteristics of large cavities around these. A discrepancy in the integrated $^{26}$Al gamma-ray flux is interpreted as an indication for excess $^{26}$Al emission from nearby, distributed in cavities that extend over major regions of the sky., Comment: 6 pages, 2 figures. Invited contribution to conference 'Nuclear Physics in Astrophysics X', Geneva Sep 2022

The presence of radioactive $^{26}$Al at 1.8 MeV reflects ongoing nucleosynthesis in the Milky Way. Diffuse emission from its decay can be measured with gamma-ray telescopes in space. The intensity, line shape, and spatial distribution of the $^{26}$Al emission allow a study of these nucleosynthesis sources. The line parameters trace massive-star feedback in the interstellar medium due to its 1~My lifetime. We aim to deepen previous studies of the $^{26}$Al emission in the Milky Way, using all gamma-ray data including single and double events as collected with SPI on INTEGRAL from 2003 until 2020. We apply improved spectral response and background as evaluated from tracing spectral details over the entire mission. The exposure for Galactic $^{26}$Al emission is enhanced using all event types measured within SPI. We re-determine the intensity of Galactic $^{26}$Al emission across the entire sky, through maximum likelihood fits of simulated and model-built sky distributions to SPI spectra for single and for double detector hits. We find an all-sky flux of (1.84$\pm$0.03$)\times$10$^{-3}$~ph~cm$^{-2}$s$^{-1}$ in the 1.809~MeV line from $^{26}$Al, determined as fitted to sky distributions from previous observations with COMPTEL. Significant emission from higher latitudes indicate an origin from nearby massive-star groups and superbubbles, also supported by a bottom-up population synthesis model. The line centroid is found at (1809.83$\pm$0.04~keV, and line broadening from source kinematics integrated over the sky is (0.62$\pm0.3$)~keV (FWHM)., Comment: 7 pages, 9 figures. Accepted for publication in Astronomy & Astrophysics

We present the Cosmological Double Radio Active Galactic Nuclei (CosmoDRAGoN) project: a large suite of simulated AGN jets in cosmological environments. These environments sample the intra-cluster media of galaxy clusters that form in cosmological smooth particle hydrodynamics (SPH) simulations, which we then use as inputs for grid-based hydrodynamic simulations of radio jets. Initially conical jets are injected with a range of jet powers, speeds (both relativistic and non-relativistic), and opening angles; we follow their collimation and propagation on scales of tens to hundreds of kiloparsecs, and calculate spatially-resolved synthetic radio spectra in post-processing. In this paper, we present a technical overview of the project, and key early science results from six representative simulations which produce radio sources with both core- (Fanaroff-Riley Type I) and edge-brightened (Fanaroff-Riley Type II) radio morphologies. Our simulations highlight the importance of accurate representation of both jets and environments for radio morphology, radio spectra, and feedback the jets provide to their surroundings., Comment: 25 pages, 20 figures. Accepted for publication in PASA

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 figures

Star clusters are known to be formed in turbulent molecular clouds. How turbulence is driven in molecular clouds and what effect this has on star formation is still unclear. We compare a simulation setup with turbulent driving everywhere in a periodic box with a setup where turbulence is only driven around the outside of the box. We analyse the resulting gas distribution, kinematics, and the population of stars that are formed from the cloud. Both setups successfully produce a turbulent velocity field with a power law structure function, the externally driven cloud has a more central, monolithic, clump, while the fully driven cloud has many smaller, more dispersed, clumps. The star formation follows the cloud morphology producing large clusters, with high star forming efficiency in the externally driven simulations and sparse individual star formation with much lower star formation efficiency in the fully driven case. We conclude that the externally driven method, which resembles a Global Hierarchical Collapse (GHC) scenario, produces star clusters that more closely match with observations., Comment: 8 pages, 6 figures, accepted by MNRAS

We present a method for applying spatially resolved adiabatic and radiative loss processes to synthetic radio emission from hydrodynamic simulations of radio sources from active galactic nuclei (AGN). Lagrangian tracer particles, each representing an ensemble of electrons, are injected into simulations and the position, grid pressure, and time since the last strong shock are recorded. These quantities are used to track the losses of the electron packet through the radio source in a manner similar to the Radio AGN in Semi-analytic Environments (RAiSE) formalism, which uses global source properties to calculate the emissivity of each particle ex-situ. Freedom in the choice of observing parameters, including redshift, is provided through the post-processing nature of this approach. We apply this framework to simulations of jets in different environments, including asymmetric ones. We find a strong dependence of radio source properties on frequency and redshift, in good agreement with observations and previous modelling work. There is a strong evolution of radio spectra with redshift due to the more prominent inverse-Compton losses at high redshift. Radio sources in denser environments have flatter spectral indices, suggesting that spectral index asymmetry may be a useful environment tracer. We simulate intermediate Mach number jets that disrupt before reaching the tip of the lobe, and find that these retain an edge-brightened Fanaroff-Riley Type II morphology, with the most prominent emission remaining near the tip of the lobes for all environments and redshifts we study., Comment: 17 pages, 15 figures. Accepted for publication in MNRAS

We have carried out relativistic three-dimensional simulations of high-power radio sources propagating into asymmetric cluster environments. We offset the environment by 0 or 1 core radii (equal to 144 kpc), and incline the jets by 0, 15, or 45{\deg} away from the environment centre. The different environment encountered by each radio lobe provides a unique opportunity to study the effect of environment on otherwise identical jets. We find that the jets become unstable towards the end of the simulations, even with a Lorentz factor of 5; they nevertheless develop typical FR II radio morphology. The jets propagating into denser environments have consistently shorter lobe lengths and brighter hotspots, while the axial ratio of the two lobes is similar. We reproduce the recently reported observational anti-correlation between lobe length asymmetry and environment asymmetry, corroborating the notion that observed large-scale radio lobe asymmetry can be driven by differences in the underlying environment., Comment: 13 pages, 11 figures. Accepted for publication in MNRAS

The cosmic evolution of the chemical elements from the Big Bang to the present time is driven by nuclear fusion reactions inside stars and stellar explosions. A cycle of matter recurrently re-processes metal-enriched stellar ejecta into the next generation of stars. The study of cosmic nucleosynthesis and of this matter cycle requires the understanding of the physics of nuclear reactions, of the conditions at which the nuclear reactions are activated inside the stars and stellar explosions, of the stellar ejection mechanisms through winds and explosions, and of the transport of the ejecta towards the next cycle, from hot plasma to cold, star-forming gas. Due to the long timescales of stellar evolution, and because of the infrequent occurrence of stellar explosions, observational studies are challenging. Due to their radioactive lifetime of million years, the 26Al and 60Fe isotopes are suitable to characterise simultaneously the processes of nuclear fusion reactions and of interstellar transport. We describe and discuss the nuclear reactions involved in the production and destruction of 26Al and 60Fe, the key characteristics of the stellar sites of their nucleosynthesis and their interstellar journey after ejection from the nucleosynthesis sites. We connect the theoretical astrophysical aspects to the variety of astronomical messengers, from stardust and cosmic-ray composition measurements, through observation of gamma rays produced by radioactivity, to material deposited in deep-sea ocean crusts and to the inferred composition of the first solids that have formed in the Solar System. We show that considering measurements of the isotopic ratio of 26Al to 60Fe eliminate some of the unknowns when interpreting astronomical results, and discuss the lessons learned from these two isotopes on cosmic chemical evolution., Comment: 51 pages, 44 figures; review from an ISSI workshop series; accepted for publication in PASA

Reticulum II (Ret II) is a satellite galaxy of the Milky Way and presents a prime target to investigate the nature of dark matter (DM) because of its high mass-to-light ratio. We evaluate a dedicated INTEGRAL observation campaign data set to obtain $\gamma$-ray fluxes from Ret II and compare those with expectations from DM. Ret II is not detected in the $\gamma$-ray band 25--8000 keV, and we derive a flux limit of $\lesssim 10^{-8}\,\mathrm{erg\,cm^{-2}\,s^{-1}}$. The previously reported 511 keV line is not seen, and we find a flux limit of $\lesssim 1.7 \times 10^{-4}\,\mathrm{ph\,cm^{-2}\,s^{-1}}$. We construct spectral models for primordial black hole (PBH) evaporation and annihilation/decay of particle DM, and subsequent annihilation of positrons produced in these processes. We exclude that the totality of DM in Ret II is made of a monochromatic distribution of PBHs of masses $\lesssim 8 \times 10^{15}\,\mathrm{g}$. Our limits on the velocity-averaged DM annihilation cross section into $e^+e^-$ are $\langle \sigma v \rangle \lesssim 5 \times 10^{-28} \left(m_{\rm DM} / \mathrm{MeV} \right)^{2.5}\,\mathrm{cm^3\,s^{-1}}$. We conclude that analysing isolated targets in the MeV $\gamma$-ray band can set strong bounds on DM properties without multi-year data sets of the entire Milky Way, and encourage follow-up observations of Ret II and other dwarf galaxies., Comment: 8 pages, 10 figures, submitted

The ARCADE 2 balloon bolometer along with a number of other instruments have detected what appears to be a radio synchrotron background at frequencies below about 3 GHz. Neither extragalactic radio sources nor diffuse Galactic emission can currently account for this finding. We use the locally measured Cosmic ray electron population, demodulated for effects of the Solar wind, and other observational constraints combined with a turbulent magnetic field model to predict the radio synchrotron emission for the Local Bubble. We find that the spectral index of the modelled radio emission is roughly consistent with the radio background. Our model can approximately reproduce the observed antenna temperatures for a mean magnetic field strength B between 3-5 nT. We argue that this would not violate observational constraints from pulsar measurements. However, the curvature in the predicted spectrum would mean that other, so far unknown sources would have to contribute below 100 MHz. Also, the magnetic energy density would then dominate over thermal and cosmic ray electron energy density, likely causing an inverse magnetic cascade with large variations of the radio emission in different sky directions as well as high polarisation. We argue that this disagrees with several observations and thus that the magnetic field is probably much lower, quite possibly limited by equipartition with the energy density in relativistic or thermal particles (B = 0.2-0.6 nT). In the latter case, we predict a contribution of the Local Bubble to the unexplained radio background at most at the per cent level., Comment: 8 pages, 2 figures, accepted by MNRAS

Radioactive 26Al is an excellent tracer for metal ejection in the Milky Way, and can provide a direct constraint on the modelling of supernova feedback in galaxy evolution. Gamma-ray observations of the 26Al decay line have found high velocities and hence require a significant fraction of the Galactic 26Al in the hot component. At the same time, meteoritic data combined with simulation results suggest that a significant amount of 26Al makes its way into stars before decay. We investigated the distribution into hot and cold channels with a simulation of a Milky-Way-like galaxy with massive-star feedback in superbubbles and with ejecta traced by 26Al. About 30-40 per cent of the ejecta remain hot, with typical cooling times of the order Gyr. 26Al traces the footpoints of a chimney-fed outflow that mixes metals turbulently into the halo of the model galaxy on a scale of at least 50 kpc. The rest diffuses into cold gas with temperatures less than about 10,000 K, and may therefore be quickly available for star formation. We discuss the robustness of the result by comparison to a simulation with a different global flow pattern. The branching ratio into hot and cold components is comparable to that of longer term average results from chemical evolution modelling of galaxies, clusters and the intracluster medium., Comment: 10 pages, 10 figures, accepted for MNRAS

Measurement and astrophysical interpretation of characteristic gamma-ray lines from nucleosynthesis was one of the prominent science goals of the INTEGRAL mission and in particular its spectrometer SPI. Emission from 26Al and from 60Fe decay lines originates from accumulated ejecta of nucleosynthesis sources, and appears diffuse in nature. 26Al and 60Fe are believed to originate mostly from massive star clusters. Gamma-ray observations open an interesting window to trace the fate and flow of nucleosynthesis ejecta, after they have left the immediate sources and their birth sites, and on their path to mix with ambient interstellar gas. The INTEGRAL 26Al emission image confirms earlier findings of clumpiness and an extent along the entire plane of the Galaxy, supporting its origin from massive-star groups. INTEGRAL spectroscopy resolved the line and found Doppler broadenings and systematic shifts from large-scale galactic rotation. But an excess velocity of ~200 km/s suggests that 26Al decays preferentially within large superbubbles that extend in forward directions between spiral arms. The detection of 26Al line emission from nearby Orion and the Eridanus superbubble supports this interpretation. Positrons from beta+ decays of 26Al and other nucleosynthesis ejecta have been found to not explain the morphology of positron annihilation gamma-rays at 511 keV that have been measured by INTEGRAL. The 60Fe signal measured by INTEGRAL is diffuse but too weak for an imaging interpretation, an origin from point-like/concentrated sources is excluded. The 60Fe/26Al ratio is constrained to a range 0.2-0.4. Beyond improving precision of these results, diffuse nucleosynthesis contributions from novae (through 22Na radioactivity) and from past neutron star mergers in our Galaxy (from r-process radioactivity) are exciting new prospects for the remaining mission extensions., Comment: 38 pages, 12 Figures. Contribution to Review Book "15 years of INTEGRAL". Accepted for publication in New Astronomy Reviews

The prospect of relativistic jets exhibiting complex morphologies as a consequence of geodetic precession has long been hypothesised. We have carried out a 3D hydrodynamics simulation study varying the precession cone angle, jet injection speed and number of turns per simulation time. Using proxies for the radio emission we project the sources with different inclinations to the line of sight to the observer. We find that a number of different precession combinations result in characteristic `X' shaped sources which are frequently observed in radio data, and some precessing jet morphologies may mimic the morphological signatures of restarting radio sources. We look at jets ranging in scale from tens to hundreds of kiloparsecs and develop tools for identifying known precession indicators of point symmetry, curvature and jet misalignment from the lobe axis and show that, based on our simulation sample of precessing and non-precessing jets, a radio source that displays any of these indicators has a 98% chance of being a precessing source., Comment: Accepted for MNRAS. A link to simulations is available at https://www.extragalactic.info/precessingjets

Star clusters are fundamental units of stellar feedback and unique tracers of their host galactic properties. In this review, we will first focus on their constituents, i.e.\ detailed insight into their stellar populations and their surrounding ionised, warm, neutral, and molecular gas. We, then, move beyond the Local Group to review star cluster populations at various evolutionary stages, and in diverse galactic environmental conditions accessible in the local Universe. At high redshift, where conditions for cluster formation and evolution are more extreme, we are only able to observe the integrated light of a handful of objects that we believe will become globular clusters. We therefore discuss how numerical and analytical methods, informed by the observed properties of cluster populations in the local Universe, are used to develop sophisticated simulations potentially capable of disentangling the genetic map of galaxy formation and assembly that is carried by globular cluster populations., Comment: Final published version. 61 pages, 10 figures. To appear in Space Science Reviews (submitted January 31, 2020; accepted May 8, 2020), topical collection Star formation

Star formation is a complex multi-scale phenomenon that is of significant importance for astrophysics in general. Stars and star formation are key pillars in observational astronomy from local star forming regions in the Milky Way up to high-redshift galaxies. From a theoretical perspective, star formation and feedback processes (radiation, winds, and supernovae) play a pivotal role in advancing our understanding of the physical processes at work, both individually and of their interactions. In this review we will give an overview of the main processes that are important for the understanding of star formation. We start with an observationally motivated view on star formation from a global perspective and outline the general paradigm of the life-cycle of molecular clouds, in which star formation is the key process to close the cycle. After that we focus on the thermal and chemical aspects in star forming regions, discuss turbulence and magnetic fields as well as gravitational forces. Finally, we review the most important stellar feedback mechanisms., Comment: 86 pages, 4 figures. To appear in Space Science Reviews, topical collection Star formation

Star clusters form in dense, hierarchically collapsing gas clouds. Bulk kinetic energy is transformed to turbulence with stars forming from cores fed by filaments. In the most compact regions, stellar feedback is least effective in removing the gas and stars may form very efficiently. These are also the regions where, in high-mass clusters, ejecta from some kind of high-mass stars are effectively captured during the formation phase of some of the low mass stars and effectively channeled into the latter to form multiple populations. Star formation epochs in star clusters are generally set by gas flows that determine the abundance of gas in the cluster. We argue that there is likely only one star formation epoch after which clusters remain essentially clear of gas by cluster winds. Collisional dynamics is important in this phase leading to core collapse, expansion and eventual dispersion of every cluster. We review recent developments in the field with a focus on theoretical work., Comment: 53 pages, 7 figures, preprint of a review to appear in Space Science Reviews Topical collection "Star formation". Changes in this version: updated comments line to to acceptance of the paper by the journal

Feedback from radio jets associated with Active Galactic Nuclei (AGN) plays a profound role in the evolution of galaxies. Kinetic power of these radio jets appears to show temporal variation, but the mechanism(s) responsible for this process are not yet clear. Recently, the LOw Frequency ARray (LOFAR) has uncovered large populations of active, remnant and restarted radio jet populations. By focusing on LOFAR data in the Lockman Hole, in this work we use the Radio AGN in Semi-Analytic Environments (RAiSE) dynamical model to present the first self-consistent modelling analysis of active, remnant and restarted radio source populations. Consistent with other recent work, our models predict that remnant radio lobes fade quickly. Any high (>10 percent) observed fraction of remnant and restarted sources therefore requires a dominant population of short-lived jets. We speculate that this could plausibly be provided by feedback-regulated accretion., Comment: 14 pages, 7 figures, accepted for publication in MNRAS. Companion paper to Jurlin et al., also in today's arXiv posting

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., Comment: 9 pages, 8 figures, accepted for Monthly Notices

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., Comment: Accepted for publication in A&A

The annihilation of cosmic positrons ($e^+$) with electrons in the interstellar medium (ISM) results in the strongest persistent gamma-ray line signal in the sky. For 50 years, this 511 keV emission has puzzled observers and theoreticians. A key issue for understanding $e^+$-astrophysics is found in cosmic-ray propagation, especially at low kinetic energies (< 10 MeV). We want to shed light on how $e^+$s propagate and the resulting morphology of the 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 source populations. By dissecting the 511 keV emission as measured with INTEGRAL/SPI, we derive spectra for individual regions in the sky. The centroid energies are converted into Doppler-shifts, representing the line-of-sight velocity along different longitudes. This results in a longitude-velocity diagram of $e^+$-annihilation. We also determine Doppler-broadenings to study annihilation conditions as they vary across the Galaxy. We find line-of-sight velocities in the 511 keV line that are consistent with zero, as well as with galactic rotation from CO measurements, and measurements of radioactive Al-26. The velocity gradient in the inner 60 deg is determined to be $4\pm6$ km/s/deg. The 511 keV line width is constant as a function of longitude at $2.43\pm0.14$ keV. The positronium fraction is found to be 1.0 along the galactic plane. The weak signals in the disk leave open the question whether $e^+$-annihilation is associated with the high velocities seen in Al-26 or rather with ordinarily rotating components of the Galaxy's ISM. We confirm previous results that $e^+$s are slowed down to the 10 eV energy scale before annihilation, and constrain bulk Doppler-broadening contributions to <1.25 keV. Consequently, the true annihilation conditions remain unclear., Comment: 8 pages, 5 main figures, 7 pages appendix with 6 auxiliary figures, accepted for publication in A&A

Gaseous halos play a key role for 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 indeed identified a galaxy stellar mass scale, $10^{11} M_\odot$ , above which the radio luminosities increase disproportionately. Here, we investigate, if radio luminosities of galaxies, especially the marked rise at galaxy masses around $10^{11} M_\odot$, can be explained with standard assumptions on jet powers, scaling between black hole-mass and galaxy mass and gaseous halos. We developed models for the radio luminosity of radio AGN in halos under infall, galactic wind and hydrostatic conditions based on observational data and theoretical constraints, and compared it to LoTSS data for a large sample of galaxies in the mass rangebetween $10^{8.5} M_\odot$ and $10^{12} M_\odot$. Assuming the same characteristic upper limit to jet powers as is known from high galaxy masses to hold at all masses, we find that the maximum radio luminosities for the hydrostatic gas halos fit well with the upper envelope of the distribution of the LOFAR data. The marked rise in radio luminosity at $10^{11} M_\odot$ is matched in our model, and is related to significant change in halo gas density around this galaxy mass, which is a consequence of the lower cooling rates at higher virial temperature. Wind and infall models overpredict the radio luminosities at small galaxy masses and have no particular steepening of the run of the radio luminosities predicted at any galaxy mass. [...], Comment: 6 pages, 3 figures, accepted by A&A, comments welcome

The 511 keV line from positron annihilation in the Galaxy was the first $\gamma$-ray line detected to originate from outside our solar system. Going into the fifth decade since the discovery, the source of positrons is still unconfirmed and remains one of the enduring mysteries in $\gamma$-ray astronomy. With a large flux of $\sim$10$^{-3}$ $\gamma$/cm$^{2}$/s, after 15 years in operation INTEGRAL/SPI has detected the 511 keV line at $>50\sigma$ and has performed high-resolution spectral studies which conclude that Galactic positrons predominantly annihilate at low energies in warm phases of the interstellar medium. The results from imaging are less certain, but show a spatial distribution with a strong concentration in the center of the Galaxy. The observed emission from the Galactic disk has low surface brightness and the scale height is poorly constrained, therefore, the shear number of annihilating positrons in our Galaxy is still not well know. Positrons produced in $\beta^+$-decay of nucleosynthesis products, such as $^{26}$Al, can account for some of the annihilation emission in the disk, but the observed spatial distribution, in particular the excess in the Galactic bulge, remains difficult to explain. Additionally, one of the largest uncertainties in these studies is the unknown distance that positrons propagate before annihilation. In this paper, we will summarize the current knowledge base of Galactic positrons, and discuss how next-generation instruments could finally provide the answers., Comment: 10 pages including 3 figures. Science White paper submitted to Astro2020