-
Magnetic buoyancy in simulated galactic discs with a realistic circum galactic medium
Authors:
Ulrich P. Steinwandel,
Marcus C. Beck,
Alexander Arth,
Klaus Dolag,
Benjamin P. Moster,
Peter Nielaba
Abstract:
We present simulations of isolated disc galaxies in a realistic environment performed with the Tree-SPMHD-Code Gadget-3. Our simulations include a spherical circum-galactic medium (CGM) surrounding the galactic disc, motivated by observations and the results of cosmological simulations. We present three galactic models with different halo masses between 10e10 Msol and 10e12 Msol, and for each we u…
▽ More
We present simulations of isolated disc galaxies in a realistic environment performed with the Tree-SPMHD-Code Gadget-3. Our simulations include a spherical circum-galactic medium (CGM) surrounding the galactic disc, motivated by observations and the results of cosmological simulations. We present three galactic models with different halo masses between 10e10 Msol and 10e12 Msol, and for each we use two different approaches to seed the magnetic field, as well as a control simulation without a magnetic field. We find that the amplification of the magnetic field in the centre of the disc leads to a biconical magnetic outflow of gas that magnetizes the CGM. This biconical magnetic outflow reduces the star formation rate (SFR) of the galaxy by roughly 40 percent compared to the simulations without magnetic fields. As the key aspect of our simulations, we find that small scale turbulent motion of the gas in the disc leads to the amplification of the magnetic field up to tens of 10e-6 G, as long as the magnetic field strength is low. For stronger magnetic fields turbulent motion does not lead to significant amplification but is replaced by an alpha-omega dynamo. The occurance of a small scale turbulent dynamo becomes apparent through the magnetic power spectrum and analysis of the field lines' curvature. In accordance with recent observations we find an anti-correlation between the spiral structure in the gas density and in the magnetic field due to a diffusion term added to the induction equation.
△ Less
Submitted 29 August, 2018;
originally announced August 2018.
-
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.
-
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.
-
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.