No Arabic abstract
The magnetic fields observed in the Milky~Way and nearby galaxies appear to be in equipartition with the turbulent, thermal, and cosmic ray energy densities, and hence are expected to be dynamically important. However, the origin of these strong magnetic fields is still unclear, and most previous attempts to simulate galaxy formation from cosmological initial conditions have ignored them altogether. Here, we analyse the magnetic fields predicted by the simulations of the Auriga Project, a set of 30 high-resolution cosmological zoom simulations of Milky~Way-like galaxies, carried out with a moving-mesh magneto-hydrodynamics code and a detailed galaxy formation physics model. We find that the magnetic fields grow exponentially at early times owing to a small-scale dynamo with an e-folding time of roughly $100,rm{Myr}$ in the center of halos until saturation occurs around $z=2-3$, when the magnetic energy density reaches about $10%$ of the turbulent energy density with a typical strength of $10-50,rm{mu G}$. In the galactic centers the ratio between magnetic and turbulent energy remains nearly constant until $z=0$. At larger radii, differential rotation in the disks leads to linear amplification that typically saturates around $z=0.5$ to $z=0$. The final radial and vertical variations of the magnetic field strength can be well described by two joint exponential profiles, and are in good agreement with observational constraints. Overall, the magnetic fields have only little effect on the global evolution of the galaxies as it takes too long to reach equipartition. We also demonstrate that our results are well converged with numerical resolution.
We introduce a dust model for cosmological simulations implemented in the moving-mesh code AREPO and present a suite of cosmological hydrodynamical zoom-in simulations to study dust formation within galactic haloes. Our model accounts for the stellar production of dust, accretion of gas-phase metals onto existing grains, destruction of dust through local supernova activity, and dust driven by winds from star-forming regions. We find that accurate stellar and active galactic nuclei feedback is needed to reproduce the observed dust-metallicity relation and that dust growth largely dominates dust destruction. Our simulations predict a dust content of the interstellar medium which is consistent with observed scaling relations at $z = 0$, including scalings between dust-to-gas ratio and metallicity, dust mass and gas mass, dust-to-gas ratio and stellar mass, and dust-to-stellar mass ratio and gas fraction. We find that roughly two-thirds of dust at $z = 0$ originated from Type II supernovae, with the contribution from asymptotic giant branch stars below 20 per cent for $z gtrsim 5$. While our suite of Milky Way-sized galaxies forms dust in good agreement with a number of key observables, it predicts a high dust-to-metal ratio in the circumgalactic medium, which motivates a more realistic treatment of thermal sputtering of grains and dust cooling channels.
We present redshift-zero synthetic observational data considering dust attenuation and dust emission for the thirty galaxies of the Auriga project, calculated with the SKIRT radiative transfer code. The post-processing procedure includes components for star-forming regions, stellar sources, and diffuse dust taking into account stochastic heating of dust grains. This allows us to obtain realistic high-resolution broadband images and fluxes from ultraviolet to sub-millimeter wavelengths. For the diffuse dust component, we consider two mechanisms for assigning dust to gas cells in the simulation. In one case, only the densest or the coldest gas cells are allowed to have dust, while in the other case this condition is relaxed to allow a larger number of dust-containing cells. The latter approach yields galaxies with a larger radial dust extent and an enhanced dust presence in the inter-spiral regions. At a global scale, we compare Auriga galaxies with observations by deriving dust scaling relations using SED fitting. At a resolved scale, we make a multi-wavelength morphological comparison with nine well-resolved spiral galaxies from the DustPedia observational database. We find that for both dust assignment methods, although the Auriga galaxies show a good overall agreement with observational dust properties, they exhibit a slightly higher specific dust mass. The multi-wavelength morphological analysis reveals a good agreement between the Auriga and the observed galaxies in the optical wavelengths. In the mid and far-infrared wavelengths, Auriga galaxies appear smaller and more centrally concentrated in comparison to their observed counterparts. We publicly release the multi-observer images and fluxes in 50 commonly used broadband filters.
We present Rotation Measures (RM) of the diffuse Galactic synchrotron emission from the Canadian Galactic Plane Survey (CGPS) and compare them to RMs of extragalactic sources in order to study the large-scale reversal in the Galactic magnetic field (GMF). Using Stokes Q, U and I measurements of the Galactic disk collected with the Synthesis Telescope at the Dominion Radio Astrophysical Observatory, we calculate RMs over an extended region of the sky, focusing on the low longitude range of the CGPS (l=52deg to l=72deg). We note the similarity in the structures traced by the compact sources and the extended emission and highlight the presence of a gradient in the RM map across an approximately diagonal line, which we identify with the well-known field reversal of the Sagittarius-Carina arm. We suggest that the orientation of this reversal is a geometric effect resulting from our location within a GMF structure arising from current sheets that are not perpendicular to the Galactic plane, as is required for a strictly radial field reversal, but that have at least some component parallel to the disk. Examples of models that fit this description are the three-dimensional dynamo-based model of Gressel et al. (2013) and a Galactic scale Parker spiral (Akasofu & Hakamada 1982), although the latter may be problematic in terms of Galactic dynamics. We emphasize the importance of constructing three-dimensional models of the GMF to account for structures like the diagonal RM gradient observed in this dataset.
The magnetic fields of spiral galaxies are so strong that they cannot be primordial. Their typical values are over one billion times higher than any value predicted for the early Universe. Explaining this immense growth and incorporating it in galaxy evolution theories is one of the long-standing challenges in astrophysics. So far, the most successful theory for the sustained growth of the galactic magnetic field is the alpha-omega dynamo. This theory predicts a characteristic dipolar or quadrupolar morphology for the galactic magnetic field, which has been observed in external galaxies. However, so far, there has been no direct demonstration of a mean-field dynamo operating in direct, multi-physics simulations of spiral galaxies. We do so in this work. We employ numerical models of isolated, star-forming spiral galaxies that include a magnetized gaseous disk, a dark matter halo, stars, and stellar feedback. Naturally, the resulting magnetic field has a complex morphology that includes a strong random component. Using a smoothing of the magnetic field on small scales, we are able to separate the mean from the turbulent component and analyze them individually. We find that a mean-field dynamo naturally occurs as a result of the dynamical evolution of the galaxy and amplifies the magnetic field by an order of magnitude over half a Gyr. Despite the highly dynamical nature of these models, the morphology of the mean component of the field is identical to analytical predictions. This result underlines the importance of the mean-field dynamo in galactic evolution. Moreover, by demonstrating the natural growth of the magnetic field in a complex galactic environment, it brings us a step closer to understanding the cosmic origin of magnetic fields.
We analyze 494 main sequence turnoff and subgiant stars from the AMBRE:HARPS survey. These stars have accurate astrometric information from textit{Gaia}/DR1, providing reliable age estimates with relative uncertainties of $pm1-2$ Gyr and allowing precise orbital determinations. The sample is split based on chemistry into a low-[Mg/Fe] sequence, which are often identified as thin disk stellar populations, and a high-[Mg/Fe] sequence, which are often associated with the thick disk. We find that the high-[Mg/Fe] chemical sequence has extended star formation for several Gyr and is coeval with the oldest stars of the low-[Mg/Fe] chemical sequence: both the low- and high-[Mg/Fe] sequences were forming stars at the same time. The high-[Mg/Fe] stellar populations are only vertically extended for the oldest, most-metal poor and highest [Mg/Fe] stars. When comparing vertical velocity dispersion for both sequences, the high-[Mg/Fe] sequence has lower velocity dispersion than the low-[Mg/Fe] sequence for stars of similar age. Identifying either group as thin or thick disk based on chemistry is misleading. The stars belonging to the high-[Mg/Fe] sequence have perigalacticons that originate in the inner disk, while the perigalacticons of stars on the low-[Mg/Fe] sequence are generally around the solar neighborhood. From the orbital properties of the stars, the high-and low-[Mg/Fe] sequences are most likely a reflection of the chemical enrichment history of the inner and outer disk populations; radial mixing causes both populations to be observed in situ at the solar position. Based on these results, we emphasize that it is important to be clear in defining what populations are being referenced when using the terms thin and thick disk, and that ideally the term thick disk should be reserved for purely geometric definitions to avoid confusion and be consistent with definitions in external galaxies.