We examine the stellar haloes of the Auriga simulations, a suite of thirty cosmological magneto-hydrodynamical high-resolution simulations of Milky Way-mass galaxies performed with the moving-mesh code arepo. We study halo global properties and radial profiles out to ∼ 150 kpc for each individual galaxy. The Auriga haloes are diverse in their masses and density profiles; mean metallicity and metallicity gradients; ages; and shapes, reflecting the stochasticity inherent in their accretion and merger histories. A comparison with observations of nearby late-type galaxies shows very good agreement between most observed and simulated halo properties. However, Auriga haloes are typically too massive. We find a connection between population gradients and mass assembly history: galaxies with few significant progenitors have more massive haloes, possess large negative halo metallicity gradients and steeper density profiles. The number of accreted galaxies, either disrupted or under disruption, that contribute 90% of the accreted halo mass ranges from 1 to 14, with a median of 6.5, and their stellar masses span over three orders of magnitude. The observed halo mass-metallicity relation is well reproduced by Auriga and is set by the stellar mass and metallicity of the dominant satellite contributors. This relationship is found not only for the accreted component but also for the total (accreted + in-situ) stellar halo. Our results highlight the potential of observable halo properties to infer the assembly history of galaxies.

The stellar disk of the Milky Way shows complex spatial and abundance structure that is central to understanding the key physical mechanisms responsible for shaping our Galaxy. In this study, we use six very high resolution cosmological zoom simulations of Milky Way-sized haloes to study the prevalence and formation of chemically distinct disc components. We find that our simulations develop a clearly bimodal distribution in the [α/Fe] -[Fe/H] plane. We find two main pathways to creating this dichotomy which operate in different regions of the galaxies: a) an early (z > 1) and intense high-[α/Fe] star formation phase in the inner region (R 5 kpc) induced by gas-rich mergers, followed by more quiescent low-[α/Fe] star formation; and b) an early phase of high-[α/Fe] star formation in the outer disc followed by a shrinking of the gas disc owing to a temporarily lowered gas accretion rate, after which disc growth resumes. In process b), a double-peaked star formation history around the time and radius of disc shrinking accentuates the dichotomy. If the early star formation phase is prolonged (rather than short and intense), chemical evolution proceeds as per process a) in the inner region, but the dichotomy is less clear. In the outer region, the dichotomy is only evident if the first intense phase of star formation covers a large enough radial range before disc shrinking occurs; otherwise, the outer disc consists of only low-[α/Fe] sequence stars. We discuss the implication that both processes occurred in the Milky Way.

Observational studies have revealed that galaxy pairs tend to have lower gas-phase metallicity than isolated galaxies. This metallicity deficiency can be caused by inflows of low-metallicity gas due to the tidal forces and gravitational torques associated with galaxy mergers, diluting the metal content of the central region. In this work we demonstrate that such metallicity dilution occurs in state-of-the-art cosmological simulations of galaxy formation. We find that the dilution is typically 0.1 dex for major mergers, and is noticeable at projected separations smaller than 40 kpc. For minor mergers the metallicity dilution is still present, even though the amplitude is significantly smaller. Consistent with previous analysis of observed galaxies we find that mergers are outliers from the fundamental metallicity relation, with deviations being larger than expected for a Gaussian distribution of residuals. Our large sample of mergers within full cosmological simulations also makes it possible to estimate how the star formation rate enhancement and gas consumption timescale behave as a function of the merger mass ratio. We confirm that strong starbursts are likely to occur in major mergers, but they can also arise in minor mergers if more than two galaxies are participating in the interaction, a scenario that has largely been ignored in previous work based on idealised isolated merger simulations.

Magnetic protection of potentially habitable planets plays a central role in determining their actual habitability and/or the chances of detecting atmospheric biosignatures. We develop here a thermal evolution model of potentially habitable Earth-like planets and super-Earths. Using up-to-date dynamo scaling laws we predict the properties of core dynamo magnetic fields and study the influence of thermal evolution on their properties. The level of magnetic protection of tidally locked and unlocked planets is estimated by combining simplified models of the planetary magnetosphere and a phenomenological description of the stellar wind. Thermal evolution introduces a strong dependence of magnetic protection on planetary mass and rotation rate. Tidally locked terrestrial planets with an Earth-like composition would have early dayside magnetospause distances between 1.5 and 4.0 R p , larger than previously estimated. Unlocked planets with periods of rotation ∼ 1 day are protected by magnetospheres extending between 3 and 8 R p . Our results are robust against variations in planetary bulk composition and uncertainties in other critical model parameters. For illustration purposes the thermal evolution and magnetic protection of the potentially habitable super-Earths GL 581d, GJ 667Cc and HD 40307g were also studied. Assuming an Earth-like composition we found that the dynamos of these planets are already extinct or close to being shut down. While GL 581d is the best protected, the protection of HD 40307g cannot be reliably estimated. GJ 667Cc, even under optimistic conditions, seems to be severely exposed to the stellar wind and, under the conditions of our model, has probably suffered massive atmospheric losses.

Citation for published item:q¡ omezD pundo eF nd qrndD oert tF tF nd wonhesiD entonel nd hiteD imon hF wF nd fustmnteD estin nd wriniD pederio nd kmorD ¤ udiger nd impsonD ghristine wF nd pringelD olker nd prenkD grlos F @PHIUA 9vessons from the eurig diss X the hunt for the wilky y9s ex situ dis is not yet overF9D wonthly xoties of the oyl estronomil oietyFD RUP @QAF ppF QUPPEQUQQF Further information on publisher's website: Additional information: Use policyThe full-text may be used and/or reproduced, and given to third parties in any format or medium, without prior permission or charge, for personal research or study, educational, or not-for-prot purposes provided that:• a full bibliographic reference is made to the original source • a link is made to the metadata record in DRO • the full-text is not changed in any way The full-text must not be sold in any format or medium without the formal permission of the copyright holders.Please consult the full DRO policy for further details. ABSTRACTWe characterize the contribution from accreted material to the galactic discs of the Auriga Project, a set of high-resolution magnetohydrodynamic cosmological simulations of late-type galaxies performed with the moving-mesh code AREPO. Our goal is to explore whether a significant accreted (or ex situ) stellar component in the Milky Way disc could be hidden within the near-circular orbit population, which is strongly dominated by stars born in situ.One-third of our models shows a significant ex situ disc but this fraction would be larger if constraints on orbital circularity were relaxed. Most of the ex situ material ( 50 per cent) comes from single massive satellites (>6 × 10 10 M ). These satellites are accreted with a wide range of infall times and inclination angles (up to 85• ). Ex situ discs are thicker, older and more metal poor than their in situ counterparts. They show a flat median age profile, which differs from the negative gradient observed in the in situ component. As a result, the likelihood of identifying an ex situ disc in samples of old stars on near-circular orbits increases towards the outskirts of the disc. We show three examples that, in addition to ex situ discs, have a strongly rotating dark matter component. Interestingly, two of these ex situ stellar discs show an orbital circularity distribution that is consistent with that of the in situ disc. Thus, they would not be detected in typical kinematic studies.