In recent observational studies, star-forming galaxies have been shown to follow a relation often dubbed the fundamental metallicity relation (FMR). This relation links the stellar mass of a galaxy with its star formation rate (SFR) and its gas-phase metallicity. Specifically, the FMR predicts that galaxies, at a given stellar mass, exhibit lower metallicities for higher SFRs. This trend is qualitatively consistent with observations of galaxy pairs, which have been robustly shown to experience increasing gas-phase metallicity dilution and enhanced star formation activity with decreasing projected separation. In this work, we show that, despite the qualitative consistency with FMR expectations, the observed O/H dilution in galaxy pairs of the Sloan Digital Sky Survey is stronger than what is predicted by the FMR. We conclude that the evolutionary phase of galaxies interacting with companions is not encoded in the FMR, and thus, mergers constitute a clearly defined population of outliers. We find that galaxies in pairs are consistent with the FMR only when their separation is larger than 110 kpc. Finally, we also quantify the local environment of the pairs using the number of galaxy neighbours within $2\, \mathrm{Mpc}$, $N_2$, and the projected separation to the second closest galaxy, $r_2$. We find that pairs are more sensitive to a second companion than to the local galaxy density, displaying less elevated SFRs with smaller values of $r_2$.

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 $\sim 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.

It is well established that the properties of supermassive black holes and their host galaxies are correlated through scaling relations. While hydrodynamical cosmological simulations have begun to account for the co-evolution of BHs and galaxies, they typically have neglected the BH spin, even though it may play an important role in modulating the growth and feedback of BHs. Here we introduce a new sub-grid model for the BH spin evolution in the moving-mesh code {\small AREPO} in order to improve the physical faithfulness of the BH modelling in galaxy formation simulations. We account for several different channels of spin evolution, in particular gas accretion through a Shakura-Sunyaev $\alpha$-disc, chaotic accretion, and BH mergers. For BH feedback, we extend the IllustrisTNG model, which considers two different BH feedback modes, a thermal quasar mode for high accretion states and a kinetic mode for low Eddington ratios, with a self-consistent accounting of spin-dependent radiative efficiencies and thus feedback strength. We find that BHs with mass $M_{\rm{bh}}\lesssim 10^{8}\, {\rm M}_\odot$ reach high spin values as they typically evolve in the coherent gas accretion regime. On the other hand, BHs with mass $M_{\rm{bh}}\gtrsim 10^{8}\, {\rm M}_\odot$ have lower spins as BH mergers become more frequent, and their accretion discs fragment due to self-gravity, inducing chaotic accretion. We also explore the hypothesis that the transition between the quasar and kinetic feedback modes is mediated by the accretion mode of the BH disc itself, i.e.~the kinetic feedback mode is activated when the disc enters the self-gravity regime. We find excellent agreement between the galaxy and BH populations for this approach and the fiducial TNG model with no spin evolution. Furthermore, our new approach alleviates a tension in the galaxy morphology -- colour relation of the original TNG model.

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 \emph{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.

The discovery of a planet around the closest star to our Sun, Proxima Centauri, represents a quantum leap in the testability of exoplanetary models. Unlike any other discovered exoplanet, models of Proxima b could be contrasted against near future telescopic observations and far future in-situ measurements. In this paper we aim at predicting the planetary radius and the magnetic properties (dynamo lifetime and magnetic dipole moment) of Proxima b analogues (solid planets with masses of $\sim 1-3\,M_\oplus$, rotation periods of several days and habitable conditions). For this purpose we build a grid of planetary models with a wide range of compositions and masses. For each point in the grid we run the planetary evolution model developed in \citet{Zuluaga2013}. Our model assumes small orbital eccentricity, negligible tidal heating and earth-like radiogenic mantle elements abundances. We devise a statistical methodology to estimate the posterior distribution of the desired planetary properties assuming simple prior distributions for the orbital inclination and bulk composition. Our model predicts that Proxima b would have a mass $1.3\leq M_{\rm p}\leq 2.3\,M_{\oplus}$ and a radius $R_{\rm p}=1.4^{+0.3}_{-0.2}\,R_{\oplus}$. In our simulations, most Proxima b analogues develop intrinsic dynamos that last for $\geq$4 Gyr (the estimated age of the host star). If alive, the dynamo of Proxima b have a dipole moment ${\cal M}_{\rm dip}>0.32^{\times 2.3}_{\div 2.9}{\cal M}_{\rm dip,\oplus}$. These results are not restricted to Proxima b but they also apply to earth-like planets having similar observed properties.

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 $[\rm \alpha/Fe]$ -- $[\rm 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-$\rm[\alpha/Fe]$ star formation phase in the inner region ($R\lesssim 5$ kpc) induced by gas-rich mergers, followed by more quiescent low-$\rm[\alpha/Fe]$ star formation; and b) an early phase of high-$\rm[\alpha/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-$\rm[\alpha/Fe]$ sequence stars. We discuss the implication that both processes occurred in the Milky Way.

We 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 ($\gtrsim 50\%$) comes from single massive satellites ($> 6 \times 10^{10}~M_{\odot}$). These satellites are accreted with a wide range of infall times and inclination angles (up to $85^{\circ}$). 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.

We present a new method to find voids in cosmological simulations based on the tidal and the velocity shear tensors definitions of the cosmic web. We use the fractional anisotropy (FA) computed from the eigenvalues of each web scheme as a void tracer. We identify voids using a watershed transform based on the local minima of the FA field without making any assumption on the shape or structure of the voids. We test the method on the Bolshoi simulation and report on the abundance and radial averaged profiles for the density, velocity and fractional anisotropy. We find that voids in the velocity shear web are smaller than voids in the tidal web, with a particular overabundance of very small voids in the inner region of filaments/sheets. We classify voids as subcompensated/overcompansated depending on the absence/presence of an overdense matter ridge in their density profile, finding that close to $65\%$ and $35\%$ of the total population are classified into each category, respectively. Finally, we find evidence for the existence of universal profiles from the radially averaged profiles for density, velocity and fractional anisotropy. This requires that the radial coordinate is normalized to the effective radius of each void. Put together, all these results show that the FA is a reliable tracer for voids, which can be used in complementarity to other existing methods and tracers.

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 Rp, larger than previously estimated. Unlocked planets with periods of rotation ~1 day are protected by magnetospheres extending between 3 and 8 Rp. 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.

Recent observations constrained the tangential velocity of M31 with respect to the Milky Way (MW) to be v_tan<34.4 km/s and the radial velocity to be in the range v_rad=-109+/- 4.4 km/s (van der Marel et al. 2012). In this study we use a large volume high resolution N-body cosmological simulation (Bolshoi) together with three constrained simulations to statistically study this kinematics in the context of the LCDM. The comparison of the ensembles of simulated pairs with the observed LG at the 1-sigma level in the uncertainties has been done with respect to the radial and tangential velocities, the reduced orbital energy (e_tot), angular momentum (l_orb) and the dimensionless spin parameter, lambda. Our main results are: (i) the preferred radial and tangential velocities for pairs in LCDM are v_rad=-80+/-20 km/s, v_tan=50+/-10 km/s, (ii) pairs around that region are 3 to 13 times more common than pairs within the observational values, (iii) 15%to 24% of LG-like pairs in LCDM have energy and angular momentum consistent with observations while (iv) 9% to 13% of pairs in the same sample show similar values in the inferred dimensionless spin parameter. It follows that within current observational uncertainties the quasi-conserved quantities that characterize the orbit of the LG, i.e. e_tot, r_orb and lambda, do not challenge the standard LCDM model, but the model is in tension with regard to the actual values of the radial and tangential velocities. This might hint to a problem of the LCDM model to reproduce the observed LG.