Source author record

Sanjib Sharma

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

49works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

49 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 "building block" 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'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 -- "C-19" ([Fe/H]$=-3.4$ dex), "Sylgr" ([Fe/H]$=-2.9$ dex) and "Phoenix" ([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'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', whereby no globular clusters are found with chemical (`metal') 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.

preprint2016arXiv

Stellar population synthesis based modelling of the Milky Way using asteroseismology of 13000 Kepler red giants

With current space-based missions it is now possible to obtain age-sensitive asteroseismic information for tens of thousands of red giants. This provides a promising opportunity to study the Galactic structure and evolution. We use asteroseismic data of red giants, observed by Kepler, to test the current theoretical framework of modelling the Galaxy based on population synthesis modeling and the use of asteroseismic scaling relations for giants. We use the open source code Galaxia to model the Milky Way and find the distribution of the masses predicted by Galaxia to be systematically offset with respect to the seismically-inferred observed masses. The Galactic model overestimates the number of low mass stars, and these stars are predominantly old and of low metallicity. Using corrections to the $Δν$ scaling relation suggested by stellar models (available for download) significantly reduces the disagreement between predicted and observed masses. For a few cases where non-seismic mass estimates are available, the corrections to $Δν$ also improve the agreement between seismic and non-seismic mass estimates. The disagreement between predictions of the Galactic model and the observations is most pronounced for stars with ${\rm [Fe/H]}<-0.5$ and ${\rm [Fe/H]}>0$ or for $T_{\rm eff}>4700$ K. Altering the star formation rate in order to suppress stars older than 10 Gyr improves the agreement for mass but leads to inconsistent color distributions. We also tested the predictions of the TRILEGAL Galactic model. However, unlike {\sl Galaxia}, it had difficulties in reproducing the photometric properties of the Kepler Input Catalog because it overestimates the number of blue stars. We conclude that either the scaling relations and/or the Galactic models need to be revised to reconcile predictions of theory with asteroseismic observations.

preprint2016arXiv

The GALAH Survey: Relative throughputs of the 2dF fibre positioner and the HERMES spectrograph from stellar targets

We present an analysis of the relative throughputs of the 3.9-metre Anglo-Australian Telescope's 2dF/HERMES system, based upon spectra acquired during the first two years of the GALAH survey. Averaged spectral fluxes of stars were compared to their photometry to determine the relative throughputs of fibres for a range of fibre position and atmospheric conditions. We find that overall the throughputs of the 771 usable fibres have been stable over the first two years of its operation. About 2.5 per cent of fibres have throughputs much lower than the average. There are also a number of yet unexplained variations between the HERMES bandpasses, and mechanically & optically linked fibre groups known as retractors or slitlets related to regions of the focal plane. These findings do not impact the science that HERMES will produce.

preprint2016arXiv

The GALAH survey: The data reduction pipeline

We present the data reduction procedures being used by the GALAH survey, carried out with the HERMES fibre-fed, multi-object spectrograph on the 3.9~m Anglo-Australian Telescope. GALAH is a unique survey, targeting 1 million stars brighter than magnitude V=14 at a resolution of 28,000 with a goal to measure the abundances of 29 elements. Such a large number of high resolution spectra necessitates the development of a reduction pipeline optimized for speed, accuracy, and consistency. We outline the design and structure of the Iraf-based reduction pipeline that we developed, specifically for GALAH, to produce fully calibrated spectra aimed for subsequent stellar atmospheric parameter estimation. The pipeline takes advantage of existing Iraf routines and other readily available software so as to be simple to maintain, testable and reliable. A radial velocity and stellar atmospheric parameter estimator code is also presented, which is used for further data analysis and yields a useful verification of the reduction quality. We have used this estimator to quantify the data quality of GALAH for fibre cross-talk level ($\lesssim0.5$%) and scattered light ($\sim5$ counts in a typical 20 minutes exposure), resolution across the field, sky spectrum properties, wavelength solution reliability (better than $1$ $\mathrm{km\ s^{-1}}$ accuracy) and radial velocity precision.

preprint2016arXiv

The K2 Ecliptic Plane Input Catalog (EPIC) and Stellar Classifications of 138,600 Targets in Campaigns 1-8

The K2 Mission uses the Kepler spacecraft to obtain high-precision photometry over ~80 day campaigns in the ecliptic plane. The Ecliptic Plane Input Catalog (EPIC) provides coordinates, photometry and kinematics based on a federation of all-sky catalogs to support target selection and target management for the K2 mission. We describe the construction of the EPIC, as well as modifications and shortcomings of the catalog. Kepler magnitudes (Kp) are shown to be accurate to ~0.1 mag for the Kepler field, and the EPIC is typically complete to Kp~17 (Kp~19 for campaigns covered by SDSS). We furthermore classify 138,600 targets in Campaigns 1-8 (~88% of the full target sample) using colors, proper motions, spectroscopy, parallaxes, and galactic population synthesis models, with typical uncertainties for G-type stars of ~3% in Teff, ~0.3 dex in log(g), ~40% in radius, ~10% in mass, and ~40% in distance. Our results show that stars targeted by K2 are dominated by K-M dwarfs (~41% of all selected targets), F-G dwarfs (~36%) and K giants (~21%), consistent with key K2 science programs to search for transiting exoplanets and galactic archeology studies using oscillating red giants. However, we find a significant variation of the fraction of cool dwarfs with galactic latitude, indicating a target selection bias due to interstellar reddening and the increased contamination by giant stars near the galactic plane. We discuss possible systematic errors in the derived stellar properties, and differences to published classifications for K2 exoplanet host stars. The EPIC is hosted at the Mikulski Archive for Space Telescopes (MAST): http://archive.stsci.edu/k2/epic/search.php.

preprint2015arXiv

First Light Results from the Hermes Spectrograph at the AAT

The High Efficiency and Resolution Multi Element Spectrograph, HERMES, is a facility-class optical spectrograph for the Anglo-Australian Telescope (AAT). It is designed primarily for Galactic Archaeology, the first major attempt to create a detailed understanding of galaxy formation and evolution by studying the history of our own galaxy, the Milky Way. The goal of the GALAH survey is to reconstruct the mass assembly history of the Milky Way through a detailed chemical abundance study of one million stars. The spectrograph is based at the AAT and is fed by the existing 2dF robotic fiber positioning system. The spectrograph uses volume phase holographic gratings to achieve a spectral resolving power of 28,000 in standard mode and also provides a high-resolution mode ranging between 40,000 and 50,000 using a slit mask. The GALAH survey requires an SNR greater than 100 for a star brightness of V ?= 14 in an exposure time of one hour. The total spectral coverage of the four channels is about 100 nm between 370 and 1000 nm for up to 392 simultaneous targets within the 2-degree field of view. HERMES has been commissioned over three runs, during bright time in October, November, and December 2013, in parallel with the beginning of the GALAH pilot survey, which started in November 2013. We present the first-light results from the commissioning run and the beginning of the GALAH survey, including performance results such as throughput and resolution, as well as instrument reliability.

preprint2014arXiv

New distances to RAVE stars

Probability density functions are determined from new stellar parameters for the distance moduli of stars for which the RAdial Velocity Experiment (RAVE) has obtained spectra with S/N>=10. Single-Gaussian fits to the pdf in distance modulus suffice for roughly half the stars, with most of the other half having satisfactory two-Gaussian representations. As expected, early-type stars rarely require more than one Gaussian. The expectation value of distance is larger than the distance implied by the expectation of distance modulus; the latter is itself larger than the distance implied by the expectation value of the parallax. Our parallaxes of Hipparcos stars agree well with the values measured by Hipparcos, so the expectation of parallax is the most reliable distance indicator. The latter are improved by taking extinction into account. The effective temperature absolute-magnitude diagram of our stars is significantly improved when these pdfs are used to make the diagram. We use the method of kinematic corrections devised by Schoenrich, Binney & Asplund to check for systematic errors for general stars and confirm that the most reliable distance indicator is the expectation of parallax. For cool dwarfs and low-gravity giants <pi> tends to be larger than the true distance by up to 30 percent. The most satisfactory distances are for dwarfs hotter than 5500 K. We compare our distances to stars in 13 open clusters with cluster distances from the literature and find excellent agreement for the dwarfs and indications that we are over-estimating distances to giants, especially in young clusters.

preprint2014arXiv

On the Shoulders of Giants: Properties of the Stellar Halo and the Milky Way Mass Distribution

Halo stars orbit within the potential of the Milky Way and hence their kinematics can be used to understand the underlying mass distribution. However, the inferred mass distribution depends sensitively upon assumptions made on the density and the velocity anisotropy profiles of the tracers. Also, there is a degeneracy between the parameters of the halo and that of the disk or bulge. Here, we decompose the Galaxy into bulge, disk and dark matter halo and then model the kinematic data of the halo BHB and K-giants from the SEGUE. Additionally, we use the gas terminal velocity curve and the Sgr A$^*$ proper motion. With $R_\odot = 8.5$kpc, our study reveals that the density of the stellar halo has a break at $17.2^{+1.1}_{-1.0}$ kpc, and an exponential cut-off in the outer parts starting at $97.7^{+15.6}_{-15.8}$kpc. Also, we find the velocity anisotropy is radially biased with $β_s= 0.4\pm{0.2}$ in the outer halo. We measure halo virial mass $M_{\text{vir}} = 0.80^{+0.31}_{-0.16} \times 10^{12} M_{\odot}$, concentration $c=21.1^{+14.8}_{-8.3}$, disk mass of $0.95^{+0.24}_{-0.30}\times10^{11} M_{\odot}$, disk scale length of $4.9^{+0.4}_{-0.4}$ kpc and bulge mass of $0.91^{+0.31}_{-0.38} \times 10^{10} M_{\odot}$. The mass of halo is found to be small and this has important consequences. The giant stars reveal that the outermost halo stars have low velocity dispersion interestingly suggesting a truncation of the stellar halo density rather than a small overall mass of the Galaxy. Our estimates of local escape velocity $v_{\rm esc} = 550.9^{+32.4}_{-22.1}$ kms$^{-1}$ and dark matter density $ρ^{\rm DM}_{\odot} = 0.0088^{+0.0024}_{-0.0018} M_{\odot} {\rm pc^{-3}} $ ($0.35^{+0.08}_{-0.07}$ GeV cm$^{-3}$) are in good agreement with recent estimates. Some of the above estimates are depended on the adopted value of $R_\odot$ and of outer power-law index of the tracer number density.

preprint2014arXiv

Pseudo Three-dimensional Maps of the Diffuse Interstellar Band at 862 nm

The diffuse interstellar bands (DIBs) are absorption lines observed in visual and near infrared spectra of stars. Understanding their origin in the interstellar medium (ISM) is one of the oldest problems in astronomical spectroscopy, as DIBs ahave been known since 1922. In a completely new approach to understanding DIBs, we combined information from nearly 500,000 stellar spectra obtained by the massive spectroscopic survey RAVE (Radial Velocity Experiment) to produce the first pseudo three-dimensional map of the strength of the DIB at 8620~Ångstroms covering the nearest 3 kiloparsecs from the Sun, and show that it follows our independently constructed spatial distribution of extinction by interstellar dust along the Galactic plane. Despite having a similar distribution in the Galactic plane, the DIB 8620 carrier has a significantly larger vertical scale height than the dust. Even if one DIB may not represent the general DIB population, our observations outline the future direction of DIB research.

preprint2014arXiv

Void asymmetries in the cosmic web: a mechanism for bulk flows

Bulk flows of galaxies moving with respect to the cosmic microwave background are well established observationally and seen in the most recent LCDM simulations. With the aid of an idealised Gadget-2 simulation, we show that void asymmetries in the cosmic web can exacerbate local bulk flows of galaxies. The Cosmicflows-2 survey, which has mapped in detail the 3D structure of the Local Universe, reveals that the Local Group resides in a "local sheet" of galaxies that borders a "local void" with a diameter of about 40 Mpc. The void is emptying out at a rate of 16 km/s/Mpc. In a co-moving frame, the Local Sheet is found to be moving away from the Local Void at ~ 260 km/s. Our model shows how asymmetric collapse due to unbalanced voids on either side of a developing sheet or wall can lead to a systematic movement of the sheet. We conjectured that asymmetries could lead to a large-scale separation of dark matter and baryons, thereby driving a dependence of galaxy properties with environment, but we do not find any evidence for this effect.

preprint2013arXiv

A Vertical Resonance Heating Model for X- or Peanut-Shaped Galactic Bulges

We explore a second order Hamiltonian vertical resonance model for X-shaped or peanut-shaped galactic bulges. The X-shape is caused by the 2:1 vertical Lindblad resonance with the bar, with two vertical oscillation periods per orbital period in the bar frame. We examine N-body simulations and find that due to the bar slowing down and disk thickening during bar buckling, the resonance and associated peanut-shape moves outward. The peanut-shape is consistent with the location of the vertical resonance, independent of whether the bar buckled or not. We estimate the resonance width from the potential m=4 Fourier component and find that the resonance is narrow, affecting orbits over a narrow range in the angular momentum distribution, dL/L ~ 0.05. As the resonance moves outward, stars originally in the mid plane are forced out of the mid plane into orbits just within the resonance separatrix. The height of the separatrix orbits, estimated from the Hamiltonian model, is approximately consistent with the peanut-shape height. The X-shape is comprised of stars in the vicinity of the resonance separatrix. The velocity distributions from the simulations illustrate that low inclination orbits are depleted within resonance. Within resonance, the vertical velocity distribution is broad, consistent with resonant heating caused by the passage of the resonance through the disk. In the Milky Way bulge we relate the azimuthally averaged mid-plane mass density near the vertical resonance to the rotation curve and bar pattern speed. At an estimated vertical resonance galactocentric radius of ~1.3 kpc, we confirm a mid-plane density of ~5x10^8 Msol/kpc^3, consistent with recently estimated mass distributions. We find that the rotation curve, bar pattern speed, 2:1 vertical resonance location, X-shape tips, and mid-plane mass density, are all self-consistent in the Milky Way galaxy bulge.

preprint2013arXiv

An empirical formula for the distribution function of a thin exponential disc

An empirical formula for a Shu distribution function that reproduces a thin disc with exponential surface density to good accuracy is presented. The formula has two free parameters that specify the functional form of the velocity dispersion. Conventionally, this requires the use of an iterative algorithm to produce the correct solution, which is computationally taxing for applications like Markov Chain Monte Carlo (MCMC) model fitting. The formula has been shown to work for flat, rising and falling rotation curves. Application of this methodology to one of the Dehnen distribution functions is also shown. Finally, an extension of this formula to reproduce velocity dispersion profiles that are an exponential function of radius is also presented. Our empirical formula should greatly aid the efficient comparison of disc models with large stellar surveys or N-body simulations.

preprint2013arXiv

Exploring the Variable Sky with LINEAR. II. Halo Structure and Substructure Traced by RR Lyrae Stars to 30 kpc

We present a sample of ~5,000 RR Lyrae stars selected from the recalibrated LINEAR dataset and detected at heliocentric distances between 5 kpc and 30 kpc over ~8,000 deg^2 of sky. The coordinates and light curve properties, such as period and Oosterhoff type, are made publicly available. We find evidence for the Oosterhoff dichotomy among field RR Lyrae stars, with the ratio of the type II and I subsamples of about 1:4. The number density distribution of halo RRab stars as a function of galactocentric distance can be described as an oblate ellipsoid with the axis ratio q=0.63 and with either a single or a double power law with a power-law index in the range -2 to -3. Using a group-finding algorithm EnLink, we detected seven candidate halo groups, only one of which is statistically spurious. Three of these groups are near globular clusters (M53/NGC 5053, M3, M13), and one is near a known halo substructure (Virgo Stellar Stream); the remaining three groups do not seem to be near any known halo substructures or globular clusters, and seem to have a higher ratio of Oosterhoff type II to Oosterhoff type I RRab stars than what is found in the halo. The extended morphology and the position (outside the tidal radius) of some of the groups near globular clusters is suggestive of tidal streams possibly originating from globular clusters. Spectroscopic followup of detected halo groups is encouraged.

preprint2013arXiv

The effect of feedback and reionization on star formation in low-mass dwarf galaxy haloes

We simulate the evolution of a 10^9 Msun dark matter halo in a cosmological setting with an adaptive-mesh refinement code as an analogue to local low luminosity dwarf irregular and dwarf spheroidal galaxies. The primary goal of our study is to investigate the roles of reionization and supernova feedback in determining the star formation histories of low mass dwarf galaxies. We include a wide range of physical effects, including metal cooling, molecular hydrogen formation and cooling, photoionization and photodissociation from a metagalactic background, a simple prescription for self-shielding, star formation, and a simple model for supernova driven energetic feedback. We carry out simulations excluding each major effect in turn. We find that reionization is primarily responsible for expelling most of the gas in our simulations, but that supernova feedback is required to disperse the dense, cold gas in the core of the halo. Moreover, we show that the timing of reionization can produce an order of magnitude difference in the final stellar mass of the system. For our full physics run with reionization at z=9, we find a stellar mass of about 10^5 Msun at z=0, and a mass-to-light ratio within the half-light radius of approximately 130 Msun/Lsun, consistent with observed low-luminosity dwarfs. However, the resulting median stellar metallicity is 0.06 Zsun, considerably larger than observed systems. In addition, we find star formation is truncated between redshifts 4 and 7, at odds with the observed late time star formation in isolated dwarf systems but in agreement with Milky Way ultrafaint dwarf spheroidals. We investigate the efficacy of energetic feedback in our simple thermal-energy driven feedback scheme, and suggest that it may still suffer from excessive radiative losses, despite reaching stellar particle masses of about 100 Msun, and a comoving spatial resolution of 11 pc.

preprint2012arXiv

Kinematics of the stellar halo and the mass distribution of the Milky Way using BHB stars

Here we present a kinematic study of the Galactic halo out to a radius of $\sim$ 60 kpc, using 4664 blue horizontal branch (BHB) stars selected from the SDSS/SEGUE survey, to determine key dynamical properties. Using a maximum likelihood analysis, we determine the velocity dispersion profiles in spherical coordinates ($σ_{r}$, $σ_θ$, $σ_ϕ$) and the anisotropy profile ($β$). The radial velocity dispersion profile ($σ_{r}$) is measured out to a galactocentric radius of $r \sim 60$ kpc, but due to the lack of proper-motion information, $σ_θ$, $σ_ϕ$ and $β$ could only be derived directly out to $r \sim25$ kpc. From a starting value of $β\approx 0.5$ in the inner parts ($9<r/\kpc<12$), the profile falls sharply in the range $r \approx 13-18$ kpc, with a minimum value of $β=-1.2$ at $r=17$ kpc, rising sharply at larger radius. In the outer parts, in the range $25<r/\kpc<56$, we predict the profile to be roughly constant with a value of $β\approx 0.5$. The newly discovered kinematic anomalies are shown not to arise from halo substructures. We also studied the anisotropy profile of simulated stellar halos formed purely by accretion and found that they cannot reproduce the sharp dip seen in the data. From the Jeans equation, we compute the stellar rotation curve ($v_{\rm circ}$) of the Galaxy out to $r \sim 25$ kpc. The mass of the Galaxy within $r \lesssim 25$ kpc is determined to be $2.1 \times 10^{11}$ $M_{\sun}$, and with a 3-component fit to $v_{\rm circ}(r)$, we determine the virial mass of the Milky Way dark matter halo to be $M_{\rm vir} = 0.9 ^{+0.4}_{-0.3} \times 10^{12}$ $M_{\sun}$ ($R_{\rm vir} = 249^{+34}_{-31}$ kpc).

preprint2012arXiv

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

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

preprint2012arXiv

Probing the Halo From the Solar Vicinity to the Outer Galaxy: Connecting Stars in Local Velocity Structures to Large-Scale Clouds

(Abridged) This paper presents the first connections made between two local features in velocity-space found in a survey of M giant stars and stellar spatial inhomogeneities on global scales. Comparison to cosmological, chemodynamical stellar halo models confirm that the M giant population is particularly sensitive to rare, recent and massive accretion events. These events can give rise to local observed velocity sequences - a signature of a small fraction of debris from a common progenitor, passing at high velocity through the survey volume, near the pericenters of their eccentric orbits. The majority of the debris is found in much larger structures, whose morphologies are more cloud-like than stream-like and which lie at the orbital apocenters. Adopting this interpretation, the full-space motions represented by the observed velocity features are derived under the assumption that the members within each sequence share a common velocity. Orbit integrations are then used to trace the past and future trajectories of these stars across the sky revealing plausible associations with large, previously-discovered, cloud-like structures. The connections made between nearby velocity structures and these distant clouds represent preliminary steps towards developing coherent maps of such giant debris systems. These maps promise to provide new insights into the origin of debris clouds, new probes of Galactic history and structure, and new constraints on the high-velocity tails of the local dark matter distribution that are essential for interpreting direct detection experiments.

preprint2012arXiv

The Lives of High Redshift Mergers

We present a comparative study of recent works on merger-timescales with dynamical friction and find a strong contrast between idealized/isolated mergers (Boylan-Kolchin et al. 2008) and mergers from a cosmological volume (Jiang et al. 2008). Our study measures the duration of mergers in a cosmological N-body simulation of dark matter, with emphasis on higher redshifts (z < 10) and a lower mass range. In our analysis we consider and compare two merger definitions; tidal disruption and coalescence. We find that the merger-time formula proposed by Jiang et al. (2008) describes our results well and conclude that cosmologically motivated merger-time formulae provide a more versatile and statistically robust approximation for practical applications such as semi-analytic/hybrid models.

preprint2011arXiv

Galaxia: a code to generate a synthetic survey of the Milky Way

We present here a fast code for creating a synthetic survey of the Milky Way. Given one or more color-magnitude bounds, a survey size and geometry, the code returns a catalog of stars in accordance with a given model of the Milky Way. The model can be specified by a set of density distributions or as an N-body realization. We provide fast and efficient algorithms for sampling both types of models. As compared to earlier sampling schemes which generate stars at specified locations along a line of sight, our scheme can generate a continuous and smooth distribution of stars over any given volume. The code is quite general and flexible and can accept input in the form of a star formation rate, age metallicity relation, age velocity dispersion relation and analytic density distribution functions. Theoretical isochrones are then used to generate a catalog of stars and support is available for a wide range of photometric bands. As a concrete example we implement the Besancon Milky Way model for the disc. For the stellar halo we employ the simulated stellar halo N-body models of Bullock & Johnston (2005). In order to sample N-body models, we present a scheme that disperses the stars spawned by an N-body particle, in such a way that the phase space density of the spawned stars is consistent with that of the N-body particles. The code is ideally suited to generating synthetic data sets that mimic near future wide area surveys such as GAIA, LSST and HERMES. As an application we study the prospect of identifying structures in the stellar halo with a simulated GAIA survey. We plan to make the code publicly available at http://galaxia.sourceforge.net.

preprint2010arXiv

Chemical signatures of the first star clusters

The chemical abundance patterns of the oldest stars in the Galaxy are expected to contain residual signatures of the first stars in the early universe. Numerous studies attempt to explain the intrinsic abundance scatter observed in some metal-poor populations in terms of chemical inhomogeneities dispersed throughout the early Galactic medium due to discrete enrichment events. Just how the complex data and models are to be interpreted with respect to "progenitor yields" remains an open question. Here we show that stochastic chemical evolution models to date have overlooked a crucial fact. Essentially all stars today are born in highly homogeneous star clusters and it is likely that this was also true at early times. When this ingredient is included, the overall scatter in the abundance plane [Fe/H] vs. [X/Fe] (C-space), where X is a nucleosynthetic element, can be much less than derived from earlier models. Moreover, for moderately flat cluster mass functions (gamma < 2), and/or for mass functions with a high mass cut-off (M_max > 10^5 M_sun), stars exhibit a high degree of clumping in C-space that can be identified even in relatively small data samples. Since stellar abundances can be modified by mass transfer in close binaries, clustered signatures are essential for deriving the yields of the first supernovae. We present a statistical test to determine whether a given set of observations exhibit such behaviour. Our initial work focusses on two dimensions in C-space, but we show that the clustering signal can be greatly enhanced by additional abundance axes. The proposed experiment will be challenging on existing 8-10m telescopes, but relatively straightforward for a multi-object echelle spectrograph mounted on a 25-40m telescope.

preprint2010arXiv

Group finding in the stellar halo using M-giants in 2MASS: An extended view of the Pisces Overdensity?

A density based hierarchical group-finding algorithm is used to identify stellar halo structures in a catalog of M-giants from the Two Micron All Sky Survey (2MASS). The intrinsic brightness of M-giant stars means that this catalog probes deep into the halo where substructures are expected to be abundant and easy to detect. Our analysis reveals 16 structures at high Galactic latitude (greater than 15 degree), of which 10 have been previously identified. Among the six new structures two could plausibly be due to masks applied to the data, one is associated with a strong extinction region and one is probably a part of the Monoceros ring. Another one originates at low latitudes, suggesting some contamination from disk stars, but also shows protrusions extending to high latitudes, implying that it could be a real feature in the stellar halo. The last remaining structure is free from the defects discussed above and hence is very likely a satellite remnant. Although the extinction in the direction of the structure is very low, the structure does match a low temperature feature in the dust maps. While this casts some doubt on its origin, the low temperature feature could plausibly be due to real dust in the structure itself. The angular position and distance of this structure encompass the Pisces overdensity traced by RR Lyraes in Stripe 82 of the Sloan Digital Sky Survey (SDSS). However, the 2MASS M-giants indicate that the structure is much more extended than what is visible with the SDSS, with the point of peak density lying just outside Stripe 82. The morphology of the structure is more like a cloud than a stream and reminiscent of that seen in simulations of satellites disrupting along highly eccentric orbits.

preprint2010arXiv

Group Finding in the Stellar Halo using Photometric Surveys: Current Sensitivity and Future Prospects

The ubiquity of substructure in the stellar halo has already been demonstrated by the SDSS and 2MASS and future surveys promise to explore the halo in ever more detail. This paper examines what can be learnt from current and future photometric-databases using group-finding techniques. We compare groups recovered from a sample of M-giants from 2MASS with those found in synthetic surveys of simulated $Λ$CDM stellar halos and demonstrate broad consistency. We also find that these recovered groups are likely to represent the majority of high-luminosity ($\log(L/L_{Sun})>6.67$) satellites accreted within the last 10 Gyr and on orbits with apocenters within 100 kpc. However the sensitivity of the M-giant survey to accretion events that were either ancient, from low-luminosity objects or those on radial orbits is limited because of the low number of stars, bias towards high-metallicity stars and the shallow depth. We examine the extent to which these limitations are addressed by current and future surveys, in particular catalogues of main-sequence turn-off (MSTO) stars from SDSS and LSST, and of RR-Lyrae stars from LSST or PanSTARRS. The MSTO surveys are more sensitive to low-luminosity events ($\log(L/L_{Sun}~5$ or less) than the 2MASS M-giant sample, while the RR-Lyrae surveys, with superior depth, are good at detecting events on highly eccentric orbits. When combined we expect these photometric surveys to provide a comprehensive picture of the last 10 Gyr of Galactic accretion. Events older than this are too phase mixed to be discovered and would require additional information such as velocity and chemical abundance.

preprint2010arXiv

The Dawning of the Stream of Aquarius in RAVE

We identify a new, nearby (0.5 < d < 10 kpc) stream in data from the RAdial Velocity Experiment (RAVE). As the majority of stars in the stream lie in the constellation of Aquarius we name it the Aquarius Stream. We identify 15 members of the stream lying between 30 < l < 75 and -70< b <-50, with heliocentric line-of-sight velocities V_los~-200 km/s. The members are outliers in the radial velocity distribution, and the overdensity is statistically significant when compared to mock samples created with both the Besançon Galaxy model and newly-developed code Galaxia. The metallicity distribution function and isochrone fit in the log g - T_eff plane suggest the stream consists of a 10 Gyr old population with [m/H]~-1.0. We explore relations to other streams and substructures, finding the stream cannot be identified with known structures: it is a new, nearby substructure in the Galaxy's halo. Using a simple dynamical model of a dissolving satellite galaxy we account for the localization of the stream. We find that the stream is dynamically young and therefore likely the debris of a recently disrupted dwarf galaxy or globular cluster. The Aquarius stream is thus a specimen of ongoing hierarchical Galaxy formation, rare for being right in the solar suburb.

preprint2008arXiv

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

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