Search

Self-enrichment is one of the leading explanations for chemical anomalies in globular clusters. In this scenario, various candidate polluter stars have been proposed to eject gas with altered chemical composition during the self-enrichment process. Most of the proposed polluters will also eject radioactive $^{26}$Al into the surroundings. Hence, any detection of $^{26}$Al in young massive star clusters (YMCs) would support the self-enrichment scenario if YMCs were indeed the progenitors of globular clusters. Observations of gamma-ray data from COMPTEL and INTEGRAL, as well as detections of $^{26}$AlF molecules by the Atacama Large Millimeter-submillimeter Array (ALMA), indicate the maturing of $^{26}$Al detection methods. Detection possibilities will be enhanced in the short- to mid-term by the upcoming launch of the Compton Spectrometer and Imager (COSI). The Square Kilometer Array (SKA) could in principle also detect radio recombination lines of the positronium formed from the decay products of $^{26}$Al. Here, we show for a sample of YMCs in the nearby Universe, where self-enrichment could plausibly take place. For some nearby galaxies, this could enhance $^{26}$Al by an order of one magnitude. Detecting $^{26}$AlF with ALMA appears feasible for many candidate self-enrichment clusters, although significant challenges remain with other detection methods. The Large Magellanic Cloud, with its YMC R136, stands out as the most promising candidate. Detecting a 1.8~MeV radioactive decay line of $^{26}$Al here would require at least 15 months of targeted observation with COSI, assuming ongoing self-enrichment in R136., Comment: 13 pages, 5 figures

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