Source author record

Marie Martig

Marie Martig appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2022arXiv

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

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.

preprint2021arXiv

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

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.

preprint2021arXiv

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

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.

preprint2021arXiv

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

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.

preprint2021arXiv

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

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.

preprint2020arXiv

Stellar populations across galaxy bars in the MUSE TIMER project

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.

preprint2016arXiv

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

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.

preprint2015arXiv

CANDELS Visual Classifications: Scheme, Data Release, and First Results

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

preprint2015arXiv

Creating S0s with major mergers: a 3D view

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

preprint2015arXiv

Formation of S0 galaxies through mergers: Explaining angular momentum and concentration change from spirals to S0s

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

preprint2015arXiv

Red giant masses and ages derived from carbon and nitrogen abundances

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

preprint2015arXiv

The Eleventh and Twelfth Data Releases of the Sloan Digital Sky Survey: Final Data from SDSS-III

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

preprint2015arXiv

Young alpha-enriched giant stars in the solar neighbourhood

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

preprint2014arXiv

Chemodynamical evolution of the Milky Way disk II: Variations with Galactic radius and height above the disk plane

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

preprint2014arXiv

Dissecting simulated disc galaxies I: the structure of mono-age populations

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

preprint2014arXiv

Dissecting simulated disc galaxies II: the age-velocity relation

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

preprint2014arXiv

IRAC Near-Infrared Features in the Outer Parts of S4G Galaxies

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

preprint2014arXiv

New Observational Constraints to Milky Way Chemodynamical models

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

preprint2014arXiv

The impact of galaxy geometry and mass evolution on the survival of star clusters

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

preprint2014arXiv

Where do galaxies end? Comparing measurement techniques of hydrodynamic-simulation galaxies' integrated properties

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

preprint2012arXiv

A diversity of progenitors and histories for isolated spiral galaxies

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

preprint2012arXiv

The two-phase formation history of spiral galaxies traced by the cosmic evolution of the bar fraction

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

preprint2010arXiv

Formation of late-type spiral galaxies: gas return from stellar populations regulates disk destruction and bulge growth

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

preprint2009arXiv

Morphological quenching of star formation: making early-type galaxies red

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

preprint2009arXiv

The thick disks of spiral galaxies as relics from gas-rich, turbulent, clumpy disks at high redshift

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

preprint2008arXiv

Numerical simulations of galaxy evolution in cosmological context

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