Search

We explore the fundamental relations governing the radial and vertical velocity dispersions of stars in the Milky Way, from combined studies of complementary surveys including GALAH, LAMOST, APOGEE, the NASA Kepler and K2 missions, and Gaia DR2. We find that different stellar samples, even though they target different tracer populations and employ a variety of age estimation techniques, follow the same set of fundamental relations. We provide the clearest evidence to date that, in addition to the well-known dependence on stellar age, the velocity dispersions of stars depend on orbital angular momentum Lz, metallicity, and height above the plane |z|, and are well described by a multiplicatively separable functional form. The dispersions have a power-law dependence on age with exponents of 0.441 ± 0.007 and 0.251 ± 0.006 for σz and σR, respectively, and the power law is valid even for the oldest stars. For the solar neighbourhood stars, the apparent break in the power law for older stars, as seen in previous studies, is due to the anticorrelation of Lz with age. The dispersions decrease with increasing Lz until we reach the Sun's orbital angular momentum, after which σz increases (implying flaring in the outer disc) while σR flattens. For a given age, the dispersions increase with decreasing metallicity, suggesting that the dispersions increase with birth radius. The dispersions also increase linearly with |z|. The same set of relations that work in the solar neighbourhood also work for stars between 3 < R/kpc < 20. Finally, the high-[α/Fe] stars follow the same relations as the low-[α/Fe] stars.

Using kinematics from Gaia and the large elemental abundance space of the second data release of the GALAH survey, we identify two new members of the Fimbulthul stellar stream, and chemically tag them to massive,multimetallic globular cluster ω Centauri. Recent analysis of the second data release of Gaia had revealed the Fimbulthul stellar stream in the halo of the Milky Way. It had been proposed that the stream is associated with the ω Cen, but this proposition relied exclusively upon the kinematics and metallicities of the stars to make the association. In thiswork,we find our two newmembers of the stream to be metal-poor stars that are enhanced in sodium and aluminium, typical of second population globular cluster stars, but not otherwise seen in field stars. Furthermore, the stars share the s-process abundance pattern seen in ω Cen, which is rare in field stars. Apart from one star within 1.5 deg of ω Cen, we find no other stars observed by GALAH spatially near ω Cen or the Fimbulthul stream that could be kinematically and chemically linked to the cluster. Chemically tagging stars in the Fimbulthul stream to ω Cen confirms the earlier work, and further links this tidal feature in the Milky Way halo to ω Cen.

We present the serendipitous discovery of the fastest main-sequence hyper-velocity star (HVS) by the Southern Stellar Stream Spectroscopic Survey (S5). The star S5-HVS1 is a ∼2.35M⊙A-type star located at a distance of ∼9 kpc from the Sun and has a heliocentric radial velocity of 1017 ± 2.7 kms-1 without any signature of velocity variability. The current 3D velocity of the star in the Galactic frame is 1755 ± 50 kms-1. When integrated backwards in time, the orbit of the star points unambiguously to the Galactic Centre, implying that S5-HVS1 was kicked away from Sgr A+ with a velocity of ∼1800 kms-1, travelled for 4.8 Myr to its current location. This is so far the only HVS confidently associated with the Galactic Centre. S5-HVS1 is also the first hyper-velocity star to provide constraints on the geometry and kinematics of the Galaxy, such as the Solar motion Vy⊙ = 246.1 ± 5.3 km s-1 or position R0 = 8.12 ± 0.23 kpc. The ejection trajectory and transit time of S5-HVS1 coincide with the orbital plane and age of the annular disc of young stars at the Galactic Centre, and thus may be linked to its formation. With the S5-HVS1 ejection velocity being almost twice the velocity of other hyper-velocity stars previously associated with the Galactic Centre, we question whether they have been generated by the same mechanism or whether the ejection velocity distribution has been constant over time.

Globular clusters are some of the oldest bound stellar structures observed in the Universe1. They are ubiquitous in large galaxies and are believed to trace intense star-formation events and the hierarchical build-up of structure2,3. 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 Sun4–6. 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 Universe7. 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 the empirical metallicity floor. The Phoenix stream thus represents the debris of the most metal-poor globular clusters 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.

Asteroseismology is a promising tool to study Galactic structure and evolution because it can probe the ages of stars. Earlier attempts comparing seismic data from the Kepler satellite with predictions from Galaxy models found that the models predicted more low-mass stars compared to the observed distribution of masses. It was unclear if the mismatch was due to inaccuracies in the Galactic models, or the unknown aspects of the selection function of the stars. Using new data from the K2 mission, which has a well-defined selection function, we find that an old metal-poor thick disc, as used in previous Galactic models, is incompatible with the asteroseismic information. We use an importance-sampling framework, which takes the selection function into account, to fit for the metallicities of a population synthesis model using spectroscopic data. We show that spectroscopic measurements of [Fe/H] and [α/Fe] elemental abundances from the GALAH survey indicate a mean metallicity of log (Z/Z☉) = −0.16 for the thick disc. Here Z is the effective solar-scaled metallicity, which is a function of [Fe/H] and [α/Fe]. With the revised disc metallicities, for the first time, the theoretically predicted distribution of seismic masses show excellent agreement with the observed distribution of masses. This indirectly verifies that the asteroseismic mass scaling relation is good to within five per cent. Assuming the asteroseismic scaling relations are correct, we estimate the mean age of the thick disc to be about 10 Gyr, in agreement with the traditional idea of an old α-enhanced thick disc.

We use the second data releases of the European Space AgencyGaia astrometric survey and the high-resolution Galactic Archaeology with HERMES (GALAH) spectroscopic survey to analyse the structure of our Galaxy's disc components. With GALAH, we separate the α-rich and α-poor discs (with respect to Fe), which are superposed in both position and velocity space, and examine their distributions in action space. We study the distribution of stars in the zVz phase plane, for both Vφ and VR, and recover the remarkable 'phase spiral' discovered by Gaia. We identify the anticipated quadrupole signature in zVz of a tilted velocity ellipsoid for stars above and below the Galactic plane. By connecting our work with earlier studies, we show that the phase spiral is likely to extend well beyond the narrow solar neighbourhood cylinder in which it was found. The phase spiral is a signature of corrugated waves that propagate through the disc, and the associated non-equilibrium phase mixing. The radially asymmetric distribution of stars involved in the phase spiral reveals that the corrugation, which is mostly confined to the α-poor disc, grows in z-amplitude with increasing radius. We present new simulations of tidal disturbance of the Galactic disc by the Sagittarius (Sgr) dwarf. The effect on the zVz phase plane lasts ≥2Gyr, but a subsequent disc crossing wipes out the coherent structure. We find that the phase spiral was excited ≤0.5 Gyr ago by an object like Sgr with total mass ∼3 × 1010 M⊙ (stripped down from ∼5 × 1010 M⊙ when it first entered the halo) passing through the plane.

The overlap between the spectroscopic Galactic Archaeology with HERMES (GALAH) survey and Gaia provides a high-dimensional chemodynamical space of unprecedented size. We present a first analysis of a subset of this overlap, of 7066 dwarf, turn-off, and sub-giant stars. These stars have spectra from the GALAH survey and high parallax precision from the Gaia DR1 Tycho-Gaia Astrometric Solution. We investigate correlations between chemical compositions, ages, and kinematics for this sample. Stellar parameters and elemental abundances are derived from the GALAH spectra with the spectral synthesis code SPECTROSCOPY MADE EASY. We determine kinematics and dynamics, including action angles, from the Gaia astrometry and GALAH radial velocities. Stellar masses and ages are determined with Bayesian isochrone matching, using our derived stellar parameters and absolute magnitudes. We report measurements of Li, C, O, Na, Mg, Al, Si, K, Ca, Sc, Ti, V, Cr, Mn, Co, Ni, Cu, Zn, Y, as well as Ba and we note that we have employed non-LTE calculations for Li, O, Al, and Fe. We show that the use of astrometric and photometric data improves the accuracy of the derived spectroscopic parameters, especially log g. Focusing our investigation on the correlations between stellar age, iron abundance [Fe/H], and mean alpha-enhancement [α/Fe] of the magnitude-selected sample, we recover the result that stars of the high-α sequence are typically older than stars in the low-α sequence, the latter spanning iron abundances of -0.7 < [Fe/H] < +0.5. While these two sequences become indistinguishable in [α/Fe] vs. [Fe/H] at the metal-rich regime, we find that age can be used to separate stars from the extended high-α and the low-α sequence even in this regime. When dissecting the sample by stellar age, we find that the old stars (>8 Gyr) have lower angular momenta Lz than the Sun, which implies that they are on eccentric orbits and originate from the inner disc. Contrary to some previous smaller scale st

Gaia DR2 has revealed new small-scale and large-scale patterns in the phase-space distribution of stars in the Milky Way. In cylindrical Galactic coordinates (R, φ, z), ridge-like structures can be seen in the (R, Vφ) plane and asymmetric arch-like structures in the (VR, Vφ) plane. We show that the ridges are also clearly present when the third dimension of the (R, Vφ) plane is represented by z, Vz, VR, [Fe/H], and [α/Fe]. The maps suggest that stars along the ridges lie preferentially close to the Galactic mid-plane (|z| < 0.2 kpc), and have metallicity and α elemental abundance similar to that of the Sun. We show that phase mixing of disrupting spiral arms can generate both the ridges and the arches. It also generates discrete groupings in orbital energy − the ridges and arches are simply surfaces of constant energy. We identify eight distinct ridges in the Gaia DR2 data: six of them have constant energy while two have constant angular momentum. Given that the signature is strongest for stars close to the plane, the presence of ridges in z and Vz suggests a coupling between planar and vertical directions. We demonstrate, using N-body simulations that such coupling can be generated both in isolated discs and in discs perturbed by an orbiting satellite like the Sagittarius dwarf galaxy.

We introduce the southern stellar stream spectroscopy survey (S5), an on-going program to map the kinematics and chemistry of stellar streams in the southern hemisphere. The initial focus of S5 has been spectroscopic observations of recently identified streams within the footprint of the dark energy survey (DES), with the eventual goal of surveying streams across the entire southern sky. Stellar streams are composed of material that has been tidally striped from dwarf galaxies and globular clusters and hence are excellent dynamical probes of the gravitational potential of the Milky Way, as well as providing a detailed snapshot of its accretion history. Observing with the 3.9 m Anglo-Australian Telescope's 2-degree-Field fibre positioner and AAOmega spectrograph, and combining the precise photometry of DES DR1 with the superb proper motions from Gaia DR2, allows us to conduct an efficient spectroscopic survey to map these stellar streams. So far S5 has mapped nine DES streams and three streams outside of DES; the former are the first spectroscopic observations of these recently discovered streams. In addition to the stream survey, we use spare fibres to undertake a Milky Way halo survey and a low-redshift galaxy survey. This paper presents an overview of the S5 program, describing the scientific motivation for the survey, target selection, observation strategy, data reduction, and survey validation. Finally, we describe early science results on stellar streams and Milky Way halo stars drawn from the survey. Updates on S5, including future public data releases, can be found at http://s5collab.github.io.

We present a study using the second data release of the GALAH survey of stellar parameters and elemental abundances of 15 pairs of stars identified by Oh et al. They identified these pairs as potentially co-moving pairs using proper motions and parallaxes from Gaia DR1. We find that 11 very wide (>1 pc) pairs of stars do in fact have similar Galactic orbits, while a further four claimed co-moving pairs are not truly co-orbiting. Eight of the 11 co-orbiting pairs have reliable stellar parameters and abundances, and we find that three of those are quite similar in their abundance patterns, while five have significant [Fe/H] differences. For the latter, this indicates that they could be co-orbiting because of the general dynamical coldness of the thin disc, or perhaps resonances induced by the Galaxy, rather than a shared formation site. Stars such as these, wide binaries, debris of past star formation episodes, and coincidental co-orbiters, are crucial for exploring the limits of chemical tagging in the Milky Way.

Context. 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. Aims. 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. Methods. 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. Results. While four of the clusters have ages between 11 and 13 Myr, the ASCC 20 cluster stands out at an age of 21 ± 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. Conclusions. 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.

The latest Gaia data release enables us to accurately identify stars that are more luminous than would be expected on the basis of their spectral type and distance. During an investigation of the 329 best solar twin candidates uncovered among the spectra acquired by the GALAH survey, we identified 64 such overluminous stars. In order to investigate their exact composition, we developed a data-driven methodology that can generate a synthetic photometric signature and spectrum of a single star. By combining multiple such synthetic stars into an unresolved binary or triple system and comparing the results to the actual photometric and spectroscopic observations, we uncovered 6 definitive triple stellar system candidates and an additional 14 potential candidates whose combined spectrum mimics the solar spectrum. Considering the volume correction factor for a magnitude-limited survey, the fraction of probable unresolved triple stars with long orbital periods is ~2 per cent. Possible orbital configurations of the candidates were investigated using the selection and observational limits. To validate the discovered multiplicity fraction, the same procedure was used to evaluate the multiplicity fraction of other stellar types.

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 first results from the K2-HERMES project, which uses the HERMES multi-object spectrograph on the Anglo-Australian Telescope to obtain R ∼ 28000 spectra of up to 360 stars in one exposure. This ongoing project aims to derive self-consistent spectroscopic parameters for about half of K2 target stars. We present complete stellar parameters and isochrone-derived masses and radii for 46 stars hosting 57 K2 candidate planets in Campaigns 1-3. Our revised host-star radii cast severe doubt on three candidate planets: EPIC 201407812.01, EPIC 203070421.01, and EPIC 202843107.01, all of which now have inferred radii well in excess of the largest known inflated Jovian planets.

The Galactic Archaeology with HERMES (GALAH) survey is a large-scale stellar spectroscopic survey of theMilkyWay, designed to deliver complementary chemical information to a large number of stars covered by the Gaia mission. We present the GALAH second public data release (GALAH DR2) containing 342 682 stars. For these stars, the GALAH collaboration provides stellar parameters and abundances for up to 23 elements to the community. Here we present the target selection, observation, data reduction, and detailed explanation of how the spectra were analysed to estimate stellar parameters and element abundances. For the stellar analysis, we have used a multistep approach. We use the physics-driven spectrum synthesis of Spectroscopy Made Easy (SME) to derive stellar labels (Teff, logg, [Fe/H], [X/Fe], vmic, vsin i, AKS) for a representative training set of stars. This information is then propagated to the whole sample with the data-driven method of The Cannon. Special care has been exercised in the spectral synthesis to only consider spectral lines that have reliable atomic input data and are little affected by blending lines. Departures from local thermodynamic equilibrium (LTE) are considered for several key elements, including Li, O, Na, Mg, Al, Si, and Fe, using 1D MARCS stellar atmosphere models. Validation tests including repeat observations, Gaia benchmark stars, open and globular clusters, and K2 asteroseismic targets lend confidence to our methods and results. Combining the GALAH DR2 catalogue with the kinematic information from Gaia will enable a wide range of Galactic Archaeology studies, with unprecedented detail, dimensionality, and scope.

The technique of chemical tagging uses the elemental abundances of stellar atmospheres to 'reconstruct' chemically homogeneous star clusters that have long since dispersed. The GALAH spectroscopic survey - which aims to observe one million stars using the Anglo- Australian Telescope - allows us to measure up to 30 elements or dimensions in the stellar chemical abundance space, many of which are not independent. How to find clustering reliably in a noisy high-dimensional space is a difficult problem that remains largely unsolved. Here, we explore t-distributed stochastic neighbour embedding (t-SNE) - which identifies an optimal mapping of a high-dimensional space into fewer dimensions - whilst conserving the original clustering information. Typically, the projection is made to a 2D space to aid recognition of clusters by eye. We show that this method is a reliable tool for chemical tagging because it can: (i) resolve clustering in chemical space alone, (ii) recover known open and globular clusters with high efficiency and low contamination, and (iii) relate field stars to known clusters. t-SNE also provides a useful visualization of a high-dimensional space. We demonstrate the method on a data set of 13 abundances measured in the spectra of 187 000 stars by the GALAH survey. We recover seven of the nine observed clusters (six globular and three open clusters) in chemical space with minimal contamination from field stars and low numbers of outliers. With chemical tagging, we also identify two Pleiades supercluster members (which we confirm kinematically), one as far as 6° - one tidal radius away from the cluster centre.

Open cluster members are coeval and share the same initial bulk chemical composition. Consequently, differences in surface abundances between members of a cluster that are at different evolutionary stages can be used to study the effects of mixing and internal chemical processing. We carry out an abundance analysis of seven elements (Li, O, Na, Mg, Al, Si, and Fe) in 66 stars belonging to the open cluster M67, based on high resolution GALAH spectra, 1D MARCS model atmospheres, and non-local thermodynamic equilibrium (non-LTE) radiative transfer. From the non-LTE analysis, we find a typical star-to-star scatter in the abundance ratios of around 0.05 dex. We find trends in the abundance ratios with effective temperature, indicating systematic differences in the surface abundances between turn-off and giant stars; these trends are more pronounced when LTE is assumed. However, trends with effective temperature remain significant for Al and Si also in non-LTE. Finally, we compare the derived abundances with prediction from stellar evolution models including effects of atomic diffusion. We find overall good agreement for the abundance patterns of dwarfs and sub-giant stars, but the abundances of cool giants are lower relative to less evolved stars than predicted by the diffusion models, in particular for Mg.

Using data from the GALAH pilot survey, we determine properties of the Galactic thin and thick discs near the solar neighbourhood. The data cover a small range of Galactocentric radius (7.9 ≲ RGC ≲ 9.5 kpc), but extend up to 4 kpc in height from the Galactic plane, and several kpc in the direction of Galactic anti-rotation (at longitude 260° ≤ ℓ ≤ 280°). This allows us to reliably measure the vertical density and abundance profiles of the chemically and kinematically defined 'thick' and 'thin' discs of the Galaxy. The thin disc (low-α population) exhibits a steep negative vertical metallicity gradient, at d[M/H]/dz = -0.18 ± 0.01 dex kpc-1, which is broadly consistent with previous studies. In contrast, its vertical α-abundance profile is almost flat, with a gradient of d[α/M]/dz = 0.008 ± 0.002 dex kpc-1. The steep vertical metallicity gradient of the low-α population is in agreement with models where radial migration has a major role in the evolution of the thin disc. The thick disc (high-α population) has a weaker vertical metallicity gradient d[M/H]/dz = -0.058 ± 0.003 dex kpc-1. The α-abundance of the thick disc is nearly constant with height, d[α/M]/dz = 0.007 ± 0.002 dex kpc-1. The negative gradient in metallicity and the small gradient in [α/M] indicate that the high-α population experienced a settling phase, but also formed prior to the onset of major Type Ia supernova enrichment. We explore the implications of the distinct α-enrichments and narrow [α/M] range of the sub-populations in the context of thick disc formation.

We present the data reduction procedures being used by the GALactic Archeology with Hermes (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 necessitate 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 (≲0.5 per cent) and scattered light (~5 counts in a typical 20 min exposure), resolution across the field, sky spectrum properties, wavelength solution reliability (better than 1 kms-1 accuracy), and radial velocity precision.

Galactic globular clusters (GCs) are now known to harbour multiple stellar populations, which are chemically distinct in many light element abundances. It is becoming increasingly clear that asymptotic giant branch (AGB) stars in GCs show different abundance distributions in light elements compared to those in the red giant branch (RGB) and other phases, skewing towards more primordial, field-star-like abundances, which we refer to as subpopulation one (SP1). As part of a larger programme targeting giants in GCs, we obtained high-resolution spectra for a sample of 106 RGB and 15 AGB stars in Messier 4 (NGC 6121) using the 2dF+HERMES facility on the Anglo-Australian Telescope. In this Letter, we report an extreme paucity of AGB stars with [Na/O] >-0.17 in M4, which contrasts with the RGB that has abundances up to [Na/O] =0.55. The AGB abundance distribution is consistent with all AGB stars being from SP1. This result appears to imply that all subpopulation two stars (SP2; Na-rich, O-poor) avoid the AGB phase. This is an unexpected result given M4's horizontal branch morphology - it does not have an extended blue horizontal branch. This is the first abundance study to be performed utilizing the HERMES spectrograph.

We present an analysis of the relative throughputs of the 3.9-m Anglo-Australian Telescope's 2dF/HERMES (High Efficiency and Resolution Multi-Element Spectrograph) system, based upon spectra acquired during the first two years of the Galactic Archaeology with HERMES 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 and 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.

The Galactic Archaeology with HERMES (GALAH) survey is a large high-resolution spectroscopic survey using the newly commissioned High Efficiency and Resolution Multi-Element Spectrograph (HERMES) on the Anglo-Australian Telescope. The HERMES spectrograph provides high-resolution (R ~ 28 000) spectra in four passbands for 392 stars simultaneously over a 2 deg field of view. The goal of the survey is to unravel the formation and evolutionary history of the Milky Way, using fossil remnants of ancient star formation events which have been disrupted and are now dispersed throughout the Galaxy. Chemical tagging seeks to identify such dispersed remnants solely from their common and unique chemical signatures; these groups are unidentifiable from their spatial, photometric or kinematic properties. To carry out chemical tagging, the GALAH survey will acquire spectra for a million stars down to V ~ 14. The HERMES spectra of FGK stars contain absorption lines from 29 elements including light proton-capture elements, α-elements, odd-Z elements, iron-peak elements and n-capture elements from the light and heavy s-process and the r-process. This paper describes the motivation and planned execution of the GALAH survey, and presents some results on the first-light performance of HERMES.