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. Its 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.
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 $mu$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.
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 m
inimum 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.
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{}3cdot{}10^{12}$ $M_{odot}$ and a virial radius of $approx{}$270 kpc. The simulations in a $Lambda$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.
