Search

We present a novel method to constrain the formation time of the Milky Way disc using the chrono-kinematic signatures of the inner Galaxy. We construct an O-rich Mira variable sample from the Gaia long-period variable catalogue to study the kinematic behaviour of stars with different ages in the inner Galaxy. From the Auriga suite of cosmological zoom-in simulations, we find that the age of the oldest stellar population with imprints of the bar in density and kinematics matches the disc spin-up epoch. This is because stars born before the spin-up show insufficient rotation and are not kinematically cold enough to be efficiently trapped by the bar. We find that the bar kinematic signature disappears for Mira variables with a period shorter than 190 d. Using the period–age relation of Mira variables, we constrain the spin-up epoch of the Milky Way to be younger than |$\sim 11{-}12$|  Gyr (redshift |$\sim 3$|⁠). We also discuss and compare our method and result to other evidence of the Milky Way spin-up epoch under the context of a realistic age uncertainty. Age uncertainty leads to an overestimation of the disc formation time when performing backward modelling. Our constrain of the spin-up epoch is independent from previous studies because it relies on the kinematics of the inner Galaxy instead of the solar vicinity. [ABSTRACT FROM AUTHOR]

Stars move away from their birthplaces over time via a process known as radial migration, which blurs chemo–kinematic relations used for reconstructing the Milky Way (MW) formation history. To understand the true time evolution of the MW, one needs to take into account the effects of this process. We show that stellar birth radii can be derived directly from the data with minimum prior assumptions on the Galactic enrichment history. This is done by first recovering the time evolution of the stellar birth metallicity gradient, |$\mathrm{ d}\mathrm{[Fe/H]}(R, \tau)/\mathrm{ d}R$|⁠ , through its inverse relation to the metallicity range as a function of age today, allowing us to place any star with age and metallicity measurements back to its birthplace, R |$_b$|⁠. Applying our method to a large high-precision data set of MW disc subgiant stars, we find a steepening of the birth metallicity gradient from 11 to 8 Gyr ago, which coincides with the time of the last massive merger, Gaia–Sausage–Enceladus (GSE). This transition appears to play a major role in shaping both the age–metallicity relation and the bimodality in the [ |$\alpha$| /Fe]–[Fe/H] plane. By dissecting the disc into mono-R |$_b$| populations, clumps in the low-[ |$\alpha$| /Fe] sequence appear, which are not seen in the total sample and coincide in time with known star-formation bursts, possibly associated with the Sagittarius Dwarf Galaxy. We estimated that the Sun was born at |$4.5\pm 0.4$| kpc from the Galactic centre. Our R |$_b$| estimates provide the missing piece needed to recover the Milky Way formation history. [ABSTRACT FROM AUTHOR]

We follow up emission line galaxies identified through the near-infrared slitless HST /WFC3 WISP survey with VLT/FORS2 optical spectroscopy. Over 4 WISP fields, we targeted 85 of 138 line emission objects at |$0.4\lt z\lt 2$| identified in WFC3 spectroscopy. Half the galaxies are fainter than |$H_{AB}=24$|  mag, and would not have been included in many well-known surveys based on broad-band magnitude selection. We confirm 95 per cent of the initial WFC3 grism redshifts in the 38 cases where we detect lines in FORS2 spectroscopy. However, for targets which exhibited a single emission line in WFC3, up to 65 per cent at |$z\lt 1.28$| did not have expected emission lines detected in FORS2 and hence may be spurious (although this false-detection rate improves to 33 per cent using the latest public WISP emission line catalogue). From the Balmer decrement, the extinction of the WISP galaxies is consistent with |$A($| H  |$\alpha)=1$|  mag. From SED fits to multiband photometry including Spitzer |$3.6\, \mu$| m, we find a median stellar mass of |$\log _{10}(M_\star /{\rm M}_{\odot })=8.94$|⁠. Our emission-line-selected galaxies tend to lie above the star-forming main sequence (i.e. higher specific star formation rates). Using [O  iii ], [O  ii ], and H β lines to derive gas-phase metallicities, we find typically sub-solar metallicities, decreasing with redshift. Our WISP galaxies lie below the |$z=0$| mass–metallicity relation, and galaxies with higher star formation rates tend to have lower metallicity. Finally, we find a strong increase with redshift of the H α rest-frame equivalent width in this emission-line selected sample, with higher |$EW_0$| galaxies having larger [O  iii ]/H β and O32 ratios on average, suggesting lower metallicity or higher ionization parameter in these extreme emission line galaxies. [ABSTRACT FROM AUTHOR]

Under the assumption of hierarchical galaxy formation, dwarf galaxies are the closest existing analogues to the high-redshift protogalaxies that merged to form the Milky Way. These low-mass systems serve as unique laboratories for studying nucleosynthetic channels given that the chemical compositions of their stars play a pivotal role in constraining their chemical enrichment history. To date, stellar abundances in dwarf galaxies have focused almost exclusively on elemental abundance ratios. While important, elemental abundances omit critical information about the isotopic composition. Here, we compute the Mg isotopic ratios of six accreted dwarf galaxy stars (low |$\alpha$|⁠) and seven Milky Way stars (high |$\alpha$|⁠) using a set of high-resolution (65 000 < R < 160 000) and high signal-to-noise ratio (⁠|$\rm {S/N} \gt 250$|⁠) optical spectra. We show, for the first time, that at a given [Fe/H] stars born in a dwarf galaxy differ in their Mg isotopic ratios from stars born in the Milky Way. However, when comparing isotopic ratios at a given [Mg/H] rather than [Fe/H], a powerful diagnostic emerges that suggests nucleosynthesis processes are consistent across different stellar environments. This universality of Mg isotopic abundances provides additional dimensionality for chemical evolution models and helps to constrain massive star nucleosynthesis across cosmic time. [ABSTRACT FROM AUTHOR]

There is no doubt that the Milky Way is a barred galaxy; however, factors that establish its prominent morphology remain largely elusive. In this work, we attempt to constrain the history of the MW by tracing the present-day parameters and evolution of a set of MW and M31 analogues from the TNG50 simulations. We find that the strength of bars at |$z=0$| correlates well not only with the stellar mass build-up but also, more crucially, with the time of onset of stellar discs. Discs of strongly barred galaxies form early (⁠|$z \gtrsim 2-3$|⁠), compared to weakly and non-barred galaxies (⁠|$z \approx 1-1.5$|⁠). Although we are cautious to draw ultimate conclusions about the governing factor of discs formation due to the complexity and correlations between different phenomena, the observed morphological diversity can be tentatively explained by a substantial variation in the gas angular momentum around proto-galaxies already at |$z\approx 3-5$|⁠ ; in such a way, early discs formed from gas with larger angular momentum. By comparing the formation time-scales of discs of barred galaxies in the TNG50 sample, we infer that the MW has a strong bar (⁠|$0.35\lt A_2\lt 0.6$|⁠), and that its stellar disc started to dominate over the spheroidal component already at |$z \approx 2$|⁠ , with a mass of |$\approx 1 \pm 0.5 \times 10^{10} M_\odot$|⁠. We conclude that the presence of a strong bar in the MW is a natural manifestation of the early formation of the stellar disc, which made possible bursty but highly efficient star formation at high redshift. [ABSTRACT FROM AUTHOR]

We present a deep X-Shooter rest-frame ultraviolet (UV) to optical spectral analysis of the heavily reddened quasar, ULASJ2315+0143 at |$z=2.566$|⁠ , known to reside in a major-merger host galaxy. The rest-frame optical is best fit by a dust-reddened quasar (E (B – V) |$^{\rm QSO}= 1.55$|⁠) with black-hole mass |$\rm log_{10}(H\beta , M{\small BH} [{\rm M}_{\odot }]) = 10.26 \pm 0.05$|⁠ , bolometric luminosity |$\rm L_{Bol}$|  = |$\rm 10^{48.16}\, erg\,\,s^{-1}$| and Eddington-scaled accretion rate log |$_{10}(\rm \lambda _{Edd}) = -0.19$|⁠. We find remarkable similarities between ULASJ2315+0143 and the high-redshift little red dots (LRDs). The rest-frame UV cannot be explained by a dusty quasar component alone and requires an additional blue component consistent with either a star-forming host galaxy or scattered AGN light. We detect broad high-ionization emission lines in the rest-UV, supporting the scattered light interpretation for the UV excess. The scattering fraction represents just 0.05  per cent of the total luminosity of ULASJ2315+0143 . Analysis of the mid-infrared SED suggests an absence of hot dust on torus-scales similar to what is observed for LRDs. The obscuring medium is therefore likely on galaxy scales. We detect narrow, blueshifted associated absorption-line systems in C  iv , N  v , Si  iv , and Si  iii. There is evidence for significant high-velocity (>1000  |$\rm km\, s^{-1}$|⁠) outflows in both the broad- and narrow-line regions as traced by C  iv and [O  iii ] emission. The kinetic power of the [O  iii ] wind is |$\dot{\epsilon }_{k}^{ion} = 10^{44.61} \rm erg\, s^{-1} \sim 0.001\, L_{Bol}$|⁠. ULASJ2315+0143 is likely in an important transition phase where star formation, black-hole accretion and multiphase gas flows are simultaneously occurring. [ABSTRACT FROM AUTHOR]

Spatial, kinematic, and orbital properties, along with ages and chemical compositions of the thin disc, thick disc, and various stellar substructures in the halo, are studied based on data from the LAMOST and Gaia surveys. The star formation in the Galactic thin and thick disc, with peak metallicities of − 0.20 and − 0.45 dex, is found to have peaked about 5.5 and 12.5 Gyr ago, respectively. The thin disc is also found to have experienced an initial star formation burst about 12.5 Gyr ago. The pro-grade population Splash and hot-disc (HD), with peak metallicity of about − 0.60 and − 0.43, are found to be about 13.03 and 12.21 Gyr old, respectively, with peak eccentricity of 0.70 and 0.35, are understood to be of in situ origin. The Gaia-Enceladus/Sausage (GE/S), Thamnos and Sequoia, with peak metallicity of about − 1.31, − 1.36, and − 1.56, are found to be about 11.66, 12.89, and 12.18 Gyr, respectively, and are understood to be remnants of dwarf galaxies merged with the Milky Way. The HD, Splash, and Thamnos have experienced chemical evolution similar to the thick disc, while GE/S, Sequoia, and Helmi stream have experienced distinct chemical enrichment of iron and α -process elements. [ABSTRACT FROM AUTHOR]

The North–South asymmetry in the number density and bulk velocity of stars in the Solar neighbourhood provides valuable insights into the formation and evolution of the Milky Way disc. Our objective is to investigate the wave-like disc oscillations of mono-age stellar populations in the Solar neighbourhood using data from Gaia Data Release 3. We have selected a comprehensive sample of main-sequence turn-off stars. The ages of these stars can be accurately determined using isochrone fitting methods. Our findings indicate that the North–South density and mean vertical velocity asymmetries remain consistent across all age groups. The uniformity of perturbations across all subsamples suggests that all populations are responding to the same external influence, which likely affects them irrespective of their age. Moreover, the fact that these perturbations appear consistently implies they could be either ongoing or recent. Regarding vertical velocity dispersions, we observe that older stars exhibit larger dispersions. [ABSTRACT FROM AUTHOR]

We present novel insights into the interplay between tidal forces and star formation in interacting galaxies before their first pericentre passage. We investigate seven close pair galaxies devoid of visible tidal disturbances, such as tails, bridges, and shells. Using integral field spectroscopy data of extended Calar Alto Legacy Integral Field Area, we unveil a previously unreported phenomenon: H  |$\alpha$| emission, a proxy for recent star formation, exhibits a significant enhancement in regions facing the companion galaxy, reaching up to 1.9 times higher flux compared to opposite directions. Notably, fainter companions within pairs display a more pronounced one-sided H  |$\alpha$| excess, exceeding the typical range observed in isolated galaxies with 2 |$\sigma$| confidence level. Furthermore, the observed H  |$\alpha$| excess in fainter companion galaxies exhibits a heightened prominence at the outer galactic regions. These findings suggest that tidal forces generated before the first pericentre passage exert a stronger influence on fainter galaxies due to their shallower potential wells by their brighter companions. This unveils a more intricate interplay between gravitational interactions and star formation history within interacting galaxies than previously understood, highlighting the need further to explore the early stages of interaction in galaxy evolution. [ABSTRACT FROM AUTHOR]

Hot, |$\log (T/\mathrm{K})$| |$\sim$| 7.5, gas was recently discovered in the Milky Way in extragalactic sightlines. In order to determine its location, here we present sightlines to Galactic X-ray binaries (XRBs) passing through the Interstellar Medium (ISM). In this pilot study, we investigate absorption features of S |$_{\mathrm{XVI}}$|⁠ , Si |$_{\mathrm{XIV}}$|⁠ , and Ne |$_{\mathrm{X}}$| in the spectra of three XRBs, namely 4U 1735–44, 4U 1820–30, and Cyg X-2, using Chandra High Energy Transmission Grating archival observations. We do not detect any of these lines. We determine the 2 |$\sigma$| upper limit for the equivalent widths of the undetected absorption lines and the column densities of the corresponding ions. We note that the 2 |$\sigma$| upper limits for S |$_{\mathrm{XVI}}$| K |$\alpha$| and Si |$_{\mathrm{XIV}}$| K |$\alpha$| are an order of magnitude smaller than those previously detected in the extragalactic sightlines. Our finding suggests that if any gas at |$\log (T/\mathrm{K})$| |$\gt\ 7$| is present in the Galactic ISM, it is unlikely to be ubiquitous. This is an important result because it implies that S |$_{\mathrm{XVI}}$|⁠ , Si |$_{\mathrm{XIV}}$|⁠ , and Ne |$_{\mathrm{X}}$| absorption detected in extragalactic sightlines is not from the ISM, but is likely from a hot gas phase in the extraplanar region beyond the ISM or in the extended Circumgalactic Medium. [ABSTRACT FROM AUTHOR]

The question of whether the Milky Way's disc extends to low metallicity has been the subject of debate for many years. We aim to address the question by employing a large sample of giant stars with radial velocities and homogeneous metallicities based on the Gaia Data Release 3 BP/RP(XP) spectra. We study the 3D velocity distribution of stars in various metallicity ranges, including the very metal-poor (VMP) regime (⁠|$\mathrm{[M/H]} < -2.0$|⁠). We find that a clear, stand-alone disc population, i.e. that with a ratio of rotational velocity to velocity dispersion |$v/\sigma > 1$|⁠ , starts to emerge only around |$\mathrm{[M/H]} \sim -1.3$|⁠ , and is not visible for |$\mathrm{[M/H]} < -1.6$|⁠. Using Gaussian mixture modelling, we show that there are two halo populations in the VMP regime: one stationary and one with a net prograde rotation of |$\sim\!\! 80\ \mathrm{km\, s}^{-1}$|⁠. In this low-metallicity range, we are able to place constraints on the contribution of a rotation-supported thick disc sub-population to a maximum of |$\sim\!\! 3$| per cent in our sample. We compare our results to previous claims of discy VMP stars in both observations and simulations and find that having a prograde halo component could explain most of these. [ABSTRACT FROM AUTHOR]

Making use of the APOGEE DR17 catalogue with high quality data for 143 509 red giant branch stars we explore the strength of different mechanisms that causes a star to radially migrate in the Milky Way stellar disc. At any position in the disc we find stars that are more metal-rich than the local interstellar medium. This is surprising and normally attributed to the migration of these stars after their formation inside their current Galactocentric radius. Such stars are prime candidates for studying the strength of different migratory processes. We specifically select two types of metal-rich stars: (i) super metal-rich stars (⁠|$\mathrm{[Fe/H]}\gt 0.2$|⁠) and (ii) stars that are more metal-rich than their local environment. For both, we explore the distribution of orbital parameters and ages as evidence of their migration history. We find that most super metal-rich stars have experienced some amount of churning as they have orbits with |$R_g\gtrsim 5 \,{\rm kpc}$|⁠. Furthermore, about half of the super metal-rich stars are on non-circular orbits (⁠|$\mathrm{ecc} \gt 0.15$|⁠) and therefore also have experienced blurring. The metallicity of young stars in our sample is generally the same as the metallicity of the interstellar medium, suggesting they have not radially migrated yet. Stars with lower metallicity than the local environment have intermediate to old ages. We further find that super metal-rich stars have approximately the same age distribution at all Galactocentric radii, which suggests that radial migration is a key mechanism responsible for the chemical compositions of stellar populations in the Milky Way. [ABSTRACT FROM AUTHOR]

The study presents a theoretical framework for understanding the role of dark matter on the stability of the galactic disc. We model the galaxy as a two-component system consisting of stars and gas in equilibrium with an external dark matter halo. We derive the equations governing the growth of perturbations and obtain a stability criterion that connects the potential of the dark matter halo and the gas fraction with the stability levels of the galaxy. We find that a two-component disc is more susceptible to the growth of gravitational instabilities than individual components, particularly as gas fractions increase. However, the external field, due to the dark matter halo, acts as a stabilizing agent and increases the net stability levels even in the presence of a cold gas component. We apply the stability criterion to models of the Milky Way, low surface brightness galaxies, and baryon-dominated cold rotating disc galaxies observed in the early universe. Our results show that the potential due to the dark matter halo plays a significant role in stabilizing nearby galaxies, such as the Milky Way, and low surface brightness galaxies, which would otherwise be prone to local gravitational instabilities. However, we find that the baryon-dominated cold disc galaxies observed in the early universe remain susceptible to the growth of local gravitational instabilities despite the stabilizing effect of the dark matter halo. [ABSTRACT FROM AUTHOR]

The role of type Ia supernovae (SN Ia), mainly the delay time distributions (DTDs) determined by the binary systems, and the yields of elements created by different explosion mechanisms, are studied by using the MulChem chemical evolution model applied to our Galaxy. We explored 15 DTDs and 12 tables of elemental yields produced by different SN Ia explosion mechanisms, doing a total of 180 models. Chemical abundances for |$\alpha$| -elements (O, Mg, Si, and Ca) and Fe derived from these models are compared with recent solar region observational data of |$\alpha$| -elements over Fe relative abundances, [X/Fe], as a function of [Fe/H] and age. A multidimensional maximum-likelihood analysis shows that 52 models are able to fit all these data sets simultaneously, considering the 1 |$\sigma$| confidence level. The combination of STROLG1 DTD from Strolger et al. (2020) and LN20181 SN Ia yields from Leung & Nomoto (2018) provides the best fit. The exponential model with very prompt events is a possible DTD, but a combination of several channels is more probable. The SN Ia yields that include MCh or Near MCh correspond to 39 (75 per cent) of the 52 best models. Regarding the DTD, 31 (60 per cent) of the 52 most probable models correspond to the SD scenario, while the remaining 21 (40 per cent) are based on the DD scenario. Our results also show that the relatively large dispersion of the observational data may be explained by the stellar migration from other radial regions, and/or perhaps a combination of DTDs and explosion channels. [ABSTRACT FROM AUTHOR]

We investigate the non-adiabatic effect of time-dependent deformations in the Milky Way (MW) halo potential on stellar streams. Specifically, we consider the MW's response to the infall of the Large Magellanic Cloud (LMC) and how this impacts our ability to recover the spherically averaged MW mass profile from observation using stream actions. Previously, action clustering methods have only been applied to static or adiabatic MW systems to constrain the properties of the host system. We use a time-evolving MW–LMC simulation described by basis function expansions. We find that for streams with realistic observational uncertainties on shorter orbital periods and without close encounters with the LMC, e.g. GD-1, the radial action distribution is sufficiently clustered to locally recover the spherical MW mass profile across the stream radial range within a |$2\sigma$| confidence interval determined using a Fisher information approach. For streams with longer orbital periods and close encounters with the LMC, e.g. Orphan–Chenab (OC), the radial action distribution disperses as the MW halo has deformed non-adiabatically. Hence, for OC streams generated in potentials that include an MW halo with any deformations, action clustering methods will fail to recover the spherical mass profile within a |$2\sigma$| uncertainty. Finally, we investigate whether the clustering of stream energies can provide similar constraints. Surprisingly, we find for OC-like streams, the recovered spherically averaged mass profiles demonstrate less sensitivity to the time-dependent deformations in the potential. [ABSTRACT FROM AUTHOR]

Exploring the marks left by galactic accretion in the Milky Way helps us understand how our Galaxy was formed. However, finding and studying accreted stars and the galaxies they came from has been challenging. This study uses a simulation from the Numerical Investigation of a Hundred Astronomical Objects project, which now includes a wider range of chemical compositions, to find better ways to spot these accreted stars. By comparing our findings with data from the GALAH spectroscopic survey, we confirm that the observationally established diagnostics of [Al/Fe] versus [Mg/Mn] also show a separation of in situ and accreted stars in the simulation, but stars from different accretion events tend to overlap in this plane even without observational uncertainties. Looking at the relationship between stellar age and linear or logarithmic abundances, such as [Fe/H], we can clearly separate different groups of these stars if the uncertainties in their chemical makeup are less than 0.15 dex and less than 20 per cent for their ages. This method shows promise for studying the history of the Milky Way and other galaxies. Our work highlights how important it is to have accurate measurements of stellar ages and chemical content. It also shows how simulations can help us understand the complex process of galaxies merging and suggest how these events might relate to the differences we see between our Galaxy's thin and thick disc stars. This study provides a way to compare theoretical models with real observations, opening new paths for research in both our own Galaxy and beyond. [ABSTRACT FROM AUTHOR]

In recent years, the study of the Milky Way has significantly advanced due to extensive spectroscopic surveys of its stars, complemented by astroseismic and astrometric data. However, it remains disjoint from recent advancements in understanding the physics of the Galactic interstellar medium (ISM). This paper introduces a new model for the chemical evolution of the Milky Way that can be constrained on stellar data, because it combines a state-of-the-art ISM model with a Milky Way stellar disc model. Utilizing a data set of red clump stars from APOGEE, known for their precise ages and metallicities, we concentrate on the last 6 billion years – a period marked by Milky Way's secular evolution. We examine the oxygen abundance in the low- |$\alpha$| disc stars relative to their ages and birth radii, validating or constraining critical ISM parameters that remain largely unexplored in extragalactic observations. The models that successfully reproduce the radius–metallicity distribution and the age–metallicity distribution of stars without violating existing ISM observations indicate a need for modest differential oxygen enrichment in Galactic outflows, meaning that the oxygen abundance of outflows is higher than the local ISM abundance, irrespective of outflow mass loading. The models also suggest somewhat elevated ISM gas velocity dispersion levels over the past 6 billion years compared to galaxies of similar mass. The extra turbulence necessary could result from energy from gas accretion onto the Galaxy, supernovae clustering in the ISM, or increased star formation efficiency per freefall time. This work provides a novel approach to constraining the Galactic ISM and outflows, leveraging the detailed insights available from contemporary Milky Way surveys. [ABSTRACT FROM AUTHOR]

We need to resolve the individual stars for binary-fraction determinations of stellar systems. Therefore, it is not possible to obtain the binary fractions for dense or distant stellar systems. We propose a method to determine the binary fraction of a dense or distant stellar system: this method is first to determine the binary-fraction variation for any two adjacent regions and then to add up those binary-fraction variations along the radial direction to obtain the binary fraction for a stellar system. The binary-fraction variation is derived using 10 binary-fraction-sensitive spectral absorption feature indices (SAFIs) and the binary-fraction variation calibrations in terms of these SAFIs. Using this method, we first present the binary-fraction variations for 21 Galactic globular clusters (GCs). By comparison, we find that they agree well with the binary fractions based on the main-sequence fiducial line method from previous studies. This verifies that the above-mentioned method is feasible. Next, we present the binary-fraction variations of 13 Galactic GCs. We give the relationships between the binary fraction and various parameters, and find that the binary fraction is negatively correlated with N HB and N RR, the binary fraction of some studies is not strongly correlated with N BS, and the number of GCs with large binary fraction is greater for the extreme blue horizontal branch population ratio. Finally, if we want to obtain a more accurate binary fraction, we suggest that spectroscopic and photometric observations are conducted at an appropriate area interval for a stellar system. [ABSTRACT FROM AUTHOR]

We present a sample of 22 massive galaxies with stellar masses >1010 M⊙ at 3 < z < 4 with deep H and K -band high resolution spectra (R  = 3500–3000) from Keck/MOSFIRE and VLT/KMOS near-infrared spectrographs. We find a large fraction have strong [O  iii ] 5007 and Hβ emission lines with large line widths (σ 100–450 km s−1). We measure the sizes of our galaxies from Hubble Space Telescope images and consider the potential kinematic scaling relations of our sample, and rule out an explanation for these broad lines in terms of galaxy-wide kinematics. Based on consideration of the [O  iii ] 5007/Hβ flux ratios, their location in the mass–excitation diagram, and the derived bolometric luminosities, we conclude that active galactic nuclei (AGN) and their narrow-line regions most likely give rise to this emission. At redshifts 3 < z < 4, we find significantly high AGN fractions in massive galaxies, ranging from 60–70 per cent for the mass range 10 < log (M ⋆/M⊙) < 11, with a lower limit 30 per cent for all galaxies within that redshift range when we apply our most stringent AGN criteria. We also find a considerably lower AGN fraction in massive quiescent galaxies, ranging from 20–30 per cent. These fractions of AGN point to the period between 3 < z < 4 being a time of heightened activity for the development of supermassive black holes in the massive end of the galaxy population and provide evidence for their role in the emergence of the first massive quenched galaxies at this epoch. [ABSTRACT FROM AUTHOR]

We derive astrophysical parameters of the open cluster NGC 1513 by means of colour indices built with new CCD  UBV (RI) KC photometry. Based on early-type members, the mean foreground reddening and total to selective extinction ratio are E (B − V) = 0.79 ± 0.09 mag and RV = 2.85 ± 0.05. Through the differential grid method, we derive the metal abundance [Fe/H] = −0.06 dex (Z = +0.013), which is consistent with the value [Fe/H] = −0.088 of the bright giant member – LAMOST 695710060. Z = +0.013 isochrone fit to the V × (B − V) colour–magnitude diagram leads to a turn-off age of 224 ± 27 Myr (thus an intermediate-age cluster), and a distance modulus of (V 0 − MV) = 10.90 ± 0.15 mag, thus implying a distance from the Sun d = 1514 ± 105 pc. Within the uncertainties, our photometric distance is consistent with the value d = 1435 ± 14 pc from the Gaia DR3 parallax. We find signs of small mass segregation through a minimum spanning tree analysis for the 190 most massive stars, together with the rather steep mass function (χ = +2.39) slope. The high core to half-light radius ratio R core/ R h = 0.82, together with the compact half-light to tidal radius ratio R h/ R t = 0.22, suggest that it is probably related to cluster-formation effects, due to little dynamical evolution, instead of driving its dynamical evolution by internal relaxation. Indeed, NGC 1513 is located in the second quadrant (ℓ = 152 |${_{.}^{\circ}}$| 59 and Galactocentric distance R GC = 9.57 kpc), which tends to minimize tidal effects by external processes and tidal disruption. Therefore, internal mass segregation effects in NGC 1513 seem to be more efficient than cluster evaporation processes. We find that NGC 1513 migrated about 0.50 kpc from its birth place. [ABSTRACT FROM AUTHOR]

We present a set of controlled hydrodynamical simulations to study the effects of strong galactic outflows on the density and temperature structures, and associated X-ray signatures, of extra-planar and circumgalactic gas. We consider three initial state models, isothermal, isentropic, and rotating cooling-flow, for the hot circumgalactic medium (CGM) into which the outflows are driven. The energy sources are either stellar winds and supernovae, or active galactic nuclei. We consider energy injection rates in the range |$10^{40} \lt \dot{E}_{\rm inj} \lt 10^{44.5}$| erg s−1, and compute the time-dependent soft X-ray (0.5–2 keV) surface brightness. For |$\dot{E}_{\rm inj} \gtrsim 10^{41} {\small --} 10^{42}$| erg s−1, with the exact threshold depending on the initial CGM state, the X-ray response is dominated by dense hot gas in the forward shock that eventually fades into the CGM as a sound wave. The shock surrounds an inner hot bubble leading to a radial flattening of the X-ray surface brightness. For lower energy injection rates, the X-ray surface brightness of the initial CGM state is almost unaffected. We present analytical approximations for the outflow shock propagation and the associated X-ray emissions. [ABSTRACT FROM AUTHOR]

The Milky Way's last significant merger, the Gaia Enceladus/Sausage (GES), is thought to have taken place between |$8\hbox{ and }11\, \mathrm{Gyr}$| ago. Recent studies in the literature suggest that the bar of the Milky Way is rather old, indicating that it formed at a similar epoch to the GES merger. We investigate the possible link between these events using one of the Auriga cosmological simulations that has salient features in common with the Milky Way, including a last significant merger with kinematic signatures resembling that of the GES. In this simulation, the GES-like merger event triggers tidal forces on the disc, gas inflows, and a burst of star formation, with the formation of a bar occurring within |$1\, \mathrm{Gyr}$| of the first pericentre. To highlight the effects of the merger, we rerun the simulation from z  = 4 with the progenitors of the GES-like galaxy removed well before the merger time. The consequence is a delay in bar formation by around |$2\, \mathrm{Gyr}$|⁠ , and this new bar forms without any significant external perturbers. We conclude that this Milky Way-like simulation shows a route to the real Milky Way's bar forming around the epoch of the GES merger due to tidal forces on its first pericentre. We explore all Auriga galaxies with GES-like merger events, and find that those with stellar mass ratios below |$10{{\, \rm per\, cent}}$| form bars within |$1\, \mathrm{Gyr}$| of the merger, while bar formation is delayed in the more massive merger scenarios. These include the four oldest bars in the simulation suite. Lastly, we note some later morphological differences between the disc of the original simulation and our rerun, in particular that the latter does not grow radially for the final |$7\, \mathrm{Gyr}$|⁠. Our study suggests that the GES may therefore be responsible for the formation of the Milky Way's bar, as well as for the build-up of its extended disc. [ABSTRACT FROM AUTHOR]

A key event in the history of the Milky Way is the formation of the bar. This event affects the subsequent structural and dynamical evolution of the entire Galaxy. When the bar formed, gas was likely rapidly funnelled to the centre of the Galaxy settling in a star-forming nuclear disc. The Milky Way bar formation can then be dated by considering the age distribution of the oldest stars in the formed nuclear stellar disc. In this highly obscured and crowded region, reliable age tracers are limited, but bright, high-amplitude Mira variables make useful age indicators as they follow a period–age relation. We fit dynamical models to the proper motions of a sample of Mira variables in the Milky Way's nuclear stellar disc region. Weak evidence for inside-out growth and both radial and vertical dynamical heating with time of the nuclear stellar disc is presented, suggesting that the nuclear stellar disc is dynamically well-mixed. Furthermore, for Mira variables around a ∼350-d period, there is a clear transition from nuclear stellar disc-dominated kinematics to background bar-bulge-dominated kinematics. Using a Mira variable period–age relation calibrated in the solar neighbourhood, this suggests the nuclear stellar disc formed in a significant burst in star formation |$(8\pm 1)\, \mathrm{Gyr}$| ago, although the data are also weakly consistent with a more gradual formation of the nuclear stellar disc at even earlier epochs. This implies a relatively early formation time for the Milky Way bar (⁠|$\gtrsim 8\, \mathrm{Gyr}$|⁠), which has implications for the growth and state of the young Milky Way and its subsequent history. [ABSTRACT FROM AUTHOR]

From the chemodynamical properties of tidal debris in the Milky Way, it has been inferred that the dwarf satellites that have been disrupted had different chemical abundances from their present-day counterparts of similar mass that survive today, specifically, they had lower [Fe/H] and higher [Mg/Fe]. Here we use the A rtemis simulations to study the relation between the chemical abundances of disrupted progenitors of MW-mass galaxies and their stellar mass, and the evolution of the stellar mass–metallicity relations (MZR) of this population with redshift. We find that these relations have significant scatter, which correlates with the accretion redshifts (z acc) of satellites, and with their cold gas fractions. We investigate the MZRs of dwarf populations accreted at different redshifts and find that they have similar slopes, and also similar with the slope of the MZR of the surviving population (≈0.32). However, the entire population of disrupted dwarfs displays a steeper MZR, with a slope of ≈0.48, which can be explained by the changes in the mass spectrum of accreted dwarf galaxies with redshift. We find strong relations between the (mass-weighted) 〈 z acc〉 of the disrupted populations and their global chemical abundances (〈[Fe/H]〉 and 〈[Mg/Fe]〉), which suggests that chemical diagnostics of disrupted dwarfs can be used to infer the types of merger histories of their hosts. For the case of the MW, our simulations predict that the bulk of the disrupted population was accreted at 〈 z acc〉 ≈ 2, in agreement with other findings. We also find that disrupted satellites form and evolve in denser environments, closer to their hosts, than their present-day counterparts. [ABSTRACT FROM AUTHOR]

In this study, we modify the semi-analytic model galacticus in order to accurately reproduce the observed properties of dwarf galaxies in the Milky Way. We find that reproducing observational determinations of the halo occupation fraction and mass–metallicity relation for dwarf galaxies requires us to include H2 cooling, an updated ultraviolet background radiation model, and to introduce a model for the metal content of the intergalactic medium. By fine-tuning various model parameters and incorporating empirical constraints, we have tailored the model to match the statistical properties of Milky Way dwarf galaxies, such as their luminosity function and size–mass relation. We have validated our modified semi-analytic framework by undertaking a comparative analysis of the resulting galaxy–halo connection. We predict a total of |$300 ^{+75} _{-99}$| satellites with an absolute V -band magnitude (MV) less than 0 within 300 kpc from our Milky Way analogues. The fraction of subhaloes that host a galaxy at least this bright drops to 50 per cent by a halo peak mass of ∼8.9 × 107 M⊙, consistent with the occupation fraction inferred from the latest observations of Milky Way satellite population. [ABSTRACT FROM AUTHOR]

Stellar chemical abundances have proved themselves a key source of information for understanding the evolution of the Milky Way, and the scale of major stellar surveys such as GALAH have massively increased the amount of chemical data available. However, progress is hampered by the level of precision in chemical abundance data as well as the visualization methods for comparing the multidimensional outputs of chemical evolution models to stellar abundance data. Machine learning methods have greatly improved the former; while the application of tree-building or phylogenetic methods borrowed from biology are beginning to show promise with the latter. Here, we analyse a sample of GALAH solar twins to address these issues. We apply The Cannon algorithm to generate a catalogue of about 40 000 solar twins with 14 high precision abundances which we use to perform a phylogenetic analysis on a selection of stars that have two different ranges of eccentricities. From our analyses, we are able to find a group with mostly stars on circular orbits and some old stars with eccentric orbits whose age–[Y/Mg] relation agrees remarkably well with the chemical clocks published by previous high precision abundance studies. Our results show the power of combining survey data with machine learning and phylogenetics to reconstruct the history of the Milky Way. [ABSTRACT FROM AUTHOR]

The warp is a well-known undulation of the Milky Way disc. Its structure has been widely studied, but only since Gaia DR2 has it been possible to reveal its kinematic signature beyond the solar neighbourhood. In this work, we present an analysis of the warp traced by Classical Cepheids by means of a Fourier decomposition of their height (Z) and, for the first time, of their vertical velocity (Vz). We find a clear but complex signal that in both variables reveals an asymmetrical warp. In Z , we find the warp to be almost symmetric in amplitude at the disc's outskirts, with the two extremes never being diametrically opposed at any radius and the line of nodes presenting a twist in the direction of stellar rotation for R > 11 kpc. For Vz , in addition to the usual m  = 1 mode, an m  = 2 mode is needed to represent the kinematic signal of the warp, reflecting its azimuthal asymmetry. The line of maximum vertical velocity is similarly twisted as the line of nodes and trails behind by ≈25°. We develop a new formalism to derive the pattern speed and change in amplitude with time |$\dot{A}$| of each Fourier mode at each radius, via a joint analysis of the Fourier decomposition in Z and Vz. By applying it to the Cepheids we find, for the m  = 1 mode, a constant pattern speed in the direction of stellar rotation of 9.2 ± 3.1 km s−1 kpc−1, a negligible |$\dot{A}$| up to R ≈ 14 kpc and a slight increase at larger radii, in agreement with previous works. [ABSTRACT FROM AUTHOR]

This study focuses on Pristine |$\_180956.78$| −294759.8 (hereafter P180956, [Fe/H] = −1.95 ± 0.02), a star selected from the Pristine Inner Galaxy Survey (PIGS), and followed-up with the recently commissioned Gemini High-resolution Optical SpecTrograph (GHOST) at the Gemini South telescope. The GHOST spectrograph's high efficiency in the blue spectral region (3700−4800 Å) enables the detection of elemental tracers of early supernovae (e.g. Al, Mn, Sr, and Eu). The star exhibits chemical signatures resembling those found in ultrafaint dwarf (UFD) systems, characterized by very low abundances of neutron-capture elements (Sr, Ba, and Eu), which are uncommon among stars in the Milky Way halo. Our analysis suggests that P180956 bears the chemical imprints of a small number (2 or 4) of low-mass hypernovae (⁠|$\sim 10{-}15{\rm \, M_\odot }$|⁠), which are needed to mostly reproduce the abundance pattern of the light-elements (e.g. [Si, Ti/Mg, Ca] ∼0.6), and one fast-rotating intermediate-mass supernova (⁠|$\sim 300{\rm \, km \ s^{-1}}$|⁠ , |$\sim 80{-}120{\rm \, M_\odot }$|⁠), which is the main channel contributing to the high [Sr/Ba] (∼+1.2). The small pericentric (⁠|$\sim 0.7{\rm \, kpc}$|⁠) and apocentric (⁠|$\sim 13{\rm \, kpc}$|⁠) distances and its orbit confined to the plane (⁠|$\lesssim 2{\rm \, kpc}$|⁠) indicate that this star was likely accreted during the early Galactic assembly phase. Its chemo-dynamical properties suggest that P180956 formed in a system similar to a UFD galaxy accreted either alone, as one of the low-mass building blocks of the proto-Galaxy, or as a satellite of Gaia–Sausage–Enceladus. The combination of Gemini's large aperture with GHOST's high efficiency and broad spectral coverage makes this new spectrograph one of the leading instruments for near-field cosmology investigations. [ABSTRACT FROM AUTHOR]

The relative enrichment of s -process to α-elements ([ s /α]) has been linked with age, providing a potentially useful avenue in exploring the Milky Way's chemical evolution. However, the age–[ s /α] relationship is non-universal, with dependencies on metallicity and current location in the Galaxy. In this work, we examine these chemical clock tracers across birth radii (⁠|${R}_\text{birth}$|⁠), recovering the inherent trends between the variables. We derive |${R}_\text{birth}$| and explore the [ s /α]–age– |${R}_\text{birth}$| relationship for 36 652 APOGEE DR17 red giant and 24 467 GALAH DR3 main-sequence turn-off and subgiant branch disc stars using [Ce/Mg], [Ba/Mg], and [Y/Mg]. We discover that the age– |$\rm [{\it s}/Mg]$| relation is strongly dependent on birth location in the Milky Way, with stars born in the inner disc having the weakest correlation. This is congruent with the Galaxy's initially weak, negative |$\rm [{\it s}/Mg]$| radial gradient, which becomes positive and steep with time. We show that the non-universal relations of chemical clocks is caused by their fundamental trends with |${R}_\text{birth}$| over time, and suggest that the tight age– |$\rm [{\it s}/Mg]$| relation obtained with solar-like stars is due to similar |${R}_\text{birth}$| for a given age. Our results are put into context with a Galactic chemical evolution model, where we demonstrate the need for data-driven nucleosynthetic yields. [ABSTRACT FROM AUTHOR]

We select 1052 469 (754 635) thin disc stars from Gaia eDR3 and LAMOST DR7 in the range of Galactocentric radius R (guiding centre radius R g) from 8 to 11 kpc to investigate the asymmetries between the North and South of the disc mid-plane. More specifically, we analyse the vertical velocity dispersion profiles (⁠|$\sigma _{v_{z}}(z$|⁠)) in different bins of R (R g) and [Fe/H]. We find troughs in the profiles of |$\sigma _{v_{z}}(z)$| located in both the North (z ∼ 0.7 kpc) and South (z ∼ −0.5 kpc) of the disc at all radial and chemical bins studied. The difference between the Northern and Southern vertical velocity dispersion profiles (⁠|$\Delta \sigma _{v_{z}}(|z|)$|⁠) shows a shift between curves of different R and R g. A similar shift exists in these North–South (NS) asymmetry profiles further divided into different [Fe/H] ranges. The sample binned with R g more clearly displays the features in the velocity dispersion profiles. The shift in the peaks of the |$\Delta \sigma _{v_{z}}$| profiles and the variation in the phase spiral shape binned by metallicity indicate the variation of the vertical potential profiles and the radial metallicity gradient. The wave-like signal in NS asymmetry of |$\sigma _{v_{z}}(z)$| largely originates from phase spiral; while the NS asymmetry profiles of [Fe/H] only display a weak wave-like feature near solar radius. We perform a test particle simulation to qualitatively reproduce the observed results. A quantitative explanation of the NS asymmetry in the metallicity profile needs careful consideration of the spiral shape and the perturbation model, and we leave this for future work. [ABSTRACT FROM AUTHOR]

We use numerical simulations to model Gaia  DR3 data with the aim of constraining the Milky Way (MW) bar and spiral structure parameters. We show that both the morphology and the velocity field in MW-like galactic disc models are strong functions of time, changing dramatically over a few tens of Myr. This suggests that by finding a good match to the observed radial velocity field, vR (x, y), we can constrain the bar-spiral orientation. Incorporating uncertainties into our models is necessary to match the data; most importantly, a heliocentric distance uncertainty above 10–15 per cent distorts the bar's shape and vR quadrupole pattern morphology, and decreases its apparent angle with respect to the Sun-Galactocentric line. An excellent match to the Gaia   DR3 vR (x, y) field is found for a simulation with a bar length Rb ≈ 3.6 kpc. We argue that the data are consistent with an MW bar as short as ∼3 kpc, for moderate strength inner disc spiral structure (A 2/ A 0 ≈ 0.25) or, alternatively, with a bar length up to ∼5.2 kpc, provided that spiral arms are quite weak (A 2/ A 0 ≈ 0.1), and is most likely in the process of disconnecting from a spiral arm. We demonstrate that the bar angle and distance uncertainty can similarly affect the match between our models and the data – a smaller bar angle (20° instead of 30°) requires smaller distance uncertainty (20 per cent instead of 30 per cent) to explain the observations. Fourier components of the face-on density distribution of our models suggest that the MW does not have strong m = 1 and/or m  = 3 spirals near the solar radius. [ABSTRACT FROM AUTHOR]

Satellite galaxies within the Milky Way's (MW's) virial radius R vir are typically devoid of cold gas due to ram pressure stripping by the MW's corona. The density of this corona is poorly constrained today and essentially unconstrained in the past, but can be estimated using ram pressure stripping. In this paper, we probe the MW's corona at z ≈ 1.6 using the Draco dwarf spheroidal galaxy. We assume that (i) Draco's orbit is determined by its interaction with the MW, whose dark matter halo we evolve in time following cosmologically motivated prescriptions, (ii) Draco's star formation was quenched by ram pressure stripping and (iii) the MW's corona is approximately smooth, spherical, and in hydrostatic equilibrium. We used Gaia proper motions to set the initial conditions and Draco's star formation history to estimate its past gas content. We found indications that Draco was stripped of its gas during the first pericentric passage. Using 3D hydrodynamical simulations at a resolution that enables us to resolve individual supernovae and assuming no tidal stripping, which we estimate to be a minor effect, we find a density of the MW corona ≥8 × 10−4 cm−3 at a radius ≈0.72 R vir. This provides evidence that the MW's corona was already in place at z ≈ 1.6 and with a higher density than today. If isothermal, this corona would have contained all the baryons expected by the cosmological baryon fraction. Extrapolating to today shows good agreement with literature constraints if feedback has removed ≲30 per cent of baryons accreted on to the halo. [ABSTRACT FROM AUTHOR]

We present a new scheme for the classification of the in-situ and accreted globular clusters (GCs). The scheme uses total energy E and z -component of the orbital angular momentum and is calibrated using the [Al/Fe] abundance ratio. We demonstrate that this classification results in two GC populations with distinct spatial, kinematic, and chemical abundance distributions. The in-situ GCs are distributed within the central 10 kpc of the Galaxy in a flattened configuration aligned with the Milky Way (MW) disc, while the accreted GCs have a wide distribution of distances and a spatial distribution close to spherical. In-situ and accreted GCs have different |$\rm [Fe/H]$| distributions with the well-known bimodality present only in the metallicity distribution of the in-situ GCs. Furthermore, the accreted and in-situ GCs are well separated in the plane of |$\rm [Al/Fe]-[Mg/Fe]$| abundance ratios and follow distinct sequences in the age– |$\rm [Fe/H]$| plane. The in-situ GCs in our classification show a clear disc spin-up signature – the increase of median V ϕ at metallicities −1.3 < [Fe/H] < −1 similar to the spin-up in the in-situ field stars. This signature signals the MW's disc formation, which occurred ≈11.7−12.7 Gyr ago (or at z ≈ 3.1−5.3) according to in-situ GC ages. In-situ GCs with metallicities of |$\rm [Fe/H]\gtrsim -1.3$| were thus born in the MW disc, while lower metallicity in-situ GCs were born during early, turbulent, pre-disc stages of the evolution of the Galaxy and are part of its Aurora stellar component. [ABSTRACT FROM AUTHOR]

We use high-resolution zoom-in simulations to study the fueling of central galaxies by filamentary and diffuse accretion at redshifts, z ≳ 2. The parent haloes were chosen with similar total masses, log (M vir/M⊙) ∼ 11.75 ± 0.05, at z  = 6, 4, and 2, in high/low overdensity environments. We analyse the kinematic and thermodynamic properties of circumgalactic medium (CGM) within few virial radii, R vir, and down to the central galaxy. Using a hybrid d-web/entropy method we mapped the gaseous filaments, and separated inflows from outflows. We find that (1) The CGM is multiphase and not in thermodynamic or dynamic equilibrium; (2) filamentary and diffuse accretion rates and densities decrease with lower redshifts, and inflow velocities decrease from |$200-300\, {\rm {km\, s}^{-1}}$| by a factor of 2; (3) temperature within the filaments increases inside R vir, faster at lower redshifts; (4) filaments show a complex structure along their spines: a core radial flow surrounded by a lower density envelope. The cores exhibit elevated densities and lower temperature, with no obvious metallicity gradient in the cross sections. Filaments also tend to separate into different infall velocity regions and split density cores, thus producing a spaghetti-type flow; (6) inside the inner |$\sim 30\, h^{-1}$|  kpc, filaments develop the Kelvin–Helmholtz instability which ablates and dissolves them, and triggers turbulence along the filaments, clearly delineating their spines; (7) finally, the galactic outflows affect mostly the inner ∼0.5 R vir ∼ 100 h −1 kpc of the CGM. [ABSTRACT FROM AUTHOR]

The kinematics, mean α-elements, and ages of the stars reveal two distinct disc populations observed in the Solar neighbourhood. Although several studies have been carried out to characterize abundances of chemical elements and kinematic identifications in these populations, there are few studies that dealt with the analysis of stellar rotation and characteristics of populations focused on low-mass stars of the F and G types. In this work, we propose a new approach to classify stellar populations from the thin and thick disc for F and G dwarf stars by relating chemical abundance and rotation (v  sin  i) of these stars in the rotational–chemical–plane (v  sin  i –[α/Fe]–[Fe/H] plot). Our results show that the two rotational-chemical sequences disentangled, high- and low-[α/Fe] components, present properties (α-enhanced, ages, kinematics, space velocity gradients, local density ratio) closely linked to the stars of the chemically and kinematically defined as thin and thick disc components. This study has an impact on calibrations that relate ages to rotation (gyro-chronology) as well as age to chemical abundances (chemochronology), as well as sheds light on the understanding of the processes of Galactic disc formation. [ABSTRACT FROM AUTHOR]

Hot subdwarf B (sdB) stars are helium core burning stars that have lost almost their entire hydrogen envelope due to binary interaction. Their assumed canonical mass of |${M_{\mathrm{sdB}}}\sim 0.47$|  M⊙ has recently been debated given a broad range found both from observations as well as from the simulations. Here, we revise and refine the mass range for sdBs derived two decades ago with the Eggleton code, using the stellar evolution code mesa , and discuss the effects of metallicity and the inclusion of core overshooting during the main sequence. We find an excellent agreement for low-mass progenitors, up to ∼2.0 M⊙. For stars more massive than ∼2.5 M⊙, we obtain a wider range of sdB masses compared to the simulations from the literature. Our mesa  models for the lower metallicity predict, on average, slightly more massive sdBs. Finally, we show the results for the sdB lifetime as a function of sdB mass and discuss the effect this might have in the comparison between simulations and observational samples. This study paves the way for reproducing the observed Galactic mass distribution of sdB binaries. [ABSTRACT FROM AUTHOR]

Based on 4098 very metal poor (VMP) stars with six-dimensional phase-space and chemical information from Gaia Data Release 3 (DR3) and Large sky Area Multi-Object fiber Spectroscopic Telescope DR9 as tracers, we apply an unsupervised machine learning algorithm, shared nearest neighbour, to identify stellar groups in the action–energy (J–E) space. We detect seven previously known mergers in local samples, including Helmi Stream, Gaia-Sausage/Enceladus (GSE), metal-weak thick disc (MWTD), Pontus, Wukong, Thamnos, and I'itoi + Sequoia + Arjuna. According to energy, we further divide GSE and Wukong into smaller parts to explore the orbital characteristics of individual fragments. Similarly, the division of Thamnos is based on action. It can be found that the apocentric distances of GSE parts of high and medium energy levels are located at |$29.5\pm 3.6$| and |$13.0\pm 2.7\ {\rm kpc}$|⁠ , respectively, which suggests that GSE could account for breaks in the density profile of the Galactic halo at both ≈30 and 15–18 kpc. The VMP stars of MWTD move along prograde orbits with larger eccentricities than those of its more metal-rich stars, which indicates that the VMP part of MWTD may be formed by accreting with dwarf galaxies. Finally, we summarize all substructures discovered in our local VMP samples. Our results provide a reference for the formation and evolution of the inner halo of the Milky Way. [ABSTRACT FROM AUTHOR]

Observational studies are finding stars believed to be relics of the earliest stages of hierarchical mass assembly of the Milky Way (i.e. proto-galaxy). In this work, we contextualize these findings by studying the masses, ages, spatial distributions, morphology, kinematics, and chemical compositions of proto-galaxy populations from the 13 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological zoom-in simulations. Our findings indicate that proto-Milky Way populations: (i) can have a stellar mass range between 1 × 108 < M⋆ < 2 × 1010 [M⊙], a virial mass range between 3 × 1010 < M⋆ < 6 × 1011 [M⊙], and be as young as 8 ≲ Age ≲ 12.8 [Gyr] (1 ≲ z ≲ 6); (ii) are pre-dominantly centrally concentrated, with |$\sim 50~{{\ \rm per\ cent}}$| of the stars contained within 5–10 kpc; (iii) on average show weak but systematic net rotation in the plane of the host's disc at z  = 0 (i.e. 0.25 ≲ 〈κ/κdisc〉 ≲ 0.8); (iv) present [α/Fe]-[Fe/H] compositions that overlap with the metal-poor tail of the host's old disc; and (v) tend to assemble slightly earlier in Local Group-like environments than in systems in isolation. Interestingly, we find that |$\sim 60~{{\ \rm per\ cent}}$| of the proto-Milky Way galaxies are comprised by 1 dominant system (1/5 ≲M⋆/M⋆, proto-MilkyWay≲ 4/5) and 4–5 lower mass systems (M⋆/M⋆, proto-MilkyWay≲ 1/10); the other |$\sim 40~{{\ \rm per\ cent}}$| are comprised by 2 dominant systems and 3–4 lower mass systems. These massive/dominant proto-Milky Way fragments can be distinguished from the lower mass ones in chemical-kinematic samples, but appear (qualitatively) indistinguishable from one another. Our results could help observational studies disentangle if the Milky Way formed from one or two dominant systems. [ABSTRACT FROM AUTHOR]

The presence of Gaia DR3 provides a large sample of stars with complete 6D information, offering a fertile ground for the exploration of stellar objects that were accreted to the Milky Way through ancient merger events. In this study, we developed a deep learning methodology to identify ex-situ stars within the Gaia DR3 catalogue. After two phases of training, our neural network (NN) model was capable of performing binary classification of stars based on input data consisting of 3D position and velocity, as well as actions. From the target sample of 27 085 748 stars, our NN model managed to identify 160 146 ex-situ stars. The metallicity distribution suggests that this ex-situ sample comprises multiple components but appears to be predominated by the Gaia-Sausage-Enceladus (GSE). We identified member stars of the Magellanic Clouds, Sagittarius, and 20 globular clusters throughout our examination. Furthermore, an extensive group of member stars from GSE, Thamnos, Sequoia, Helmi streams, Wukong, and Pontus were meticulously selected, constituting an ideal sample for the comprehensive study of substructures. Finally, we conducted a preliminary estimation to determine the proportions of ex-situ stars in the thin disc, thick disc, and halo, which resulted in percentages of 0.1 per cent, 1.6 per cent, and 63.2 per cent, respectively. As the vertical height from the Galactic disc and distance from the Galactic centre increased, there was a corresponding upward trend in the ex-situ fraction of the target sample. [ABSTRACT FROM AUTHOR]

We employ a sample of 135 873 RR Lyrae stars (RRLs) with precise photometric-metallicity and distance estimates from the newly calibrated P –ϕ31– R 21–[Fe/H] and Gaia G band P – R 21–[Fe/H] absolute magnitude–metallicity relations of Li et al. combined with available proper motions from Gaia EDR3, and 6955 systemic radial velocities from Gaia DR3 and other sources, in order to explore the chemistry and kinematics of the halo of the Milky Way (MW). This sample is ideally suited for characterization of the inner- and outer-halo populations of the stellar halo, free from the bias associated with spectroscopically selected probes, and for estimation of their relative contributions as a function of Galactocentric distance. The results of a Gaussian mixture model analysis of these contributions are broadly consistent with other observational studies of the halo, and with expectations from recent MW simulation studies. We apply the hdbscan clustering method to the specific energies and cylindrical actions (E , Jr, Jϕ, J z ), identifying 97 dynamically tagged groups (DTGs) of RRLs, and explore their associations with recognized substructures of the MW. The precise photometric-distance determinations (relative distance errors on the order of 5 per cent or better), and the resulting high-quality determination of dynamical parameters, yield highly statistically significant (low) dispersions of [Fe/H] for the stellar members of the DTGs compared to random draws from the full sample, indicating that they share common star-formation and chemical histories, influenced by their birth environments. [ABSTRACT FROM AUTHOR]

Pristine_183.6849 + 04.8619 (P1836849) is an extremely metal-poor ([Fe/H] = −3.3 ± 0.1) star on a prograde orbit confined to the Galactic disc. Such stars are rare and may have their origins in protogalactic fragments that formed the early Milky Way, in low-mass satellites accreted later, or forming in situ in the Galactic plane. Here, we present a chemo-dynamical analysis of the spectral features between 3700−11 000 Å from a high-resolution spectrum taken during Science Verification of the new Gemini High-resolution Optical SpecTrograph. Spectral features for many chemical elements are analysed (Mg, Al, Si, Ca, Sc, Ti, Cr, Mn, Fe, Ni), and valuable upper limits are determined for others (C, Na, Sr, Ba). This main sequence star exhibits several rare chemical signatures, including (i) extremely low metallicity for a star in the Galactic disc, (ii) very low abundances of the light α-elements (Na, Mg, Si) compared to other metal-poor stars, and (iii) unusually large abundances of Cr and Mn, where [Cr, Mn/Fe]NLTE > +0.5. A comparison to theoretical yields from supernova models suggests that two low-mass Population III objects (one 10 M⊙ supernova and one 17 M⊙ hypernova) can reproduce the abundance pattern well (reduced χ2 < 1). When this star is compared to other extremely metal-poor stars on quasi-circular, prograde planar orbits, differences in both chemistry and kinematics imply there is little evidence for a common origin. The unique chemistry of P1836849 is discussed in terms of the earliest stages in the formation of the Milky Way. [ABSTRACT FROM AUTHOR]

The variations of elemental abundance and their ratios along the Galactocentric radius result from the chemical evolution of the Milky Way discs. The |$\rm ^{12}C/^{13}C$| ratio in particular is often used as a proxy to determine other isotopic ratios, such as |$\rm ^{16}O/^{18}O$| and |$\rm ^{14}N/^{15}N$|⁠. Measurements of |$\rm ^{12}CN$| and |$\rm ^{13}CN$| (or |$\rm C^{15}N$|⁠) – with their optical depths corrected via their hyper-fine structure lines – have traditionally been exploited to constrain the Galactocentric gradients of the CNO isotopic ratios. Such methods typically make several simplifying assumptions (e.g. a filling factor of unity, the Rayleigh–Jeans approximation, and the neglect of the cosmic microwave background) while adopting a single average gas phase. However, these simplifications introduce significant biases to the measured |$\rm ^{12}C/^{13}C$| and |$\rm ^{14}N/^{15}N$|⁠. We demonstrate that exploiting the optically thin satellite lines of 12CN constitutes a more reliable new method to derive |$\rm ^{12}C/^{13}C$| and |$\rm ^{14}N/^{15}N$| from CN isotopologues. We apply this satellite-line method to new IRAM 30-m observations of 12CN, 13CN, and C15N N  = 1 → 0 towards 15 metal-poor molecular clouds in the Galactic outer disc (R gc > 12 kpc), supplemented by data from the literature. After updating their Galactocentric distances, we find that |$\rm ^{12}C/^{13}C$| and |$\rm ^{14}N/^{15}N$| gradients are in good agreement with those derived using independent optically thin molecular tracers, even in regions with the lowest metallicities. We therefore recommend using optically thin tracers for Galactic and extragalactic CNO isotopic measurements, which avoids the biases associated with the traditional method. [ABSTRACT FROM AUTHOR]

The iMaNGA project uses a forward-modelling approach to compare the predictions of cosmological simulations with observations from SDSS-IV/MaNGA. We investigate the dependency of age and metallicity radial gradients on galaxy morphology, stellar mass, stellar surface mass density (Σ*), and environment. The key of our analysis is that observational biases affecting the interpretation of MaNGA data are emulated in the theoretical iMaNGA sample. The simulations reproduce the observed global stellar population scaling relations with positive correlations between galaxy mass and age/metallicity quite well and also produce younger stellar populations in late-type in agreement with observations. We do find interesting discrepancies, though, that can inform the physics and further development of the simulations. Ages of spiral galaxies and low-mass ellipticals are overestimated by about 2–4 Gyr. Radial metallicity gradients are steeper in iMaNGA than in MaNGA, a discrepancy most prominent in spiral and lenticular galaxies. Also, the observed steepening of metallicity gradients with increasing galaxy mass is not well matched by the simulations. We find that the theoretical radial profiles of surface mass density Σ* are steeper than in observations except for the most massive galaxies. In both MaNGA and iMaNGA [Z/H] correlates with Σ*, however, the simulations systematically predict lower [Z/H] by almost a factor of 2 at any Σ*. Most interestingly, for galaxies with stellar mass log  M * ≤ 10.80 M⊙, the MaNGA data reveal a positive correlation between galaxy radius and [Z/H] at fixed Σ*, which is not recovered in iMaNGA. Finally, the dependence on environmental density is negligible in both the theoretical iMaNGA and the observed MaNGA data. [ABSTRACT FROM AUTHOR]

An enormous number of compact binary systems, spanning from stellar to supermassive levels, emit substantial gravitational waves during their final evolutionary stages, thereby creating a stochastic gravitational wave background (SGWB). We calculate the merger rates of stellar compact binaries and massive black hole binaries using a semi-analytic galaxy formation model – Galaxy Assembly with Binary Evolution (GABE) – in a unified and self-consistent approach, followed by an estimation of the multiband SGWB contributed by those systems. We find that the amplitudes of the principal peaks of the SGWB energy density are within one order of magnitude ΩGW ∼ 10−9–10−8. This SGWB could easily be detected by the Square Kilometre Array (SKA), as well as by the planned interferometric detectors, such as the Einstein Telescope and the Laser Interferometer Space Antenna (LISA). The energy density of this background varies as ΩGW ∝  f 2/3 in the SKA band. The shape of the SGWB spectrum in the frequency range ∼[10−4,1] Hz could allow the LISA to distinguish the black hole seed models. The amplitude of the SGWB from merging stellar binary black holes at ∼100 Hz is approximately 10 and 100 times greater than those from merging binary neutron stars and neutron-star–black-hole mergers, respectively. Note that, since the cosmic star formation rate density predicted by GABE is somewhat lower than the observational results by ∼0.2 dex at z < ∼2, the amplitude of the SGWB in the frequency range ∼[1, 104] Hz may be underestimated by a similar factor at most. [ABSTRACT FROM AUTHOR]

We study the kinematics of stars both at their formation and today within 14 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological zoom-in simulations. We quantify the relative importance of cosmological disc settling and post-formation dynamical heating. We identify three eras: a Pre-Disc Era (typically ≳ 8 Gyr ago), when stars formed on dispersion-dominated orbits; an Early-Disc Era (≈8–4 Gyr ago), when stars started to form on rotation-dominated orbits but with high velocity dispersion, σform; and a Late-Disc Era (≲ 4 Gyr ago), when stars formed with low σform. σform increased with time during the Pre-Disc Era, peaking ≈8 Gyr ago, then decreased throughout the Early-Disc Era as the disc settled and remained low throughout the Late-Disc Era. By contrast, the dispersion measured today, σnow, increases monotonically with age because of stronger post-formation heating for Pre-Disc stars. Importantly, most of σnow was in place at formation, not added post-formation, for stars younger than ≈10 Gyr. We compare the evolution of the three velocity components: at all times, σR, form > σϕ, form > σZ, form. Post-formation heating primarily increased σR at ages ≲ 4 Gyr but acted nearly isotropically for older stars. The kinematics of young stars in FIRE-2 broadly agree with the range observed across the MW, M31, M33, and PHANGS-MUSE galaxies. The lookback time that the disc began to settle correlates with its dynamical state today: earlier-settling galaxies currently form colder discs. Including stellar cosmic-ray feedback does not significantly change disc rotational support at fixed stellar mass. [ABSTRACT FROM AUTHOR]

Recent discoveries of a supervirial hot phase of the Milky Way circumgalactic medium (CGM) have launched new questions regarding the multiphase structure of the CGM around the Galaxy. We use 1.05 Ms of archival Chandra /high-energy transmission grating observations to characterize highly ionized metal absorption at z  = 0 along the line of sight of the quasar NGC 3783. We detect two distinct temperature phases with T |$_1 = 5.83^{+0.15}_{-0.07}$| K, warm–hot virial temperature, and T |$_2=6.61^{+0.12}_{-0.06}$| K, hot supervirial temperature. The supervirial hot phase coexisting with the warm–hot virial phase has been detected in absorption along only two other sightlines and in one stacking analysis. There is scatter in temperature of the hot as well as warm–hot gas. Similar to previous observations, we detect supersolar abundance ratios of metals in the hot phase, with a Ne/O ratio 2σ above solar mixtures. These new detections continue the mystery of the mechanism behind the supervirial hot phase, but provide evidence that this is a true property of the CGM rather than an isolated observation. The supervirial CGM could hold the key to understanding the physical and chemical history of the Milky Way. [ABSTRACT FROM AUTHOR]

We investigate the formation (spin-up) of galactic discs in the artemis simulations of Milky Way (MW)-mass galaxies. In almost all galaxies, discs spin up at higher [Fe/H] than the MW. Those galaxies that contain an analogue of the Gaia Sausage-Enceladus (GSE) spin up at a lower average metallicity than those without. We identify six galaxies with spin-up metallicity similar to that of the MW, which formed their discs ∼8–11 Gyr ago. Five of these experience a merger similar to the GSE. The spin-up times correlate with the halo masses at early times: galaxies with early spin-up have larger virial masses at a lookback time t L = 12 Gyr. The fraction of stars accreted from outside the host galaxy is smaller in galaxies with earlier spin-ups. Accreted fractions small enough to be comparable to the MW are only found in galaxies with the earliest disc formation and large initial virial masses (M 200c ≈ 2 × 1011 M⊙ at t L = 12 Gyr). We find that discs form when the halo's virial mass reaches a threshold of M 200c ≈ (6 ± 3) × 1011 M⊙, independent of the spin-up time. However, the failure to form a disc in other galaxies appears to be instead related to mergers at early times. We also find that discs form when the central potential is not particularly steep. Our results indicate that the MW assembled its mass and formed its disc earlier than the average galaxy of a similar mass. [ABSTRACT FROM AUTHOR]

We present a comprehensive study of the molecular gas properties of 17 Type 2 quasars at z < 0.2 from the Quasar Feedback Survey (L |$_{\rm [O~{\small III}]}$| > 1042.1  |$\rm ergs^{-1}$|⁠), selected by their high [O  iii ] luminosities and displaying a large diversity of radio jet properties, but dominated by LIRG-like galaxies. With these data, we are able to investigate the impact of AGN and AGN feedback mechanisms on the global molecular interstellar medium. Using Atacama Pathfinder EXperiment and ALMA ACA observations, we measure the total molecular gas content using the CO(1-0) emission and homogeneously sample the carbon monoxide (CO) spectral line energy distributions, observing CO transitions (J up  = 1, 2, 3, 6, 7). We observe high r 21 ratios (r 21  =  L 'CO(2-1)/ L 'CO(1-0)) with a median r 21  = 1.06, similar to local (U)LIRGs (with r 21 ∼ 1) and higher than normal star-forming galaxies (with r 21 ∼ 0.65). Despite the high r 21 values, for the seven targets with the required data, we find low excitation in CO(6-5) & CO(7-6) (r 61 and r 62 < 0.6 in all but one target), unlike high-redshift quasars in the literature, which are far more luminous and show higher line ratios. The ionized gas traced by [O  iii ] exhibits systematically higher velocities than the molecular gas traced by CO. We conclude that any effects of quasar feedback (e.g. via outflows and radio jets) do not have a significant instantaneous impact on the global molecular gas content and excitation and we suggest that it only occurs on more localized scales. [ABSTRACT FROM AUTHOR]

Individual abundances in the Milky Way disc record stellar birth properties [e.g. age, birth radius (R birth)], and capture the diversity of the star-forming environments over time. Assuming an analytical relationship between ([Fe/H], [α/Fe]), and R birth, we examine the distributions of individual abundances [X/Fe] of elements C, O, Mg, Si, Ca (α), Al (odd-z), Mn (iron-peak), Y, and Ba (neutron-capture) for stars in the Milky Way. We want to understand how these elements might differentiate environments across the disc. We assign tracks of R birth in the [α/Fe] versus [Fe/H] plane as informed by expectations from simulations for ∼59 000 GALAH stars in the solar neighborhood (R ∼ 7−9 kpc) which also have inferred ages. Our formalism for R birth shows that older stars (∼10 Gyrs) have an R birth distribution with smaller mean values (i.e. |$\bar{R}_{\mbox{birth}} \sim 5\pm 0.8$| kpc) compared to younger stars (∼6 Gyrs; |$\bar{R}_{\mbox{birth}} \sim 10\pm 1.5$| kpc), for a given [Fe/H], consistent with inside–out growth. The α-, odd-z, and iron-peak element abundances decrease as a function of R birth, whereas the neutron-capture abundances increase. The R birth–[Fe/H] gradient we measure is steeper compared to the present-day gradient (−0.066 dex kpc−1 versus −0.058 dex kpc−1), which we also find true for R birth–[X/Fe] gradients. These results (i) showcase the feasibility of relating the birth radius of stars to their element abundances, (ii) demonstrate that the Milky Way abundance gradients across R birth have evolved to be shallower over time, and (iii) offer an observational comparison to element abundance distributions in hydrodynamical simulations. [ABSTRACT FROM AUTHOR]

We explore the influence of the Milky Way galaxy's chemical evolution on the formation, structure, and habitability of the Solar system. Using a multizone Galactic chemical evolution (GCE) model, we successfully reproduce key observational constraints, including the age–metallicity ([Fe/H]) relation, metallicity distribution functions, abundance gradients, and [X/Fe] ratio trends for critical elements involved in planetary mineralogy, including C, O, Mg, and Si. Our GCE model suggests that the Sun formed in the inner Galactic disc, R birth,⊙ ≈ 5 kpc. We also combined a stoichiometric model with the GCE model to examine the temporal evolution and spatial distribution of planet building blocks within the Milky Way galaxy, revealing trends in the condensed mass fraction (f cond), iron-to-silicon mass fraction (f iron), and water mass fraction (f water) over time and towards the inner Galactic disc regions. Specifically, our model predicts a higher f cond in the protoplanetary disc within the inner regions of the Milky Way galaxy, as well as an increased f iron and a decreased f water in the inner regions. Based on these findings, we discuss the potential impact of the Sun's birth location on the overall structure and habitability of the Solar system. [ABSTRACT FROM AUTHOR]