Search

We use numerical simulations to model Gaia DR3 data with the aim of constraining the Milky Way bar and spiral structure parameters. We show that both the morphology and the velocity field in Milky Way-like galactic disc models are strong functions of time, changing dramatically over a few tens of Myr. This suggests that by finding a good match to the observed radial velocity field, v_R(x,y), we can constrain the bar-spiral orientation. Incorporating uncertainties into our models is necessary to match the data; most importantly, a heliocentric distance uncertainty above 10-15% distorts the bar's shape and v_R quadrupole pattern morphology, and decreases its apparent angle with respect to the Sun-Galactocentric line. An excellent match to the Gaia DR3 v_R(x,y) field is found for a simulation with a bar length R_b\approx3.6 kpc. We argue that the data are consistent with a MW bar as short as ~3 kpc, for moderate strength inner disc spiral structure (A_2/A_0\approx0.25) or, alternatively, with a bar length up to ~5.2 kpc, provided that spiral arms are quite weak (A_2/A_0\approx0.1), and is most likely in the process of disconnecting from a spiral arm. We demonstrate that the bar angle and distance uncertainty can similarly affect the match between our models and the data - a smaller bar angle (20 deg instead of 30 deg) requires smaller distance uncertainty (20% instead of 30%) to explain the observations. Fourier components of the face-on density distribution of our models suggest that the MW does not have strong m=1 and/or m=3 spirals near the solar radius., Comment: 17p, 12 figs. Accepted to MNRAS on Dec 22, 2023

We present GECKOS (Generalising Edge-on galaxies and their Chemical bimodalities, Kinematics, and Outflows out to Solar environments), a new ESO VLT/MUSE large program. The main aim of GECKOS is to reveal the variation in key physical processes of disk formation by connecting Galactic Archaeology with integral field spectroscopic observations of nearby galaxies. Edge-on galaxies are ideal for this task: they allow us to disentangle the assembly history imprinted in thick disks and provide the greatest insights into outflows. The GECKOS sample of 35 nearby edge-on disk galaxies is designed to trace the assembly histories and properties of galaxies across a large range of star formation rates, bulge-to-total ratios, and boxy and non-boxy bulges. GECKOS will deliver spatially resolved measurements of stellar abundances, ages, and kinematics, as well as ionised gas metallicities, ionisation parameters, pressure, and inflow and outflow kinematics; all key parameters for building a complete chemodynamical picture of disk galaxies. With these data, we aim to extend Galactic analysis methods to the wider galaxy population, reaping the benefits of detailed Milky Way studies, while probing the diverse mechanisms of galaxy evolution., Comment: 6 pages, 2 figures, Proceedings of IAUS 377, eds. F. Tabatabaei, B. Barbuy, and Y. Ting

Warps are observed in a large fraction of disc galaxies, and can be due to a large number of different processes. Some of these processes might also cause vertical heating and flaring. Using a sample of galaxies simulated in their cosmological context, we study the connection between warping and disc heating. We analyse the vertical stellar density structure within warped stellar discs, and monitor the evolution of the scale-heights of the mono-age populations and the geometrical thin and thick disc during the warp's lifetime. We also compare the overall thickness and the vertical velocity dispersion in the disc before and after the warp. We find that for warps made of pre-existing stellar particles shifted off-plane, the scale-heights do not change within the disc's warped region: discs tilt rigidly. For warps made of off-plane new stellar material (either born in-situ or accreted), the warped region of the disc is not well described by a double $\mathrm{sech^2}$ density profile. Yet, once the warp is gone, the thin and thick disc structure is recovered, with their scale-heights following the same trends as in the region that was never warped. Finally, we find that the overall thickness and vertical velocity dispersion do not increase during a warp, regardless of the warp's origin. This holds even for warps triggered by interactions with satellites, which cause disc heating but before the warp forms. Our findings suggest that the vertical structure of galaxies does not hold any memory of past warps., Comment: 10 pages, 7 figures, 1 table. Submitted to MNRAS

We introduce BAYES-LOSVD, a novel implementation of the non-parametric extraction of line-of-sight velocity distributions (LOSVDs) in galaxies. We employ bayesian inference to obtain robust LOSVDs and associated uncertainties. Our method relies on principal component analysis to reduce the dimensionality of the base of templates required for the extraction and thus increase the performance of the code. In addition, we implement several options to regularise the output solutions. Our tests, conducted on mock spectra, confirm the ability of our approach to model a wide range of LOSVD shapes, overcoming limitations of the most widely used parametric methods (e.g. Gauss-Hermite expansion). We present examples of LOSVD extractions for real galaxies with known peculiar LOSVD shapes, i.e. NGC4371, IC0719 and NGC4550, using MUSE and SAURON integral-field unit (IFU) data. Our implementation can also handle data from other popular IFU surveys (e.g. ATLAS3D, CALIFA, MaNGA, SAMI). Details of the code and relevant documentation are freely available to the community in the dedicated repositories., Comment: 13 pages, 7 figures. Accepted for publication in Astronomy & Astrophysics. Public repository with the code can be found at: https://github.com/jfalconbarroso/BAYES-LOSVD

Studying the stellar kinematics of galaxies is a key tool in the reconstruction of their evolution. However, the current measurements of the stellar kinematics are complicated by several factors, including dust extinction and the presence of multiple stellar populations. We use integral field spectroscopic data of four galaxies from the TIMER survey to explore and compare the kinematics measured in different spectral regions that are sensitive to distinct stellar populations. We derive the line-of-sight velocity and velocity dispersion of both a young (<2 Gyr) and an old stellar population from the spectral regions around the H$\beta$ line and the Ca II Triplet. In addition we obtain colour excess, mean age, and metallicity. We report a correlation of the colour excess with the difference in the kinematic parameters of the H$\beta$ line and the Ca II Triplet range, which are dominated by young and old stellar populations, respectively. Young stellar populations, located primarily in nuclear rings, have higher velocity dispersions than old ones. These differences in the rings are typically 10 km/s in velocity dispersion, but up to a mean value of 24 km/s in the most extreme case. Trends with age exist in the nuclear rings but are less significant than those with dust extinction. We report different degrees of correlation of these trends among the galaxies in the sample, which are related to the size of the Voronoi bins in their rings. No clear trends for the difference of line-of-sight velocity are observed. The absence of these trends can be explained as a consequence of the masking process of the H$\beta$ line during the kinematic extraction, as confirmed by dedicated simulations. Our study demonstrates that kinematic differences caused by different stellar populations can be identified in the central regions of nearby galaxies even from intermediate resolution spectroscopy., Comment: 22 pages, 24 figures

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., Comment: 25 p., 19 figures, Accepted to MNRAS

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., Comment: 11 pages, 5 figures, accepted for publication in MNRAS

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., Comment: 16 pages, 14 figures, accepted in A&A (v2 after language editing)

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., Comment: 13 pages, 8 figure; Movie (figure 8) https://youtu.be/vCei7k1bVNE

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., Comment: 24 pages, 18 figures

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., Comment: 17 p., 13 fig., accepted MNRAS

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., Comment: 13 p., 7 Fig., Accepted to MNRAS

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., Comment: 22 pages, 15 figures

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., Comment: Accepted for publication in Astronomy and Astrophysics (A&A), 17 pages, 18 Figures

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., Comment: 14 p., 9 fig., Accepted for publication by MNRAS after minor revisions

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., Comment: Published in MNRAS, 05 January 2018. 17 pages, 9 figures, 4 tables

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., Comment: 15 pages, 16 figures, MNRAS accepted

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., Comment: 17 pages, 9 figures, accepted for publication in Astronomical Notes

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., Comment: 16 pages, 4 main figures main, 7 supplementary figures, 2 supplementary tables; accepted for publication in Nature Astronomy

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., Comment: 8 p., 8 fig., replaced with accepted version (minor changes)

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., Comment: 8 pages, 3 figures, submitted to A&A, version in reply to first referee report

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 $

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., Comment: Interesting new Fig.2. To appear in Astronomische Nachrichten, special issue "Reconstruction the Milky Way's History: Spectroscopic surveys, Asteroseismology and Chemo-dynamical models", Guest Editors C. Chiappini, J. Montalb\'an, and M. Steffen, AN 2016 (in press)"

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., Comment: Submitted to ApJ 13 October 2015

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., Comment: Letter in press in Astronomy and Astrophysics (5 pages, 4 figures + 2 pages Appendix)

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., Comment: 5 pages, 2 figures, Accepted for publication in ApJ Letters

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., Comment: 6 pages, 3 figures & 1 table. Accepted for publication in ApJL

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., Comment: 25 pages (double-column), 17 Figures, under review at A&A

We determine the intrinsic, 3-dimensional shape distribution of star-forming galaxies at 0 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., Comment: Published in ApJ Letters

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., Comment: 9 pages, 5 figures, Accepted to ApJL. Interpretation section expanded, different Fig.5, typos corrected, references added. Conclusions unchanged

[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., Comment: 27 pages, 18 figures, exact match to published version

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., Comment: 18 Pages, 13 Figures, A&A Accepted with minor revisions: a paragraph added in Sec.10 discussing the effect of resolution. Higher resolution version can be found at www.ivanminchev.com

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., Comment: 4 pages, Proceedings of the conference "Galaxies in isolation: exploring Nature vs Nurture" (May 2009, Granada, Spain)

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., Comment: 22 pages, 28 figures, A&A accepted. High resolution version available at http://aramis.obspm.fr/~paola/SFR_frequency/

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., Comment: 10 pages, 5 figures, submitted to PRD