Search

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

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

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

Early on, extragalactic radio sources have pointed to a cosmologically evolving Universe. They were also an important piece of evidence for the existence of supermassive black holes, now thought to be a key component of galaxies. The observation that the power of radio sources increases with redshift, whereas the cosmological assembly of mass proceeds vice versa means that radio sources have their strongest impact in the early Universe. Our simulations suggest that radio sources heat hot halo gas, boost star formation in disc galaxies and other cold gas in the vicinity, possibly filaments, by a surround and squash mechanism. They might cause gaseous outflows in connection with stellar feedback. This might be an important mode of star formation for forming massive galaxies. Analysis of the jet-environment interaction may provide insights into black-hole physics and jet formation, e.g., rotational energy extraction (Blandford-Znajek) or how frequent black-hole binaries or multiple systems are. The former relates to fundamental questions about the nature of black holes. The latter is expected from hierarchical cosmology. Extragalactic radio sources thus continue to corroborate the cosmological picture and lead the way towards new, exciting discoveries., Comment: Invited review given at the international conference "The Physical Universe" in Nagpur, India, 26 Feb - 1 March 2018. Minor style adaptation and a clarification in Sect. 4.2 compared to the version printed in the conference book

Superclusters are the largest, observed matter density structures in the Universe. Recently Chon et al.(2013) presented the first supercluster catalogue constructed with a well-defined selection function based on the X-ray flux-limited cluster survey, REFLEX II. For the construction of the sample we proposed a concept to find the large objects with a minimum overdensity such that most of their mass will collapse in the future. The main goal of the paper is to provide support for our concept using simulations that we can, on the basis of our observational sample of X-ray clusters, construct a supercluster sample defined by a certain minimum overdensity, and to test how superclusters trace the underlying dark matter distribution. Our results confirm that an overdensity in the number of clusters is tightly correlated with an overdensity of the dark matter distribution. This enables us to define superclusters such that most of the mass will collapse in the future and to get first-order mass estimates of superclusters on the basis of the properties of the member clusters. We also show that in this context the ratio of the cluster number density and dark matter mass density is consistent with the theoretically expected cluster bias. Our previous work provided evidence that superclusters are a special environment for density structures of the dark matter to grow differently from the field as characterised by the X-ray luminosity function. Here we confirm for the first time that this originates from a top-heavy mass function at high statistical significance provided by a Kolmogorov-Smirnov test. We also find in close agreement with observations that the superclusters occupy only a small volume of few percent while they contain more than half of the clusters in the present day Universe., Comment: 11 pages, 11 figures, accepted for publication in Astronomy & Astrophysics

We perform hydrodynamical simulations of a young galactic disc embedded in a hot gaseous halo using parameters typical for Lyman break galaxies (LBGs). We take into account the (static) gravitational potentials due to a dark matter halo, a stellar bulge and a disc of stars and gas. Star formation is treated by a local Kennicutt-Schmidt law. We simplify the structure of the interstellar medium by restricting the computational domain to a 25th of the full azimuthal angle, effectively assuming large-scale axisymmetry and neglecting any effects of spiral structure, and focus on the large-scale ISM drivers, the superbubbles. Supernovae are triggered randomly and have preset event sizes of several tens to hundreds. We further investigate different halo gas pressures and energy injection methods. Many of our simulated galaxies, but not all, develop bipolar outflows. We characterise the strength of the outflow by mass and energy outflow rates, and investigate the effect of changes to the details of the model. We find that supernovae are more effective if comprised into larger superbubbles. The weight and the pressure of the halo gas is able to quench galactic outflows. The wind emerges from a series of superbubbles in regions where a critical star formation density is exceeded. The superbubbles expand into the gaseous halo at slightly supersonic speed, producing radiative shock waves with similar characteristics as the absorptions systems observed around LBGs., Comment: 20 pages, 21 figures, accepted by MNRAS

The self-enrichment scenario for globular clusters (GC) requires large amounts of residual gas after the initial formation of the first stellar generation. Recently, Krause et al. (2012) found that supernovae may not be able to expel that gas, as required to explain their present day gas-free state, and suggested that a sudden accretion on to the dark remnants, at a stage when type II supernovae have ceased, may plausibly lead to fast gas expulsion. Here, we explore the consequences of these results for the self-enrichment scenario via fast rotating massive stars (FRMS). We analyse the interaction of FRMS with the intra-cluster medium (ICM), in particular where, when and how the second generation of stars may form. From the results, we develop a timeline of the first approximately 40 Myr of GC evolution. The results of Paper I imply three phases during which the ICM is in a fundamentally different state, namely the wind bubble phase (lasting 3.5 to 8.8 Myr), the supernova phase (lasting 26.2 to 31.5 Myr), and the dark remnant accretion phase (lasting 0.1 to 4 Myr): (i) Quickly after the first generation massive stars have formed, stellar wind bubbles compress the ICM into thin filaments. No stars may form in the normal way during this phase, due to the high Lyman-Werner flux density. If the first generation massive stars have however equatorial ejections, as we proposed in the FRMS scenario, accretion may resume in the shadow of the equatorial ejecta. The second generation stars may then form due to gravitational instability in these discs that are fed by both the FRMS ejecta and pristine gas. [...], Comment: 15 pages, 6 figures, A&A accepted

The hypervelocity OB stars in the Milky Way Galaxy were ejected from the central regions some 10-100 million years ago. We argue that these stars, {as well as many more abundant bound OB stars in the innermost few parsecs,} were generated by the interactions of an AGN jet from the central black hole with a dense molecular cloud. Considerations of the associated energy and momentum injection have broader implications for the possible origin of the Fermi bubbles and for the enrichment of the intergalactic medium., Comment: 4 pages, 1 figure. Astronomy and Astrophysics Letters, in press

We study the evolution of isolated and merging interstellar bubbles of three stars (25, 32 and 60 M\odot) in a homogeneous background medium with a density of 10 mp /ccm via 3D-hydrodynamic simulations with standard ISM thermodynamics (optically thin radiative cooling and photo-electric heating) and time dependent energy and mass input according to stellar evolutionary tracks. We vary the position of the three stars relative to each other to compare the energy response for cases of isolated, merging and initially cospatial bubbles. Due to mainly the Vishniac instability, our simulated bubbles develop thick shells and filamentary internal structures in column density. The shell widths reach tens of per cent of the outer bubble radius, which compares favourably to observations. More energy is retained in the ISM for more closely packed groups, by up to a factor of three and typically a factor of two for intermediate times after the first supernova. Once the superbubble is established, different positions of the contained stars make only a minor difference to the energy tracks. For our case of three massive stars, the energy deposition varies only very little for distances up to about 30 pc between the stars. Energy injected by supernovae is entirely dissipated in a superbubble on a timescale of about 1 Myr, which increases slightly with the superbubble size at the time of the explosion. The Vishniac instability may be responsible for the broadening of the shells of interstellar bubbles. Massive star winds are significant energetically due to their - in the long run - more efficient, steady energy injection and because they evacuate the space around the massive stars. For larger scale simulations, the feedback effect of close groups of stars or clusters may be subsumed into one effective energy input with insignificant loss of energy accuracy., Comment: 14 pages, 15 figures, accepted by A&A, shortened, different code test

Hydrodynamic stability has been a longstanding issue for the cloud model of the broad line region in active galactic nuclei. We argue that the clouds may be gravitationally bound to the supermassive black hole. If true, stabilisation by thermal pressure alone becomes even more difficult. We further argue that if magnetic fields should be present in such clouds at a level that could affect the stability properties, they need to be strong enough to compete with the radiation pressure on the cloud. This would imply magnetic field values of a few Gauss for a sample of Active Galactic Nuclei we draw from the literature. We then investigate the effect of several magnetic configurations on cloud stability in axi-symmetric magnetohydrodynamic simulations. For a purely azimuthal magnetic field which provides the dominant pressure support, the cloud first gets compressed by the opposing radiative and gravitational forces. The pressure inside the cloud then increases, and it expands vertically. Kelvin-Helmholtz and column density instability lead to a filamentary fragmentation of the cloud. This radiative dispersion continues until the cloud is shredded down to the resolution level. For a helical magnetic field configuration, a much more stable cloud core survives with a stationary density histogram which takes the form of a power law. Our simulated clouds develop sub-Alfvenic internal motions on the level of a few hundred km/s., Comment: 16 pages, 11 figures, accepted by MNRAS, Figure 6 updated due to inconsistent linestyles, corrected mistake in Alfven speed formula, some very minor language corrections. The definitive version is available at www.blackwell- synergy.com

Cygnus A harbours the nearest powerful radio jet of an Fanaroff-Riley (FR) class II radio galaxy in a galaxy cluster where the interaction of the jet with the intracluster medium (ICM) can be studied in detail. We use a large set of Chandra archival data, VLA and new LOFAR observations to shed new light on the interaction of the jets with the ICM. We identify an X-ray cavity in the distribution of the X-ray emitting plasma in the region south of the Cyg A nucleus which has lower pressure than the surrounding medium. The LOFAR and VLA radio observations show that the cavity is filled with synchrotron emitting plasma. The spectral age and the buoyancy time of the cavity indicates an age at least as large as the current Cyg A jets and not much larger than twice this time. We suggest that this cavity was created in a previous active phase of Cyg A when the energy output of the Active Galactic Nucleus (AGN) was about two orders of magnitude less than today., Comment: Letter submitted on 4 May 2012 to A&A, 4 pages, 4 figures

In the case of an initially conical jet, we study the relation between jet collimation by the external pressure and large-scale morphology. We first consider the important length-scales in the problem, and then carry out axisymmetric hydrodynamic simulations that include, for certain parameters, all these length-scales. We find three important scales related to the collimation region: (i) where the sideways ram-pressure equals the external pressure, (ii) where the jet density equals the ambient density, and (iii) where the forward ram-pressure falls below the ambient pressure. These scales are set by the external Mach-number and opening angle of the jet. We demonstrate that the relative magnitudes of these scales determine the collimation, Mach-number, density and morphology of the large scale jet. Based on analysis of the shock structure, we reproduce successfully the morphology of Fanaroff-Riley (FR) class I and II radio sources. Within the framework of the model, an FR I radio source must have a large intrinsic opening angle. Entrainment of ambient gas might also be important. We also show that all FR I sources with radio lobes or similar features must have had an earlier FR II phase., Comment: 14 pages, 10 figures, accepted by MNRAS, same as previous version

A combination of three-dimensional (3D) magnetohydrodynamics (MHD) and synthetic numerical simulations are presented to follow the evolution of a randomly magnetised plasma that models the intra-cluster medium (ICM), under the isolated effects of powerful, light, hypersonic and bipolar Fanaroff-Riley class II (FR II) jets. We prescribe the cluster magnetic field (CMF) as a Gaussian random field with a Kolmogorov-like energy spectrum. Both the power of the jets and the viewing angle that is used for the synthetic Rotation Measure (RM) observations are investigated. We find the model radio sources introduce and amplify fluctuations on the RM statistical properties which we analyse as a function of time as well as the viewing angle. The average RM and the RM standard deviation are increased by the action of the jets. Energetics, RM statistics and magnetic power spectral analysis consistently show that the effects also correlate with the jets' power, and that the lightest, fastest jets produce the strongest changes in their environment. We see jets distort and amplify the CMFs especially near the edges of the lobes and the jets' heads. This process leads to a flattening of the RM structure functions at scales comparable to the source size. The edge features we find are similar to ones observed in Hydra A. The results show that jet-produced RM enhancements are more apparent in quasars than in radio galaxies. Globally, jets tend to enhance the RM standard deviation which may lead to overestimations of the CMFs' strength by about 70%. This study means to serve as a pathfinder for the SKA, EVLA and LOFAR to follow the evolution of cosmic magnetic fields., Comment: Accepted for publication in the MNRAS. 21 pages, 15 figures

Radio observations of Fanaroff-Riley class II sources often show correlations between the synchrotron emission and the linear-polarimetric distributions. Magnetic position vectors seem to align with the projected emission of both the radio jets and the sources' edges. Using statistics we study such relation as well as its unknown time evolution via synthetic polarisation maps of model FR II sources formed in 3D-MHD numerical simulations of bipolar, hypersonic and weakly magnetised jets. The magnetic field is initially random with a Kolmogorov power spectrum, everywhere. We investigate the structure and evolution of magnetic fields in the sources as a function of the power of jets and the observational viewing angle. Our synthetic polarisation maps agree with observations, showing B-field vectors which are predominantly aligned with the jet axis, and show that magnetic fields inside sources are shaped by the jets' backflow. Polarimetry is found to correlate with time, the viewing angle and the jet-to-ambient density contrast. The magnetic structure inside thin elongated sources is more uniform than inside more spherical ones. We see jets increase the magnetic energy in cocoons in proportion to the jet velocity and the cocoon width. Filaments in the synthetic emission maps suggest turbulence develops in evolved sources., Comment: Accepted for publication in the MNRAS. 21 pages, 11 figures

We investigate the evolution of dusty gas clouds falling into the centre of an active Seyfert nucleus. Two-dimensional high-resolution radiation hydrodynamics simulations are performed to study the fate of single clouds and the interaction between two clouds approaching the Active Galactic Nucleus. We find three distinct phases of the evolution of the cloud: (i) formation of a lenticular shape with dense inner rim caused by the interaction of gravity and radiation pressure (the lense phase), (ii) formation of a clumpy sickle-shaped structure as the result of a converging flow (the clumpy sickle phase) and (iii) a filamentary phase caused by a rapidly varying optical depth along the sickle. Depending on the column density of the cloud, it will either be pushed outwards or its central (highest column density) parts move inwards, while there is always some material pushed outwards by radiation pressure effects. The general dynamical evolution of the cloud can approximately be described by a simple analytical model., Comment: 13 pages, 18 figures, accepted by MNRAS

3D-MHD numerical simulations of bipolar, hypersonic, weakly magnetized jets and synthetic synchrotron observations are presented to study the structure and evolution of magnetic fields in FR II radio sources. The magnetic field setup in the jet is initially random. The power of the jets as well as the observational viewing angle are investigated. We find that synthetic polarization maps agree with observations and show that magnetic fields inside the sources are shaped by the jets' backflow. Polarimetry statistics correlates with time, the viewing angle and the jet-to-ambient density contrast. The magnetic structure inside thin elongated sources is more uniform than for ones with fatter cocoons. Jets increase the magnetic energy in cocoons, in proportion to the jet velocity. Both, filaments in synthetic emission maps and 3D magnetic power spectra suggest that turbulence develops in evolved sources., Comment: 2 pages, 1 figure, accepted for publication in the Proceedings of IAU Symposium 275, "Jets at all Scales", Buenos Aires, Argentina, 13-17 September 2010, eds. G.E. Romero, R.A. Sunyaev and T. Belloni

We investigate the global dynamic stability of spherical clouds in the Broad Line Region (BLR) of Active Galactic Nuclei (AGN), exposed to radial radiation pressure, gravity of the central black hole (BH), and centrifugal forces assuming the clouds adapt their size according to the local pressure. We consider both, isotropic and anisotropic light sources. In both cases, stable orbits exist also for very sub-Keplerian rotation for which the radiation pressure contributes substantially to the force budget. We demonstrate that highly eccentric, very sub-Keplerian stable orbits may be found. This gives further support for the model of Marconi et al. 2008, who pointed out that black hole masses might be significantly underestimated if radiation pressure is neglected. That model improved the agreement between black hole masses derived in certain active galaxies based on BLR dynamics, and black hole masses derived by other means in other galaxies by inclusion of a luminosity dependent term. For anisotropic illumination, energy is conserved for averages over long time intervals, only, but not for individual orbits. This leads to Rosetta orbits that are systematically less extended in the direction of maximum radiation force. Initially isotropic relatively low column density systems would therefore turn into a disk when an anisotropic AGN is switched on., Comment: 7 pages, 4 fig.s, accepted by MNRAS, some text was missing - corrected

We investigate the potential of the Square Kilometre Array (SKA) for measuring the magnetic fields in clusters of galaxies via Faraday rotation of background polarised sources. [...] We find that about 10 per cent of the sky is covered by a significant extragalactic Faraday screen. Most of it has rotation measures between 10 and 100 rad/m/m. We argue that the cluster centres should have up to about 5000 rad/m/m. We show that the proposed mid frequency aperture array of the SKA as well as the lowest band of the SKA dish array are well suited to make measurements for most of these rotation measure values, typically requiring a signal-to-noise of ten. We calculate the spacing of sources forming a grid for the purpose of measuring foreground rotation measures: it reaches a spacing of 36 arcsec for a 100 hour SKA observation per field. We also calculate the statistics for background RM measurements in clusters of galaxies. We find that a first phase of the SKA would allow us to take stacking experiments out to high redshifts (>1), and provide improved magnetic field structure measurements for individual nearby clusters. The full SKA aperture array would be able to make very detailed magnetic field structure measurements of clusters with more than 100 background sources per cluster up to a redshift of 0.5 and more than 1000 background sources per cluster for nearby clusters, and could for reasonable assumptions about future measurements of electron densities in high redshift clusters constrain the power law index for the magnetic field evolution to better than dm=0.4, if the magnetic field in clusters should follow B ~ (1+z)^m., Comment: 12 pages, 10 figures, 3 tables, accepted by MNRAS, minor correction to eq (5)

At high redshift, extragalactic jets are associated with extended emission line regions. We present global and local hydrodynamic simulations of the interaction of jets with their environment relevant to this phenomenon. We determine the fraction of kinetic energy and momentum of the jet that should appear in the emission line gas. From momentum considerations, we argue that the gas seen in outward motion in many examples must to a large part originate from within the galaxy. The host galaxies should have been massive, gas rich and heavily star forming galaxies. The expelled gas is not expected to return to the galaxy quickly. The observed kinetic energy in the emission line gas constrains the jet powers to be greater than 10^47 erg/s., Comment: 12 pages, 5 figures