Source author record

Ralph Schönrich

Ralph Schönrich 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

15works
4topics
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

15 published item(s)

preprint2022arXiv

Age distribution of stars in boxy/peanut/X-shaped bulges formed without bar buckling

Some barred galaxies, including the Milky Way, host a boxy/peanut/X-shaped bulge (BPX-shaped bulge). Previous studies suggested that the BPX-shaped bulge can either be developed by bar buckling or by vertical inner Lindblad resonance (vILR) heating without buckling. In this paper, we study the observable consequence of an BPX-shaped bulge built up quickly after bar formation via vILR heating without buckling, using an N-body/hydrodynamics simulation of an isolated Milky Way-like galaxy. We found that the BPX-shaped bulge is dominated by stars born prior to bar formation. This is because the bar suppresses star formation, except for the nuclear stellar disc (NSD) region and its tips. The stars formed near the bar ends have higher Jacobi energy, and when these stars lose their angular momentum, their non-circular energy increases to conserve Jacobi energy. This prevents them from reaching the vILR to be heated to the BPX region. By contrast, the NSD forms after the bar formation. From this simulation and general considerations, we expect that the age distributions of the NSD and BPX-shaped bulge formed without bar buckling do not overlap each other. Then, the transition age between these components betrays the formation time of the bar, and is testable in future observations of the Milky Way and extra-galactic barred galaxies.

preprint2022arXiv

Oscillating dynamical friction on galactic bars by trapped dark matter

The dynamic evolution of galactic bars in standard $Λ$CDM models is dominated by angular momentum loss to the dark matter haloes via dynamical friction. Traditional approximations to dynamical friction are formulated using the linearized collisionless Boltzmann equation and have been shown to be valid in the fast limit, i.e. for rapidly slowing bars. However, the linear assumption breaks down within a few dynamical periods for typical slowly evolving bars, which trap a significant amount of disc stars and dark matter in resonances. Recent observations of the Galactic bar imply this slow regime at the main bar resonances. We formulate the time-dependent dynamical friction in the slow limit and explore its mechanism in the general slow regime with test-particle simulations. Here, angular momentum exchange is dominated by resonantly trapped orbits which slowly librate around the resonances. In typical equilibrium haloes, the initial phase-space density within the trapped zone is higher at lower angular momentum. Since the libration frequency falls towards the separatrix, this density contrast winds up into a phase-space spiral, resulting in a dynamical friction that oscillates with $\sim$Gyr periods and damps over secular timescales. We quantify the long-term behaviour of this torque with secular perturbation theory, and predict two observable consequences: i) The phase-space spirals may be detectable in the stellar disc where the number of windings encodes the age of the bar. ii) The torque causes weak oscillations in the bar's pattern speed, overlaying the overall slowdown -- while not discussed, this feature is visible in previous simulations.

preprint2022arXiv

The K2 Galactic Archaeology Program Data Release 3: Age-abundance patterns in C1-C8 and C10-C18

We present the third and final data release of the K2 Galactic Archaeology Program (K2 GAP) for Campaigns C1-C8 and C10-C18. We provide asteroseismic radius and mass coefficients, $κ_R$ and $κ_M$, for $\sim 19,000$ red giant stars, which translate directly to radius and mass given a temperature. As such, K2 GAP DR3 represents the largest asteroseismic sample in the literature to date. K2 GAP DR3 stellar parameters are calibrated to be on an absolute parallactic scale based on Gaia DR2, with red giant branch and red clump evolutionary state classifications provided via a machine-learning approach. Combining these stellar parameters with GALAH DR3 spectroscopy, we determine asteroseismic ages with precisions of $\sim 20-30\%$ and compare age-abundance relations to Galactic chemical evolution models among both low- and high-$α$ populations for $α$, light, iron-peak, and neutron-capture elements. We confirm recent indications in the literature of both increased Ba production at late Galactic times, as well as significant contribution to r-process enrichment from prompt sources associated with, e.g., core-collapse supernovae. With an eye toward other Galactic archaeology applications, we characterize K2 GAP DR3 uncertainties and completeness using injection tests, suggesting K2 GAP DR3 is largely unbiased in mass/age and with uncertainties of $2.9\%\,(\rm{stat.})\,\pm0.1\%\,(\rm{syst.})$ & $6.7\%\,(\rm{stat.})\,\pm0.3\%\,(\rm{syst.})$ in $κ_R$ & $κ_M$ for red giant branch stars and $4.7\%\,(\rm{stat.})\,\pm0.3\%\,(\rm{syst.})$ & $11\%\,(\rm{stat.})\,\pm0.9\%\,(\rm{syst.})$ for red clump stars. We also identify percent-level asteroseismic systematics, which are likely related to the time baseline of the underlying data, and which therefore should be considered in TESS asteroseismic analysis.

preprint2021arXiv

Capture of interstellar objects II: by the Solar system

Capture of interstellar objects (ISOs) into the Solar system is dominated by ISOs with asymptotic incoming speeds $v_\infty<4\,$km\,s$^{-1}$. The capture rate is proportional to the ISO phase-space density in the Solar vicinity and does not vary along the Sun's Galactic orbit, i.e.\ is not enhanced during a passage through a cloud of ISOs (in contrast to previous suggestions). Most bound orbits crossing those of Jupiter and Saturn are fully mixed with unbound phase space, implying that they hold the same ISO phase-space density. Assuming an interstellar number density $n_{iso}\sim0.1\,$au$^{-3}$, we estimate that in 1000 years the planets capture $\sim2$ ISOs (while $\sim17$ fall into the Sun), resulting in a population of $\sim8$ captured ISOs within 5\,au of the Sun at any time, less than the number of visiting ISOs passing through the same volume on hyperbolic orbits. In terms of phase-space volume, capture onto and ejection from the Solar system are equal, such that on average ISOs will not remain captive at $a\lesssim2000\,$au for extensive periods.

preprint2016arXiv

Age velocity dispersion relations and heating histories in disc galaxies

We analyse the heating of stellar discs by non axisymmetric structures and giant molecular clouds (GMCs) in N-body simulations of growing disc galaxies. The analysis resolves long-standing discrepancies between models and data by demonstrating the importance of distinguishing between measured age-velocity dispersion relations (AVRs) and the heating histories of the stars that make up the AVR. We fit both AVRs and heating histories with formulae proportional to t^beta and determine the exponents beta_R and beta_z derived from in-plane and vertical AVRs and ~beta_R and ~beta_z from heating histories. Values of beta_z are in almost all simulations larger than values of ~beta_z, whereas values of beta_R are similar to or mildly larger than values of ~beta_R. Moreover, values of beta_z (~beta_z) are generally larger than values of beta_R (~beta_R). The dominant cause of these relations is the decline over the life of the disc in importance of GMCs as heating agents relative to spiral structure and the bar. We examine how age errors and biases in solar neighbourhood surveys influence the measured AVR: they tend to decrease beta values by smearing out ages and thus measured dispersions. We compare AVRs and velocity ellipsoid shapes sigma_z/sigma_R from simulations to Solar neighbourhood data. We conclude that for the expected disc mass and dark halo structure, combined GMC and spiral/bar heating can explain the AVR of the Galactic thin disc. Strong departures of the disc mass or the dark halo structure from expectation spoil fits to the data.

preprint2016arXiv

The quiescent phase of galactic disc growth

We perform a series of controlled N-body simulations of growing disc galaxies within non-growing, live dark matter haloes of varying mass and concentration. Our initial conditions include either a low-mass disc or a compact bulge. New stellar particles are continuously added on near-circular orbits to the existing disc, so spiral structure is continuously excited. To study the effect of combined spiral and giant molecular cloud (GMC) heating on the discs we introduce massive, short-lived particles that sample a GMC mass function. An isothermal gas component is introduced for a subset of the models. We perform a resolution study and vary parameters governing the GMC population, the histories of star formation and radial scale growth. Models with GMCs and standard values for the disc mass and halo density provide the right level of self-gravity to explain the age velocity dispersion relation of the Solar neighbourhood (Snhd). GMC heating generates remarkably exponential vertical profiles with scaleheights that are radially constant and agree with observations of galactic thin discs. GMCs are also capable of significantly delaying bar formation. The amount of spiral induced radial migration agrees with what is required for the metallicity distribution of the Snhd. However, in our standard models the outward migrating populations are not hot enough vertically to create thick discs. Thick discs can form in models with high baryon fractions, but the corresponding bars are too long, the young stellar populations too hot and the discs flare considerably.

preprint2015arXiv

Kinematic Detection of the Galactic Nuclear Disc

We report the detection of the Galactic nuclear disc in line-of-sight kinematics of stars, measured with infrared spectroscopy from APOGEE. This stellar component of the nuclear disc has an extent and rotation velocity V ~ 120kms comparable to the gas disc in the central molecular zone. The current data suggest that this disc is kinematically cool and has a small vertical extent of order 50pc. The stellar kinematics suggest a truncation radius/steep decline of the stellar disc at a galactocentric radius R ~ 150pc, and provide tentative evidence for an overdensity at the position of the ring found in the molecular gas disc.

preprint2015arXiv

Origin of the high v_los feature in the Galactic bar

We analyse a controlled N-body + smoothed particle hydrodynamics simulation of a growing disc galaxy within a non-growing, live dark halo. The disc is continuously fed with gas and star particles on near-circular orbits and develops a bar comparable in size to the one of the Milky Way (MW). We extract line of sight velocity v_los distributions from the model and compare it to data recently obtained from the APOGEE survey which show distinct high velocity features around v_los ~ 200 km/s. With an APOGEE like selection function, but without any scaling nor adjustment, we find v_los distributions very similar to those in APOGEE. The stars that make up the high v_los features at positive longitudes l are preferentially young bar stars (age <~ 2-3 Gyr) which move away from us along the rear side of the bar. At negative l, we find the corresponding low v_los feature from stars moving towards us. At l>10 degrees the highest v_los stars are a mixture of bar and background disc stars which complicates the interpretation of observations. The main peak in v_los is dominated by fore- and background stars. At a given time, ~40-50 per cent of high v_los stars occupy x_1 like orbits, but a significant fraction are on more complex orbits. The observed feature is likely due to a population of dynamically cool, young stars formed from gas just outside the bar and subsequently captured by the growing bar. The high v_los features disappear at high latitudes |b|>~2 degrees which explains the non-detection of such features in other surveys.

preprint2015arXiv

Rings and Radial Waves in the Disk of the Milky Way

We show that in the anticenter region, between Galactic longitudes of $110^\circ<l<229^\circ$, there is an oscillating asymmetry in the main sequence star counts on either side of the Galactic plane using data from the Sloan Digital Sky Survey. This asymmetry oscillates from more stars in the north at distances of about 2 kpc from the Sun to more stars in the south at 4-6 kpc from the Sun to more stars in the north at distances of 8-10 kpc from the Sun. We also see evidence that there are more stars in the south at distances of 12-16 kpc from the Sun. The three more distant asymmetries form roughly concentric rings around the Galactic center, opening in the direction of the Milky Way's spiral arms. The northern ring, 9 kpc from the Sun, is easily identified with the previously discovered Monoceros Ring. Parts of the southern ring at 14 kpc from the Sun (which we call the TriAnd Ring) have previously been identified as related to the Monoceros Ring and others have been called the Triangulum Andromeda Overdensity. The two nearer oscillations are approximated by a toy model in which the disk plane is offset by of the order 100 pc up and then down at different radii. We also show that the disk is not azimuthally symmetric around the Galactic anticenter and that there could be a correspondence between our observed oscillations and the spiral structure of the Galaxy. Our observations suggest that the TriAnd and Monoceros Rings (which extend to at least 25 kpc from the Galactic center) are primarily the result of disk oscillations.

preprint2014arXiv

Fundamental stellar parameters and metallicities from Bayesian spectroscopy. Application to low- and high-resolution spectra

We present a unified framework to derive fundamental stellar parameters by combining all available observational and theoretical information for a star. The algorithm relies on the method of Bayesian inference, which for the first time directly integrates the spectroscopic analysis pipeline based on the global spectrum synthesis and allows for comprehensive and objective error calculations given the priors. Arbitrary input datasets can be included into our analysis and other stellar quantities, in addition to stellar age, effective temperature, surface gravity, and metallicity, can be computed on demand. We lay out the mathematical framework of the method and apply it to several observational datasets, including high- and low-resolution spectra (UVES, NARVAL, HARPS, SDSS/SEGUE). We find that simpler approximations for the spectroscopic PDF, which are inherent to past Bayesian approaches, lead to deviations of several standard deviations and unreliable errors on the same data. By its flexibility and the simultaneous analysis of multiple independent measurements for a star, it will be ideal to analyse and cross-calibrate the large ongoing and forthcoming surveys, like Gaia-ESO, SDSS, Gaia and LSST.

preprint2014arXiv

Testing the Asteroseismic Mass Scale Using Metal-Poor Stars Characterized with APOGEE and Kepler

Fundamental stellar properties, such as mass, radius, and age, can be inferred using asteroseismology. Cool stars with convective envelopes have turbulent motions that can stochastically drive and damp pulsations. The properties of the oscillation frequency power spectrum can be tied to mass and radius through solar-scaled asteroseismic relations. Stellar properties derived using these scaling relations need verification over a range of metallicities. Because the age and mass of halo stars are well-constrained by astrophysical priors, they provide an independent, empirical check on asteroseismic mass estimates in the low-metallicity regime. We identify nine metal-poor red giants (including six stars that are kinematically associated with the halo) from a sample observed by both the Kepler space telescope and the Sloan Digital Sky Survey-III APOGEE spectroscopic survey. We compare masses inferred using asteroseismology to those expected for halo and thick-disk stars. Although our sample is small, standard scaling relations, combined with asteroseismic parameters from the APOKASC Catalog, produce masses that are systematically higher (<ΔM>=0.17+/-0.05 Msun) than astrophysical expectations. The magnitude of the mass discrepancy is reduced by known theoretical corrections to the measured large frequency separation scaling relationship. Using alternative methods for measuring asteroseismic parameters induces systematic shifts at the 0.04 Msun level. We also compare published asteroseismic analyses with scaling relationship masses to examine the impact of using the frequency of maximum power as a constraint. Upcoming APOKASC observations will provide a larger sample of ~100 metal-poor stars, important for detailed asteroseismic characterization of Galactic stellar populations.

preprint2013arXiv

A new calibration for the Blue Horizontal Branch

We suggest a simple analytic approximation for magnitudes and hence distances of Blue Horizontal Branch (BHB) stars in Sloan colours. Precedent formulations do not offer a simple closed formula, nor do they cover the full dependences, e.g. on metallicity. Further, using BHB star samples from the Sloan Digital Sky Survey (SDSS) we validate our distance calibration directly on field stars instead of globular clusters and assess the performance of other available distance calibrations. Our method to statistically measure distances is sufficiently accurate to measure the colour and metallicity dependence on our sample and can be applied to other sets of stars or filters.

preprint2013arXiv

Rotational signature of the Milky Way stellar halo

We measure the rotation of the Milky Way stellar halo on two samples of Blue Horizontal Branch (BHB) field halo stars from the Sloan Digital Sky Survey (SDSS) with four different methods. The two samples comprise 1582 and 2563 stars respectively and reach out to ~50 kpc in galactocentric distance. Two of the methods to measure rotation rely exclusively on line-of-sight velocities, namely the popular double power-law model and a direct estimate of the de-projected l.o.s. velocity. The other two techniques use the full 3D motions: the radial velocity based rotation estimator of Schönrich, Binney & Asplund (2012) and a simple 3D azimuthal velocity mean. In this context we a) critique the popular model and b) assess the reliability of the estimators. All four methods agree on a weakly prograde or non-rotating halo. Further, we observe no duality in the rotation of sub-samples with different metallicities or at different radii. We trace the rotation gradient across metallicity measured by Deason et al. (2011) on a similar sample of BHB stars back to the inclusion of regions in the apparent magnitude-surface gravity plane known to be contaminated. In the spectroscopically selected sample of Xue et al. (2011), we flag ~500 hot metal-poor stars for their peculiar kinematics w.r.t. to both their cooler metal-poor counter-parts and to the metal-rich stars in the same sample. They show a seemingly retrograde behaviour in line-of-sight velocities, which is not confirmed by the 3D estimators. Their anomalous vertical motion hints at either a pipeline problem or a stream-like component rather than a smooth retrograde population.

preprint2012arXiv

Galactic Rotation and Solar Motion from Stellar Kinematics

I present three methods to determine the distance to the Galactic centre R0, the solar azimuthal velocity in the Galactic rest frame Vg,\odot and hence the local circular speed Vc at R0. These simple, model-independent strategies reduce the set of assumptions to near axisymmetry of the disc and are designed for kinematically hot stars, which are less affected by spiral arms and other effects. The first two methods use the position-dependent rotational streaming in the heliocentric radial velocities U. The resulting rotation estimate θ from U velocities does not depend on Vg,\odot. The first approach compares this with rotation from the galactic azimuthal velocities to constrain Vg,\odot at an assumed R0. Both Vg,\odot and R0 can be determined using the proper motion of Sgr A\ast as a second constraint. The second strategy makes use of θ being roughly proportional to R0. Therefore a wrong R0 can be detected by an unphysical trend of Vg,\odot with the intrinsic rotation of different populations. From these two strategies I estimate R0 = (8.27 \pm 0.29) kpc and Vg,\odot = (250 \pm 9) kms-1 for a stellar sample from SEGUE, or respectively Vc = (238 \pm 9) kms-1. The result is consistent with the third estimator, where I use the angle of the mean motion of stars, which should follow the geometry of the Galactic disc. This method also gives the Solar radial motion with high accuracy. The rotation effect on U velocities must not be neglected when measuring the Solar radial velocity U\odot. It biases U\odot in any extended sample that is lop-sided in position angle α by of order 10 kms-1. Combining different methods I find U\odot \sim 14 kms-1, moderately higher than previous results from the Geneva-Copenhagen Survey.

preprint2012arXiv

Radial Flows and Angular Momentum Conservation in Galactic Chemical Evolution

We study the effects of radial flows on Galactic chemical evolution. A simple analytic scheme is developed prescribing the coupling of infall from the intergalactic medium and radial flows within the disc based on angular momentum conservation. We show that model parameters are tightly constrained by the observed [Fe/H]-abundance gradient in the Galactic disc. By this comparison the average rotational velocity of the onfalling material can be constrained to 0.7 < v/V_c < 0.75, or respectively ~ 160 km/s when assuming a constant disc circular velocity of V_c = 220 km/s. We test the robustness of this value against the influence of other processes. For a very simple model of inside-out formation this value changes only by Δv/V_c ~ 0.1, i.e. ~ 20 km/s, and significantly less on more realistic scenarios, showing that inside-out formation does not alone explain the abundance gradient. Effects of other uncertain parameters, e.g. star formation history and star formation efficiency have very small impact. Other drivers of inflow beyond our explicit modelling are assessed by adding a fixed inflow across the whole disc. The churning amplitude only mildly affects the results mostly by slightly flattening the metallicity gradient in the inner disc. A new process causing radial gas flows due to the ejection of material by stars moving on non-circular orbits is studied and seems to contribute negligibly to the total flows. We further show that gaseous outer discs cannot be the main source feeding the persistent star formation in the inner regions by a direct inflow.