-
On the effect of galactic outflows in cosmological simulations of disc galaxies
Authors:
Milena Valentini,
Giuseppe Murante,
Stefano Borgani,
Pierluigi Monaco,
Alessandro Bressan,
Alexander M. Beck
Abstract:
We investigate the impact of galactic outflow modelling on the formation and evolution of a disc galaxy, by performing a suite of cosmological simulations with zoomed-in initial conditions of a Milky Way-sized halo. We verify how sensitive the general properties of the simulated galaxy are to the way in which stellar feedback triggered outflows are implemented, keeping initial conditions, simulati…
▽ More
We investigate the impact of galactic outflow modelling on the formation and evolution of a disc galaxy, by performing a suite of cosmological simulations with zoomed-in initial conditions of a Milky Way-sized halo. We verify how sensitive the general properties of the simulated galaxy are to the way in which stellar feedback triggered outflows are implemented, keeping initial conditions, simulation code and star formation (SF) model all fixed. We present simulations that are based on a version of the GADGET3 code where our sub-resolution model is coupled with an advanced implementation of Smoothed Particle Hydrodynamics that ensures a more accurate fluid sampling and an improved description of gas mixing and hydrodynamical instabilities. We quantify the strong interplay between the adopted hydrodynamic scheme and the sub-resolution model describing SF and feedback. We consider four different galactic outflow models, including the one introduced by Dalla Vecchia and Schaye (2012) and a scheme that is inspired by the Springel and Hernquist (2003) model. We find that the sub-resolution prescriptions adopted to generate galactic outflows are the main shaping factor of the stellar disc component at low redshift. The key requirement that a feedback model must have to be successful in producing a disc-dominated galaxy is the ability to regulate the high-redshift SF (responsible for the formation of the bulge component), the cosmological infall of gas from the large-scale environment, and gas fall-back within the galactic radius at low redshift, in order to avoid a too high SF rate at $z=0$.
△ Less
Submitted 31 December, 2017; v1 submitted 29 May, 2017;
originally announced May 2017.
-
Simulations of the Galaxy Cluster CIZA J2242.8+5301 I: Thermal Model and Shock Properties
Authors:
J. M. F. Donnert,
A. M. Beck,
K. Dolag,
H. J. A. Röttgering
Abstract:
The giant radio relic in CIZA J2242.8+5301 is likely evidence of a Mpc sized shock in a massive merging galaxy cluster. However, the exact shock properties are still not clearly determined. In particular, the Mach number derived from the integrated radio spectrum exceeds the Mach number derived from the X-ray temperature jump by a factor of two. We present here a numerical study, aiming for a mode…
▽ More
The giant radio relic in CIZA J2242.8+5301 is likely evidence of a Mpc sized shock in a massive merging galaxy cluster. However, the exact shock properties are still not clearly determined. In particular, the Mach number derived from the integrated radio spectrum exceeds the Mach number derived from the X-ray temperature jump by a factor of two. We present here a numerical study, aiming for a model that is consistent with the majority of observations of this galaxy cluster. We first show that in the northern shock upstream X-ray temperature and radio data are consistent with each other. We then derive progenitor masses for the system using standard density profiles, X-ray properties and the assumption of hydrostatic equilibrium. We find a class of models that is roughly consistent with weak lensing data, radio data and some of the X-ray data. Assuming a cool-core versus non-cool-core merger, we find a fiducial model with a total mass of $1.6 \times 10^{15}\,M_\odot$, a mass ratio of 1.76 and a Mach number that is consistent with estimates from the radio spectrum. We are not able to match X-ray derived Mach numbers, because even low mass models over-predict the X-ray derived shock speeds. We argue that deep X-ray observations of CIZA J2242.8+5301 will be able to test our model and potentially reconcile X-ray and radio derived Mach numbers in relics.
△ Less
Submitted 8 April, 2017; v1 submitted 16 March, 2017;
originally announced March 2017.
-
The history of chemical enrichment in the intracluster medium from cosmological simulations
Authors:
V. Biffi,
S. Planelles,
S. Borgani,
D. Fabjan,
E. Rasia,
G. Murante,
L. Tornatore,
K. Dolag,
G. L. Granato,
M. Gaspari,
A. M. Beck
Abstract:
The distribution of metals in the intracluster medium (ICM) of galaxy clusters provides valuable information on their formation and evolution, on the connection with the cosmic star formation and on the effects of different gas processes. By analyzing a sample of simulated galaxy clusters, we study the chemical enrichment of the ICM, its evolution, and its relation with the physical processes incl…
▽ More
The distribution of metals in the intracluster medium (ICM) of galaxy clusters provides valuable information on their formation and evolution, on the connection with the cosmic star formation and on the effects of different gas processes. By analyzing a sample of simulated galaxy clusters, we study the chemical enrichment of the ICM, its evolution, and its relation with the physical processes included in the simulation and with the thermal properties of the core. These simulations, consisting of re-simulations of 29 Lagrangian regions performed with an upgraded version of the SPH GADGET-3 code, have been run including two different sets of baryonic physics: one accounts for radiative cooling, star formation, metal enrichment and supernova (SN) feedback, and the other one further includes the effects of feedback from active galactic nuclei (AGN). In agreement with observations, we find an anti-correlation between entropy and metallicity in cluster cores, and similar radial distributions of heavy-element abundances and abundance ratios out to large cluster-centric distances (~R180). In the outskirts, namely outside of ~0.2R180, we find a remarkably homogeneous metallicity distribution, with almost flat profiles of the elements produced by either SNIa or SNII. We investigated the origin of this phenomenon and discovered that it is due to the widespread displacement of metal-rich gas by early (z>2-3) AGN powerful bursts, acting on small high-redshift haloes. Our results also indicate that the intrinsic metallicity of the hot gas for this sample is on average consistent with no evolution between z=2 and z=0, across the entire radial range.
△ Less
Submitted 27 January, 2017;
originally announced January 2017.
-
Pressure of the hot gas in simulations of galaxy clusters
Authors:
S. Planelles,
D. Fabjan,
S. Borgani,
G. Murante,
E. Rasia,
V. Biffi,
N. Truong,
C. Ragone-Figueroa,
G. L. Granato,
K. Dolag,
E. Pierpaoli,
A. M. Beck,
Lisa K. Steinborn,
M. Gaspari
Abstract:
We analyze the radial pressure profiles, the ICM clumping factor and the Sunyaev-Zel'dovich (SZ) scaling relations of a sample of simulated galaxy clusters and groups identified in a set of hydrodynamical simulations based on an updated version of the TreePM-SPH GADGET-3 code. Three different sets of simulations are performed: the first assumes non-radiative physics, the others include, among othe…
▽ More
We analyze the radial pressure profiles, the ICM clumping factor and the Sunyaev-Zel'dovich (SZ) scaling relations of a sample of simulated galaxy clusters and groups identified in a set of hydrodynamical simulations based on an updated version of the TreePM-SPH GADGET-3 code. Three different sets of simulations are performed: the first assumes non-radiative physics, the others include, among other processes, AGN and/or stellar feedback. Our results are analyzed as a function of redshift, ICM physics, cluster mass and cluster cool-coreness or dynamical state. In general, the mean pressure profiles obtained for our sample of groups and clusters show a good agreement with X-ray and SZ observations. Simulated cool-core (CC) and non-cool-core (NCC) clusters also show a good match with real data. We obtain in all cases a small (if any) redshift evolution of the pressure profiles of massive clusters, at least back to z=1. We find that the clumpiness of gas density and pressure increases with the distance from the cluster center and with the dynamical activity. The inclusion of AGN feedback in our simulations generates values for the gas clumping ($\sqrt C_ρ\sim 1.2$ at $R_{200}$) in good agreement with recent observational estimates. The simulated $Y_{SZ}-M$ scaling relations are in good accordance with several observed samples, especially for massive clusters. As for the scatter of these relations, we obtain a clear dependence on the cluster dynamical state, whereas this distinction is not so evident when looking at the subsamples of CC and NCC clusters.
△ Less
Submitted 21 December, 2016;
originally announced December 2016.
-
nIFTy galaxy cluster simulations V: Investigation of the Cluster Infall Region
Authors:
Jake Arthur,
Frazer R. Pearce,
Meghan E. Gray,
Pascal J. Elahi,
Alexander Knebe,
Alexander M. Beck,
Weiguang Cui,
Daniel Cunnama,
Romeel Davé,
Sean February,
Shuiyao Huang,
Neal Katz,
Scott T. Kay,
Ian G. McCarthy,
Giuseppe Murante,
Valentin Perret,
Chris Power,
Ewald Puchwein,
Alexandro Saro,
Federico Sembolini,
Romain Teyssier,
Gustavo Yepes
Abstract:
We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $Λ$CDM galaxy cluster studied by Elahi et al. (2016) at z=0. The $1.1\times10^{15}h^{-1}\text{M}_{\odot}$ galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and subgrid schemes. All models completed a…
▽ More
We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $Λ$CDM galaxy cluster studied by Elahi et al. (2016) at z=0. The $1.1\times10^{15}h^{-1}\text{M}_{\odot}$ galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and subgrid schemes. All models completed a dark-matter only, non-radiative and full-physics run from the same initial conditions. The simulations contain dark matter and gas with mass resolution $m_{\text{DM}}=9.01\times 10^8h^{-1}\text{M}_{\odot}$ and $m_{\text{gas}}=1.9\times 10^8h^{-1}\text{M}_{\odot}$ respectively. We find that the synthetic cluster is surrounded by clear filamentary structures that contain ~60% of haloes in the infall region with mass ~$10^{12.5} - 10^{14} h^{-1}\text{M}_{\odot}$, including 2-3 group-sized haloes ($> 10^{13}h^{-1}\text{M}_{\odot}$). However, we find that only ~10% of objects in the infall region are subhaloes residing in haloes, which may suggest that there is not much ongoing preprocessing occurring in the infall region at z=0. By examining the baryonic content contained within the haloes, we also show that the code-to-code scatter in stellar fraction across all halo masses is typically ~2 orders of magnitude between the two most extreme cases, and this is predominantly due to the differences in subgrid schemes and calibration procedures that each model uses. Models that do not include AGN feedback typically produce too high stellar fractions compared to observations by at least ~1 order of magnitude.
△ Less
Submitted 23 September, 2016;
originally announced September 2016.
-
Cosmological hydrodynamical simulations of galaxy clusters: X-ray scaling relations and their evolution
Authors:
N. Truong,
E. Rasia,
P. Mazzotta,
S. Planelles,
V. Biffi,
D. Fabjan,
A. M. Beck,
S. Borgani,
K. Dolag,
M. Gaspari,
G. L. Granato,
G. Murante,
C. Ragone-Figueroa,
L. K. Steinborn
Abstract:
We analyse cosmological hydrodynamical simulations of galaxy clusters to study the X-ray scaling relations between total masses and observable quantities such as X-ray luminosity, gas mass, X-ray temperature, and $Y_{X}$. Three sets of simulations are performed with an improved version of the smoothed particle hydrodynamics GADGET-3 code. These consider the following: non-radiative gas, star forma…
▽ More
We analyse cosmological hydrodynamical simulations of galaxy clusters to study the X-ray scaling relations between total masses and observable quantities such as X-ray luminosity, gas mass, X-ray temperature, and $Y_{X}$. Three sets of simulations are performed with an improved version of the smoothed particle hydrodynamics GADGET-3 code. These consider the following: non-radiative gas, star formation and stellar feedback, and the addition of feedback by active galactic nuclei (AGN). We select clusters with $M_{500} > 10^{14} M_{\odot} E(z)^{-1}$, mimicking the typical selection of Sunyaev-Zeldovich samples. This permits to have a mass range large enough to enable robust fitting of the relations even at $z \sim 2$. The results of the analysis show a general agreement with observations. The values of the slope of the mass-gas mass and mass-temperature relations at $z=2$ are 10 per cent lower with respect to $z=0$ due to the applied mass selection, in the former case, and to the effect of early merger in the latter. We investigate the impact of the slope variation on the study of the evolution of the normalization. We conclude that cosmological studies through scaling relations should be limited to the redshift range $z=0-1$, where we find that the slope, the scatter, and the covariance matrix of the relations are stable. The scaling between mass and $Y_X$ is confirmed to be the most robust relation, being almost independent of the gas physics. At higher redshifts, the scaling relations are sensitive to the inclusion of AGNs which influences low-mass systems. The detailed study of these objects will be crucial to evaluate the AGN effect on the ICM.
△ Less
Submitted 25 January, 2018; v1 submitted 30 June, 2016;
originally announced July 2016.
-
On the nature of hydrostatic equilibrium in galaxy clusters
Authors:
V. Biffi,
S. Borgani,
G. Murante,
E. Rasia,
S. Planelles,
G. L. Granato,
C. Ragone-Figueroa,
A. M. Beck,
M. Gaspari,
K. Dolag
Abstract:
In this paper we investigate the level of hydrostatic equilibrium (HE) in the intra-cluster medium of simulated galaxy clusters, extracted from state-of-the-art cosmological hydrodynamical simulations performed with the Smoothed-Particle-Hydrodynamic code GADGET-3. These simulations include several physical processes, among which stellar and AGN feedback, and have been performed with an improved v…
▽ More
In this paper we investigate the level of hydrostatic equilibrium (HE) in the intra-cluster medium of simulated galaxy clusters, extracted from state-of-the-art cosmological hydrodynamical simulations performed with the Smoothed-Particle-Hydrodynamic code GADGET-3. These simulations include several physical processes, among which stellar and AGN feedback, and have been performed with an improved version of the code that allows for a better description of hydrodynamical instabilities and gas mixing processes. Evaluating the radial balance between the gravitational and hydrodynamical forces, via the gas accelerations generated, we effectively examine the level of HE in every object of the sample, its dependence on the radial distance from the center and on the classification of the cluster in terms of either cool-coreness or dynamical state. We find an average deviation of 10-20% out to the virial radius, with no evident distinction between cool-core and non-cool-core clusters. Instead, we observe a clear separation between regular and disturbed systems, with a more significant deviation from HE for the disturbed objects. The investigation of the bias between the hydrostatic estimate and the total gravitating mass indicates that, on average, this traces very well the deviation from HE, even though individual cases show a more complex picture. Typically, in the radial ranges where mass bias and deviation from HE are substantially different, the gas is characterized by a significant amount of random motions (>~30 per cent), relative to thermal ones. As a general result, the HE-deviation and mass bias, at given interesting distance from the cluster center, are not very sensitive to the temperature inhomogeneities in the gas.
△ Less
Submitted 7 June, 2016;
originally announced June 2016.
-
nIFTy Galaxy Cluster simulations IV: Quantifying the Influence of Baryons on Halo Properties
Authors:
Weiguang Cui,
Chris Power,
Alexander Knebe,
Scott T. Kay,
Federico Sembolini,
Pascal J. Elahi,
Gustavo Yepes,
Frazer Pearce,
Daniel Cunnama,
Alexander M. Beck,
Claudio Dalla Vecchia,
Romeel Davé,
Sean February,
Shuiyao Huang,
Alex Hobbs,
Neal Katz,
Ian G. McCarthy,
Giuseppe Murante,
Valentin Perret,
Ewald Puchwein,
Justin I. Read,
Alexandro Saro,
Romain Teyssier,
Robert J. Thacker
Abstract:
Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a…
▽ More
Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency (<= 20%) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R200. However, we see larger differences for quantities within R2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include AGN feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers.
△ Less
Submitted 22 February, 2016;
originally announced February 2016.
-
A supernova scenario for magnetic fields and rotation measures in galaxies
Authors:
Alexander M. Beck,
Klaus Dolag
Abstract:
We present a model for the seeding and evolution of magnetic fields in galaxies by supernovae (SN). SN explosions during galaxy assembly provide seed fields, which are subsequently amplified by compression, shear flows and random motions. Our model explains the origin of microG magnetic fields within galactic structures. We implement our model in the MHD version of the cosmological simulation code…
▽ More
We present a model for the seeding and evolution of magnetic fields in galaxies by supernovae (SN). SN explosions during galaxy assembly provide seed fields, which are subsequently amplified by compression, shear flows and random motions. Our model explains the origin of microG magnetic fields within galactic structures. We implement our model in the MHD version of the cosmological simulation code Gadget-3 and couple it with a multi-phase description of the interstellar medium. We perform simulations of Milky Way-like galactic halo formation and analyze the distribution and strength of the magnetic field. We investigate the intrinsic rotation measure (RM) evolution and find RM values exceeding 1000 rad/m*m at high redshifts and RM values around 10 rad/m*m at present-day. We compare our simulations to a limited set of observational data points and find encouraging similarities. In our model, galactic magnetic fields are a natural consequence of the very basic processes of star formation and galaxy assembly.
△ Less
Submitted 27 November, 2015;
originally announced December 2015.
-
nIFTY galaxy cluster simulations III: The Similarity & Diversity of Galaxies & Subhaloes
Authors:
Pascal J. Elahi,
Alexander Knebe,
Frazer R. Pearce,
Chris Power,
Gustavo Yepes,
Weiguang Cui,
Daniel Cunnama,
Scott T. Kay,
Federico Sembolini,
Alexander M. Beck,
Romeel Davé,
Sean February,
Shuiyao Huang,
Neal Katz,
Ian G. McCarthy,
Giuseppe Murante,
Valentin Perret,
Ewald Puchwein,
Alexandro Saro,
Romain Teyssier
Abstract:
We examine subhaloes and galaxies residing in a simulated LCDM galaxy cluster ($M^{\rm crit}_{200}=1.1\times10^{15}M_\odot/h$) produced by hydrodynamical codes ranging from classic Smooth Particle Hydrodynamics (SPH), newer SPH codes, adaptive and moving mesh codes. These codes use subgrid models to capture galaxy formation physics. We compare how well these codes reproduce the same subhaloes/gala…
▽ More
We examine subhaloes and galaxies residing in a simulated LCDM galaxy cluster ($M^{\rm crit}_{200}=1.1\times10^{15}M_\odot/h$) produced by hydrodynamical codes ranging from classic Smooth Particle Hydrodynamics (SPH), newer SPH codes, adaptive and moving mesh codes. These codes use subgrid models to capture galaxy formation physics. We compare how well these codes reproduce the same subhaloes/galaxies in gravity only, non-radiative hydrodynamics and full feedback physics runs by looking at the overall subhalo/galaxy distribution and on an individual objects basis. We find the subhalo population is reproduced to within $\lesssim10\%$ for both dark matter only and non-radiative runs, with individual objects showing code-to-code scatter of $\lesssim0.1$ dex, although the gas in non-radiative simulations shows significant scatter. Including feedback physics significantly increases the diversity. Subhalo mass and $V_{max}$ distributions vary by $\approx20\%$. The galaxy populations also show striking code-to-code variations. Although the Tully-Fisher relation is similar in almost all codes, the number of galaxies with $10^{9}M_\odot/h\lesssim M_*\lesssim 10^{12}M_\odot/h$ can differ by a factor of 4. Individual galaxies show code-to-code scatter of $\sim0.5$ dex in stellar mass. Moreover, strong systematic differences exist, with some codes producing galaxies $70\%$ smaller than others. The diversity partially arises from the inclusion/absence of AGN feedback. Our results combined with our companion papers demonstrate that subgrid physics is not just subject to fine-tuning, but the complexity of building galaxies in all environments remains a challenge. We argue even basic galaxy properties, such as the stellar mass to halo mass, should be treated with errors bars of $\sim0.2-0.4$ dex.
△ Less
Submitted 10 February, 2016; v1 submitted 25 November, 2015;
originally announced November 2015.
-
Geometrical on-the-fly shock detection in SPH
Authors:
Alexander M. Beck,
K. Dolag,
Julius M. F. Donnert
Abstract:
We present an on-the-fly geometrical approach for shock detection and Mach number calculation in simulations employing smoothed particle hydrodynamics (SPH). We utilize pressure gradients to select shock candidates and define up- and downstream positions. We obtain hydrodynamical states in the up- and downstream regimes with a series of normal and inverted kernel weightings parallel and perpendicu…
▽ More
We present an on-the-fly geometrical approach for shock detection and Mach number calculation in simulations employing smoothed particle hydrodynamics (SPH). We utilize pressure gradients to select shock candidates and define up- and downstream positions. We obtain hydrodynamical states in the up- and downstream regimes with a series of normal and inverted kernel weightings parallel and perpendicular to the shock normals. Our on-the-fly geometrical Mach detector incorporates well within the SPH formalism and has low computational cost.
We implement our Mach detector into the simulation code GADGET and alongside many SPH improvements. We test our shock finder in a sequence of shock-tube tests with successively increasing Mach numbers exceeding by far the typical values inside galaxy clusters. For all shocks, we resolve the shocks well and the correct Mach numbers are assigned. An application to a strong magnetized shock-tube gives stable results in full magnetohydrodynamic set-ups. We simulate a merger of two idealized galaxy clusters and study the shock front. Shock structures within the merging clusters as well as the cluster shock are well-captured by our algorithm and assigned correct Mach numbers.
△ Less
Submitted 26 February, 2016; v1 submitted 23 November, 2015;
originally announced November 2015.
-
nIFTy galaxy cluster simulations II: radiative models
Authors:
Federico Sembolini,
Pascal Jahan Elahi,
Frazer R. Pearce,
Chris Power,
Alexander Knebe,
Scott T. Kay,
Weiguang Cui,
Gustavo Yepes,
Alexander M. Beck,
Stefano Borgani,
Daniel Cunnama,
Romeel Davé,
Sean February,
Shuiyao Huang,
Neal Katz,
Ian G. McCarthy,
Giuseppe Murante,
Richard D. A. Newton,
Valentin Perret,
Alexandro Saro,
Joop Schaye,
Romain Teyssier
Abstract:
We have simulated the formation of a massive galaxy cluster (M$_{200}^{\rm crit}$ = 1.1$\times$10$^{15}h^{-1}M_{\odot}$) in a $Λ$CDM universe using 10 different codes (RAMSES, 2 incarnations of AREPO and 7 of GADGET), modeling hydrodynamics with full radiative subgrid physics. These codes include Smoothed-Particle Hydrodynamics (SPH), spanning traditional and advanced SPH schemes, adaptive mesh an…
▽ More
We have simulated the formation of a massive galaxy cluster (M$_{200}^{\rm crit}$ = 1.1$\times$10$^{15}h^{-1}M_{\odot}$) in a $Λ$CDM universe using 10 different codes (RAMSES, 2 incarnations of AREPO and 7 of GADGET), modeling hydrodynamics with full radiative subgrid physics. These codes include Smoothed-Particle Hydrodynamics (SPH), spanning traditional and advanced SPH schemes, adaptive mesh and moving mesh codes. Our goal is to study the consistency between simulated clusters modeled with different radiative physical implementations - such as cooling, star formation and AGN feedback. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. We find that, with respect to non-radiative simulations, dark matter is more centrally concentrated, the extent not simply depending on the presence/absence of AGN feedback. The scatter in global quantities is substantially higher than for non-radiative runs. Intriguingly, adding radiative physics seems to have washed away the marked code-based differences present in the entropy profile seen for non-radiative simulations in Sembolini et al. (2015): radiative physics + classic SPH can produce entropy cores. Furthermore, the inclusion/absence of AGN feedback is not the dividing line -as in the case of describing the stellar content- for whether a code produces an unrealistic temperature inversion and a falling central entropy profile. However, AGN feedback does strongly affect the overall stellar distribution, limiting the effect of overcooling and reducing sensibly the stellar fraction.
△ Less
Submitted 11 November, 2015;
originally announced November 2015.
-
Neutral hydrogen in galaxy clusters: impact of AGN feedback and implications for intensity mapping
Authors:
Francisco Villaescusa-Navarro,
Susana Planelles,
Stefano Borgani,
Matteo Viel,
Elena Rasia,
Giuseppe Murante,
Klaus Dolag,
Lisa K. Steinborn,
Veronica Biffi,
Alexander M. Beck,
Cinthia Ragone-Figueroa
Abstract:
By means of zoom-in hydrodynamic simulations we quantify the amount of neutral hydrogen (HI) hosted by groups and clusters of galaxies. Our simulations, which are based on an improved formulation of smoothed particle hydrodynamics (SPH), include radiative cooling, star formation, metal enrichment and supernova feedback, and can be split in two different groups, depending on whether feedback from a…
▽ More
By means of zoom-in hydrodynamic simulations we quantify the amount of neutral hydrogen (HI) hosted by groups and clusters of galaxies. Our simulations, which are based on an improved formulation of smoothed particle hydrodynamics (SPH), include radiative cooling, star formation, metal enrichment and supernova feedback, and can be split in two different groups, depending on whether feedback from active galactic nuclei (AGN) is turned on or off. Simulations are analyzed to account for HI self-shielding and the presence of molecular hydrogen. We find that the mass in neutral hydrogen of dark matter halos monotonically increases with the halo mass and can be well described by a power-law of the form $M_{\rm HI}(M,z)\propto M^{3/4}$. Our results point out that AGN feedback reduces both the total halo mass and its HI mass, although it is more efficient in removing HI. We conclude that AGN feedback reduces the neutral hydrogen mass of a given halo by $\sim50\%$, with a weak dependence on halo mass and redshift. The spatial distribution of neutral hydrogen within halos is also affected by AGN feedback, whose effect is to decrease the fraction of HI that resides in the halo inner regions. By extrapolating our results to halos not resolved in our simulations we derive astrophysical implications from the measurements of $Ω_{\rm HI}(z)$: halos with circular velocities larger than $\sim25~{\rm km/s}$ are needed to host HI in order to reproduce observations. We find that only the model with AGN feedback is capable of reproducing the value of $Ω_{\rm HI}b_{\rm HI}$ derived from available 21cm intensity mapping observations.
△ Less
Submitted 8 December, 2015; v1 submitted 14 October, 2015;
originally announced October 2015.
-
Cool Core Clusters from Cosmological Simulations
Authors:
E. Rasia,
S. Borgani,
G. Murante,
S. Planelles,
A. M. Beck,
V. Biffi,
C. Ragone-Figueroa,
G. L. Granato,
L. K. Steinborn,
K. Dolag
Abstract:
We present results obtained from a set of cosmological hydrodynamic simulations of galaxy clusters, aimed at comparing predictions with observational data on the diversity between cool-core (CC) and non-cool-core (NCC) clusters. Our simulations include the effects of stellar and AGN feedback and are based on an improved version of the smoothed particle hydrodynamics code GADGET-3, which ameliorate…
▽ More
We present results obtained from a set of cosmological hydrodynamic simulations of galaxy clusters, aimed at comparing predictions with observational data on the diversity between cool-core (CC) and non-cool-core (NCC) clusters. Our simulations include the effects of stellar and AGN feedback and are based on an improved version of the smoothed particle hydrodynamics code GADGET-3, which ameliorates gas mixing and better captures gas-dynamical instabilities by including a suitable artificial thermal diffusion. In this Letter, we focus our analysis on the entropy profiles, the primary diagnostic we used to classify the degree of cool-coreness of clusters, and on the iron profiles. In keeping with observations, our simulated clusters display a variety of behaviors in entropy profiles: they range from steadily decreasing profiles at small radii, characteristic of cool-core systems, to nearly flat core isentropic profiles, characteristic of non-cool-core systems. Using observational criteria to distinguish between the two classes of objects, we find that they occur in similar proportions in both simulations and in observations. Furthermore, we also find that simulated cool-core clusters have profiles of iron abundance that are steeper than those of NCC clusters, which is also in agreement with observational results. We show that the capability of our simulations to generate a realistic cool-core structure in the cluster population is due to AGN feedback and artificial thermal diffusion: their combined action allows us to naturally distribute the energy extracted from super-massive black holes and to compensate for the radiative losses of low-entropy gas with short cooling time residing in the cluster core.
△ Less
Submitted 31 October, 2015; v1 submitted 14 September, 2015;
originally announced September 2015.
-
nIFTy galaxy cluster simulations I: dark matter & non-radiative models
Authors:
Federico Sembolini,
Gustavo Yepes,
Frazer R. Pearce,
Alexander Knebe,
Scott T. Kay,
Chris Power,
Weiguang Cui,
Alexander M. Beck,
Stefano Borgani,
Claudio Dalla Vecchia,
Romeel Davé,
Pascal Jahan Elahi,
Sean February,
Shuiyao Huang,
Alex Hobbs,
Neal Katz,
Erwin Lau,
Ian G. McCarthy,
Giuseppe Murante,
Daisuke Nagai,
Kaylea Nelson,
Richard D. A. Newton,
Ewald Puchwein,
Justin I. Read,
Alexandro Saro
, et al. (2 additional authors not shown)
Abstract:
We have simulated the formation of a galaxy cluster in a $Λ$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanc…
▽ More
We have simulated the formation of a galaxy cluster in a $Λ$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanced smoothed-particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ART and AREPO form extended entropy cores in the gas with rising central gas temperatures. Those codes employing traditional SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid based methods.
△ Less
Submitted 20 March, 2015;
originally announced March 2015.
-
Connecting Angular Momentum and Galactic Dynamics: The complex Interplay between Spin, Mass, and Morphology
Authors:
A. F. Teklu,
R. -S. Remus,
K. Dolag,
A. M. Beck,
A. Burkert,
A. S. Schmidt,
F. Schulze,
L. K. Steinborn
Abstract:
The evolution and distribution of the angular momentum of dark matter (DM) halos have been discussed in several studies over the past decades. In particular, the idea arose that angular momentum conservation should allow to infer the total angular momentum of the entire DM halo from measuring the angular momentum of the baryonic component, which is populating the center of the halo, especially for…
▽ More
The evolution and distribution of the angular momentum of dark matter (DM) halos have been discussed in several studies over the past decades. In particular, the idea arose that angular momentum conservation should allow to infer the total angular momentum of the entire DM halo from measuring the angular momentum of the baryonic component, which is populating the center of the halo, especially for disk galaxies. To test this idea and to understand the connection between the angular momentum of the DM halo and its galaxy, we use the Magneticum simulations. We successfully produce populations of spheroidal and disk galaxies self-consistently. Thus, we are able to study the dependence of galactic properties on their morphology. We find that (1) the specific angular momentum of stars in disk and spheroidal galaxies as a function of their stellar mass compares well with observational results; (2) the specific angular momentum of the stars in disk galaxies is slightly smaller compared to the specific angular momentum of the cold gas, in good agreement with observations; (3) simulations including the baryonic component show a dichotomy in the specific stellar angular momentum distribution when splitting the galaxies according to their morphological type (this dichotomy can also be seen in the spin parameter, where disk galaxies populate halos with slightly larger spin compared to spheroidal galaxies); (4) disk galaxies preferentially populate halos in which the angular momentum vector of the DM component in the central part shows a better alignment to the angular momentum vector of the entire halo; and (5) the specific angular momentum of the cold gas in disk galaxies is approximately 40 percent smaller than the specific angular momentum of the total DM halo and shows a significant scatter.
△ Less
Submitted 28 September, 2015; v1 submitted 11 March, 2015;
originally announced March 2015.
-
An improved SPH scheme for cosmological simulations
Authors:
A. M. Beck,
G. Murante,
A. Arth,
R. -S. Remus,
A. F. Teklu,
J. M. F. Donnert,
S. Planelles,
M. C. Beck,
P. Foerster,
M. Imgrund,
K. Dolag,
S. Borgani
Abstract:
We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement techniques in the latest instalment of the GADGET code. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-depend…
▽ More
We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement techniques in the latest instalment of the GADGET code. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-dependent scheme for artificial viscosity, which includes a high-order gradient computation and shear flow limiter. Additionally, we include a novel prescription for time-dependent artificial conduction, which corrects for gravitationally induced pressure gradients and largely improves the SPH performance in capturing the development of gas-dynamical instabilities. We extensively test our new implementation in a wide range of hydrodynamical standard tests including weak and strong shocks as well as shear flows, turbulent spectra, gas mixing, hydrostatic equilibria and self-gravitating gas clouds. We jointly employ all modifications; however, when necessary we study the performance of individual code modules. We approximate hydrodynamical states more accurately and with significantly less noise than standard SPH. Furthermore, the new implementation promotes the mixing of entropy between different fluid phases, also within cosmological simulations. Finally, we study the performance of the hydrodynamical solver in the context of radiative galaxy formation and non-radiative galaxy cluster formation. We find galactic disks to be colder, thinner and more extended and our results on galaxy clusters show entropy cores instead of steadily declining entropy profiles. In summary, we demonstrate that our improved SPH implementation overcomes most of the undesirable limitations of standard SPH, thus becoming the core of an efficient code for large cosmological simulations.
△ Less
Submitted 20 October, 2015; v1 submitted 25 February, 2015;
originally announced February 2015.
-
Anisotropic thermal conduction in galaxy clusters with MHD in Gadget
Authors:
Alexander Arth,
Klaus Dolag,
Alexander M. Beck,
Margarita Petkova,
Harald Lesch
Abstract:
We present an implementation of thermal conduction including the anisotropic effects of magnetic fields for SPH. The anisotropic thermal conduction is mainly proceeding parallel to magnetic fields and suppressed perpendicular to the fields. We derive the SPH formalism for the anisotropic heat transport and solve the corresponding equation with an implicit conjugate gradient scheme. We discuss seve…
▽ More
We present an implementation of thermal conduction including the anisotropic effects of magnetic fields for SPH. The anisotropic thermal conduction is mainly proceeding parallel to magnetic fields and suppressed perpendicular to the fields. We derive the SPH formalism for the anisotropic heat transport and solve the corresponding equation with an implicit conjugate gradient scheme. We discuss several issues of unphysical heat transport in the cases of extreme ansiotropies or unmagnetized regions and present possible numerical workarounds. We implement our algorithm into the GADGET code and study its behaviour in several test cases. In general, we reproduce the analytical solutions of our idealised test problems, and obtain good results in cosmological simulations of galaxy cluster formations. Within galaxy clusters, the anisotropic conduction produces a net heat transport similar to an isotropic Spitzer conduction model with an efficiency of one per cent. In contrast to isotropic conduction our new formalism allows small-scale structure in the temperature distribution to remain stable, because of their decoupling caused by magnetic field lines. Compared to observations, isotropic conduction with more than 10 per cent of the Spitzer value leads to an oversmoothed temperature distribution within clusters, while the results obtained with anisotropic thermal conduction reproduce the observed temperature fluctuations well. A proper treatment of heat transport is crucial especially in the outskirts of clusters and also in high density regions. It's connection to the local dynamical state of the cluster also might contribute to the observed bimodal distribution of cool core and non cool core clusters. Our new scheme significantly advances the modelling of thermal conduction in numerical simulations and overall gives better results compared to observations.
△ Less
Submitted 15 September, 2017; v1 submitted 19 December, 2014;
originally announced December 2014.
-
Constraints on the distribution and energetics of fast radio bursts using cosmological hydrodynamic simulations
Authors:
K. Dolag,
B. M. Gaensler,
A. M. Beck,
M. C. Beck
Abstract:
We present constraints on the origins of fast radio bursts (FRBs) using large cosmological simulations. We calculate contributions to FRB dispersion measures (DMs) from the Milky Way, from the local Universe, from cosmological large-scale structure, and from potential FRB host galaxies, and then compare these simulations to the DMs of observed FRBs. We find that the Milky Way contribution has prev…
▽ More
We present constraints on the origins of fast radio bursts (FRBs) using large cosmological simulations. We calculate contributions to FRB dispersion measures (DMs) from the Milky Way, from the local Universe, from cosmological large-scale structure, and from potential FRB host galaxies, and then compare these simulations to the DMs of observed FRBs. We find that the Milky Way contribution has previously been underestimated by a factor of ~2, and that the foreground-subtracted DMs are consistent with a cosmological origin, corresponding to a source population observable to a maximum redshift z~0.6-0.9. We consider models for the spatial distribution of FRBs in which they are randomly distributed in the Universe, track the star-formation rate of their host galaxies, track total stellar mass, or require a central supermassive black hole. Current data do not discriminate between these possibilities, but the predicted DM distributions for different models will differ considerably once we begin detecting FRBs at higher DMs and higher redshifts. We additionally consider the distribution of FRB fluences, and show that the observations are consistent with FRBs being standard candles, each burst producing the same radiated isotropic energy. The data imply a constant isotropic burst energy of ~7e40 erg if FRBs are embedded in host galaxies, or ~9e40 erg if FRBs are randomly distributed. These energies are 10-100 times larger than had previously been inferred. Within the constraints of the available small sample of data, our analysis favours FRB mechanisms for which the isotropic radiated energy has a narrow distribution in excess of 1e40 erg.
△ Less
Submitted 22 May, 2015; v1 submitted 15 December, 2014;
originally announced December 2014.
-
New constraints on modelling the random magnetic field of the MW
Authors:
Marcus C. Beck,
Alexander M. Beck,
Rainer Beck,
Klaus Dolag,
Andrew W. Strong,
Peter Nielaba
Abstract:
We extend the description of the isotropic and anisotropic random component of the small-scale magnetic field within the existing magnetic field model of the Milky Way from Jansson&Farrar, by including random realizations of the small-scale component. Using a magnetic-field power spectrum with Gaussian random fields, the NE2001 model for the thermal electrons and the Galactic cosmic-ray electron d…
▽ More
We extend the description of the isotropic and anisotropic random component of the small-scale magnetic field within the existing magnetic field model of the Milky Way from Jansson&Farrar, by including random realizations of the small-scale component. Using a magnetic-field power spectrum with Gaussian random fields, the NE2001 model for the thermal electrons and the Galactic cosmic-ray electron distribution from the current GALPROP model we derive full-sky maps for the total and polarized synchrotron intensity as well as the Faraday rotation-measure distribution. While previous work assumed that small-scale fluctuations average out along the line-of-sight or which only computed ensemble averages of random fields, we show that these fluctuations need to be carefully taken into account. Comparing with observational data we obtain not only good agreement with 408 MHz total and WMAP7 22 GHz polarized intensity emission maps, but also an improved agreement with Galactic foreground rotation-measure maps and power spectra, whose amplitude and shape strongly depend on the parameters of the random field. We demonstrate that a correlation length of ~220pc (50pc being a 5σlower limit) is needed to match the slope of the observed power spectrum of Galactic foreground rotation-measure maps. Using multiple realizations allows us also to infer errors on individual observables. We find that previously-used amplitudes for random and anisotropic random magnetic field components need to be rescaled by factors of ~0.3 and 0.6 to account for the new small-scale contributions. Our model predicts a rotation measure of -2.8\pm7.1 rad/m^2 and 4.4\pm11.0 rad/m^2 for the north and south Galactic poles respectively, in good agreement with observations. Applying our model to deflections of ultra-high-energy cosmic rays we infer a mean deflection of ~3.5\pm1.1 degree for 60 EeV protons arriving from CenA.
△ Less
Submitted 1 June, 2016; v1 submitted 17 September, 2014;
originally announced September 2014.
-
Strong magnetic fields and large rotation measures in protogalaxies by supernova seeding
Authors:
Alexander M. Beck,
Klaus Dolag,
Harald Lesch,
Philipp P. Kronberg
Abstract:
We present a model for the seeding and evolution of magnetic fields in protogalaxies. Supernova (SN) explosions during the assembly of a protogalaxy provide magnetic seed fields, which are subsequently amplified by compression, shear flows and random motions. We implement the model into the MHD version of the cosmological N-body / SPH simulation code GADGET and we couple the magnetic seeding direc…
▽ More
We present a model for the seeding and evolution of magnetic fields in protogalaxies. Supernova (SN) explosions during the assembly of a protogalaxy provide magnetic seed fields, which are subsequently amplified by compression, shear flows and random motions. We implement the model into the MHD version of the cosmological N-body / SPH simulation code GADGET and we couple the magnetic seeding directly to the underlying multi-phase description of star formation. We perform simulations of Milky Way-like galactic halo formation using a standard LCDM cosmology and analyse the strength and distribution of the subsequent evolving magnetic field. A dipole-shape divergence-free magnetic field is injected at a rate of 10^{-9}G / Gyr within starforming regions, given typical dimensions and magnetic field strengths in canonical SN remnants. Subsequently, the magnetic field strength increases exponentially on timescales of a few ten million years. At redshift z=0, the entire galactic halo is magnetized and the field amplitude is of the order of a few $μ$G in the center of the halo, and 10^{-9} G at the virial radius. Additionally, we analyse the intrinsic rotation measure (RM) of the forming galactic halo over redshift. The mean halo intrinsic RM peaks between redshifts z=4 and z=2 and reaches absolute values around 1000 rad m^{-2}. While the halo virializes towards redshift z=0, the intrinsic RM values decline to a mean value below 10 rad m^{-2}. At high redshifts, the distribution of individual starforming, and thus magnetized regions is widespread. In our model for the evolution of galactic magnetic fields, the seed magnetic field amplitude and distribution is no longer a free parameter, but determined self-consistently by the star formation process occuring during the formation of cosmic structures.
△ Less
Submitted 15 August, 2013;
originally announced August 2013.
-
On the magnetic fields in voids
Authors:
Alexander M. Beck,
Michal Hanasz,
Harald Lesch,
Rhea-Silvia Remus,
Federico A. Stasyszyn
Abstract:
We study the possible magnetization of cosmic voids by void galaxies. Recently, observations revealed isolated starforming galaxies within the voids. Furthermore, a major fraction of a voids volume is expected to be filled with magnetic fields of a minimum strength of about $10^{-15}$ G on Mpc scales. We estimate the transport of magnetic energy by cosmic rays (CR) from the void galaxies into the…
▽ More
We study the possible magnetization of cosmic voids by void galaxies. Recently, observations revealed isolated starforming galaxies within the voids. Furthermore, a major fraction of a voids volume is expected to be filled with magnetic fields of a minimum strength of about $10^{-15}$ G on Mpc scales. We estimate the transport of magnetic energy by cosmic rays (CR) from the void galaxies into the voids. We assume that CRs and winds are able to leave small isolated void galaxies shortly after they assembled, and then propagate within the voids. For a typical void, we estimate the magnetic field strength and volume filling factor depending on its void galaxy population and possible contributions of strong active galactic nuclei (AGN) which border the voids. We argue that the lower limit on the void magnetic field can be recovered, if a small fraction of the magnetic energy contained in the void galaxies or void bordering AGNs is distributed within the voids.
△ Less
Submitted 31 October, 2012;
originally announced October 2012.
-
Synthetic X-ray and radio maps for two different models of Stephan's Quintet
Authors:
Annette Geng,
Alexander M. Beck,
Klaus Dolag,
Florian Bürzle,
Marcus C. Beck,
Hanna Kotarba,
Peter Nielaba
Abstract:
We present simulations of the compact galaxy group Stephan's Quintet (SQ) including magnetic fields, performed with the N-body/smoothed particle hydrodynamics (SPH) code \textsc{Gadget}. The simulations include radiative cooling, star formation and supernova feedback. Magnetohydrodynamics (MHD) is implemented using the standard smoothed particle magnetohydrodynamics (SPMHD) method. We adapt two di…
▽ More
We present simulations of the compact galaxy group Stephan's Quintet (SQ) including magnetic fields, performed with the N-body/smoothed particle hydrodynamics (SPH) code \textsc{Gadget}. The simulations include radiative cooling, star formation and supernova feedback. Magnetohydrodynamics (MHD) is implemented using the standard smoothed particle magnetohydrodynamics (SPMHD) method. We adapt two different initial models for SQ based on Renaud et al. and Hwang et al., both including four galaxies (NGC 7319, NGC 7320c, NGC 7318a and NGC 7318b). Additionally, the galaxies are embedded in a magnetized, low density intergalactic medium (IGM). The ambient IGM has an initial magnetic field of $10^{-9}$ G and the four progenitor discs have initial magnetic fields of $10^{-9} - 10^{-7}$ G. We investigate the morphology, regions of star formation, temperature, X-ray emission, magnetic field structure and radio emission within the two different SQ models. In general, the enhancement and propagation of the studied gaseous properties (temperature, X-ray emission, magnetic field strength and synchrotron intensity) is more efficient for the SQ model based on Renaud et al., whose galaxies are more massive, whereas the less massive SQ model based on Hwang et al. shows generally similar effects but with smaller efficiency. We show that the large shock found in observations of SQ is most likely the result of a collision of the galaxy NGC 7318b with the IGM. This large group-wide shock is clearly visible in the X-ray emission and synchrotron intensity within the simulations of both SQ models. The order of magnitude of the observed synchrotron emission within the shock front is slightly better reproduced by the SQ model based on Renaud et al., whereas the distribution and structure of the synchrotron emission is better reproduced by the SQ model based on Hwang et al..
△ Less
Submitted 14 August, 2012; v1 submitted 6 June, 2012;
originally announced June 2012.
-
A divergence-cleaning scheme for cosmological SPMHD simulations
Authors:
F. A. Stasyszyn,
K. Dolag,
A. M. Beck
Abstract:
In magnetohydrodynamics (MHD), the magnetic field is evolved by the induction equation and coupled to the gas dynamics by the Lorentz force. We perform numerical smoothed particle magnetohydrodynamics (Spmhd) simulations and study the influence of a numerical magnetic divergence. For instabilities arising from divergence B related errors, we find the hyperbolic/parabolic cleaning scheme suggested…
▽ More
In magnetohydrodynamics (MHD), the magnetic field is evolved by the induction equation and coupled to the gas dynamics by the Lorentz force. We perform numerical smoothed particle magnetohydrodynamics (Spmhd) simulations and study the influence of a numerical magnetic divergence. For instabilities arising from divergence B related errors, we find the hyperbolic/parabolic cleaning scheme suggested by Dedner et al. 2002 to give good results and prevent numerical artifacts from growing. Additionally, we demonstrate that certain current Spmhd implementations of magnetic field regularizations give rise to unphysical instabilities in long-time simulations. We also find this effect when employing Euler potentials (divergenceless by definition), which are not able to follow the winding-up process of magnetic field lines properly. Furthermore, we present cosmological simulations of galaxy cluster formation at extremely high resolution including the evolution of magnetic fields. We show synthetic Faraday rotation maps and derive structure functions to compare them with observations. Comparing all the simulations with and without divergence cleaning, we are able to confirm the results of previous simulations performed with the standard implementation of MHD in Spmhd at normal resolution. However, at extremely high resolution, a cleaning scheme is needed to prevent the growth of numerical errors at small scales.
△ Less
Submitted 18 May, 2012;
originally announced May 2012.
-
Origin of strong magnetic fields in Milky-Way like galactic haloes
Authors:
Alexander M. Beck,
Harald Lesch,
Klaus Dolag,
Hanna Kotarba,
Annette Geng,
Federico A. Stasyszyn
Abstract:
An analytical model predicting the growth rates, the absolute growth times and the saturation values of the magnetic field strength within galactic haloes is presented. The analytical results are compared to cosmological MHD simulations of Milky-Way like galactic halo formation performed with the N-body / \textsc{Spmhd} code \textsc{Gadget}. The halo has a mass of $\approx{}3\cdot{}10^{12}$…
▽ More
An analytical model predicting the growth rates, the absolute growth times and the saturation values of the magnetic field strength within galactic haloes is presented. The analytical results are compared to cosmological MHD simulations of Milky-Way like galactic halo formation performed with the N-body / \textsc{Spmhd} code \textsc{Gadget}. The halo has a mass of $\approx{}3\cdot{}10^{12}$ $M_{\odot}$ and a virial radius of $\approx{}$270 kpc. The simulations in a $Λ$CDM cosmology also include radiative cooling, star formation, supernova feedback and the description of non-ideal MHD. A primordial magnetic seed field ranging from $10^{-10}$ to $10^{-34}$ G in strength agglomerates together with the gas within filaments and protohaloes. There, it is amplified within a couple of hundred million years up to equipartition with the corresponding turbulent energy. The magnetic field strength increases by turbulent small-scale dynamo action. The turbulence is generated by the gravitational collapse and by supernova feedback. Subsequently, a series of halo mergers leads to shock waves and amplification processes magnetizing the surrounding gas within a few billion years. At first, the magnetic energy grows on small scales and then self-organizes to larger scales. Magnetic field strengths of $\approx{}10^{-6}$ G are reached in the center of the halo and drop to $\approx{}10^{-9}$ G in the IGM. Analyzing the saturation levels and growth rates, the model is able to describe the process of magnetic amplification notably well and confirms the results of the simulations.
△ Less
Submitted 15 February, 2012;
originally announced February 2012.
