Recent determinations of the radial distributions of mono-metallicity populations (MMPs, i.e., stars in narrow bins in [Fe/H] within wider [$\alpha$/Fe] ranges) by the SDSS-III/APOGEE DR12 survey cast doubts on the classical thin - thick disk dichotomy. The analysis of these observations lead to the non-$[\alpha$/Fe] enhanced populations splitting into MMPs with different surface densities according to their [Fe/H]. By contrast, $[\alpha$/Fe] enhanced (i.e., old) populations show an homogeneous behaviour. We analyze these results in the wider context of disk formation within non-isolated halos embedded in the Cosmic Web, resulting in a two-phase mass assembly. By performing hydrodynamical simulations in the context of the $\rm \Lambda CDM$ model, we have found that the two phases of halo mass assembly (an early, fast phase, followed by a slow one, with low mass assembly rates) are very relevant to determine the radial structure of MMP distributions, while radial mixing has only a secondary role, depending on the coeval dynamical and/or destabilizing events. Indeed, while the frequent dynamical violent events occuring at high redshift remove metallicity gradients, and imply efficient stellar mixing, the relatively quiescent dynamics after the transition keeps [Fe/H] gaseous gradients and prevents newly formed stars to suffer from strong radial mixing. By linking the two-component disk concept with the two-phase halo mass assembly scenario, our results set halo virialization (the event marking the transition from the fast to the slow phases) as the separating event marking periods characterized by different physical conditions under which thick and thin disk stars were born.

We present a detailed analysis of the local evolution of 206 Lagrangian Volumes (LVs) selected at high redshift around galaxy seeds, identified in a large-volume $\Lambda$ cold dark matter ($\Lambda$CDM) hydrodynamical simulation. The LVs have a mass range of $1 - 1500 \times 10^{10} M_\odot$. We follow the dynamical evolution of the density field inside these initially spherical LVs from $z=10$ up to $z_{\rm low} = 0.05$, witnessing highly non-linear, anisotropic mass rearrangements within them, leading to the emergence of the local cosmic web (CW). These mass arrangements have been analysed in terms of the reduced inertia tensor $I_{ij}^r$, focusing on the evolution of the principal axes of inertia and their corresponding eigendirections, and paying particular attention to the times when the evolution of these two structural elements declines. In addition, mass and component effects along this process have also been investigated. We have found that deformations are led by dark matter dynamics and they transform most of the initially spherical LVs into prolate shapes, i.e. filamentary structures. An analysis of the individual freezing-out time distributions for shapes and eigendirections shows that first most of the LVs fix their three axes of symmetry (like a skeleton) early on, while accretion flows towards them still continue. Very remarkably, we have found that more massive LVs fix their skeleton earlier on than less massive ones. We briefly discuss the astrophysical implications our findings could have, including the galaxy mass-morphology relation and the effects on the galaxy-galaxy merger parameter space, among others.

We analyze and compare the bulges of a sample of L* spiral galaxies in hydrodynamical simulations in a cosmological context, using two different codes, P-DEVA and GASOLINE. The codes regulate star formation in very different ways, with P-DEVA simulations inputing low star formation efficiency under the assumption that feedback occurs on subgrid scales, while the GASOLINE simulations have feedback which drives large scale outflows. In all cases, the marked knee-shape in mass aggregation tracks, corresponding to the transition from an early phase of rapid mass assembly to a later slower one, separates the properties of two populations within the simulated bulges. The bulges analyzed show an important early starburst resulting from the collapse-like fast phase of mass assembly, followed by a second phase with lower star formation, driven by a variety of processes such as disk instabilities and/or mergers. Classifying bulge stellar particles identified at z=0 into old and young according to these two phases, we found bulge stellar sub-populations with distinct kinematics, shapes, stellar ages and metal contents. The young components are more oblate, generally smaller, more rotationally supported, with higher metallicity and less alpha-element enhanced than the old ones. These results are consistent with the current observational status of bulges, and provide an explanation for some apparently paradoxical observations, such as bulge rejuvenation and metal-content gradients observed. Our results suggest that bulges of L* galaxies will generically have two bulge populations which can be likened to classical and pseudo-bulges, with differences being in the relative proportions of the two, which may vary due to galaxy mass and specific mass accretion and merger histories.

We present a detailed analysis of the global and fine structure of four middle-mass disc galaxies obtained from simulations in a $\Lambda$CDM scenario. These objects have photometric D/T ratios in good agreement with those observed for late-type spirals, as well as kinematic properties in agreement with the observational Tully-Fisher relation. We identify the different dynamical components at z=0 on the basis of both orbital parameters and the binding energy of stars in the galaxy. In this way, we recognize a slowly rotating centrally concentrated spheroid, and two disc components supported by rotation: a thin disc with stars in nearly circular orbits, and a thick disc with orbital parameters transitional between the thin disc and the spheroid. The spheroidal component is composed mainly by old, metal-poor and {\alpha}-enhanced stars. The distribution of metals in this component shows, however, a clear bimodality with a low-metallicity peak, which could be related to a classical bulge, and a high-metallicity peak, which could be related to a pseudo-bulge. The thin disc appears in our simulations as the youngest and most metal-rich component. The radial distribution of ages and colours in this component are U-shaped: the new stars are forming in the inner regions, and then migrate through secular processes. Finally, we also find a thick disc containing about 16% of the total stellar mass and with properties that are intermediate between those of the thin disc and the spheroid. Its low-metallicity stars are {\alpha}-enhanced when compared to thin disc stars of the same metallicity. The structural parameters (e.g., the scale height) of the simulated thick discs suggest that such a component could result from the combination of different thickening mechanisms that include merger-driven processes, but also long-lived internal perturbations of the thin disc. [Abridged]

The mass assembly and star formation histories of massive galaxies identified at low redshift z in different cosmological hydrodynamical simulations, have been studied through a detailed follow-up backwards in time of their constituent mass elements (sampled by particles) of different types. Then, the configurations they depict at progressively higher zs have been analysed. The analyses show that these histories share common generic patterns, irrespective of particular circumstances. In any case, the results we have found are different depending on the particle type. The most outstanding differences follow. We have found that by z ~ 3.5 - 6, mass elements identified as stellar particles at z=0 exhibit a gaseous cosmic-web-like morphology with scales of ~ 1 physical Mpc, where the densest mass elements have already turned into stars by z ~ 6. These settings are in fact the densest pieces of the cosmic web, where no hot particles show up, and dynamically organized as a hierarchy of flow convergence regions, that is, attraction basins for mass flows. On the other hand, mass elements identified at the diffuse hot coronae surrounding massive galaxies at z = 0, do not display a clear web-like morphology at any z. Diffuse gas is heated when flow convergence regions go through contractive deformations, and most of it keeps hot and with low density along the evolution. To shed light on the physical foundations of the behaviour our analyses show up, as well as on their possible observational implications, these patterns have been confronted with some generic properties of singular flows as described by the adhesion model. We have found that these common patterns simulations show can be interpreted as a consequence of flow properties, that, moreover, could explain different generic observational results on massive galaxies or their samples. We briefly discuss some of them.[Abridged]

Relaxed, massive galactic objects have been identified at redshifts z = 4;5; and 6 in hydrodynamical simulations run in a large cosmological volume. This allowed us to analyze the assembly patterns of the high mass end of the galaxy distribution at these high zs, by focusing on their structural and dynamical properties. Our simulations indicate that massive objects at high redshift already follow certain scaling relations. These relations define virial planes at the halo scale, whereas at the galactic scale they define intrinsic dynamical planes that are, however, tilted relative to the virial plane. Therefore, we predict that massive galaxies must lie on fundamental planes from their formation. We briefly discuss the physical origin of the tilt in terms the physical processes underlying massive galaxy formation at high z, in the context of a two-phase galaxy formation scenario. Specifically, we have found that it lies on the different behavior of the gravitationally heated gas as compared with cold gas previously involved in caustic formation, and the mass dependence of the energy available to heat the gas.

We present SPH cosmological simulations of the formation of three disk galaxies with a detailed treatment of chemical evolution and cooling. The resulting galaxies have properties compatible with observations: relatively high disk-to-total ratios, thin stellar disks and good agreement with the Tully-Fisher and the luminosity-size relations. They present a break in the luminosity profile at 3.0 +- 0.5 disk scale lengths, while showing an exponential mass profile without any apparent breaks, in line with recent observational results. Since the stellar mass profile is exponential, only differences in the stellar populations can be the cause of the luminosity break. Although we find a cutoff for the star formation rate imposed by a density threshold in our star formation model, it does not coincide with the luminosity break and is located at 4.3 +- 0.4 disk scale lengths, with star formation going on between both radii. The color profiles and the age profiles are "U-shaped", with the minimum for both profiles located approximately at the break radius. The SFR to stellar mass ratio increases until the break, explaining the coincidence of the break with the minimum of the age profile. Beyond the break we find a steep decline in the gas density and, consequently, a decline in the SFR and redder colors. We show that most stars (64-78%) in the outer disk originate in the inner disk and afterwards migrate there. Such stellar migrations are likely the main origin of the U-shaped age profile and, therefore, of the luminosity break.

We describe an SPH model for chemical enrichment and radiative cooling in cosmological simulations of structure formation. This model includes: i) the delayed gas restitution from stars by means of a probabilistic approach designed to reduce the statistical noise and, hence, to allow for the study of the inner chemical structure of objects with moderately high numbers of particles; ii) the full dependence of metal production on the detailed chemical composition of stellar particles by using, for the first time in SPH codes, the Qij matrix formalism that relates each nucleosynthetic product to its sources; and iii) the full dependence of radiative cooling on the detailed chemical composition of gas particles, achieved through a fast algorithm using a new metallicity parameter zeta(T) that gives the weight of each element on the total cooling function. The resolution effects and the results obtained from this SPH chemical model have been tested by comparing its predictions in different problems with known theoretical solutions. We also present some preliminary results on the chemical properties of elliptical galaxies found in self-consistent cosmological simulations. Such simulations show that the above zeta-cooling method is important to prevent an overestimation of the metallicity-dependent cooling rate, whereas the Qij formalism is important to prevent a significant underestimation of the [alpha/Fe] ratio in simulated galaxy-like objects.