-
Axisymmetric Schwarzschild models of an isothermal axisymmetric mock dwarf spheroidal galaxy
Authors:
Jorrit H. J. Hagen,
Amina Helmi,
Maarten A. Breddels
Abstract:
We test the ability of Schwarzschild's orbit superposition method in measuring the mass content, scale radius and shape of a flattened dwarf spheroidal galaxy. Until now, most dynamical model efforts have assumed that dwarf spheroidal galaxies and their host halos are spherical. We use an Evans model (1993) to construct an isothermal mock galaxy whose properties somewhat resemble those of the Scul…
▽ More
We test the ability of Schwarzschild's orbit superposition method in measuring the mass content, scale radius and shape of a flattened dwarf spheroidal galaxy. Until now, most dynamical model efforts have assumed that dwarf spheroidal galaxies and their host halos are spherical. We use an Evans model (1993) to construct an isothermal mock galaxy whose properties somewhat resemble those of the Sculptor dwarf spheroidal galaxy. This mock galaxy contains flattened luminous and dark matter components, resulting in a logarithmic profile for the global potential. We have tested how well our Schwarzschild method could constrain the characteristic parameters of the system for different sample sizes, and also if the functional form of the potential was unknown. When assuming the true functional form of the potential, the Schwarzschild modelling technique is able to provide an accurate and precise measurement of the characteristic mass parameter of the system and reproduces well the light distribution and the stellar kinematics of our mock galaxy. When assuming a different functional form for the potential, such as a flattened NFW profile, we also constrain the mass and scale radius to their expected values. However in both cases, we find that the flattening parameter remains largely unconstrained. This is likely because the information content of the velocity dispersion on the geometric shape of the potential is too small, since $σ$ is constant across our mock dSph. Our results using Schwarzschild's method indicate that the mass enclosed can be derived reliably, even if the flattening parameter is unknown, and already for samples containing 2000 line-of-sight radial velocities, such as those currently available. Further applications of the method to more general distribution functions of flattened systems are needed to establish how well the flattening of dSph dark halos can be determined.
△ Less
Submitted 29 June, 2019;
originally announced July 2019.
-
Gaia Data Release 2: Observational Hertzsprung-Russell diagrams
Authors:
Gaia Collaboration,
C. Babusiaux,
F. van Leeuwen,
M. A. Barstow,
C. Jordi,
A. Vallenari,
D. Bossini,
A. Bressan,
T. Cantat-Gaudin,
M. van Leeuwen,
A. G. A. Brown,
T. Prusti,
J. H. J. de Bruijne,
C. A. L. Bailer-Jones,
M. Biermann,
D. W. Evans,
L. Eyer,
F. Jansen,
S. A. Klioner,
U. Lammers,
L. Lindegren,
X. Luri,
F. Mignard,
C. Panem,
D. Pourbaix
, et al. (428 additional authors not shown)
Abstract:
We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the select…
▽ More
We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.
△ Less
Submitted 13 August, 2018; v1 submitted 25 April, 2018;
originally announced April 2018.
-
Gaia Data Release 2: Catalogue validation
Authors:
F. Arenou,
X. Luri,
C. Babusiaux,
C. Fabricius,
A. Helmi,
T. Muraveva,
A. C. Robin,
F. Spoto,
A. Vallenari,
T. Antoja,
T. Cantat-Gaudin,
C. Jordi,
N. Leclerc,
C. Reylé,
M. Romero-Gómez,
I-C. Shih,
S. Soria,
C. Barache,
D. Bossini,
A. Bragaglia,
M. A. Breddels,
M. Fabrizio,
S. Lambert,
P. M. Marrese,
D. Massari
, et al. (14 additional authors not shown)
Abstract:
The second Gaia data release (DR2), contains very precise astrometric and photometric properties for more than one billion sources, astrophysical parameters for dozens of millions, radial velocities for millions, variability information for half a million of stellar sources and orbits for thousands of solar system objects. Before the Catalogue publication, these data have undergone dedicated valid…
▽ More
The second Gaia data release (DR2), contains very precise astrometric and photometric properties for more than one billion sources, astrophysical parameters for dozens of millions, radial velocities for millions, variability information for half a million of stellar sources and orbits for thousands of solar system objects. Before the Catalogue publication, these data have undergone dedicated validation processes. The goal of this paper is to describe the validation results in terms of completeness, accuracy and precision of the various Gaia DR2 data. The validation processes include a systematic analysis of the Catalogue content to detect anomalies, either individual errors or statistical properties, using statistical analysis, and comparisons to external data or to models. Although the astrometric, photometric and spectroscopic data are of unprecedented quality and quantity, it is shown that the data cannot be used without a dedicated attention to the limitations described here, in the Catalogue documentation and in accompanying papers. A particular emphasis is put on the caveats for the statistical use of the data in scientific exploitation.
△ Less
Submitted 25 April, 2018;
originally announced April 2018.
-
Vaex: Big Data exploration in the era of Gaia
Authors:
Maarten A. Breddels,
Jovan Veljanoski
Abstract:
We present a new Python library called vaex, to handle extremely large tabular datasets, such as astronomical catalogues like the Gaia catalogue, N-body simulations or any other regular datasets which can be structured in rows and columns. Fast computations of statistics on regular N-dimensional grids allows analysis and visualization in the order of a billion rows per second. We use streaming alg…
▽ More
We present a new Python library called vaex, to handle extremely large tabular datasets, such as astronomical catalogues like the Gaia catalogue, N-body simulations or any other regular datasets which can be structured in rows and columns. Fast computations of statistics on regular N-dimensional grids allows analysis and visualization in the order of a billion rows per second. We use streaming algorithms, memory mapped files and a zero memory copy policy to allow exploration of datasets larger than memory, e.g. out-of-core algorithms. Vaex allows arbitrary (mathematical) transformations using normal Python expressions and (a subset of) numpy functions which are lazily evaluated and computed when needed in small chunks, which avoids wasting of RAM. Boolean expressions (which are also lazily evaluated) can be used to explore subsets of the data, which we call selections. Vaex uses a similar DataFrame API as Pandas, a very popular library, which helps migration from Pandas. Visualization is one of the key points of vaex, and is done using binned statistics in 1d (e.g. histogram), in 2d (e.g. 2d histograms with colormapping) and 3d (using volume rendering). Vaex is split in in several packages: vaex-core for the computational part, vaex-viz for visualization mostly based on matplotlib, vaex-jupyter for visualization in the Jupyter notebook/lab based in IPyWidgets, vaex-server for the (optional) client-server communication, vaex-ui for the Qt based interface, vaex-hdf5 for hdf5 based memory mapped storage, vaex-astro for astronomy related selections, transformations and memory mapped (column based) fits storage. Vaex is open source and available under MIT license on github, documentation and other information can be found on the main website: https://vaex.io, https://docs.vaex.io or https://github.com/maartenbreddels/vaex
△ Less
Submitted 8 January, 2018;
originally announced January 2018.
-
3D motions in the Sculptor dwarf galaxy as a glimpse of a new era
Authors:
D. Massari,
M. A. Breddels,
A. Helmi,
L. Posti,
A. G. A Brown,
E. Tolstoy
Abstract:
The 3D motions of stars in small galaxies beyond our own are minute and yet they are crucial for our understanding of the nature of gravity and dark matter. Even for the dwarf galaxy Sculptor which is one of the best studied systems and inferred to be strongly dark matter dominated, there are conflicting reports on its mean motion around the Milky Way and the 3D internal motions of its stars have…
▽ More
The 3D motions of stars in small galaxies beyond our own are minute and yet they are crucial for our understanding of the nature of gravity and dark matter. Even for the dwarf galaxy Sculptor which is one of the best studied systems and inferred to be strongly dark matter dominated, there are conflicting reports on its mean motion around the Milky Way and the 3D internal motions of its stars have never been measured. Here we report, based on data from the Gaia space mission and the Hubble Space Telescope, a new precise measurement of Sculptor's mean proper motion. From this we deduce that Sculptor is currently at its closest approach to the Milky Way and moving on an elongated high-inclination orbit that takes it much farther away than previously thought. For the first time we are also able to measure the internal motions of stars in Sculptor. We find $σ_{R}=11.5 \pm 4.3$ km/s and $σ_{T}=8.5\pm3.2$ km/s along the projected radial and tangential directions, implying that the stars in our sample move preferentially on radial orbits as quantified by the anisotropy parameter, which we find to be $β\sim 0.86^{+0.12}_{-0.83}$ at a location beyond the core radius. Taken at face value such a high radial anisotropy requires abandoning conventional models for the mass distribution in Sculptor. Our sample is dominated by metal-rich stars and for these we find $β^{MR} \sim 0.95^{+0.04}_{-0.27}$, a value consistent with multi-component models where Sculptor is embedded in a cuspy dark halo as expected for cold dark matter.
△ Less
Submitted 24 November, 2017;
originally announced November 2017.
-
Gaia Data Release 1. Testing the parallaxes with local Cepheids and RR Lyrae stars
Authors:
Gaia Collaboration,
G. Clementini,
L. Eyer,
V. Ripepi,
M. Marconi,
T. Muraveva,
A. Garofalo,
L. M. Sarro,
M. Palmer,
X. Luri,
R. Molinaro,
L. Rimoldini,
L. Szabados,
I. Musella,
R. I. Anderson,
T. Prusti,
J. H. J. de Bruijne,
A. G. A. Brown,
A. Vallenari,
C. Babusiaux,
C. A. L. Bailer-Jones,
U. Bastian,
M. Biermann,
D. W. Evans,
F. Jansen
, et al. (566 additional authors not shown)
Abstract:
Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by…
▽ More
Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.
△ Less
Submitted 1 May, 2017;
originally announced May 2017.
-
Gaia Data Release 1. Open cluster astrometry: performance, limitations, and future prospects
Authors:
Gaia Collaboration,
F. van Leeuwen,
A. Vallenari,
C. Jordi,
L. Lindegren,
U. Bastian,
T. Prusti,
J. H. J. de Bruijne,
A. G. A. Brown,
C. Babusiaux,
C. A. L. Bailer-Jones,
M. Biermann,
D. W. Evans,
L. Eyer,
F. Jansen,
S. A. Klioner,
U. Lammers,
X. Luri,
F. Mignard,
C. Panem,
D. Pourbaix,
S. Randich,
P. Sartoretti,
H. I. Siddiqui,
C. Soubiran
, et al. (567 additional authors not shown)
Abstract:
Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the ast…
▽ More
Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.
△ Less
Submitted 3 March, 2017;
originally announced March 2017.
-
Interactive (statistical) visualisation and exploration of a billion objects with Vaex
Authors:
Maarten A. Breddels
Abstract:
With new catalogues arriving such as the Gaia DR1, containing more than a billion objects, new methods of handling and visualizing these data volumes are needed. In visualization, one problem is that the number of datapoints can become so large, that a scatter plot becomes cluttered. Another problem is that with over a billion objects, only a few cpu cycles are available per object if one wants to…
▽ More
With new catalogues arriving such as the Gaia DR1, containing more than a billion objects, new methods of handling and visualizing these data volumes are needed. In visualization, one problem is that the number of datapoints can become so large, that a scatter plot becomes cluttered. Another problem is that with over a billion objects, only a few cpu cycles are available per object if one wants to process them within a second, making traditional methods by rendering glyphs not viable. Instead, we show that by calculating statistics on a regular (N-dimensional) grid, visualizations of a billion objects can be done within a second on a modern desktop computer. This is achieved using memory mapping of hdf5 files together with a simple binning algorithm, which are part of a Python library called vaex. This enables efficient exploration or large datasets interactively, making science exploration of large catalogues feasible. Vaex is a Python library, which also integrates well in the Jupyter/Numpy/Astropy/matplotlib stack. Build on top of this is the vaex application, which allows for interactive exploration and visualization. The motivation for developing vaex is the catalogue of the Gaia satellite, however, vaex can also be used on SPH or N-body simulations, any other (future) catalogues such as SDSS, Pan-STARRS, LSST, WISE, 2MASS, etc. or other tabular data. The homepage for vaex is http://vaex.astro.rug.nl.
△ Less
Submitted 12 December, 2016;
originally announced December 2016.
-
A box full of chocolates: The rich structure of the nearby stellar halo revealed by Gaia and RAVE
Authors:
Amina Helmi,
Jovan Veljanoski,
Maarten A. Breddels,
Hao Tian,
Laura V. Sales
Abstract:
The hierarchical structure formation model predicts that stellar halos should form, at least partly, via mergers. If this was a predominant formation channel for the Milky Way's halo, imprints of this merger history in the form of moving groups or streams should exist also in the vicinity of the Sun. Here we study the kinematics of halo stars in the Solar neighbourhood using the very recent first…
▽ More
The hierarchical structure formation model predicts that stellar halos should form, at least partly, via mergers. If this was a predominant formation channel for the Milky Way's halo, imprints of this merger history in the form of moving groups or streams should exist also in the vicinity of the Sun. Here we study the kinematics of halo stars in the Solar neighbourhood using the very recent first data release from the Gaia mission, and in particular the TGAS dataset, in combination with data from the RAVE survey. Our aim is to determine the amount of substructure present in the phase-space distribution of halo stars that could be linked to merger debris. To characterise kinematic substructure, we measure the velocity correlation function in our sample of halo (low metallicity) stars. We also study the distribution of these stars in the space of energy and two components of the angular momentum, in what we call "Integrals of Motion" space. The velocity correlation function reveals substructure in the form of an excess of pairs of stars with similar velocities, well above that expected for a smooth distribution. Comparison to cosmological simulations of the formation of stellar halos indicate that the levels found are consistent with the Galactic halo having been built fully via accretion. Similarly, the distribution of stars in the space of "Integrals of motion" is highly complex. A strikingly high fraction (between 58% and upto 73%) of the stars that are somewhat less bound than the Sun are on (highly) retrograde orbits. A simple comparison to Milky Way-mass galaxies in cosmological hydrodynamical simulations suggests that less than 1% have such prominently retrograde outer halos. We also identify several other statistically significant structures in "Integrals of Motion" space that could potentially be related to merger events.
△ Less
Submitted 23 December, 2016; v1 submitted 1 November, 2016;
originally announced November 2016.
-
Matching the dark matter profiles of dSph galaxies with those of simulated satellites: a two parameter comparison
Authors:
Maarten A. Breddels,
Carlos Vera-Ciro,
Amina Helmi
Abstract:
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 clo…
▽ More
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.
△ Less
Submitted 27 July, 2015; v1 submitted 14 July, 2015;
originally announced July 2015.
-
Complexity on dwarf galaxies scale: A bimodal distribution function in Sculptor
Authors:
Maarten A. Breddels,
Amina Helmi
Abstract:
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 bimo…
▽ More
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.
△ Less
Submitted 16 September, 2014; v1 submitted 23 April, 2014;
originally announced April 2014.
-
An analytic distribution function for a massless cored stellar system in a cuspy dark matter halo
Authors:
Maarten A. Breddels,
Amina Helmi
Abstract:
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…
▽ More
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.
△ Less
Submitted 7 August, 2013;
originally announced August 2013.
-
Model comparison of the dark matter profiles of Fornax, Sculptor, Carina and Sextans
Authors:
M. A. Breddels,
A. Helmi
Abstract:
We use orbit based dynamical models to fit the 2nd and 4th moments of the line of sight velocity distributions of the Fornax, Sculptor, Carina and Sextans dwarf spheroidal galaxies. Our goal is to compare dark matter profile models of these four systems using Bayesian evidence. We consider NFW, Einasto and several cored profiles for their dark halos and present the probability distribution functio…
▽ More
We use orbit based dynamical models to fit the 2nd and 4th moments of the line of sight velocity distributions of the Fornax, Sculptor, Carina and Sextans dwarf spheroidal galaxies. Our goal is to compare dark matter profile models of these four systems using Bayesian evidence. We consider NFW, Einasto and several cored profiles for their dark halos and present the probability distribution functions of the model parameters. When considering each system separately, we find there is no preference for one of these specific parametric density profiles. However, the combined evidence shows that it is unlikely that all galaxies are embedded in the same type of cored profiles of the form ρ_{DM} \propto 1/(1 + r^2)^{β/2}, where β= 3, 4. For each galaxy, we also obtain an almost model independent, and therefore accurate, measurement of the logarithmic slope of the dark matter density distribution at a radius ~ r_{-3}, i.e. where the logarithmic slope of the stellar density profile is -3. This slope ranges from ~ -1.4 for Fornax to ~ -1.1 for Sextans, both at ~ 1 kpc. All our best fit models essentially have the same mass distribution over a large range in radius (from just below r_{-3} to the last measured data point). This remarkable finding likely implies much stronger constraints on the characteristics that subhalos extracted from cosmological simulations should have in order to host the dSph galaxies around the Milky Way.
△ Less
Submitted 10 April, 2013;
originally announced April 2013.
-
Orbit-based dynamical models of the Sculptor dSph galaxy
Authors:
Maarten A. Breddels,
A. Helmi,
R. C. E. van den Bosch,
G. van de Ven,
G. Battaglia
Abstract:
We have developed spherically symmetric dynamical models of dwarf spheroidal galaxies using Schwarzschild's orbit superposition 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 (e.g. velocity anisotropy). This method is thus less prone to biases introduced by assumptions in…
▽ More
We have developed spherically symmetric dynamical models of dwarf spheroidal galaxies using Schwarzschild's orbit superposition 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 (e.g. velocity anisotropy). This method is thus less prone to biases introduced by assumptions in comparison to the more commonly used Jeans modelling, and it allows us to derive the dark matter content in a robust way. Here we present our results for the Sculptor dwarf spheroidal galaxy, after testing our methods on mock data sets. We fit both the second and fourth velocity moment profile to break the mass-anisotropy degeneracy. For an NFW dark matter halo profile, we find that the mass of Sculptor within 1 kpc is M_1kpc = 1.03 \pm 0.07 x 10^8 Msol, and that its velocity anisotropy profile is tangentially biased and nearly constant with radius. The preferred concentration (c ~ 15) is low for its dark matter mass but consistent within the scatter found in N-body cosmological simulations. When we let the value of the central logarithmic slope α vary, we find that the best-fit model has α = 0, although an NFW cusp or shallower is consistent at 1 σ confidence level. On the other hand, very cuspy density profiles with logarithmic central slopes α < -1.5 are strongly disfavoured for Sculptor.
△ Less
Submitted 6 August, 2013; v1 submitted 21 May, 2012;
originally announced May 2012.
-
Not too big, not too small: the dark halos of the dwarf spheroidals in the Milky Way
Authors:
Carlos A. Vera-Ciro,
Amina Helmi,
Else Starkenburg,
Maarten A. Breddels
Abstract:
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…
▽ More
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.
△ Less
Submitted 1 November, 2012; v1 submitted 27 February, 2012;
originally announced February 2012.
-
Schwarzschild models of the Sculptor dSph galaxy
Authors:
Maarten A. Breddels,
Amina Helmi,
R. C. E. van den Bosch,
G. van de Ven,
G. Battaglia
Abstract:
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…
▽ More
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.
△ Less
Submitted 21 October, 2011;
originally announced October 2011.
-
Testing formation mechanisms of the Milky Way's thick disc with RAVE
Authors:
Michelle Wilson,
Amina Helmi,
H. L. Morrison,
M. A. Breddels,
O. Bienayme,
J. Binney,
J. Bland-Hawthorn,
R. Campbell,
K. C. Freeman,
J. P. Fulbright,
B. K. Gibson,
G. Gilmore,
E. K. Grebel,
U. Munari,
J. F. Navarro,
Q. A. Parker,
W. Reid,
G. Seabroke,
A. Siebert,
A. Siviero,
M. Steinmetz,
M. E. K. Williams,
R. F. G. Wyse,
T. Zwitter
Abstract:
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 o…
▽ More
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.
△ Less
Submitted 10 September, 2010;
originally announced September 2010.
-
Distance determination for RAVE stars using stellar models II: Most likely values assuming a standard stellar evolution scenario
Authors:
T. Zwitter,
G. Matijevič,
M. A. Breddels,
M. C. Smith,
A. Helmi,
U. Munari,
O. Bienaymé,
J. Binney,
J. Bland-Hawthorn,
C. Boeche,
A. G. A. Brown,
R. Campbell,
K. C. Freeman,
J. Fulbright,
B. Gibson,
G. Gilmore,
E. K. Grebel,
J. F. Navarro,
Q. A. Parker,
G. M. Seabroke,
A. Siebert,
A. Siviero,
M. Steinmetz,
F. G. Watson,
M. Williams
, et al. (1 additional authors not shown)
Abstract:
The RAdial Velocity Experiment (RAVE) is a spectroscopic survey of the Milky Way. We use the subsample of spectra with spectroscopically determined values of stellar parameters to determine the distances to these stars. The list currently contains 235,064 high quality spectra which show no peculiarities and belong to 210,872 different stars. The numbers will grow as the RAVE survey progresses. The…
▽ More
The RAdial Velocity Experiment (RAVE) is a spectroscopic survey of the Milky Way. We use the subsample of spectra with spectroscopically determined values of stellar parameters to determine the distances to these stars. The list currently contains 235,064 high quality spectra which show no peculiarities and belong to 210,872 different stars. The numbers will grow as the RAVE survey progresses. The public version of the catalog will be made available through the CDS services along with the ongoing RAVE public data releases.
The distances are determined with a method based on the work by Breddels et al.~(2010). Here we assume that the star undergoes a standard stellar evolution and that its spectrum shows no peculiarities. The refinements include: the use of either of the three isochrone sets, a better account of the stellar ages and masses, use of more realistic errors of stellar parameter values, and application to a larger dataset. The derived distances of both dwarfs and giants match within ~21% to the astrometric distances of Hipparcos stars and to the distances of observed members of open and globular clusters. Multiple observations of a fraction of RAVE stars show that repeatability of the derived distances is even better, with half of the objects showing a distance scatter of \simlt 11%.
RAVE dwarfs are ~300 pc from the Sun, and giants are at distances of 1 to 2 kpc, and up to 10 kpc. This places the RAVE dataset between the more local Geneva-Copenhagen survey and the more distant and fainter SDSS sample. As such it is ideal to address some of the fundamental questions of Galactic structure and evolution in the pre-Gaia era. Individual applications are left to separate papers, here we show that the full 6-dimensional information on position and velocity is accurate enough to discuss the vertical structure and kinematic properties of the thin and thick disks.
△ Less
Submitted 13 October, 2010; v1 submitted 26 July, 2010;
originally announced July 2010.
-
Distance determination for RAVE stars using stellar models
Authors:
M. A. Breddels,
M. C. Smith,
A. Helmi,
O. Bienayme,
J. Binney,
J. Bland-Hawthorn,
C. Boeche,
B. C. M. Burnett,
R. Campbell,
K. C. Freeman,
B. Gibson,
G. Gilmore,
E. K. Grebel,
U. Munari,
J. F. Navarro,
Q. A. Parker,
G. M. Seabroke,
A. Siebert,
A. Siviero,
M. Steinmetz,
F. G. Watson,
M. Williams,
R. F. G. Wyse,
T. Zwitter
Abstract:
(Abridged) Aims:We develop a method for deriving distances from spectroscopic data and obtaining full 6D phase-space coordinates for the RAVE survey's second data release. Methods: We used stellar models combined with atmospheric properties from RAVE (Teff, logg and [Fe/H]) and (J-Ks) photometry from archival sources to derive absolute magnitudes. We are able to derive the full 6D phase-space co…
▽ More
(Abridged) Aims:We develop a method for deriving distances from spectroscopic data and obtaining full 6D phase-space coordinates for the RAVE survey's second data release. Methods: We used stellar models combined with atmospheric properties from RAVE (Teff, logg and [Fe/H]) and (J-Ks) photometry from archival sources to derive absolute magnitudes. We are able to derive the full 6D phase-space coordinates for a large sample of RAVE stars. This method is tested with artificial data, Hipparcos trigonometric parallaxes and observations of the open cluster M67. Results: When we applied our method to a set of 16 146 stars, we found that 25% (4 037) of the stars have relative (statistical) distance errors of < 35%, while 50% (8 073) and 75% (12 110) have relative (statistical) errors smaller than 45% and 50%, respectively. Our various tests show that we can reliably estimate distances for main-sequence stars, but there is an indication of potential systematic problems with giant stars. For the main-sequence star sample (defined as those with log(g) > 4), 25% (1 744) have relative distance errors < 31%, while 50% (3 488) and 75% (5 231) have relative errors smaller than 36% and 42%, respectively. Our full dataset shows the expected decrease in the metallicity of stars as a function of distance from the Galactic plane. The known kinematic substructures in the U and V velocity components of nearby dwarf stars are apparent in our dataset, confirming the accuracy of our data and the reliability of our technique. We provide independent measurements of the orientation of the UV velocity ellipsoid and of the solar motion, and they are in very good agreement with previous work. Conclusions: The distance catalogue for the RAVE second data release is available at http://www.astro.rug.nl/~rave
△ Less
Submitted 3 March, 2010;
originally announced March 2010.
