Researcher profile

Martin Krause

Martin Krause contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
28works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2025arXiv

The origin of the spectral versus dynamical age discrepancy in radio galaxies

We investigate the effect of turbulent magnetic fields on the observed spectral properties of synchrotron radio emission in large-scale radio galaxy lobes. We use three-dimensional relativistic magnetohydrodynamic simulations of fast, high-powered jets to study the structure of the lobe magnetic fields and how this structure affects the radio spectrum of the lobes. It has previously been argued that lobe ages inferred from radio spectra underestimate the true ages of radio galaxies due to re-acceleration of electrons in the lobe, mixing of electron populations, or the presence of turbulent magnetic fields in the lobes. We find that the spectral ages with and without accounting for the lobe magnetic field structure are consistent with each other, suggesting that mixing of radiating populations of different ages is the primary cause of the underestimation of radio lobe ages. By accounting for the structure of lobe magnetic fields, we find greater spectral steepening in the equatorial regions of the lobe. We demonstrate that the assumptions of the continuous injection, Jaffe-Perola, and Tribble models for radio lobe spectra do not hold in our simulations, and we show that young particles with high magnetic field strengths are the dominant contributors to the overall radio lobe spectrum.

preprint2020arXiv

The duty cycle of radio galaxies revealed by LOFAR: remnant and restarted radio source populations in the Lockman Hole

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.

preprint2014arXiv

Are there any first-generation stars in globular clusters today?

Several models compete to explain the abundance properties of stellar populations in globular clusters. One of the main constraints is the present-day ratio of first- and second-generation stars that are currently identified based on their sodium content. We propose an alternative interpretation of the observed sodium distribution, and suggest that stars with low sodium abundance that are counted as members of the first stellar generation could actually be second-generation stars. We compute the number ratio of second-generation stars along the Na distribution following the fast rotating massive star model using the same constraints from the well-documented case of NGC 6752 as in our previous developments. We reproduce the typical percentage of low-sodium stars usually classified as first-generation stars by invoking only secondary star formation from material ejected by massive stars and mixed with original globular cluster material in proportions that account for the Li-Na anti-correlation in this cluster. Globular clusters could be totally devoid of first-generation low-mass stars today. This can be tested with the determination of the carbon isotopic ratio and nitrogen abundance in turn-off globular cluster stars. Consequences and related issues are briefly discussed.

preprint2014arXiv

Characterising superclusters with the galaxy cluster distribution

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.

preprint2014arXiv

Early 56Ni decay γ-rays from SN2014J suggest an unusual explosion

Type-Ia supernovae result from binary systems that include a carbon-oxygen white dwarf, and these thermonuclear explosions typically produce 0.5 M_solar of radioactive 56Ni. The 56Ni is commonly believed to be buried deeply in the expanding supernova cloud. Surprisingly, in SN2014J we detected the lines at 158 and 812 keV from 56Ni decay (τ~8.8 days) earlier than the expected several-week time scale, only ~20 days after the explosion, and with flux levels corresponding to roughly 10% of the total expected amount of 56Ni. Some mechanism must break the spherical symmetry of the supernova, and at the same time create a major amount of 56Ni at the outskirts. A plausible explanation is that a belt of helium from the companion star is accreted by the white dwarf, where this material explodes and then triggers the supernova event.

preprint2014arXiv

Feedback by massive stars and the emergence of superbubbles II. X-ray properties

In a previous paper we investigated the energy transfer of massive stars to the interstellar medium as a function of time and the geometrical configuration of three massive stars via 3D-mesh-refining hydrodynamics simulations, following the complete evolution of the massive stars and their supernovae except non-thermal processes . We analysed our ISM simulation results with the help of spectra for plasma temperatures between 0.1 and 10 keV and computed the spectral evolution and the spatio-temporal distribution of the hot gas. Results. Despite significant input of high temperature gas from supernovae and fast stellar winds, the resulting thermal X-ray spectra are generally very soft, with most of the emission well below 1 keV. We show that this is due to mixing triggered by resolved hydrodynamic instabilities. Supernovae enhance the X-ray luminosity of a superbubble by 1-2 orders of magnitude for a time span of about 0.1 Myr; longer if a supernova occurs in a larger superbubble and shorter in higher energy bands. Peak superbubble luminosities of the order of 10^{36} erg/s are reproduced well. The strong decay of the X-ray luminosity is due to bubble expansion, hydrodynamic instabilities related to the acceleration of the superbubble's shell thanks to the sudden energy input, and subsequent mixing. We also find global oscillations of our simulated superbubbles, which produce spatial variations of the X-ray spectrum, similar to what we see in the Orion-Eridanus cavity. We calculated the fraction of energy emitted in X-rays and find that with a value of a few times 10^{-4}, it is about a factor of ten below the measurements for nearby galaxies.

preprint2014arXiv

SN2014J gamma-rays from the 56Ni decay chain

The measurement of gamma-ray lines from the decay chain of 56Ni provides unique information about the explosion in supernovae. The 56Ni freshly-produced in the supernova powers the optical light curve, as it emits gamma-rays upon its radioactive decay to 56Co and then 56Fe. Gamma-ray lines from 56Co decay are expected to become directly visible through the overlying white dwarf material several weeks after the explosion, as they progressively penetrate the overlying material of the supernova envelope, diluted as it expands. The lines are expected to be Doppler-shifted or broadened from the kinematics of the 56Ni ejecta. With the SPI spectrometer on INTEGRAL and using an improved instrumental background method, we detect the two main lines from 56Co decay at 847 and 1238 keV from SN2014J at 3.3 Mpc, significantly Doppler-broadened, and at intensities (3.65+/-1.21) 10^-4 and (2.27+/-0.69) 10^-4 ph cm^-2s^-1, respectively, at brightness maximum. We measure their rise towards a maximum after about 60-100 days and decline thereafter. The intensity ratio of the two lines is found consistent with expectations from 56Co decay (0.62+/-0.28 at brightness maximum, expected is 0.68). We find that the broad lines seen in the late, gamma-ray transparent phase are not representative for the early gamma-ray emission, and rather notice the emission spectrum to be complex and irregular until the supernova is fully transparent to gamma-rays, with progressive uncovering of the bulk of 56Ni. We infer that the explosion morphology is not spherically symmetric, both in the distribution of 56Ni and of the unburnt material which occults the 56Co emission. Comparing light curves from different plausible models, the resulting 56Ni mass is determined as 0.49+/-0.09 M_sol.

preprint2013arXiv

Galactic winds - How to launch galactic outflows in typical Lyman-break galaxies

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.

preprint2013arXiv

Kinematics of massive star ejecta in the Milky Way as traced by 26Al

Massive stars form in groups and their winds and supernova explosions create superbubbles up to kpc in size. Their ejecta are important for the dynamics of the interstellar medium and chemical evolution models. However, ejecta kinematics and the characteristic scales in space and time are rather unexplored beyond ~10 ka. Through measurement of radioactive 26Al with its decay time constant of ~10^6 years, we aim to trace the kinematics of massive-star and supernova ejecta over million-year time scales. Our goal is to identify the mixing time scale and the spatio-kinematics of such ejecta from the pc to kpc scale. We use the SPI spectrometer on INTEGRAL and its observations along the Galactic ridge to trace the line shape of the 1808.63 keV gamma-ray line from 26Al decay. We determine line centroids and compare these to Doppler shift expectations from large-scale systematic rotation around the Galaxy's center, as observed in other Galactic objects. We measure the radial velocities of gas traced by 26Al, averaged over the line of sight, as a function of Galactic longitude. We find substantially higher velocities than expected from Galactic rotation, the average bulk velocity being ~200 km s^-1 larger than the Galactic-rotation prediction. The observed radial velocity spread implies a Doppler-broadening of the gamma-ray line that is consistent with our measurements of the overall line width. We can reproduce the observed characteristics with 26Al sources located along the inner spiral arms, when we add a global blow-out into the forward direction away from arms, such as expected from massive stars offset towards the spiral-arm leading edge. The superbubbles-halo connection implies angular-momentum transfer in the disk-halo system and also radial gas flows.

preprint2013arXiv

SPI Measurements of Nucleosynthesis Gamma-Rays

Studies based on the gamma-ray lines from radioactive decay of unstable isotopes produced in massive-star and supernova nucleosynthesis have been among INTEGRAL's prominent science achievements. 26Al has become a tool to study specific source regions, such as massive-star groups and associations in nearby regions which can be discriminated from the galactic-plane background, and the inner Galaxy where Doppler shifted lines add to the astronomical information. 60Fe is co-produced by the sources of 26Al, and the isotopic ratio from their nucleosynthesis encodes stellar-structure information. Here we summarize latest results using the accumulated multi-year database of observations, point to the relevant publications, and discuss their astrophysical implications.

preprint2013arXiv

Superbubble dynamics in globular cluster infancy II. Consequences for secondary star formation in the context of self-enrichment via fast rotating massive stars

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. [...]

preprint2012arXiv

A new connection between the opening angle and the large-scale morphology of extragalactic radio sources

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.

preprint2012arXiv

Discovery of an X-ray cavity near the radio lobes of Cygnus A indicating previous AGN activity

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.

preprint2012arXiv

Feedback by massive stars and the emergence of superbubbles I. Energy efficiency & Vishniac instabilities

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.

preprint2012arXiv

Jet interactions with a giant molecular cloud in the Galactic centre and ejection of hypervelocity stars

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.

preprint2012arXiv

Jet-induced star formation in gas-rich galaxies

Feedback from active galactic nuclei (AGN) has become a major component in simulations of galaxy evolution, in particular for massive galaxies. AGN jets have been shown to provide a large amount of energy and are capable of quenching cooling flows. Their impact on the host galaxy, however, is still not understood. Subgrid models of AGN activity in a galaxy evolution context so far have been mostly focused on the quenching of star formation. To shed more light on the actual physics of the "radio mode" part of AGN activity, we have performed simulations of the interaction of a powerful AGN jet with the massive gaseous disc (10^11 solar masses) of a high-redshift galaxy. We spatially resolve both the jet and the clumpy, multi-phase interstellar medium (ISM) and include an explicit star formation model in the simulation. Following the system over more than 10^7 years, we find that the jet activity excavates the central region, but overall causes a significant change to the shape of the density probability distribution function and hence the star formation rate due to the formation of a blast wave with strong compression and cooling in the ISM. This results in a ring- or disc-shaped population of young stars. At later times, the increase in star formation rate also occurs in the disc regions further out since the jet cocoon pressurizes the ISM. The total mass of the additionally formed stars may be up to 10^10 solar masses for one duty cycle. We discuss the details of this jet-induced star formation (positive feedback) and its potential consequences for galaxy evolution and observable signatures.

preprint2012arXiv

Magnetohydrodynamic stability of broad line region clouds

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.

preprint2012arXiv

Superbubble dynamics in globular cluster infancy I. How do globular clusters first lose their cold gas?

The picture of the early evolution of globular clusters has been significantly revised in recent years. Current scenarios require at least two generations of stars of which the first generation (1G), and therefore also the protocluster cloud, has been much more massive than the currently predominating second generation (2G). Fast gas expulsion is thought to unbind the majority of the 1G stars. Gas expulsion is also mandatory to remove metal-enriched supernova ejecta, which are not found in the 2G stars. It has long been thought that the supernovae themselves are the agent of the gas expulsion, based on crude energetics arguments. Here, we assume that gas expulsion happens via the formation of a superbubble, and describe the kinematics by a thin-shell model. We find that supernova- driven shells are destroyed by the Rayleigh-Taylor instability before they reach escape speed for all but perhaps the least massive and most extended clusters. More power is required to expel the gas, which might plausibly be provided by a coherent onset of accretion onto the stellar remnants. The resulting kpc-sized bubbles might be observable in Faraday rotation maps with the planned Square Kilometre Array radio telescope against polarised background radio lobes if a globular cluster would happen to form in front of such a radio lobe.

preprint2011arXiv

Interaction of Fanaroff-Riley class II jets with a magnetised intra-cluster medium

We present 3-D MHD and synthetic numerical simulations to follow the evolution of randomly magnetized intra-cluster medium plasma under the effects of powerful, light, hypersonic and bipolar jets. We prescribe the cluster magnetic field (CMF) as a Gaussian random field with power law energy spectrum tuned to the expectation for Kolmogorov turbulence. We investigate the power of jets and the viewing angle used for the synthetic Rotation Measure (RM) observations. We find the model radio sources introduce and amplify fluctuations on the RM statistical properties; the average RM and the RM standard deviation are increased by the action of the jets. This may lead to overestimations of the CMFs' strength up to 70%. The effect correlates with the jet power. Jets distort and amplify CMFs especially near the edges of the lobes and the jets' heads. Thus the RM structure functions are flattened at scales comparable to the source size. Jet-produced RM enhancements depend on the orientation of the jet axis to the line of sight and are thus more apparent in quasars than in radio galaxies.

preprint2011arXiv

Interaction of Fanaroff-Riley class II radio jets with a randomly magnetised intra-cluster medium

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.

preprint2011arXiv

Radiation feedback on dusty clouds during Seyfert activity

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.

preprint2011arXiv

Three-dimensional magnetohydrodynamic simulations of the evolution of magnetic fields in Fanaroff-Riley class II radio sources

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.

preprint2010arXiv

3D-MHD simulations of the evolution of magnetic fields in FR II radio sources

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.

preprint2010arXiv

Asymmetries in Extragalactic Double Radio Sources: Clues from 3D Simulations of Jet - Disc Interaction

Observational and theoretical studies of extragalactic radio sources have suggested that an inhomogeneous environment may be responsible for observed arm length asymmetries of jets and the properties of extended emission line regions in high redshift radio galaxies. We perform 3D hydrodynamic simulations of the interaction of a powerful extragalactic bipolar jet with a disc-shaped clumpy interstellar medium of log-normal density distribution and analyze the asymmetry. Furthermore, we compute the relation between jet asymmetry and the ISM properties by means of Monte Carlo simulations based on a 1D propagation model for the jet through the dense medium. We find that the properties of the ISM can be related to a probability distribution of jet arm length asymmetries: Disc density and height are found to have the largest effect on the asymmetry for realistic parameter ranges, while the Fourier energy spectrum of the ISM and turbulent Mach number only have a smaller effect. The hydrodynamic simulations show that asymmetries generally may be even larger than expected from the 1D model due to the complex interaction of the jet and its bow shock with gaseous clumps, which goes much beyond simple energy disposal. From our results, observed asymmetries of medium-sized local radio galaxies may be explained by gas masses of 10^9 to 10^10 solar masses in massive elliptical galaxies. Furthermore, the simulations provide a theoretical basis for the observed correlation that emission line nebulae are generally found to be brighter on the side of the shorter lobe in high redshift radio galaxies (McCarthy et al. 1991). This interaction of jets with the cold gas phase suggests that star formation in evolving high redshift galaxies may be affected considerably by jet activity.

preprint2010arXiv

Stability of Cloud Orbits in the Broad Line Region of Active Galactic Nuclei

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.

preprint2009arXiv

Measurements of the Cosmological Evolution of Magnetic Fields with the Square Kilometre Array

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.

preprint2009arXiv

Physics and fate of jet related emission line regions

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.