Showing 1–50 of 87 results for author: Martig, M

Gaia DR3 data consistent with a short bar connected to a spiral arm

Authors: E. Vislosky, I. Minchev, S. Khoperskov, M. Martig, T. Buck, T. Hilmi, B. Ratcliffe, J. Bland-Hawthorn, A. C. Quillen, M. Steinmetz, R. de Jong

Abstract: 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 ca… ▽ More 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. △ Less

Submitted 4 January, 2024; v1 submitted 6 December, 2023; originally announced December 2023.

Comments: 17p, 12 figs. Accepted to MNRAS on Dec 22, 2023

Did the Gaia Enceladus/Sausage merger form the Milky Way's bar?

Authors: Alex Merrow, Robert J. J. Grand, Francesca Fragkoudi, Marie Martig

Abstract: The Milky Way's last significant merger, the Gaia Enceladus/Sausage (GES), is thought to have taken place between 8-11 Gyr ago. Recent studies in the literature suggest that the bar of the Milky Way is rather old, indicating that it formed at a similar epoch to the GES merger. We investigate the possible link between these events using one of the Auriga cosmological simulations which has salient f… ▽ More The Milky Way's last significant merger, the Gaia Enceladus/Sausage (GES), is thought to have taken place between 8-11 Gyr ago. Recent studies in the literature suggest that the bar of the Milky Way is rather old, indicating that it formed at a similar epoch to the GES merger. We investigate the possible link between these events using one of the Auriga cosmological simulations which has salient features in common with the Milky Way, including a last significant merger with kinematic signatures resembling that of the GES. In this simulation, the GES-like merger event triggers tidal forces on the disc, gas inflows and a burst of star formation, with the formation of a bar occuring within 1 Gyr of the first pericentre. To highlight the effects of the merger, we rerun the simulation from z=4 with the progenitors of the GES-like galaxy removed well before the merger time. The consequence is a delay in bar formation by around 2 Gyr, and this new bar forms without any significant external perturbers. We conclude that this Milky Way-like simulation shows a route to the real Milky Way's bar forming around the epoch of the GES merger due to tidal forces on its first pericentre. We explore all Auriga galaxies with GES-like merger events, and find that those with stellar mass ratios below 10% form bars within 1 Gyr of the merger, while bar formation is delayed in the more massive merger scenarios. These include the 4 oldest bars in the simulation suite. Lastly, we note some later morphological differences between the disc of the original simulation and our rerun, in particular that the latter does not grow radially for the final 7 Gyr. Our study suggests that the GES may therefore be responsible for the formation of the Milky Way's bar, as well as for the build-up of its extended disc. △ Less

Submitted 28 March, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: 14 pages, 9 figures, submitted to MNRAS

Stellar populations and origin of thick disks in AURIGA simulations

Authors: Francesca Pinna, Daniel Walo-Martín, Robert J. J. Grand, Marie Martig, Francesca Fragkoudi, Facundo A. Gómez, Federico Marinacci, Rüdiger Pakmor

Abstract: The origin of thick disks and their evolutionary connection with thin disks are still a matter of debate. We provide new insights into this topic by connecting the stellar populations of thick disks at redshift $z=0$ with their past formation and growth, in 24 Milky Way-mass galaxies from the AURIGA zoom-in cosmological simulations. We projected each galaxy edge on, and decomposed it morphological… ▽ More The origin of thick disks and their evolutionary connection with thin disks are still a matter of debate. We provide new insights into this topic by connecting the stellar populations of thick disks at redshift $z=0$ with their past formation and growth, in 24 Milky Way-mass galaxies from the AURIGA zoom-in cosmological simulations. We projected each galaxy edge on, and decomposed it morphologically into two disk components, in order to define geometrically the thin and the thick disks as usually done in observations. We produced age, metallicity and [Mg/Fe] edge-on maps. We quantified the impact of satellite mergers by mapping the distribution of ex-situ stars. Thick disks are on average $\sim 3$~Gyr older, $\sim 0.25$~dex more metal poor and $\sim 0.06$~dex more [Mg/Fe]-enhanced than thin disks. Their average ages range from $\sim 6$ to $\sim 9$~Gyr, metallicities from $\sim -0.15$ to $\sim 0.1$~dex, and [Mg/Fe] from $\sim 0.12$ to $\sim 0.16$~dex. These properties are the result of an early initial in-situ formation, followed by a later growth driven by the combination of direct accretion of stars, some in-situ star formation fueled by mergers, and dynamical heating of stars. The balance between these processes varies from galaxy to galaxy. Mergers play a key role in the mass assembly of thick disks, contributing an average accreted mass fraction of $\sim 22$\% in the analyzed thick-disk dominated regions. In two galaxies, about half of the geometric thick-disk mass was directly accreted. While primordial thick disks form at high redshift in all galaxies, young metal-rich thin disks, with much lower [Mg/Fe] abundances, start to form later but at different times (higher or lower redshift) depending on the galaxy. We conclude that thick disks result from the interplay of external processes with the internal evolution of the galaxy. △ Less

Submitted 2 April, 2024; v1 submitted 22 November, 2023; originally announced November 2023.

Comments: 57 pages, 43 figures, 3 tables; Published in Astronomy & Astrophysics

Journal ref: A&A, 683 (2024) A236

The Milky Way in Context: Building an integral-field spectrograph data cube of the Galaxy

Authors: Zixian Wang, Michael R. Hayden, Sanjib Sharma, Jesse van de Sande, Joss Bland-Hawthorn, Sam Vaughan, Marie Martig, Francesca Pinna

Abstract: The Milky Way (MW) is by far the best-studied galaxy and has been regarded as an ideal laboratory for understanding galaxy evolution. However, direct comparisons of Galactic and extra-galactic observations are marred by many challenges, including selection effects and differences in observations and methodology. In this study, we present a novel code GalCraft to address these challenges by generat… ▽ More The Milky Way (MW) is by far the best-studied galaxy and has been regarded as an ideal laboratory for understanding galaxy evolution. However, direct comparisons of Galactic and extra-galactic observations are marred by many challenges, including selection effects and differences in observations and methodology. In this study, we present a novel code GalCraft to address these challenges by generating mock integral-field spectrograph data cubes of the MW using simple stellar population models and a mock stellar catalog of the Galaxy derived from E-Galaxia. The data products are in the same format as external galaxies, allowing for direct comparisons. We investigate the ability of pPXF to recover kinematics and stellar population properties for an edge-on mock observation of the MW. We confirm that pPXF can distinguish kinematic and stellar population differences between thin and thick disks. However, pPXF struggles to recover star formation history, where the SFR is overestimated in the ranges between 2-4 and 12-14 Gyr compared to the expected values. This is likely due to the template age spacing, pPXF regularization algorithm, and spectral similarities in old population templates. Furthermore, we find systematic offsets in the recovered kinematics, potentially due to insufficient spectral resolution and the variation of line-of-sight velocity with [M/H] and age through a line-of-sight. With future higher resolution and multi-[$α$/Fe] SSP templates, GalCraft will be useful to identify key signatures such as [$α$/Fe]-[M/H] distribution at different $R$ and $|z|$ and potentially measure radial migration and kinematic heating efficiency to study detailed chemodynamical evolution of MW-like galaxies. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: 27 pages, 27 figures (8 figs in appendix), submitted to MNRAS. Comments are welcome

GECKOS: Turning galaxy evolution on its side with deep observations of edge-on galaxies

Authors: J. van de Sande, A. Fraser-McKelvie, D. B. Fisher, M. Martig, M. R. Hayden, the GECKOS Survey collaboration

Abstract: 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… ▽ More 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. △ Less

Submitted 31 May, 2023; originally announced June 2023.

Comments: 6 pages, 2 figures, Proceedings of IAUS 377, eds. F. Tabatabaei, B. Barbuy, and Y. Ting

The vertical structure of the spiral galaxy NGC 3501: first stages of the formation of a thin metal-rich disc

Authors: Natascha Sattler, Francesca Pinna, Nadine Neumayer, Jesus Falcón-Barroso, Marie Martig, Dimitri A. Gadotti, Glenn van de Ven, Ivan Minchev

Abstract: We trace the evolution of the edge-on spiral galaxy NGC 3501, making use of its stellar populations extracted from deep integral-field spectroscopy MUSE observations. We present stellar kinematic and population maps, as well as the star formation history, of the south-western half of the galaxy. The derived maps of the stellar line-of-sight velocity and velocity dispersion are quite regular, show… ▽ More We trace the evolution of the edge-on spiral galaxy NGC 3501, making use of its stellar populations extracted from deep integral-field spectroscopy MUSE observations. We present stellar kinematic and population maps, as well as the star formation history, of the south-western half of the galaxy. The derived maps of the stellar line-of-sight velocity and velocity dispersion are quite regular, show disc-like rotation, and no other structural component of the galaxy. However, maps of the stellar populations exhibit structures in the mass-weighted and light-weighted age, total metallicity and [Mg/Fe] abundance. These maps indicate that NGC 3501 is a young galaxy, consisting mostly of stars with ages between 2 to 8 Gyr. Also, they show a thicker more extended structure that is metal-poor and $α$-rich, and another inner metal-rich and $α$-poor one with smaller radial extension. While previous studies revealed that NGC 3501 shows only one morphological disc component in its vertical structure, we divided the galaxy into two regions: an inner metal-rich midplane and a metal-poor thicker envelope. Comparing the star formation history of the inner thinner metal-rich disc and the thicker metal-poor disc, we see that the metal-rich component evolved more steadily, while the metal-poor one experienced several bursts of star formation. We propose this spiral galaxy is being observed in an early evolutionary phase, with a thicker disc already in place and an inner thin disc in an early formation stage. So we are probably witnessing the birth of a future massive thin disc, continuously growing embedded in a preexisting thicker disc. △ Less

Submitted 23 January, 2023; originally announced January 2023.

Comments: 15 pages, 12 figures. Accepted for publication in MNRAS

A new method for age-dating the formation of bars in disc galaxies: The TIMER view on NGC1433's old bar and the inside-out growth of its nuclear disc

Authors: Camila de Sá-Freitas, Francesca Fragkoudi, Dimitri A. Gadotti, Jesús Falcón-Barroso, Adrian Bittner, Patricia Sánchez-Blázquez, Glenn van de Ven, Rebekka Bieri, Lodovico Coccato, Paula Coelho, Katja Fahrion, Geraldo Gonçalves, Taehyun Kim, Adriana de Lorenzo-Cáceres, Marie Martig, Ignacio Martín-Navarro, Jairo Mendez-Abreu, Justus Neumann, Miguel Querejeta

Abstract: The epoch in which galactic discs settle is a major benchmark to test models of galaxy formation and evolution but is as yet largely unknown. Once discs settle and become self-gravitating enough, stellar bars are able to form; therefore, determining the ages of bars can shed light on the epoch of disc settling, and on the onset of secular evolution. Nevertheless, until now, timing when the bar for… ▽ More The epoch in which galactic discs settle is a major benchmark to test models of galaxy formation and evolution but is as yet largely unknown. Once discs settle and become self-gravitating enough, stellar bars are able to form; therefore, determining the ages of bars can shed light on the epoch of disc settling, and on the onset of secular evolution. Nevertheless, until now, timing when the bar formed has proven challenging. In this work, we present a new methodology for obtaining the bar age, using the star formation history of nuclear discs. Nuclear discs are rotation-supported structures, built by gas pushed to the centre via bar-induced torques, and their formation is thus coincident with bar formation. In particular, we use integral field spectroscopic (IFS) data from the TIMER survey to disentangle the star formation history of the nuclear disc from that of the underlying main disc, which enables us to more accurately determine when the nuclear disc forms. We demonstrate the methodology on the galaxy NGC 1433 -- which we find to host an old bar that is $8.0^{+1.6}_{-1.1}\rm{(sys)}^{+0.2}_{-0.5}\rm{(stat)}$ Gyr old -- and describe a number of tests carried out on both the observational data and numerical simulations. In addition, we present evidence that the nuclear disc of NGC 1433 grows in accordance with an inside-out formation scenario. This methodology is applicable to high-resolution IFS data of barred galaxies with nuclear discs, making it ideally suited for the TIMER survey sample. In the future we will thus be able to determine the bar age for a large sample of galaxies, shedding light on the epoch of disc settling and bar formation. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: Accepted for Publication by A&A, 18 pages, 14 figures

Investigating the Effect of Galaxy Interactions on Star Formation at 0.5<z<3.0

Authors: Ekta A. Shah, Jeyhan S. Kartaltepe, Christina T. Magagnoli, Isabella G. Cox, Caleb T. Wetherell, Brittany N. Vanderhoof, Kevin C. Cooke, Antonello Calabro, Nima Chartab, Christopher J. Conselice, Darren J. Croton, Alexander de la Vega, Nimish P. Hathi, Olivier Ilbert, Hanae Inami, Dale D. Kocevski, Anton M. Koekemoer, Brian C. Lemaux, Lori Lubin, Kameswara Bharadwaj Mantha, Stefano Marchesi, Marie Martig, Jorge Moreno, Belen Alcalde Pampliega, David R. Patton , et al. (2 additional authors not shown)

Abstract: Observations and simulations of interacting galaxies and mergers in the local universe have shown that interactions can significantly enhance the star formation rates (SFR) and fueling of Active Galactic Nuclei (AGN). However, at higher redshift, some simulations suggest that the level of star formation enhancement induced by interactions is lower due to the higher gas fractions and already increa… ▽ More Observations and simulations of interacting galaxies and mergers in the local universe have shown that interactions can significantly enhance the star formation rates (SFR) and fueling of Active Galactic Nuclei (AGN). However, at higher redshift, some simulations suggest that the level of star formation enhancement induced by interactions is lower due to the higher gas fractions and already increased SFRs in these galaxies. To test this, we measure the SFR enhancement in a total of 2351 (1327) massive ($M_*>10^{10}M_\odot$) major ($1<M_1/M_2<4$) spectroscopic galaxy pairs at 0.5<z<3.0 with $ΔV <5000$ km s$^{-1}$ (1000 km s$^{-1}$) and projected separation <150 kpc selected from the extensive spectroscopic coverage in the COSMOS and CANDELS fields. We find that the highest level of SFR enhancement is a factor of 1.23$^{+0.08}_{-0.09}$ in the closest projected separation bin (<25 kpc) relative to a stellar mass-, redshift-, and environment-matched control sample of isolated galaxies. We find that the level of SFR enhancement is a factor of $\sim1.5$ higher at 0.5<z<1 than at 1<z<3 in the closest projected separation bin. Among a sample of visually identified mergers, we find an enhancement of a factor of 1.86$^{+0.29}_{-0.18}$ for coalesced systems. For this visually identified sample, we see a clear trend of increased SFR enhancement with decreasing projected separation (2.40$^{+0.62}_{-0.37}$ vs.\ 1.58$^{+0.29}_{-0.20}$ for 0.5<z<1.6 and 1.6<z<3.0, respectively). The SFR enhancement seen in our interactions and mergers are all lower than the level seen in local samples at the same separation, suggesting that the level of interaction-induced star formation evolves significantly over this time period. △ Less

Submitted 30 September, 2022; originally announced September 2022.

Comments: 23 pages, 13 figures, Accepted for publication in ApJ

Local variations of the Stellar Velocity Ellipsoid-II: the effect of the bar in the inner regions of Auriga galaxies

Authors: Daniel Walo-Martín, Francesca Pinna, Robert J. J. Grand, Isabel Pérez, Jesús Falcón-Barroso, Francesca Fragkoudi, Marie Martig

Abstract: Theoretical works have shown that off-plane motions of bars can heat stars in the vertical direction during buckling but is not clear how do they affect the rest of components of the Stellar Velocity Ellipsoid (SVE). We study the 2D spatial distribution of the vertical, $σ_{z}$, azimuthal, $σ_φ$ and radial, $σ_{r}$ velocity dispersions in the inner regions of Auriga galaxies, a set of high-resolut… ▽ More Theoretical works have shown that off-plane motions of bars can heat stars in the vertical direction during buckling but is not clear how do they affect the rest of components of the Stellar Velocity Ellipsoid (SVE). We study the 2D spatial distribution of the vertical, $σ_{z}$, azimuthal, $σ_φ$ and radial, $σ_{r}$ velocity dispersions in the inner regions of Auriga galaxies, a set of high-resolution magneto-hydrodynamical cosmological zoom-in simulations, to unveil the influence of the bar on the stellar kinematics. $σ_{z}$ and $σ_φ$ maps exhibit non-axisymmetric features that closely match the bar light distribution with low $σ$ regions along the bar major axis and high values in the perpendicular direction. On the other hand, $σ_{r}$ velocity dispersion maps present more axisymmetric distributions. We show that isophotal profile differences best capture the impact of the bar on the three SVE components providing strong correlations with bar morphology proxies although there is no relation with individual $σ$. Time evolution analysis shows that these differences are a consequence of the bar formation and that they tightly coevolve with the strength of the bar. We discuss the presence of different behaviours of $σ_{z}$ and its connection with observations. This work helps us understand the intrinsic $σ$ distribution and motivates the use of isophotal profiles as a mean to quantify the effect of bars. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: 13 pages, 6 Figures. Accepted for publications in MNRAS

The physics governing the upper truncation mass of the globular cluster mass function

Authors: Meghan E. Hughes, Joel L. Pfeffer, Nate Bastian, Marie Martig, J. M. Diederik Kruijssen, Robert A. Crain, Marta Reina-Campos, Sebastian Trujillo-Gomez

Abstract: The mass function of globular cluster (GC) populations is a fundamental observable that encodes the physical conditions under which these massive stellar clusters formed and evolved. The high-mass end of star cluster mass functions are commonly described using a Schechter function, with an exponential truncation mass $M_{c,*}$. For the GC mass functions in the Virgo galaxy cluster, this truncation… ▽ More The mass function of globular cluster (GC) populations is a fundamental observable that encodes the physical conditions under which these massive stellar clusters formed and evolved. The high-mass end of star cluster mass functions are commonly described using a Schechter function, with an exponential truncation mass $M_{c,*}$. For the GC mass functions in the Virgo galaxy cluster, this truncation mass increases with galaxy mass ($M_{*}$). In this paper we fit Schechter mass functions to the GCs in the most massive galaxy group ($M_{\mathrm{200}} = 5.14 \times 10^{13} M_{\odot}$) in the E-MOSAICS simulations. The fiducial cluster formation model in E-MOSAICS reproduces the observed trend of $M_{c,*}$ with $M_{*}$ for the Virgo cluster. We therefore examine the origin of the relation by fitting $M_{c,*}$ as a function of galaxy mass, with and without accounting for mass loss by two-body relaxation, tidal shocks and/or dynamical friction. In the absence of these mass-loss mechanisms, the $M_{c,*}$-$M_{*}$ relation is flat above $M_* > 10^{10} M_{\odot}$. It is therefore the disruption of high-mass GCs in galaxies with $M_{*}\sim 10^{10} M_{\odot}$ that lowers the $M_{c,*}$ in these galaxies. High-mass GCs are able to survive in more massive galaxies, since there are more mergers to facilitate their redistribution to less-dense environments. The $M_{c,*}-M_*$ relation is therefore a consequence of both the formation conditions of massive star clusters and their environmentally-dependent disruption mechanisms. △ Less

Submitted 3 December, 2021; originally announced December 2021.

Comments: Accepted to MNRAS

No memory of past warps in the vertical density structure of galaxies

Authors: J. García de la Cruz, M. Martig, I. Minchev

Abstract: 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 monito… ▽ More 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. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: 10 pages, 7 figures, 1 table. Submitted to MNRAS

NGC 5746: formation history of a massive disc-dominated galaxy

Authors: Marie Martig, Francesca Pinna, Jesús Falcón-Barroso, Dimitri Gadotti, Bernd Husemann, Ivan Minchev, Justus Neumann, Tomás Ruiz-Lara, Glenn van de Ven

Abstract: The existence of massive galaxies lacking a classical bulge has often been proposed as a challenge to $Λ$CDM. However, recent simulations propose that a fraction of massive disc galaxies might have had very quiescent merger histories, and also that mergers do not necessarily build classical bulges. We test these ideas with deep MUSE observations of NGC 5746, a massive ($\sim 10^{11}$ M$_\odot$) ed… ▽ More The existence of massive galaxies lacking a classical bulge has often been proposed as a challenge to $Λ$CDM. However, recent simulations propose that a fraction of massive disc galaxies might have had very quiescent merger histories, and also that mergers do not necessarily build classical bulges. We test these ideas with deep MUSE observations of NGC 5746, a massive ($\sim 10^{11}$ M$_\odot$) edge-on disc galaxy with no classical bulge. We analyse its stellar kinematics and stellar populations, and infer that a massive and extended disc formed very early: 80% of the galaxy's stellar mass formed more than 10 Gyr ago. Most of the thick disc and the bar formed during that early phase. The bar drove gas towards the center and triggered the formation of the nuclear disc followed by the growth of a boxy/peanut-shaped bulge. Around $\sim$ 8 Gyr ago, a $\sim$1:10 merger happened, possibly on a low-inclination orbit. The satellite did not cause significant vertical heating, did not contribute to the growth of a classical bulge, and did not destroy the bar and the nuclear disc. It was however an important event for the galaxy: by depositing its stars throughout the whole galaxy it contributed $\sim 30$% of accreted stars to the thick disc. NGC 5746 thus did not completely escape mergers, but the only relatively recent significant merger did not damage the galaxy and did not create a classical bulge. Future observations will reveal if this is representative of the formation histories of massive disc galaxies. △ Less

Submitted 20 September, 2021; originally announced September 2021.

Comments: Accepted for publication in MNRAS. 22 pages, 21 figures (including appendix)

Local variations of the Stellar Velocity Ellipsoid-I: the disc of galaxies in the Auriga simulations

Authors: Daniel Walo-Martín, Isabel Pérez, Robert J. J. Grand, Jesús Falcón-Barroso, Francesca Pinna, Marie Martig

Abstract: The connection between the Stellar Velocity Ellipsoid (SVE) and the dynamical evolution of galaxies has been a matter of debate in the last years and there is no clear consensus whether different heating agents (e.g. spiral arms, giant molecular clouds, bars and mergers) leave clear detectable signatures in the present day kinematics. Most of these results are based on a single and global SVE and… ▽ More The connection between the Stellar Velocity Ellipsoid (SVE) and the dynamical evolution of galaxies has been a matter of debate in the last years and there is no clear consensus whether different heating agents (e.g. spiral arms, giant molecular clouds, bars and mergers) leave clear detectable signatures in the present day kinematics. Most of these results are based on a single and global SVE and have not taken into account that these agents do not necessarily equally affect all regions of the stellar disc.We study the 2D spatial distribution of the SVE across the stellar discs of Auriga galaxies, a set of high resolution magneto-hydrodynamical cosmological zoom-in simulations, to unveil the connection between local and global kinematic properties in the disc region. We find very similar, global, $σ_{z}/σ_{r}$= 0.80$\pm$ 0.08 values for galaxies of different Hubble types. This shows that the global properties of the SVE at z=0 are not a good indicator of the heating and cooling events experienced by galaxies. We also find that similar $σ_{z}/σ_{r}$radial profiles are obtained through different combinations of $σ_{z}$ and $σ_{r}$ trends: at a local level, the vertical and radial components can evolve differently, leading to similar $σ_{z}/σ_{r}$ profiles at z=0. By contrast, the 2D spatial distribution of the SVE varies a lot more from galaxy to galaxy. Present day features in the SVE spatial distribution may be associated with specific interactions such as fly-by encounters or the accretion of low mass satellites even in the cases when the global SVE is not affected. The stellar populations decomposition reveals that young stellar populations present colder and less isotropic SVEs and more complex 2D distributions than their older and hotter counterparts. △ Less

Submitted 27 July, 2021; v1 submitted 8 June, 2021; originally announced June 2021.

Comments: 13 pages, 9 Figures, Figure 1 compares simulations and observations. Accepted for publication in MNRAS Title Updated

What to expect when using globular clusters as tracers of the total mass distribution in Milky Way-mass galaxies

Authors: Meghan E. Hughes, Prashin Jethwa, Michael Hilker, Glenn van de Ven, Marie Martig, Joel L. Pfeffer, Nate Bastian, J. M. Diederik Kruijssen, Sebastian Trujillo-Gomez, Marta Reina-Campos, Robert A. Crain

Abstract: Dynamical models allow us to connect the motion of a set of tracers to the underlying gravitational potential, and thus to the total (luminous and dark) matter distribution. They are particularly useful for understanding the mass and spatial distribution of dark matter (DM) in a galaxy. Globular clusters (GCs) are an ideal tracer population in dynamical models, since they are bright and can be fou… ▽ More Dynamical models allow us to connect the motion of a set of tracers to the underlying gravitational potential, and thus to the total (luminous and dark) matter distribution. They are particularly useful for understanding the mass and spatial distribution of dark matter (DM) in a galaxy. Globular clusters (GCs) are an ideal tracer population in dynamical models, since they are bright and can be found far out into the halo of galaxies. We aim to test how well Jeans-Anisotropic-MGE (JAM) models using GCs (positions and line-of-sight velocities) as tracers can constrain the mass and radial distribution of DM halos. For this, we use the E-MOSAICS suite of 25 zoom-in simulations of L* galaxies. We find that the DM halo properties are reasonably well recovered by the JAM models. There is, however, a strong correlation between how well we recover the mass and the radial distribution of the DM and the number of GCs in the galaxy: the constraints get exponentially worse with fewer GCs, and at least 150 GCs are needed in order to guarantee that the JAM model will perform well. We find that while the data quality (uncertainty on the radial velocities) can be important, the number of GCs is the dominant factor in terms of the accuracy and precision of the measurements. This work shows promising results for these models to be used in extragalactic systems with a sample of more than 150 GCs. △ Less

Submitted 20 January, 2021; originally announced January 2021.

Comments: 18 pages, 13 figures. Accepted for publication in MNRAS

On the Flaring of Thick Disc of Galaxies: Insights from Simulations

Authors: Joaquín García de la Cruz, Marie Martig, Ivan Minchev, Philip James

Abstract: Using simulated galaxies in their cosmological context, we analyse how the flaring of mono-age populations (MAPs) influences the flaring and the age structure of geometrically-defined thick discs. We also explore under which circumstances the geometric thin and thick discs are meaningfully distinct components, or are part of a single continuous structure as in the Milky Way. We find that flat thic… ▽ More Using simulated galaxies in their cosmological context, we analyse how the flaring of mono-age populations (MAPs) influences the flaring and the age structure of geometrically-defined thick discs. We also explore under which circumstances the geometric thin and thick discs are meaningfully distinct components, or are part of a single continuous structure as in the Milky Way. We find that flat thick discs are created when MAPs barely flare or have low surface density at the radius where they start flaring. When looking at the vertical distribution of MAPs, these galaxies show a continuous thin/thick structure. They also have radial age gradients and tend to have quiescent merger histories. Those characteristics are consistent with what is observed in the Milky Way. Flared thick discs, on the other hand, are created when the MAPs that flare have a high surface density at the radius where they start flaring. The thick discs' scale-heights can either be dominated by multiple MAPs or just a few, depending on the mass and scale-height distribution of the MAPs. In a large fraction of these galaxies, thin and thick discs are clearly distinct structures. Finally, flared thick discs have diverse radial age gradients and merger histories, with galaxies that are more massive or that have undergone massive mergers showing flatter age radial gradients in their thick disc. △ Less

Submitted 28 January, 2021; v1 submitted 4 December, 2020; originally announced December 2020.

Comments: 15 pages, 14 figures. Accepted to MNRAS

BAYES-LOSVD: a bayesian framework for non-parametric extraction of the line-of-sight velocity distribution of galaxies

Authors: J. Falcon-Barroso, M. Martig

Abstract: 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… ▽ More 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. △ Less

Submitted 24 November, 2020; originally announced November 2020.

Comments: 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

Journal ref: A&A 646, A31 (2021)

The kinematics of young and old stellar populations in nuclear rings of MUSE TIMER galaxies

Authors: D. Rosado-Belza, J. Falcón-Barroso, J. H. Knapen, A. Bittner, D. A. Gadotti, J. Neumann, A. de Lorenzo-Cáceres, J. Méndez-Abreu, M. Querejeta, I. Martín-Navarro, P. Sánchez-Blázquez, P. R. T. Coelho, M. Martig, G. van de Ven, T. Kim

Abstract: 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 i… ▽ More 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$β$ 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$β$ 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$β$ 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. △ Less

Submitted 27 October, 2020; v1 submitted 22 October, 2020; originally announced October 2020.

Comments: 22 pages, 24 figures

Journal ref: A&A 644, A116 (2020)

Investigating the Effect of Galaxy Interactions on AGN Enhancement at $0.5<z<3.0$

Authors: Ekta A. Shah, Jeyhan S. Kartaltepe, Christina T. Magagnoli, Isabella G. Cox, Caleb T. Wetherell, Brittany N. Vanderhoof, Antonello Calabro, Nima Chartab, Christopher J. Conselice, Darren J. Croton, Jennifer Donley, Laura de Groot, Alexander de la Vega, Nimish P. Hathi, Olivier Ilbert, Hanae Inami, Dale D. Kocevski, Anton M. Koekemoer, Brian C. Lemaux, Kameswara Bharadwaj Mantha, Stefano Marchesi, Marie Martig, Daniel C. Masters, Elizabeth J. McGrath, Daniel H. McIntosh , et al. (8 additional authors not shown)

Abstract: Galaxy interactions and mergers are thought to play an important role in the evolution of galaxies. Studies in the nearby universe show a higher AGN fraction in interacting and merging galaxies than their isolated counterparts, indicating that such interactions are important contributors to black hole growth. To investigate the evolution of this role at higher redshifts, we have compiled the large… ▽ More Galaxy interactions and mergers are thought to play an important role in the evolution of galaxies. Studies in the nearby universe show a higher AGN fraction in interacting and merging galaxies than their isolated counterparts, indicating that such interactions are important contributors to black hole growth. To investigate the evolution of this role at higher redshifts, we have compiled the largest known sample of major spectroscopic galaxy pairs (2381 with $ΔV <5000$ km s$^{-1}$) at $0.5<z<3.0$ from observations in the COSMOS and CANDELS surveys. We identify X-ray and IR AGN among this kinematic pair sample, a visually identified sample of mergers and interactions, and a mass-, redshift-, and environment-matched control sample for each in order to calculate AGN fractions and the level of AGN enhancement as a function of relative velocity, redshift, and X-ray luminosity. While we see a slight increase in AGN fraction with decreasing projected separation, overall, we find no significant enhancement relative to the control sample at any separation. In the closest projected separation bin ($<25$ kpc, $ΔV <1000$ km s$^{-1}$), we find enhancements of a factor of 0.94$^{+0.21}_{-0.16}$ and 1.00$^{+0.58}_{-0.31}$ for X-ray and IR-selected AGN, respectively. While we conclude that galaxy interactions do not significantly enhance AGN activity on average over $0.5<z<3.0$ at these separations, given the errors and the small sample size at the closest projected separations, our results would be consistent with the presence of low-level AGN enhancement. △ Less

Submitted 6 October, 2020; originally announced October 2020.

Comments: Accepted for publication in ApJ, 28 pages, 15 figures

Kinematic signatures of nuclear discs and bar-driven secular evolution in nearby galaxies of the MUSE TIMER project

Authors: Dimitri A. Gadotti, Adrian Bittner, Jesus Falcon-Barroso, Jairo Mendez-Abreu, Taehyun Kim, Francesca Fragkoudi, Adriana de Lorenzo-Caceres, Ryan Leaman, Justus Neumann, Miguel Querejeta, Patricia Sanchez-Blazquez, Marie Martig, Ignacio Martin-Navarro, Isabel Perez, Marja K. Seidel, Glenn van de Ven

Abstract: The central regions of disc galaxies hold clues to the processes that dominate their formation and evolution. The TIMER project has obtained high signal-to-noise and spatial resolution integral-field spectroscopy data of the inner few kpc of 21 nearby massive barred galaxies, allowing studies of the stellar kinematics with unprecedented spatial resolution. We confirm theoretical predictions of the… ▽ More The central regions of disc galaxies hold clues to the processes that dominate their formation and evolution. The TIMER project has obtained high signal-to-noise and spatial resolution integral-field spectroscopy data of the inner few kpc of 21 nearby massive barred galaxies, allowing studies of the stellar kinematics with unprecedented spatial resolution. We confirm theoretical predictions of the effects of bars on stellar kinematics, and identify box/peanuts through kinematic signatures in mildly and moderately inclined galaxies, finding a lower limit to the fraction of massive barred galaxies with box/peanuts at ~62%. Further, we provide kinematic evidence of the connection between barlenses, box/peanuts and bars. We establish the presence of nuclear discs in 19 galaxies and show that their kinematics are characterised by near-circular orbits with low pressure support, and are consistent with the bar-driven secular evolution picture for their formation. In fact, we show that these nuclear discs have, in the region where they dominate, larger rotational support than the underlying main galaxy disc. We define a kinematic radius for the nuclear discs and show that it relates to bar radius, ellipticity and strength, and bar-to-total ratio. Comparing our results with photometric studies, we find that state-of-the-art galaxy image decompositions are able to discern nuclear discs from classical bulges, if the images employed have enough physical spatial resolution. In fact, we show that nuclear discs are typically identified in such image decompositions as photometric bulges with (near-)exponential profiles. However, we find that the presence of composite bulges (galaxies hosting both a classical bulge and a nuclear disc) can often be unnoticed in studies based on photometry alone, and suggest a more stringent threshold to the Sersic index to identify galaxies with pure classical bulges. △ Less

Submitted 3 September, 2020; originally announced September 2020.

Comments: Accepted for publication in A&A; 24 pages, 9 figures, 3 tables, including 2 appendices. A full resolution PDF is available at http://www.sc.eso.org/~dgadotti/kinematics.pdf and FITS files with the derived (spatially-resolved) stellar kinematics (v, sigma, h3 and h4) for all TIMER galaxies are available at https://www.muse-timer.org/data

Journal ref: A&A 643, A14 (2020)

Weighing stars from birth to death: mass determination methods across the HRD

Authors: Aldo Serenelli, Achim Weiss, Conny Aerts, George C. Angelou, David Baroch, Nate Bastian, Paul G. Beck, Maria Bergemann, Joachim M. Bestenlehner, Ian Czekala, Nancy Elias-Rosa, Ana Escorza, Vincent Van Eylen, Diane K. Feuillet, Davide Gandolfi, Mark Gieles, Leo Girardi, Yveline Lebreton, Nicolas Lodieu, Marie Martig, Marcelo M. Miller Bertolami, Joey S. G. Mombarg, Juan Carlos Morales, Andres Moya, Benard Nsamba , et al. (9 additional authors not shown)

Abstract: 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 approac… ▽ More 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]\,\msun$, $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. △ Less

Submitted 9 April, 2021; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: Invited review article for The Astronomy and Astrophysics Review. 146 pages, 16 figures, 11 tables. Accepted version by the Journal. It includes summary figure of accuracy/precision of methods for mass ranges and summary table for individual methods

Stellar populations across galaxy bars in the MUSE TIMER project

Authors: Justus Neumann, Francesca Fragkoudi, Isabel Pérez, Dimitri A. Gadotti, Jesús Falcón-Barroso, Patricia Sánchez-Blázquez, Adrian Bittner, Bernd Husemann, Facundo A. Gómez, Robert J. J. Grand, Charlotte E. Donohoe-Keyes, Taehyun Kim, Adriana de Lorenzo-Cáceres, Marie Martig, Jairo Méndez-Abreu, Rüdiger Pakmor, Marja K. Seidel, Glenn van de Ven

Abstract: 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 stella… ▽ More 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$α$ 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. △ Less

Submitted 19 March, 2020; originally announced March 2020.

Comments: Accepted for publication in Astronomy and Astrophysics; 14 pages, 14 figures, 3 tables (including appendix: 29 pages, 29 figures)

Journal ref: A&A 637, A56 (2020)

Fluctuations in galactic bar parameters due to bar-spiral interaction

Authors: T. Hilmi, I. Minchev, T. Buck, M. Martig, A. C. Quillen, G. Monari, B. Famaey, R. S. de Jong, C. F. P. Laporte, J. Read, J. L. Sanders, M. Steinmetz, C. Wegg

Abstract: 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%, correlati… ▽ More 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. △ Less

Submitted 30 June, 2020; v1 submitted 11 March, 2020; originally announced March 2020.

Comments: 25 p., 19 figures, Accepted to MNRAS

The loss of the intra-cluster medium in globular clusters

Authors: W. Chantereau, P. Biernacki, M. Martig, N. Bastian, M. Salaris, R. Teyssier

Abstract: 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 simu… ▽ More 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. △ Less

Submitted 6 February, 2020; originally announced February 2020.

Comments: 11 pages, 5 figures, accepted for publication in MNRAS

The [α/Fe]-[Fe/H] relation in the E-MOSAICS simulations: its connection to the birth place of globular clusters and the fraction of globular cluster field stars in the bulge

Authors: Meghan E. Hughes, Joel L. Pfeffer, Marie Martig, Marta Reina-Campos, Nate Bastian, Robert A. Crain, J. M. Diederik Kruijssen

Abstract: The α-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 [α/Fe]-[Fe/H] distribution of field stars and GCs in 25 Milky Way-mass galaxies. The [α/Fe]-[Fe/H] distribution go GCs largely follows that of t… ▽ More The α-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 [α/Fe]-[Fe/H] distribution of field stars and GCs in 25 Milky Way-mass galaxies. The [α/Fe]-[Fe/H] distribution go GCs largely follows that of the field stars and can also therefore be used as tracers of the [α/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 [α/Fe] at fixed [Fe/H]. Therefore, if a GC is observed to have low [α/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 [α/Fe]-[Fe/H] distribution, with a notable subset of galaxies exhibiting bimodal distributions, in which the high [α/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. △ Less

Submitted 5 December, 2019; v1 submitted 3 December, 2019; originally announced December 2019.

Probing the merger history of red early-type galaxies with their faint stellar substructures

Authors: B. Mancillas, P. -A. Duc, F. Combes, F. Bournaud, E. Emsellem, M. Martig, L. Michel-Dansac

Abstract: 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 perturbation… ▽ More 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. △ Less

Submitted 25 October, 2019; v1 submitted 16 September, 2019; originally announced September 2019.

Comments: 16 pages, 14 figures, accepted in A&A (v2 after language editing)

Journal ref: A&A 632, A122 (2019)

Redistribution of Stars and Gas in the Star Formation Deserts of Barred Galaxies

Authors: C. E. Donohoe-Keyes, M. Martig, P. A. James, K. Kraljic

Abstract: 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 distribu… ▽ More 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. △ Less

Submitted 29 August, 2019; originally announced August 2019.

Comments: 13 pages, 8 figure; Movie (figure 8) https://youtu.be/vCei7k1bVNE

Survival of molecular gas in a stellar feedback-driven outflow witnessed with the MUSE TIMER project and ALMA

Authors: Ryan Leaman, Francesca Fragkoudi, Miguel Querejeta, Gigi Y. C. Leung, Dimitri A. Gadotti, Bernd Husemann, Jesus Falcon-Barroso, Patricia Sanchez-Blazquez, Glenn van de Ven, Taehyun Kim, Paula Coelho, Mariya Lyubenova, Adriana de Lorenzo-Caceres, Marie Martig, Inma Martinez-Valpuesta, Justus Neumann, Isabel Perez, Marja Seidel

Abstract: 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… ▽ More 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. △ Less

Submitted 30 July, 2019; originally announced July 2019.

Comments: Accepted for publication in MNRAS. 27 pages, 17 figures

The Fornax 3D project: Thick disks in a cluster environment

Authors: F. Pinna, J. Falcón-Barroso, M. Martig, L. Coccato, E. M. Corsini, P. T. de Zeeuw, D. A. Gadotti, E. Iodice, R. Leaman, M. Lyubenova, I. Martín-Navarro, L. Morelli, M. Sarzi, G. van de Ven, S. Viaene, R. M. McDermid

Abstract: 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 (SF… ▽ More 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. △ Less

Submitted 2 April, 2019; originally announced April 2019.

Comments: 24 pages, 18 figures

Journal ref: A&A 625, A95 (2019)

Kinematics with Gaia DR2: The Force of a Dwarf

Authors: I. Carrillo, I. Minchev, M. Steinmetz, G. Monari, C. F. P. Laporte, F. Anders, A. B. A. Queiroz, C. Chiappini, A. Khalatyan, M. Martig, P. McMillan, B. X. Santiago, K. Youakim

Abstract: 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_φ$, $V_R$ and $V_z$) with Galactic radius, azimuth, and distance from the disc mid-plane. We confirm recent wor… ▽ More 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_φ$, $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<φ<45^\circ$ ahead of the Sun-Galactic centre line, to $-16$ km s$^{-1}$ kpc$^{-1}$ at an azimuth of $-45^\circ<φ<-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. △ Less

Submitted 23 August, 2019; v1 submitted 4 March, 2019; originally announced March 2019.

Comments: 17 p., 13 fig., accepted MNRAS

Yule-Simpson's paradox in Galactic Archaeology

Authors: I. Minchev, G. Matijevic, D. W. Hogg, G. Guiglion, M. Steinmetz, F. Anders, C. Chiappini, M. Martig, A. Queiroz, C. Scannapieco

Abstract: 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 abunda… ▽ More 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. △ Less

Submitted 1 May, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

Comments: 13 p., 7 Fig., Accepted to MNRAS

Dynamical heating across the Milky Way disc using APOGEE and $\it{Gaia}$

Authors: J. Ted Mackereth, Jo Bovy, Henry W. Leung, Ricardo P. Schiavon, Wilma H. Trick, William J. Chaplin, Katia Cunha, Diane K. Feuillet, Steven R. Majewski, Marie Martig, Andrea Miglio, David Nidever, Marc H. Pinsonneault, Victor Silva Aguirre, Jennifer Sobeck, Jamie Tayar, Gail Zasowski

Abstract: 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{[α/Fe]}$ discs between $4 \lesssim R \lesssim 13$ kpc and $|z| \lesssim 2$ kpc using 65,719 stars in common between APOGEE D… ▽ More 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{[α/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{[α/Fe]}$ discs display markedly different age--velocity-dispersion relations (AVRs) and shapes $σ_z/σ_R$. The high $\mathrm{[α/Fe]}$ disc has roughly flat AVRs and constant $σ_z/σ_R = 0.64\pm 0.04$, whereas the low $\mathrm{[α/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{[α/Fe]}$ disc component's flat AVRs and constant $σ_z/σ_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 $σ_z/σ_R$ across the disc indicate that low $\mathrm{[α/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. △ Less

Submitted 30 May, 2019; v1 submitted 14 January, 2019; originally announced January 2019.

Comments: Accepted for publication in MNRAS - The revised version has various updates to the text following suggestions from the referee, but general results remain the same. Code is available https://github.com/jmackereth/monoage-velocity-dispersion, and the APOGEE DR14 age catalogue can be found alongside a paper summary http://www.astro.ljmu.ac.uk/~astjmack/dynamical-heating.html

The Fornax 3D project: unveiling the thick disk origin in FCC 170: signs of accretion?

Authors: F. Pinna, J. Falcón-Barroso, M. Martig, M. Sarzi, L. Coccato, E. Iodice, E. M. Corsini, P. T. de Zeeuw, D. A. Gadotti, R. Leaman, M. Lyubenova, R. M. McDermid, I. Minchev, L. Morelli, G. van de Ven, S. Viaene

Abstract: 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… ▽ More 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. △ Less

Submitted 14 January, 2019; originally announced January 2019.

Comments: 22 pages, 15 figures

Journal ref: A&A 623, A19 (2019)

2D chemical evolution model: the impact of galactic disc asymmetries on azimuthal chemical abundance variations

Authors: E. Spitoni, G. Cescutti, I. Minchev, F. Matteucci, V. Silva Aguirre, M. Martig, G. Bono, C. Chiappini

Abstract: 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 variation… ▽ More 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. △ Less

Submitted 24 June, 2019; v1 submitted 27 November, 2018; originally announced November 2018.

Comments: Accepted for publication in Astronomy and Astrophysics (A&A), 17 pages, 18 Figures

Journal ref: A&A 628, A38 (2019)

Fossil stellar streams and their globular cluster populations in the E-MOSAICS simulations

Authors: Meghan E. Hughes, Joel Pfeffer, Marie Martig, Nate Bastian, Robert A. Crain, J. M. Diederik Kruijssen, Marta Reina-Campos

Abstract: 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… ▽ More 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$. △ Less

Submitted 23 October, 2018; originally announced October 2018.

Comments: Accepted to MNRAS

Time Inference with MUSE in Extragalactic Rings (TIMER): Properties of the Survey and High-Level Data Products

Authors: Dimitri A. Gadotti, Patricia Sánchez-Blázquez, Jesús Falcón-Barroso, Bernd Husemann, Marja K. Seidel, Isabel Pérez, Adriana de Lorenzo-Cáceres, Inma Martinez-Valpuesta, Francesca Fragkoudi, Gigi Leung, Glenn van de Ven, Ryan Leaman, Paula Coelho, Marie Martig, Taehyun Kim, Justus Neumann, Miguel Querejeta

Abstract: 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… ▽ More 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. △ Less

Submitted 2 October, 2018; originally announced October 2018.

Comments: Accepted for publication in MNRAS; 24 pages including 15 colour figures, two tables and an appendix; a version with high resolution figures can be found at http://www.sc.eso.org/~dgadotti/paper_1.pdf

Gas accretion in Milky Way-like galaxies: temporal and radial dependencies

Authors: Sebastian E. Nuza, Cecilia Scannapieco, Cristina Chiappini, Thiago C. Junqueira, Ivan Minchev, Marie Martig

Abstract: 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… ▽ More 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 $Λ$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. △ Less

Submitted 11 December, 2018; v1 submitted 16 May, 2018; originally announced May 2018.

Comments: Replaced to match published version. References and table added; 20 pages, 13 figures, 2 tables. Gas accretion timescale vs. galactocentric radius for simulated MW candidates shown in Table B1

Journal ref: Mon. Not. Roy. Astron. Soc., 482, 3089-3108, 2019

How does the stellar disk of the Milky Way get its gas?

Authors: Sebastian E. Nuza, Cristina Chiappini, Cecilia Scannapieco, Ivan Minchev, Marie Martig, Thiago C. Junqueira

Abstract: 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… ▽ More 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. △ Less

Submitted 3 May, 2018; originally announced May 2018.

Comments: 4 pages, 2 figures. To appear in the Proceedings of IAU Symposium 334,"Rediscovering our Galaxy", Eds. C. Chiappini, I. Minchev, E. Starkenburg & M. Valentini, Cambridge University Press, Cambridge

Estimating stellar birth radii and the time evolution of the Milky Way's ISM metallicity gradient

Authors: I. Minchev, F. Anders, A. Recio-Blanco, C. Chiappini, P. de Laverny, A. Queiroz, M. Steinmetz, V. Adibekyan, I. Carrillo, G. Cescutti, G. Guiglion, M. Hayden, R. S. de Jong, G. Kordopatis, S. R. Majewski, M. Martig, B. X. Santiago

Abstract: 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 an… ▽ More 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. △ Less

Submitted 26 July, 2018; v1 submitted 18 April, 2018; originally announced April 2018.

Comments: 14 p., 9 fig., Accepted for publication by MNRAS after minor revisions

Revisiting the Stellar Velocity Ellipsoid - Hubble type relation: observations versus simulations

Authors: F. Pinna, J. Falcón-Barroso, M. Martig, I. Martínez-Valpuesta, J. Méndez-Abreu, G. van de Ven, R. Leaman, M. Lyubenova

Abstract: 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 e… ▽ More 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 $σ_z/σ_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. △ Less

Submitted 10 January, 2018; originally announced January 2018.

Comments: Published in MNRAS, 05 January 2018. 17 pages, 9 figures, 4 tables

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 475, Issue 2, 1 April 2018, Pages 2697 2712

The stellar orbit distribution in present-day galaxies inferred from the CALIFA survey

Authors: Ling Zhu, Glenn van de Ven, Remco van den Bosch, Hans-Walter Rix, Mariya Lyubenova, Jesús Falcón-Barroso, Marie Martig, Shude Mao, Dandan Xu, Yunpeng Jin, Aura Obreja, Robert J. J. Grand, Aaron A. Dutton, Andrea V. Maccio, Facundo A. Gómez, Jakob C. Walcher, Rubén García-Benito, Stefano Zibetti, Sebastian F. Sánchez

Abstract: 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 $λ_z \simeq1$ and the random-motion dominated o… ▽ More 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 $λ_z \simeq1$ and the random-motion dominated orbits with low circularity $λ_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(λ_z~|~M_\star)$, and its volume-averaged total distribution, $p(λ_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 λ_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. △ Less

Submitted 7 January, 2018; v1 submitted 17 November, 2017; originally announced November 2017.

Comments: Published in Nature Astronomy, 1 January 2018; doi:10.1038/s41550-017-0348-1; 22pages, 8 figures, 2 tables

Confirming chemical clocks: asteroseismic age dissection of the Milky Way disk(s)

Authors: V. Silva Aguirre, M. Bojsen-Hansen, D. Slumstrup, L. Casagrande, D. Kawata, I. Ciuca, R. Handberg, M. N. Lund, J. R. Mosumgaard, D. Huber, J. A. Johnson, M. H. Pinsonneault, A. M. Serenelli, D. Stello, J. Tayar, J. C. Bird, S. Cassisi, M. Hon, M. Martig, P. E. Nissen, H. W. Rix, R. Schönrich, C. Sahlholdt, W. H. Trick, J. Yu

Abstract: 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 c… ▽ More 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-$α$ populations, which also show distinct velocity dispersions in the $V$ and $W$ components. We find no tight correlation between age and metallicity nor [$α$/Fe] for the high-$α$ 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 [$α$/Fe] independent of time. Our findings show that the kinematic properties of young $α$-rich stars are consistent with the rest of the high-$α$ population and different from the low-$α$ 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. △ Less

Submitted 15 January, 2018; v1 submitted 26 October, 2017; originally announced October 2017.

Comments: 15 pages, 16 figures, MNRAS accepted

A unified model for age-velocity dispersion relations in Local Group galaxies: Disentangling ISM turbulence and latent dynamical heating

Authors: Ryan Leaman, J. Trevor Mendel, Emily Wisnioski, Alyson M. Brooks, Michael A. Beasley, Else Starkenburg, Marie Martig, Giuseppina Battaglia, Charlotte Christensen, Andrew A. Cole, T. J. L. de Boer, Drew Wills

Abstract: 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 $σ$ 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 o… ▽ More 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 $σ$ 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 $σ\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. △ Less

Submitted 3 October, 2017; originally announced October 2017.

Comments: 21 pages, 10 figures. accepted for publication in MNRAS

Orbital decomposition of CALIFA spiral galaxies

Authors: Ling Zhu, Remco van den Bosch, Glenn van de Ven, Mariya Lyubenova, Jesús Falcón-Barroso, Sharon E. Meidt, Marie Martig, Juntai Shen, Zhao-Yu Li, Akin Yildirim, C. Jakob Walcher, Sebastian F. Sanchez

Abstract: 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σ$ statistical error of $\sim 10\%$. The lum… ▽ More 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σ$ 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 $σ_r/σ_t$ and $σ_z/σ_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. △ Less

Submitted 19 September, 2017; originally announced September 2017.

Comments: accepted by MNRAS

PLATO as it is: a legacy mission for Galactic archaeology

Authors: A. Miglio, C. Chiappini, B. Mosser, G. R. Davies, K. Freeman, L. Girardi, P. Jofre, D. Kawata, B. M. Rendle, M. Valentini, L. Casagrande, W. J. Chaplin, G. Gilmore, K. Hawkins, B. Holl, T. Appourchaux, K. Belkacem, D. Bossini, K. Brogaard, M. -J. Goupil, J. Montalban, A. Noels, F. Anders, T. Rodrigues, G. Piotto , et al. (80 additional authors not shown)

Abstract: 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 s… ▽ More 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. △ Less

Submitted 7 July, 2017; v1 submitted 12 June, 2017; originally announced June 2017.

Comments: 17 pages, 9 figures, accepted for publication in Astronomical Notes

The unexpectedly large dust and gas content of quiescent galaxies at z>1.4

Authors: R. Gobat, E. Daddi, G. Magdis, F. Bournaud, M. Sargent, M. Martig, S. Jin, A. Finoguenov, M. Béthermin, H. S. Hwang, A. Renzini, G. W. Wilson, I. Aretxaga, M. Yun, V. Strazzullo, F. Valentino

Abstract: 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 o… ▽ More 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. △ Less

Submitted 15 January, 2018; v1 submitted 6 March, 2017; originally announced March 2017.

Comments: 16 pages, 4 main figures main, 7 supplementary figures, 2 supplementary tables; accepted for publication in Nature Astronomy

A radial age gradient in the geometrically thick disk of the Milky Way

Authors: Marie Martig, Ivan Minchev, Melissa Ness, Morgan Fouesneau, Hans-Walter Rix

Abstract: 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 exte… ▽ More 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. △ Less

Submitted 5 September, 2016; originally announced September 2016.

Comments: 6 pages, 4 figures. Accepted for publication in ApJ

The relationship between mono-abundance and mono-age stellar populations in the Milky Way disk

Authors: I. Minchev, M. Steinmetz, C. Chiappini, M. Martig, F. Anders, G. Matijevic, R. S. de Jong

Abstract: 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 equiv… ▽ More 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. △ Less

Submitted 22 November, 2016; v1 submitted 16 August, 2016; originally announced August 2016.

Comments: 8 p., 8 fig., replaced with accepted version (minor changes)

Self-similarity in the chemical evolution of galaxies and the delay time distribution of SNe Ia

Authors: C. J. Walcher, R. M. Yates, I. Minchev, C. Chiappini, M. Bergemann, G. Bruzual, S. Charlot, P. R. T. Coelho, A. Gallazzi, M. Martig

Abstract: 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 α-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 f… ▽ More 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 α-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-[α/Fe] relation as the manifestation of averaging of different stellar populations with varying chemical evolution histories. △ Less

Submitted 30 June, 2016; originally announced July 2016.

Comments: 8 pages, 3 figures, submitted to A&A, version in reply to first referee report

Journal ref: A&A 594, A61 (2016)

Creating lenticular galaxies with major mergers

Authors: Miguel Querejeta, M. Carmen Eliche-Moral, Trinidad Tapia, Alejandro Borlaff, Glenn van de Ven, Mariya Lyubenova, Marie Martig, Jesús Falcón-Barroso, Jairo Méndez-Abreu, Jaime Zamorano, Jesús Gallego

Abstract: 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,… ▽ More 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). △ Less

Submitted 6 May, 2016; originally announced May 2016.

Comments: 3 pages, 1 figure, Proceedings of IAU Symposium 321, "Formation and evolution of galaxy outskirts", Eds. A. Gil de Paz, J. C. Lee & J. H. Knapen, Cambridge University Press, Cambridge

Galactic Archaeology with asteroseismology and spectroscopy: Red giants observed by CoRoT and APOGEE

Authors: F. Anders, C. Chiappini, T. S. Rodrigues, A. Miglio, J. Montalbán, B. Mosser, L. Girardi, M. Valentini, A. Noels, T. Morel, J. A. Johnson, M. Schultheis, F. Baudin, R. de Assis Peralta, S. Hekker, N. Themeßl, T. Kallinger, R. A. García, S. Mathur, A. Baglin, B. X. Santiago, M. Martig, I. Minchev, M. Steinmetz, L. N. da Costa , et al. (20 additional authors not shown)

Abstract: 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 pa… ▽ More 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 [$α$/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 [$α$/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. △ Less

Submitted 17 August, 2016; v1 submitted 26 April, 2016; originally announced April 2016.

Comments: 16 pages + references and appendix, 17 figures, accepted for publication in A&A. Data described in Appendix B will be published through the CDS upon publication

Journal ref: A&A 597, A30 (2017)