Source author record

Matthias Steinmetz

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

74works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

74 published item(s)

preprint2022arXiv

Cold and Hot gas distribution around the Milky-Way-M31 system in the HESTIA simulations

Recent observations have revealed remarkable insights into the gas reservoir in the circumgalactic medium (CGM) of galaxy haloes. In this paper, we characterise the gas in the vicinity of Milky Way and Andromeda analogues in the HESTIA (High resolution Environmental Simulations of The Immediate Area) suite of constrained Local Group (LG) simulations. The HESTIA suite comprise of a set of three high-resolution {\sc arepo}-based simulations of the LG, run using the Auriga galaxy formation model. For this paper, we focus only on the $z = 0$ simulation datasets and generate mock skymaps along with a power spectrum analysis to show that the distributions of ions tracing low-temperature gas (HI and SiIII) are more clumpy in comparison to warmer gas tracers (OVI, OVII and OVIII). We compare to the spectroscopic CGM observations of M31 and low-redshift galaxies. HESTIA under-produces the column densities of the M31 observations, but the simulations are consistent with the observations of low-redshift galaxies. A possible explanation for these findings is that the spectroscopic observations of M31 are contaminated by gas residing in the CGM of the Milky Way.

preprint2022arXiv

Hermeian haloes: Field haloes that interacted with both the Milky Way and M31

The Local Group is a unique environment in which to study the astrophysics of galaxy formation. The proximity of the Milky Way and M31 enhances the frequency of interactions of the low-mass halo population with more massive dark matter haloes, which increases their concentrations and strips them of gas and other material. Some low-mass haloes pass through the haloes of the Milky Way or M31 and are either ejected into the field or exchanged between the two primary hosts. We use high resolution gas-dynamical simulations to describe a new class of field haloes that passed through the haloes of both the Milky Way and M31 at early times and are almost twice as concentrated as field haloes that do not interact with the primary pair. These 'Hermeian' haloes are distributed anisotropically at larger distances from the Local Group barycentre than the primary haloes and appear to cluster along the line connecting the Milky Way and M31. Hermeian haloes facilitate the exchange of dark matter, gas, and stars between the Milky Way and M31 and can enhance the star formation rate of the gas in the primary haloes during their interactions with them. We also show that some Hermeian haloes can host galaxies that, because they are embedded in haloes that are more concentrated than regular field haloes, are promising targets for indirect dark matter searches beyond the Milky Way virial radius and can produce signals that are competitive with those of some dwarf galaxies. Hermeian galaxies in the Local Group should be detectable by forthcoming wide-field imaging surveys.

preprint2022arXiv

The MUSE Hubble Ultra Deep Field Survey XVI. The angular momentum of low-mass star-forming galaxies. A cautionary tale and insights from TNG50

We investigate the specific angular momentum (sAM) $ j(<r)$ profiles of intermediate redshift ($0.4<z<1.4$) star-forming galaxies (SFGs) in the relatively unexplored regime of low masses (down to $M_\star\sim 10^8$M$_{\odot}$), and small sizes (down to $R_{\rm e}\sim 1.5$ kpc) and characterize the sAM scaling relation and its redshift evolution. We have developed a 3D methodology to constrain sAM profiles of the star-forming gas using a forward modeling approach with \galpak{} that incorporates the effects of beam smearing, yielding the intrinsic morpho-kinematic properties even with limited spatial resolution data. Using mock observations from the TNG50 simulation, we find that our 3D methodology robustly recovers the star formation rate (SFR)-weighted $j(<r)$ profiles down to low effective signal-to-noise ratio (SNR) of $\gtrapprox3$. We applied our methodology blindly to a sample of 494 \OII{}-selected SFGs in the MUSE Ultra Deep Field (UDF) 9~arcmin$^2$ mosaic data, covering the unexplored $8<\log M_*/$M$_{\odot}<9$ mass range. We find that the (SFR-weighted) sAM relation follows $j\propto M_\star^α$ with an index $α$ varying from $α=0.3$ to $α=0.5$, from $\log M_\star/$M$_{\odot}=8$ to $\log M_*/$M$_{\odot}=10.5$. The UDF sample supports a redshift evolution consistent with the $(1+z)^{-0.5}$ expectation from a Universe in expansion. The scatter of the sAM sequence is a strong function of the dynamical state with $\log j|_{M_*}\propto 0.65 \times \log(V_{\rm max}/σ)$ where $σ$ is the velocity dispersion at $2 R_{\rm e}$. In TNG50, SFGs also form a $j-M_{\star}-(V/σ)$ plane but it correlates more with galaxy size than with morphological parameters. Our results suggest that SFGs might experience a dynamical transformation before their morphological transformation to becoming passive via either merging or secular evolution.

preprint2021arXiv

Searching for Light in the Darkness: Bounds on ALP Dark Matter with the optical MUSE-Faint survey

We use MUSE spectroscopic observations of the dwarf spheroidal galaxy Leo T between 470 and 935 nm to search for radiative decays of axion like particles (ALPs). Under the assumption that ALPs constitute the dark matter component of the Leo T halo, we derive bounds on the effective ALP-two-photon coupling. We improve existing limits by more than one order of magnitude in the ALP mass range 2.7-5.3 eV.

preprint2021arXiv

The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) Survey Design, Reductions, and Detections

We describe the survey design, calibration, commissioning, and emission-line detection algorithms for the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The goal of HETDEX is to measure the redshifts of over a million Ly$α$ emitting galaxies between 1.88<z<3.52, in a 540 deg^2 area encompassing a co-moving volume of 10.9 Gpc^3. No pre-selection of targets is involved; instead the HETDEX measurements are accomplished via a spectroscopic survey using a suite of wide-field integral field units distributed over the focal plane of the telescope. This survey measures the Hubble expansion parameter and angular diameter distance, with a final expected accuracy of better than 1%. We detail the project's observational strategy, reduction pipeline, source detection, and catalog generation, and present initial results for science verification in the COSMOS, Extended Groth Strip, and GOODS-N fields. We demonstrate that our data reach the required specifications in throughput, astrometric accuracy, flux limit, and object detection, with the end products being a catalog of emission-line sources, their object classifications, and flux-calibrated spectra.

preprint2021arXiv

The MUSE Extremely Deep Field: Evidence for SFR-induced cores in dark-matter dominated galaxies at z=1

Disc-halo decomposition on rotationally supported star-forming galaxies (SFGs) at $z>1$ are often limited to massive galaxies ($M_\star>10^{10}~M_\odot$) and rely on either deep Integral Field Spectroscopy data or stacking analyses. We present a study of the dark matter (DM) content of nine $z\approx1$ SFGs selected Using the brightest [OII] emitters in the deepest Multi-Unit Spectrograph Explorer (MUSE) field to date, namely the 140hr MUSE Extremely Deep Field, we perform disk-halo decompositions on 9 low-mass SFGs (with $10^{8.5}<M_\star<10^{10.5}~M_\odot$) using a novel 3D modeling approach, which together with the exquisite S/N allows us to measure individual rotation curves to $3\times R_e$. The disk-halo decomposition includes a stellar, DM, gas, and occasionally a bulge component. The DM component primarily uses the generalized $α,β,γ$ profile or a Navarro-Frenk-White (NFW) profile. The disk stellar masses $M_\star$ obtained from the [OII] disk-halo decomposition agree with the values inferred from the spectral energy distributions. While the rotation curves show diverse shapes, ranging from rising to declining at large radii, the DM fractions within the half-light radius $f_{\rm DM}(<R_e)$ are found to be 60\% to 95\%, extending to lower masses (densities) recent results on massive SFGs with $M_\star>10^{10}~M_\odot$. The DM halos show constant surface densities of $\sim100~M_\odot$ pc$^{-2}$. Half of the sample shows a strong preference for cored over cuspy DM profiles. The presence of DM cores appears to be related to galaxies with stellar-to-halo mass $\log M_\star/M_{\rm vir}\approx-2.5$. In addition, the cuspiness of the DM profiles is found to be a strong function of the recent star-formation activity. Both of these results are interpreted as evidence for feedback-induced core formation in the Cold Dark Matter context.

preprint2020arXiv

A LOFAR Observation of Ionospheric Scintillation from Two Simultaneous Travelling Ionospheric Disturbances

This paper presents the results from one of the first observations of ionospheric scintillation taken using the Low-Frequency Array (LOFAR). The observation was of the strong natural radio source Cas A, taken overnight on 18-19 August 2013, and exhibited moderately strong scattering effects in dynamic spectra of intensity received across an observing bandwidth of 10-80MHz. Delay-Doppler spectra (the 2-D FFT of the dynamic spectrum) from the first hour of observation showed two discrete parabolic arcs, one with a steep curvature and the other shallow, which can be used to provide estimates of the distance to, and velocity of, the scattering plasma. A cross-correlation analysis of data received by the dense array of stations in the LOFAR "core" reveals two different velocities in the scintillation pattern: a primary velocity of ~30m/s with a north-west to south-east direction, associated with the steep parabolic arc and a scattering altitude in the F-region or higher, and a secondary velocity of ~110m/s with a north-east to south-west direction, associated with the shallow arc and a scattering altitude in the D-region. Geomagnetic activity was low in the mid-latitudes at the time, but a weak sub-storm at high latitudes reached its peak at the start of the observation. An analysis of Global Navigation Satellite Systems (GNSS) and ionosonde data from the time reveals a larger-scale travelling ionospheric disturbance (TID), possibly the result of the high-latitude activity, travelling in the north-west to south-east direction, and, simultaneously, a smaller--scale TID travelling in a north-east to south-west direction, which could be associated with atmospheric gravity wave activity. The LOFAR observation shows scattering from both TIDs, at different altitudes and propagating in different directions. To the best of our knowledge this is the first time that such a phenomenon has been reported.

preprint2020arXiv

The $R$-Process Alliance: Fourth Data Release from the Search for $r$-Process-Enhanced Stars in the Galactic Halo

This compilation is the fourth data release from the $R$-Process Alliance (RPA) search for $r$-process-enhanced stars, and the second release based on "snapshot" high-resolution ($R \sim 30,000$) spectra collected with the du Pont 2.5m Telescope. In this data release, we propose a new delineation between the $r$-I and $r$-II stellar classes at $\mathrm{[Eu/Fe]} = +0.7$, instead of the empirically chosen $\mathrm{[Eu/Fe]} = +1.0$ level previously in use, based on statistical tests of the complete set of RPA data released to date. We also statistically justify the minimum level of [Eu/Fe] for definition of the $r$-I stars, [Eu/Fe] $> +0.3$. Redefining the separation between $r$-I and $r$-II stars will aid in analysis of the possible progenitors of these two classes of stars and whether these signatures arise from separate astrophysical sources at all. Applying this redefinition to previous RPA data, the number of identified $r$-II and $r$-I stars changes to 51 and 121, respectively, from the initial set of data releases published thus far. In this data release, we identify 21 new $r$-II, 111 new $r$-I (plus three re-identified), and 7 new (plus one re-identified) limited-$r$ stars out of a total of 232 target stars, resulting in a total sample of 72 new $r$-II stars, 232 new $r$-I stars, and 42 new limited-$r$ stars identified by the RPA to date.

preprint2020arXiv

The Hestia project: simulations of the Local Group

We present the Hestia simulation suite: High-resolutions Environmental Simulations of The Immediate Area, a set of cosmological simulations of the Local Group. Initial conditions constrained by the observed peculiar velocity of nearby galaxies are employed to accurately simulate the local cosmography. Halo pairs that resemble the Local Group are found in low resolutions constrained, dark matter only simulations, and selected for higher resolution magneto hydrodynamic simulation using the Arepo code. Baryonic physics follows the Auriga model of galaxy formation. The simulations contain a high resolution region of 3-5 Mpc in radius from the Local Group midpoint embedded in the correct cosmographic landscape. Within this region a simulated Local Group consisting of a Milky Way and Andromeda like galaxy forms, whose description is in excellent agreement with observations. The simulated Local Group galaxies resemble the Milky Way and Andromeda in terms of their halo mass, mass ratio, stellar disc mass, morphology separation, relative velocity, rotation curves, bulge-disc morphology, satellite galaxy stellar mass function, satellite radial distribution and in some cases, the presence of a Magellanic cloud like object. Because these simulations properly model the Local Group in their cosmographic context, they provide a testing ground for questions where environment is thought to play an important role.

preprint2020arXiv

The Sixteenth Data Release of the Sloan Digital Sky Surveys: First Release from the APOGEE-2 Southern Survey and Full Release of eBOSS Spectra

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

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.

preprint2016arXiv

Chemical separation of disc components using RAVE

We present evidence from the RAdial Velocity Experiment (RAVE) survey of chemically separated, kinematically distinct disc components in the solar neighbourhood. We apply probabilistic chemical selection criteria to separate our sample into $α$-low (`thin disc') and $α$-high (`thick disc') sequences. Using newly derived distances, which will be utilized in the upcoming RAVE DR5, we explore the kinematic trends as a function of metallicity for each of the disc components. For our $α$-low disc, we find a negative trend in the mean rotational velocity ($V_{\mathrmϕ}$) as a function of iron abundance ([Fe/H]). We measure a positive gradient $\partial V_{\mathrmϕ}$/$\partial$[Fe/H] for the $α$-high disc, consistent with results from high-resolution surveys. We also find differences between the $α$-low and $α$-high discs in all three components of velocity dispersion. We discuss the implications of an $α$-low, metal-rich population originating from the inner Galaxy, where the orbits of these stars have been significantly altered by radial mixing mechanisms in order to bring them into the solar neighbourhood. The probabilistic separation we propose can be extended to other data sets for which the accuracy in [$α$/Fe] is not sufficient to disentangle the chemical disc components a priori. For such datasets which will also have significant overlap with Gaia DR1, we can therefore make full use of the improved parallax and proper motion data as it becomes available to investigate kinematic trends in these chemical disc components.

preprint2016arXiv

Constrained Local UniversE Simulations: A Local Group Factory

Near field cosmology is practiced by studying the Local Group (LG) and its neighbourhood. The present paper describes a framework for simulating the near field on the computer. Assuming the LCDM model as a prior and applying the Bayesian tools of the Wiener filter (WF) and constrained realizations of Gaussian fields to the Cosmicflows-2 (CF2) survey of peculiar velocities, constrained simulations of our cosmic environment are performed. The aim of these simulations is to reproduce the LG and its local environment. Our main result is that the LG is likely a robust outcome of the LCDM scenario when subjected to the constraint derived from CF2 data, emerging in an environment akin to the observed one. Three levels of criteria are used to define the simulated LGs. At the base level, pairs of halos must obey specific isolation, mass and separation criteria. At the second level the orbital angular momentum and energy are constrained and on the third one the phase of the orbit is constrained. Out of the 300 constrained simulations 146 LGs obey the first set of criteria, 51 the second and 6 the third. The robustness of our LG factory enables the construction of a large ensemble of simulated LGs. Suitable candidates for high resolution hydrodynamical simulations of the LG can be drawn from this ensemble, which can be used to perform comprehensive studies of the formation of the LG

preprint2016arXiv

MusE GAs FLOw and Wind (MEGAFLOW) I: First MUSE results on background quasars

The physical properties of galactic winds are one of the keys to understand galaxy formation and evolution. These properties can be constrained thanks to background quasar lines of sight (LOS) passing near star-forming galaxies (SFGs). We present the first results of the MusE GAs FLOw and Wind (MEGAFLOW) survey obtained of 2 quasar fields which have 8 MgII absorbers of which 3 have rest-equivalent width greater than 0.8 Å. With the new Multi Unit Spectroscopic Explorer (MUSE) spectrograph on the Very Large Telescope (VLT), we detect 6 (75$\%$) MgII host galaxy candidates withing a radius of 30 arcsec from the quasar LOS. Out of these 6 galaxy--quasar pairs, from geometrical arguments, one is likely probing galactic outflows, two are classified as "ambiguous", two are likely probing extended gaseous disks and one pair seems to be a merger. We focus on the wind$-$pair and constrain the outflow using a high resolution quasar spectra from Ultraviolet and Visual Echelle Spectrograph (UVES). Assuming the metal absorption to be due to gas flowing out of the detected galaxy through a cone along the minor axis, we find outflow velocities of the order of $\approx$ 150 km/s (i.e. smaller than the escape velocity) with a loading factor, $η=\dot M_{\rm out}/$SFR, of $\approx$ 0.7. We see evidence for an open conical flow, with a low-density inner core. In the future, MUSE will provide us with about 80 multiple galaxy$-$quasar pairs in two dozen fields.

preprint2016arXiv

SPIDERS: the spectroscopic follow-up of X-ray selected clusters of galaxies in SDSS-IV

SPIDERS (The SPectroscopic IDentification of eROSITA Sources) is a program dedicated to the homogeneous and complete spectroscopic follow-up of X-ray AGN and galaxy clusters over a large area ($\sim$7500 deg$^2$) of the extragalactic sky. SPIDERS is part of the SDSS-IV project, together with the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) and the Time-Domain Spectroscopic Survey (TDSS). This paper describes the largest project within SPIDERS before the launch of eROSITA: an optical spectroscopic survey of X-ray selected, massive ($\sim 10^{14}$ to $10^{15}~M_{\odot}$) galaxy clusters discovered in ROSAT and XMM-Newton imaging. The immediate aim is to determine precise ($Δ_z \sim 0.001$) redshifts for 4,000-5,000 of these systems out to $z \sim 0.6$. The scientific goal of the program is precision cosmology, using clusters as probes of large-scale structure in the expanding Universe. We present the cluster samples, target selection algorithms and observation strategies. We demonstrate the efficiency of selecting targets using a combination of SDSS imaging data, a robust red-sequence finder and a dedicated prioritization scheme. We describe a set of algorithms and work-flow developed to collate spectra and assign cluster membership, and to deliver catalogues of spectroscopically confirmed clusters. We discuss the relevance of line-of-sight velocity dispersion estimators for the richer systems. We illustrate our techniques by constructing a catalogue of 230 spectroscopically validated clusters ($0.031 < z < 0.658$), found in pilot observations. We discuss two potential science applications of the SPIDERS sample: the study of the X-ray luminosity-velocity dispersion ($L_X-σ$) relation and the building of stacked phase-space diagrams.

preprint2016arXiv

The LOFAR Solar Imaging Pipeline and the LOFAR Solar Data Center

LOFAR is a new and sensitive radio interferometer that can be used for dynamic high-resolution imaging spectroscopy at low radio frequencies from 10 to 90 and 110 to 250 MHz. Here we describe its usage for observations of the Sun and in particular of solar radio bursts. We also describe the processing, archiving and accessing of solar LOFAR data, which is accomplished via the LOFAR Solar Imaging Pipeline and the LOFAR Solar Data Center.

preprint2015arXiv

Chemical Cartography with APOGEE: Metallicity Distribution Functions and the Chemical Structure of the Milky Way Disk

Using a sample of 69,919 red giants from the SDSS-III/APOGEE Data Release 12, we measure the distribution of stars in the [$α$/Fe] vs. [Fe/H] plane and the metallicity distribution functions (MDF) across an unprecedented volume of the Milky Way disk, with radius $3<R<15$ kpc and height $|z|<2$ kpc. Stars in the inner disk ($R<5$ kpc) lie along a single track in [$α$/Fe] vs. [Fe/H], starting with $α$-enhanced, metal-poor stars and ending at [$α$/Fe]$\sim0$ and [Fe/H]$\sim+0.4$. At larger radii we find two distinct sequences in [$α$/Fe] vs. [Fe/H] space, with a roughly solar-$α$ sequence that spans a decade in metallicity and a high-$α$ sequence that merges with the low-$α$ sequence at super-solar [Fe/H]. The location of the high-$α$ sequence is nearly constant across the disk, however there are very few high-$α$ stars at $R>11$ kpc. The peak of the midplane MDF shifts to lower metallicity at larger $R$, reflecting the Galactic metallicity gradient. Most strikingly, the shape of the midplane MDF changes systematically with radius, with a negatively skewed distribution at $3<R<7$ kpc, to a roughly Gaussian distribution at the solar annulus, to a positively skewed shape in the outer Galaxy. For stars with $|z|>1$ kpc or [$α$/Fe]$>0.18$, the MDF shows little dependence on $R$. The positive skewness of the outer disk MDF may be a signature of radial migration; we show that blurring of stellar populations by orbital eccentricities is not enough to explain the reversal of MDF shape but a simple model of radial migration can do so.

preprint2015arXiv

Cosmicflows Constrained Local UniversE Simulations

This paper combines observational datasets and cosmological simulations to generate realistic numerical replicas of the nearby Universe. These latter are excellent laboratories for studies of the non-linear process of structure formation in our neighborhood. With measurements of radial peculiar velocities in the Local Universe (cosmicflows-2) and a newly developed technique, we produce Constrained Local UniversE Simulations (CLUES). To assess the quality of these constrained simulations, we compare them with random simulations as well as with local observations. The cosmic variance, defined as the mean one-sigma scatter of cell-to-cell comparison between two fields, is significantly smaller for the constrained simulations than for the random simulations. Within the inner part of the box where most of the constraints are, the scatter is smaller by a factor 2 to 3 on a 5 Mpc/h scale with respect to that found for random simulations. This one-sigma scatter obtained when comparing the simulated and the observation-reconstructed velocity fields is only 104 +/- 4 km/s i.e. the linear theory threshold. These two results demonstrate that these simulations are in agreement with each other and with the observations of our neighborhood. For the first time, simulations constrained with observational radial peculiar velocities resemble the Local Universe up to a distance of 150 Mpc/h on a scale of a few tens of megaparsecs. When focusing on the inner part of the box, the resemblance with our cosmic neighborhood extends to a few megaparsecs (< 5 Mpc/h). The simulations provide a proper Large Scale environment for studies of the formation of nearby objects.

preprint2015arXiv

Mergers and the outside-in formation of dwarf spheroidals

We use a cosmological simulation of the formation of the Local Group to explore the origin of age and metallicity gradients in dwarf spheroidal galaxies. We find that a number of simulated dwarfs form "outside-in", with an old, metal-poor population that surrounds a younger, more concentrated metal-rich component, reminiscent of dwarf spheroidals like Sculptor or Sextans. We focus on a few examples where stars form in two populations distinct in age in order to elucidate the origin of these gradients. The spatial distributions of the two components reflect their diverse origin; the old stellar component is assembled through mergers, but the young population forms largely in situ. The older component results from a first episode of star formation that begins early but is quickly shut off by the combined effects of stellar feedback and reionization. The younger component forms when a late accretion event adds gas and reignites star formation. The effect of mergers is to disperse the old stellar population, increasing their radius and decreasing their central density relative to the young population. We argue that dwarf-dwarf mergers offer a plausible scenario for the formation of systems with multiple distinct populations and, more generally, for the origin of age and metallicity gradients in dwarf spheroidals.

preprint2015arXiv

Planes of satellite galaxies and the cosmic web

Recent observational studies have demonstrated that the majority of satellite galaxies tend to orbit their hosts on highly flattened, vast, possibly co-rotating planes. Two nearly parallel planes of satellites have been confirmed around the M31 galaxy and around the Centaurus A galaxy, while the Milky Way also sports a plane of satellites. It has been argued that such an alignment of satellites on vast planes is unexpected in the standard (ΛCDM) model of cosmology if not even in contradiction to its generic predictions. Guided by ΛCDM numerical simulations, which suggest that satellites are channeled towards hosts along the axis of the slowest collapse as dictated by the ambient velocity shear tensor, we re-examine the planes of local satellites systems within the framework of the local shear tensor derived from the Cosmicflows-2 dataset. The analysis reveals that the Local Group and Centaurus A reside in a filament stretched by the Virgo cluster and compressed by the expansion of the Local Void. Four out of five thin planes of satellite galaxies are indeed closely aligned with the axis of compression induced by the Local Void. Being the less massive system, the moderate misalignment of the Milky Way's satellite plane can likely be ascribed to its greater susceptibility to tidal torques, as suggested by numerical simulations. The alignment of satellite systems in the local universe with the ambient shear field is thus in general agreement with predictions of the ΛCDM model.

preprint2015arXiv

Spectro-photometric distances to stars: a general-purpose Bayesian approach

We developed a code that estimates distances to stars using measured spectroscopic and photometric quantities. We employ a Bayesian approach to build the probability distribution function over stellar evolutionary models given these data, delivering estimates of model parameters for each star individually. The code was first tested on simulations, successfully recovering input distances to mock stars with <1% bias.The method-intrinsic random distance uncertainties for typical spectroscopic survey measurements amount to around 10% for dwarf stars and 20\% for giants, and are most sensitive to the quality of $\log g$ measurements. The code was validated by comparing our distance estimates to parallax measurements from the Hipparcos mission for nearby stars (< 300 pc), to asteroseismic distances of CoRoT red giant stars, and to known distances of well-studied open and globular clusters. The external comparisons confirm that our distances are subject to very small systematic biases with respect to the fundamental Hipparcos scale (+0.4 % for dwarfs, and +1.6% for giants). The typical random distance scatter is 18% for dwarfs, and 26% for giants. For the CoRoT-APOGEE sample, the typical random distance scatter is ~15%, both for the nearby and farther data. Our distances are systematically larger than the CoRoT ones by about +9%, which can mostly be attributed to the different choice of priors. The comparison to known distances of star clusters from SEGUE and APOGEE has led to significant systematic differences for many cluster stars, but with opposite signs, and with substantial scatter. Finally, we tested our distances against those previously determined for a high-quality sample of giant stars from the RAVE survey, again finding a small systematic trend of +5% and an rms scatter of 30%.

preprint2015arXiv

The alignment of galaxy spin with the shear field in observations

Tidal torque theory suggests that galaxies gain angular momentum in the linear stage of structure formation. Such a theory predicts alignments between the spin of haloes and tidal shear field. However, non-linear evolution and angular momentum acquisition may alter this prediction significantly. In this paper, we use a reconstruction of the cosmic shear field from observed peculiar velocities combined with spin axes extracted from galaxies within $115\, \mathrm{Mpc} $ ($\sim8000 \, {\mathrm {km}}{\mathrm s}^{-1}$) from 2MRS catalog, to test whether or not galaxies appear aligned with principal axes of shear field. Although linear reconstructions of the tidal field have looked at similar issues, this is the first such study to examine galaxy alignments with velocity-shear field. Ellipticals in the 2MRS sample, show a statistically significant alignment with two of the principal axes of the shear field. In general, elliptical galaxies have their short axis aligned with the axis of greatest compression and perpendicular to the axis of slowest compression. Spiral galaxies show no signal. Such an alignment is significantly strengthened when considering only those galaxies that are used in velocity field reconstruction. When examining such a subsample, a weak alignment with the axis of greatest compression emerges for spiral galaxies as well. This result indicates that although velocity field reconstructions still rely on fairly noisy and sparse data, the underlying alignment with shear field is strong enough to be visible even when small numbers of galaxies are considered - especially if those galaxies are used as constraints in the reconstruction.

preprint2015arXiv

The Effect of Environment on Milky Way-mass galaxies in a Constrained Simulation of the Local Group

In this letter we present, for the first time, a study of star formation rate, gas fraction and galaxy morphology of a constrained simulation of the Milky Way (MW) and Andromeda (M31) galaxies, compared to other MW-mass galaxies. By combining with unconstrained simulations we cover a sufficient volume to compare these galaxies environmental densities ranging from the field to that of the Local Group (LG). This is particularly relevant as it has been shown that, quite generally, galaxy properties depend intimately upon their environment, most prominently when galaxies in clusters are compared to those in the field. For galaxies in loose groups such as the LG, however, environmental effects have been less clear. We consider the galaxy's environmental density in spheres of 1200 kpc (comoving) and find that whilst environment does not appear to directly affect morphology, there is a positive trend with star formation rates. This enhancement in star formation occurs systematically for galaxies in higher density environments, regardless whether they are part of the LG or in filaments. Our simulations suggest that the richer environment at Mpc-scales may help replenish the star-forming gas, allowing higher specific star formation rates in galaxies such as the MW.

preprint2015arXiv

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

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

preprint2015arXiv

The Milky Way and Andromeda galaxies in a constrained hydrodynamical simulation: morphological evolution

We study the two main constituent galaxies of a constrained simulation of the Local Group as candidates for the Milky Way (MW) and Andromeda (M31). We focus on the formation of the stellar discs and its relation to the formation of the group as a rich system with two massive galaxies, and investigate the effects of mergers and accretion as drivers of morphological transformations. We use a state-of-the-art hydrodynamical code which includes star formation, feedback and chemical enrichment to carry out our study. We run two simulations, where we include or neglect the effects of radiation pressure from stars, to investigate the impact of this process on the morphologies and star formation rates of the simulated galaxies. We find that the simulated M31 and MW have different formation histories, even though both inhabit, at z=0, the same environment. These differences directly translate into and explain variations in their star formation rates, in-situ fractions and final morphologies. The M31 candidate has an active merger history, as a result of which its stellar disc is unable to survive unaffected until the present time. In contrast, the MW candidate has a smoother history with no major mergers at late times, and forms a disc that grows steadily; at z=0 the simulated MW has an extended, rotationally-supported disc which is dominant over the bulge. Our two feedback implementations predict similar evolution of the galaxies and their discs, although some variations are detected, the most important of which is the formation time of the discs: in the model with weaker/stronger feedback the discs form earlier/later. In summary, by comparing the formation histories of the two galaxies, we conclude that the particular merger/accretion history of a galaxy rather than its environment at the LG-scales is the main driver of the formation and subsequent growth or destruction of galaxy discs.

preprint2015arXiv

Unveiling the counter-rotating nature of the kinematically distinct core in NGC5813 with MUSE

MUSE observations of NGC5813 reveal a complex structure in the velocity dispersion map, previously hinted by SAURON observations. The structure is reminiscent of velocity dispersion maps of galaxies comprising two counter-rotating discs, and may explain the existence of the kinematically distinct core (KDC). Further evidence for two counter-rotating components comes from the analysis of the higher moments of the stellar line-of-sight velocity distributions and fitting MUSE spectra with two separate Gaussian line-of-sight velocity distributions. The emission-line kinematics show evidence of being linked to the present cooling flows and the buoyant cavities seen in X-rays. We detect ionised gas in a nuclear disc-like structure, oriented like the KDC, which is, however, not directly related to the KDC. We build an axisymmetric Schwarzschild dynamical model, which shows that the MUSE kinematics can be reproduced well with two counter-rotating orbit families, characterised by relatively low angular momentum components, but clearly separated in integral phase space and with radially varying contributions. The model indicates that the counter-rotating components in NGC5813 are not thin discs, but dynamically hot structures. Our findings give further evidence that KDCs in massive galaxies should not necessarily be considered as structurally or dynamically decoupled regions, but as the outcomes of the mixing of different orbital families, where the balance in the distribution of mass of the orbital families is crucial. We discuss the formation of the KDC in NGC5813 within the framework of gas accretion, binary mergers and formation of turbulent thick discs from cold streams at high redshift.

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

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

Overview of the SDSS-IV MaNGA Survey: Mapping Nearby Galaxies at Apache Point Observatory

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

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 distribution of gas in the Local Group from constrained cosmological simulations: the case for Andromeda and the Milky Way galaxies

We study the gas distribution in the Milky Way and Andromeda using a constrained cosmological simulation of the Local Group (LG) within the context of the CLUES (Constrained Local UniversE Simulations) project. We analyse the properties of gas in the simulated galaxies at $z=0$ for three different phases: `cold', `hot' and HI, and compare our results with observations. The amount of material in the hot halo ($M_{hot}\approx 4-5\times10^{10}\,$M$_{\odot}$), and the cold ($M_{cold}(r\lesssim10\,$kpc$)\approx10^{8}\,$M$_{\odot}$) and HI ($M_{HI}(r\lesssim50\,$kpc$)\approx 3-4\times10^8\,$M$_{\odot}$) components display a reasonable agreement with observations. We also compute the accretion/ejection rates together with the HI (radial and all-sky) covering fractions. The integrated HI accretion rate within $r=50\,$kpc gives $\sim$$0.2-0.3\,$M$_{\odot}\,$yr$^{-1}$, i.e. close to that obtained from high-velocity clouds in the Milky Way. We find that the global accretion rate is dominated by hot material, although ionized gas with $T\lesssim10^5\,$K can contribute significantly too. The $net$ accretion rates of $all$ material at the virial radii are $6-8\,$M$_{\odot}\,$yr$^{-1}$. At $z=0$, we find a significant gas excess between the two galaxies, as compared to any other direction, resulting from the overlap of their gaseous haloes. In our simulation, the gas excess first occurs at $z\sim1$, as a consequence of the kinematical evolution of the LG.

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 Imprint of Reionization on the Star Formation Histories of Dwarf Galaxies

We explore the impact of cosmic reionization on nearby isolated dwarf galaxies using a compilation of SFHs estimated from deep HST data and a cosmological hydrodynamical simulation of the Local Group. The nearby dwarfs show a wide diversity of star formation histories; from ancient systems that have largely completed their star formation $\sim 10$ Gyr ago to young dwarfs that have formed the majority of their stars in the past $\sim 5$ Gyr to two-component systems characterized by the overlap of comparable numbers of old and young stars. Taken as an ensemble, star formation in nearby dwarfs dips to lower-than-average rates at intermediate times ($4<t$/Gyr $<8$), a feature that we trace in the simulation to the effects of cosmic reionization. Reionization heats the gas and drives it out of the shallow potential wells of low mass halos, affecting especially those below a sharp mass threshold that corresponds to a virial temperature of $\sim 2 \times 10^4 $ $\mathrm{K}$ at $z_{\rm reion}$. The loss of baryons leads to a sharp decline in the star forming activity of early-collapsing systems, which, compounded by feedback from early star formation, empties halos of gas and leaves behind systems where a single old stellar component prevails. In halos below the threshold at $z_{\rm reion}$, reionization heating leads to a delay in the onset of star formation that lasts until the halo grows massive enough to allow some of the remaining gas to cool and form stars. Young stellar components therefore dominate in dwarfs whose halos assemble late and thus form few stars before reionization. Two-component systems may be traced to late mergers of individual examples of the two aforementioned cases. The relative dearth of intermediate-age stars in nearby dwarfs might thus be the clearest signature yet identified of the imprint of cosmic reionization on the star formation history of dwarf galaxies.

preprint2014arXiv

The Tenth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Apache Point Observatory Galactic Evolution Experiment

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

preprint2013arXiv

Cold versus Warm Dark Matter simulations of a galaxy group

The differences between cold (CDM) and warm (WDM) dark matter in the formation of a group of galaxies is examined by running two identical simulations where in the WDM case the initial power spectrum has been altered to mimic a 1keV dark matter particle. The CDM initial conditions were constrained to reproduce at z = 0 the correct local environment within which a "Local Group" (LG) of galaxies may form. Two significant differences between the two simulations are found. While in the CDM case a group of galaxies that resembles the real LG forms, the WDM run fails to reproduce a viable LG, instead forming a diffuse group which is still expanding at z = 0. This is surprising since, due to the suppression of small scale power in its power spectrum, WDM is naively expected to only affect the collapse of small haloes and not necessarily the dynamics on a scale of a group of galaxies. Furthermore the concentration of baryons in halo center's is greater in CDM than in WDM and the properties of the disks differ.

preprint2013arXiv

Counterrotating Stars in Simulated Galaxy Disks

Counterrotating stars in disk galaxies are a puzzling dynamical feature whose origin has been ascribed to either satellite accretion events or to disk instabilities triggered by deviations from axisymmetry. We use a cosmological simulation of the formation of a disk galaxy to show that counterrotating stellar disk components may arise naturally in hierarchically-clustering scenarios even in the absence of merging. The simulated disk galaxy consists of two coplanar, overlapping stellar components with opposite spins: an inner counterrotating bar-like structure made up mostly of old stars surrounded by an extended, rotationally-supported disk of younger stars. The opposite-spin components originate from material accreted from two distinct filamentary structures which at turn around, when their net spin is acquired, intersect delineating a "V"-like structure. Each filament torques the other in opposite directions; the filament that first drains into the galaxy forms the inner counterrotating bar, while material accreted from the other filament forms the outer disk. Mergers do not play a substantial role and most stars in the galaxy are formed in situ; only 9% of all stars are contributed by accretion events. The formation scenario we describe here implies a significant age difference between the co- and counterrotating components, which may be used to discriminate between competing scenarios for the origin of counterrotating stars in disk galaxies.

preprint2013arXiv

The clustering of galaxies at z~0.5 in the SDSS-III Data Release 9 BOSS-CMASS sample: a test for the LCDM cosmology

We present results on the clustering of 282,068 galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) sample of massive galaxies with redshifts 0.4<z<0.7 which is part of the Sloan Digital Sky Survey III project. Our results cover a large range of scales from ~0.5 to ~90 Mpc/h. We compare these estimates with the expectations of the flat LCDM cosmological model with parameters compatible with WMAP7 data. We use the MultiDark cosmological simulation together with a simple halo abundance matching technique, to estimate galaxy correlation functions, power spectra, abundance of subhaloes and galaxy biases. We find that the LCDM model gives a reasonable description to the observed correlation functions at z~0.5, which is a remarkably good agreement considering that the model, once matched to the observed abundance of BOSS galaxies, does not have any free parameters. However, we find a deviation (>~10%) in the correlation functions for scales less than ~1 Mpc/h and ~10-40 Mpc/h. A more realistic abundance matching model and better statistics from upcoming observations are needed to clarify the situation. We also estimate that about 12% of the "galaxies" in the abundance-matched sample are satellites inhabiting central haloes with mass M>~1e14 M_sun/h. Using the MultiDark simulation we also study the real space halo bias b(r) of the matched catalogue finding that b=2.00+/-0.07 at large scales, consistent with the one obtained using the measured BOSS projected correlation function. Furthermore, the linear large-scale bias depends on the number density n of the abundance-matched sample as b=-0.048-(0.594+/-0.02)*log(n/(h/Mpc)^3). Extrapolating these results to BAO scales we measure a scale-dependent damping of the acoustic signal produced by non-linear evolution that leads to ~2-4% dips at ~3 sigma level for wavenumbers k>~0.1 h/Mpc in the linear large-scale bias.

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

Dwarf Galaxies and the Cosmic Web

We use a cosmological simulation of the formation of the Local Group of Galaxies to identify a mechanism that enables the removal of baryons from low-mass halos without appealing to feedback or reionization. As the Local Group forms, matter bound to it develops a network of filaments and pancakes. This moving web of gas and dark matter drifts and sweeps a large volume, overtaking many halos in the process. The dark matter content of these halos is unaffected but their gas can be efficiently removed by ram-pressure. The loss of gas is especially pronounced in low-mass halos due to their lower binding energy and has a dramatic effect on the star formation history of affected systems. This "cosmic web stripping" may help to explain the scarcity of dwarf galaxies compared with the numerous low-mass halos expected in ΛCDM and the large diversity of star formation histories and morphologies characteristic of faint galaxies. Although our results are based on a single high-resolution simulation, it is likely that the hydrodynamical interaction of dwarf galaxies with the cosmic web is a crucial ingredient so far missing from galaxy formation models.

preprint2012arXiv

On the origin of the angular momentum properties of gas and dark matter in galactic halos and its implications

We perform a set of non-radiative hydrodynamical simulations of merging spherical halos in order to understand the angular momentum (AM) properties of the galactic halos seen in cosmological simulations. The universal shape of AM distributions seen in simulations is found to be generically produced as a result of mergers. The universal shape is such that it has an excess of low AM material and hence cannot explain the exponential structure of disk galaxies. A resolution to this is suggested by the spatial distribution of low AM material which is found to be in the centre and a conical region close to the axis of rotation. A mechanism that preferentially discards the material in the centre and prevents the material along the poles from falling onto the disc is proposed as a solution. We implement a simple geometric criteria for selective removal of low AM material and show that in order for 90% of halos to host exponential discs one has to reject at least 40% of material. Next, we explore the physical mechanisms responsible for distributing the AM within the halo during a merger. For dark matter there is an inside-out transfer of AM, whereas for gas there is an outside-in transfer, which is due to differences between collisionless and gas dynamics. We also explain the apparent high spin of dark matter halos undergoing mergers and show that a criteria stricter than what is currently used, would be required to detect such unrelaxed halos. Finally, we demonstrate that the misalignment of AM between gas and dark matter only occurs when the intrinsic spins of the merging halos are not aligned with the orbital AM of the system. The self-misalignment (orientation of AM when measured in radial shells not being constant), which could be the cause of warps and anomalous rotation in disks galaxies, also occurs under similar conditions.

preprint2012arXiv

The cosmic web and the orientation of angular momenta

We use a 64$h^{-1}$Mpc dark matter (DM) only cosmological simulation to examine the large scale orientation of haloes and substructures with respect the cosmic web. A web classification scheme based on the velocity shear tensor is used to assign to each halo in the simulation a web type: knot, filament, sheet or void. Using $\sim10^6$ haloes that span ~3 orders of magnitude in mass the orientation of the halo's spin and the orbital angular momentum of subhaloes with respect to the eigenvectors of the shear tensor is examined. We find that the orbital angular momentum of subhaloes tends to align with the intermediate eigenvector of the velocity shear tensor for all haloes in knots, filaments and sheets. This result indicates that the kinematics of substructures located deep within the virialized regions of a halo is determined by its infall which in turn is determined by the large scale velocity shear, a surprising result given the virilaized nature of haloes. The non-random nature of subhalo accretion is thus imprinted on the angular momentum measured at z = 0. We also find that haloes' spin axis is aligned with the third eigenvector of the velocity shear tensor in filaments and sheets: the halo spin axis points along filaments and lies in the plane of cosmic sheets.

preprint2012arXiv

The Galactic thin and thick disk

The ongoing large spectroscopic surveys of the Milky Way such as SEGUE and RAVE have enabled us to take a fresh look at the structure of the Galactic thin and thick disks, and how their structure fits within the framework of structure formation via hierarchical clustering. In this article I will summarize some recent results mainly based in the RAVE survey with respect to the structure of the Galactic disks, their possible origin as well as indications of substructure and asymmetries.

preprint2012arXiv

The Ninth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Baryon Oscillation Spectroscopic Survey

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

preprint2012arXiv

The velocity shear tensor: tracer of halo alignment

The alignment of DM halos and the surrounding large scale structure (LSS) is examined in the context of the cosmic web. Halo spin, shape and the orbital angular momentum of subhaloes is investigated relative to the LSS using the eigenvectors of the velocity shear tensor evaluated on a grid with a scale of 1 Mpc/h, deep within the non-linear regime. Knots, filaments, sheets and voids are associated with regions that are collapsing along 3, 2, 1 or 0 principal directions simultaneously. Each halo is tagged with a web classification (i.e. knot halo, filament halo, etc) according to the nature of the collapse at the halo's position. The full distribution of shear eigenvalues is found to be substantially different from that tagged to haloes, indicating that the observed velocity shear is significantly biased. We find that larger mass haloes live in regions where the shear is more isotropic, namely the expansion or collapse is more spherical. A correlation is found between the halo's shape and the eigenvectors of the shear tensor, with the longest (shortest) axis of the halo's shape being aligned with the slowest (fastest) collapsing eigenvector. This correlation is web independent, suggesting that the velocity shear is a fundamental tracer of the halo alignment. A similar result is found for the alignment of halo spin with the cosmic web. It has been shown that high mass haloes exhibit a spin flip with respect to the LSS: we find the mass at which this spin flip occurs is web dependent and not universal as suggested previously. Although weaker than haloes, subhalo orbits too exhibit an alignment with the LSS, providing a possible insight into the highly correlated co-rotation of the Milky Way's satellite system. The present study suggests that the velocity shear tensor constitutes the natural framework for studying the directional properties of the non-linear LSS and of halos and galaxies.

preprint2012arXiv

Tracing the origin of the Aquarius Stream: I

We present an abundance analysis of six member stars of the recently discovered Aquarius stream, in an attempt to ascertain whether this halo stream is real and, if so, to understand its origin. The mean metallicities of the six stars have a dispersion of only 0.10 dex, indicating that they are part of a chemically coherent structure. We then investigate whether the stream represents the debris of a disrupted dwarf galaxy or a disrupted globular cluster. The [Ni/Fe] - [Na/Fe] plane provides a good diagnostic: globular cluster stars and dwarf spheroidal galaxy stars are well separated in this plane, and the Aquarius stream stars lie unambiguously in the globular cluster region. The Aquarius stream stars also lie on the distinct [Na/Fe] - [O/Fe] and [Mg/Fe] - [Al/Fe] relations delineated by Galactic globular cluster stars. Spectroscopic parameters for the six Aquarius stars show that they are tightly confined to a 12 Gyr, [Fe/H] = -1.0, alpha-enhanced isochrone, consistent with their identification as globular cluster debris. We present evidence that the Aquarius stream may continue through the disk and out into the northern halo. Our results indicate a globular cluster origin for the Aquarius stream, and demonstrate the potential for chemical tagging to identify the origins of Galactic substructures.

preprint2011arXiv

A Sagittarius-Induced Origin for the Monoceros Ring

The Monoceros ring is a collection of stars in nearly-circular orbits at roughly 18 kpc from the Galactic center. It may have originated (i) as the response of the disc to perturbations excited by satellite companions or (ii) from the tidal debris of a disrupted dwarf galaxy. The metallicity of Monoceros stars differs from that of disc stars at comparable Galactocentric distances, an observation that disfavours the first scenario. On the other hand, circular orbits are difficult to accommodate in the tidal-disruption scenario, since it requires a satellite which at the time of disruption was itself in a nearly circular orbit. Such satellite could not have formed at the location of the ring and, given its low mass, dynamical friction is unlikely to have played a major role in its orbital evolution. We search cosmological simulations for low-mass satellites in nearly-circular orbits and find that they result, almost invariably, from orbital changes induced by collisions with more massive satellites: the radius of the circular orbit thus traces the galactocentric distance of the collision. Interestingly, the Sagittarius dwarf, one of the most luminous satellites of the Milky Way, is in a polar orbit that crosses the Galactic plane at roughly the same Galactocentric distance as Monoceros. We use idealized simulations to demonstrate that an encounter with Sagittarius might well have led to the circularization and subsequent tidal demise of the progenitor of the Monoceros ring.

preprint2011arXiv

Disc Heating: Comparing the Milky Way with Cosmological Simulations

We present the analysis of a suite of simulations run with different particle-and grid-based cosmological hydrodynamical codes and compare them with observational data of the Milky Way. This is the first study to make comparisons of properties of galaxies simulated with particle and grid-based codes. Our analysis indicates that there is broad agreement between these different modelling techniques. We study the velocity dispersion - age relation for disc stars at z=0 and find that four of the simulations are more consistent with observations by Holmberg et al. (2008) in which the stellar disc appears to undergo continual/secular heating. Two other simulations are in better agreement with the Quillen & Garnett (2001) observations that suggest a "saturation" in the heating profile for young stars in the disc. None of the simulations have thin discs as old as that of the Milky Way. We also analyse the kinematics of disc stars at the time of their birth for different epochs in the galaxies' evolution and find that in some simulations old stars are born cold within the disc and are subsequently heated, while other simulations possess old stellar populations which are born relatively hot. The models which are in better agreement with observations of the Milky Way's stellar disc undergo significantly lower minor-merger/assembly activity after the last major merger - i.e. once the disc has formed. All of the simulations are significantly "hotter" than the Milky Way disc; on top of the effects of mergers, we find a "floor" in the dispersion that is related to the underlying treatment of the heating and cooling of the interstellar medium, and the low density threshold which such codes use for star formation. This finding has important implications for all studies of disc heating that use hydrodynamical codes.

preprint2011arXiv

Kinematics in galactic tidal tails - A source for Hypervelocity stars?

[abridged] In this work we study in detail the kinematics of tidal debris stars to investigate the implications of the new scenario that the observed sample of Hypervelocity stars could partly originate from a dwarf-host galaxy collision. We use a suite of N-body simulations following the encounter of a satellite galaxy with its Milky Way-type host galaxy to gather statistics on the properties of stripped-off stars. We study especially the orbital energy distribution of this population. We quantify the typical pattern in angular and phase space formed by the debris stars. We further develop a simple stripping model predicting the kinematics of stripped-off stars. We show that the distribution of orbital energies in the tidal debris has a typical form which can be described quite accurately by a simple function. Based on this we develop a method to predict the energy distribution which allows us to evaluate the significance and the implications of high velocity stars in satellite tidal debris. Generally tidal collisions of satellite galaxy produce stars which escape into the intragalactic space even if the satellite itself is on a bound orbit. The main parameters determining the maximum energy kick a tidal debris star can get is the initial mass of the satellite and only to a lower extent its orbit. Main contributors to an unbound stellar population created in this way are massive satellites (M_sat > 10^9 M_sun). We thus expect intragalactic stars to have a metallicity higher than the surviving satellite population of the Milky Way. However, the probability that the observed HVS population is significantly contaminated by tidal debris stars appears small in the light of our results.

preprint2011arXiv

Metal-Poor Lithium-Rich Giants in the Radial Velocity Experiment Survey

We report the discovery of eight lithium-rich field giants found in a high resolution spectroscopic sample of over 700 metal-poor stars ([Fe/H]<-0.5) selected from the RAVE survey. The majority of the Li-rich giants in our sample are very metal-poor ([Fe/H]<-1.9), and have a Li abundance (in the form of 7Li), A(Li)=log(n(Li)/n(H))+12, between 2.30 and 3.63, well above the typical upper red giant branch limit, A(Li)<0.5, while two stars, with A(Li)~1.7-1.8, show similar lithium abundances to normal giants at the same gravity. We further included two metal-poor, Li-rich globular cluster giants in our sample, namely the previously discovered M3-IV101 and newly discovered (in this work) M68-A96. This comprises the largest sample of metal-poor Li-rich giants to date. We performed a detailed abundance analysis of all stars, finding that the majority our sample stars have elemental abundances similar to that of Li-normal halo giants. Although the evolutionary phase of each Li-rich giant cannot be definitively determined, the Li-rich phase is likely connected to extra mixing at the red giant branch bump or early asymptotic giant branch that triggers cool bottom processing in which the bottom of the outer convective envelope is connected to the H-burning shell in the star. The surface of a star becomes Li-enhanced as 7Be (which burns to 7Li) is transported to the stellar surface via the Cameron-Fowler mechanism. We discuss and discriminate among several models for the extra mixing that can cause Li-production, given the detailed abundances of the Li-rich giants in our sample.

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

SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

preprint2011arXiv

The Eighth Data Release of the Sloan Digital Sky Survey: First Data from SDSS-III

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

preprint2011arXiv

The MultiDark Database: Release of the Bolshoi and MultiDark Cosmological Simulations

We present the online MultiDark Database -- a Virtual Observatory-oriented, relational database for hosting various cosmological simulations. The data is accessible via an SQL (Structured Query Language) query interface, which also allows users to directly pose scientific questions, as shown in a number of examples in this paper. Further examples for the usage of the database are given in its extensive online documentation (www.multidark.org). The database is based on the same technology as the Millennium Database, a fact that will greatly facilitate the usage of both suites of cosmological simulations. The first release of the MultiDark Database hosts two 8.6 billion particle cosmological N-body simulations: the Bolshoi (250/h Mpc simulation box, 1/h kpc resolution) and MultiDark Run1 simulation (MDR1, or BigBolshoi, 1000/h Mpc simulation box, 7/h kpc resolution). The extraction methods for halos/subhalos from the raw simulation data, and how this data is structured in the database are explained in this paper. With the first data release, users get full access to halo/subhalo catalogs, various profiles of the halos at redshifts z=0-15, and raw dark matter data for one time-step of the Bolshoi and four time-steps of the MultiDark simulation. Later releases will also include galaxy mock catalogs and additional merging trees for both simulations as well as new large volume simulations with high resolution. This project is further proof of the viability to store and present complex data using relational database technology. We encourage other simulators to publish their results in a similar manner.

preprint2010arXiv

AstroGrid-D: Grid Technology for Astronomical Science

We present status and results of AstroGrid-D, a joint effort of astrophysicists and computer scientists to employ grid technology for scientific applications. AstroGrid-D provides access to a network of distributed machines with a set of commands as well as software interfaces. It allows simple use of computer and storage facilities and to schedule or monitor compute tasks and data management. It is based on the Globus Toolkit middleware (GT4). Chapter 1 describes the context which led to the demand for advanced software solutions in Astrophysics, and we state the goals of the project. We then present characteristic astrophysical applications that have been implemented on AstroGrid-D in chapter 2. We describe simulations of different complexity, compute-intensive calculations running on multiple sites, and advanced applications for specific scientific purposes, such as a connection to robotic telescopes. We can show from these examples how grid execution improves e.g. the scientific workflow. Chapter 3 explains the software tools and services that we adapted or newly developed. Section 3.1 is focused on the administrative aspects of the infrastructure, to manage users and monitor activity. Section 3.2 characterises the central components of our architecture: The AstroGrid-D information service to collect and store metadata, a file management system, the data management system, and a job manager for automatic submission of compute tasks. We summarise the successfully established infrastructure in chapter 4, concluding with our future plans to establish AstroGrid-D as a platform of modern e-Astronomy.

preprint2010arXiv

eROSITA on SRG

eROSITA (extended ROentgen Survey with an Imaging Telescope Array) is the core instrument on the Russian Spektrum-Roentgen-Gamma (SRG) mission which is scheduled for launch in late 2012. eROSITA is fully approved and funded by the German Space Agency DLR and the Max-Planck-Society.The design driving science is the detection of 50-100 thousands Clusters of Galaxies up to redshift z~1.3 in order to study the large scale structure in the Universe and test cosmological models, especially Dark Energy. This will be accomplished by an all-sky survey lasting for four years plus a phase of pointed observations. eROSITA consists of seven Wolter-I telescope modules, each equipped with 54 Wolter-I shells having an outer diameter of 360 mm. This would provide and effective area at 1.5 keV of ~1500 cm2 and an on axis PSF HEW of 15" which would provide an effective angular resolution of 25"-30". In the focus of each mirror module, a fast frame-store pn-CCD will provide a field of view of 1deg in diameter for an active FOV of ~0.83 deg2. At the time of writing the instrument development is currently in phase C/D.

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 preferred direction of infalling satellite galaxies in the Local Group

Using a high resolution DM simulation of the Local Group (LG), conducted within the framework of the Constrained Local UniversE Simulation (CLUES) project, we investigate the nature of how satellites of the MW and M31 are accreted. Satellites of these galaxies are accreted anisotropically, entering the virial radius of their hosts from specific "spots" with respect to the large scale structure. Furthermore, material which is tidally stripped from these satellites is also, at z=0, distributed anisotropically and is characterized by an ellipsoidal sub-volume embedded in the halo. The angular pattern created by the locus of satellite infall points and the projected stripped DM is investigated within a coordinate system determined by the location of the LG companion and the simulated Virgo cluster across concentric shells ranging from 0.1 to 5 r_vir. Remarkably, the principal axis of the ellipsoidal sub-volume shows a coherent alignment extending from well within the halo to a few r_vir. A spherical harmonics transform applied to the angular distributions confirms the visual impression: namely, the angular distributions of both the satellites entry points and stripped DM is dominated by the l=2 quadrupole term, whose major principal axis is aligned across the shells considered. It follows that the structure of the outer (r>0.5 r_vir) regions of the main halos is closely related to the cosmic web within which it is embedded. Given the hypothesis that a large fraction of the stellar halo of the Milky Way has been accreted from satellites, the present results can be applied to the stellar halo of the MW and M31. We predict that the remnants of tidally stripped satellites should be embedded in streams of material composed of dark matter and stars. The present results can therefore shed light on the existence of satellites embedded within larger streams of matter, such as the Segue 2 satellite.

preprint2009arXiv

Galaxy-Induced Transformation of Dark Matter Halos

We use N-body/gasdynamical LambdaCDM cosmological simulations to examine the effect of the assembly of a central galaxy on the shape and mass profile of its dark halo. Two series of simulations are compared; one that follows only the evolution of the dark matter component and a second one where a baryonic component is added. These simulations include radiative cooling but neglect star formation and feedback, leading most baryons to collect at the halo center in a disk which is too small and too massive when compared with typical spiral. This unrealistic model allows us, nevertheless, to gauge the maximum effect that galaxies may have in transforming their dark halos. We find that the shape of the halo becomes more axisymmetric: halos are transformed from triaxial into essentially oblate systems, with well-aligned isopotential contours of roughly constant flattening (c/a ~ 0.85). Halos always contract as a result of galaxy assembly, but the effect is substantially less pronounced than predicted by the "adiabatic contraction" hypothesis. The reduced contraction helps to reconcile LambdaCDM halos with constraints on the dark matter content inside the solar circle and should alleviate the long-standing difficulty of matching simultaneously the scaling properties of galaxy disks and the luminosity function. The halo contraction is also less pronounced than found in earlier simulations, a disagreement that suggests that halo contraction is not solely a function of the initial and final distribution of baryons. Not only how much baryonic mass has been deposited at the center of a halo matters, but also the mode of its deposition. It might prove impossible to predict the halo response without a detailed understanding of a galaxy's assembly history. (Abriged)

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

The Angular Momentum Problem in Cosmological Simulations of Disk Galaxy Formation

We conduct a systematic study of the angular momentum problem in numerical simulations of disk galaxy formation. We investigate the role of numerical resolution using a semi-cosmological setup which combines an efficient use of the number of particles in an isolated halo while preserving the hierarchical build-up of the disk through the merging of clumps. We perform the same simulation varying the resolution over 4 orders of magnitude. Independent on the level of resolution, the loss of angular momentum stays the same and can be tied to dynamical friction during the build-up phase. This is confirmed in a cosmological simulation. We also perform simulations including star formation and star formation and supernova feedback. While the former has no influence on the angular momentum problem, the latter reduces the loss to a level potentially in agreement with observations. This is achieved through a suppression of early star formation and therefore the formation of a large, slowly rotating bulge. We conclude that feedback is a critical component to achieve realistic disk galaxies in cosmological simulations. Numerical resolution is important, but by itself not capable of solving the angular momentum problem.

preprint2008arXiv

Multi-Dimensional Density Estimation and Phase Space Structure Of Dark Matter Halos

We present a method to numerically estimate the densities of a discretely sampled data based on binary space partitioning tree. We start with a root node containing all the particles and then recursively divide each node into two nodes each containing roughly equal number of particles,until each of the nodes contains only one particle. The volume of such a leaf node provides an estimate of the local density. We implement an entropy-based node splitting criterion that results in a significant improvement in the estimation of densities compared to earlier work. The method is completely metric free and can be applied to arbitrary number of dimensions. We apply this method to determine the phase space densities of dark matter halos obtained from cosmological N-body simulations. We find that contrary to earlier studies, the volume distribution function $v(f)$ of phase space density $f$ does not have a constant slope but rather a small hump at high phase space densities. We demonstrate that a model in which a halo is made up by a superposition of Hernquist spheres is not capable in explaining the shape of $v(f)$ vs $f$ relation, whereas a model which takes into account the contribution of the main halo separately roughly reproduces the behavior as seen in simulations. The use of the presented method is not limited to calculation of phase space densities, but can be used as a general-purpose data-mining tool and due to its speed and accuracy it is ideally suited for analysis of large multidimensional data sets.

preprint2000arXiv

Dark Halo and Disk Galaxy Scaling Laws in Hierarchical Universes

We use cosmological N-body/gasdynamical simulations that include star formation and feedback to examine the proposal that scaling laws between the total luminosity, rotation speed, and angular momentum of disk galaxies reflect analogous correlations between the structural parameters of their surrounding dark matter halos. The numerical experiments follow the formation of galaxy-sized halos in two Cold Dark Matter dominated universes: the standard Omega=1 CDM scenario and the currently popular LCDM model. We find that the slope and scatter of the I-band Tully-Fisher relation are well reproduced in the simulations, although not, as proposed in recent work, as a result of the cosmological equivalence between halo mass and circular velocity: large systematic variations in the fraction of baryons that collapse to form galaxies and in the ratio between halo and disk circular velocities are observed in our numerical experiments. The Tully-Fisher slope and scatter are recovered in this model as a direct result of the dynamical response of the halo to the assembly of the luminous component of the galaxy. We conclude that models that neglect the self-gravity of the disk and its influence on the detailed structure of the halo cannot be used to derive meaningful estimates of the scatter or slope of the Tully-Fisher relation. Our models fail, however, to match the zero-point of the Tully-Fisher relation, as well as that of the relation linking disk rotation speed and angular momentum. These failures can be traced, respectively, to the excessive central concentration of dark halos formed in the Cold Dark Matter cosmogonies we explore and to the formation of galaxy disks as the final outcome of a sequence of merger events. (abridged)

preprint1998arXiv

The cosmological origin of the Tully-Fisher relation

We use high-resolution cosmological simulations that include the effects of gasdynamics and star formation to investigate the origin of the Tully-Fisher relation in the standard Cold Dark Matter cosmogony. Luminosities are computed for each model galaxy using their full star formation histories and the latest spectrophotometric models. We find that at z=0 the stellar mass of model galaxies is proportional to the total baryonic mass within the virial radius of their surrounding halos. Circular velocity then correlates tightly with the total luminosity of the galaxy, reflecting the equivalence between mass and circular velocity of systems identified in a cosmological context. The slope of the relation steepens slightly from the red to the blue bandpasses, and is in fairly good agreement with observations. Its scatter is small, decreasing from \~0.45 mag in the U-band to ~0.34 mag in the K-band. The particular cosmological model we explore here seems unable to account for the zero-point of the correlation. Model galaxies are too faint at z=0 (by about two magnitudes) if the circular velocity at the edge of the luminous galaxy is used as an estimator of the rotation speed. The Tully-Fisher relation is brighter in the past, by about ~0.7 magnitudes in the B-band at z=1, at odds with recent observations of z~1 galaxies. We conclude that the slope and tightness of the Tully-Fisher relation can be naturally explained in hierarchical models but that its normalization and evolution depend strongly on the star formation algorithm chosen and on the cosmological parameters that determine the universal baryon fraction and the time of assembly of galaxies of different mass.

preprint1996arXiv

A Comparison of X-ray and Strong Lensing Properties of Simulated X-ray Clusters

We use gas-dynamical simulations of galaxy clusters to compare their X-ray and strong lensing properties. Special emphasis is laid on mass estimates. The cluster masses range between 6 x 10^14 solar masses and 4 x 10^15 solar masses, and they are examined at redshifts between 1 and 0. We compute the X-ray emission of the intracluster gas by thermal bremsstrahlung, add background contamination, and mimic imaging and spectral observations with current X-ray telescopes. Although the beta model routinely provides excellent fits to the X-ray emission profiles, the derived masses are typically biased low because of the restricted range of radii within which the fit can be done. For beta values of ~ 2/3, which is the average in our numerically simulated sample, the mass is typically underestimated by ~ 40 per cent. The masses of clusters which exhibit pronounced substructure are often substantially underestimated. We suggest that the ratio between peak temperature and emission-weighted average cluster temperature may provide a good indicator for ongoing merging and, therefore, for unreliable mass estimates. X-ray mass estimates are substantially improved if we fit a King density profile rather than the beta model to the X-ray emission, thereby dropping the degree of freedom associated with beta. Clusters selected for their strong lensing properties are typically dynamically more active than typical clusters. Bulk flows in the intracluster gas contain a larger than average fraction of the internal energy of the gas in such objects, hence the measured gas temperatures are biased low. The bulk of the optical depth for arc formation is contributed by clusters with intermediate rather than high X-ray luminosity. Arcs occur predominantly in clusters which exhibit substructure and are not in an equilibrium state. Finally we explain why the

preprint1996arXiv

Non-equilibrium effects on line-of-sight size estimates of QSO absorption systems

Estimates of the linear extent of heavy-element absorption systems along the line-of-sight to a QSO often assume that the cloud is photoionized and that the temperature takes the equilibrium value where photo-heating balances line cooling. We show that rather small deviations from this photoionization equilibrium temperature caused by additional heating processes will lead to an overestimate of the neutral hydrogen fraction and thus to an underestimate of the thickness of the absorber by about two orders of magnitude. Such temperature deviations are indicated both by observations and numerical simulations. This interpretation reconciles the discrepancy between the rather small extent of heavy-element-absorption systems parallel to the line-of-sight obtained from a standard photoionization analysis and the much larger transverse sizes estimates inferred from the observation of common absorption in the spectra of close quasar pairs.

preprint1996arXiv

QSO Metal Absorption Systems at High Redshift and the Signature of Hierarchical Galaxy Formation

In hierarchical cosmogonies at redshift three the matter content of a typical present-day galaxy is dispersed over several individual clumps embedded in sheet-like structures, often aligned along filaments. We have used hydrodynamical simulations to investigate the spatial distribution and absorption properties of metal enriched gas in such regions of ongoing galaxy formation. Strong CII and SiIV absorption occurs for lines-of-sight passing the densest regions close to the center of the protogalactic clumps, while CIV is a good tracer of the prominent filamentary structures and OVI becomes the strongest absorption feature for lines-of-sight through low-density regions far away from fully-collapsed objects. Observed column density ratios can be well reproduced if a mean metallicity [Z/H] = -2.5, relative abundances as found in metal-poor stars, a UV background with intensity J_{-22} = 3 at the Lyman limit, and either a power law spectrum (alpha = -1.5) or the spectral shape proposed by Haardt & Madau (1996) are assumed. The observed scatter in [C/H] is about a magnitude larger than that in the simulations suggesting inhomogeneous metal distribution. The detailed small and large scale kinematics observed in CIV and HI are also well reproduced.

preprint1996arXiv

Two--body heating in numerical galaxy formation experiments

We show that discreteness effects related to classical two-body relaxation produce spurious heating of the gaseous component in numerical simulations of galaxy formation. A simple analytic model demonstrates that this artificial heating will dominate radiative cooling in any simulation where the mass of an individual dark matter particle exceeds a certain critical value. This maximum mass depends only on the cooling function of the gas, on the fraction of the material in gaseous form, and (weakly) on typical temperatures in the gas. It is comparable to, or smaller than, the dark matter particle masses employed in most published simulations of cosmological hydrodynamics and galaxy formation. Any simulation which violates this constraint will be unable to follow cooling flows, although catastrophic cooling of gas may still occur in regions with very short cooling times. We use a series of N--body/smoothed particle hydrodynamics simulations to explore this effect. In simulations which neglect radiative cooling, two--body heating causes a gradual expansion of the gas component. When radiative effects are included, we find that gas cooling is almost completely suppressed for dark matter particle masses above our limit. Although our test simulations use smoothed particle hydrodynamics, similar effects, and a similar critical mass, are expected in any simulation where the dark matter is represented by discrete particles.

preprint1995arXiv

Galaxy Formation and Large-Scale Bias

We outline a simple approach to understanding the physical origin of bias in the distribution of galaxies relative to that of dark matter. The first step is to specify how collapsed, virialized halos of dark matter trace the overall matter distribution. The next step is to make a connection between halos and the luminous galaxies we observe. We appeal to the results of semi-analytic models of galaxy formation that are tuned to fit the observed luminosity functions of local groups and clusters. We have also used a high-resolution N-body simulation of a cold dark matter (CDM) universe to study the bias relation in more detail. The differences between the galaxy and dark matter distributions are quantified using a number of different clustering statistics. We arrive at the following general conclusions: 1) A comparison of the galaxy and dark matter density fields shows that linear biasing is a good description on large scales. 2) The bias factor b depends on the shape and normalization of the power spectrum. The lower the normalization, the larger the bias. More bias is obtained for spectra with more power on large scales. For "realistic" models, b ranges from 1 to 2.5. 3) Galaxies of different luminosity or morphology have different bias factors. 4) The scale dependence of the bias factor is weak.

preprint1995arXiv

GRAPESPH: Cosmological SPH simulations with the special purpose hardware GRAPE

A combined N--body/SPH code is presented which benefits from the high speed of the special purpose hardware GRAPE (GRAvity PipE). Besides gravitational forces, GRAPE also returns the list of neighbours and can, therefore, be used to speed up the hydrodynamical part, too. After the interaction list has been passed, density, pressure forces, propagation and interpolation of particles etc. are calculated on the front end, a 50 MHz SUN SPARC 10. In order to combine SPH and GRAPE, possible limitations due to the hardware design of GRAPE are carefully analyzed and modifications compared to current SPH codes are discussed. The resulting code, GRAPESPH, is similarly flexible as TREESPH. It seems especially well suited to investigate the formation of individual objects in a large scale structure environment, like eg. galaxies or clusters. The total performance is at least half as good as treesph on a CRAY and for most applications it seems to be even better. The CPU time per time step is only slightly dependent on the clustering state. GRAPESPH, therefore, provides a very attractive alternative to the use of supercomputers in cosmology.

preprint1995arXiv

Mergers and the Formation of Disk Galaxies in Hierarchically Clustering Universes

New simulations are presented which investigate the formation of smaller groups of galaxies in a CDM like universe. The simulations follow the evolution of dark matter and gas and are performed with an SPH code adapted for the special-purpose hardware GRAPE. The mass resolution in the baryonic component is 5 10^6 Msol, the spatial resolution is 1.5 kpc. A sufficiently large volume is simulated in order to accurately include the tidal field and the mass inflow. The mass resolution enables us to resolve the galaxies formed in the simulation. Individual properties of the galaxies, each consisting of several thousand particles, can also be investigated. It turns out, that in the center of the dark haloes the gas accumulates to form a rotationally supported disk. However, in comparison with observations, disks which form in numerical simulations are too concentrated: Due to the merging of substructures, angular momentum is transported from the gas to the dark halo, and far too much gas is accumulated at very small radii. It is argued that changing the cosmogony (e.g. to low Omega) or adding photoionisation due to an external UV background do not provide an easy fix to the angular momentum problem. Feedback processes caused by supernovae or stellar winds seem to be the most likely solution to this problem.

preprint1994arXiv

On the spin parameter of dark-matter haloes

The study by White (1984) on the growth of angular momentum in dark haloes is extended towards a more detailed investigation of the spin parameter $λ\equiv L\sqrt{E}/{G M^{2.5}}$. Starting from the Zel'dovich approximation to structure formation, a dark halo is approximated by a homogeneous ellipsoid with the inertial tensor of the (highly irregular) Lagrangian region $Υ$ from which the dark halo forms. Within this approximation, an expression for the spin parameter can be derived, which depends on the geometry of $Υ$, the cosmological density parameter $Ω_0$, the overdensity of the dark halo, and the tidal torque exerted on it. For Gaussian random fields, this expression can be evaluated statistically. As a result, we derive a probability distribution of the spin parameter which gives $λ\simeq0.07^{+0.04}_{-0.05}$, consistent with numerical investigations. This probability distribution steeply rises with increasing spin parameter, reaching its maximum at $λ\simeq0.025$. The 10 (50,90) percentile values are $λ=0.02$ (0.05,0.11, respectively). There is a weak anticorrelation of the spin parameter with the peak height $ν$ of the density fluctuation field $λ\propto ν^{-0.29}$. The dependence on $Ω_0$ and the variance $σ$ of the density-contrast field is very weak; there is only a marginal tendency for the spin parameter to be slightly larger for late-forming objects in an open universe. Due to the weak dependence on $σ$, our results should be quite generally applicable and independent on

preprint1994arXiv

The Formation of Disk Galaxies in a Cosmological Context: Structure and Kinematics

We present results concerning the internal structure and kinematics of disk galaxies formed in cosmologically motivated simulations. The calculations include dark matter, gas dynamics, radiative cooling, star formation, supernova feedback and metal enrichment. The initial model is a rigidly rotating overdense sphere with a mass of about 8 10^11 Msol which is perturbed by small scale fluctuations according to a biased CDM power spectrum. Converging, Jeans unstable and rapidly cooling regions are allowed to form stars. Via supernovae, metal enriched gas is returned to the interstellar medium. {}From these initial conditions a galaxy forms which shows the main properties of spiral galaxies: a rotationally supported exponential disk which consists of young stars with about solar metallicity, a slowly rotating halo of old metal poor stars, a bulge of old metal rich stars and a slowly rotating extended halo of dark matter. Bulge, stellar and dark halo are supported by an anisotropic velocity dispersion and have a de Vaucouleurs surface density profile. The flattening of the dark and stellar halo is too large to be explained by rotation only. Whether the flattening of the bulge is caused by an anisotropic velocity dispersion or by its rotation cannot be answered, because of the limited numerical resolution due to gravitational softening. The velocity dispersion and the thickness of the stellar disk increase with the age of the stars. Considering only the young stellar component, the disk is cold (sigma=20 km/sec) and thin (z <1 kpc). The dynamical formation process ends after about 4\,Gyr, when