Source author record

Amina Helmi

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

ResearcherUnclaimed source record

astro-ph.GA astro-ph.CO astro-ph.SR astro-ph astro-ph.IM astro-ph.HE

Catalog footprint

What is connected

71works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2025arXiv

Chronology of our Galaxy from Gaia colour-magnitude diagram fitting (ChronoGal) -- III. Age and metallicity distribution of Gaia-Sausage-Enceladus stars near the Sun

Context. Gaia-Sausage-Enceladus is considered the last major merger that contributed to the formation of the Milky Way. Its remnants dominate the nearby accreted stellar halo of the Milky Way. Aim. We aim to characterise the star formation history of Gaia-Sausage-Enceladus through the age and metallicity of its stellar populations. Methods. From Gaia DR3 data, we dynamically define three Gaia-Sausage-Enceladus samples with different criteria and possible degrees of contamination from other substructures in the halo. Then, we derive the stellar age and metallicity distributions using the CMDfit.Gaia package. Results. We identify three main populations of stars and a fourth smaller one following an almost linear age-[M/H] relation. The three oldest populations correspond to the bulk of the star formation that lasted for, at least, $\sim$3-4 Gyr and ended about 10 Gyr ago, its metallicities ranging from $-$1.7 to $-$0.8. We categorise these populations into two main epochs: the evolution of GSE in isolation and the merger event. This separation finds independent support from the age-metallicty relation of GSE globular clusters (Aguado-Agelet et al., subm.). The fourth population is younger and more metal-rich, at $\sim$8.5 Gyr and [M/H]$\sim-0.4$; its link to GSE is unclear.

preprint2025arXiv

Chronology of our Galaxy from Gaia colour-magnitude diagram fitting (ChronoGal) II. Unveiling the formation and evolution of the kinematically selected Thick and Thin Discs

Understanding the formation and evolution of the Milky Way's thin and thick discs is crucial to galaxy formation studies. We derive age and metallicity distributions of the kinematic thick and thin discs using the CMDft.Gaia pipeline and Gaia DR3 data within 250 pc of the Sun, covering 1 kpc in height. Our results show that the kinematic thick disc is mostly older than 10 Gyr, undergoing three main metallicity enrichment episodes: (1) over 12 Gyr ago, peaking at [M/H] $\sim$ -0.5 dex, (2) $\sim$11 Gyr ago, rapidly increasing to solar [M/H] and spanning [$α$/Fe] from 0.3 to solar, and (3) just over 10 Gyr ago, reaching supersolar metallicities. Meanwhile, the kinematic thin disc began forming $\sim$10 Gyr ago, just as thick disc star formation ended, characterized by supersolar metallicities and low [$α$/Fe]. This transition coincides with the Milky Way's last major merger: Gaia-Sausage Enceladus (GSE). We also identify a subset of kinematic thin disc stars older than 10 Gyr with high/intermediate [$α$/Fe], indicating a transition phase. The age-metallicity relation of the thin disc suggests overlapping star formation episodes and radial mixing in the solar neighborhood, with the greatest spread $\sim$6 Gyr ago. Additionally, we detect an isolated thick disc star formation event at solar metallicity, coinciding with Sagittarius' first pericenter passage. These findings provide precise age-metallicity distributions and star formation rates, offering key insights for chemical evolution models and cosmological simulations.

preprint2022arXiv

Characterization and dynamics of the peculiar stream Jhelum, A tentative role for the Sagittarius dwarf galaxy

Stellar streams are a promising tool to study the Milky Way's dark matter subhalo population, as interactions with subhalos could produce substructure in streams. However, other possible causes for substructure first need to be well understood. Here we studied the kinematics and the unusual morphology of the stellar stream Jhelum. Using a combination of ground-based photometry and Gaia EDR3 astrometry, we characterized the morphology of Jhelum. We combined this new data with radial velocities from the literature to perform orbit integrations of the stream in static Galactic potentials. We also carried out N-body simulations in the presence of the Sagittarius dwarf galaxy. The new data reveal a previously unreported tertiary component in the stream, as well as several gaps and a kink-like feature in its narrow component. We find that for a range of realistic Galactic potentials, no single orbit is able to reproduce Jhelum's radial velocity data entirely. A generic property of the orbital solutions is that they share a similar orbital plane to Sagittarius and this leads to repeated encounters with the stream. Using N-body simulations that include a massive Sagittarius, we explored its effect on Jhelum, and we show that these encounters can qualitatively reproduce the narrow and broad components in Jhelum, as well as create a tertiary component in some cases. We also find evidence that such encounters can result in an apparent increase in the velocity dispersion of the stream by a factor up to four due to overlapping narrow and broad components. Our findings suggest that the Jhelum stream is even more complex than once thought; however, its morphology and kinematics can tentatively be explained via the interactions with Sagittarius. In this scenario, the formation of Jhelum's narrow and broad components occurs naturally, yet some of the smaller gap-like features remain to be explained.

preprint2022arXiv

Substructure in the stellar halo near the Sun. I. Data-driven clustering in Integrals of Motion space

Aims: Develop a data-driven and statistically based method for finding such clumps in Integrals of Motion space for nearby halo stars and evaluating their significance robustly. Methods: We use data from Gaia EDR3 extended with radial velocities from ground-based spectroscopic surveys to construct a sample of halo stars within 2.5 kpc from the Sun. We apply a hierarchical clustering method that uses the single linkage algorithm in a 3D space defined by the commonly used integrals of motion energy $E$, together with two components of the angular momentum, $L_z$ and $L_\perp$. To evaluate the statistical significance of the clusters found, we compare the density within an ellipsoidal region centered on the cluster to that of random sets with similar global dynamical properties. We pick out the signal at the location of their maximum statistical significance in the hierarchical tree. We estimate the proximity of a star to the cluster center using the Mahalanobis distance. We also apply the HDBSCAN clustering algorithm in velocity space. Results: Our procedure identifies 67 highly significant clusters ($ > 3σ$), containing 12\% of the sources in our halo set, and in total 232 subgroups or individual streams in velocity space. In total, 13.8\% of the stars in our data set can be confidently associated to a significant cluster based on their Mahalanobis distance. Inspection of our data set reveals a complex web of relationships between the significant clusters, suggesting that they can be tentatively grouped into at least 6 main structures, many of which can be associated to previously identified halo substructures, and a number of independent substructures. This preliminary conclusion is further explored in an accompanying paper by Ruiz-Lara et al., where we also characterize the substructures in terms of their stellar populations. Conclusions: We find... (abridged version)

preprint2022arXiv

Substructure in the stellar halo near the Sun. II. Characterisation of independent structures

In Lövdal et al, we presented a data-driven method for clustering in Integrals of Motion space and applied it to a large sample of nearby halo stars with 6D phase-space information. We identified a large number of clusters, many of which could tentatively be merged into larger groups. Our goal is to establish the reality of the clusters through a combined study of their stellar populations to gain more insights into the accretion history of the Milky Way. We develop a procedure that quantifies the similarity of clusters based on KS tests using their metallicity distribution functions, and an isochrone fitting method to determine their average age, which is also used to compare the distribution of stars in the Colour-Absolute magnitude diagram. This allows us to group clusters into substructures, and to compare substructures with one another. The clusters identified are merged into 12 extended substructures, while 8 small clusters remain as such. The large substructures include the previously known Gaia-Enceladus, Helmi streams, Sequoia, and Thamnos 1 and 2. We identify overdensities associated with the hot thick disc and hosting a metal-poor population. Especially notable is our largest substructure which, although peaking at the metallicity characteristic of the thick disk has a well populated metal-poor component, and dynamics in-between hot thick disc and halo. We identify additional debris in the region occupied by Sequoia with distinct kinematics, likely remnants of three different accretion events with progenitors of similar mass. We also identify different trends of [Mg/Fe] vs [Fe/H] for the various substructures confirming our dissection of the nearby halo. At least 20\% of the halo near the Sun is associated to substructures. When comparing their global properties, we note that those substructures on retrograde orbits are not only more metal-poor on average but also older.

preprint2021arXiv

Age Determination of Galaxy Merger Remnant Stars using Asteroseismology

The Milky Way was shaped by the mergers with several galaxies in the past. We search for remnant stars that were born in these foreign galaxies and assess their ages in an effort to put upper limits on the merger times and thereby better understand the evolutionary history of our Galaxy. Using 6D-phase space information from Gaia eDR3 and chemical information from APOGEE DR16, we kinematically and chemically select $23$ red giant stars belonging to former dwarf galaxies that merged with the Milky Way. With added asteroseismology from Kepler and K2, we determine the ages of the $23$ ex-situ stars and $55$ in-situ stars with great precision. We find that all the ex-situ stars are consistent with being older than $8$ Gyr. While it is not possible to associate all the stars with a specific dwarf galaxy we classify eight of them as Gaia-Enceladus/Sausage stars, which is one of the most massive mergers in our Galaxy's history. We determine their mean age to be $9.5^{+1.2}_{-1.3}$ Gyr consistent with a merger time of $8$-$10$ Gyr ago. The rest of the stars are possibly associated with Kraken, Thamnos, Sequoia, or another extragalactic progenitor. The age determination of ex-situ stars paves the way to more accurately pinning down when the merger events occurred and hence provide tight constraints useful for simulating how these events unfolded.

preprint2021arXiv

Galactic potential constraints from clustering in action space of combined stellar stream data

Stream stars removed by tides from their progenitor satellite galaxy or globular cluster act as a group of test particles on neighboring orbits, probing the gravitational field of the Milky Way. While constraints from individual streams have been shown to be susceptible to biases, combining several streams from orbits with various distances reduces these biases. We fit a common gravitational potential to multiple stellar streams simultaneously by maximizing the clustering of the stream stars in action space. We apply this technique to members of the GD-1, Pal 5, Orphan and Helmi streams, exploiting both the individual and combined data sets. We describe the Galactic potential with a Stäckel model, and vary up to five parameters simultaneously. We find that we can only constrain the enclosed mass, and that the strongest constraints come from the GD-1, Pal 5 and Orphan streams whose combined data set yields $M(< 20\ \mathrm{kpc}) = 2.96^{+0.25}_{-0.26} \times 10^{11} \ M_{\odot}$. When including the Helmi stream in the data set, the mass uncertainty increases to $M(< 20\ \mathrm{kpc}) = 3.12^{+3.21}_{-0.46} \times 10^{11} \ M_{\odot}$.

preprint2021arXiv

Substructures, Resonances and debris streams

The local stellar halo of the Milky Way contains the debris from several past accretion events. Here we study in detail the structure and properties of nearby debris associated with the Helmi streams, originally identified as an overdensity in integrals of motion space. We use 6D phase-space information from Gaia EDR3 combined with spectroscopic surveys, and we analyse the orbits and frequencies of the stars in the streams using various Galactic potentials. We also explore how the Helmi streams constrain the flattening q, of the Galactic dark matter halo. We find that the streams are split into substructures in integrals of motion space, most notably into two clumps in angular momentum space. The clumps have consistent metallicity distributions and stellar populations, supporting a common progeny. In all the realistic Galactic potentials explored, the Helmi streams stars depict a diffuse distribution close to Omega_z /Omega_R ~0.7. At the same time, the reason for the substructure in angular momentum space appears to be a Omega_z : Omega_phi resonance close to the 1:1. This resonance is exactly the 1:1 in the case that the (density) flattening of the dark halo is q = 1.2. For this halo shape the substructure in angular momenta is also long lasting. Our findings suggest that the structure of the Galactic potential leaves a clear imprint on the properties of phase-mixed debris streams.

preprint2020arXiv

Cosmological insights into the assembly of the radial and compact stellar halo of the Milky Way

Recent studies using Gaia DR2 have identified a massive merger in the history of the Milky Way (MW) whose debris is markedly radial and counterrotating. This event, known as the Gaia-Enceladus/Gaia-Sausage (GE/GS), is also hypothesized to have built the majority of the inner stellar halo. We use the cosmological hydrodynamic simulation Illustris to place this merger in the context of galaxy assembly within $Λ$CDM. From $\sim$150 MW analogs, $\sim 80 \%$ have experienced at least one merger of similar mass and infall time as GE/GS. Within this sample, 37 have debris as radial as that of the GE/GS, which we dub the Ancient Radial Mergers (ARMs). Counterrotation is not rare among ARMs, with $43 \%$ having $>40 \%$ of their debris in counterrotating orbits. However, the compactness inferred for the GE/GS debris, given its large $β$ and its substantial contribution to the stellar halo, is more difficult to reproduce. The median radius of ARM debris is r$_{*,deb}\simeq 45$kpc, while GE/GS is thought to be mostly contained within $r\sim 30$ kpc. For most MW analogs, a few mergers are required to build the inner stellar halo, and ARM debris only accounts for $\sim 12 \%$ of inner accreted stars. Encouragingly, we find one ARM that is both compact and dominates the inner halo of its central, making it our best GE/GS analog. Interestingly, this merger deposits a significant number of stars (M$_*\simeq1.5 \times 10^9 M_\odot$) in the outer halo, suggesting that an undiscovered section of GE/GS may await detection.

preprint2020arXiv

Frequencies, chaos and resonances: a study of orbital parameters of nearby thick disc and halo stars

We study the distribution of nearby thick disc and halo stars in subspaces defined by their characteristic orbital parameters. Our aim is to establish the origin of the structure reported in particular, in the $R_{\rm max}-z_{\rm max}$ space. To this end, we compute the orbital parameters and frequencies of stars for a generic and for a Stäckel Milky Way potential. We find that for both the thick disc and halo populations very similar prominent substructures are apparent for the generic Galactic potential, while no substructure is seen for the Stäckel model. This indicates that the origin of these features is not merger-related, but due to non-integrability of the generic potential. This conclusion is strengthened by our frequency analysis of the orbits of stars, which reveals the presence of prominent resonances, with $\sim 30\%$ of the halo stars associated to resonance families. In fact, the stars in resonances define the substructures seen in the spaces of characteristic orbital parameters. Furthermore, we find that some stars in our sample and in debris streams are on the same resonance as the Sagittarius dwarf, suggesting this system has also influenced the distribution of stars in the Galactic thick disc and halo components. Our study constitutes a step towards disentangling the imprint of merger debris from substructures driven by internal dynamics. Given their prominence, these resonant-driven overdensities could potentially be useful to constrain the exact form of the Galactic potential.

preprint2020arXiv

Streams, substructures and the early history of the Milky Way

The advent of Gaia's 2nd data release in combination with large spectroscopic surveys are revolutionizing our understanding of the Galaxy. Thanks to these and the knowledge accumulated thus far, a more mature picture of the evolution of the early Milky Way is emerging: * Two of the traditional Galactic components, i.e. the stellar halo and the thick disk, appear to be intimately linked: stars with halo-like kinematics originate in similar proportions, from a "heated" (thick) disk and from debris from a system named Gaia-Enceladus. Gaia-Enceladus was the last big merger event experienced by the Milky Way and probably completed around 10 Gyr ago. The puffed-up stars now present in the halo as a consequence of the merger have thus exposed the existence of a disk component at z ~ 1.8. * The Helmi streams, Sequoia, and Thamnos are amongst the newly uncovered or better characterized merger events. Knowledge of their progenitor's properties, star formation and chemical histories is still incomplete. * Debris' from different objects often overlap in phase-space. A task for the next years will be to use spectroscopic surveys for chemical labelling and to disentangle events from one another using dimensions other than only phase-space, metallicity or [alpha/Fe]. * These surveys will also provide line-of-sight velocities missing for faint stars and more accurate distance determinations for distant objects. The resulting samples of stars will cover a much wider volume of the Galaxy allowing, for example, linking kinematic substructures in the inner halo to spatial overdensities in the outer halo. * All the results obtained so far are in-line with expectations of current cosmological models. Yet, tailored hydrodynamical simulations as well as "constrained" cosmological simulations are needed to push our knowledge of the assembly of the Milky Way back to the earliest times. [abridged]

preprint2020arXiv

The Sixth Data Release of the Radial Velocity Experiment (RAVE) -- I: Survey Description, Spectra and Radial Velocities

The Radial Velocity Experiment (RAVE) is a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A). The 6th and final data release (DR6 or FDR) is based on 518387 observations of 451783 unique stars. RAVE observations were taken between 12 April 2003 and 4 April 2013. Here we present the genesis, setup and data reduction of RAVE as well as wavelength-calibrated and flux-normalized spectra and error spectra for all observations in RAVE DR6. Furthermore, we present derived spectral classification and radial velocities for the RAVE targets, complemented by cross matches with Gaia DR2 and other relevant catalogs. A comparison between internal error estimates, variances derived from stars with more than one observing epoch and a comparison with radial velocities of Gaia DR2 reveals consistently that 68% of the objects have a velocity accuracy better than 1.4 km/s, while 95% of the objects have radial velocities better than 4.0 km/s. Stellar atmospheric parameters, abundances and distances are presented in subsequent publication. The data can be accessed via the RAVE Web (http://rave-survey.org) or the Vizier database.

preprint2020arXiv

The Sixth Data Release of the Radial Velocity Experiment (RAVE) -- II: Stellar Atmospheric Parameters, Chemical Abundances and Distances

We present part 2 of the 6th and final Data Release (DR6 or FDR) of the Radial Velocity Experiment (RAVE), a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A) and span the complete time frame from the start of RAVE observations on 12 April 2003 to their completion on 4 April 2013. In the second of two publications, we present the data products derived from 518387 observations of 451783 unique stars using a suite of advanced reduction pipelines focussing on stellar atmospheric parameters, in particular purely spectroscopically derived stellar atmospheric parameters (Teff, log(g), and the overall metallicity), enhanced stellar atmospheric parameters inferred via a Bayesian pipeline using Gaia DR2 astrometric priors, and asteroseismically calibrated stellar atmospheric parameters for giant stars based on asteroseismic observations for 699 K2 stars. In addition, we provide abundances of the elements Fe, Al, and Ni, as well as an overall [alpha/Fe] ratio obtained using a new pipeline based on the GAUGUIN optimization method that is able to deal with variable signal-to-noise ratios. The RAVE DR6 catalogs are cross matched with relevant astrometric and photometric catalogs, and are complemented by orbital parameters and effective temperatures based on the infrared flux method. The data can be accessed via the RAVE Web site (http://rave-survey.org) or the Vizier database.

preprint2019arXiv

Stellar 3-D kinematics in the Draco dwarf spheroidal galaxy

Aims. We present the first three-dimensional internal motions for individual stars in the Draco dwarf spheroidal galaxy. Methods. By combining first-epoch $Hubble$ $Space$ $Telescope$ observations and second-epoch $Gaia$ Data Release 2 positions, we measured the proper motions of $149$ sources in the direction of Draco. We determined the line-of-sight velocities for a sub-sample of $81$ red giant branch stars using medium resolution spectra acquired with the DEIMOS spectrograph at the Keck II telescope. Altogether, this resulted in a final sample of $45$ Draco members with high-precision and accurate 3D motions, which we present as a table in this paper. Results. Based on this high-quality dataset, we determined the velocity dispersions at a projected distance of $\sim120$ pc from the centre of Draco to be $σ_{R} =11.0^{+2.1}_{-1.5}$ km/s, $σ_{T}=9.9^{+2.3}_{-3.1}$ km/s and $σ_{LOS}=9.0^{+1.1}_{-1.1}$ km/s in the projected radial, tangential, and line-of-sight directions. This results in a velocity anisotropy $β=0.25^{+0.47}_{-1.38}$ at $r \gtrsim120$ pc. Tighter constraints may be obtained using the spherical Jeans equations and assuming constant anisotropy and Navarro-Frenk-White (NFW) mass profiles, also based on the assumption that the 3D velocity dispersion should be lower than $\approx 1/3$ of the escape velocity of the system. In this case, we constrain the maximum circular velocity $V_{max}$ of Draco to be in the range of $10.2-17.0$ km/s. The corresponding mass range is in good agreement with previous estimates based on line-of-sight velocities only. Conclusions. Our Jeans modelling supports the case for a cuspy dark matter profile in this galaxy. Firmer conclusions may be drawn by applying more sophisticated models to this dataset and with new datasets from upcoming $Gaia$ releases.

preprint2016arXiv

Building Blocks of the Milky Way's Accreted Spheroid

In the $Λ$CDM model of structure formation, a stellar spheroid grows by the assembly of smaller galaxies, the so-called building blocks. Combining the Munich-Groningen semi-analytical model of galaxy formation with the high resolution Aquarius simulations of dark matter haloes, we study the assembly history of the stellar spheroids of six Milky Way-mass galaxies, focussing on building block properties such as mass, age and metallicity. These properties are compared to those of the surviving satellites in the same models. We find that the building blocks have higher star formation rates on average, and this is especially the case for the more massive objects. At high redshift these dominate in star formation over the satellites, whose star formation timescales are longer on average. These differences ought to result in a larger $α$-element enhancement from Type II supernovae in the building blocks (compared to the satellites) by the time Type Ia supernovae would start to enrich them in iron, explaining the observational trends. Interestingly, there are some variations in the star formation timescales of the building blocks amongst the simulated haloes, indicating that [$α$/Fe] abundances in spheroids of other galaxies might differ from those in our own Milky Way.

preprint2016arXiv

Dark influences III. Structural characterization of minor mergers of dwarf galaxies with dark satellites

In the current concordance cosmology small halos are expected to be completely dark and can significantly perturb low-mass galaxies during minor merger interactions. These interactions may well contribute to the diversity of the dwarf galaxy population. Dwarf galaxies in the field are often observed to have peculiarities in their structure, morphology, and kinematics, as well as strong bursts of star formation without apparent cause. We aim to characterize the signatures of minor mergers of dwarf galaxies with dark satellites to aid their observational identification. We explore and quantify a variety of structural, morphological, and kinematic indicators of merging dwarf galaxies and their remnants using a suite of hydrodynamical simulations. The most sensitive indicators of mergers with dark satellites are large asymmetries in the gaseous and stellar distributions, enhanced central surface brightness and starbursts, and velocity offsets and misalignments between the cold gas and stellar components. In general, merging systems span a wide range of values of the most commonly used indicators, while isolated objects tend to have more confined values. Interestingly, we find in our simulations that a significantly off-centered burst of star formation can pinpoint the location of the dark satellite. Observational systems with such characteristics are perhaps the most promising for unveiling the presence of the hitherto, missing satellites.

preprint2016arXiv

Stream-subhalo interactions in the Aquarius simulations

We perform the first self-consistent measurement of the rate of interactions between stellar tidal streams created by disrupting satellites and dark subhalos in a cosmological simulation of a Milky-Way-mass galaxy. Using a retagged version of the Aquarius A dark-matter-only simulation, we selected 18 streams of tagged star particles that appear thin at the present day and followed them from the point their progenitors accrete onto the main halo, recording in each snapshot the characteristics of all dark-matter subhalos passing within several distance thresholds of any tagged star particle in each stream. We considered distance thresholds corresponding to constant impact parameters (1, 2, and 5 kpc), as well as those proportional to the region of influence of each subhalo (one and two times its half-mass radius $r_{1/2}$). We then measured the age and present-day, phase-unwrapped length of each stream in order to compute the interaction rate in different mass bins and for different thresholds, and compared these to analytic predictions from the literature. We measure a median rate of $1.5^{+3.0}_{-1.1}\ (9.1^{+17.5}_{-7.1},\ 61.8^{+211}_{-40.6})$ interactions within 1 (2, 5) kpc of the stream per 10 kpc of stream length per 10 Gyr. Resolution effects (both time and particle number) affect these estimated rates by lowering them.

preprint2016arXiv

The kinematics of globular clusters systems in the outer halos of the Aquarius simulations

Stellar halos and globular cluster (GC) systems contain valuable information regarding the assembly history of their host galaxies. Motivated by the detection of a significant rotation signal in the outer halo GC system of M31, we investigate the likelihood of detecting such a rotation signal in projection, using cosmological simulations. To this end we select subsets of tagged particles in the halos of the Aquarius simulations to represent mock GC systems, and analyse their kinematics. We find that GC systems can exhibit a non-negligible rotation signal provided the associated stellar halo also has a net angular momentum. The ability to detect this rotation signal is highly dependent on the viewing perspective, and the probability of seeing a signal larger than that measured in M31 ranges from 10% to 90% for the different halos in the Aquarius suite. High values are found from a perspective such that the projected angular momentum of the GC system is within 40 deg of the rotation axis determined via the projected positions and line-of-sight velocities of the GCs. Furthermore, the true 3D angular momentum of the outer stellar halo is relatively well aligned, within 35 deg, with that of the mock GC systems. We argue that the net angular momentum in the mock GC systems arises naturally when the majority of the material is accreted from a preferred direction, namely along the dominant dark matter filament of the large-scale structure that the halos are embedded in. This, together with the favourable edge-on view of M31's disk suggests that it is not a coincidence that a large rotation signal has been measured for its outer halo GC system.

preprint2016arXiv

The time evolution of gaps in tidal streams

We model the time evolution of gaps in tidal streams caused by the impact of a dark matter subhalo, while both orbit a spherical gravitational potential. To this end, we make use of the simple behaviour of orbits in action-angle space. A gap effectively results from the divergence of two nearby orbits whose initial phase-space separation is, for very cold thin streams, largely given by the impulse induced by the subhalo. We find that in a spherical potential the size of a gap increases linearly with time for sufficiently long timescales. We have derived an analytic expression that shows how the growth rate depends on the mass of the perturbing subhalo, its scale and its relative velocity with respect to the stream. We have verified these scalings using N-body simulations and find excellent agreement. For example, a subhalo of mass 10^8 Msun directly impacting a very cold thin stream on an inclined orbit can induce a gap that may reach a size of several tens of kpc after a few Gyr. The gap size fluctuates importantly with phase on the orbit, and it is largest close to pericentre. This indicates that it may not be fully straightforward to invert the spectrum of gaps present in a stream to recover the mass spectrum of the subhalos.

preprint2015arXiv

Action-space clustering of tidal streams to infer the Galactic potential

We present a new method for constraining the Milky Way halo gravitational potential by simultaneously fitting multiple tidal streams. This method requires full three-dimensional positions and velocities for all stars to be fit, but does not require identification of any specific stream or determination of stream membership for any star. We exploit the principle that the action distribution of stream stars is most clustered when the potential used to calculate the actions is closest to the true potential. Clustering is quantified with the Kullback-Leibler Divergence (KLD), which also provides conditional uncertainties for our parameter estimates. We show, for toy Gaia-like data in a spherical isochrone potential, that maximizing the KLD of the action distribution relative to a smoother distribution recovers the true values of the potential parameters. The precision depends on the observational errors and the number of streams in the sample; using KIII giants as tracers, we measure the enclosed mass at the average radius of the sample stars accurate to 3% and precise to 20-40%. Recovery of the scale radius is precise to 25%, and is biased 50% high by the small galactocentric distance range of stars in our mock sample (1-25 kpc, or about three scale radii, with mean 6.5 kpc). About 15 streams, with at least 100 stars per stream, are needed to obtain upper and lower bounds on the enclosed mass and scale radius when observational errors are taken into account; 20-25 streams are required to stabilize the size of the confidence interval. If radial velocities are provided for stars out to 100 kpc (10 scale radii), all parameters can be determined with 10% accuracy and 20% precision (1.3% accuracy in the case of the enclosed mass), underlining the need for ground-based spectroscopic follow-up to complete the radial velocity catalog for faint halo stars observed by Gaia.

preprint2015arXiv

Dark influences II: gas and star formation in minor mergers of dwarf galaxies with dark satellites

It has been proposed that mergers induce starbursts and lead to important morphological changes in galaxies. Most studies so far have focused on large galaxies, but dwarfs might also experience such events, since the halo mass function is scale-free in the concordance cosmological model. Notably, because of their low mass, most of their interactions will be with dark satellites. In this paper we follow the evolution of gas-rich disky dwarf galaxies as they experience a minor merger with a dark satellite. We aim to characterize the effects of such an interaction on the dwarf's star formation, morphology, and kinematical properties. We performed a suite of carefully set-up hydrodynamical simulations of dwarf galaxies that include dark matter, gas, and stars merging with a satellite consisting solely of dark matter. For the host system we vary the gas fraction, disk size and thickness, halo mass, and concentration, while we explore different masses, concentrations, and orbits for the satellite. We find that the interactions cause strong starbursts of both short and long duration in the dwarfs. Their star formation rates increase by factors of a few to 10 or more. They are strongest for systems with extended gas disks and high gas fractions merging with a high-concentration satellite on a planar, radial orbit. In contrast to analogous simulations of Milky Way-mass galaxies, many of the systems experience strong morphological changes and become spheroidal even in the presence of significant amounts of gas. The simulated systems compare remarkably well with the observational properties of a large selection of irregular dwarf galaxies and blue compact dwarfs. This implies that mergers with dark satellites might well be happening but not be fully evident, and may thus play a role in the diversity of the dwarf galaxy population.

preprint2015arXiv

Matching the dark matter profiles of dSph galaxies with those of simulated satellites: a two parameter comparison

We compare the dark matter halos' structural parameters derived for four Milky Way dwarf spheroidal galaxies to those of subhalos found in cosmological $N$-body simulations. We confirm that estimates of the mass at a single fixed radius are fully consistent with the observations. However, when a second structural parameter such as the logarithmic slope of the dark halo density profile measured close to the half-light radius is included in the comparison, we find little to no overlap between the satellites and the subhalos. Typically the right mass subhalos have steeper profiles at these radii than measurements of the dSph suggest. Using energy arguments we explore if it is possible to solve this discrepancy by invoking baryonic effects. Assuming that feedback from supernovae can lead to a reshaping of the halos, we compute the required efficiency and find entirely plausible values for a significant fraction of the subhalos and even as low as 0.1%. This implies that care must be taken not to exaggerate the effect of supernovae feedback as this could make the halos too shallow. These results could be used to calibrate and possibly constrain feedback recipes in hydrodynamical simulations.

preprint2015arXiv

The evolution of streams in a time-dependent potential

We study the evolution of streams in a time-dependent spherical gravitational potential. Our goal is to establish what are the imprints of this time evolution on the properties of streams as well as their observability. To this end, we have performed a suite of numerical experiments for a host system that doubles its mass during the integration time and for a variety of initial conditions. In these experiments we found that the most striking imprint is a misalignment of 10 degrees in the angular location of the apocentres of the streams compared to the static case (and to the orbit of the centre of mass), which only becomes apparent for sufficiently long streams. We have also developed an analytic model using action-angle variables which allows us to explain this behaviour and to identify the most important signature of time evolution, namely a difference in the slope defined by the distribution of particles along a stream in frequency and in angle space. Although a difference in slope can arise when the present-day potential is not correctly modelled, this shortcoming can be by-passed because in this case, streams are no longer straight lines in angle space, but depict a wiggly appearance and an implausible energy gradient. The difference in slope due to time-evolution is small, typically $\sim10^{-2}$ and its amplitude depends on the growth rate of the potential, but nonetheless we find that it could be observable if accurate full-space information for nearby long streams is available.

preprint2015arXiv

What can Gaia proper motions tell us about Milky Way dwarf galaxies?

We present a proper-motion study on models of the dwarf spheroidal galaxy Sculptor, based on the predicted proper-motion accuracy of Gaia measurements. Gaia will measure proper motions of several hundreds of stars for a Sculptor-like system. Even with an uncertainty on the proper motion of order 1.5 times the size of an individual proper-motion value of ~10 mas/century, we find that it is possible to recover Sculptor's systemic proper motion at its distance of 79 kpc.

preprint2014arXiv

A new fitting-function to describe the time evolution of a galaxy's gravitational potential

We present a new simple functional form to model the evolution of a spherical mass distribution in a cosmological context. Two parameters control the growth of the system and this is modelled using a redshift dependent exponential for the scale mass and scale radius. In this new model, systems form inside out and the mass of a given shell can be made to never decrease, as generally expected. This feature makes it more suitable for studying the smooth growth of galactic potentials or cosmological halos than other parametrizations often used in the literature. This is further confirmed through a comparison to the growth of dark matter halos in the Aquarius simulations.

preprint2014arXiv

Binary white dwarfs in the halo of the Milky Way

Aims: We study single and binary white dwarfs in the inner halo of the Milky Way in order to learn more about the conditions under which the population of halo stars was born, such as the initial mass function (IMF), the star formation history, or the binary fraction. Methods: We simulate the evolution of low-metallicity halo stars at distances up to ~ 3 kpc using the binary population synthesis code SeBa. We use two different white dwarf cooling models to predict the present-day luminosities of halo white dwarfs. We determine the white dwarf luminosity functions (WDLFs) for eight different halo models and compare these with the observed halo WDLF of white dwarfs in the SuperCOSMOS Sky Survey. Furthermore, we predict the properties of binary white dwarfs in the halo and determine the number of halo white dwarfs that is expected to be observed with the Gaia satellite. Results: By comparing the WDLFs, we find that a standard IMF matches the observations more accurately than a top-heavy one, but the difference with a bottom-heavy IMF is small. A burst of star formation 13 Gyr ago fits slightly better than a star formation burst 10 Gyr ago and also slightly better than continuous star formation $10-13$ Gyrs ago. Gaia will be the first instrument to constrain the bright end of the field halo WDLF, where contributions from binary WDs are considerable. Many of these will have He cores, of which a handful have atypical surface gravities ($\log g < 6$) and reach luminosities $\log(L/L_\odot) > 0$ in our standard model for WD cooling. These so called pre-WDs, if observed, can help us to constrain white dwarf cooling models and might teach us something about the fraction of halo stars that reside in binaries.

preprint2014arXiv

Complexity on dwarf galaxies scale: A bimodal distribution function in Sculptor

In previous work we have presented Schwarzschild models of the Sculptor dSph, demonstrating that this system could be embedded in dark matter halos that are either cusped or cored. Here we show that the non-parametric distribution function recovered through Schwarschild's method is bimodal in energy and angular momentum space for all best fitting mass models explored. We demonstrate that this bimodality is directly related to the two components known to be present in Sculptor through stellar populations analysis, although our method is purely dynamical in nature and does not use this prior information. It therefore constitutes independent confirmation of the existence of two physically distinct dynamical components in Sculptor and suggests a rather complex assembly history for this dwarf galaxy.

preprint2014arXiv

Dark influences: imprints of dark satellites on dwarf galaxies

In the context of the current $Λ$CDM cosmological model small dark matter haloes are abundant and satellites of dwarf galaxies are expected to be predominantly dark. Since low mass galaxies have smaller baryon fractions interactions with these satellites may leave particularly dramatic imprints. We uncover the influence of the most massive of these dark satellites on disky dwarf galaxies and the possible dynamical and morphological transformations that result from these interactions. We use a suite of carefully set-up, controlled simulations of isolated dwarf galaxies. The primary dwarf galaxies have solely a stellar disk in the dark matter halo and the secundaries are completely devoid of baryons. We vary the disk mass, halo concentration, initial disk thickness and inclination of the satellite orbit. The disky dwarf galaxies are heated and disrupted due to the minor merger event, more extremely for higher satellite over disk mass ratios, and the morphology and kinematics are significantly altered. Moreover, for less concentrated haloes the minor merger can completely destroy the disk leaving a low-luminosity spheroidal-like galaxy instead. We conclude that dwarf galaxies are very much susceptible to being disturbed by dark galaxies and that even a minor merger event can significantly disrupt and alter the structure and kinematics of a dwarf galaxy. This process may be seen as a new channel for the formation of dwarf spheroidal galaxies.

preprint2014arXiv

Elemental Abundances in Milky Way-like Galaxies from a Hierarchical Galaxy Formation Model

We develop a new method to account for the finite lifetimes of stars and trace individual abundances within a semi-analytic model of galaxy formation. At variance with previous methods, based on the storage of the (binned) past star formation history of model galaxies, our method projects the information about the metals produced by each simple stellar population (SSP) in the future. Using this approach, an accurate accounting of the timings and properties of the individual SSPs composing model galaxies is possible. We analyse the dependence of our chemical model on various ingredients, and apply it to six simulated haloes of roughly Milky Way mass and with no massive close neighbour at z=0. For all models considered, the [Fe/H] distributions of the stars in the disc component are in good agreement with Milky Way data, while for the spheroid component (whose formation we model only through mergers) these are offset low with respect to observational measurements for the Milky Way bulge. This is a consequence of narrow star formation histories, with relatively low rates of star formation. The slow recycling of gas and energy from supernovae in our chemical model has important consequences on the predicted star formation rates, which are systematically lower than the corresponding rates in the same physical model but with an instantaneous recycling approximation. The halo that resembles most our Galaxy in terms of its global properties also reproduces the observed relation between the average metallicity and luminosity of the Milky Way satellites, albeit with a slightly steeper slope.

preprint2014arXiv

Galactic kinematics and dynamics from RAVE stars

We analyse the kinematics of ~400000 RAVE stars. We split the sample into hot and cold dwarfs, red-clump and non-clump giants. The kinematics of the clump giants are consistent with being identical with those of non-clump giants. We fit Gaussian velocity ellipsoids to the meridional-plane components of velocity of each star class and give formulae from which the shape and orientation of the velocity ellipsoid can be determined at any location. The data are consistent with the giants and the cool dwarfs sharing the same velocity ellipsoids; sigma_z rises from 21 kms in the plane to sim 55 kms at |z|=2 kpc, while sigma_r rises from 37 kms to 82 kms. At (R,z) the longest axis of one of these velocity ellipsoids is inclined to the Galactic plane by an angle ~0.8 arctan(z/R). We use a novel formula to obtain precise fits to the highly non-Gaussian distributions of v_phi components. We compare the observed velocity distributions with the predictions of a dynamical model fitted to the velocities of stars that lie within ~150 pc of the Sun and star counts towards the Galactic pole. The model accurately reproduces the non-Gaussian nature of the v_r and v_z distributions and provides excellent fits to the data for v_z at all locations. The model v_phi distributions for the cool dwarfs fit the data extremely well, while those for the hot dwarfs have displacements to low v_phi that grow with |z| from very small values near the plane. At |z|>0.5 kpc, the theoretical v_phi distributions for giants show a deficit of stars with large v_phi and the model v_r distributions are too narrow. Systematically over-estimating distances by 20 per cent introduces asymmetry into the model v_r and v_z distributions near the plane and but significantly improves the fits to the data at |z|>0.5 kpc. The quality of the fits lends credence to the assumed, disc-dominated, gravitational potential.

preprint2014arXiv

More pieces of the puzzle: Chemistry and substructures in the Galactic thick disk

We present a study of the chemical abundances of Solar neighbourhood stars associated to dynamical structures in the Milky Way's (thick) disk. These stars were identified as overdensity in the eccentricity range 0.3< ecc < 0.5 in the Copenhagen-Geneva Survey by Helmi et al. (2006). We find that the stars with these dynamical characteristics do not constitute a homogeneous population. A relatively sharp transition in dynamical and chemical properties appears to occur at a metallicity of [Fe/H] ~ -0.4. Stars with [Fe/H] > -0.4 have mostly lower eccentricities, smaller vertical velocity dispersions, are alpha-enhanced and define a rather narrow sequence in [alpha/Fe] vs [Fe/H], clearly distinct from that of the thin disk. Stars with [Fe/H] < -0.4 have a range of eccentricities, are hotter vertically, and depict a larger spread in [alpha/Fe]. We have also found tentative evidence of substructure possibly associated to the disruption of a metal-rich star cluster. The differences between these populations of stars is also present in e.g. [Zn/Fe], [Ni/Fe] and [SmII/Fe], suggesting a real physical distinction.

preprint2014arXiv

New distances to RAVE stars

Probability density functions are determined from new stellar parameters for the distance moduli of stars for which the RAdial Velocity Experiment (RAVE) has obtained spectra with S/N>=10. Single-Gaussian fits to the pdf in distance modulus suffice for roughly half the stars, with most of the other half having satisfactory two-Gaussian representations. As expected, early-type stars rarely require more than one Gaussian. The expectation value of distance is larger than the distance implied by the expectation of distance modulus; the latter is itself larger than the distance implied by the expectation value of the parallax. Our parallaxes of Hipparcos stars agree well with the values measured by Hipparcos, so the expectation of parallax is the most reliable distance indicator. The latter are improved by taking extinction into account. The effective temperature absolute-magnitude diagram of our stars is significantly improved when these pdfs are used to make the diagram. We use the method of kinematic corrections devised by Schoenrich, Binney & Asplund to check for systematic errors for general stars and confirm that the most reliable distance indicator is the expectation of parallax. For cool dwarfs and low-gravity giants <pi> tends to be larger than the true distance by up to 30 percent. The most satisfactory distances are for dwarfs hotter than 5500 K. We compare our distances to stars in 13 open clusters with cluster distances from the literature and find excellent agreement for the dwarfs and indications that we are over-estimating distances to giants, especially in young clusters.

preprint2014arXiv

Pseudo Three-dimensional Maps of the Diffuse Interstellar Band at 862 nm

The diffuse interstellar bands (DIBs) are absorption lines observed in visual and near infrared spectra of stars. Understanding their origin in the interstellar medium (ISM) is one of the oldest problems in astronomical spectroscopy, as DIBs ahave been known since 1922. In a completely new approach to understanding DIBs, we combined information from nearly 500,000 stellar spectra obtained by the massive spectroscopic survey RAVE (Radial Velocity Experiment) to produce the first pseudo three-dimensional map of the strength of the DIB at 8620~Ångstroms covering the nearest 3 kiloparsecs from the Sun, and show that it follows our independently constructed spatial distribution of extinction by interstellar dust along the Galactic plane. Despite having a similar distribution in the Galactic plane, the DIB 8620 carrier has a significantly larger vertical scale height than the dust. Even if one DIB may not represent the general DIB population, our observations outline the future direction of DIB research.

preprint2014arXiv

The Galactic bar and the large scale velocity gradients in the Galactic disk

Aims: We investigate whether the cylindrical (galactocentric) radial velocity gradient of ~ -3 km/s/kpc, directed radially from the Galactic center and recently observed in the stars of the solar neighborhood with the RAVE survey, can be explained by the resonant effects of the bar near the solar neighborhood. Methods: We compared the results of test particle simulations of the Milky Way with a potential that includes a rotating bar with observations from the RAVE survey. To this end we applied the RAVE selection function to the simulations and convolved these with the characteristic RAVE errors. We explored different "solar neighborhoods" in the simulations, as well as different bar models Results: We find that the bar induces a negative radial velocity gradient at every height from the Galactic plane, outside the outer Lindblad resonance and for angles from the long axis of the bar compatible with the current estimates. The selection function and errors do not wash away the gradient, but often make it steeper, especially near the Galactic plane, because this is where the RAVE survey is less radially extended. No gradient in the vertical velocity ispresent in our simulations, from which we may conclude that this cannot be induced by the bar.

preprint2014arXiv

The shape of dark matter subhalos in the Aquarius simulations

We analyze the Aquarius simulations to characterize the shape of dark matter halos with peak circular velocity in the range 8<Vmax<200 km/s, and perform a convergence study using the various Aquarius resolution levels. For the converged objects, we determine the principal axis (a<b<c) of the normalized inertia tensor as a function of radius. We find that the triaxiality of field halos is an increasing function of halo mass, so that the smallest halos in our sample are ~40-50% rounder than Milky Way-like objects at the radius where the circular velocity peaks, rmax. We find that the distribution of subhalo axis ratios is consistent with that of field halos of comparable Vmax. Inner and outer contours within each object are well aligned, with the major axis preferentially pointing in the radial direction for subhalos closest to the center of their host halo. We also analyze the dynamical structure of subhalos likely to host luminous satellites comparable to the classical dwarf spheroidals in the Local Group. These halos have axis ratios that increase with radius, and which are mildly triaxial with <b/a>~0.75 and <c/a>~0.60 at r~1 kpc. Their velocity ellipsoid become strongly tangentially biased in the outskirts as a consequence of tidal stripping.

preprint2013arXiv

3D test particle simulations of the Galactic disks. The kinematical effects of the bar

Aims: To study the imprints of a rotating bar on the kinematics of stars in the thin and thick disks throughout the Galaxy. Methods: We perform test particle numerical simulations of the thin and thick disks in a 3D Galactic potential that includes a halo, a bulge, thin and thick disks, and a Ferrers bar. We analyze the resulting velocity distributions of populations corresponding to both disks, for different positions in the Galaxy and for different structural parameters of the bar. Results: We find that the velocity distributions of the disks are affected by the bar, and that strong transient effects are present for approximately 10 bar rotations after this is introduced adiabatically. On long (more realistic) timescales, the effects of the bar are strong on the kinematics of thin disk stars, and weaker on those in the thick disk, but in any case significant. Furthermore, we find that it is possible to trace the imprints of the bar also vertically and at least up to z ~ 1 kpc for the thin disk and z ~ 2 kpc for the thick disk.

preprint2013arXiv

An analytic distribution function for a massless cored stellar system in a cuspy dark matter halo

We demonstrate the existence of distribution functions that can be used to represent spherical massless cored stellar systems embedded in cuspy dark matter halos with constant mildly tangential velocity anisotropy. In particular, we derive analytically the functional form of the distribution function for a Plummer stellar sphere in a Hernquist dark halo, for β_0 = -0.5 and for different degrees of embedding. This particular example satisfies the condition that the central logarithmic slope of the light profile γ_0 > 2 β_0. Our models have velocity dispersion profiles similar to those observed in nearby dwarf spheroidal galaxies. Hence they can be used to generate initial conditions for a variety of problems, including N-body simulations that may represent dwarf galaxies in the Local Group.

preprint2013arXiv

An analytical phase-space model for tidal caustics

The class of tidal features around galaxies known variously as "shells" or "umbrellas" comprises debris that has arisen from high-mass-ratio mergers with low impact parameter; the nearly radial orbits of the debris give rise to a unique morphology, a universal density profile, and a tight correlation between positions and velocities of the material. As such they are accessible to analytical treatment, and can provide a relatively clean system for probing the gravitational potential of the host galaxy. In this work we present a simple analytical model that describes the density profile, phase-space distribution, and geometry of a shell, and whose parameters are directly related to physical characteristics of the interacting galaxies. The model makes three assumptions: that their orbit is radial, that the potential of the host is spherical at the shell radii, and that the satellite galaxy had a Maxwellian velocity distribution. We quantify the error introduced by the first two assumptions and show that selecting shells by their appearance on the sky is a sufficient basis to assume that these simplifications are valid. We further demonstrate that (1) given only an image of a shell, the radial gravitational force at the shell edge and the phase-space density of the satellite are jointly constrained, (2) that combining the image with measurements of either point line-of-sight velocities or integrated spectra will yield an independent estimate of the gravitational force at a shell, and (3) that an independent measurement of this force is obtained for each shell observed around a given galaxy, potentially enabling a determination of the galactic mass distribution.

preprint2013arXiv

Constraints on the shape of the Milky Way dark matter halo from the Sagittarius stream

We propose a new model for the dark matter halo of the Milky Way that fits the properties of the stellar stream associated with the Sagittarius dwarf galaxy. Our dark halo is oblate with q_z = 0.9 for r < 10 kpc, and can be made to follow the Law & Majewski model at larger radii. However, we find that the dynamical perturbations induced by the Large Magellanic Cloud on the orbit of Sgr cannot be neglected when modeling its streams. When taken into account, this leads us to constrain the Galaxy's outer halo shape to have minor-to-major axis ratio (c/a)_Φ= 0.8 and intermediate-to-major axis ratio (b/a)_Φ= 0.9, in good agreement with cosmological expectations.

preprint2013arXiv

Internal kinematics and dynamical models of dwarf spheroidal galaxies around the Milky Way

We review our current understanding of the internal dynamical properties of the dwarf spheroidal galaxies surrounding the Milky Way. These are the most dark matter dominated galaxies, and as such may be considered ideal laboratories to test the current concordance cosmological model, and in particular provide constraints on the nature of the dominant form of dark matter. We discuss the latest observations of the kinematics of stars in these systems, and how these may be used to derive their mass distribution. We tour through the various dynamical techniques used, with emphasis on the complementarity and limitations, and discuss what the results imply also in the context of cosmological models. Finally we provide an outlook on exciting developments in this field.

preprint2013arXiv

Streams in the Aquarius stellar haloes

We use the very high resolution, fully cosmological simulations from the Aquarius project, coupled to a semi-analytical model of galaxy formation, to study the phase-space distribution of halo stars in "solar neighbourhood"-like volumes. We find that this distribution is very rich in substructure in the form of stellar streams for all five stellar haloes we have analysed. These streams can be easily identified in velocity space, as well as in spaces of pseudo-conserved quantities such as E vs. Lz. In our best-resolved local volumes, the number of identified streams ranges from ~ 300 to 600, in very good agreement with previous analytical predictions, even in the presence of chaotic mixing. The fraction of particles linked to (massive) stellar streams in these volumes can be as large as 84%. The number of identified streams is found to decrease as a power-law with galactocentric radius. We show that the strongest limitation to the quantification of substructure in our poorest-resolved local volumes is particle resolution rather than strong diffusion due to chaotic mixing.

preprint2013arXiv

The satellites of the Milky Way - Insights from semi-analytic modelling in a LambdaCDM cosmology

We combine the six high-resolution Aquarius dark matter simulations with a semi-analytic galaxy formation model to investigate the properties of the satellites of Milky Way-like galaxies. We find good correspondence with the observed luminosity function, luminosity-metallicity relation and radial distribution of the Milky Way satellites. The star formation histories of the dwarf galaxies in our model vary widely, in accordance with what is seen observationally. Ram-pressure stripping of hot gas from the satellites leaves a clear imprint of the environment on the characteristics of a dwarf galaxy. We find that the fraction of satellites dominated by old populations of stars matches observations well. However, the internal metallicity distributions of the model satellites appear to be narrower than observed. This may indicate limitations in our treatment of chemical enrichment, which is based on the instantaneous recycling approximation. Our model works best if the dark matter halo of the Milky Way has a mass of ~8 x 10^11 Msun, in agreement with the lower estimates from observations. The galaxy that resembles the Milky Way the most also has the best matching satellite luminosity function, although it does not contain an object as bright as the SMC or LMC. Compared to other semi-analytic models and abundance matching relations we find that central galaxies reside in less massive haloes, but the halo mass-stellar mass relation in our model is consistent both with hydrodynamical simulations and with recent observations.

preprint2012arXiv

4MOST - 4-metre Multi-Object Spectroscopic Telescope

The 4MOST consortium is currently halfway through a Conceptual Design study for ESO with the aim to develop a wide-field (>3 square degree, goal >5 square degree), high-multiplex (>1500 fibres, goal 3000 fibres) spectroscopic survey facility for an ESO 4m-class telescope (VISTA). 4MOST will run permanently on the telescope to perform a 5 year public survey yielding more than 20 million spectra at resolution R~5000 (λ=390-1000 nm) and more than 2 million spectra at R~20,000 (395-456.5 nm & 587-673 nm). The 4MOST design is especially intended to complement three key all-sky, space-based observatories of prime European interest: Gaia, eROSITA and Euclid. Initial design and performance estimates for the wide-field corrector concepts are presented. We consider two fibre positioner concepts, a well-known Phi-Theta system and a new R-Theta concept with a large patrol area. The spectrographs are fixed configuration two-arm spectrographs, with dedicated spectrographs for the high- and low-resolution. A full facility simulator is being developed to guide trade-off decisions regarding the optimal field-of-view, number of fibres needed, and the relative fraction of high-to-low resolution fibres. Mock catalogues with template spectra from seven Design Reference Surveys are simulated to verify the science requirements of 4MOST. The 4MOST consortium aims to deliver the full 4MOST facility by the end of 2018 and start delivering high-level data products for both consortium and ESO community targets a year later with yearly increments.

preprint2012arXiv

Dark satellites and the morphology of dwarf galaxies

One of the strongest predictions of the LambdaCDM cosmological model is the presence of dark satellites orbiting all types of galaxies. We focus here on the dynamical effects of such satellites on disky dwarf galaxies, and demonstrate that these encounters can be dramatic. Although mergers with M_sat > M_d are not very common, because of the lower baryonic content they occur much more frequently on the dwarf scale than for L_*-galaxies. As an example, we present a numerical simulation of a 20% (virial) mass ratio merger between a dark satellite and a disky dwarf (akin to the Fornax dwarf galaxy in luminosity) that shows that the merger remnant has a spheroidal morphology. We conclude that perturbations by dark satellites provide a plausible path for the formation of dSph systems and also could trigger starbursts in gas rich dwarf galaxies. Therefore the transition from disky to the often amorphous, irregular, or spheroidal morphologies of dwarfs could be a natural consequence of the dynamical heating of hitherto unobservable dark satellites.

preprint2012arXiv

Not too big, not too small: the dark halos of the dwarf spheroidals in the Milky Way

We present a new analysis of the Aquarius simulations done in combination with a semi-analytic galaxy formation model. Our goal is to establish whether the subhalos present in LCDM simulations of Milky Way-like systems could host the dwarf spheroidal (dSph) satellites of our Galaxy. Our analysis shows that, contrary to what has been assumed in most previous work, the mass profiles of subhalos are generally not well fit by NFW models but that Einasto profiles are preferred. We find that for shape parameters alpha = 0.2 - 0.5 and Vmax = 10 - 30 km/s there is very good correspondence with the observational constraints obtained for the nine brightest dSph of the Milky Way. However, to explain the internal dynamics of these systems as well as the number of objects of a given circular velocity the total mass of the Milky Way should be ~ 8x10^11 Msun, a value that is in agreement with many recent determinations, and at the low mass end of the range explored by the Aquarius simulations. Our simulations show important scatter in the number of bright satellites, even when the Aquarius Milky Way-like hosts are scaled to a common mass, and we find no evidence for a missing population of massive subhalos in the Galaxy. This conclusion is also supported when we examine the dynamics of the satellites of M31.

preprint2012arXiv

Nucleosynthesis and the Inhomogeneous Chemical Evolution of the Carina Dwarf Galaxy

The detailed abundances of 23 elements in nine bright RGB stars in the Carina dSph are presented based on high resolution spectra gathered at the VLT and Magellan telescopes. A spherical model atmospheres analysis is applied using standard methods to spectra ranging from 380 to 680 nm. The stars in this analysis range from -2.9 < [Fe/H] < -1.3, and adopting the ages determined by Lemasle et al. (2012), we are able to examine the chemical evolution of Carina's old and intermediate-aged populations. One of the main results from this work is the evidence for inhomogeneous mixing in Carina; a large dispersion in [Mg/Fe] indicates poor mixing in the old population, an offset in the [alpha/Fe] ratios between the old and intermediate-aged populations (when examined with previously published results) suggests that the second star formation event occurred in alpha-enriched gas, and one star, Car-612, seems to have formed in a pocket enhanced in SN Ia/II products. This latter star provides the first direct link between the formation of stars with enhanced SN Ia/II ratios in dwarf galaxies to those found in the outer Galactic halo (Ivans et al. 2003). Another important result is the potential evidence for SN II driven winds. We show that the very metal-poor stars in Carina have not been enhanced in AGB or SN Ia products, and therefore their very low ratios of [Sr/Ba] suggests the loss of contributions from the early SNe II. Low ratios of [Na/Fe], [Mn/Fe], and [Cr/Fe] in two of these stars support this scenario, with additional evidence from the low [Zn/Fe] upper limit for one star. It is interesting that the chemistry of the metal-poor stars in Carina is not similar to those in the Galaxy, most of the other dSphs, or the UFDs, and suggests that Carina may be at the critical mass where some chemical enrichment events are lost through SN II driven winds.

preprint2011arXiv

Clues to the "Magellanic Galaxy" from Cosmological Simulations

We use cosmological simulations from the Aquarius Project to study the orbital history of the Large Magellanic Cloud (LMC) and its potential association with other satellites of the Milky Way (MW). We search for dynamical analogs to the LMC and find a subhalo that matches the LMC position and velocity at either of its two most recent pericentric passages. This suggests that the LMC is not necessarily on its first approach to the MW, provided that the virial mass of the Milky Way is M_200 ~1.8e12 M_sun. The simulation results yield specific predictions for the position and velocity of systems associated with the LMC prior to infall. If on first approach, most should lie close to the LMC because the Galactic tidal field has not yet had enough time to disperse them. If on second approach, the list of potential associates increases substantially. Interestingly, our analysis rules out an LMC association for Draco and Ursa Minor, two of the dwarf spheroidals suggested by Lynden-Bell & Lynden-Bell to form part of the "Magellanic Ghostly Stream". Our results also indicate that the direction of the orbital angular momentum is a powerful test of LMC association. This test, however, requires precise proper motions, which are unavailable for most MW satellites. Of the 4 satellites with published proper motions, only the Small Magellanic Cloud is clearly associated with the LMC. Taken at face value, the proper motions of Carina, Fornax and Sculptor rule them out as potential associates, but this conclusion should be revisited when better data become available. The dearth of satellites clearly associated with the Clouds might be solved by wide-field imaging surveys that target its surroundings, a region that may prove a fertile hunting ground for faint, previously unnoticed MW satellites.

preprint2011arXiv

Constraining the Milky Way potential using the dynamical kinematic substructures

We present a method to constrain the potential of the non-axisymmetric components of the Galaxy using the kinematics of stars in the solar neighborhood. The basic premise is that dynamical substructures in phase-space (i.e. due to the bar and/or spiral arms) are associated with families of periodic or irregular orbits, which may be easily identified in orbital frequency space. We use the "observed" positions and velocities of stars as initial conditions for orbital integrations in a variety of gravitational potentials. We then compute their characteristic frequencies, and study the structure present in the frequency maps. We find that the distribution of dynamical substructures in velocity- and frequency-space is best preserved when the integrations are performed in the "true" gravitational potential.

preprint2011arXiv

Observational Properties of the Metal-Poor Thick Disk of the Milky Way Galaxy and Insights into Its Origins

We have undertaken the study of the elemental abundances and kinematic properties of a metal-poor sample of candidate thick-disk stars selected from the RAVE spectroscopic survey of bright stars to differentiate among the present scenarios of the formation of the thick disk. In this paper, we report on a sample of 214 red giant branch, 31 red clump/horizontal branch, and 74 main-sequence/sub-giant branch metal-poor stars, which serves to augment our previous sample of only giant stars. We find that the thick disk [alpha/Fe] ratios are enhanced, and have little variation (<0.1 dex), in agreement with our previous study. The augmented sample further allows, for the first time, investigation of the gradients in the metal-poor thick disk. For stars with [Fe/H] < -1.2, the thick disk shows very small gradients, <0.03 +/- 0.02 dex/kpc, in alpha-enhancement, while we find a +0.01 +/- 0.04 dex/kpc radial gradient and a -0.09 +/- 0.05 dex/kpc vertical gradient in iron abundance. In addition, we show that the peak of the distribution of orbital eccentricities for our sample agrees better with models in which the stars that comprise the thick disk were formed primarily in the Galaxy, with direct accretion of stars contributing little. Our results thus disfavor direct accretion of stars from dwarf galaxies into the thick disk as a major contributor to the thick disk population, but cannot discriminate between alternative models for the thick disk, such as those that invoke high-redshift (gas-rich) mergers, heating of a pre-existing thin stellar disk by a minor merger, or efficient radial migration of stars.

preprint2011arXiv

Schwarzschild models of the Sculptor dSph galaxy

We have developed a spherically symmetric dynamical model of a dwarf spheroidal galaxy using the Schwarzschild method. This type of modelling yields constraints both on the total mass distribution (e.g. enclosed mass and scale radius) as well as on the orbital structure of the system modelled (e.g. velocity anisotropy). Therefore not only can we derive the dark matter content of these systems, but also explore possible formation scenarios. Here we present preliminary results for the Sculptor dSph. We find that the mass of Sculptor within 1kpc is 8.5\times10^(7\pm0.05) M\odot, its anisotropy profile is tangentially biased and slightly more isotropic near the center. For an NFW profile, the preferred concentration (~15) is compatible with cosmological models. Very cuspy density profiles (steeper than NFW) are strongly disfavoured for Sculptor.

preprint2011arXiv

Substructure in the stellar halos of the Aquarius simulations

We characterize substructure in the simulated stellar halos of Cooper et al. (2010) which were formed by the disruption of satellite galaxies within the cosmological N-body simulations of galactic halos of the Aquarius Project. These stellar halos exhibit a wealth of tidal features: broad overdensities and very narrow faint streams akin to those observed around the Milky Way. The substructures are distributed anisotropically on the sky, a characteristic that should become apparent in the next generation of photometric surveys. The normalized RMS of the density of stars on the sky appears to be systematically larger for our halos compared to the value estimated for the Milky Way from main sequence turn-off stars in the Sloan Digital Sky Survey. We show that this is likely to be due in part to contamination by faint QSOs and redder main sequence stars, and might suggest that ~10% of the Milky Way halo stars have formed in-situ.

preprint2011arXiv

The population of Milky Way satellites in the LambdaCDM cosmology

We present a model for the satellites of the Milky Way in which galaxy formation is followed using semi-analytic techniques applied to the six high-resolution N-body simulations of galactic halos of the Aquarius project. The model, calculated using the Galform code, incorporates improved treatments of the relevant physics in the LambdaCDM cosmogony, particularly a self-consistent calculation of reionization by UV photons emitted by the forming galaxy population, including the progenitors of the central galaxy. Along the merger tree of each halo, the model calculates gas cooling (by Compton scattering off cosmic microwave background photons, molecular hydrogen and atomic processes), gas heating (from hydrogen photoionization and supernova energy), star formation and evolution. The evolution of the intergalactic medium is followed simultaneously with that of the galaxies. Star formation in the more massive progenitor subhalos is suppressed primarily by supernova feedback, while for smaller subhalos it is suppressed primarily by photoionization due to external and internal sources. The model is constrained to match a wide range of properties of the present day galaxy population as a whole, but at high redshift it requires an escape fraction of UV photons near unity in order completely to reionize the universe by redshift z ~ 8. In the most successful model the local sources photoionize the pre-galactic region completely by z ~ 10. In addition to the luminosity function of Milky Way satellites, the model matches their observed luminosity-metallicity relation, their radial distribution and the inferred values of the mass within 300 pc, which in the models increase slowly but significantly with luminosity. There is a large variation in satellite properties from halo to halo, with the luminosity function, for example, varying by a factor of ~ 2 among the six simulations.

preprint2011arXiv

The Shape of Dark Matter Haloes in the Aquarius Simulations: Evolution and Memory

We use the high resolution cosmological N-body simulations from the Aquarius project to investigate in detail the mechanisms that determine the shape of Milky Way-type dark matter haloes. We find that, when measured at the instantaneous virial radius, the shape of individual haloes changes with time, evolving from a typically prolate configuration at early stages to a more triaxial/oblate geometry at the present day. This evolution in halo shape correlates well with the distribution of the infalling material: prolate configurations arise when haloes are fed through narrow filaments, which characterizes the early epochs of halo assembly, whereas triaxial/oblate configurations result as the accretion turns more isotropic at later times. Interestingly, at redshift z=0, clear imprints of the past history of each halo are recorded in their shapes at different radii, which also exhibit a variation from prolate in the inner regions to triaxial/oblate in the outskirts. Provided that the Aquarius haloes are fair representatives of Milky Way-like 10^12 Msun objects, we conclude that the shape of such dark matter haloes is a complex, time-dependent property, with each radial shell retaining memory of the conditions at the time of collapse.

preprint2010arXiv

Design drivers for a wide-field multi-object spectrograph for the William Herschel Telescope

Wide-field multi-object spectroscopy is a high priority for European astronomy over the next decade. Most 8-10m telescopes have a small field of view, making 4-m class telescopes a particularly attractive option for wide-field instruments. We present a science case and design drivers for a wide-field multi-object spectrograph (MOS) with integral field units for the 4.2-m William Herschel Telescope (WHT) on La Palma. The instrument intends to take advantage of a future prime-focus corrector and atmospheric-dispersion corrector that will deliver a field of view 2 deg in diameter, with good throughput from 370 to 1,000 nm. The science programs cluster into three groups needing three different resolving powers R: (1) high-precision radial-velocities for Gaia-related Milky Way dynamics, cosmological redshift surveys, and galaxy evolution studies (R = 5,000), (2) galaxy disk velocity dispersions (R = 10,000) and (3) high-precision stellar element abundances for Milky Way archaeology (R = 20,000). The multiplex requirements of the different science cases range from a few hundred to a few thousand, and a range of fibre-positioner technologies are considered. Several options for the spectrograph are discussed, building in part on published design studies for E-ELT spectrographs. Indeed, a WHT MOS will not only efficiently deliver data for exploitation of important imaging surveys planned for the coming decade, but will also serve as a test-bed to optimize the design of MOS instruments for the future E-ELT.

preprint2010arXiv

On the identification of merger debris in the {\it Gaia} Era

We model the formation of the Galactic stellar halo via the accretion of satellite galaxies onto a time-dependent semi-cosmological galactic potential. Our goal is to characterize the substructure left by these accretion events in a close manner to what may be possible with the {\it Gaia} mission. We have created a synthetic {\it Gaia} Solar Neighbourhood catalogue by convolving the 6D phase-space coordinates of stellar particles from our disrupted satellites with the latest estimates of the {\it Gaia} measurement errors, and included realistic background contamination due to the Galactic disc(s) and bulge. We find that, even after accounting for the expected observational errors, the resulting phase-space is full of substructure. We are able to successfully isolate roughly 50% of the different satellites contributing to the `Solar Neighbourhood' by applying the Mean-Shift clustering algorithm in energy-angular momentum space. Furthermore, a Fourier analysis of the space of orbital frequencies allows us to obtain accurate estimates of time since accretion for approximately 30% of the recovered satellites.

preprint2010arXiv

On the nature of the Milky Way satellites

We combine a series of high-resolution simulations with semi-analytic galaxy formation models to follow the evolution of a system resembling the Milky Way and its satellites. The semi-analytic model is based on that developed for the Millennium Simulation, and successfully reproduces the properties of galaxies on large scales, as well as those of the Milky Way. In this model, we are able to reproduce the luminosity function of the satellites around the Milky Way by preventing cooling in haloes with Vvir < 16.7 km/s (i.e. the atomic hydrogen cooling limit) and including the impact of the reionization of the Universe. The physical properties of our model satellites (e.g. mean metallicities, ages, half-light radii and mass-to-light ratios) are in good agreement with the latest observational measurements. We do not find a strong dependence upon the particular implementation of supernova feedback, but a scheme which is more efficient in galaxies embedded in smaller haloes, i.e. shallower potential wells, gives better agreement with the properties of the ultra-faint satellites. Our model predicts that the brightest satellites are associated with the most massive subhaloes, are accreted later (z $\lta$ 1), and have extended star formation histories, with only 1 per cent of their stars made by the end of the reionization. On the other hand, the faintest satellites were accreted early, are dominated by stars with age > 10 Gyr, and a few of them formed most of their stars before the reionization was complete. Objects with luminosities comparable to those of the classical MW satellites are associated with dark matter subhaloes with a peak circular velocity $\gta$ 10 km/s, in agreement with the latest constraints.

preprint2010arXiv

Testing formation mechanisms of the Milky Way's thick disc with RAVE

We study the eccentricity distribution of a thick disc sample of stars observed in the Radial Velocity Experiment (RAVE) and compare it to that expected in four simulations of thick disc formation in the literature (accretion of satellites, heating of a primordial thin disc during a merger, radial migration, and gas-rich mergers), as compiled by Sales et al. (2009). We find that the distribution of our sample is peaked at low eccentricities and falls off smoothly and rather steeply to high eccentricities. This distribution is fairly robust to changes in distances, thin disc contamination, and the particular thick disc sample used. Our results are inconsistent with what is expected for the pure accretion simulation, since we find that the dynamics of local thick disc stars implies that the majority must have formed "in situ". Of the remaining models explored, the eccentricity distribution of our stars appears to be most consistent with the gas-rich merger case.

preprint2010arXiv

The Dawning of the Stream of Aquarius in RAVE

We identify a new, nearby (0.5 < d < 10 kpc) stream in data from the RAdial Velocity Experiment (RAVE). As the majority of stars in the stream lie in the constellation of Aquarius we name it the Aquarius Stream. We identify 15 members of the stream lying between 30 < l < 75 and -70< b <-50, with heliocentric line-of-sight velocities V_los~-200 km/s. The members are outliers in the radial velocity distribution, and the overdensity is statistically significant when compared to mock samples created with both the Besançon Galaxy model and newly-developed code Galaxia. The metallicity distribution function and isochrone fit in the log g - T_eff plane suggest the stream consists of a 10 Gyr old population with [m/H]~-1.0. We explore relations to other streams and substructures, finding the stream cannot be identified with known structures: it is a new, nearby substructure in the Galaxy's halo. Using a simple dynamical model of a dissolving satellite galaxy we account for the localization of the stream. We find that the stream is dynamically young and therefore likely the debris of a recently disrupted dwarf galaxy or globular cluster. The Aquarius stream is thus a specimen of ongoing hierarchical Galaxy formation, rare for being right in the solar suburb.

preprint2010arXiv

The NIR Ca II triplet at low metallicity - Searching for extremely low-metallicity stars in classical dwarf galaxies

The NIR Ca II triplet absorption lines have proven to be an important tool for quantitative spectroscopy of individual red giant branch stars in the Local Group, providing a better understanding of metallicities of stars in the Milky Way and dwarf galaxies and thereby an opportunity to constrain their chemical evolution processes. An interesting puzzle in this field is the significant lack of extremely metal-poor stars, below [Fe/H]=-3, found in classical dwarf galaxies around the Milky Way using this technique. The question arises whether these stars are really absent, or if the empirical Ca II triplet method used to study these systems is biased in the low-metallicity regime. Here we present results of synthetic spectral analysis of the Ca II triplet, that is focused on a better understanding of spectroscopic measurements of low-metallicity giant stars. Our results start to deviate strongly from the widely-used and linear empirical calibrations at [Fe/H]<-2. We provide a new calibration for Ca II triplet studies which is valid for -0.5<[Fe/H]<-4. We subsequently apply this new calibration to current data sets and suggest that the classical dwarf galaxies are not so devoid of extremely low-metallicity stars as was previously thought.

preprint2010arXiv

Thick-disk evolution induced by the growth of an embedded thin disk

We perform collisionless N-body simulations to investigate the evolution of the structural and kinematical properties of simulated thick disks induced by the growth of an embedded thin disk. The thick disks used in the present study originate from cosmologically-common 5:1 encounters between initially-thin primary disk galaxies and infalling satellites. The growing thin disks are modeled as static gravitational potentials and we explore a variety of growing-disk parameters that are likely to influence the response of thick disks. We find that the final thick-disk properties depend strongly on the total mass and radial scale-length of the growing thin disk, and much less sensitively on its growth timescale and vertical scale-height as well as the initial sense of thick-disk rotation. Overall, the growth of an embedded thin disk can cause a substantial contraction in both the radial and vertical direction, resulting in a significant decrease in the scale-lengths and scale-heights of thick disks. Kinematically, a growing thin disk can induce a notable increase in the mean rotation and velocity dispersions of thick-disk stars. We conclude that the reformation of a thin disk via gas accretion may play a significant role in setting the structure and kinematics of thick disks, and thus it is an important ingredient in models of thick-disk formation.

preprint2009arXiv

Mapping the Galactic Halo VIII: Quantifying substructure

We have measured the amount of kinematic substructure in the Galactic halo using the final data set from the Spaghetti project, a pencil-beam high latitude sky survey. Our sample contains 101 photometrically selected and spectroscopically confirmed giants with accurate distance, radial velocity and metallicity information. We have developed a new clustering estimator: the "4distance" measure, which when applied to our data set leads to the identification of 1 group and 7 pairs of clumped stars. The group, with 6 members, can confidently be matched to tidal debris of the Sagittarius dwarf galaxy. Two pairs match the properties of known Virgo structures. Using models of the disruption of Sagittarius in Galactic potentials with different degrees of dark halo flattening, we show that this favors a spherical or prolate halo shape, as demonstrated by Newberg et al. (2007) using SDSS data. One additional pair can be linked to older Sagittarius debris. We find that 20% of the stars in the Spaghetti data set are in substructures. From comparison with random data sets we derive a very conservative lower limit of 10% to the amount of substructure in the halo. However, comparison to numerical simulations shows that our results are also consistent with a halo entirely built up from disrupted satellites, provided the dominating features are relatively broad due to early merging or relatively heavy progenitor satellites.

preprint2009arXiv

On the identification of substructure in phase-space using orbital frequencies

We study the evolution of satellite debris to establish the most suitable space to identify past merger events. We confirm that the space of orbital frequencies is very promising in this respect. In frequency space individual streams can be easily identified, and their separation provides a direct measurement of the time of accretion. We are able to show for a few idealised gravitational potentials that these features are preserved also in systems that have evolved strongly in time. Furthermore, this time evolution is imprinted in the distribution of streams in frequency space. We have also tested the power of the orbital frequencies in a fully self-consistent (live) N-body simulation of the merger between a disk galaxy and a massive satellite. Even in this case streams can be easily identified and the time of accretion of the satellite can be accurately estimated.

preprint2009arXiv

Orbital Eccentricity as a probe of Thick Disk Formation Scenarios

We study the orbital properties of stars in four (published) simulations of thick disks formed by: i) accretion from disrupted satellites, ii) heating of a pre-existing thin disk by a minor merger, iii) radial migration and iv) gas rich mergers. We find that the distribution of orbital eccentricities are predicted to be different for each model: a prominent peak at low eccentricity is expected for the heating, migration and gas-rich merging scenarios, while the eccentricity distribution is broader and shifted towards higher values for the accretion model. These differences can be traced back to whether the bulk of the stars in each case is formed 'in-situ' or is 'accreted', and are robust to the peculiarities of each model. A simple test based on the eccentricity distribution of nearby thick disk stars may thus help elucidate the dominant formation mechanism of the Galactic thick disk.

preprint2009arXiv

Simulations of minor mergers. II. The phase-space structure of thick discs

We analyse the phase-space structure of simulated thick discs that are the result of a significant merger between a disc galaxy and a satellite. Our main goal is to establish what would be the characteristic imprints of a merger origin for the Galactic thick disc. We find that the spatial distribution predicted for thick disc stars is asymmetric, seemingly in agreement with recent observations of the Milky Way thick disc. Near the Sun, the accreted stars are expected to rotate more slowly, to have broad velocity distributions, and to occupy preferentially the wings of the line-of-sight velocity distributions. The majority of the stars in our model thick discs have low eccentricity orbits (in clear reference to the pre-existing heated disc) which gives rise to a characteristic (sinusoidal) pattern for their line of sight velocities as function of galactic longitude. The z-component of the angular momentum of thick disc stars provides a clear discriminant between stars from the pre-existing disc and those from the satellite, particularly at large radii. These results are robust against the particular choices of initial conditions made in our simulations, and thus provide clean tests of the disc heating via a minor merger scenario for the formation of thick discs.

preprint2009arXiv

The Diversity and Similarity of Simulated Cold Dark Matter Halos

We study the mass, velocity dispersion, and anisotropy profiles of $Λ$CDM halos using a suite of N-body simulations of unprecedented numerical resolution (the {\it Aquarius Project}). Our analysis confirms a number of results claimed by earlier work, and clarifies a few issues where conflicting claims may be found in the recent literature. The spherically-averaged density profile becomes progressively shallower inwards and, at the innermost resolved radius, the logarithmic slope is $γ\equiv -$d$\lnρ/$d$\ln r \simlt 1$. Asymptotic inner slopes as steep as the recently claimed $ρ\propto r^{-1.2}$ are clearly ruled out. The radial dependence of $γ$ is well approximated by a power-law, $γ\propto r^α$ (the Einasto profile). The shape parameter, $α$, varies slightly but significantly from halo to halo, implying that the mass profiles of $Λ$CDM halos are not strictly universal: different halos cannot, in general, be rescaled to look identical. Departures from similarity are also seen in velocity dispersion profiles and correlate with those in density profiles so as to preserve a power-law form for the spherically averaged pseudo-phase-space density, $ρ/σ^3\propto r^{-1.875}$. Our conclusions are reliable down to radii below 0.4% of the virial radius, providing well-defined predictions for halo structure when baryonic effects are neglected, and thus an instructive theoretical template against which the modifications induced by the baryonic components of real galaxies can be judged.

preprint2009arXiv

The Effect of Tidal Stripping on Composite Stellar Populations in Dwarf Spheroidal Galaxies

We use N-body simulations to study the effects of tides on the kinematical structure of satellite galaxies orbiting a Milky Way-like potential. Our work is motivated by observations of dwarf spheroidal galaxies in the Local Group, for which often a distinction is possible between a cold centrally concentrated metal rich and a hot, extended metal poor population. We find that an important attenuation of the initial differences in the distribution of the two stellar components occurs for orbits with small pericentric radii (r_per < 20 kpc). This is mainly due to: i) the loss of the gravitational support provided by the dark matter component after tidal stripping takes place, which forces a re-configuration of the luminous components, and ii) tides preferentially affect the more extended stellar component, leading to a net decrease in its velocity dispersion as a response for the mass loss, which thus shrinks the kinematical gap. We apply these ideas to the Sculptor and Carina dwarf spheroidals. Differences in their orbits might help to explain, under the assumption of similar initial configurations, why in the former a clear kinematical separation between metal poor and metal rich stars is apparent, while in Carina this segregation is significantly more subtle.

preprint2008arXiv

Group infall of substructures on to a Milky Way-like dark halo

We report the discovery that substructures/subhaloes of a galaxy-size halo tend to fall in together in groups in cosmological simulations, something that may explain the oddity of the MW satellite distribution. The original clustering at the time of infall is still discernible in the angular momenta of the subhaloes even for events which took place up to eight Gyrs ago, $z \sim 1$. This phenomenon appears to be rather common since at least 1/3 of the present-day subhaloes have fallen in groups in our simulations. Hence, this may well explain the Lynden-Bell & Lynden-Bell ghostly streams. We have also found that the probability of building up a flattened distribution similar to the MW satellites is as high as $\sim 80%$ if the MW satellites were from only one group and $\sim 20%$ when five groups are involved. Therefore, we conclude that the `peculiar' distribution of satellites around the MW can be expected with the CDM structure formation theory. This non-random assignment of satellites to subhaloes implies an environmental dependence on whether these low-mass objects are able to form stars, possibly related to the nature of reionization in the early Universe.

preprint2008arXiv

The Aquarius Project: the subhalos of galactic halos

We have performed the largest ever particle simulation of a Milky Way-sized dark matter halo, and present the most comprehensive convergence study for an individual dark matter halo carried out thus far. We have also simulated a sample of 6 ultra-highly resolved Milky-way sized halos, allowing us to estimate the halo-to-halo scatter in substructure statistics. In our largest simulation, we resolve nearly 300,000 gravitationally bound subhalos within the virialized region of the halo. Simulations of the same object differing in mass resolution by factors up to 1800 accurately reproduce the largest subhalos with the same mass, maximum circular velocity and position, and yield good convergence for the abundance and internal properties of dark matter substructures. We detect up to four generations of subhalos within subhalos, but contrary to recent claims, we find less substructure in subhalos than in the main halo when regions of equal mean overdensity are compared. The overall substructure mass fraction is much lower in subhalos than in the main halo. Extrapolating the main halo's subhalo mass spectrum down to an Earth mass, we predict the mass fraction in substructure to be well below 3% within 100 kpc, and to be below 0.1% within the Solar Circle. The inner density profiles of subhalos show no sign of converging to a fixed asymptotic slope and are well fit by gently curving profiles of Einasto form. The mean concentrations of isolated halos are accurately described by the fitting formula of Neto et al. down to maximum circular velocities of 1.5 km/s, an extrapolation over some 5 orders of magnitude in mass. However, at equal maximum circular velocity, subhalos are more concentrated than field halos, with a characteristic density that is typically ~2.6 times larger and increases towards the halo centre.

preprint2008arXiv

The Arcturus Moving Group: Its Place in the Galaxy

The Arcturus moving group is a well-populated example of phase space substructure within the thick disk of our Galaxy. Traditionally regarded as the remains of a dissolved open cluster, it has recently been suggested to be a remnant of a satellite accreted by our Galaxy. To investigate the origin of the group we undertook a high-resolution spectroscopic abundance study of Arcturus group members and candidates. Examining abundance of Fe, Mg, Ca, Ti, Cr, Ni, Zn, Ce, Nd, Sm and Gd for 134 stars we found that the group is chemically similar to disk stars and does not exhibit a clear chemical homogeneity. Furthermore, we confirm the existence of the Arcturus group using the Nordstroem et al. (2004), Schuster et al. (2006) and RAdial Velocity Experiment (RAVE) surveys (Steinmetz et al. 2006). The origin of the group still remains unresolved: the chemical results are consistent with a dynamical origin but do not entirely rule out a merger one. Certainly, the Arcturus group provides a challenge to our understanding of the nature and origin of the Galaxy's thick disk.

preprint2008arXiv

The Galaxy and its stellar halo - insights from a hybrid cosmological approach

We use a series of high-resolution N-body simulations of a `Milky-Way' halo, coupled to semi-analytic techniques, to study the formation of our own Galaxy and of its stellar halo. Our model Milky Way galaxy is a relatively young system whose physical properties are in quite good agreement with observational determinations. In our model, the stellar halo is mainly formed from a few massive satellites accreted early on during the galaxy's lifetime. The stars in the halo do not exhibit any metallicity gradient, but higher metallicity stars are more centrally concentrated than stars with lower abundances. This is due to the fact that the most massive satellites contributing to the stellar halo are also more metal rich, and dynamical friction drags them closer to the inner regions of the host halo.

preprint1999arXiv

Debris streams in the solar neighbourhood as relicts from the formation of the Milky Way

It is now generally believed that galaxies were built up through gravitational amplification of primordial fluctuations and the subsequent merging of smaller precursor structures. The stars of the structures that assembled to form the Milky Way should now make up much or all of its bulge and halo, in which case one hopes to find "fossil" evidence for those precursor structures in the present distribution of halo stars. Confirmation that this process is continuing came with the discovery of the Sagittarius dwarf galaxy, which is being disrupted by the Milky Way, but direct evidence that this process provided the bulk of the Milky Way's population of old stars has so far been lacking. Here we show that about ten per cent of the metal--poor stars in the halo of the Milky Way, outside the radius of the Sun's orbit, come from a single coherent structure that was disrupted during or soon after the Galaxy's formation. This object had a highly inclined orbit about the Milky Way at a maximum distance of $\sim$ 16 kpc, and it probably resembled the Fornax and Sagittarius dwarf spheroidal galaxies.

Amina Helmi

What is connected

Connect this record

See the researcher in context

Building this map preview

71 published item(s)

Chronology of our Galaxy from Gaia colour-magnitude diagram fitting (ChronoGal) -- III. Age and metallicity distribution of Gaia-Sausage-Enceladus stars near the Sun

Chronology of our Galaxy from Gaia colour-magnitude diagram fitting (ChronoGal) II. Unveiling the formation and evolution of the kinematically selected Thick and Thin Discs

Characterization and dynamics of the peculiar stream Jhelum, A tentative role for the Sagittarius dwarf galaxy

Substructure in the stellar halo near the Sun. I. Data-driven clustering in Integrals of Motion space

Substructure in the stellar halo near the Sun. II. Characterisation of independent structures

Age Determination of Galaxy Merger Remnant Stars using Asteroseismology

Galactic potential constraints from clustering in action space of combined stellar stream data

Substructures, Resonances and debris streams

Cosmological insights into the assembly of the radial and compact stellar halo of the Milky Way

Frequencies, chaos and resonances: a study of orbital parameters of nearby thick disc and halo stars

Streams, substructures and the early history of the Milky Way

The Sixth Data Release of the Radial Velocity Experiment (RAVE) -- I: Survey Description, Spectra and Radial Velocities

The Sixth Data Release of the Radial Velocity Experiment (RAVE) -- II: Stellar Atmospheric Parameters, Chemical Abundances and Distances

Stellar 3-D kinematics in the Draco dwarf spheroidal galaxy

Building Blocks of the Milky Way's Accreted Spheroid

Dark influences III. Structural characterization of minor mergers of dwarf galaxies with dark satellites

Stream-subhalo interactions in the Aquarius simulations

The kinematics of globular clusters systems in the outer halos of the Aquarius simulations

The time evolution of gaps in tidal streams

Action-space clustering of tidal streams to infer the Galactic potential

Dark influences II: gas and star formation in minor mergers of dwarf galaxies with dark satellites

Matching the dark matter profiles of dSph galaxies with those of simulated satellites: a two parameter comparison

The evolution of streams in a time-dependent potential

What can Gaia proper motions tell us about Milky Way dwarf galaxies?

A new fitting-function to describe the time evolution of a galaxy's gravitational potential

Binary white dwarfs in the halo of the Milky Way

Complexity on dwarf galaxies scale: A bimodal distribution function in Sculptor

Dark influences: imprints of dark satellites on dwarf galaxies

Elemental Abundances in Milky Way-like Galaxies from a Hierarchical Galaxy Formation Model

Galactic kinematics and dynamics from RAVE stars

More pieces of the puzzle: Chemistry and substructures in the Galactic thick disk

New distances to RAVE stars

Pseudo Three-dimensional Maps of the Diffuse Interstellar Band at 862 nm

The Galactic bar and the large scale velocity gradients in the Galactic disk

The shape of dark matter subhalos in the Aquarius simulations

3D test particle simulations of the Galactic disks. The kinematical effects of the bar

An analytic distribution function for a massless cored stellar system in a cuspy dark matter halo

An analytical phase-space model for tidal caustics

Constraints on the shape of the Milky Way dark matter halo from the Sagittarius stream

Internal kinematics and dynamical models of dwarf spheroidal galaxies around the Milky Way

Streams in the Aquarius stellar haloes

The satellites of the Milky Way - Insights from semi-analytic modelling in a LambdaCDM cosmology

4MOST - 4-metre Multi-Object Spectroscopic Telescope

Dark satellites and the morphology of dwarf galaxies

Not too big, not too small: the dark halos of the dwarf spheroidals in the Milky Way

Nucleosynthesis and the Inhomogeneous Chemical Evolution of the Carina Dwarf Galaxy

Clues to the "Magellanic Galaxy" from Cosmological Simulations

Constraining the Milky Way potential using the dynamical kinematic substructures

Observational Properties of the Metal-Poor Thick Disk of the Milky Way Galaxy and Insights into Its Origins

Schwarzschild models of the Sculptor dSph galaxy

Substructure in the stellar halos of the Aquarius simulations

The population of Milky Way satellites in the LambdaCDM cosmology

The Shape of Dark Matter Haloes in the Aquarius Simulations: Evolution and Memory

Design drivers for a wide-field multi-object spectrograph for the William Herschel Telescope

On the identification of merger debris in the {\it Gaia} Era

On the nature of the Milky Way satellites

Testing formation mechanisms of the Milky Way's thick disc with RAVE

The Dawning of the Stream of Aquarius in RAVE

The NIR Ca II triplet at low metallicity - Searching for extremely low-metallicity stars in classical dwarf galaxies

Thick-disk evolution induced by the growth of an embedded thin disk

Mapping the Galactic Halo VIII: Quantifying substructure

On the identification of substructure in phase-space using orbital frequencies

Orbital Eccentricity as a probe of Thick Disk Formation Scenarios

Simulations of minor mergers. II. The phase-space structure of thick discs

The Diversity and Similarity of Simulated Cold Dark Matter Halos

The Effect of Tidal Stripping on Composite Stellar Populations in Dwarf Spheroidal Galaxies

Group infall of substructures on to a Milky Way-like dark halo

The Aquarius Project: the subhalos of galactic halos

The Arcturus Moving Group: Its Place in the Galaxy

The Galaxy and its stellar halo - insights from a hybrid cosmological approach

Debris streams in the solar neighbourhood as relicts from the formation of the Milky Way