-
Pattern speed evolution of barred galaxies in TNG50
Authors:
Marcin Semczuk,
Walter Dehnen,
Ralph Schönrich,
E. Athanassoula
Abstract:
Galactic bars are found in the majority of disc galaxies. They rotate nearly rigidly with an angular frequency called pattern speed. Previous idealised simulations have shown that bar pattern speed decreases with time due to dynamical friction exerted by the dark matter halo, while cold gas can reduce or even reverse this trend. We want to understand how different galaxy properties affect the evol…
▽ More
Galactic bars are found in the majority of disc galaxies. They rotate nearly rigidly with an angular frequency called pattern speed. Previous idealised simulations have shown that bar pattern speed decreases with time due to dynamical friction exerted by the dark matter halo, while cold gas can reduce or even reverse this trend. We want to understand how different galaxy properties affect the evolution of the bar pattern speed in more realistic situations, including ongoing star formation, mass infall, AGN feedback and galaxy interactions. We used the high-resolution run TNG50-1 of the magnetohydrodynamical cosmological simulations suite IllustrisTNG to trace the pattern speed of simulated bars and see how it depends on various galaxy properties. Simulated bars with initially high pattern speed and a subsequent rapid slowdown are more likely found in more massive galaxies. Lower mass galaxies, on the other hand, preferentially host bars that start at relatively low pattern speeds and retain the same value until the end of the simulation. More massive barred galaxies are also more affected by the AGN feedback, which removes (or heats up) the cold gas that could have prevented the slowdown. We find that bars grow and strengthen with slowdown, in agreement with higher resolution simulations. We find that strong correlations between the bar slowdown rate and galaxy mass weaken considerably when we use dimensionless measures to quantify the slowdown. In TNG50, the AGN feedback prescription amplifies the mass dependence. Turned around, this provides an interesting statistic to constrain subgrid physics by bar growth/slowing.
△ Less
Submitted 15 July, 2024;
originally announced July 2024.
-
Comment on 'a slightly oblate dark matter halo revealed by a retrograde precessing Galactic disk warp' by Huang et al
Authors:
Walter Dehnen,
Ralph Schönrich,
Ronald Drimmel,
Žofia Chrobáková,
Eloisa Poggio,
Marcin Semczuk
Abstract:
Huang et al. (2024) measured the derivative of the phase $φ_{\mathrm{w}}$ of the Galactic warp traced by classical Cepheids with respect to their age $τ$ and interpreted it as the warp precession rate $ω\equiv\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}t=-\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}τ$. This interpretation is unfounded: young stars follow trajectories close to those of their parental gas and trace t…
▽ More
Huang et al. (2024) measured the derivative of the phase $φ_{\mathrm{w}}$ of the Galactic warp traced by classical Cepheids with respect to their age $τ$ and interpreted it as the warp precession rate $ω\equiv\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}t=-\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}τ$. This interpretation is unfounded: young stars follow trajectories close to those of their parental gas and trace the instantaneous gas warp, not its shape at their time of birth: $φ_{\mathrm{w}}$ should hardly depend on Cepheid age. We show that the measured $\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}τ>0$ is consistent with an omitted-variable bias from neglecting the natural twist $\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}R$ of the warp and the $R$-$τ$ correlation for Cepheids (originating from the Galactic metallicity gradient and the Cepheid metallicity-age correlation).
△ Less
Submitted 8 July, 2024;
originally announced July 2024.
-
Radial Wave in the Galactic Disk: New Clues to Discriminate Different Perturbations
Authors:
Chengye Cao,
Zhao-Yu Li,
Ralph Schönrich,
Teresa Antoja
Abstract:
Decoding the key dynamical processes that shape the Galactic disk structure is crucial for reconstructing the Milky Way's evolution history. The second Gaia data release unveils a novel wave pattern in the $L_Z-\langle V_R\rangle$ space, but its formation mechanism remains elusive due to the intricate nature of involved perturbations and the challenges in disentangling their effects. Utilizing the…
▽ More
Decoding the key dynamical processes that shape the Galactic disk structure is crucial for reconstructing the Milky Way's evolution history. The second Gaia data release unveils a novel wave pattern in the $L_Z-\langle V_R\rangle$ space, but its formation mechanism remains elusive due to the intricate nature of involved perturbations and the challenges in disentangling their effects. Utilizing the latest Gaia DR3 data, we find that the $L_Z-\langle V_R\rangle$ wave systematically shifts towards lower $L_Z$ for dynamically hotter stars. The amplitude of this phase shift between stars of different dynamical hotness ($ΔL_Z$) peaks around $\mathrm{2300\,km\,s^{-1}\,kpc}$. To differentiate the role of different perturbations, we perform three sets of test particle simulations, wherein a satellite galaxy, corotating transient spiral arms, and a bar plus the corotating transient spiral arms act as the sole perturber, respectively. Under the satellite impact, the phase shift amplitude decreases towards higher $L_Z$, which we interpret through a toy model of radial phase mixing. While the corotating transient spiral arms do not generate an azimuthally universal phase shift variation pattern, combining the bar and spirals generates a characteristic $ΔL_Z$ peak at 2:1 Outer Lindblad Resonance, qualitatively resembling the observed feature. Therefore, the $L_Z-\langle V_R\rangle$ is more likely of internal origin. Furthermore, linking the $ΔL_Z$ peak to the 2:1 Lindblad resonance offers a novel approach to estimating the pattern speed of the Galactic Bar, supporting the long/slow bar model.
△ Less
Submitted 22 March, 2024;
originally announced March 2024.
-
Chemical Evolution in Galactic Nuclear Discs
Authors:
Jennifer K S Friske,
Ralph Schönrich
Abstract:
Nuclear Stellar Discs have been observed in the vast majority of barred disc galaxies including the Milky Way. Their intense star formation is sustained by the intense gas inflows driven by their surrounding bars and frequently supports a large-scale galactic fountain. Despite their central role in galaxy evolution, their chemical evolution remains largely unexplored. Here we present the first sys…
▽ More
Nuclear Stellar Discs have been observed in the vast majority of barred disc galaxies including the Milky Way. Their intense star formation is sustained by the intense gas inflows driven by their surrounding bars and frequently supports a large-scale galactic fountain. Despite their central role in galaxy evolution, their chemical evolution remains largely unexplored. Here we present the first systematic, multizonal modelling of the chemical evolution of these nuclear stellar discs. We argue that the chemical composition of nuclear stellar discs is best understood relative to the bar tips from which their gas is drawn. We make predictions of the detailed abundance profiles within the nuclear stellar disc under different accretion scenarios from the galactic bar and show by which observable differences they can be decided. We show that with their difference to normal disc star formation, nuclear discs offer a unique laboratory to break parameter degeneracies in chemical evolution models. This allows us to identify the effects of the main gas parameters and disentangle them from the global enrichment history. We show this for the example of the ejection fractions from the outer and nuclear disk.
△ Less
Submitted 3 October, 2023;
originally announced October 2023.
-
Roman Early-Definition Astrophysics Survey Opportunity: Galactic Roman Infrared Plane Survey (GRIPS)
Authors:
Roberta Paladini,
Catherine Zucker,
Robert Benjamin,
David Nataf,
Dante Minniti,
Gail Zasowski,
Joshua Peek,
Sean Carey,
Lori Allen,
Javier Alonso-Garcia,
Joao Alves,
Friederich Anders,
Evangelie Athanassoula,
Timothy C. Beers,
Jonathan Bird,
Joss Bland-Hwathorn,
Anthony Brown,
Sven Buder,
Luca Casagrande,
Andrew Casey,
Santi Cassisi,
Marcio Catelan,
Ranga-Ram Chary,
Andre-Nicolas Chene,
David Ciardi
, et al. (45 additional authors not shown)
Abstract:
A wide-field near-infrared survey of the Galactic disk and bulge/bar(s) is supported by a large representation of the community of Galactic astronomers. The combination of sensitivity, angular resolution and large field of view make Roman uniquely able to study the crowded and highly extincted lines of sight in the Galactic plane. A ~1000 deg2 survey of the bulge and inner Galactic disk would yiel…
▽ More
A wide-field near-infrared survey of the Galactic disk and bulge/bar(s) is supported by a large representation of the community of Galactic astronomers. The combination of sensitivity, angular resolution and large field of view make Roman uniquely able to study the crowded and highly extincted lines of sight in the Galactic plane. A ~1000 deg2 survey of the bulge and inner Galactic disk would yield an impressive dataset of ~120 billion sources and map the structure of our Galaxy. The effort would foster subsequent expansions in numerous dimensions (spatial, depth, wavelengths, epochs). Importantly, the survey would benefit from early defintion by the community, namely because the Galactic disk is a complex environment, and different science goals will require trade offs.
△ Less
Submitted 14 July, 2023;
originally announced July 2023.
-
JASMINE: Near-Infrared Astrometry and Time Series Photometry Science
Authors:
Daisuke Kawata,
Hajime Kawahara,
Naoteru Gouda,
Nathan J. Secrest,
Ryouhei Kano,
Hirokazu Kataza,
Naoki Isobe,
Ryou Ohsawa,
Fumihiko Usui,
Yoshiyuki Yamada,
Alister W. Graham,
Alex R. Pettitt,
Hideki Asada,
Junichi Baba,
Kenji Bekki,
Bryan N. Dorland,
Michiko Fujii,
Akihiko Fukui,
Kohei Hattori,
Teruyuki Hirano,
Takafumi Kamizuka,
Shingo Kashima,
Norita Kawanaka,
Yui Kawashima,
Sergei A. Klioner
, et al. (64 additional authors not shown)
Abstract:
Japan Astrometry Satellite Mission for INfrared Exploration (JASMINE) is a planned M-class science space mission by the Institute of Space and Astronautical Science, the Japan Aerospace Exploration Agency. JASMINE has two main science goals. One is the Galactic archaeology with Galactic Center Survey, which aims to reveal the Milky Way's central core structure and formation history from Gaia-level…
▽ More
Japan Astrometry Satellite Mission for INfrared Exploration (JASMINE) is a planned M-class science space mission by the Institute of Space and Astronautical Science, the Japan Aerospace Exploration Agency. JASMINE has two main science goals. One is the Galactic archaeology with Galactic Center Survey, which aims to reveal the Milky Way's central core structure and formation history from Gaia-level (~25 $μ$as) astrometry in the Near-Infrared (NIR) Hw-band (1.0-1.6 $μ$m). The other is the Exoplanet Survey, which aims to discover transiting Earth-like exoplanets in the habitable zone from NIR time-series photometry of M dwarfs when the Galactic center is not accessible. We introduce the mission, review many science objectives, and present the instrument concept. JASMINE will be the first dedicated NIR astrometry space mission and provide precise astrometric information of the stars in the Galactic center, taking advantage of the significantly lower extinction in the NIR. The precise astrometry is obtained by taking many short-exposure images. Hence, the JASMINE Galactic center survey data will be valuable for studies of exoplanet transits, asteroseismology, variable stars and microlensing studies, including discovery of (intermediate mass) black holes. We highlight a swath of such potential science, and also describe synergies with other missions.
△ Less
Submitted 4 March, 2024; v1 submitted 11 July, 2023;
originally announced July 2023.
-
A twisted and precessing Cepheid warp in the outer Milky Way disc
Authors:
Walter Dehnen,
Marcin Semczuk,
Ralph Schönrich
Abstract:
We examine the Galactic warp in a sample of all classical Cepheids with Gaia DR3 radial velocity. In each radial bin, we determine (1) the inclined plane normal to the mean orbital angular momentum of the stars and (2) that best fitting their positions. We find no warping inside $R\approx 11$ kpc; for larger $R$ the disc is increasingly inclined, reaching $i\sim 3^{\circ}$ at $R \ge 14$ kpc. With…
▽ More
We examine the Galactic warp in a sample of all classical Cepheids with Gaia DR3 radial velocity. In each radial bin, we determine (1) the inclined plane normal to the mean orbital angular momentum of the stars and (2) that best fitting their positions. We find no warping inside $R\approx 11$ kpc; for larger $R$ the disc is increasingly inclined, reaching $i\sim 3^{\circ}$ at $R \ge 14$ kpc. With larger $R$ the azimuth of the warp's ascending node shifts from $\varphi_{\mathrm{lon}}\approx-15^\circ$ at 11 kpc by about $14^{\circ}$/kpc in the direction of Galactic rotation, implying a leading spiral of nodes, the general behaviour of warped galaxies. From the method of fitting planes to the positions we also obtain $\dot{\varphi}_{\mathrm{lon}}$ and find prograde precession of $\dot{\varphi}_{\mathrm{lon}} \sim 12$ km/s/kpc at 12 kpc decreasing to $\sim 6$ km/s/kpc at 14 kpc and beyond. This would unwind the leading spiral of nodes in $\sim 100$ Myr, suggesting that our instantaneous measurements of $\dot{\varphi}_{\mathrm{lon}}$ reflect transient behaviour. This is consistent with existing simulations, which show oscillations in $\dot{\varphi}_{\mathrm{lon}}$ overlaying a long-term retrograde differential precession which generates the leading spiral of nodes.
△ Less
Submitted 16 May, 2023;
originally announced May 2023.
-
Measuring bar pattern speeds from single simulation snapshots
Authors:
Walter Dehnen,
Marcin Semczuk,
Ralph Schoenrich
Abstract:
We describe methods to measure simultaneously the orientation angle $ψ$ and pattern speed $Ω$ from single snapshots of simulated barred galaxies. Unlike previous attempts, our approach is unbiased, precise, and consistent in the sense that $ψ=\intΩ\mathrm{d}t$. It can be extended to obtain the rate and axis of rotation, i.e. the vector $\boldsymbolΩ$. We provide computer code implementing our meth…
▽ More
We describe methods to measure simultaneously the orientation angle $ψ$ and pattern speed $Ω$ from single snapshots of simulated barred galaxies. Unlike previous attempts, our approach is unbiased, precise, and consistent in the sense that $ψ=\intΩ\mathrm{d}t$. It can be extended to obtain the rate and axis of rotation, i.e. the vector $\boldsymbolΩ$. We provide computer code implementing our method.
△ Less
Submitted 7 November, 2022; v1 submitted 1 November, 2022;
originally announced November 2022.
-
Spiral arms and the angular momentum gap in Milky Way Cepheids
Authors:
Marcin Semczuk,
Walter Dehnen,
Ralph Schoenrich,
E. Athanassoula
Abstract:
The angular-momentum distribution of classical Cepheids in the outer Milky Way disc is bi-modal with a gap at $L_\mathrm{gap}=2950\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}$, corresponding to $R=13\,\mathrm{kpc}$, while no similar feature has been found in the general population of disc stars. We show that star formation in multiple spiral arm segments at the same azimuth leads to such multi-mod…
▽ More
The angular-momentum distribution of classical Cepheids in the outer Milky Way disc is bi-modal with a gap at $L_\mathrm{gap}=2950\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}$, corresponding to $R=13\,\mathrm{kpc}$, while no similar feature has been found in the general population of disc stars. We show that star formation in multiple spiral arm segments at the same azimuth leads to such multi-modality which quickly dissolves and only shows in young stars. Unlike other explanations, such as a 1:1 orbital resonance with the Galactic bar, this also accounts for the observed steepening of the stellar warp at $L_\mathrm{gap}$, since the adjacent spiral arms represent different parts of the warped gas disc. In this scenario the gap is clearly present only in young stars, as observed, while most purely stellar dynamical origins would affect all disc populations, including older disc stars.
△ Less
Submitted 21 October, 2022;
originally announced October 2022.
-
Age-morphology dependence of the Milky Way boxy/peanut bulge seen in Mira variables
Authors:
Marcin Semczuk,
Walter Dehnen,
Ralph Schoenrich,
E. Athanassoula
Abstract:
We analyse the distribution of Mira variable stars in the central region of the Milky Way. We find that with increasing period, i.e. decreasing age, the Miras shift towards negative Galactic longitudes $\ell$. Comparing to a cosmological zoom simulation of a barred galaxy, we find that this shift with age can be explained by an age-morphology dependence of the boxy peanut/bulge. Owing to a combina…
▽ More
We analyse the distribution of Mira variable stars in the central region of the Milky Way. We find that with increasing period, i.e. decreasing age, the Miras shift towards negative Galactic longitudes $\ell$. Comparing to a cosmological zoom simulation of a barred galaxy, we find that this shift with age can be explained by an age-morphology dependence of the boxy peanut/bulge. Owing to a combination of projection effects and the limitation of the range of Galactic longitudes, the near hump at $\ell>0$ is more truncated for younger populations, and the far hump at $\ell<0$ dominates the observed distributions.
△ Less
Submitted 9 June, 2022;
originally announced June 2022.
-
Capture of interstellar objects II: by the Solar system
Authors:
Walter Dehnen,
Thomas O. Hands,
Ralph Schönrich
Abstract:
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 cr…
▽ More
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.
△ Less
Submitted 14 December, 2021;
originally announced December 2021.
-
Metallicity-Suppressed Collapsars Cannot be the Dominant r-Process Source in the Milky Way
Authors:
Jack Fraser,
Ralph Schönrich
Abstract:
We develop a high-performance analytical model of Galactic Chemical Evolution, which accounts for delay time distributions and lock-up of stellar yields in a thermal-phased ISM. The model is capable of searching, for the first time, through the high-dimensional parameter space associated with the r-process enrichment of the Milky Way by its possible sources: Neutron Star Mergers and Collapsar even…
▽ More
We develop a high-performance analytical model of Galactic Chemical Evolution, which accounts for delay time distributions and lock-up of stellar yields in a thermal-phased ISM. The model is capable of searching, for the first time, through the high-dimensional parameter space associated with the r-process enrichment of the Milky Way by its possible sources: Neutron Star Mergers and Collapsar events. Their differing formation mechanisms give these two processes different time dependencies, a property which has frequently been used to argue in favour of collapsars as the dominant r-process source. However, we show that even with large degrees of freedom in the allowed thermal, structural, and chemical properties of the galaxy, large regions of parameter space are in strong tension with the data. In particular, whilst we are able to find models in which neutron star mergers produce the majority of r-process material, the data rule out all models with dominant collapsar yields. With no other identified source, we conclude that Neutron Star Mergers must be the dominant contributors to the modern Milky Way r-process budget.
△ Less
Submitted 22 November, 2021;
originally announced November 2021.
-
Oscillating dynamical friction on galactic bars by trapped dark matter
Authors:
Rimpei Chiba,
Ralph Schönrich
Abstract:
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 with…
▽ More
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.
△ Less
Submitted 28 April, 2022; v1 submitted 22 September, 2021;
originally announced September 2021.
-
The K2 Galactic Archaeology Program Data Release 3: Age-abundance patterns in C1-C8 and C10-C18
Authors:
Joel C. Zinn,
Dennis Stello,
Yvonne Elsworth,
Rafael A. García,
Thomas Kallinger,
Savita Mathur,
Benoît Mosser,
Marc Hon,
Lisa Bugnet,
Caitlin Jones,
Claudia Reyes,
Sanjib Sharma,
Ralph Schönrich,
Jack T. Warfield,
Rodrigo Luger,
Andrew Vanderburg,
Chiaki Kobayashi,
Marc H. Pinsonneault,
Jennifer A. Johnson,
Daniel Huber,
Sven Buder,
Meridith Joyce,
Joss Bland-Hawthorn,
Luca Casagrande,
Geraint F. Lewis
, et al. (6 additional authors not shown)
Abstract:
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 s…
▽ More
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.
△ Less
Submitted 10 March, 2022; v1 submitted 11 August, 2021;
originally announced August 2021.
-
The small boxy/peanut structure of the Milky Way traced by old stars
Authors:
Marcin Semczuk,
Walter Dehnen,
Ralph Schoenrich,
E. Athanassoula
Abstract:
We analyse the positions of RR Lyrae stars in the central region of the Milky Way. In addition to the overall bar shape detected previously, we find evidence for a peanut shaped structure, in form of overdensities near $\ell=-2$ deg and $1$ deg at $b\sim-3.5$ deg. The corresponding physical distance between the two peaks of the peanut is $\sim0.7$ kpc, significantly shorter than that found from ne…
▽ More
We analyse the positions of RR Lyrae stars in the central region of the Milky Way. In addition to the overall bar shape detected previously, we find evidence for a peanut shaped structure, in form of overdensities near $\ell=-2$ deg and $1$ deg at $b\sim-3.5$ deg. The corresponding physical distance between the two peaks of the peanut is $\sim0.7$ kpc, significantly shorter than that found from near-IR images (3.3 kpc) and red-clump stars. Qualitatively this is expected from `fractionation' of bar orbits, which we demonstrate to be matched in a simulation of an inside-out growing disc subsequently forming a bar.
△ Less
Submitted 11 November, 2021; v1 submitted 22 July, 2021;
originally announced July 2021.
-
Age distribution of stars in boxy/peanut/X-shaped bulges formed without bar buckling
Authors:
Junichi Baba,
Daisuke Kawata,
Ralph Schönrich
Abstract:
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 witho…
▽ More
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.
△ Less
Submitted 28 February, 2022; v1 submitted 19 April, 2021;
originally announced April 2021.
-
Tree-ring structure of Galactic bar resonance
Authors:
Rimpei Chiba,
Ralph Schönrich
Abstract:
Galaxy models have long predicted that galactic bars slow down by losing angular momentum to their postulated dark haloes. When the bar slows down, resonance sweeps radially outwards through the galactic disc while growing in volume, thereby sequentially capturing new stars at its surface/separatrix. Since trapped stars conserve their action of libration, which measures the relative distance to th…
▽ More
Galaxy models have long predicted that galactic bars slow down by losing angular momentum to their postulated dark haloes. When the bar slows down, resonance sweeps radially outwards through the galactic disc while growing in volume, thereby sequentially capturing new stars at its surface/separatrix. Since trapped stars conserve their action of libration, which measures the relative distance to the resonance centre, the order of capturing is preserved: the surface of the resonance is dominated by stars captured recently at large radius, while the core of the resonance is occupied by stars trapped early at small radius. The slow-down of the bar thus results in a rising mean metallicity of trapped stars from the surface towards the centre of the resonance as the Galaxy's metallicity declines towards large radii. This argument, when applied to Solar neighbourhood stars, allows a novel precision measurement of the bar's current pattern speed $Ω_p = 35.5 \pm 0.8$ km/s/kpc, placing the corotation radius at $R_{CR} = 6.6 \pm 0.2$ kpc. With this pattern speed, the corotation resonance precisely fits the Hercules stream in agreement with kinematics. Beyond corroborating the slow bar theory, this measurement manifests the deceleration of the bar of more than 24% since its formation and thus the angular momentum transfer to the dark halo by dynamical friction. The measurement therefore supports the existence of a standard dark-matter halo rather than alternative models of gravity.
△ Less
Submitted 15 April, 2021; v1 submitted 16 February, 2021;
originally announced February 2021.
-
Galactic Bar Resonances Inferred from Kinematically Hot Stars in Gaia EDR3
Authors:
Daisuke Kawata,
Junichi Baba,
Jason A. S. Hunt,
Ralph Schönrich,
Ioana Ciucǎ,
Jennifer Friske,
George Seabroke,
Mark Cropper
Abstract:
Using a numerical simulation of an isolated barred disc galaxy, we first demonstrate that the resonances of the inner bar structure induce more prominent features in the action space distribution for the kinematically hotter stars, which are less sensitive to the local perturbation, such as the transient spiral arms. Then, we analyse the action distribution for the kinematically hotter stars selec…
▽ More
Using a numerical simulation of an isolated barred disc galaxy, we first demonstrate that the resonances of the inner bar structure induce more prominent features in the action space distribution for the kinematically hotter stars, which are less sensitive to the local perturbation, such as the transient spiral arms. Then, we analyse the action distribution for the kinematically hotter stars selected from the Gaia EDR3 data as the stars with higher values of radial and vertical actions. We find several resonance features, including two new features, in the angular momentum distribution similar to what are seen in our numerical simulations. We show that the bar pattern speeds of about $Ω_{\rm bar}\sim34$~km~s$^{-1}$~kpc$^{-1}$ and 42~km~s$^{-1}$~kpc$^{-1}$ explain all these features equally well. The resonance features we find correspond to the inner 4:1, co-rotation, outer 4:1, outer Lindblad and outer 4:3 (co-rotation, outer 4:1, outer Lindblad, outer 4:3 and outer 1:1) resonances, when $Ω_{\rm bar}\sim34$ (42) km~s$^{-1}$~kpc$^{-1}$ is assumed.
△ Less
Submitted 8 September, 2021; v1 submitted 10 December, 2020;
originally announced December 2020.
-
The K2 Galactic Archaeology Program Data Release 2: Asteroseismic results from campaigns 4, 6, & 7
Authors:
Joel C. Zinn,
Dennis Stello,
Yvonne Elsworth,
Rafael A. García,
Thomas Kallinger,
Savita Mathur,
Benoît Mosser,
Lisa Bugnet,
Caitlin Jones,
Marc Hon,
Sanjib Sharma,
Ralph Schönrich,
Jack T. Warfield,
Rodrigo Luger,
Marc H. Pinsonneault,
Jennifer A. Johnson,
Daniel Huber,
Victor Silva Aguirre,
William J. Chaplin,
Guy R. Davies,
Andrea Miglio
Abstract:
Studies of Galactic structure and evolution have benefitted enormously from Gaia kinematic information, though additional, intrinsic stellar parameters like age are required to best constrain Galactic models. Asteroseismology is the most precise method of providing such information for field star populations $\textit{en masse}$, but existing samples for the most part have been limited to a few nar…
▽ More
Studies of Galactic structure and evolution have benefitted enormously from Gaia kinematic information, though additional, intrinsic stellar parameters like age are required to best constrain Galactic models. Asteroseismology is the most precise method of providing such information for field star populations $\textit{en masse}$, but existing samples for the most part have been limited to a few narrow fields of view by the CoRoT and Kepler missions. In an effort to provide well-characterized stellar parameters across a wide range in Galactic position, we present the second data release of red giant asteroseismic parameters for the K2 Galactic Archaeology Program (GAP). We provide $ν_{\mathrm{max}}$ and $Δν$ based on six independent pipeline analyses; first-ascent red giant branch (RGB) and red clump (RC) evolutionary state classifications from machine learning; and ready-to-use radius & mass coefficients, $κ_R$ & $κ_M$, which, when appropriately multiplied by a solar-scaled effective temperature factor, yield physical stellar radii and masses. In total, we report 4395 radius and mass coefficients, with typical uncertainties of $3.3\% \mathrm{\ (stat.)} \pm 1\% \mathrm{\ (syst.)}$ for $κ_R$ and $7.7\% \mathrm{\ (stat.)} \pm 2\% \mathrm{\ (syst.)}$ for $κ_M$ among RGB stars, and $5.0\% \mathrm{\ (stat.)} \pm 1\% \mathrm{\ (syst.)}$ for $κ_R$ and $10.5\% \mathrm{\ (stat.)} \pm 2\% \mathrm{\ (syst.)}$ for $κ_M$ among RC stars. We verify that the sample is nearly complete -- except for a dearth of stars with $ν_{\mathrm{max}} \lesssim 10-20μ$Hz -- by comparing to Galactic models and visual inspection. Our asteroseismic radii agree with radii derived from Gaia Data Release 2 parallaxes to within $2.2 \pm 0.3\%$ for RGB stars and $2.0 \pm 0.6\%$ for RC stars.
△ Less
Submitted 7 December, 2020;
originally announced December 2020.
-
Jeans modelling of the Milky Way's nuclear stellar disc
Authors:
Mattia C. Sormani,
John Magorrian,
Francisco Nogueras-Lara,
Nadine Neumayer,
Ralph Schoenrich,
Ralf S. Klessen,
Alessandra Mastrobuono-Battisti
Abstract:
The nuclear stellar disc (NSD) is a flattened stellar structure that dominates the gravitational potential of the Milky Way at Galactocentric radii $30 \lesssim R \lesssim 300{\, \rm pc}$. In this paper, we construct axisymmetric Jeans dynamical models of the NSD based on previous photometric studies and we fit them to line-of-sight kinematic data of APOGEE and SiO maser stars. We find that (i) th…
▽ More
The nuclear stellar disc (NSD) is a flattened stellar structure that dominates the gravitational potential of the Milky Way at Galactocentric radii $30 \lesssim R \lesssim 300{\, \rm pc}$. In this paper, we construct axisymmetric Jeans dynamical models of the NSD based on previous photometric studies and we fit them to line-of-sight kinematic data of APOGEE and SiO maser stars. We find that (i) the NSD mass is lower but consistent with the mass independently determined from photometry by Launhardt et al. (2002). Our fiducial model has a mass contained within spherical radius $r=100{\, \rm pc}$ of $M(r<100{\, \rm pc}) = 3.9 \pm 1 \times 10^8 {\, \rm M_\odot}$ and a total mass of $M_{\rm NSD} = 6.9 \pm 2 \times 10^8 {\, \rm M_\odot}$. (ii) The NSD might be the first example of a vertically biased disc, i.e. with ratio between the vertical and radial velocity dispersion $σ_z/σ_R>1$. Observations and theoretical models of the star-forming molecular gas in the central molecular zone suggest that large vertical oscillations may be already imprinted at stellar birth. However, the finding $σ_z/σ_R > 1$ depends on a drop in the velocity dispersion in the innermost few tens of parsecs, on our assumption that the NSD is axisymmetric, and that the available (extinction corrected) stellar samples broadly trace the underlying light and mass distributions, all of which need to be established by future observations and/or modelling. (iii) We provide the most accurate rotation curve to date for the innermost $500 {\, \rm pc}$ of our Galaxy.
△ Less
Submitted 9 September, 2020; v1 submitted 13 July, 2020;
originally announced July 2020.
-
Mapping the Galactic disk with the LAMOST and Gaia Red clump sample: I: precise distances, masses, ages and 3D velocities of $\sim$ 140000 red clump stars
Authors:
Yang Huang,
Ralph Schonrich,
Huawei Zhang,
Yaqian Wu,
Bingqiu Chen,
Haifeng Wang,
Maosheng Xiang,
Chun Wang,
Haibo Yuan,
Xinyi Li,
Weixiang Sun,
Ji Li,
Xiaowei Liu
Abstract:
We present a sample of $\sim$ 140,000 primary red clump (RC) stars of spectral signal-to-noise ratios higher than 20 from the LAMOST Galactic spectroscopic surveys, selected based on their positions in the metallicity-dependent effective temperature--surface gravity and color--metallicity diagrams, supervised by high-quality $Kepler$ asteroseismology data. The stellar masses and ages of those star…
▽ More
We present a sample of $\sim$ 140,000 primary red clump (RC) stars of spectral signal-to-noise ratios higher than 20 from the LAMOST Galactic spectroscopic surveys, selected based on their positions in the metallicity-dependent effective temperature--surface gravity and color--metallicity diagrams, supervised by high-quality $Kepler$ asteroseismology data. The stellar masses and ages of those stars are further determined from the LAMOST spectra, using the Kernel Principal Component Analysis method, trained with thousands of RCs in the LAMOST-$Kepler$ fields with accurate asteroseismic mass measurements. The purity and completeness of our primary RC sample are generally higher than 80 per cent. For the mass and age, a variety of tests show typical uncertainties of 15 and 30 per cent, respectively. Using over ten thousand primary RCs with accurate distance measurements from the parallaxes of Gaia DR2, we re-calibrate the $K_{\rm s}$ absolute magnitudes of primary RCs by, for the first time, considering both the metallicity and age dependencies. With the the new calibration, distances are derived for all the primary RCs, with a typical uncertainty of 5--10 per cent, even better than the values yielded by the Gaia parallax measurements for stars beyond 3--4 kpc. The sample covers a significant volume of the Galactic disk of $4 \leq R \leq 16$ kpc, $|Z| \leq 5$ kpc, and $-20 \leq φ\leq 50^{\circ}$. Stellar atmospheric parameters, line-of-sight velocities and elemental abundances derived from the LAMOST spectra and proper motions of Gaia DR2 are also provided for the sample stars. Finally, the selection function of the sample is carefully evaluated in the color-magnitude plane for different sky areas. The sample is publicly available.
△ Less
Submitted 5 June, 2020; v1 submitted 4 June, 2020;
originally announced June 2020.
-
Resonance sweeping by a decelerating Galactic bar
Authors:
Rimpei Chiba,
Jennifer K. S. Friske,
Ralph Schönrich
Abstract:
We provide the first quantitative evidence for the deceleration of the Galactic bar from local stellar kinematics in agreement with dynamical friction by a typical dark matter halo. The kinematic response of the stellar disk to a decelerating bar is studied using secular perturbation theory and test particle simulations. We show that the velocity distribution at any point in the disk affected by a…
▽ More
We provide the first quantitative evidence for the deceleration of the Galactic bar from local stellar kinematics in agreement with dynamical friction by a typical dark matter halo. The kinematic response of the stellar disk to a decelerating bar is studied using secular perturbation theory and test particle simulations. We show that the velocity distribution at any point in the disk affected by a naturally slowing bar is qualitatively different from that perturbed by a steadily rotating bar with the same current pattern speed $Ω_{\rm p}$ and amplitude. When the bar slows down, its resonances sweep through phase space, trapping and dragging along a portion of previously free orbits. This enhances occupation on resonances, but also changes the distribution of stars within the resonance. Due to the accumulation of orbits near the boundary of the resonance, the decelerating bar model reproduces with its corotation resonance the offset and strength of the Hercules stream in the local $v_R$-$v_\varphi$ plane and the double-peaked structure of mean $v_R$ in the $L_z$-$\varphi$ plane. At resonances other than the corotation, resonant dragging by a slowing bar is associated with a continuing increase in radial action, leading to multiple resonance ridges in the action plane as identified in the Gaia data. This work shows models using a constant bar pattern speed likely lead to qualitatively wrong conclusions. Most importantly we provide a quantitative estimate of the current slowing rate of the bar $dΩ_{\rm p}/dt = (-4.5 \pm 1.4)~{\rm km} {\rm s}^{-1} {\rm kpc}^{-1} {\rm Gyr}^{-1}$ with additional systematic uncertainty arising from unmodeled impacts of e.g. spiral arms.
△ Less
Submitted 14 November, 2020; v1 submitted 9 December, 2019;
originally announced December 2019.
-
The Tilt of the Local Velocity Ellipsoid as Seen by Gaia
Authors:
A. Everall,
N. W. Evans,
V. Belokurov,
R. Schönrich
Abstract:
The Gaia Radial Velocity Spectrometer (RVS) provides a sample of 7,224,631 stars with full six-dimensional phase space information. Bayesian distances of these stars are available from the catalogue of Schönrich et al. (2019). We exploit this to map out the behaviour of the velocity ellipsoid within 5 kpc of the Sun. We find that the tilt of the disc-dominated RVS sample is accurately described by…
▽ More
The Gaia Radial Velocity Spectrometer (RVS) provides a sample of 7,224,631 stars with full six-dimensional phase space information. Bayesian distances of these stars are available from the catalogue of Schönrich et al. (2019). We exploit this to map out the behaviour of the velocity ellipsoid within 5 kpc of the Sun. We find that the tilt of the disc-dominated RVS sample is accurately described by the relation $α= (0.952 \pm 0.007)\arctan (|z|/R)$, where ($R,z$) are cylindrical polar coordinates. This corresponds to velocity ellipsoids close to spherical alignment (for which the normalising constant would be unity) and pointing towards the Galactic centre. Flattening of the tilt of the velocity ellipsoids is enhanced close to the plane and Galactic centre, whilst at high elevations far from the Galactic center the population is consistent with exact spherical alignment. Using the LAMOST catalogue cross-matched with Gaia DR2, we construct thin disc and halo samples of reasonable purity based on metallicity. We find that the tilt of thin disc stars straddles $α= (0.909-1.038)\arctan (|z|/R)$, and of halo stars straddles $α= (0.927-1.063)\arctan (|z|/R)$. We caution against the use of reciprocal parallax for distances in studies of the tilt, as this can lead to serious artefacts.
△ Less
Submitted 7 August, 2019; v1 submitted 17 April, 2019;
originally announced April 2019.
-
More than just a wrinkle: A wave-like pattern in radial velocity vs. angular momentum from Gaia Data
Authors:
Jennifer Friske,
Ralph Schönrich
Abstract:
We present a newly found wave--like pattern in mean Galactocentric radial velocity U_g vs. guiding centre radius R_g or angular momentum L_z of stars in the RV subsample of Gaia DR2. The short-wave pattern has a wavelength of order 1.2 kpc in R_g or 285 kpc km/s in L_z. The pattern shows only weak changes with Galactocentric radius R and little change in strength in particular with the vertical en…
▽ More
We present a newly found wave--like pattern in mean Galactocentric radial velocity U_g vs. guiding centre radius R_g or angular momentum L_z of stars in the RV subsample of Gaia DR2. The short-wave pattern has a wavelength of order 1.2 kpc in R_g or 285 kpc km/s in L_z. The pattern shows only weak changes with Galactocentric radius R and little change in strength in particular with the vertical energy E_z of the stars or respectively the distance to the Galactic plane |z|. The pattern is to first order symmetric around the plane, i.e. has no significant odd terms in z. There is weak phase shift with the pattern moving towards slightly lower L_z (i.e. trailing) with |z| and E_z. However, we observe a highly significant phase shift in Galactic azimuth phi, which is different for different peaks. The peak around L_z ~ 2100 kpc km/s only shows a weak change with galactic azimuth, while the rest of the pattern shows a clearly detectable shift of d L_z /d phi = (200 +/- 22) kpc km/s/rad. If we consider all peaks to belong to the same pattern, this would suggest a wavenumber m = 4. We further find that the wave-like pattern in U_g appears to be related to the W vs. L_z pattern detected in Gaia DR1. A comparison of the U_g-L_z wave pattern with changes of U_g vs. R, which have been previously discussed, suggests that the latter can be understood as just the U_g-L_z pattern washed out by blurring (i.e. orbital excursions around their guiding-centre) of disc stars.
△ Less
Submitted 17 October, 2019; v1 submitted 25 February, 2019;
originally announced February 2019.
-
Distances and parallax bias in Gaia DR2
Authors:
Ralph Schönrich,
Paul McMillan,
Laurent Eyer
Abstract:
We derive Bayesian distances for all stars in the RV sample of Gaia DR2, and use the statistical method of Schoenrich, Binney & Asplund(2012) to validate the distances and test the Gaia parallaxes. In contrast to other methods, which rely on special sources, our method directly tests the distances to all stars in our sample. We find clear evidence for a near-linear trend of distance bias f with di…
▽ More
We derive Bayesian distances for all stars in the RV sample of Gaia DR2, and use the statistical method of Schoenrich, Binney & Asplund(2012) to validate the distances and test the Gaia parallaxes. In contrast to other methods, which rely on special sources, our method directly tests the distances to all stars in our sample. We find clear evidence for a near-linear trend of distance bias f with distance s, proving a parallax offset delta p. On average, we find delta p = -0.054 mas (parallaxes in Gaia DR2 need to be increased) when accounting for the parallax uncertainty under-estimate in the Gaia set (delta p = -0.048 mas on the raw parallax errors) with negligible formal error and a systematic uncertainty of about 0.006 mas. The value is in concordance with results from asteroseismic measurements, but differs from the much lower bias found on quasar samples. We further use our method to compile a comprehensive set of quality cuts in colour, apparent magnitude, and astrometric parameters. Last, we find that for this sample delta p appears to strongly depend on the parallax uncertainty sigma_p (when including the additional 0.043 mas) with a statistical confidence far in excess of 10σand a proportionality factor close to 1, though the dependence varies somewhat with sigma_p. Correcting for the sigma_p dependence also resolves otherwise unexplained correlations of the offset with the number of observation periods n_{vis} and ecliptic latitude. Every study using Gaia DR2 parallaxes/distances should investigate the sensitivity of their results on the parallax biases described here and - for fainter samples - in the DR2 astrometry paper.
△ Less
Submitted 29 May, 2019; v1 submitted 6 February, 2019;
originally announced February 2019.
-
The chemical evolution of r-process elements from neutron star mergers: the role of a 2-phase interstellar medium
Authors:
Ralph Schönrich,
David H. Weinberg
Abstract:
Neutron star mergers (NM) are a plausible source of heavy r-process elements such as Europium, but previous chemical evolution models have either failed to reproduce the observed Europium trends for Milky Way thick disc stars (with [Fe/H] ~ -1) or have done so only by adopting unrealistically short merger timescales. Using analytic arguments and numerical simulations, we demonstrate that models wi…
▽ More
Neutron star mergers (NM) are a plausible source of heavy r-process elements such as Europium, but previous chemical evolution models have either failed to reproduce the observed Europium trends for Milky Way thick disc stars (with [Fe/H] ~ -1) or have done so only by adopting unrealistically short merger timescales. Using analytic arguments and numerical simulations, we demonstrate that models with a single-phase interstellar medium (ISM) and metallicity-independent yields cannot reproduce observations showing [Eu/alpha] > 0 or [Eu/Fe] > [alpha/Fe] for alpha-elements such as Mg and Si. However, this problem is easily resolved if we allow for a 2-phase ISM, with hot-phase cooling times τ_{cool} of order 1 Gyr and a larger fraction of NM yields injected directly into the cold star-forming phase relative to alpha-element yields from core collapse supernovae (ccSNe). We find good agreement with observations in models with a cold phase injection ratio f_{c,NM}/f_{c,ccSN} of order 2, and a characteristic merger timescale τ_NM=150 Myr. We show that the observed super-solar [Eu/alpha] at intermediate metallicities implies that a significant fraction of Eu originates from NM or another source besides ccSNe, and that these non-ccSN yields are preferentially deposited in the star-forming phase of the ISM at early times.
△ Less
Submitted 28 May, 2019; v1 submitted 28 January, 2019;
originally announced January 2019.
-
The origin of the Gaia phase-plane spiral
Authors:
James Binney,
Ralph Schoenrich
Abstract:
A simple model is presented of the formation of the spiral the (z,v_z) phase plane of solar-neighbourhood stars that was recently discovered in Gaia data. The key is that the frequency Omega_z at which stars oscillate vertically depends on angular momentum about the z axis in addition to the amplitude of the star's vertical oscillations. Spirals should form in both <v_phi> and <v_R> whenever a mas…
▽ More
A simple model is presented of the formation of the spiral the (z,v_z) phase plane of solar-neighbourhood stars that was recently discovered in Gaia data. The key is that the frequency Omega_z at which stars oscillate vertically depends on angular momentum about the z axis in addition to the amplitude of the star's vertical oscillations. Spirals should form in both <v_phi> and <v_R> whenever a massive substructure, such as the Sgr dwarf galaxy, passes through the Galactic plane. The model yields similar spirals to those observed in both <v_phi> and <v_R>. The primary driver is the component of the tidal force that lies in the plane. We investigate the longevity of the spirals and the mass of the substructure, but the approximations inherent in the model make quantitative results unreliable. The work relies heavily on a self-consistent, multi-component model of our Galaxy produced by the AGAMA package for f(J) modelling.
△ Less
Submitted 25 July, 2018;
originally announced July 2018.
-
On the kinematic signature of the Galactic warp as revealed by the LAMOST-TGAS data
Authors:
Yang Huang,
R. Schönrich,
Xiaowei Liu,
Bingqiu Chen,
Huawei Zhang,
Haibo Yuan,
Maosheng Xiang,
Chun Wang,
Zhijia Tian
Abstract:
Using a sample of about 123,000 stars with accurate 3D velocity measurements from the LAMOST-TGAS data, we confirm the kinematic signature of the Galactic warp found by Schönrich & Dehnen recently. The data reveal a clear trend of increasing mean vertical velocity $\overline{V_{z}}$ as a function of absolute vertical angular momentum $L_{z}$ and azimuthal velocity $V_φ$ for guiding center radius…
▽ More
Using a sample of about 123,000 stars with accurate 3D velocity measurements from the LAMOST-TGAS data, we confirm the kinematic signature of the Galactic warp found by Schönrich & Dehnen recently. The data reveal a clear trend of increasing mean vertical velocity $\overline{V_{z}}$ as a function of absolute vertical angular momentum $L_{z}$ and azimuthal velocity $V_φ$ for guiding center radius $R_{g}$ between 6.0 and 10.5 kpc. The trend is consistent with a large-scale Galactic warp. Similar to Schönrich & Dehnen, we also find a wave-like pattern of $\overline{V_{z}}$ versus $L_{z}$ with an amplitude of $\sim 0.9$ km s$^{-1}$ on a scale of $\sim 2.0$ kpc, which could arise from bending waves or a winding warp. Finally, we confirm a prominent, localized peak in $\overline{V_z}$ near $L_z \sim 2150$ kpc km s$^{-1}$ (corresponding to $R_{g} \sim 9$ kpc and $V_φ \sim 255$ km s$^{-1}$). The additional line-of-sight velocity information from LAMOST reveals that stars in this feature have a large, inward radial velocity of $V_{R} \sim -13.33 \pm 0.59$ km s$^{-1}$ and a small radial velocity dispersion $σ_{R} \sim 25.27 \pm 0.89$ km s$^{-1}$, suggesting that a stellar stream gives rise to this feature.
△ Less
Submitted 8 July, 2018;
originally announced July 2018.
-
Two chemically similar stellar overdensities on opposite sides of the plane of the Galaxy
Authors:
Maria Bergemann,
Branimir Sesar,
Judith G. Cohen,
Aldo M. Serenelli,
Allyson Sheffield,
Ting S. Li,
Luca Casagrande,
Kathryn V. Johnston,
Chervin F. P. Laporte,
Adrian M. Price-Whelan,
Ralph Schoenrich,
Andrew Gould
Abstract:
Our Galaxy is thought to have undergone an active evolutionary history dominated by star formation, the accretion of cold gas, and, in particular, mergers up to 10 gigayear ago. The stellar halo reveals rich fossil evidence of these interactions in the form of stellar streams, substructures, and chemically distinct stellar components. The impact of dwarf galaxy mergers on the content and morpholog…
▽ More
Our Galaxy is thought to have undergone an active evolutionary history dominated by star formation, the accretion of cold gas, and, in particular, mergers up to 10 gigayear ago. The stellar halo reveals rich fossil evidence of these interactions in the form of stellar streams, substructures, and chemically distinct stellar components. The impact of dwarf galaxy mergers on the content and morphology of the Galactic disk is still being explored. Recent studies have identified kinematically distinct stellar substructures and moving groups, which may have extragalactic origin. However, there is mounting evidence that stellar overdensities at the outer disk/halo interface could have been caused by the interaction of a dwarf galaxy with the disk. Here we report detailed spectroscopic analysis of 14 stars drawn from two stellar overdensities, each lying about 5 kiloparsecs above and below the Galactic plane - locations suggestive of association with the stellar halo. However, we find that the chemical compositions of these stars are almost identical, both within and between these groups, and closely match the abundance patterns of the Milky Way disk stars. This study hence provides compelling evidence that these stars originate from the disk and the overdensities they are part of were created by tidal interactions of the disk with passing or merging dwarf galaxies.
△ Less
Submitted 1 March, 2018;
originally announced March 2018.
-
Warp, Waves, and Wrinkles in the Milky Way
Authors:
Ralph Schönrich,
Walter Dehnen
Abstract:
We derive unbiased distance estimates for the Gaia-TGAS dataset by correcting for the bias due to the distance dependence of the selection function, which we measure directly from the data. From these distances and proper motions, we estimate the vertical and azimuthal velocities, $W$ and $V_φ$, and angular momentum $L_z$ for stars in the Galactic centre and anti-centre directions. The resulting m…
▽ More
We derive unbiased distance estimates for the Gaia-TGAS dataset by correcting for the bias due to the distance dependence of the selection function, which we measure directly from the data. From these distances and proper motions, we estimate the vertical and azimuthal velocities, $W$ and $V_φ$, and angular momentum $L_z$ for stars in the Galactic centre and anti-centre directions. The resulting mean vertical motion $\overline{W}$ shows a linear increase with both $V_φ$ and $L_z$ at $10 σ$ significance. Such a trend is expected from and consistent with the known Galactic warp. This signal extends to stars with guiding centre radii $R_g<R_0$, placing the onset of the warp at $R\lesssim7{\rm kpc}$. At equally high significance, we detect a previously unknown wave-like pattern of $\overline{W}$ over guiding centre $R_g$ with amplitude $\sim1{\rm kms}^{-1}$ and wavelength $\sim2.5{\rm kpc}$. This pattern is present in both the centre and anti-centre directions, consistent with a winding (corrugated) warp or bending wave, likely related to known features in the outer disc (TriAnd and Monoceros over-densities), and may be caused by the interaction with the Sgr dwarf galaxy $\sim1{\rm Gyr}$ ago. The only significant deviation from this simple fit is a stream-like feature near $R_g\sim9{\rm kpc}$ ($|L_z|\sim2150{\rm kpckms}^{-1}$).
△ Less
Submitted 19 June, 2018; v1 submitted 18 December, 2017;
originally announced December 2017.
-
Confirming chemical clocks: asteroseismic age dissection of the Milky Way disk(s)
Authors:
V. Silva Aguirre,
M. Bojsen-Hansen,
D. Slumstrup,
L. Casagrande,
D. Kawata,
I. Ciuca,
R. Handberg,
M. N. Lund,
J. R. Mosumgaard,
D. Huber,
J. A. Johnson,
M. H. Pinsonneault,
A. M. Serenelli,
D. Stello,
J. Tayar,
J. C. Bird,
S. Cassisi,
M. Hon,
M. Martig,
P. E. Nissen,
H. W. Rix,
R. Schönrich,
C. Sahlholdt,
W. H. Trick,
J. Yu
Abstract:
Investigations of the origin and evolution of the Milky Way disk have long relied on chemical and kinematic identification of its components to reconstruct our Galactic past. Difficulties in determining precise stellar ages have restricted most studies to small samples, normally confined to the solar neighbourhood. Here we break this impasse with the help of asteroseismic inference and perform a c…
▽ More
Investigations of the origin and evolution of the Milky Way disk have long relied on chemical and kinematic identification of its components to reconstruct our Galactic past. Difficulties in determining precise stellar ages have restricted most studies to small samples, normally confined to the solar neighbourhood. Here we break this impasse with the help of asteroseismic inference and perform a chronology of the evolution of the disk throughout the age of the Galaxy. We chemically dissect the Milky Way disk population using a sample of red giant stars spanning out to 2~kpc in the solar annulus observed by the {\it Kepler} satellite, with the added dimension of asteroseismic ages. Our results reveal a clear difference in age between the low- and high-$α$ populations, which also show distinct velocity dispersions in the $V$ and $W$ components. We find no tight correlation between age and metallicity nor [$α$/Fe] for the high-$α$ disk stars. Our results indicate that this component formed over a period of more than 2~Gyr with a wide range of [M/H] and [$α$/Fe] independent of time. Our findings show that the kinematic properties of young $α$-rich stars are consistent with the rest of the high-$α$ population and different from the low-$α$ stars of similar age, rendering support to their origin being old stars that went through a mass transfer or stellar merger event, making them appear younger, instead of migration of truly young stars formed close to the Galactic bar.
△ Less
Submitted 15 January, 2018; v1 submitted 26 October, 2017;
originally announced October 2017.
-
Migration and kinematics in growing disc galaxies with thin and thick discs
Authors:
Michael Aumer,
James Binney,
Ralph Schönrich
Abstract:
We analyse disc heating and radial migration in N-body models of growing disc galaxies with thick and thin discs. Similar to thin-disc-only models, galaxies with appropriate non-axisymmetric structures reproduce observational constraints on radial disc heating in and migration to the Solar Neighbourhood (Snhd). The presence of thick discs can suppress non-axisymmetries and thus higher baryonic-to-…
▽ More
We analyse disc heating and radial migration in N-body models of growing disc galaxies with thick and thin discs. Similar to thin-disc-only models, galaxies with appropriate non-axisymmetric structures reproduce observational constraints on radial disc heating in and migration to the Solar Neighbourhood (Snhd). The presence of thick discs can suppress non-axisymmetries and thus higher baryonic-to-dark matter fractions are required than in models that only have a thin disc. Models that are baryon-dominated to roughly the Solar radius R_0 are favoured, in agreement with data for the Milky Way. For inside-out growing discs, today's thick-disc stars at R_0 are dominated by outwards migrators. Whether outwards migrators are vertically hotter than non-migrators depends on the radial gradient of the thick disc vertical velocity dispersion. There is an effective upper boundary in angular momentum that thick disc stars born in the centre of a galaxy can reach by migration, which explains the fading of the high-alpha sequence outside R_0. Our models compare well to Snhd kinematics from RAVE-TGAS. For such comparisons it is important to take into account the azimuthal variation of kinematics at R ~ R_0 and biases from survey selection functions. The vertical heating of thin disc stars by giant molecular clouds is only mildly affected by the presence of thick discs. Our models predict higher vertical velocity dispersions for the oldest stars than found in the Snhd age-velocity dispersion relation, possibly because of measurement uncertainties or an underestimation of the number of old cold stars in our models.
△ Less
Submitted 13 June, 2017;
originally announced June 2017.
-
Assessing distances and consistency of kinematics in Gaia/TGAS
Authors:
Ralph Schönrich,
Michael Aumer
Abstract:
We apply the statistical methods by Schoenrich, Binney & Asplund to assess the quality of distances and kinematics in the RAVE-TGAS and LAMOST-TGAS samples of Solar neighbourhood stars. These methods yield a nominal distance accuracy of 1-2%. Other than common tests on parallax accuracy, they directly test distance estimations including the effects of distance priors. We show how to construct thes…
▽ More
We apply the statistical methods by Schoenrich, Binney & Asplund to assess the quality of distances and kinematics in the RAVE-TGAS and LAMOST-TGAS samples of Solar neighbourhood stars. These methods yield a nominal distance accuracy of 1-2%. Other than common tests on parallax accuracy, they directly test distance estimations including the effects of distance priors. We show how to construct these priors including the survey selection functions (SSFs) directly from the data. We demonstrate that neglecting the SSFs causes severe distance biases. Due to the decline of the SSFs in distance, the simple 1/parallax estimate only mildly underestimates distances. We test the accuracy of measured line-of-sight velocities (v_los) by binning the samples in the nominal v_los uncertainties. We find: a) the LAMOST v_los have a ~ -5 km/s offset; b) the average LAMOST measurement error for v_los is ~7 km/s, significantly smaller than, and nearly uncorrelated with the nominal LAMOST estimates. The RAVE sample shows either a moderate distance underestimate, or an unaccounted source of v_los dispersion (e_v) from measurement errors and binary stars. For a subsample of suspected binary stars in RAVE, our methods indicate significant distance underestimates. Separating a sample in metallicity or kinematics to select thick-disc/halo stars, discriminates between distance bias and e_v. For LAMOST, this separation yields consistency with pure v_los measurement errors. We find an anomaly near longitude l~(300+/-60)deg and distance s~(0.32+/-0.03)kpc on both sides of the galactic plane, which could be explained by either a localised distance error or a breathing mode.
△ Less
Submitted 30 August, 2017; v1 submitted 5 April, 2017;
originally announced April 2017.
-
Non-local thermodynamic equilibrium stellar spectroscopy with 1D and 3D models - II. Chemical properties of the Galactic metal-poor disc and the halo
Authors:
Maria Bergemann,
Remo Collet,
Ralph Schoenrich,
Rene Andrae,
Mikhail Kovalev,
Gregory Ruchti,
Camilla J. Hansen,
Zazralt Magic
Abstract:
From exploratory studies and theoretical expectations it is known that simplifying approximations in spectroscopic analysis (LTE, 1D) lead to systematic biases of stellar parameters and abundances. These biases depend strongly on surface gravity, temperature, and, in particular, for LTE vs. non-LTE (NLTE) on metallicity of the stars. Here we analyse the [Mg/Fe] and [Fe/H] plane of a sample of 326…
▽ More
From exploratory studies and theoretical expectations it is known that simplifying approximations in spectroscopic analysis (LTE, 1D) lead to systematic biases of stellar parameters and abundances. These biases depend strongly on surface gravity, temperature, and, in particular, for LTE vs. non-LTE (NLTE) on metallicity of the stars. Here we analyse the [Mg/Fe] and [Fe/H] plane of a sample of 326 stars, comparing LTE and NLTE results obtained using 1D hydrostatic models and averaged <3D> models. We show that compared to the <3D>NLTE benchmark, all other three methods display increasing biases towards lower metallicities, resulting in false trends of [Mg/Fe] against [Fe/H], which have profound implications for interpretations by chemical evolution models. In our best <3D> NLTE model, the halo and disc stars show a clearer behaviour in the [Mg/Fe] - [Fe/H] plane, from the knee in abundance space down to the lowest metallicities. Our sample has a large fraction of thick disc stars and this population extends down to at least [Fe/H] ~ -1.6 dex, further than previously proven. The thick disc stars display a constant [Mg/Fe] ~ 0.3 dex, with a small intrinsic dispersion in [Mg/Fe] that suggests that a fast SN Ia channel is not relevant for the disc formation. The halo stars reach higher [Mg/Fe] ratios and display a net trend of [Mg/Fe] at low metallicities, paired with a large dispersion in [Mg/Fe]. These indicate the diverse origin of halo stars from accreted low-mass systems to stochastic/inhomogeneous chemical evolution in the Galactic halo.
△ Less
Submitted 16 August, 2017; v1 submitted 21 December, 2016;
originally announced December 2016.
-
The Thirteenth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-IV Survey MApping Nearby Galaxies at Apache Point Observatory
Authors:
SDSS Collaboration,
Franco D. Albareti,
Carlos Allende Prieto,
Andres Almeida,
Friedrich Anders,
Scott Anderson,
Brett H. Andrews,
Alfonso Aragon-Salamanca,
Maria Argudo-Fernandez,
Eric Armengaud,
Eric Aubourg,
Vladimir Avila-Reese,
Carles Badenes,
Stephen Bailey,
Beatriz Barbuy,
Kat Barger,
Jorge Barrera-Ballesteros,
Curtis Bartosz,
Sarbani Basu,
Dominic Bates,
Giuseppina Battaglia,
Falk Baumgarten,
Julien Baur,
Julian Bautista,
Timothy C. Beers
, et al. (314 additional authors not shown)
Abstract:
The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases,…
▽ More
The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.
△ Less
Submitted 25 September, 2017; v1 submitted 5 August, 2016;
originally announced August 2016.
-
Age velocity dispersion relations and heating histories in disc galaxies
Authors:
Michael Aumer,
James Binney,
Ralph Schönrich
Abstract:
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 f…
▽ More
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.
△ Less
Submitted 7 July, 2016;
originally announced July 2016.
-
The Gaia-ESO survey: Hydrogen lines in red giants directly trace stellar mass
Authors:
Maria Bergemann,
Aldo Serenelli,
Ralph Schoenrich,
Greg Ruchti,
Andreas Korn,
Saskia Hekker,
Mikhail Kovalev,
Lyudmila Mashonkina,
Gerry Gilmore,
Sofia Randich,
Martin Asplund,
Hans-Walter Rix,
Andrew R. Casey,
Paula Jofre,
Elena Pancino,
Alejandra Recio-Blanco,
Patrick de Laverny,
Rodolfo Smiljanic,
Grazina Tautvaisiene,
Amelia Bayo,
Jim Lewis,
Sergey Koposov,
Anna Hourihane,
Clare Worley,
Lorenzo Morbidelli
, et al. (5 additional authors not shown)
Abstract:
Red giant stars are perhaps the most important type of stars for Galactic and extra-galactic archaeology: they are luminous, occur in all stellar populations, and their surface temperatures allow precise abundance determinations for many different chemical elements. Yet, the full star formation and enrichment history of a galaxy can be traced directly only if two key observables can be determined…
▽ More
Red giant stars are perhaps the most important type of stars for Galactic and extra-galactic archaeology: they are luminous, occur in all stellar populations, and their surface temperatures allow precise abundance determinations for many different chemical elements. Yet, the full star formation and enrichment history of a galaxy can be traced directly only if two key observables can be determined for large stellar samples - age and chemical composition. While spectroscopy is a powerful method to analyse the detailed abundances of stars, stellar ages are the "missing link in the chain", since they are not a direct observable. However, spectroscopy should be able to estimate stellar masses, which for red giants directly infer ages provided their chemical composition is known.
Here we establish a new empirical relation between the shape of the hydrogen line in the observed spectra of red giants and stellar mass determined from asteroseismology. The relation allows to determine stellar masses and ages with the accuracy of 10-15%. The method can be used with confidence for stars in the following range of stellar parameters: 4000 < Teff < 5000 K, 0.5 < log g < 3.5, -2.0 < [Fe/H] < 0.3, and luminosities log L/LSun < 2.5. Our analysis provides observational evidence that the Halpha spectral characteristics of red giant stars are tightly correlated with their mass and therefore their age. We also show that the method samples well all stellar populations with ages above 1 Gyr. Targeting bright giants, the method allows to obtain simultaneous age and chemical abundance information far deeper than would be possible with asteroseismology, extending the possible survey volume to remote regions of the Milky Way and even to neighbouring galaxies like Andromeda or the Magellanic Clouds already with present instrumentation, like VLT and Keck facilities.
△ Less
Submitted 17 June, 2016;
originally announced June 2016.
-
Chemical tagging with APOGEE: Discovery of a large population of N-rich stars in the inner Galaxy
Authors:
Ricardo P. Schiavon,
Olga Zamora,
Ricardo Carrera,
Sara Lucatello,
A. C. Robin,
Melissa Ness,
Sarah L. Martell,
Verne V. Smith,
D. A. Garcia Hernandez,
Arturo Manchado,
Ralph Schoenrich,
Nate Bastian,
Cristina Chiappini,
Matthew Shetrone,
J. Ted Mackereth,
Rob A. Williams,
Szabolcs Meszaros,
Carlos Allende Prieto,
Friedrich Anders,
Dmitry Bizyaev,
Timothy C. Beers,
S. Drew Chojnowski,
Katia Cunha,
Courtney Epstein,
Peter M. Frinchaboy
, et al. (20 additional authors not shown)
Abstract:
Formation of globular clusters (GCs), the Galactic bulge, or galaxy bulges in general, are important unsolved problems in Galactic astronomy. Homogeneous infrared observations of large samples of stars belonging to GCs and the Galactic bulge field are one of the best ways to study these problems. We report the discovery by APOGEE of a population of field stars in the inner Galaxy with abundances o…
▽ More
Formation of globular clusters (GCs), the Galactic bulge, or galaxy bulges in general, are important unsolved problems in Galactic astronomy. Homogeneous infrared observations of large samples of stars belonging to GCs and the Galactic bulge field are one of the best ways to study these problems. We report the discovery by APOGEE of a population of field stars in the inner Galaxy with abundances of N, C, and Al that are typically found in GC stars. The newly discovered stars have high [N/Fe], which is correlated with [Al/Fe] and anti-correlated with [C/Fe]. They are homogeneously distributed across, and kinematically indistinguishable from, other field stars in the same volume. Their metallicity distribution is seemingly unimodal, peaking at [Fe/H]~-1, thus being in disagreement with that of the Galactic GC system. Our results can be understood in terms of different scenarios. N-rich stars could be former members of dissolved GCs, in which case the mass in destroyed GCs exceeds that of the surviving GC system by a factor of ~8. In that scenario, the total mass contained in so-called "first-generation" stars cannot be larger than that in "second-generation" stars by more than a factor of ~9 and was certainly smaller. Conversely, our results may imply the absence of a mandatory genetic link between "second generation" stars and GCs. Last, but not least, N-rich stars could be the oldest stars in the Galaxy, the by-products of chemical enrichment by the first stellar generations formed in the heart of the Galaxy.
△ Less
Submitted 17 June, 2016;
originally announced June 2016.
-
Understanding inverse metallicity gradients in galactic discs as a consequence of inside-out formation
Authors:
Ralph Schönrich,
Paul McMillan
Abstract:
The early stages of a galaxy's evolution leave an imprint on its metallicity distribution. We discuss the origins and evolution of radial metallicity gradients in discs of spiral galaxies using an analytical chemical evolution model. We explain how radial metallicity gradients in stellar populations are determined by three factors: the radial metallicity profile of the star-forming ISM, radial cha…
▽ More
The early stages of a galaxy's evolution leave an imprint on its metallicity distribution. We discuss the origins and evolution of radial metallicity gradients in discs of spiral galaxies using an analytical chemical evolution model. We explain how radial metallicity gradients in stellar populations are determined by three factors: the radial metallicity profile of the star-forming ISM, radial changes in the star-formation history (in particular inside-out formation), and radial mixing of stars. Under reasonable assumptions, inside-out formation steepens the negative ISM metallicity gradient, but contributes positively to the stellar metallicity gradient, up to inverting the metallicity profile to a positive d[Fe/H]/dR. This reconciles steep negative d[Fe/H]/dR in some high redshift galaxies to generally flatter gradients in local observations.
We discuss the evidence for inverse radial metallicity gradients (positive d[X/H]/dR) at high redshifts and the inverse relationship between azimuthal velocity and the metallicity (positive dV_φ/d[Fe/H]) of stars for the Milky Way's thick disc. The former can be achieved by high central gas-loss rates and re-distribution processes, e.g. re-accretion of enriched material in conjunction with inside-out formation, and near-disc galactic fountaining. For the Milky Way thick disc, we show that the positive dV_φ/d[Fe/H] correlation points to comparable timescales for inside-out formation, initial metal enrichment and SNIa enrichment. We argue that the original ISM metallicity gradient could be inferred with better data from the high-metallicity tail of the alpha enhanced population. Including inside-out formation in our models changes the local vertical metallicity gradient by about -0.2dex/kpc, in line with local measurements.
△ Less
Submitted 5 April, 2017; v1 submitted 8 May, 2016;
originally announced May 2016.
-
Inflow, Outflow, Yields, and Stellar Population Mixing in Chemical Evolution Models
Authors:
Brett H. Andrews,
David H. Weinberg,
Ralph Schönrich,
Jennifer A. Johnson
Abstract:
Chemical evolution models are powerful tools for interpreting stellar abundance surveys and understanding galaxy evolution. However, their predictions depend heavily on the treatment of inflow, outflow, star formation efficiency (SFE), the stellar initial mass function, the Type Ia supernova delay time distribution, stellar yields, and stellar population mixing. Using flexCE, a flexible one-zone c…
▽ More
Chemical evolution models are powerful tools for interpreting stellar abundance surveys and understanding galaxy evolution. However, their predictions depend heavily on the treatment of inflow, outflow, star formation efficiency (SFE), the stellar initial mass function, the Type Ia supernova delay time distribution, stellar yields, and stellar population mixing. Using flexCE, a flexible one-zone chemical evolution code, we investigate the effects of and trade-offs between parameters. Two critical parameters are SFE and the outflow mass-loading parameter, which shift the knee in [O/Fe]-[Fe/H] and the equilibrium abundances that the simulations asymptotically approach, respectively. One-zone models with simple star formation histories follow narrow tracks in [O/Fe]-[Fe/H] unlike the observed bimodality (separate high-alpha and low-alpha sequences) in this plane. A mix of one-zone models with inflow timescale and outflow mass-loading parameter variations, motivated by the inside-out galaxy formation scenario with radial mixing, reproduces the two sequences better than a one-zone model with two infall epochs. We present [X/Fe]-[Fe/H] tracks for 20 elements assuming three different supernova yield models and find some significant discrepancies with solar neighborhood observations, especially for elements with strongly metallicity-dependent yields. We apply principal component abundance analysis (PCAA) to the simulations and existing data to reveal the main correlations amongst abundances and quantify their contributions to variation in abundance space. For the stellar population mixing scenario, the abundances of alpha-elements and elements with metallicity-dependent yields dominate the first and second principal components, respectively, and collectively explain 99% of the variance in the model. flexCE is a python package available at https://github.com/bretthandrews/flexCE.
△ Less
Submitted 21 December, 2016; v1 submitted 28 April, 2016;
originally announced April 2016.
-
The quiescent phase of galactic disc growth
Authors:
Michael Aumer,
James Binney,
Ralph Schönrich
Abstract:
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 s…
▽ More
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.
△ Less
Submitted 1 April, 2016;
originally announced April 2016.
-
Measuring the vertical age structure of the Galactic disc using asteroseismology and SAGA
Authors:
L. Casagrande,
V. Silva Aguirre,
K. J. Schlesinger,
D. Stello,
D. Huber,
A. M. Serenelli,
R. Schoenrich,
S. Cassisi,
A. Pietrinferni,
S. Hodgkin,
A. P. Milone,
S. Feltzing,
M. Asplund
Abstract:
The existence of a vertical age gradient in the Milky Way disc has been indirectly known for long. Here, we measure it directly for the first time with seismic ages, using red giants observed by Kepler. We use Stroemgren photometry to gauge the selection function of asteroseismic targets, and derive colour and magnitude limits where giants with measured oscillations are representative of the under…
▽ More
The existence of a vertical age gradient in the Milky Way disc has been indirectly known for long. Here, we measure it directly for the first time with seismic ages, using red giants observed by Kepler. We use Stroemgren photometry to gauge the selection function of asteroseismic targets, and derive colour and magnitude limits where giants with measured oscillations are representative of the underlying population in the field. Limits in the 2MASS system are also derived. We lay out a method to assess and correct for target selection effects independent of Galaxy models. We find that low mass, i.e. old red giants dominate at increasing Galactic heights, whereas closer to the Galactic plane they exhibit a wide range of ages and metallicities. Parametrizing this as a vertical gradient returns approximately 4 Gyr/kpc for the disc we probe, although with a large dispersion of ages at all heights. The ages of stars show a smooth distribution over the last 10 Gyr, consistent with a mostly quiescent evolution for the Milky Way disc since a redshift of about 2. We also find a flat age-metallicity relation for disc stars. Finally, we show how to use secondary clump stars to estimate the present-day intrinsic metallicity spread, and suggest using their number count as a new proxy for tracing the ageing of the disc. This work highlights the power of asteroseismology for Galactic studies; however, we also emphasize the need for better constraints on stellar mass-loss, which is a major source of systematic age uncertainties in red giant stars.
△ Less
Submitted 5 October, 2015;
originally announced October 2015.
-
Kinematic Detection of the Galactic Nuclear Disc
Authors:
Ralph Schönrich,
Michael Aumer,
Stuart E. Sale
Abstract:
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…
▽ More
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.
△ Less
Submitted 5 October, 2015; v1 submitted 9 July, 2015;
originally announced July 2015.
-
Origin of the high v_los feature in the Galactic bar
Authors:
Michael Aumer,
Ralph Schönrich
Abstract:
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 obtaine…
▽ More
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.
△ Less
Submitted 5 October, 2015; v1 submitted 3 July, 2015;
originally announced July 2015.
-
Rings and Radial Waves in the Disk of the Milky Way
Authors:
Yan Xu,
Heidi Jo Newberg,
Jeffrey L. Carlin,
Chao Liu,
Licai Deng,
Jing Li,
Ralph Schönrich,
Brian Yanny
Abstract:
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…
▽ More
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.
△ Less
Submitted 1 March, 2015;
originally announced March 2015.
-
The Vertical Metallicity Gradient of the Milky Way Disk: Transitions in [a/Fe] Populations
Authors:
Katharine J. Schlesinger,
Jennifer A. Johnson,
Constance M. Rockosi,
Young Sun Lee,
Timothy C. Beers,
Paul Harding,
Carlos Allende Prieto,
Jonathan C. Bird,
Ralph Schoenrich,
Brian Yanny,
Donald P. Schneider,
Benjamin A. Weaver,
Jon Brinkmann
Abstract:
Using G dwarfs from the Sloan Extension for Galactic Understanding and Exploration (SEGUE) survey, we have determined a vertical metallicity gradient over a large volume of the Milky Way's disk, and examined how this gradient varies for different [a/Fe] subsamples. This sample contains over 40,000 stars with low-resolution spectroscopy over 144 lines of sight. We employ the SEGUE Stellar Parameter…
▽ More
Using G dwarfs from the Sloan Extension for Galactic Understanding and Exploration (SEGUE) survey, we have determined a vertical metallicity gradient over a large volume of the Milky Way's disk, and examined how this gradient varies for different [a/Fe] subsamples. This sample contains over 40,000 stars with low-resolution spectroscopy over 144 lines of sight. We employ the SEGUE Stellar Parameter Pipeline (SSPP) to obtain estimates of effective temperature, surface gravity, [Fe/H], and [a/Fe] for each star and extract multiple volume-complete subsamples of approximately 1000 stars each. Based on the survey's consistent target-selection algorithm, we adjust each subsample to determine an unbiased picture of the disk in [Fe/H] and [a/Fe]; consequently, each individual star represents the properties of many. The SEGUE sample allows us to constrain the vertical metallicity gradient for a large number of stars over a significant volume of the disk, between ~0.3 and 1.6 kpc from the Galactic plane, and examine the in situ structure, in contrast to previous analyses which are more limited in scope. This work does not pre-suppose a disk structure, whether composed of a single complex population or a distinct thin and thick disk component. The metallicity gradient is -0.243 +0.039 -0.053 dex/kpc for the sample as a whole, which we compare to various literature results. Each [a/Fe] subsample dominates at a different range of heights above the plane of the Galaxy, which is exhibited in the gradient found in the sample as a whole. Stars over a limited range in [a/Fe] show little change in median [Fe/H] with height. If we associate [a/Fe] with age, our consistent vertical metallicity gradients with [a/Fe] suggest that stars formed in different epochs exhibit comparable vertical structure, implying similar star-formation processes and evolution.
△ Less
Submitted 26 May, 2014;
originally announced May 2014.
-
Stroemgren survey for Asteroseismology and Galactic Archaeology: let the SAGA begin
Authors:
L. Casagrande,
V. Silva Aguirre,
D. Stello,
D. Huber,
A. M. Serenelli,
S. Cassisi,
A. Dotter,
A. P. Milone,
S. Hodgkin,
A. F. Marino,
M. N. Lund,
A. Pietrinferni,
M. Asplund,
S. Feltzing,
C. Flynn,
F. Grundahl,
P. E. Nissen,
R. Schoenrich,
K. J. Schlesinger,
W. Wang
Abstract:
Asteroseismology has the capability of precisely determining stellar properties which would otherwise be inaccessible, such as radii, masses and thus ages of stars. When coupling this information with classical determinations of stellar parameters, such as metallicities, effective temperatures and angular diameters, powerful new diagnostics for Galactic studies can be obtained. The ongoing Stroemg…
▽ More
Asteroseismology has the capability of precisely determining stellar properties which would otherwise be inaccessible, such as radii, masses and thus ages of stars. When coupling this information with classical determinations of stellar parameters, such as metallicities, effective temperatures and angular diameters, powerful new diagnostics for Galactic studies can be obtained. The ongoing Stroemgren survey for Asteroseismology and Galactic Archaeology (SAGA) has the goal of transforming the Kepler field into a new benchmark for Galactic studies, similarly to the solar neighborhood. Here we present first results from a stripe centred at Galactic longitude 74deg and covering latitude from about 8 to 20deg, which includes almost 1000 K-giants with seismic information and the benchmark open cluster NGC 6819. We describe the coupling of classical and seismic parameters, the accuracy as well as the caveats of the derived effective temperatures, metallicities, distances, surface gravities, masses, and radii. Confidence in the achieved precision is corroborated by the detection of the first and secondary clump in a population of field stars with a ratio of 2 to 1, and by the negligible scatter in the seismic distances among NGC 6819 member stars. An assessment of the reliability of stellar parameters in the Kepler Input Catalogue is also performed, and the impact of our results for population studies in the Milky Way is discussed, along with the importance of an all-sky Stroemgren survey.
△ Less
Submitted 11 March, 2014;
originally announced March 2014.
-
Testing the Asteroseismic Mass Scale Using Metal-Poor Stars Characterized with APOGEE and Kepler
Authors:
Courtney R. Epstein,
Yvonne P. Elsworth,
Jennifer A. Johnson,
Matthew Shetrone,
Benoît Mosser,
Saskia Hekker,
Jamie Tayar,
Paul Harding,
Marc Pinsonneault,
Víctor Silva Aguirre,
Sarbani Basu,
Timothy C. Beers,
Dmitry Bizyaev,
Timothy R. Bedding,
William J. Chaplin,
Peter M. Frinchaboy,
Rafael A. García,
Ana E. García Pérez,
Fred R. Hearty,
Daniel Huber,
Inese I. Ivans,
Steven R. Majewski,
Savita Mathur,
David Nidever,
Aldo Serenelli
, et al. (7 additional authors not shown)
Abstract:
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 rel…
▽ More
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.
△ Less
Submitted 7 March, 2014;
originally announced March 2014.
-
Does SEGUE/SDSS indicate a dual Galactic halo?
Authors:
Ralph Schoenrich,
Martin Asplund,
Luca Casagrande
Abstract:
We re-examine recent claims of observational evidence for a dual Galactic halo in SEGUE/SDSS data, and trace them back to improper error treatment and neglect of selection effects. In particular, the detection of a vertical abundance gradient in the halo can be explained as a metallicity bias in distance. A similar bias, and the impact of disk contamination, affect the sample of blue horizontal br…
▽ More
We re-examine recent claims of observational evidence for a dual Galactic halo in SEGUE/SDSS data, and trace them back to improper error treatment and neglect of selection effects. In particular, the detection of a vertical abundance gradient in the halo can be explained as a metallicity bias in distance. A similar bias, and the impact of disk contamination, affect the sample of blue horizontal branch stars. These examples highlight why non-volume complete samples require forward-modelling from theoretical models or extensive bias-corrections. We also show how observational uncertainties produce the specific non-Gaussianity in the observed azimuthal velocity distribution of halo stars, which can be erroneously identified as two Gaussian components. A single kinematic component yields an excellent fit to the observed data, when we model the measurement process including distance uncertainties. Furthermore, we show that sample differences in proper motion space are the direct consequence of kinematic cuts, and are enhanced when distance estimates are less accurate. Thus, their presence is neither a proof for a separate population, nor a measure of reliability for the applied distances. We conclude that currently there is no evidence from SEGUE/SDSS that would favour a dual Galactic halo over a single halo full of substructure.
△ Less
Submitted 4 March, 2014;
originally announced March 2014.
-
Fundamental stellar parameters and metallicities from Bayesian spectroscopy. Application to low- and high-resolution spectra
Authors:
Ralph Schönrich,
Maria Bergemann
Abstract:
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 t…
▽ More
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.
△ Less
Submitted 11 June, 2014; v1 submitted 21 November, 2013;
originally announced November 2013.