Magnetized Interstellar Medium Research Articles

Why only a small fraction of quasars are radio loud?

AbstractIt is still a mystery why only a small fraction of quasars contain relativistic jets. A strong magnetic field is a necessary ingredient for jet formation. Gas falls from the Bondi radius RB nearly freely to the circularization radius Rc, and a thin accretion disk is formed within Rc We suggest that the external weak magnetic field threading interstellar medium is substantially enhanced in this region, and the magnetic field at Rc can be sufficiently strong to drive outflows from the disk if the angular velocity of the gas is low at RB. In this case, the magnetic field is efficiently dragged in the disk, because most angular momentum of the disk is removed by the outflows that lead to a significantly high radial velocity. The strong magnetic field formed in this way may accelerate jets in the region near the black hole, either by the Blandford-Payne or/and Blandford-Znajek mechanisms. If the angular velocity of the circumnuclear gas is low, the field advection in the thin disk is inefficient, and it will appear as a radio-quiet (RQ) quasar.

Galactic neutral hydrogen and the magnetic ISM foreground

AbstractThe interstellar medium is suffused with magnetic fields, which inform the shape of structures in the diffuse gas. Recent high-dynamic range observations of Galactic neutral hydrogen, combined with novel data analysis techniques, have revealed a deep link between the morphology of neutral gas and the ambient magnetic field. At the same time, an observational revolution is underway in low-frequency radio polarimetry, driven in part by the need to characterize foregrounds to the cosmological 21-cm signal. A new generation of experiments, capable of high angular and Faraday depth resolution, are revealing complex filamentary structures in diffuse polarization. The relationship between filamentary structures observed in radio-polarimetric data and those observed in atomic hydrogen is not yet well understood. Multiwavelength observations will enable new insights into the magnetic interstellar medium across phases.

Structure distribution and turbulence in self-consistently supernova-driven ISM of multiphase magnetized galactic discs

Context. Galaxy evolution and star formation are two multi-scale problems tightly linked to each other. Aims. We aim to describe simultaneously the large-scale evolution widely induced by the feedback processes and the details of the gas dynamics that controls the star formation process through gravitational collapse. This is a necessary step in understanding the interstellar cycle, which triggers galaxy evolution. Methods. We performed a set of three-dimensional high-resolution numerical simulations of a turbulent, self-gravitating and magnetized interstellar medium within a 1 kpc stratified box with supernova feedback correlated with star-forming regions. In particular, we focussed on the role played by the magnetic field and the feedback on the galactic vertical structure, the star formation rate (SFR) and the flow dynamics. For this purpose we have varied their respective intensities. We extracted properties of the dense clouds arising from the turbulent motions and compute power spectra of various quantities. Results. Using a distribution of supernovae sufficiently correlated with the dense gas, we find that supernova explosions can reproduce the observed SFR, particularly if the magnetic field is on the order of a few μG. The vertical structure, which results from a dynamical and an energy equilibrium is well reproduced by a simple analytical model, which allows us to roughly estimate the efficiency of the supernovae in driving the turbulence in the disc to be rather low, of the order of 1.5%. Strong magnetic fields may help to increase this efficiency by a factor of between two and three. To characterize the flow we compute the power spectra of various quantities in 3D but also in 2D in order to account for the stratification of the galactic disc. We find that within our setup, the compressive modes tend to dominate in the equatorial plane, while at about one scale height above it, solenoidal modes become dominant. We measured the angle between the magnetic and velocity fields and we conclude that they tend to be well aligned particularly at high magnetization and lower feedback. Finally, the dense structures present scaling relations that are reminiscent of the observational ones. The virial parameter is typically larger than 10 and shows a large spread of masses below 1000 M⊙. For masses larger than 104M⊙, its value tends to a few. Conclusions. Using a relatively simple scheme for the supernova feedback, which is self-consistently proportional to the SFR and spatially correlated to the star formation process, we reproduce a stratified galactic disc that presents reasonable scale height, SFR as well as a cloud distribution with characteristics close to the observed ones.

Modelling and simulation of large-scale polarized dust emission over the southern Galactic cap using the GASS Hi data

The Planck survey has quantified polarized Galactic foregrounds and established that they are a main limiting factor in the quest for the cosmic microwave background B-mode signal induced by primordial gravitational waves during cosmic inflation. Accurate separation of the Galactic foregrounds therefore binds this quest to our understanding of the magnetized interstellar medium. The two most relevant empirical results from analysis of Planck data are line of sight depolarization arising from fluctuations of the Galactic magnetic field orientation and alignment of filamentary dust structures with the magnetic field at high Galactic latitude. Furthermore, Planck and H I emission data in combination indicate that most of the filamentary dust structures are in the cold neutral medium. The goal of this paper is to test whether these salient observational results, taken together, can account fully for the statistical properties of the dust polarization over a selected low column density region comprising 34% of the southern Galactic cap (b ≤ −30°). To do this, we construct a dust model that incorporates H I column density maps as tracers of the dust intensity structures and a phenomenological description of the Galactic magnetic field. By adjusting the parameters of the dust model, we were able to reproduce the Planck dust observations at 353GHz in the selected region. Realistic simulations of the polarized dust emission enabled by such a dust model are useful for testing the accuracy of component separation methods, studying non-Gaussianity, and constraining the amount of decorrelation with frequency.

Dust-polarization Maps and Interstellar Turbulence

Abstract Perhaps the most intriguing result of Planck’s dust-polarization measurements is the observation that the power in the E-mode polarization is twice that in the B mode, as opposed to pre-Planck expectations of roughly equal dust powers in the E and B modes. Here we show how the E- and B-mode powers depend on the detailed properties of the fluctuations in the magnetized interstellar medium (ISM). These fluctuations can be decomposed into slow, fast, and Alfvén magnetohydrodynamic (MHD) waves, which comprise a complete basis that can be used to describe linear fluctuations of a magnetized fluid. They can alternatively be decomposed in terms of one longitudinal and two transverse components of a fluid-displacement field. The intensity (T) and E- and B-mode amplitudes induced by each of the three types of waves, in both decompositions, are then calculated. To illustrate how these tools can be applied, we consider a toy model of the ISM in which dust traces a single component of plasma, and obtain the EE/BB ratio and TE correlation for several ansatzes for the power in slow/fast/Alfvén waves and in longitudinal/transverse waves. Although our model may be too simplistic to properly describe the nonlinear structure of interstellar magnetic fields, we find that the observed EE/BB ratio (and its scale invariance) and positive TE correlation—as well as the observed power-law index for the angular spectrum of these fluctuations—are not easily accommodated in terms of simple MHD turbulence prescriptions for the expected powers in slow, fast, and Alfvén waves. We speculate that the ∼0.1–30 pc length scales probed by these dust-polarization measurements are not described by MHD turbulence, but rather probe the large-scale physics that drives ISM turbulence. We find that a slightly anisotropic spectrum of random fluid displacements produces and a positive TE cross-correlation. Furthermore, we find that large EE/BB and positive TE are due primarily to longitudinal, rather than transverse, modes in the random-displacement field, providing, perhaps, some clue to the mechanism that stirs the ISM. Future investigations involving the spatial dependence of the EE/BB ratio, TE correlation, and local departures from statistical isotropy in dust-polarization maps, as well as further tests of some of the assumptions in this analysis, are outlined. This work may also aid in the improvement of foreground-separation techniques for studies of cosmic microwave background polarization.

Evolutionary Description of Giant Molecular Cloud Mass Functions on Galactic Disks

Abstract Recent radio observations show that giant molecular cloud (GMC) mass functions noticeably vary across galactic disks. High-resolution magnetohydrodynamics simulations show that multiple episodes of compression are required for creating a molecular cloud in the magnetized interstellar medium. In this article, we formulate the evolution equation for the GMC mass function to reproduce the observed profiles, for which multiple compressions are driven by a network of expanding shells due to H ii regions and supernova remnants. We introduce the cloud–cloud collision (CCC) terms in the evolution equation in contrast to previous work (Inutsuka et al.). The computed time evolution suggests that the GMC mass function slope is governed by the ratio of GMC formation timescale to its dispersal timescale, and that the CCC effect is limited only in the massive end of the mass function. In addition, we identify a gas resurrection channel that allows the gas dispersed by massive stars to regenerate GMC populations or to accrete onto pre-existing GMCs. Our results show that almost all of the dispersed gas contributes to the mass growth of pre-existing GMCs in arm regions whereas less than 60% contributes in inter-arm regions. Our results also predict that GMC mass functions have a single power-law exponent in the mass range <105.5 (where represents the solar mass), which is well characterized by GMC self-growth and dispersal timescales. Measurement of the GMC mass function slope provides a powerful method to constrain those GMC timescales and the gas resurrecting factor in various environments across galactic disks.

Planckintermediate results

We study the statistical properties of interstellar dust polarization at high Galactic latitude, using the Stokes parameter Planck maps at 353 GHz. Our aim is to advance the understanding of the magnetized interstellar medium (ISM), and to provide a model of the polarized dust foreground for cosmic microwave background component-separation procedures. Focusing on the southern Galactic cap, we examine the statistical distributions of the polarization fraction ($p$) and angle ($\psi$) to characterize the ordered and turbulent components of the Galactic magnetic field (GMF) in the solar neighbourhood. We relate patterns at large angular scales in polarization to the orientation of the mean (ordered) GMF towards Galactic coordinates $(l_0,b_0)=(70^\circ \pm 5^\circ,24^\circ \pm 5^\circ)$. The histogram of $p$ shows a wide dispersion up to 25 %. The histogram of $\psi$ has a standard deviation of $12^\circ$ about the regular pattern expected from the ordered GMF. We use these histograms to build a phenomenological model of the turbulent component of the GMF, assuming a uniform effective polarization fraction ($p_0$) of dust emission. To model the Stokes parameters, we approximate the integration along the line of sight (LOS) as a sum over a set of $N$ independent polarization layers, in each of which the turbulent component of the GMF is obtained from Gaussian realizations of a power-law power spectrum. We are able to reproduce the observed $p$ and $\psi$ distributions using: a $p_0$ value of (26 $\pm$ 3)%; a ratio of 0.9 $\pm$ 0.1 between the strengths of the turbulent and mean components of the GMF; and a small value of $N$. We relate the polarization layers to the density structure and to the correlation length of the GMF along the LOS.

Planckintermediate results

The role of the magnetic field in the formation of the filamentary structures observed in the interstellar medium (ISM) is a debated topic. The Planck all-sky maps of linearly polarized emission from dust at 353GHz provide the required combination of imaging and statistics to study the correlation between the structures of the Galactic magnetic field and of interstellar matter, both in the diffuse ISM and in molecular clouds. The data reveal structures, or ridges, in the intensity map with counterparts in the Stokes Q and/or U maps. We focus on structures at intermediate and high Galactic latitudes with column density from $10^{20}$ to $10^{22}$ cm$^{-2}$. We measure the magnetic field orientation on the plane of the sky from the polarization data, and present an algorithm to estimate the orientation of the ridges from the dust intensity map. We use analytical models to account for projection effects. Comparing polarization angles on and off the structures, we estimate the mean ratio between the strengths of the turbulent and mean components of the magnetic field to be between 0.6 and 1.0, with a preferred value of 0.8. We find that the ridges are preferentially aligned with the magnetic field measured on the structures. This trend becomes more striking for increasing polarization fraction and decreasing column density. We interpret the increase of alignment with polarization fraction as a consequence of projections effects. The decrease of alignment for high column density is not due to a loss of correlation between the structures and the geometry of the magnetic field. In molecular complexes, we observe structures perpendicular to the magnetic field, which cannot be accounted for by projection effects. We discuss our results in the context of models and MHD simulations, which describe the formation of structures in the magnetized ISM.

COSMIC RAYS CAN DRIVE STRONG OUTFLOWS FROM GAS-RICH HIGH-REDSHIFT DISK GALAXIES

We present simulations of the magnetized interstellar medium (ISM) in models of massive star forming (40 Msun / yr) disk galaxies with high gas surface densities (~100 Msun / pc^2) similar to observed star forming high-redshift disks. We assume that type II supernovae deposit 10 per cent of their energy into the ISM as cosmic rays and neglect the additional deposition of thermal energy or momentum. With a typical Galactic diffusion coefficient for CRs (3e28 cm^2 / s) we demonstrate that this process alone can trigger the local formation of a strong low density galactic wind maintaining vertically open field lines. Driven by the additional pressure gradient of the relativistic fluid the wind speed can exceed 1000 km/s, much higher than the escape velocity of the galaxy. The global mass loading, i.e. the ratio of the gas mass leaving the galactic disk in a wind to the star formation rate becomes of order unity once the system has settled into an equilibrium. We conclude that relativistic particles accelerated in supernova remnants alone provide a natural and efficient mechanism to trigger winds similar to observed mass-loaded galactic winds in high-redshift galaxies. These winds also help explaining the low efficiencies for the conversion of gas into stars in galaxies as well as the early enrichment of the intergalactic medium with metals. This mechanism can be at least of similar importance than the traditionally considered momentum feedback from massive stars and thermal and kinetic feedback from supernova explosions.

THE MAXIMUM ENERGY OF ACCELERATED PARTICLES IN RELATIVISTIC COLLISIONLESS SHOCKS

The afterglow emission from gamma-ray bursts (GRBs) is usually interpreted as synchrotron radiation from electrons accelerated at the GRB external shock, that propagates with relativistic velocities into the magnetized interstellar medium. By means of multi-dimensional particle-in-cell simulations, we investigate the acceleration performance of weakly magnetized relativistic shocks, in the magnetization range 0<sigma<1e-1. The pre-shock magnetic field is orthogonal to the flow, as generically expected for relativistic shocks. We find that relativistic perpendicular shocks propagating in electron-positron plasmas are efficient particle accelerators if the magnetization is sigma<1e-3. For electron-ion plasmas, the transition to efficient acceleration occurs for sigma<3e-5. Here, the acceleration process proceeds similarly for the two species, since the electrons enter the shock nearly in equipartition with the ions, as a result of strong pre-heating in the self-generated upstream turbulence. In both electron-positron and electron-ion shocks, we find that the maximum energy of the accelerated particles scales in time as t^(1/2). This scaling is shallower than the so-called (and commonly assumed) Bohm limit, and it naturally results from the small-scale nature of the Weibel turbulence generated in the shock layer. In magnetized plasmas, the energy of the accelerated particles increases until it reaches a saturation value that scales with the magnetization as sigma^(-1/4). Further energization is prevented by the fact that the self-generated turbulence is confined within a finite region of thickness proportional to sigma^(-1/2) around the shock. Our results can provide physically-grounded inputs for models of non-thermal emission from a variety of astrophysical sources, with particular relevance to GRB afterglows.

Clumping of ejecta and accelerated cosmic rays in the evolution of type Ia SNRs

AbstractWe investigate the role played by initial clumping of ejecta and by efficient acceleration of cosmic rays (CRs) in determining the density structure of the post-shock region of a Type Ia supernova remnant (SNR) through detailed 3D MHD modeling. Our model describes the expansion of a SNR through a magnetized interstellar medium (ISM), including the initial clumping of ejecta and the effects on shock dynamics due to back-reaction of accelerated CRs. The model predictions are compared to the observations of SN 1006. We found that the back-reaction of accelerated CRs alone cannot reproduce the observed separation between the forward shock (FS) and the contact discontinuity (CD) unless the energy losses through CR acceleration and escape are very large and independent of the obliquity angle. On the contrary, the clumping of ejecta can naturally reproduce the observed small separation and the occurrence of protrusions observed in SN 1006, even without the need of accelerated CRs. We conclude that FS-CD separation is a probe of the ejecta structure at the time of explosion rather than a probe of the efficiency of CR acceleration in young SNRs.

THE SINO-GERMAN λ6 CM POLARIZATION SURVEY OF THE GALACTIC PLANE: A SUMMARY

We have finished the λ6 cm polarization survey of the Galactic plane using the Urumqi 25 m radio telescope. It covers 10° ≤ l ≤ 230° in Galactic longitude and |b| ≤ 5° in Galactic latitude. The new polarization maps not only reveal new properties of the diffuse magnetized interstellar medium, but also are very useful for studying individual objects such as HII regions, which may act as Faraday screens with strong regular magnetic fields inside, and supernova remnants for their polarization properties and spectra. The high sensitivity of the survey enables us to discover two new SNRs G178.2–4.2 and G25.3–2.1 and a number of HII regions.

ROLE OF EJECTA CLUMPING AND BACK-REACTION OF ACCELERATED COSMIC RAYS IN THE EVOLUTION OF TYPE Ia SUPERNOVA REMNANTS

We investigate the role played by initial clumping of ejecta and by efficient acceleration of cosmic rays (CRs) in determining the density structure of the post-shock region of a Type Ia supernova remnant (SNR) through detailed 3D MHD modeling. Our model describes the expansion of a SNR through a magnetized interstellar medium (ISM), including the initial clumping of ejecta and the effects on shock dynamics due to back-reaction of accelerated CRs. The model predictions are compared to the observations of SN 1006. We found that the back-reaction of accelerated CRs alone cannot reproduce the observed separation between the forward shock (FS) and the contact discontinuity (CD) unless the energy losses through CR acceleration and escape are very large and independent of the obliquity angle. On the contrary, the clumping of ejecta can naturally reproduce the observed small separation and the occurrence of protrusions observed in SN 1006, even without the need of accelerated CRs. We conclude that FS-CD separation is a probe of the ejecta structure at the time of explosion rather than a probe of the efficiency of CR acceleration in young SNRs.

THE EFFECT OF ANISOTROPIC CONDUCTION ON THE THERMAL INSTABILITY IN THE INTERSTELLAR MEDIUM

Thermal instability (TI) can strongly affect the structure and dynamics of the interstellar medium (ISM) in the Milky Way and other disk galaxies. Thermal conduction plays an important role in the TI by stabilizing small scales and limiting the size of the smallest condensates. In the magnetized ISM, however, heat is conducted anisotropically (primarily along magnetic field lines). We investigate the effects of anisotropic thermal conduction on the nonlinear regime of the TI by performing two-dimensional magnetohydrodynamic simulations. We present models with magnetic fields of different initial geometries and strengths, and compare them to hydrodynamic models with isotropic conduction. We find anisotropic conduction does not significantly alter the overall density and temperature statistics in the saturated state of the TI. However, it can strongly affect the shapes and sizes of cold clouds formed by the TI. For example, for uniform initial fields long filaments of cold gas are produced that are reminiscent of some observed HI clouds. For initially tangled fields, such filaments are not produced. We also show that anisotropic conduction suppresses turbulence generated by evaporative flows from the surfaces of cold blobs, which may have implications for mechanisms for driving turbulence in the ISM.

A Sino-Germanλ6 cm polarization survey of the Galactic plane

Linear polarization of diffuse Galactic emission is a signature of magnetic fields in the interstellar medium of our Galaxy. Observations at high frequencies are less affected by Faraday depolarization than those at lower frequencies and are able to detect polarized emission from more distant Galactic regions. We attempt to perform a sensitive survey of the polarized emission from the Galactic disk at 6 cm wavelength. We made polarization observations of the Galactic plane using the Urumqi 25-m telescope at 6 cm covering the area of 60deg< l <129deg and |b|<5deg. Missing large-scale structures in polarization were restored by extrapolation of the WMAP polarization data. We present the 6cm total intensity and linear polarization maps of the surveyed region. We identify two new extended HII regions G98.3-1.6 and G119.6+0.4 in this region. Numerous polarized patches and depolarization structures are visible in the polarization maps. Depolarization along the periphery of a few HII complexes was detected and can be explained by a Faraday screen model. We discuss some prominent depolarization HII regions, which have regular magnetic fields of several $\mu$G. Structure functions of U, Q and PI images of the entire 6cm survey region of 10deg< l < 230deg exhibit much larger fluctuation power towards the inner Galaxy, which suggests a higher turbulence in the arm regions of the inner Galaxy. The Sino-German 6cm survey reveals new properties of the diffuse magnetized interstellar medium. The survey is also very useful for studying individual objects such as HII regions, which act as Faraday screens and have high rotation measures and therefore strong regular magnetic fields inside the regions.

Effects of non-uniform interstellar magnetic field on synchrotron X-ray and inverse-Comptonγ-ray morphology of supernova remnants

Observations of SNRs in X-ray and gamma-ray bands promise to contribute with important information in our understanding on the nature of galactic cosmic rays. The analysis of SNRs images collected in different energy bands requires the support of theoretical modeling of synchrotron and inverse Compton (IC) emission. We develop a numerical code (REMLIGHT) to synthesize, from MHD simulations, the synchrotron radio, X-ray and IC gamma-ray emission from SNRs expanding in non-uniform interstellar medium (ISM) and/or non-uniform interstellar magnetic field (ISMF). As a first application, the code is used to investigate the effects of non-uniform ISMF on the SNR morphology in the non-thermal X-ray and gamma-ray bands. We perform 3D MHD simulations of a spherical SNR shock expanding through a magnetized ISM with a gradient of ambient magnetic field strength. The model includes an approximate treatment of upstream magnetic field amplification and the effect of shock modification due to back reaction of accelerated cosmic rays. From the simulations, we synthesize the synchrotron radio, X-ray and IC gamma-ray emission with REMLIGHT, making different assumptions about the details of acceleration and injection of relativistic electrons. A gradient of the ambient magnetic field strength induces asymmetric morphologies in radio, X-ray and gamma-ray bands independently from the model of electron injection if the gradient has a component perpendicular to the line-of-sight. The degree of asymmetry of the remnant morphology depends on the details of the electron injection and acceleration and is different in the radio, X-ray, and gamma-ray bands. The non-thermal X-ray morphology is the most sensitive to the gradient, showing the highest degree of asymmetry. The IC gamma-ray emission is weakly sensitive to the non-uniform ISMF, the degree of asymmetry of the SNR morphology being the lowest in this band.

The haloes of planetary nebulae in the mid-infrared: evidence for interaction with the interstellar medium

The motion of planetary nebulae through the interstellar medium (ISM) is thought to lead to a variety of observational consequences, including the formation of bright rims, deformation and fragmentation of the shells, and a shift of the central stars away from the geometric centres of the envelopes. These and other characteristics have been noted through imaging in the visual wavelength regime. We report further observations of such shells taken in the mid-infrared (MIR), acquired through programmes of Infrared Array Camera imaging undertaken using the SpitzerSpace Telescope. NGC 2440 and NGC 6629 are shown to possess likely interacting haloes, together with ram-pressure-stripped material to one side of their shells. Similarly, the outer haloes of NGC 3242 and NGC 6772 appear to have been fragmented through Rayleigh–Taylor (RT) instabilities, leading to a possible flow of ISM material towards the inner portions of their envelopes. If this interpretation is correct, then it would suggest that NGC 3242 is moving towards the NE, a suggestion which is also supported through the presence of a 60 μm tail extending in the opposite direction, and curved bands of Hα emission in the direction of motion – components which may arise through RT instabilities in the magnetized ISM. NGC 2438 possesses strong scalloping at the outer limits of its asymptotic giant branch (AGB) halo, probably reflecting RT instabilities at the nebular/ISM interface We also note that the interior structure of the source has been interpreted in terms of a recombining shell, a hypothesis which may not be consistent with the central star luminosities. Finally, we point out that two of the rims (and likely shock interfaces) appear to have a distinct signature in the MIR, whereby relative levels of 8.0 μm emission are reduced. This may imply that the grain emission agents are depleted in the post-shock AGB regimes.

GLOBAL GALACTIC DYNAMO DRIVEN BY COSMIC RAYS AND EXPLODING MAGNETIZED STARS

We report first results of first global galactic-scale CR-MHD simulations of cosmic-ray-driven dynamo. We investigate the dynamics of magnetized interstellar medium (ISM), which is dynamically coupled with the cosmic-ray (CR) gas. We assume that exploding stars deposit small-scale, randomly oriented, dipolar magnetic fields into the differentially rotating ISM, together with a portion of cosmic rays, accelerated in supernova shocks. We conduct numerical simulations with the aid of a new parallel MHD code PIERNIK. We find that the initial magnetization of galactic disks by exploding magnetized stars forms a favourable conditions for the cosmic-ray-driven dynamo. We demonstrate that dipolar magnetic fields supplied on small SN-remnant scales, can be amplified exponentially, by the CR-driven dynamo, to the present equipartition values, and transformed simultaneously to large galactic-scales. The resulting magnetic field structure in an evolved galaxy appears spiral in the face-on view and reveals the so called X-shaped structure in the edge-on view.

MAGNETOHYDRODYNAMIC SIMULATIONS OF DISK GALAXY FORMATION: THE MAGNETIZATION OF THE COLD AND WARM MEDIUM

We modeled the formation and early evolution of disk galaxies with a magnetized interstellar medium using magnetohydrodynamic (MHD) adaptive mesh refinement simulations. For a 1010 M☉ halo with initial Navarro–Frenk–White dark matter and gas profiles, we impose a uniform 10−9 G magnetic field and follow its collapse, disk formation and evolution up to 1 Gyr. We find that the initial magnetic fields are quickly amplified by the differentially rotating turbulent disk with the amplification rate roughly one e-folding per orbit. After the initial rapid amplification lasting ∼500 Myr, subsequent field amplification appears self-regulated. The field strengths in the self-regulated regime have similar strength as the observed fields in the Milky Way galaxy both in the warm and the cold H i phases. Since supernova explosions, which we neglected in our current model, are likely to further amplify the magnetic field, our calculation suggests that Milky Way strength magnetic field might be common in high redshift disk galaxies. After saturation, highly magnetized material also begins to form above and below the disk, which may affect subsequent galaxy evolution, especially mergers, significantly. The global azimuthal magnetic fields reverse at different radii and the amplitude declines as a function of radius of the disk. We also find that magnetic force can provide further support in the cold gas and lead to a decline of the amount of cold gas at high density which may lead to a decline in the star formation rate.

Simulating polarized Galactic synchrotron emission at all frequencies

<i>Context. <i/>Galactic synchrotron emission, rotation measure (RM), and the deflection of ultra-high-energy-cosmic-rays (UHECR) permit detailed studies of the Galactic magnetized inter-stellar-medium (ISM). The synchrotron emission has also to be characterized properly to enable its distinction from other astrophysically interesting signals such as the CMB.<i>Aims. <i/>We present a publicly available code called hammurabi for generating mock, polarized observations of Galactic synchrotron emission with telescopes such as LOFAR, SKA, Planck, and WMAP, based on model inputs for the Galactic magnetic field (GMF), the cosmic-ray density distribution, and the thermal electron density. We also present mock UHECR deflection measure (UDM) maps based on model inputs for the GMF. In future, when UHECR sources are identified, this will allow us to use UDM as a probe of GMF in a way similar to how polarized radio sources enable us to define rotation measures.<i>Methods. <i/>To demonstrate the code's abilities, mock observations are compared with real data. This allows us to constrain the input parameters of our simulations with a focus on large-scale magnetic field properties.<i>Results. <i/>Magnetized ISM models in the literature appear to be unable to reproduce any additional observational data not included in their design.<i>Conclusions. <i/>As expected, attempts to model the synchrotron, UHECR-deflection, and RM input parameters indicate that any additional observational data set significantly increases the constraints on the models. The hammurabi code addresses this by allowing one to perform simulations of several different data sets simultaneously, providing a more reliable constraint of the magnetized ISM.