Researcher profile

Sanjib Sharma

Sanjib Sharma contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
26works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

26 published item(s)

preprint2022arXiv

Discovery of post-mass-transfer helium-burning red giants using asteroseismology

A star expands to become a red giant when it has fused all the hydrogen in its core into helium. If the star is in a binary system, its envelope can overflow onto its companion or be ejected into space, leaving a hot core and potentially forming a subdwarf-B star. However, most red giants that have partially transferred envelopes in this way remain cool on the surface and are almost indistinguishable from those that have not. Among $\sim$7000 helium-burning red giants observed by NASA's Kepler mission, we use asteroseismology to identify two classes of stars that must have undergone dramatic mass loss, presumably due to stripping in binary interactions. The first class comprises about 7 underluminous stars with smaller helium-burning cores than their single-star counterparts. Theoretical models show that these small cores imply the stars had much larger masses when ascending the red giant branch. The second class consists of 32 red giants with masses down to 0.5 M$_\odot$, whose implied ages would exceed the age of the universe had no mass loss occurred. The numbers are consistent with binary statistics, and our results open up new possibilities to study the evolution of post-mass-transfer binary systems.

preprint2022arXiv

Flood hazard model calibration using multiresolution model output

Riverine floods pose a considerable risk to many communities. Improving flood hazard projections has the potential to inform the design and implementation of flood risk management strategies. Current flood hazard projections are uncertain, especially due to uncertain model parameters. Calibration methods use observations to quantify model parameter uncertainty. With limited computational resources, researchers typically calibrate models using either relatively few expensive model runs at high spatial resolutions or many cheaper runs at lower spatial resolutions. This leads to an open question: Is it possible to effectively combine information from the high and low resolution model runs? We propose a Bayesian emulation-calibration approach that assimilates model outputs and observations at multiple resolutions. As a case study for a riverine community in Pennsylvania, we demonstrate our approach using the LISFLOOD-FP flood hazard model. The multiresolution approach results in improved parameter inference over the single resolution approach in multiple scenarios. Results vary based on the parameter values and the number of available models runs. Our method is general and can be used to calibrate other high dimensional computer models to improve projections.

preprint2022arXiv

Reliable stellar abundances of individual stars with the MUSE integral-field spectrograph

We present a novel approach to deriving stellar labels for stars observed in MUSE fields making use of data-driven machine learning methods. Taking advantage of the comparable spectral properties (resolution, wavelength coverage) of the LAMOST and MUSE instruments, we adopt the Data-Driven Payne (DD-Payne) model used on LAMOST observations and apply it to stars observed in MUSE fields. Remarkably, in spite of instrumental differences, according to the cross-validation of 27 LAMOST-MUSE common stars, we are able to determine stellar labels with precision better than 75K in $T_{\rm eff}$, 0.15 dex in $\log g$, and 0.1 dex in abundances of [Fe/H], [Mg/Fe], [Si/Fe], [Ti/Fe], [C/Fe], [Ni/Fe] and [Cr/Fe] for current MUSE observations over a parameter range of 3800<$T_{\rm eff}$<7000 K, -1.5<[Fe/H]<0.5 dex. To date, MUSE has been used to target 13,000 fields across the southern sky since it was first commissioned six years ago and it is unique in its ability to study dense star fields such as globular clusters or the Milky Way bulge. Our method will enable the automated determination of stellar parameters for all stars in these fields. Additionally, it opens the door for applications to data collected by other spectrographs having resolution similar to LAMOST. With the upcoming BlueMUSE and MAVIS, we will gain access to a whole new range of chemical abundances with higher precision, especially critical s-process elements such as [Y/Fe] and [Ba/Fe] that provide key age diagnostics for stellar targets.

preprint2022arXiv

The GALAH Survey: A New Sample of Extremely Metal-Poor Stars Using A Machine Learning Classification Algorithm

Extremely Metal-Poor (EMP) stars provide a valuable probe of early chemical enrichment in the Milky Way. Here we leverage a large sample of $\sim600,000$ high-resolution stellar spectra from the GALAH survey plus a machine learning algorithm to find 54 candidates with estimated [Fe/H]~$\leq$~-3.0, 6 of which have [Fe/H]~$\leq$~-3.5. Our sample includes $\sim 20 \%$ main sequence EMP candidates, unusually high for \emp surveys. We find the magnitude-limited metallicity distribution function of our sample is consistent with previous work that used more complex selection criteria. The method we present has significant potential for application to the next generation of massive stellar spectroscopic surveys, which will expand the available spectroscopic data well into the millions of stars.

preprint2022arXiv

The GALAH Survey: Chemical tagging and chrono-chemodynamics of accreted halo stars with GALAH+ DR3 and $Gaia$ eDR3

Since the advent of $Gaia$ astrometry, it is possible to identify massive accreted systems within the Galaxy through their unique dynamical signatures. One such system, $Gaia$-Sausage-Enceladus (GSE), appears to be an early &#34;building block&#34; given its virial mass $> 10^{10}\,\mathrm{M_\odot}$ at infall ($z\sim1-3$). In order to separate the progenitor population from the background stars, we investigate its chemical properties with up to 30 element abundances from the GALAH+ Survey Data Release 3 (DR3). To inform our choice of elements for purely chemically selecting accreted stars, we analyse 4164 stars with low-$α$ abundances and halo kinematics. These are most different to the Milky Way stars for abundances of Mg, Si, Na, Al, Mn, Fe, Ni, and Cu. Based on the significance of abundance differences and detection rates, we apply Gaussian mixture models to various element abundance combinations. We find the most populated and least contaminated component, which we confirm to represent GSE, contains 1049 stars selected via [Na/Fe] vs. [Mg/Mn] in GALAH+ DR3. We provide tables of our selections and report the chrono-chemodynamical properties (age, chemistry, and dynamics). Through a previously reported clean dynamical selection of GSE stars, including $30 < \sqrt{J_R~/~\mathrm{kpc\,km\,s^{-1}}} < 55$, we can characterise an unprecedented 24 abundances of this structure with GALAH+ DR3. Our chemical selection allows us to prevent circular reasoning and characterise the dynamical properties of the GSE, for example mean $\sqrt{J_R~/~\mathrm{kpc\,km\,s^{-1}}} = 26_{-14}^{+9}$. We find only $(29\pm1)\%$ of the GSE stars within the clean dynamical selection region. Our methodology will improve future studies of accreted structures and their importance for the formation of the Milky Way.

preprint2022arXiv

The Global Dynamical Atlas of the Milky Way mergers: Constraints from Gaia EDR3 based orbits of globular clusters, stellar streams and satellite galaxies

The Milky Way halo was predominantly formed by the merging of numerous progenitor galaxies. However, our knowledge of this process is still incomplete, especially in regard to the total number of mergers, their global dynamical properties and their contribution to the stellar population of the Galactic halo. Here, we uncover the Milky Way mergers by detecting groupings of globular clusters, stellar streams and satellite galaxies in action ($\mathbf{J}$) space. While actions fully characterize the orbits, we additionally use the redundant information on their energy ($\textit{E}$) to enhance the contrast between groupings. For this endeavour, we use $\textit{Gaia}$ EDR3 based measurements of $170$ globular clusters, $41$ streams and $46$ satellites to derive their $\mathbf{J}$ and $\textit{E}$. To detect groups, we use the $\texttt{ENLINK}$ software, coupled with a statistical procedure that accounts for the observed phase-space uncertainties of these objects. We detect a total of $N=6$ groups, including the previously known mergers $\textit{Sagittarius}$, $\textit{Cetus}$, $\textit{Gaia-Sausage/Enceladus}$, $\textit{LMS-1/Wukong}$, $\textit{Arjuna/Sequoia/I&#39;itoi}$ and one new merger that we call $\textit{Pontus}$. All of these mergers, together, comprise $62$ objects ($\approx 25\%$ of our sample). We discuss their members, orbital properties and metallicity distributions. We find that the three most metal-poor streams of our Galaxy -- &#34;C-19&#34; ([Fe/H]$=-3.4$ dex), &#34;Sylgr&#34; ([Fe/H]$=-2.9$ dex) and &#34;Phoenix&#34; ([Fe/H]$=-2.7$ dex) -- are associated with $\textit{LMS-1/Wukong}$; showing it to be the most metal-poor merger. The global dynamical atlas of Milky Way mergers that we present here provides a present-day reference for galaxy formation models.

preprint2022arXiv

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

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

preprint2022arXiv

Trade-offs in the design and communication of flood-risk information

There is an increasingly urgent need to develop knowledge and practices to manage climate risks. For example, flood-risk information can inform household decisions such as purchasing a home or flood insurance. However, flood-risk estimates are deeply uncertain, meaning that they are subject to sizeable disagreement. Available flood-risk estimates provide inconsistent and incomplete information and pose communication challenges. The effects of different choices of design and communication options can create confusion in decision-making processes. The climate services literature includes insights into desirable features for producing information that is credible and relevant. Using examples of riverine (fluvial) flood-risk information products and studies in the United States, we assess how existing risk characterizations integrate desirable features outlined in the climate services literature. Improved characterization and communication of decision-relevant (and often deep) uncertainties, including those arising from human decisions, is a crucial next step. We argue that producing relevant flood-risk information requires applying principles of open science and co-production.

preprint2021arXiv

A safety factor approach to designing urban infrastructure for dynamic conditions

Current approaches to design flood-sensitive infrastructure typically assume a stationary rainfall distribution and neglect many uncertainties. These assumptions are inconsistent with observations that suggest intensifying extreme precipitation events and the uncertainties surrounding projections of the coupled natural-human systems. Here we demonstrate a safety factor approach to designing urban infrastructure in a changing climate. Our results show that assuming climate stationarity and neglecting deep uncertainties can drastically underestimate flood risks and lead to poor infrastructure design choices. We find that climate uncertainty dominates the socioeconomic and engineering uncertainties that impact the hydraulic reliability in stormwater drainage systems. We quantify the upfront costs needed to achieve higher hydraulic reliability and robustness against the deep uncertainties surrounding projections of rainfall, surface runoff characteristics, and infrastructure lifetime. Depending on the location, we find that adding safety factors of 1.4 to 1.7 to the standard stormwater pipe design guidance produces robust performance to the considered deep uncertainties. The insights gained from this study highlight the need for updating traditional engineering design strategies to improve infrastructure reliability under socioeconomic and environmental changes.

preprint2021arXiv

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

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

preprint2021arXiv

The GALAH survey: tracing the Galactic disk with Open Clusters

Open clusters are unique tracers of the history of our own Galaxy&#39;s disk. According to our membership analysis based on \textit{Gaia} astrometry, out of the 226 potential clusters falling in the footprint of GALAH or APOGEE, we find that 205 have secure members that were observed by at least one of the survey. Furthermore, members of 134 clusters have high-quality spectroscopic data that we use to determine their chemical composition. We leverage this information to study the chemical distribution throughout the Galactic disk of 21 elements, from C to Eu. The radial metallicity gradient obtained from our analysis is $-$0.076$\pm$0.009 dex kpc$^{-1}$, which is in agreement with previous works based on smaller samples. Furthermore, the gradient in the [Fe/H] - guiding radius (r$_{\rm guid}$) plane is $-$0.073$\pm$0.008 dex kpc$^{-1}$. We show consistently that open clusters trace the distribution of chemical elements throughout the Galactic disk differently than field stars. In particular, at given radius, open clusters show an age-metallicity relation that has less scatter than field stars. As such scatter is often interpreted as an effect of radial migration, we suggest that these differences are due to the physical selection effect imposed by our Galaxy: clusters that would have migrated significantly also had higher chances to get destroyed. Finally, our results reveal trends in the [X/Fe]$-$r$_{\rm guid}$$-$age space, which are important to understand production rates of different elements as a function of space and time.

preprint2021arXiv

The K2 Galactic Archaeology Program: Overview, target selection and survey properties

K2 was a community-driven NASA mission where all targets were proposed through guest observer programs. Among those, the K2 Galactic Archaeology Program (K2GAP) was devoted to measuring asteroseismic signals from giant stars to inform studies of the Galaxy, with about 25% of the observed K2 targets being allocated to this program. Here, we provide an overview of this program. We discuss in detail the target selection procedure and provide a python code that implements the selection function, see https://github.com/sanjibs/k2gap . Additionally, we discuss the detection completeness of the asteroseismic parameters $ν_{\rm max}$ and $Δν$. Broadly speaking, the targets were selected based on 2MASS color $J-Ks > 0.5$, with finely tuned adjustments for each campaign. Making use of the selection function we compare the observed distribution of asteroseismic masses to theoretical predictions. The median asteroseismic mass is higher by about 4% compared to predictions. Additionally, the number of seismic detections is on average 14% lower than expected. We provide a selection-function-matched mock catalog of stars based on a synthetic model of the Galaxy for the community to be used in subsequent analyses of the K2GAP data set, see http://www.physics.usyd.edu.au/k2gap .

preprint2020arXiv

A Gaia early DR3 mock stellar catalog: Galactic prior and selection function

We present a mock stellar catalog, matching in volume, depth and data model the content of the planned Gaia early data release 3 (Gaia EDR3). We have generated our catalog (GeDR3mock) using galaxia, a tool to sample stars from an underlying Milky Way (MW) model or from N-body data. We used an updated Besançon Galactic model together with the latest PARSEC stellar evolutionary tracks, now also including white dwarfs. We added the Magellanic clouds and realistic open clusters with internal rotation. We empirically modelled uncertainties based on Gaia DR2 (GDR2) and scaled them according to the longer baseline in Gaia EDR3. The apparent magnitudes were reddened according to a new selection of 3D extinction maps. To help with the Gaia selection function we provide all-sky magnitude limit maps in G and BP for a few relevant GDR2 subsets together with the routines to produce these maps for user-defined subsets. We supplement the catalog with photometry and extinctions in non-Gaia bands. The catalog is available in the Virtual Observatory and can be queried just like the actual Gaia EDR3 will be. We highlight a few capabilities of the Astronomy Data Query Language (ADQL) with educative catalog queries. We use the data extracted from those queries to compare GeDR3mock to GDR2, which emphasises the importance of adding observational noise to the mock data. Since the underlying truth, e.g. stellar parameters, is know in GeDR3mock, it can be used to construct priors as well as mock data tests for parameter estimation. All code, models and data used to produce GeDR3mock are linked and contained in galaxia_wrap, a python package, representing a fast galactic forward model, able to project MW models and N-body data into realistic Gaia observables.

preprint2020arXiv

Abundances in the Milky Way across five nucleosynthetic channels from 4 million LAMOST stars

Large stellar surveys are revealing the chemodynamical structure of the Galaxy across a vast spatial extent. However, the many millions of low-resolution spectra observed to date are yet to be fully exploited. We employ The Cannon, a data-driven approach to estimating abundances, to obtain detailed abundances from low-resolution (R = 1800) LAMOST spectra, using the GALAH survey as our reference. We deliver five (for dwarfs) or six (for giants) estimated abundances representing five different nucleosynthetic channels, for 3.9 million stars, to a precision of 0.05 - 0.23 dex. Using wide binary pairs, we demonstrate that our abundance estimates provide chemical discriminating power beyond metallicity alone. We show the coverage of our catalogue with radial, azimuthal and dynamical abundance maps, and examine the neutron capture abundances across the disk and halo, which indicate different origins for the in-situ and accreted halo populations. LAMOST has near-complete Gaia coverage and provides an unprecedented perspective on chemistry across the Milky Way.

preprint2020arXiv

K2-HERMES II. Planet-candidate properties from K2 Campaigns 1-13

Accurate and precise radius estimates of transiting exoplanets are critical for understanding their compositions and formation mechanisms. To know the planet, we must know the host star in as much detail as possible. We present complete results for planet-candidate hosts from the K2-HERMES survey, which uses the HERMES multi-object spectrograph on the Anglo-Australian Telescope to obtain R$\sim$28,000 spectra for more than 30,000 K2 stars. We present complete host-star parameters and planet-candidate radii for 224 K2 candidate planets from C1-C13. Our results cast severe doubt on 30 K2 candidates, as we derive unphysically large radii, larger than 2 $R_{Jup}$. This work highlights the importance of obtaining accurate, precise, and self-consistent stellar parameters for ongoing large planet search programs - something that will only become more important in the coming years, as TESS begins to deliver its own harvest of exoplanets.

preprint2020arXiv

Milky Way Tomography with the SkyMapper Southern Survey. II. Photometric Re-calibration of SMSS DR2

We apply the spectroscopy-based stellar-color regression (SCR) method to perform an accurate photometric re-calibration of the second data release from the SkyMapper Southern Survey (SMSS DR2). From comparison with a sample of over 200,000 dwarf stars with stellar atmospheric parameters taken from GALAH+ DR3 and with accurate, homogeneous photometry from $Gaia$ DR2, zero-point offsets are detected in the original photometric catalog of SMSS DR2, in particular for the gravity- and metallicity-sensitive $uv$ bands. For $uv$ bands, the zero-point offsets are close to zero at very low extinction, and then steadily increase with $E (B - V)$, reaching as large as 0.174 and 0.134 mag respectively, at $E (B - V) \sim 0.5$ mag. These offsets largely arise from the adopted dust term in the transformations used by SMSS DR2 to construct photometric calibrators from the ATLAS reference catalog. For the $gr$ bands, the zero-point offsets exhibit negligible variations with SFD $E(B - V )$, due to their tiny coefficients on the dust term in the transformation. Our study also reveals small, but significant, spatial variations of the zero-point offsets in all $uvgr$ bands. External checks using Strömgren photometry, WD loci and the SDSS Stripe 82 standard-star catalog independently confirm the zero-points found by our revised SCR method.

preprint2020arXiv

Synthetic Gaia surveys from the FIRE cosmological simulations of Milky Way-mass galaxies

With Gaia Data Release 2, the astronomical community is entering a new era of multidimensional surveys of the Milky Way. This new phase-space view of our Galaxy demands new tools for comparing observations to simulations of Milky-Way-mass galaxies in a cosmological context, to test the physics of both dark matter and galaxy formation. We present ananke, a framework for generating synthetic phase-space surveys from high-resolution baryonic simulations, and use it to generate a suite of synthetic surveys resembling Gaia DR2 in data structure, magnitude limits, and observational errors. We use three cosmological simulations of Milky-Way-mass galaxies from the Latte suite of the Feedback In Realistic Environments (FIRE) project, which feature self-consistent clustering of star formation in dense molecular clouds and thin stellar/gaseous disks in live cosmological halos with satellite dwarf galaxies and stellar halos. We select three solar viewpoints from each simulation to generate nine synthetic Gaia-like surveys. We sample synthetic stars by assuming each star particle (of mass 7070 $M_{\odot}$) represents a single stellar population. At each viewpoint, we compute dust extinction from the simulated gas metallicity distribution and apply a simple error model to produce a synthetic Gaia-like survey that includes both observational properties and a pointer to the generating star particle. We provide the complete simulation snapshot at $z = 0$ for each simulated galaxy. We describe data access points, the data model, and plans for future upgrades. These synthetic surveys provide a tool for the scientific community to test analysis methods and interpret Gaia data.

preprint2020arXiv

Testing the intrinsic scatter of the asteroseismic scaling relations with Kepler red giants

Asteroseismic scaling relations are often used to derive stellar masses and radii, particulaly for stellar, exoplanet, and Galactic studies. It is therefore important that their precisions are known. Here we measure the intrinsic scatter of the underlying seismic scaling relations for $Δν$ and $ν_{\rm max}$, using two sharp features that are formed in the H--R diagram (or related diagrams) by the red giant populations. These features are the edge near the zero-age core-helium-burning phase, and the strong clustering of stars at the so-called red giant branch bump. The broadening of those features is determined by factors including the intrinsic scatter of the scaling relations themselves, and therefore it is capable of imposing constraints on them. We modelled Kepler stars with a Galaxia synthetic population, upon which we applied the intrinsic scatter of the scaling relations to match the degree of sharpness seen in the observation. We found that the random errors from measuring $Δν$ and $ν_{\rm max}$ provide the dominating scatter that blurs the features. As a consequence, we conclude that the scaling relations have intrinsic scatter of $\sim0.5\%$ ($Δν$), $\sim1.1\%$ ($ν_{\rm max}$), $\sim1.7\%$ ($M$) and $\sim0.4\%$ ($R$), for the SYD pipeline measured $Δν$ and $ν_{\rm max}$. This confirms that the scaling relations are very powerful tools. In addition, we show that standard evolution models fail to predict some of the structures in the observed population of both the HeB and RGB stars. Further stellar model improvements are needed to reproduce the exact distributions.

preprint2020arXiv

The GALAH Survey: A new constraint on cosmological lithium and Galactic lithium evolution from warm dwarf stars

Lithium depletion and enrichment in the cosmos is not yet well understood. To help tighten constraints on stellar and Galactic evolution models, we present the largest high-resolution analysis of Li abundances A(Li) to date, with results for over 100 000 GALAH field stars spanning effective temperatures $5900\,\mathrm{K} \lesssim \rm{T_{eff}} \lesssim7000\,\mathrm{K}$ and metallicities $-3 \lesssim \rm[Fe/H] \lesssim +0.5$. We separated these stars into two groups, on the warm and cool side of the so-called Li-dip, a localised region of the Kiel diagram wherein lithium is severely depleted. We discovered that stars in these two groups show similar trends in the A(Li)-[Fe/H] plane, but with a roughly constant offset in A(Li) of 0.4 dex, the warm group having higher Li abundances. At $\rm[Fe/H]\gtrsim-0.5$, a significant increasing in Li abundance with increasing metallicity is evident in both groups, signalling the onset of significant Galactic production. At lower metallicity, stars in the cool group sit on the Spite plateau, showing a reduced lithium of around 0.4 dex relative to the primordial value predicted from Big Bang nucleosynthesis (BBN). However, stars in the warm group between [Fe/H] = -1.0 and -0.5, form an elevated plateau that is largely consistent with the BBN prediction. This may indicate that these stars in fact preserve the primordial Li produced in the early Universe.

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

The tidal remnant of an unusually metal-poor globular cluster

Globular clusters are some of the oldest bound stellar structures observed in the Universe. They are ubiquitous in large galaxies and are believed to trace intense star formation events and the hierarchical build-up of structure. Observations of globular clusters in the Milky Way, and a wide variety of other galaxies, have found evidence for a `metallicity floor&#39;, whereby no globular clusters are found with chemical (`metal&#39;) abundances below approximately 0.3 to 0.4 per cent of that of the Sun. The existence of this metallicity floor may reflect a minimum mass and a maximum redshift for surviving globular clusters to form, both critical components for understanding the build-up of mass in the universe. Here we report measurements from the Southern Stellar Streams Spectroscopic Survey of the spatially thin, dynamically cold Phoenix stellar stream in the halo of the Milky Way. The properties of the Phoenix stream are consistent with it being the tidally disrupted remains of a globular cluster. However, its metal abundance ([Fe/H] = -2.7) is substantially below that of the empirical metallicity floor. The Phoenix stream thus represents the debris of the most metal-poor globular cluster discovered so far, and its progenitor is distinct from the present-day globular cluster population in the local Universe. Its existence implies that globular clusters below the metallicity floor have probably existed, but were destroyed during Galactic evolution.

preprint2019arXiv

A data-driven model of nucleosynthesis with chemical tagging in a lower-dimensional latent space

Chemical tagging seeks to identify unique star formation sites from present-day stellar abundances. Previous techniques have treated each abundance dimension as being statistically independent, despite theoretical expectations that many elements can be produced by more than one nucleosynthetic process. In this work we introduce a data-driven model of nucleosynthesis where a set of latent factors (e.g., nucleosynthetic yields) contribute to all stars with different scores, and clustering (e.g., chemical tagging) is modelled by a mixture of multivariate Gaussians in a lower-dimensional latent space. We use an exact method to simultaneously estimate the factor scores for each star, the partial assignment of each star to each cluster, and the latent factors common to all stars, even in the presence of missing data entries. We use an information-theoretic Bayesian principle to estimate the number of latent factors and clusters. Using the second Galah data release we find that six latent factors are preferred to explain N = 2,566 stars with 17 chemical abundances. We identify the rapid- and slow-neutron capture processes, as well as latent factors consistent with Fe-peak and α-element production, and another where K and Zn dominate. When we consider N ~ 160,000 stars with missing abundances we find another 7 factors, as well as 16 components in latent space. Despite these components showing separation in chemistry that is explained through different yield contributions, none show significant structure in their positions or motions. We argue that more data, and joint priors on cluster membership that are constrained by dynamical models, are necessary to realise chemical tagging at a galactic-scale. We release software that allows for model parameters to be optimised in seconds given a fixed number of latent factors, components, and $10^7$ abundance measurements.

preprint2019arXiv

Discovery of a 21 Myr old stellar population in the Orion complex

The Orion complex is arguably the most studied star-forming region in the Galaxy. While stars are still being born in the Orion nebula, the oldest part was believed to be no more than 13 Myr old. In order to study the full hierarchy of star formation across the Orion complex, we perform a clustering analysis of the Ori OB1a region using new stellar surveys and derive robust ages for each identified stellar aggregate. We use Gaia DR2 parameters supplemented with radial velocities from the GALAH and APOGEE surveys to perform clustering of the Ori OB1a association. Five overdensities are resolved in a six-dimensional parameter space (positions, distance, proper motions, and radial velocity). Most correspond to previously known structures (ASCC 16, 25 Orionis, ASCC 20, ASCC 21). We use Gaia DR2, Pan-STARRS1 and 2MASS photometry to fit isochrones to the colour-magnitude diagrams of the identified clusters. The ages of the clusters can thus be measured with ~10% precision. While four of the clusters have ages between 11 and 13 Myr, the ASCC 20 cluster stands out at an age of 21 $\pm$ 3 Myr. This is significantly greater than the age of any previously known component of the Orion complex. To some degree, all clusters overlap in at least one of the six phase-space dimensions. We argue that the formation history of the Orion complex, and its relation to the Gould belt, must be reconsidered. A significant challenge in reconstructing the history of the Ori OB1a association is to understand the impact of the newly discovered 21 Myr old population on the younger parts of the complex, including their formation.

preprint2019arXiv

The GALAH Survey: Chemodynamics of the Solar Neighbourhood

We present the chemodynamic structure of the solar neighbourhood using 62 814 stars within a 500 pc sphere of the Sun observed by GALAH and with astrometric parameters from Gaia DR2. We measure the velocity dispersion for all three components (vertical, radial, and tangential) and find that it varies smoothly with [Fe/H] and [$α$/Fe] for each component. The vertical component is especially clean, with $σ_{v_z}$ increasing from a low of $8$ km s$^{-1}$ at solar-[$α$/Fe] and [Fe/H] to a high of more than 50 km s$^{-1}$ for more metal-poor and [$α$/Fe] enhanced populations. We find no evidence of a large decrease in the velocity dispersion of the higher-[$α$/Fe] populations as claimed in analysis prior to Gaia DR2, although the trend of increasing velocity dispersion with [$α$/Fe] for the same metallicity does significantly flatten at high-[$α$/Fe]. The eccentricity distribution for local stars varies most strongly as a function of [$α$/Fe], where stars with [$α$/Fe]$<0.1$ dex having generally circular orbits ($e<0.15$), while the median eccentricity increases rapidly for more [$α$/Fe]-enhanced stellar populations up to $e\sim0.35$. These [$α$/Fe]-enhanced populations have guiding radii consistent with origins in the inner Galaxy. Of the stars with metallicities much higher than the local ISM ([Fe/H]>0.1 dex), we find that more than 70\% have $e<0.2$ and are likely observed in the solar neighbourhood through churning/migration rather than blurring effects, as the epicyclic motion for these stars is not large enough to reach the radii at which they were likely born based on their metallicity.

preprint2019arXiv

The GALAH Survey: Temporal Chemical Enrichment of the Galactic Disk

We present isochrone ages and initial bulk metallicities ($\rm [Fe/H]_{bulk}$, by accounting for diffusion) of 163,722 stars from the GALAH Data Release 2, mainly composed of main sequence turn-off stars and subgiants ($\rm 7000 K>T_{eff}>4000 K$ and $\rm log g>3$ dex). The local age-metallicity relationship (AMR) is nearly flat but with significant scatter at all ages; the scatter is even higher when considering the observed surface abundances. After correcting for selection effects, the AMR appear to have intrinsic structures indicative of two star formation events, which we speculate are connected to the thin and thick disks in the solar neighborhood. We also present abundance ratio trends for 16 elements as a function of age, across different $\rm [Fe/H]_{bulk}$ bins. In general, we find the trends in terms of [X/Fe] vs age from our far larger sample to be compatible with studies based on small ($\sim$ 100 stars) samples of solar twins but we now extend it to both sub- and super-solar metallicities. The $α$-elements show differing behaviour: the hydrostatic $α$-elements O and Mg show a steady decline with time for all metallicities while the explosive $α$-elements Si, Ca and Ti are nearly constant during the thin disk epoch (ages $\lessapprox $ 12 Gyr). The s-process elements Y and Ba show increasing [X/Fe] with time while the r-process element Eu have the opposite trend, thus favouring a primary production from sources with a short time-delay such as core-collapse supernovae over long-delay events such as neutron star mergers.