We study the late-time evolution of the central regions of two Milky Way-like simulations of galaxies formed in a cosmological context, one hosting a fast bar and the other a slow one. We find that bar length, R_b, measurements fluctuate on a dynamical timescale by up to 100%, depending on the spiral structure strength and measurement threshold. The bar amplitude oscillates by about 15%, correlating with R_b. The Tremaine-Weinberg-method estimates of the bars' instantaneous pattern speeds show variations around the mean of up to ~20%, typically anti-correlating with the bar length and strength. Through power spectrum analyses, we establish that these bar pulsations, with a period in the range ~60-200 Myr, result from its interaction with multiple spiral modes, which are coupled with the bar. Because of the presence of odd spiral modes, the two bar halves typically do not connect at exactly the same time to a spiral arm, and their individual lengths can be significantly offset. We estimated that in about 50% of bar measurements in Milky Way-mass external galaxies, the bar lengths of SBab type galaxies are overestimated by ~15% and those of SBbc types by ~55%. Consequently, bars longer than their corotation radius reported in the literature, dubbed "ultra-fast bars", may simply correspond to the largest biases. Given that the Scutum-Centaurus arm is likely connected to the near half of the Milky Way bar, recent direct measurements may be overestimating its length by 1-1.5 kpc, while its present pattern speed may be 5-10 km/s/kpc smaller than its time-averaged value.

The mass of a star is the most fundamental parameter for its structure, evolution, and final fate. It is particularly important for any kind of stellar archaeology and characterization of exoplanets. There exists a variety of methods in astronomy to estimate or determine it. In this review we present a significant number of such methods, beginning with the most direct and model-independent approach using detached eclipsing binaries. We then move to more indirect and model-dependent methods, such as the quite commonly used isochrone or stellar track fitting. The arrival of quantitative asteroseismology has opened a completely new approach to determine stellar masses and to complement and improve the accuracy of other methods. We include methods for different evolutionary stages, from the pre-main sequence to evolved (super)giants and final remnants. For all methods uncertainties and restrictions will be discussed. We provide lists of altogether more than 200 benchmark stars with relative mass accuracies between $[0.3,2]\%$ for the covered mass range of $M\in [0.1,16]\,M_\odot$, $75\%$ of which are stars burning hydrogen in their core and the other $25\%$ covering all other evolved stages. We close with a recommendation how to combine various methods to arrive at a "mass-ladder" for stars.

Stellar populations in barred galaxies save an imprint of the influence of the bar on the host galaxy's evolution. We present a detailed analysis of star formation histories (SFHs) and chemical enrichment of stellar populations in nine nearby barred galaxies from the TIMER project. We use integral field observations with the MUSE instrument to derive unprecedented spatially resolved maps of stellar ages, metallicities, [Mg/Fe] abundances and SFHs, as well as H$\alpha$ as a tracer of ongoing star formation. We find a characteristic V-shaped signature in the SFH perpendicular to the bar major axis which supports the scenario where intermediate age stars ($\sim 2$-$6\ \mathrm{Gyr}$) are trapped on more elongated orbits shaping a thinner part of the bar, while older stars ($> 8\ \mathrm{Gyr}$) are trapped on less elongated orbits shaping a rounder and thicker part of the bar. We compare our data to state-of-the-art cosmological magneto-hydrodynamical simulations of barred galaxies and show that such V-shaped SFHs arise naturally due to the dynamical influence of the bar on stellar populations with different ages and kinematic properties. Additionally, we find an excess of very young stars ($< 2\ \mathrm{Gyr}$) on the edges of the bars, predominantly on the leading side, confirming typical star formation patterns in bars. Furthermore, mass-weighted age and metallicity gradients are slightly shallower along the bar than in the disc likely due to orbital mixing in the bar. Finally, we find that bars are mostly more metal-rich and less [Mg/Fe]-enhanced than the surrounding discs. We interpret this as a signature that the bar quenches star formation in the inner region of discs, usually referred to as star formation deserts. We discuss these results and their implications on two different scenarios of bar formation and evolution.

Stars in globular clusters (GCs) lose a non negligible amount of mass during their post-main sequence evolution. This material is then expected to build up a substantial intra-cluster medium (ICM) within the GC. However, the observed gas content in GCs is a couple of orders of magnitude below these expectations. Here we follow the evolution of this stellar wind material through hydrodynamical simulations to attempt to reconcile theoretical predictions with observations. We test different mechanisms proposed in the literature to clear out the gas such as ram-pressure stripping by the motion of the GC in the Galactic halo medium and ionisation by UV sources. We use the code ramses to run 3D hydrodynamical simulations to study for the first time the ICM evolution within discretised multi-mass GC models including stellar winds and full radiative transfer. We find that the inclusion of both ram-pressure and ionisation is mandatory to explain why only a very low amount of ionised gas is observed in the core of GCs. The same mechanisms operating in ancient GCs that clear the gas could also be efficient at younger ages, meaning that young GCs would not be able to retain gas and form multiple generations of stars as assumed in many models to explain "multiple populations". However, this rapid clearing of gas is consistent with observations of young massive clusters.

The {\alpha}-element abundances of the globular cluster (GC) and field star populations of galaxies encode information about the formation of each of these components. We use the E-MOSAICS cosmological simulations of ~L* galaxies and their GCs to investigate the [{\alpha}/Fe]-[Fe/H] distribution of field stars and GCs in 25 Milky Way-mass galaxies. The [{\alpha}/Fe]-[Fe/H] distribution go GCs largely follows that of the field stars and can also therefore be used as tracers of the [{\alpha}/Fe]-[Fe/H] evolution of the galaxy. Due to the difference in their star formation histories, GCs associated with stellar streams (i.e. which have recently been accreted) have systematically lower [{\alpha}/Fe] at fixed [Fe/H]. Therefore, if a GC is observed to have low [{\alpha}/Fe] for its [Fe/H] there is an increased probability that this GC was accreted recently alongside a dwarf galaxy. There is a wide range of shapes for the field star [{\alpha}/Fe]-[Fe/H] distribution, with a notable subset of galaxies exhibiting bimodal distributions, in which the high [{\alpha}/Fe] sequence is mostly comprised of stars in the bulge, a high fraction of which are from disrupted GCs. We calculate the contribution of disrupted GCs to the bulge component of the 25 simulated galaxies and find values between 0.3-14 per cent, where this fraction correlates with the galaxy's formation time. The upper range of these fractions is compatible with observationally-inferred measurements for the Milky Way, suggesting that in this respect the Milky Way is not typical of L* galaxies, having experienced a phase of unusually rapid growth at early times.

Several deep observations such as those carried out at the Canada-France-Hawaii Telescope (CFHT) have revealed prominent Low Surface Brightness (LSB) fine structures that change the apparent morphology of galaxies. Previous photometry surveys have developed observational techniques which exploit the diffuse light detected in the external regions of galaxies. In these studies the outer perturbations have been identified and classified like tidal tails, stellar streams, and shells. These structures are tracers of interacting and merging events and they keep a memory of the mass assembly of galaxies. Cosmological numerical simulations are needed to estimate their visibility time-scale (among other properties) in order to reconstruct the past merger history of galaxies. In the present work, we analyze a hydrodynamical cosmological simulation to build a comprehensive interpretation of the properties of fine structures. We make a census of several types of LSB fine structures by visual inspection of individual snapshots at various time. We reconstruct the evolution of the number of fine structures detected around an early-type galaxy and we compare with the merger history of the galaxy. We find that most of fine structures are associated with major and intermediate mass merger events. Their survival time scale ranges between 0.7 and 4 Gyr. Shells and streams remain visible for a longer time than tidal tails. These estimates of survival times provide clues to interpret the shape and frequency of fine structures observed in deep images in terms of mass assembly. We find that the detectability of stellar streams is the most sensitive to the surface brightness limit. We see 2-3 times more streams with a surface brightness cut of 33 mag arcsec$^{-2}$ than with 29 mag arcsec$^{-2}$. The detection of shells display a strong dependence on the projection angle.

Bars strongly influence the distribution of gas and stars within the central regions of their host galaxies. This is particularly pronounced in the star formation desert (SFD) which is defined as two symmetrical regions either side of the bar that show a deficit in young stars. Previous studies proposed that, if star formation is truncated because of the influence of the bar, then the age distribution of stars within the SFD could be used to determine the epoch of bar formation. To test this, we study the properties of SFDs in 6 galaxies from zoom-in cosmological re-simulations. Age maps reveal old regions on both sides of the bars, with a lack of stars younger than 10 Myr, confirming the SFD phenomenon. Local star formation is truncated in the SFDs because after the bar forms, gas in these regions is removed on 1 Gyr timescales. However, the overall age distribution of stars in the SFD does not show a sharp truncation after bar formation but rather a gradual downturn in comparison to that of the bar. This more subtle signature may still give information on bar formation epochs in observed galaxies, but the interpretation will be more difficult than originally hoped. The gradual drop in the SFD age distribution, instead of a truncation, is due to radial migration of stars born in the disk. The SFD is thus one of the only regions where an uncontaminated sample of stars only affected by radial migration can be studied.

We use Gaia DR2 astrometric and line-of-sight velocity information combined with two sets of distances obtained with a Bayesian inference method to study the 3D velocity distribution in the Milky Way disc. We search for variations in all Galactocentric cylindrical velocity components ($V_{\phi}$, $V_R$ and $V_z$) with Galactic radius, azimuth, and distance from the disc mid-plane. We confirm recent work showing that bulk vertical motions in the $R\text{-}z$ plane are consistent with a combination of breathing and bending modes. In the $x\text{-}y$ plane, we show that, although the amplitudes change, the structure produced by these modes is mostly invariant as a function of distance from the plane. Comparing to two different Galactic disc models, we demonstrate that the observed patterns can drastically change in short time intervals, showing the complexity of understanding the origin of vertical perturbations. A strong radial $V_R$ gradient was identified in the inner disc, transitioning smoothly from $16$ km s$^{-1}$ kpc$^{-1}$ at an azimuth of $30^\circ<\phi<45^\circ$ ahead of the Sun-Galactic centre line, to $-16$ km s$^{-1}$ kpc$^{-1}$ at an azimuth of $-45^\circ<\phi<-30^\circ$ lagging the solar azimuth. We use a simulation with no significant recent mergers to show that exactly the opposite trend is expected from a barred potential, but overestimated distances can flip this trend to match the data. Alternatively, using an $N$-body simulation of the Sagittarius dwarf-Milky Way interaction, we demonstrate that a major recent perturbation is necessary to reproduce the observations. Such an impact may have strongly perturbed the existing bar or even triggered its formation in the last $1\text{-}2$ Gyr.

Stellar feedback plays a significant role in modulating star formation, redistributing metals, and shaping the baryonic and dark structure of galaxies -- however, the efficiency of its energy deposition to the interstellar medium is challenging to constrain observationally. Here we leverage HST and ALMA imaging of a molecular gas and dust shell ($M_{H2} \sim 2\times 10^{5} ~{\rm M}_{\odot}$) in an outflow from the nuclear star forming ring of the galaxy NGC 3351, to serve as a boundary condition for a dynamical and energetic analysis of the outflowing ionised gas seen in our MUSE TIMER survey. We use \texttt{STARBURST99} models and prescriptions for feedback from simulations to demonstrate that the observed star formation energetics can reproduce the ionised and molecular gas dynamics -- provided a dominant component of the momentum injection comes from direct photon pressure from young stars, on top of supernovae, photoionisation heating and stellar winds. The mechanical energy budget from these sources is comparable to low luminosity AGN, suggesting that stellar feedback can be a relevant driver of bulk gas motions in galaxy centres - although here $\lesssim 10^{-3}$ of the ionized gas mass is escaping the galaxy. We test several scenarios for the survival/formation of the cold gas in the outflow, including in-situ condensation and cooling. Interestingly, the geometry of the molecular gas shell, observed magnetic field strengths and emission line diagnostics are consistent with a scenario where magnetic field lines aided survival of the dusty ISM as it was initially launched (with mass loading factor $\lesssim 1$) from the ring by stellar feedback. This system's unique feedback driven morphology can hopefully serve as a useful litmus test for feedback prescriptions in magnetohydrodynamical galaxy simulations.

Galactic disc chemical evolution models generally ignore azimuthal surface density variation that can introduce chemical abundance azimuthal gradients. Recent observations, however, have revealed chemical abundance changes with azimuth in the gas and stellar components of both the Milky Way and external galaxies. To quantify the effects of spiral arm density fluctuations on the azimuthal variations of the oxygen and iron abundances in disc galaxies. We develop a new 2D galactic disc chemical evolution model, capable of following not just radial but also azimuthal inhomogeneities. The density fluctuations resulting from a Milky Way-like N-body disc formation simulation produce azimuthal variations in the oxygen abundance gradients of the order of 0.1 dex. Moreover, in agreement with the most recent observations in external galaxies, the azimuthal variations are more evident in the outer galactic regions. Using a simple analytical model, we show that the largest fluctuations with azimuth result near the spiral structure corotation resonance, where the relative speed between spiral and gaseous disc is the slowest. In conclusion we provided a new 2D chemical evolution model capable of following azimuthal density variations. Density fluctuations extracted from a Milky Way-like dynamical model lead to a scatter in the azimuthal variations of the oxygen abundance gradient in agreement with observations in external galaxies. We interpret the presence of azimuthal scatter at all radii by the presence of multiple spiral modes moving at different pattern speeds, as found in both observations and numerical simulations.

The kinematics of the Milky Way disc as a function of age are well measured at the solar radius, but have not been studied over a wider range of Galactocentric radii. Here, we measure the kinematics of mono-age, mono-$\mathrm{[Fe/H]}$ populations in the low and high $\mathrm{[\alpha/Fe]}$ discs between $4 \lesssim R \lesssim 13$ kpc and $|z| \lesssim 2$ kpc using 65,719 stars in common between APOGEE DR14 and $\it{Gaia}$ DR2 for which we estimate ages using a Bayesian neural network model trained on asteroseismic ages. We determine the vertical and radial velocity dispersions, finding that the low and high $\mathrm{[\alpha/Fe]}$ discs display markedly different age--velocity-dispersion relations (AVRs) and shapes $\sigma_z/\sigma_R$. The high $\mathrm{[\alpha/Fe]}$ disc has roughly flat AVRs and constant $\sigma_z/\sigma_R = 0.64\pm 0.04$, whereas the low $\mathrm{[\alpha/Fe]}$ disc has large variations in this ratio which positively correlate with the mean orbital radius of the population at fixed age. The high $\mathrm{[\alpha/Fe]}$ disc component's flat AVRs and constant $\sigma_z/\sigma_R$ clearly indicates an entirely different heating history. Outer disc populations also have flatter radial AVRs than those in the inner disc, likely due to the waning effect of spiral arms. Our detailed measurements of AVRs and $\sigma_z/\sigma_R$ across the disc indicate that low $\mathrm{[\alpha/Fe]}$, inner disc ($R \lesssim 10\,\mathrm{kpc}$) stellar populations are likely dynamically heated by both giant molecular clouds and spiral arms, while the observed trends for outer disc populations require a significant contribution from another heating mechanism such as satellite perturbations. We also find that outer disc populations have slightly positive mean vertical and radial velocities, likely because they are part of the warped disc.

Simpson's paradox, or Yule-Simpson effect, arises when a trend appears in different subsets of data but disappears or reverses when these subsets are combined. We describe here seven cases of this phenomenon for chemo-kinematical relations believed to constrain the Milky Way disk formation and evolution. We show that interpreting trends in relations, such as the radial and vertical chemical abundance gradients, the age-metallicity relation, and the metallicity-rotational velocity relation (MVR), can lead to conflicting conclusions about the Galaxy past if analyses marginalize over stellar age and/or birth radius. It is demonstrated that the MVR in RAVE giants is consistent with being always strongly negative, when narrow bins of [Mg/Fe] are considered. This is directly related to the negative radial metallicity gradients of stars grouped by common age (mono-age populations) due to the inside out disk formation. The effect of the asymmetric drift can then give rise to a positive MVR trend in high-[alpha/Fe] stars, with a slope dependent on a given survey's selection function and observational uncertainties. We also study the variation of lithium abundance, A(Li), with [Fe/H] of AMBRE:HARPS dwarfs. A strong reversal in the positive A(Li)-[Fe/H] trend of the total sample is found for mono-age populations, flattening for younger groups of stars. Dissecting by birth radius shows strengthening in the positive A(Li)-[Fe/H] trend, shifting to higher [Fe/H] with decreasing birth radius; these observational results suggest new constraints on chemical evolution models. This work highlights the necessity for precise age estimates for large stellar samples covering wide spatial regions.

We used deep MUSE observations to perform a stellar-kinematic and population analysis of FCC 153 and FCC 177, two edge-on S0 galaxies in the Fornax cluster. The geometrical definition of the different structural components of these two galaxies allows us to describe the nature of their thick disks. These are both old, relatively metal poor and [Mg/Fe]-enhanced, and their star formation history (SFH) reveals a minor younger component whose chemical properties suggest its later accretion. Moreover, the outer regions of these geometrically defined thick disks show higher values of metallicity and lower values of [Mg/Fe]. These stars probably formed in the thin-disk region and they were dynamically heated to form the flares present in these two galaxies. We propose different formation scenarios for the three populations of these thick disks: in-situ formation, accretion and disk heating. A clear distinction in age is found between the metal poor and [Mg/Fe]-enhanced thick disks (old, $\sim 12-13$ Gyr), and the metal rich and less [Mg/Fe]-enhanced thin disks (young, $\sim 4-5$ Gyr). These two galaxies show signs of relatively recent star formation in their thin disks and nuclear regions. While the thin disks show more continuous SFHs, the nuclei display a rather bursty SFH. These two galaxies are located outside of the densest region of the Fornax cluster where FCC 170 resides. This other edge-on S0 galaxy was studied by \citet{Pinna2019}. We compare and discuss our results with this previous study. The differences between these three galaxies, at different distances from the cluster center, suggest that the environment can have a strong effect on the galaxy evolutionary path.

We present and discuss the stellar kinematics and populations of the S0 galaxy FCC 170 (NGC 1381) in the Fornax cluster, using deep MUSE data from the Fornax 3D survey. We show the maps of the first four moments of the stellar line-of-sight velocity distribution and of the mass-weighted mean stellar age, metallicity and [Mg/Fe] abundance ratio. The high-quality MUSE stellar kinematic measurements unveil the structure of this massive galaxy: a nuclear disk, a bar seen as a boxy bulge with a clear higher-velocity-dispersion X shape, a fast-rotating and flaring thin disk and a slower rotating thick disk. Whereas their overall old age makes it difficult to discuss differences in the formation epoch between these components, we find a clear-cut distinction between metal-rich and less [Mg/Fe]-enhanced populations in the thin-disk, boxy-bulge and nuclear disk, and more metal-poor and [Mg/Fe]-enhanced stars in the thick disk. Located in the densest region of the Fornax cluster, where signs of tidal stripping have been recently found, the evolution of FCC 170 might have been seriously affected by its environment. We discuss the possibility of its "pre-processing" in a subgroup before falling into the present-day cluster, which would have shaped this galaxy a long time ago. The thick disk displays a composite star formation history, as a significant fraction of younger stars co-exist with the main older thick-disk population. The former sub-population is characterized by even lower-metallicity and higher-[Mg/Fe] values, suggesting that these stars formed later and faster in a less chemically evolved satellite, which was subsequently accreted. Finally, we discuss evidence that metal-rich and less [Mg/Fe]-enhanced stars were brought in the outer parts of the thick disk by the flaring of the thin disk.

One of the fundamental assumptions of chemical evolution models (CEMs) of the Milky Way (MW) and other spirals is that higher gas accretion rates are expected in the past, and in the inner regions of the Galaxy. This leads to the so-called `inside-out disc formation scenario'. Yet, these are probably the most unconstrained inputs of such models. In the present paper, we aim at investigating these main assumptions by studying how gas is accreted in four simulated MW-like galaxies assembled within the $\Lambda$CDM scenario. The galaxies were obtained using two different simulation techniques, cosmological setups and initial conditions. Two of them are MW candidates corresponding to the chemodynamical model of Minchev et al. (2013, 2014) (known as MCM) and the Local Group cosmological simulation of Nuza et al. (2014). We investigate vertical and radial gas accretion on to galaxy discs as a function of cosmic time and disc radius. We find that accretion in the MW-like galaxies seem to happen in two distinct phases, namely: an early, more violent period; followed by a subsequent, slowly declining phase. Our simulations seem to give support to the assumption that the amount of gas incorporated into the MW disc exponentially decreases with time, leading to current net accretion rates of $0.6-1\,$M$_\odot\,$yr$^{-1}$. In particular, accretion timescales on to the simulated thin-disc-like structures are within $\sim5-7\,$Gyr, consistent with expectations from CEMs. Moreover, our simulated MW discs are assembled from the inside-out with gas in the inner disc regions accreted in shorter timescales than in external ones, in qualitative agreement with CEMs of the Galaxy. However, this type of growth is not general to all galaxies and it is intimately linked to their particular merger and gas accretion history.

Stellar haloes encode a fossil record of a galaxy's accretion history, generally in the form of structures of low surface brightness, such as stellar streams. While their low surface brightness makes it challenging to determine their age, metallicity, kinematics and spatial structure, the infalling galaxies also deposit globular clusters (GCs) in the halo, which are bright and therefore easier to observe and characterise. To understand how GCs associated with stellar streams can be used to estimate the stellar mass and the infall time of their parent galaxy, we examine a subset of 15 simulations of galaxies and their star clusters from the E-MOSAICS project. E-MOSAICS is a suite of hydrodynamical simulations incorporating a sub-grid model for GC formation and evolution. We find that more massive accreted galaxies typically contribute younger and more metal rich GCs. This lower age results from a more extended cluster formation history in more massive galaxies. In addition, at fixed stellar mass, galaxies that are accreted later host younger clusters, because they can continue to form GCs without being subjected to environmental influences for longer. This explains the large range of ages observed for clusters associated with the Sagittarius dwarf galaxy in the halo of the Milky Way compared to clusters which are thought to have formed in satellites accreted early in the Milky Way's formation history. Using the ages of the GCs associated with the Sagittarius dwarf, we estimate a virial radius crossing lookback time (infall time) of $9.3 \pm 1.8 Gyr$.

The Time Inference with MUSE in Extragalactic Rings (TIMER) project is a survey with the VLT-MUSE integral-field spectrograph of 24 nearby barred galaxies with prominent central structures (e.g., nuclear rings or inner discs). The main goals of the project are: (i) estimating the cosmic epoch when discs of galaxies settle, leading to the formation of bars; (ii) testing the hypothesis whereby discs in more massive galaxies are assembled first; and (iii) characterising the history of external gas accretion in disc galaxies. We present details on the sample selection, observations, data reduction, and derivation of high-level data products, including stellar kinematics, ages and metallicities. We also derive star formation histories and physical properties and kinematics of ionised gas. We illustrate how this dataset can be used for a plethora of scientific applications, e.g., stellar feedback, outflows, nuclear and primary bars, stellar migration and chemical enrichment, and the gaseous and stellar dynamics of nuclear spiral arms, barlenses, box/peanuts and bulges. Amongst our first results - based on a few selected galaxies -, we show that the dynamics of nuclear rings and inner discs is consistent with the picture in which they are formed by bars, that the central few hundred parsecs in massive disc galaxies tend to show a pronounced peak in stellar metallicity, and that nuclear rings can efficiently prevent star formation in this region. Finally, we present evidence that star-bursting nuclear rings can be fed with low-metallicity gas from low-mass companions.

We present a semi-empirical, largely model-independent approach for estimating Galactic birth radii, r_birth, for Milky Way disk stars. The technique relies on the justifiable assumption that a negative radial metallicity gradient in the interstellar medium (ISM) existed for most of the disk lifetime. Stars are projected back to their birth positions according to the observationally derived age and [Fe/H] with no kinematical information required. Applying our approach to the AMBRE:HARPS and HARPS-GTO local samples, we show that we can constrain the ISM metallicity evolution with Galactic radius and cosmic time, [Fe/H]_ISM(r, t), by requiring a physically meaningful r_birth distribution. We find that the data are consistent with an ISM radial metallicity gradient that flattens with time from ~-0.15 dex/kpc at the beginning of disk formation, to its measured present-day value (-0.07 dex/kpc). We present several chemo-kinematical relations in terms of mono-r_birth populations. One remarkable result is that the kinematically hottest stars would have been born locally or in the outer disk, consistent with thick disk formation from the nested flares of mono-age populations and predictions from cosmological simulations. This phenomenon can be also seen in the observed age-velocity dispersion relation, in that its upper boundary is dominated by stars born at larger radii. We also find that the flatness of the local age-metallicity relation (AMR) is the result of the superposition of the AMRs of mono-r_birth populations, each with a well-defined negative slope. The solar birth radius is estimated to be 7.3+-0.6 kpc, for a current Galactocentric radius of 8 kpc.

In chemodynamical evolution models it is usually assumed that the Milky Way galaxy forms from the inside-out implying that gas inflows onto the disk decrease with galactocentric distance. Similarly, to reproduce differences between chemical abundances of the thick disk and bulge with respect to those of the thin disk, higher accretion fluxes at early times are postulated. By using a suite of Milky Way-like galaxies extracted from cosmological simulations, we investigate the accretion of gas on the simulated stellar disks during their whole evolution. In general, we find that the picture outlined above holds, although the detailed behavior depends on the assembly history of the Galaxy and the complexities inherent to the physics of galaxy formation.

Investigations of the origin and evolution of the Milky Way disk have long relied on chemical and kinematic identification of its components to reconstruct our Galactic past. Difficulties in determining precise stellar ages have restricted most studies to small samples, normally confined to the solar neighbourhood. Here we break this impasse with the help of asteroseismic inference and perform a chronology of the evolution of the disk throughout the age of the Galaxy. We chemically dissect the Milky Way disk population using a sample of red giant stars spanning out to 2~kpc in the solar annulus observed by the {\it Kepler} satellite, with the added dimension of asteroseismic ages. Our results reveal a clear difference in age between the low- and high-$\alpha$ populations, which also show distinct velocity dispersions in the $V$ and $W$ components. We find no tight correlation between age and metallicity nor [$\alpha$/Fe] for the high-$\alpha$ disk stars. Our results indicate that this component formed over a period of more than 2~Gyr with a wide range of [M/H] and [$\alpha$/Fe] independent of time. Our findings show that the kinematic properties of young $\alpha$-rich stars are consistent with the rest of the high-$\alpha$ population and different from the low-$\alpha$ stars of similar age, rendering support to their origin being old stars that went through a mass transfer or stellar merger event, making them appear younger, instead of migration of truly young stars formed close to the Galactic bar.

Early type galaxies (ETG) contain most of the stars present in the local Universe and, above a stellar mass of ~5e10 Msun, vastly outnumber spiral galaxies like the Milky Way. These massive spheroidal galaxies have, in the present day, very little gas or dust, and their stellar populations have been evolving passively for over 10 billion years. The physical mechanisms that led to the termination of star formation in these galaxies and depletion of their interstellar medium remain largely conjectural. In particular, there are currently no direct measurements of the amount of residual gas that might be still present in newly quiescent spheroids at high redshift. Here we show that quiescent ETGs at z~1.8, close to their epoch of quenching, contained 2-3 orders of magnitude more dust at fixed stellar mass than local ETGs. This implies the presence of substantial amounts of gas (5-10%), which was however consumed less efficiently than in more active galaxies, probably due to their spheroidal morphology, and consistently with our simulations. This lower star formation efficiency, and an extended hot gas halo possibly maintained by persistent feedback from an active galactic nucleus (AGN), combine to keep ETGs mostly passive throughout cosmic time.

The stellar velocity ellipsoid (SVE) in galaxies can provide important information on the processes that participate in the dynamical heating of their disc components (e.g. giant molecular clouds, mergers, spiral density waves, bars). Earlier findings suggested a strong relation between the shape of the disc SVE and Hubble type, with later-type galaxies displaying more anisotropic ellipsoids and early-types being more isotropic. In this paper, we revisit the strength of this relation using an exhaustive compilation of observational results from the literature on this issue. We find no clear correlation between the shape of the disc SVE and morphological type, and show that galaxies with the same Hubble type display a wide range of vertical-to-radial velocity dispersion ratios. The points are distributed around a mean value and scatter of $\sigma_z/\sigma_R=0.7\pm 0.2$. With the aid of numerical simulations, we argue that different mechanisms might influence the shape of the SVE in the same manner and that the same process (e.g. mergers) does not have the same impact in all the galaxies. The complexity of the observational picture is confirmed by these simulations, which suggest that the vertical-to-radial axis ratio of the SVE is not a good indicator of the main source of disc heating. Our analysis of those simulations also indicates that the observed shape of the disc SVE may be affected by several processes simultaneously and that the signatures of some of them (e.g. mergers) fade over time.

Galaxy formation entails the hierarchical assembly of mass, along with the condensation of baryons and the ensuing, self-regulating star formation. The stars form a collisionless system whose orbit distribution retains dynamical memory that can constrain a galaxy's formation history. The ordered-rotation dominated orbits with near maximum circularity $\lambda_z \simeq1$ and the random-motion dominated orbits with low circularity $\lambda_z \simeq0$ are called kinematically cold and kinematically hot, respectively. The fraction of stars on `cold' orbits, compared to the fraction of stars on `hot' orbits, speaks directly to the quiescence or violence of the galaxies' formation histories. Here we present such orbit distributions, derived from stellar kinematic maps via orbit-based modelling for a well defined, large sample of 300 nearby galaxies. The sample, drawn from the CALIFA survey, includes the main morphological galaxy types and spans the total stellar mass range from $10^{8.7}$ to $10^{11.9}$ solar masses. Our analysis derives the orbit-circularity distribution as a function of galaxy mass, $p(\lambda_z~|~M_\star)$, and its volume-averaged total distribution, $p(\lambda_z)$. We find that across most of the considered mass range and across morphological types, there are more stars on `warm' orbits defined as $0.25\le \lambda_z \le 0.8$ than on either `cold' or `hot' orbits. This orbit-based "Hubble diagram" provides a benchmark for galaxy formation simulations in a cosmological context.

We analyze age-velocity dispersion relations (AVRs) from kinematics of individual stars in eight Local Group galaxies ranging in mass from Carina ($M_{*} \sim 10^{6}$) to M31 ($M_{*} \sim 10^{11}$). Observationally the $\sigma$ vs. stellar age trends can be interpreted as dynamical heating of the stars by GMCs, bars/spiral arms, or merging subhalos; alternatively the stars could have simply been born out of a more turbulent ISM at high redshift and retain that larger velocity dispersion till present day - consistent with recent IFU studies. To ascertain the dominant mechanism and better understand the impact of instabilities and feedback, we develop models based on observed SFHs of these Local Group galaxies in order to create an evolutionary formalism which describes the ISM velocity dispersion due to a galaxy's evolving gas fraction. These empirical models relax the common assumption that the stars are born from gas which has constant velocity dispersion at all redshifts. Using only the observed SFHs as input, the ISM velocity dispersion and a mid-plane scattering model fits the observed AVRs of low mass galaxies without fine tuning. Higher mass galaxies above $M_{vir} > 10^{11}$ need a larger contribution from latent dynamical heating processes (for example minor mergers), in excess of the ISM model. Using the SFHs we also find that supernovae feedback does not appear to be a dominant driver of the gas velocity dispersion compared to gravitational instabilities - at least for dispersions $\sigma \gtrsim 25$ km/s. Together our results point to stars being born with a velocity dispersion close to that of the gas at the time of their formation, with latent dynamical heating operating with a galaxy mass-dependent efficiency. These semi-empirical relations may help constrain the efficiency of feedback and its impact on the physics of disk settling in galaxy formation simulations.

Schwarzschild orbit-based dynamical models are widely used to uncover the internal dynamics of early-type galaxies and globular clusters. Here we present for the first time the Schwarzschild models of late-type galaxies: an SBb galaxy NGC 4210 and an S0 galaxy NGC 6278 from the CALIFA survey. The mass profiles within $2\,R_e$ are constrained well with $1\sigma$ statistical error of $\sim 10\%$. The luminous and dark mass can be disentangled with uncertainties of $\sim 20\%$ and $\sim 50\%$ respectively. From $R_e$ to $2\,R_e$, the dark matter fraction increases from $14\pm10\%$ to $18\pm10\%$ for NGC 4210 and from $15\pm10\%$ to $30\pm20\%$ for NGC 6278. The velocity anisotropy profiles of both $\sigma_r/\sigma_t$ and $\sigma_z/\sigma_R$ are well constrained. The inferred internal orbital distributions reveal clear substructures. The orbits are naturally separated into three components: a cold component with near circular orbits, a hot component with near radial orbits, and a warm component in between. The photometrically-identified exponential disks are predominantly made up of cold orbits only beyond $\sim 1\,R_e$, while they are constructed mainly with the warm orbits inside. Our dynamical hot components are concentrated in the inner regions, similar to the photometrically-identified bulges. The reliability of the results, especially the orbit distribution, are verified by applying the model to mock data.

Deciphering the assembly history of the Milky Way is a formidable task, which becomes possible only if one can produce high-resolution chrono-chemo-kinematical maps of the Galaxy. Data from large-scale astrometric and spectroscopic surveys will soon provide us with a well-defined view of the current chemo-kinematical structure of the Milky Way, but will only enable a blurred view on the temporal sequence that led to the present-day Galaxy. As demonstrated by the (ongoing) exploitation of data from the pioneering photometric missions CoRoT, Kepler, and K2, asteroseismology provides the way forward: solar-like oscillating giants are excellent evolutionary clocks thanks to the availability of seismic constraints on their mass and to the tight age-initial-mass relation they adhere to. In this paper we identify five key outstanding questions relating to the formation and evolution of the Milky Way that will need precise and accurate ages for large samples of stars to be addressed, and we identify the requirements in terms of number of targets and the precision on the stellar properties that are needed to tackle such questions. By quantifying the asteroseismic yields expected from PLATO for red-giant stars, we demonstrate that these requirements are within the capabilities of the current instrument design, provided that observations are sufficiently long to identify the evolutionary state and allow robust and precise determination of acoustic-mode frequencies. This will allow us to harvest data of sufficient quality to reach a 10% precision in age. This is a fundamental pre-requisite to then reach the more ambitious goal of a similar level of accuracy, which will only be possible if we have to hand a careful appraisal of systematic uncertainties on age deriving from our limited understanding of stellar physics, a goal which conveniently falls within the main aims of PLATO's core science.

The radial profiles of stars in disc galaxies are observed to be either purely exponential (Type-I), truncated (Type-II) or anti-truncated (Type-III) exponentials. Controlled formation simulations of isolated galaxies can reproduce all of these profile types by varying a single parameter, the initial halo spin. In this paper we examine these simulations in more detail in an effort to identify the physical mechanism that leads to the formation of Type-III profiles. The stars in the anti-truncated outskirts of such discs are now on eccentric orbits, but were born on near-circular orbits at much smaller radii. We show that, and explain how, they were driven to the outskirts via non-linear interactions with a strong and long-lived central bar, which greatly boosted their semi-major axis but also their eccentricity. While bars have been known to cause radial heating and outward migration to stellar orbits, we link this effect to the formation of Type-III profiles. This predicts that the anti-truncated parts of galaxies have unusual kinematics for disc-like stellar configurations: high radial velocity dispersions and slow net rotation. Whether such discs exist in nature, can be tested by future observations.

Studying the Milky Way disk structure using stars in narrow bins of [Fe/H] and [alpha/Fe] has recently been proposed as a powerful method to understand the Galactic thick and thin disk formation. It has been assumed so far that these mono-abundance populations (MAPs) are also coeval, or mono-age, populations. Here we study this relationship for a Milky Way chemo-dynamical model and show that equivalence between MAPs and mono-age populations exists only for the high-[alpha/Fe] tail, where the chemical evolution curves of different Galactic radii are far apart. At lower [alpha/Fe]-values a MAP is composed of stars with a range in ages, even for small observational uncertainties and a small MAP bin size. Due to the disk inside-out formation, for these MAPs younger stars are typically located at larger radii, which results in negative radial age gradients that can be as large as 2 Gyr/kpc. Positive radial age gradients can result for MAPs at the lowest [alpha/Fe] and highest [Fe/H] end. Such variations with age prevent the simple interpretation of observations for which accurate ages are not available. Studying the variation with radius of the stellar surface density and scale-height in our model, we find good agreement to recent analyses of the APOGEE red-clump (RC) sample when 1-4 Gyr old stars dominate (as expected for the RC). Our results suggest that the APOGEE data are consistent with a Milky Way model for which mono-age populations flare for all ages. We propose observational tests for the validity of our predictions and argue that using accurate age measurements, such as from asteroseismology, is crucial for putting constraints on the Galactic formation and evolution.

In the Milky Way, the thick disk can be defined using individual stellar abundances, kinematics, or age; or geometrically, as stars high above the mid-plane. In nearby galaxies, where only a geometric definition can be used, thick disks appear to have large radial scale-lengths, and their red colors suggest that they are uniformly old. The Milky Way's geometrically thick disk is also radially extended, but it is far from chemically uniform: alpha-enhanced stars are confined within the inner Galaxy. In simulated galaxies, where old stars are centrally concentrated, geometrically thick disks are radially extended, too. Younger stellar populations flare in the simulated disks' outer regions, bringing those stars high above the mid-plane. The resulting geometrically thick disks therefore show a radial age gradient, from old in their central regions to younger in their outskirts. Based on our age estimates for a large sample of giant stars in the APOGEE survey, we can now test this scenario for the Milky Way. We find that the geometrically-defined thick disk in the Milky Way has indeed a strong radial age gradient: the median age for red clump stars goes from ~9 Gyr in the inner disk to 5 Gyr in the outer disk. We propose that at least some nearby galaxies could also have thick disks that are not uniformly old, and that geometrically thick disks might be complex structures resulting from different formation mechanisms in their inner and outer parts.

With the advent of the space missions CoRoT and Kepler, it has become feasible to determine precise asteroseismic masses and ages for large samples of red-giant stars. In this paper, we present the CoRoGEE dataset -- obtained from CoRoT lightcurves for 606 red giant stars in two fields of the Galactic disc which have been co-observed for an ancillary project of APOGEE. We have used the Bayesian parameter estimation code PARAM to calculate distances, extinctions, masses, and ages for these stars in a homogeneous analysis, resulting in relative statistical uncertainties of $\sim2\%$ in distance, $\sim4\%$ in radius, $\sim9\%$ in mass and $\sim25\%$ in age. We also assess systematic age uncertainties due to different input physics and mass loss. We discuss the correlation between ages and chemical abundance patterns of field stars over a large radial range of the Milky Way's disc (5 kpc $<R_{\rm Gal}<$ 14 kpc), focussing on the [$\alpha$/Fe]-[Fe/H]-age plane in five radial bins of the Galactic disc. We find an overall agreement with the expectations of chemical-evolution models computed before the present data were available, especially for the outer regions. However, our data also indicate that a significant fraction of stars now observed near and beyond the Solar Neighbourhood migrated from inner regions. Mock CoRoGEE observations of a chemo-dynamical Milky Way disc model show that the number of high-metallicity stars in the outer disc is too high to be accounted for even by the strong radial mixing present in the model. The mock observations also reveal that the age distribution of the [$\alpha$/Fe]-enhanced sequence in the CoRoGEE inner-disc field is much broader than expected from a combination of radial mixing and observational errors. We suggest that a thick disc/bulge component that formed stars for more than 3 Gyr may account for these discrepancies.

Recent improvements in the age dating of stellar populations and single stars allow us to study the ages and abundance of stars and galaxies with unprecedented accuracy. We here compare the relation between age and \alpha-element abundances for stars in the solar neighborhood to that of local, early-type galaxies. We find both relations to be very similar. Both fall into two regimes with a flat slope for ages younger than ~9 Gyr and a steeper slope for ages older than that value. This quantitative similarity seems surprising, given the different types of galaxies and scales involved. For the sample of early-type galaxies we also show that the data are inconsistent with literature delay time distributions of either single or double Gaussian shape. The data are consistent with a power law delay time distribution. We thus confirm that the delay time distribution inferred for the Milky Way from chemical evolution arguments also must apply to massive early-type galaxies. We also offer a tentative explanation for the seeming universality of the age-[\alpha/Fe] relation as the manifestation of averaging of different stellar populations with varying chemical evolution histories.

Lenticular galaxies (S0s) represent the majority of early-type galaxies in the local Universe, but their formation channels are still poorly understood. While galaxy mergers are obvious pathways to suppress star formation and increase bulge sizes, the marked parallelism between spiral and lenticular galaxies (e.g. photometric bulge-disc coupling) seemed to rule out a potential merger origin. Here, we summarise our recent work in which we have shown, through N-body numerical simulations, that disc-dominated lenticulars can emerge from major mergers of spiral galaxies, in good agreement with observational photometric scaling relations. Moreover, we show that mergers simultaneously increase the light concentration and reduce the angular momentum relative to their spiral progenitors. This explains the mismatch in angular momentum and concentration between spirals and lenticulars recently revealed by CALIFA observations, which is hard to reconcile with simple fading mechanisms (e.g. ram-pressure stripping).

The main goal of Galactic Archaeology is understanding the formation and evolution of the basic Galactic components. This requires sophisticated chemo-dynamical modeling, where disc asymmetries (e.g., perturbations from the central bar, spirals arms, and infalling satellites) and non-equilibrium processes are taken into account self-consistently. Here we discuss the current status of Galactic chemo-dynamical modeling and focus on a recent hybrid technique, which helps circumvent traditional problems with chemical enrichment and star formation encountered in fully self-consistent cosmological simulations. We show that this model can account for a number of chemo-kinematic relations in the Milky Way. In addition, we demonstrate that (1) our model matches well the observed age-[alpha/Fe] and age-[Fe/H] relations and (2) that the scatter in the age-[Fe/H] relation cannot be simply explained by blurring (stars on apo- and pericenters visiting the solar vicinity) but significant radial migration (stars born elsewhere but ending up at the solar vicinity today because of a change in guiding radius) is needed. We emphasize the importance of accurate stellar ages, such as those obtained through asteroseismology by the CoRoT and Kepler missions, for breaking the degeneracy among different Galactic evolution scenarios.

The mass of a star is arguably its most fundamental parameter. For red giant stars, tracers luminous enough to be observed across the Galaxy, mass implies a stellar evolution age. It has proven to be extremely difficult to infer ages and masses directly from red giant spectra using existing methods. From the KEPLER and APOGEE surveys, samples of several thousand stars exist with high-quality spectra and asteroseismic masses. Here we show that from these data we can build a data-driven spectral model using The Cannon, which can determine stellar masses to $\sim$ 0.07 dex from APOGEE DR12 spectra of red giants; these imply age estimates accurate to $\sim$ 0.2 dex (40 percent). We show that The Cannon constrains these ages foremost from spectral regions with CN absorption lines, elements whose surface abundances reflect mass-dependent dredge-up. We deliver an unprecedented catalog of 80,000 giants (including 20,000 red-clump stars) with mass and age estimates, spanning the entire disk (from the Galactic center to R $\sim$ 20 kpc). We show that the age information in the spectra is not simply a corollary of the birth-material abundances [Fe/H] and [$\alpha$/Fe], and that even within a mono-abundance population of stars, there are age variations that vary sensibly with Galactic position. Such stellar age constraints across the Milky Way open up new avenues in Galactic archeology.

We show that the masses of red giant stars can be well predicted from their photospheric carbon and nitrogen abundances, in conjunction with their spectroscopic stellar labels log g, Teff, and [Fe/H]. This is qualitatively expected from mass-dependent post main sequence evolution. We here establish an empirical relation between these quantities by drawing on 1,475 red giants with asteroseismic mass estimates from Kepler that also have spectroscopic labels from APOGEE DR12. We assess the accuracy of our model, and find that it predicts stellar masses with fractional r.m.s. errors of about 14% (typically 0.2 Msun). From these masses, we derive ages with r.m.s errors of 40%. This empirical model allows us for the first time to make age determinations (in the range 1-13 Gyr) for vast numbers of giant stars across the Galaxy. We apply our model to 52,000 stars in APOGEE DR12, for which no direct mass and age information was previously available. We find that these estimates highlight the vertical age structure of the Milky Way disk, and that the relation of age with [alpha/M] and metallicity is broadly consistent with established expectations based on detailed studies of the solar neighbourhood.

A number of simulators have argued that major mergers can sometimes preserve discs (e.g. Springel & Hernquist 2005), but the possibility that they could explain the emergence of lenticular galaxies (S0s) has been generally neglected. In fact, observations of S0s reveal a strong structural coupling between their bulges and discs, which seems difficult to reconcile with the idea that they come from major mergers. However, in Querejeta et al. (2015a) we have used N-body simulations of binary mergers to show that, under favourable conditions, discs are first destroyed but soon regrow out of the leftover debris, matching observational photometric scaling relations (e.g. Laurikainen et al. 2010). Additionally, in Querejeta et al. (2015b) we have shown how the merger scenario agrees with the recent discovery that S0s and most spirals are not compatible in an angular momentum--concentration plane. This important result from CALIFA constitutes a serious objection to the idea that spirals transform into S0s mainly by fading (e.g. via ram-pressure stripping, as that would not explain the observed simultaneous change in $\lambda_\mathrm{Re}$ and concentration), but our simulations of major mergers do explain that mismatch. From such a 3D comparison we conclude that mergers must be a relevant process in the build-up of the current population of S0s.

We have undertaken an ambitious program to visually classify all galaxies in the five CANDELS fields down to H<24.5 involving the dedicated efforts of 65 individual classifiers. Once completed, we expect to have detailed morphological classifications for over 50,000 galaxies up to z<4 over all the fields. Here, we present our detailed visual classification scheme, which was designed to cover a wide range of CANDELS science goals. This scheme includes the basic Hubble sequence types, but also includes a detailed look at mergers and interactions, the clumpiness of galaxies, $k$-corrections, and a variety of other structural properties. In this paper, we focus on the first field to be completed -- GOODS-S. The wide area coverage spanning the full field includes 7634 galaxies that have been classified by at least three different people. In the deep area of the field, 2534 galaxies have been classified by at least five different people at three different depths. With this paper, we release to the public all of the visual classifications in GOODS-S along with the GUI that we developed to classify galaxies. We find that the level of agreement among classifiers is good and depends on both the galaxy magnitude and the galaxy type, with disks showing the highest level of agreement and irregulars the lowest. A comparison of our classifications with the Sersic index and rest-frame colors shows a clear separation between disk and spheroid populations. Finally, we explore morphological k-corrections between the V-band and H-band observations and find that a small fraction (84 galaxies in total) are classified as being very different between these two bands. These galaxies typically have very clumpy and extended morphology or are very faint in the V-band.

We derive age constraints for 1639 red giants in the APOKASC sample for which seismic parameters from Kepler, as well as effective temperatures, metallicities and [alpha/Fe] values from APOGEE DR12 are available. We investigate the relation between age and chemical abundances for these stars, using a simple and robust approach to obtain ages. We first derive stellar masses using standard seismic scaling relations, then determine the maximum possible age for each star as function of its mass and metallicity, independently of its evolutionary stage. While the overall trend between maximum age and chemical abundances is a declining fraction of young stars with increasing [alpha/Fe], at least 14 out of 241 stars with [alpha/Fe]>0.13 are younger than 6 Gyr. Five stars with [alpha/Fe]>0.2 have ages below 4 Gyr. We examine the effect of modifications in the standard seismic scaling relations, as well as the effect of very low helium fractions, but these changes are not enough to make these stars as old as usually expected for alpha-rich stars (i.e., ages greater than 8-9 Gyr). Such unusual alpha-rich young stars have also been detected by other surveys, but defy simple explanations in a galaxy evolution context.

The CALIFA team has recently found that the stellar angular momentum and concentration of late-type spiral galaxies are incompatible with those of lenticular galaxies (S0s), concluding that fading alone cannot satisfactorily explain the evolution from spirals into S0s. Here we explore whether major mergers can provide an alternative way to transform spirals into S0s by analysing the spiral-spiral major mergers from the GalMer database that lead to realistic, relaxed S0-like galaxies. We find that the change in stellar angular momentum and concentration can explain the differences in the $\lambda_\mathrm{Re}$--$R_{90}/R_{50}$ plane found by the CALIFA team. Major mergers thus offer a feasible explanation for the transformation of spirals into S0s.

Gas and dust in star-forming galaxies at the peak epoch of galaxy assembly are presently the topic of intense study, but little is known about the interstellar medium (ISM) of distant, passively evolving galaxies. We report on a deep 3 mm-band search with IRAM/PdBI for molecular (H$_2$) gas in a massive ($M_{\star}{\sim}6{\times}10^{11}M_{\odot}$) elliptical galaxy at z=1.4277, the first observation of this kind ever attempted. We place a 3$\sigma$ upper limit of 0.32 Jy km/s on the flux of the CO($J$=2$\rightarrow$1) line or $L'_{\rm CO}$$<$8.8$\times$10$^{9}$ K km/s pc$^2$, assuming a disk-like CO-morphology and a circular velocity scaling with the stellar velocity dispersion as in local early-type galaxies (ETGs). This translates to an H$_2$ mass of $<$3.9$\times$10$^{10}$($\alpha_{\rm CO}$/4.4)$M_{\odot}$ or a gas fraction of $\lesssim$6% assuming a Salpeter initial mass function (IMF) and an ISM dominated by H$_2$, as observed in many local, high-mass ellipticals. This low value approaches that of local ETGs, suggesting that the low star formation activity in massive, high-z passive galaxies reflects a true dearth of gas and a lesser role for inhibitive mechanisms like morphological quenching.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

Recent spectroscopic observations in the Milky Way suggest that the chemically defined thick disk (stars with high [alpha/Fe] ratios and thus old) has a significantly smaller scale-length than the thin disk. This is in apparent contradiction with observations of external edge-on galaxies, where the thin and thick components have comparable scale-lengths. Moreover, while observed disks do not flare (scale-height does not increase with radius), numerical simulations suggest that disk flaring is unavoidable, resulting from both environmental effects and secular evolution. Here we address these problems by studying two different suites of simulated galactic disks formed in the cosmological context. We show that the scale-heights of coeval populations always increase with radius. However, the total population can be decomposed morphologically into thin and thick disks, which do not flare. We relate this to the disk inside-out formation, where younger populations have increasingly larger scale-lengths and flare at progressively larger radii. In this new picture, thick disks are composed of the imbedded flares of mono-age stellar populations. Assuming that disks form inside out, we predict that morphologically defined thick disks must show a decrease in age (or [alpha/Fe] ratios) with radius and that coeval populations should always flare. This also explains the observed inversion in the metallicity and [alpha/Fe] gradients for stars away from the disk midplane in the Milky Way. The results of this work are directly linked to, and can be seen as evidence of, inside-out disk growth.

We report the discovery of a group of apparently young CoRoT red-giant stars exhibiting enhanced [alpha/Fe] abundance ratios (as determined from APOGEE spectra) with respect to Solar values. Their existence is not explained by standard chemical evolution models of the Milky Way, and shows that the chemical-enrichment history of the Galactic disc is more complex. We find similar stars in previously published samples for which isochrone-ages could be robustly obtained, although in smaller relative numbers, which could explain why these stars have not received prior attention. The young [alpha/Fe]-rich stars are much more numerous in the CoRoT-APOGEE (CoRoGEE) inner-field sample than in any other high-resolution sample available at present, as only CoRoGEE can explore the inner-disc regions and provide ages for its field stars. The kinematic properties of the young [$\alpha$/Fe]-rich stars are not clearly thick-disc like, despite their rather large distances from the Galactic mid-plane. Our tentative interpretation of these and previous intriguing observations in the Milky Way is that these stars were formed close to the end of the Galactic bar, near corotation -- a region where gas can be kept inert for longer times, compared to other regions shocked more frequently by the passage of spiral arms. Moreover, that is where the mass return from older inner-disc stellar generations should be maximal (according to an inside-out disc-formation scenario), further diluting the in-situ gas. Other possibilities to explain these observations (e.g., a recent gas-accretion event) are also discussed.

The radial density profiles of stellar galaxy discs can be well approximated as an exponential. Compared to this canonical form, however, the profiles in the majority of disc galaxies show downward or upward breaks at large radii. Currently, there is no coherent explanation in a galaxy formation context of the radial profile per se, along with the two types of profile breaks. Using a set of controlled hydrodynamic simulations of disc galaxy formation, we find a correlation between the host halo's initial angular momentum and the resulting radial profile of the stellar disc: galaxies that live in haloes with a low spin parameter {\lambda} <~ 0.03 show an up-bending break in their disc density profiles, while galaxies in haloes of higher angular momentum show a down-bending break. We find that the case of pure exponential profiles ({\lambda} ~ 0.035) coincides with the peak of the spin parameter distribution from cosmological simulations. Our simulations not only imply an explanation of the observed behaviours, but also suggest that the physical origin of this effect is related to the amount of radial redistribution of stellar mass, which is anti-correlated with {\lambda}.

[Abridge] In the first paper of this series (paper I) we presented a new approach for studying the chemo-odynamical evolution in disk galaxies, focusing on the Milky Way. Here we extend these results to different distances from the Galactic center, looking for variations of observables that can be related to on-going and future spectroscopic surveys. By separating the effects of kinematic heating and radial migration, we show that migration is much more important, even for the oldest and hottest stellar population. The distributions of stellar birth guiding radii and final guiding radii (signifying contamination from migration and heating, respectively) widen with increasing distance from the Galactic center. As a result, the slope in the age-metallicity relation flattens significantly at Galactic radii larger than solar. The radial metallicity and [Mg/Fe] gradients in our model show significant variations with height above the plane due to changes in the mixture of stellar ages. An inversion in the radial metallicity gradient is found from negative to weakly positive (at r<10 kpc), and from positive to negative for the [Mg/Fe] gradient, with increasing distance from the disk plane. We relate this to the combined effect of (i) the predominance of young stars close to the disk plane and old stars away from it, (ii) the more concentrated older stellar component, and (iii) the flaring of mono-age disk populations. We also investigate the effect of recycled gas flows on the mean [Fe/H] and find that in the region 4<r<12 kpc the introduced errors are less than 0.05-0.1 dex, related to the fact that inward and outward flows mostly cancel in that radial range. We show that radial migration cannot compete with the inside-out formation of the disk, exposed by the more centrally concentrated older disk populations, and consistent with recent observations.

We study the relationship between the structure and star-formation rate (SFR) of X-ray selected low and moderate luminosity active galactic nuclei (AGNs) in the two Chandra Deep Fields, using Hubble Space Telescope imaging from the Cosmic Assembly Near Infrared Extragalactic Legacy Survey (CANDELS) and deep far-infrared maps from the PEP+GOODS-Herschel survey. We derive detailed distributions of structural parameters and FIR luminosities from carefully constructed control samples of galaxies, which we then compare to those of the AGNs. At z~1, AGNs show slightly diskier light profiles than massive inactive (non-AGN) galaxies, as well as modestly higher levels of gross galaxy disturbance (as measured by visual signatures of interactions and clumpy structure). In contrast, at z~2, AGNs show similar levels of galaxy disturbance as inactive galaxies, but display a red central light enhancement, which may arise due to a more pronounced bulge in AGN hosts or due to extinguished nuclear light. We undertake a number of tests of these alternatives, but our results do not strongly favour one interpretation over the other. The mean SFR and its distribution among AGNs and inactive galaxies are similar at z>1.5. At z<1, however, clear and significant enhancements are seen in the SFRs of AGNs with bulge-dominated light profiles. These trends suggest an evolution in the relation between nuclear activity and host properties with redshift, towards a minor role for mergers and interactions at z>1.5.

Galactic Archaeology, i.e. the use of chemo-dynamical information for stellar samples covering large portions of the Milky Way to infer the dominant processes involved in its formation and evolution, is now a powerful method thanks to the large recently completed and ongoing spectroscopic surveys. It is now important to ask the right questions when analyzing and interpreting the information contained in these rich datasets. To this aim, we have developed a chemodynamical model for the Milky Way that provides quantitative predictions to be compared with the chemo-kinematical properties extracted from the stellar spectra. Three key parameters are needed to make the comparison between data and model predictions useful in order to advance in the field, namely: precise proper-motions, distances and ages. The uncertainties involved in the estimate of ages and distances for field stars are currently the main obstacles in the Galactic Archaeology method. Two important developments might change this situation in the near future: asteroseismology and the now launched Gaia. When combined with the large datasets from surveys like RAVE, SEGUE, LAMOST, Gaia-ESO, APOGEE , HERMES and the future 4MOST we will have the basic ingredients for the reconstruction of the MW history in hands. In the light of these observational advances, the development of detailed chemo-dynamical models tailored to the Milky Way is urgently needed in the field. Here we show the steps we have taken, both in terms of data analysis and modelling. The examples shown here illustrate how powerful can the Galactic Archaeology method become once ages and distances are known with better precision than what is currently feasible.

Using the suite of high-resolution zoom re-simulations of individual haloes by Martig et al., and the large-scale simulation \emph{MassiveBlack-II}, we examine the differences in measured galaxy properties from techniques with various aperture definitions of where galaxies end. We perform techniques popular in the literature and present a new technique of our own, where the aperture radius is based on the baryonic mass profiles of simulated (sub)haloes. For the average Milky-Way-mass system, we find the two most popular techniques in the literature return differences of order 30 per cent for stellar mass, a factor of 3 for gas mass, 40 per cent for star formation rate, and factors of several for gas accretion and ejection rates. Individual cases can show variations greater than this, with the severity dependent on the concentration of a given system. The average difference in integrated properties for a more general galaxy population are not as striking, but are still significant for stellar and gas mass, especially for optical-limit apertures. The large differences that can occur are problematic for comparing results from various publications. We stress the importance of both defining and justifying a technique choice and discourage using popular apertures that use an exact fraction of the virial radius, due to the unignorable variation in galaxy-to-(sub)halo size. Finally, we note that technique choice does not greatly affect simulated galaxies from lying within the scatter of observed scaling relations, but it can alter the derived best-fit slope for the Kennicutt-Schmidt relation.

We determine the intrinsic, 3-dimensional shape distribution of star-forming galaxies at 0<z<2.5, as inferred from their observed projected axis ratios. In the present-day universe star-forming galaxies of all masses 1e9 - 1e11 Msol are predominantly thin, nearly oblate disks, in line with previous studies. We now extend this to higher redshifts, and find that among massive galaxies (M* > 1e10 Msol) disks are the most common geometric shape at all z < 2. Lower-mass galaxies at z>1 possess a broad range of geometric shapes: the fraction of elongated (prolate) galaxies increases toward higher redshifts and lower masses. Galaxies with stellar mass 1e9 Msol (1e10 Msol) are a mix of roughly equal numbers of elongated and disk galaxies at z~1 (z~2). This suggests that galaxies in this mass range do not yet have disks that are sustained over many orbital periods, implying that galaxies with present-day stellar mass comparable to that of the Milky Way typically first formed such sustained stellar disks at redshift z~1.5-2. Combined with constraints on the evolution of the star formation rate density and the distribution of star formation over galaxies with different masses, our findings imply that, averaged over cosmic time, the majority of stars formed in disks.

We present a catalogue and images of visually detected features, such as asymmetries, extensions, warps, shells, tidal tails, polar rings, and obvious signs of mergers or interactions, in the faint outer regions (at and outside of R_25) of nearby galaxies. This catalogue can be used in future quantitative studies that examine galaxy evolution due to internal and external factors. We are able to reliably detect outer region features down to a brightness level of 0.03 MJy/sr per pixel at 3.6 microns in the Spitzer Survey of Stellar Structure in Galaxies (S4G). We also tabulate companion galaxies. We find asymmetries in the outer isophotes in 22+/-1 per cent of the sample. The asymmetry fraction does not correlate with galaxy classification as an interacting galaxy or merger remnant, or with the presence of companions. We also compare the detected features to similar features in galaxies taken from cosmological zoom re-simulations. The simulated images have a higher fraction (33 per cent) of outer disc asymmetries, which may be due to selection effects and an uncertain star formation threshold in the models. The asymmetries may have either an internal (e.g., lopsidedness due to dark halo asymmetry) or external origin.

We study seven simulated disc galaxies, three with a quiescent merger history, and four with mergers in their last 9 Gyr of evolution. We compare their structure at z=0 by decomposing them into "mono-age populations" (MAPs) of stars within 500 Myr age bins. All studied galaxies undergo a phase of merging activity at high redshift, so that stars older than 9 Gyr are found in a centrally concentrated component, while younger stars are mostly found in discs. We find that most MAPs have simple exponential radial and vertical density profiles, with a scale-height that typically increases with age. Because a large range of merger histories can create populations with simple structures, this suggests that the simplicity of the structure of mono-abundance populations observed in the Milky Way by Bovy et al. (2012b,c) is not necessarily a direct indicator of a quiescent history for the Milky Way. Similarly, the anti-correlation between scale-length and scale-height does not necessarily imply a merger-free history. However, mergers produce discontinuities between thin and thick disc components, and jumps in the age-velocity relation. The absence of a structural discontinuity between thin and thick disc observed in the Milky Way would seem to be a good indicator that no merger with a mass ratio larger than 1:15-1:10 occurred in the last 9 Gyr. Mergers at higher redshift might nevertheless be necessary to produce the thickest, hottest components of the Milky Way's disc.

We study the relation between stellar ages and vertical velocity dispersion (the age-velocity relation, or AVR) in a sample of seven simulated disc galaxies. In our simulations, the shape of the AVR for stars younger than 9 Gyr depends strongly on the merger history at low redshift, with even 1:10 - 1:15 mergers being able to create jumps in the AVR (although these jumps might not be detectable if the errors on stellar ages are on the order of 30%). For galaxies with a quiescent history at low redshift, we find that the vertical velocity dispersion rises smoothly for ages up to 8-9 Gyr, following a power law with a slope of ~0.5, similar to what is observed in the solar neighbourhood by the Geneva-Copenhagen Survey. For these galaxies, we show that the slope of the AVR is not imprinted at birth, but is the result of subsequent heating. By contrast, in all our simulations, the oldest stars form a significantly different population, with a high velocity dispersion. These stars are usually born kinematically hot in a turbulent phase of intense mergers at high redshift, and also include some stars accreted from satellites. This maximum in velocity dispersion is strongly decreased when age errors are included, suggesting that observations can easily miss such a jump with the current accuracy of age measurements.

Direct N-body simulations of globular clusters in a realistic Milky Way-like potential are carried out using the code NBODY6 to determine the impact of the host galaxy disk mass and geometry on the survival of star clusters. A relationship between disk mass and star cluster dissolution timescale is derived. These N-body models show that doubling the mass of the disk from 5x10^10 solar masses to 10x10^10 solar masses halves the dissolution time of a satellite star cluster orbiting the host galaxy at 6 kpc from the galactic center. Different geometries in a disk of identical mass can determine either the survival or dissolution of a star cluster orbiting within the inner 6 kpc of the galactic center. Furthermore, disk geometry has measurable effects on the mass loss of star clusters up to 15 kpc from the galactic center. N-body simulations performed with a fine output time step show that at each disk crossing the outer layers of star clusters experience an increase in velocity dispersion of ~5% of the average velocity dispersion in the outer section of star clusters. This leads to an enhancement of mass-loss -- a clearly discernable effect of disk shocking. By running models with different inclinations we determine that star clusters with an orbit perpendicular to the Galactic plane have larger mass loss rates than both clusters evolving in the Galactic plane or in an inclined orbit.

Despite the recent advancements in the field of galaxy formation and evolution, fully self-consistent simulations are still unable to make the detailed predictions necessary for the planned and ongoing large spectroscopic and photometry surveys of the Milky Way disc. These difficulties arise from the very uncertain nature of sub-grid physical energy feedback within models, affecting both star formation rates and chemical enrichment. To avoid these problems, we have introduced a new approach which consists of fusing disc chemical evolution models with compatible numerical simulations. We demonstrate the power of this method by showing that a range of observational results can be explained by our new model. We show that due to radial migration from mergers at high redshift and the central bar at later times, a sizable fraction of old metal-poor, high-[alpha/Fe] stars can reach the solar vicinity. This naturally accounts for a number of recent observations related to both the thin and thick discs, despite the fact that we use thin-disc chemistry only. Within the framework of our model, the MW thick disc has emerged naturally from (i) stars born with high velocity dispersions at high redshift, (ii) stars migrating from the inner disc very early on due to strong merger activity, and (iii) further radial migration driven by the bar and spirals at later times. A significant fraction of old stars with thick-disc characteristics could have been born near the solar radius.

The velocity dispersions of stars near the Sun are known to increase with stellar age, but age can be difficult to determine so a proxy like the abundance of alpha elements (e.g., Mg) with respect to iron, [alpha/Fe], is used. Here we report an unexpected behavior found in the velocity dispersion of a sample of giant stars from the RAdial Velocity Experiment (RAVE) survey with high quality chemical and kinematical information, in that it decreases strongly for stars with [Mg/Fe] > 0.4 dex (i.e., those that formed in the first Gyr of the Galaxy's life). These findings can be explained by perturbations from massive mergers in the early Universe, which have affected more strongly the outer parts of the disc, and the subsequent radial migration of stars with cooler kinematics from the inner disc. Similar reversed trends in velocity dispersion are also found for different metallicity subpopulations. Our results suggest that the Milky Way disc merger history can be recovered by relating the observed chemo-kinematic relations to the properties of past merger events.

[Abridged] In this first paper of this series, we present a new approach for studying the chemo-dynamical evolution in disk galaxies, which consists of fusing disk chemical evolution models with compatible numerical simulations of galactic disks. This method avoids known star formation and chemical enrichment problems encountered in simulations. Here we focus on the Milky Way, by using a detailed thin-disk chemical evolution model and a simulation in the cosmological context, with dynamical properties close to those of our Galaxy. We show that, due to radial migration from mergers at high redshift and the central bar at later times, a sizable fraction of old metal-poor high-[alpha/Fe] stars reaches the solar vicinity. This naturally accounts for a number of recent observations related to both the thin and thick disks, despite the fact that we use thin-disk chemistry only. Although significant radial mixing is present, the slope in the age-metallicity relation is only weakly affected, with a scatter compatible with recent observational work. While we find a smooth density distribution in the [O/Fe]-[Fe/H] plane, we can recover the observed discontinuity by selecting particles according to kinematic criteria used in high-resolution samples to define the thin and thick disks. We outline a new method for estimating the birth place of the Sun and predict that the most likely radius lies in the range 4.4 < r < 7.7 kpc (for a location at r = 8 kpc). A new, unifying model for the Milky Way thick disk is offered, where both mergers and radial migration play a role at different stages of the disk evolution. We show that in the absence of early-on massive mergers the vertical velocity dispersion of the oldest stars is underestimated by a factor of ~2 compared with observations. We can, therefore, argue that the Milky Way thick disk is unlikely to have been formed through a quiescent disk evolution.

We study the global efficiency of star formation in high resolution hydrodynamical simulations of gas discs embedded in isolated early-type and spiral galaxies. Despite using a universal local law to form stars in the simulations, we find that the early-type galaxies are offset from the spirals on the large-scale Kennicutt relation, and form stars 2 to 5 times less efficiently. This offset is in agreement with previous results on morphological quenching: gas discs are more stable against star formation when embedded in early-type galaxies due to the lower disc self-gravity and increased shear. As a result, these gas discs do not fragment into dense clumps and do not reach as high densities as in the spiral galaxies. Even if some molecular gas is present, the fraction of very dense gas (above 10^4 cm-3) is significantly reduced, which explains the overall lower star formation efficiency. We also analyse a sample of local early-type and spiral galaxies, measuring their CO and HI surface densities and their star formation rates as determined by their non-stellar 8um emission. As predicted by the simulations, we find that the early-type galaxies are offset from the Kennicutt relation compared to the spirals, with a twice lower efficiency. Finally, we validate our approach by performing a direct comparison between models and observations. We run a simulation designed to mimic the stellar and gaseous properties of NGC524, a lenticular galaxy, and find a gas disc structure and global star formation rate in good agreement with the observations. Morphological quenching thus seems to be a robust mechanism, and is also consistent with other observations of a reduced star formation efficiency in early-type galaxies in the COLD GASS survey. This lower efficiency of star formation is not enough to explain the formation of the whole Red Sequence, but can contribute to the reddening of some galaxies.

Non-axisymmetric components, such as spirals and central bars, play a major role in shaping galactic discs. An important aspect of the disc secular evolution driven by these perturbers is the radial migration of stars. It has been suggested recently that migration can populate a thick-disc component from inner-disc stars with high vertical energies. Since this has never been demonstrated in simulations, we study in detail the effect of radial migration on the disc velocity dispersion and disc thickness, by separating simulated stars into migrators and non-migrators. We apply this method to three isolated barred Tree-SPH N-body galaxies with strong radial migration. Contrary to expectations, we find that as stellar samples migrate, on the average, their velocity dispersion change (by as much as 50%) in such a way as to approximately match the non-migrating population at the radius at which they arrive. We show that, in fact, migrators suppress heating in parts of the disc. To confirm the validity of our findings, we also apply our technique to three cosmological re-simulations, which use a completely different simulation scheme and, remarkably, find very similar results. We believe the inability of migration to thicken discs is a fundamental property of internal disc evolution, irrespective of the migration mechanism at work. We explain this with the approximate conservation of the (average) vertical and radial actions rather than the energy. This "action mixing" can be used to constrain the migration rate in the Milky Way: estimates of the average vertical action in observations for different populations of stars should reveal flattening with radius for older groups of stars.

We analyze a suite of 33 cosmological simulations of the evolution of Milky Way-mass galaxies in low-density environments. Our sample spans a broad range of Hubble types at z=0, from nearly bulgeless disks to bulge-dominated galaxies. Despite the fact that a large fraction of the bulge is typically in place by z=1, we find no significant correlation between the morphology at z=1 and at z=0. The z=1 progenitors of disk galaxies span a range of morphologies, including smooth disks, unstable disks, interacting galaxies and bulge-dominated systems. By z=0.5, spiral arms and bars are largely in place and the progenitor morphology is correlated with the final morphology. We next focus on late-type galaxies with a bulge-to-total ratio B/T<0.3 at z=0. These show a correlation between B/T at z=0 and the mass ratio of the largest merger at z<2, as well as with the gas accretion rate at z>1. We find that the galaxies with the lowest B/T tend to have a quiet baryon input history, with no major mergers at z<2, and with a low and constant gas accretion rate that keeps a stable angular-momentum direction. More violent merger or gas accretion histories lead to galaxies with more prominent bulges. Most disk galaxies have a bulge Sersic index n<2. The galaxies with the highest bulge Sersic index tend to have histories of intense gas accretion and disk instability rather than active mergers.

We study the evolution of galactic bars and the link with disk and spheroid formation in a sample of zoom-in cosmological simulations. Our simulation sample focuses on galaxies with present-day stellar masses in the 10^10-10^11 Msun range, in field and loose group environments, with a broad variety of mass growth histories. In our models, bars are almost absent from the progenitors of present-day spirals at z>1.5, and they remain rare and generally too weak to be observable down to z~1. After this characteristic epoch, the fractions of observable and strong bars raise rapidly, bars being present in 80% of spiral galaxies and easily observable in two thirds of these at z<0.5. This is quantitatively consistent with the redshift evolution of the observed bar fraction. Our models predict that the decrease in the bar fraction with increasing redshift should continue with a fraction of observable bars <10-15% in disk galaxies at z>1. Our models also predict later bar formation in lower-mass galaxies, in agreement with existing data. We find that the characteristic epoch of bar formation, namely redshift z~0.8-1, corresponds to the epoch at which today's spirals acquire their disk-dominated morphology. At higher redshift, disks tend to be rapidly destroyed by mergers and gravitational instabilities and rarely develop significant bars. The bar formation epoch corresponds to the transition between an early "violent" phase of spiral galaxy formation at z>1 and a late "secular" phase at z<0.8. In the secular phase, the presence of bars substantially contributes to the growth of the bulge, but the bulge mass budget remains statistically dominated by the contribution of mergers, interactions and disk instabilities at high redshift. Early bars at z>1 are often short-lived, while most of the bars formed at z<1 persist down to z=0, late cosmological gas infall being necessary to maintain some of them.

Spiral galaxies have most of their stellar mass in a large rotating disk, and only a modest fraction in a central spheroidal bulge. This poses a major challenge for cosmological models of galaxy formation. Galaxies form at the centre of dark matter halos through a combination of hierarchical merging and gas accretion along cold streams, and should rapidly grow their bulge through mergers and instabilities. Cosmological simulations predict galaxies to have most of their mass in the central bulge, and therefore an angular momentum much below the observed level, except in dwarf galaxies. We propose that the continuous return of fresh gas by stellar populations over cosmic times could solve this issue. A population of stars formed at a given instant typically returns half of its initial mass in the form of gas over 10 billion years, and the process is not dominated by rapid supernovae explosions but by the long-term mass-loss from low- and intermediate-mass stars. Using simulations of galaxy formation, we show that this recycling of gas can strongly affect the structural evolution of massive galaxies, potentially solving the bulge fraction issue: we find that the bulge-to-disk ratio of a massive galaxy can be divided by a factor of 3. The continuous recycling of baryons through star formation and stellar mass loss helps the growth of disks and their survival to interactions and mergers. Instead of forming only early-type, spheroid-dominated galaxies (S0 and ellipticals), the standard cosmological model can then successfully account for massive late-type, disk-dominated spiral galaxies (Sb-Sc).

We point out a natural mechanism for quenching of star formation in early-type galaxies. It automatically links the color of a galaxy with its morphology and does not require gas consumption, removal or termination of gas supply. Given that star formation takes place in gravitationally unstable gas disks, it can be quenched when a disk becomes stable against fragmentation to bound clumps. This can result from the growth of a stellar spheroid, for instance by mergers. We present the concept of morphological quenching (MQ) using standard disk instability analysis, and demonstrate its natural occurrence in a cosmological simulation using an efficient zoom-in technique. We show that the transition from a stellar disk to a spheroid can be sufficient to stabilize the gas disk, quench star formation, and turn an early-type galaxy red and dead while gas accretion continues. The turbulence necessary for disk stability can be stirred up by sheared perturbations within the disk in the absence of bound star-forming clumps. While gas stripping processes are limited to dense groups and clusters, and other quenching mechanisms like AGN feedback, virial shock heating and gravitational heating, are limited to halos more massive than 10^12 Mo, the MQ can explain the appearance of red ellipticals even in less massive halos and in the field. The dense gas disks observed in some of today's red ellipticals may be the relics of this mechanism, whereas red galaxies with quenched gas disks are expected to be more frequent at high redshift.

The formation of thick stellar disks in spiral galaxies is studied. Simulations of gas-rich young galaxies show formation of internal clumps by gravitational instabilities, clump coalescence into a bulge, and disk thickening by strong stellar scattering. The bulge and thick disks of modern galaxies may form this way. Simulations of minor mergers make thick disks too, but there is an important difference. Thick disks made by internal processes have a constant scale height with galactocentric radius, but thick disks made by mergers flare. The difference arises because in the first case, perpendicular forcing and disk-gravity resistance are both proportional to the disk column density, so the resulting scale height is independent of this density. In the case of mergers, perpendicular forcing is independent of the column density and the low density regions get thicker; the resulting flaring is inconsistent with observations. Late-stage gas accretion and thin disk growth are shown to preserve the constant scale heights of thick disks formed by internal evolution. These results reinforce the idea that disk galaxies accrete most of their mass smoothly and acquire their structure by internal processes, in particular through turbulent and clumpy phases at high redshift.

Many evidences exist for a connection between galaxy interactions and induced star formation. However, a large range of responses of galaxies to tidal interactions is found, both in observations and in numerical simulations. We will discuss some recent results obtained analysing a large sample (~ 1000) of simulations of interacting pairs and their agreement with the most recent observational works.

We investigate the intensity enhancement and the duration of starburst episodes, triggered by major galaxy interactions and mergers. To this aim, we analyze two large statistical datasets of numerical simulations. These have been obtained using two independent and different numerical techniques to model baryonic and dark matter evolution, that are extensively compared for the first time. One is a Tree-SPH code, the other one is a grid-based N-body sticky-particles code. We show that, at low redshift, galaxy interactions and mergers in general trigger only moderate star formation enhancements. Strong starbursts where the star formation rate is increased by a factor larger than 5 are rare and found only in about 15% of major galaxy interactions and mergers. Merger-driven starbursts are also rather short-lived, with a typical duration of the activity of a few 10^8 yr. These conclusions are found to be robust, independent from the numerical techniques and star formation models. At higher redshifts where galaxies contain more gas, gas inflow-induced starbursts are neither stronger neither longer than their local counterparts. In turn, the formation of massive gas clumps, results of local Jeans instability that can occur spontaneously in gas-rich disks or be indirectly favored by galaxy interactions, could play a more important role in determining the duration and intensity of star formation episodes.

Large volume cosmological simulations succeed in reproducing the large-scale structure of the Universe. However, they lack resolution and may not take into account all relevant physical processes to test if the detail properties of galaxies can be explained by the CDM paradigm. On the other hand, galaxy-scale simulations could resolve this in a robust way but do not usually include a realistic cosmological context. To study galaxy evolution in cosmological context, we use a new method that consists in coupling cosmological simulations and galactic scale simulations. For this, we record merger and gas accretion histories from cosmological simulations and re-simulate at very high resolution the evolution of baryons and dark matter within the virial radius of a target galaxy. This allows us for example to better take into account gas evolution and associated star formation, to finely study the internal evolution of galaxies and their disks in a realistic cosmological context. We aim at obtaining a statistical view on galaxy evolution from z = 2 to 0, and we present here the first results of the study: we mainly stress the importance of taking into account gas accretion along filaments to understand galaxy evolution.

General relativity as one the pillar of modern cosmology has to be thoroughly tested if we want to achieve an accurate cosmology. We present the results from such a test on cosmological scales using cosmic shear and galaxy clustering measurements. We parametrize potential deviation from general relativity as a modification to the cosmological Poisson equation. We consider two models relevant either for some linearized theory of massive gravity or for the physics of extra-dimensions. We use the latest observations from the CFHTLS-Wide survey and the SDSS survey to set our constraints. We do not find any deviation from general relativity on scales between 0.04 and 10 Mpc. We derive constraints on the graviton mass in a restricted class of model.

Star formation in galaxies is for a part driven by galaxy mergers. At low redshift, star formation activity is low in high-density environments like groups and clusters, and the star formation activity of galaxies increases with their isolation. This star formation-density relation is observed to be reversed at z~1, which is not explained by theoretical models so far. We study the influence of the tidal field of a galaxy group or cluster on the star formation activity of merging galaxies, using N-body simulations including gas dynamics and star formation. We find that the merger-driven star formation is significantly more active in the vicinity of such cosmological structures compared to mergers in the field. The large-scale tidal field can thus enhance the activity of galaxies in dense cosmic structures, and should be particularly efficient at high redshift before quenching processes take effect in the densest regions.