Search

The cold neutral medium (CNM) is where neutral atomic hydrogen (H i) is converted into molecular clouds, so the structure and kinematics of the CNM are key drivers of galaxy evolution. Here we provide new constraints on the vertical distribution of the CNM using the recently developed kinematic_scaleheight software package and a large catalog of sensitive H i absorption observations. We estimate the thickness of the CNM in the solar neighborhood to be σ z ∼ 50–90 pc, assuming a Gaussian vertical distribution. This is a factor of ∼2 smaller than typically assumed, indicating that the thickness of the CNM in the solar neighborhood is similar to that found in the inner Galaxy, consistent with recent simulation results. If we consider only structures with H i optical depths τ > 0.1 or column densities N (H i) > 1019.5 cm−2, which recent work suggests are thresholds for molecule formation, we find σ z ∼ 50 pc. Meanwhile, for structures with τ < 0.1 or column densities N (H i) < 1019.5 cm−2, we find σ z ∼ 120 pc. These thicknesses are similar to those derived for the thin- and thick-disk molecular cloud populations traced by CO emission, possibly suggesting that cold H i and CO are well mixed. Approximately 20% of CNM structures are identified as outliers, with kinematics that are not well explained by Galactic rotation. We show that some of these CNM structures—perhaps representing intermediate-velocity clouds—are associated with the Local Bubble wall. We compare our results to recent observations and simulations, and we discuss their implications for the multiphase structure of the Milky Way's interstellar medium. [ABSTRACT FROM AUTHOR]

Stars in the Galactic disk are born on cold, nearly circular orbits with small vertical excursions. After their birth, their orbits evolve, driven by small- or large-scale perturbations in the Galactic disk's gravitational potential. Here, we study the vertical motions of young stars over their first few orbital periods, using a sample of OBA stars from Gaia E/DR3, which includes radial velocities and ages τ from LAMOST. We construct a parametric model for the time evolution of the stellar orbits' mean vertical actions J z as a function of Galactocentric radius, R GC. Accounting for data uncertainties, we use Markov Chain Monte Carlo analysis in annuli of Galactocentric radius (R GC) to constrain the model parameters. Our best-fit model shows a remarkably linear increase of vertical actions with age across all Galactocentric radii examined. Orbital heating by random scattering could offer a straightforward interpretation for this trend. However, various other dynamical aspects of the Galactic disk, such as stars being born in a warped disk, might offer alternative explanations that could be tested in the future. [ABSTRACT FROM AUTHOR]

The Milky Way is widely considered to exhibit features of a rotational bar or quadrupole bar. In either case, the feature of the resonance of the Galactic bar should be present in the properties of the chemistry and kinematics, over a large area of the disk. With a sample of over 170,000 red clump stars from LAMOST-APOGEE data, we attempt to detect the chemical and kinematic signatures of the resonances of the Galactic bar, within 4.0 ≤ R ≤ 15.0 kpc and ∣ Z ∣ ≤ 3.0 kpc. The measurement of the Δ[Fe/H]/Δ∣ Z ∣– R with subtracted global profiles trends shows that the thin and thick disks values are Cor_Δ[Fe/H]/Δ∣ Z ∣ = 0.010 sin (1.598 R + 2.551) and Cor_Δ[Fe/H]/Δ∣ Z ∣ = 0.006 sin (1.258 R − 0.019), respectively. The analysis of the tilt angle of the velocity ellipsoid indicates that the thin and thick disks are accurately described as α = α _0 arctan ( Z / R ), with ${\alpha }_{0}=0.198\,\sin $ (0.853 R + 1.982) + 0.630 and ${\alpha }_{0}=0.220\,\sin $ (0.884 R + 2.012) + 0.679 for thin and thick disks, respectively. These periodic oscillations in Cor_Δ[Fe/H]/Δ∣ Z ∣ and α _0 with R appear in both thin and thick disks are the most likely chemical and kinematic signatures of the resonance of the Galactic bar. The difference in the phase of the functions of the fitted periodic oscillations for the thin and thick disks may be related to the presence of a second Galactic bar.

The Parkes 20 cm multibeam pulsar surveys have discovered nearly half of the known pulsars and revealed many distant pulsars with high dispersion measures. Using a sample of 1301 pulsars from these surveys, we have explored the spatial distribution and birth rate of normal pulsars. The pulsar distances used to calculate the pulsar surface density are estimated from the YMW16 electron-density model. When estimating the impact of the Galactic background radiation on our survey, we projected pulsars in the galaxy onto the Galactic plane, assuming that the flux density distribution of pulsars is uniform in all directions, and utilized the most up-to-date background temperature map. We also used an up-to-date version of the ATNF Pulsar Catalogue to model the distribution of pulsar flux densities at 1400 MHz. We derive an improved radial distribution for the pulsar surface density projected onto the Galactic plane, which has a maximum value at ∼4 kpc from the Galactic center. We also derive the local surface density and birth rate of pulsars, obtaining 47 ± 5 kpc ^−2 and ∼4.7 ± 0.5 kpc ^−2 Myr ^−1 , respectively. For the total number of potentially detectable pulsars in the galaxy, we obtain (1.1 ± 0.2) × 10 ^4 and (1.1 ± 0.2)×10 ^5 before and after applying the Tauris & Manchester beaming correction model. The radial distribution function is used to estimate the proportion of pulsars in each spiral arm and the Galactic center.

The phase space, as coded by kinematic parameters and chemical abundances, is crucial for understanding the formation of the Galactic disk. Using red giant stars from the Galactic thin disk with [Fe/H] > − 0.8 and low- α ratios identified in LAMOST-Gaia, we detect numerous ridges and undulations in the R – V _ϕ diagram coded by median V _R . Strikingly, the slope of these features changes from −22 km s ^−1 kpc ^−1 to −8 km s ^−1 kpc ^−1 at R ∼ 11.5 kpc. Accordingly, the R _g – V _ϕ plane, also coded by median V _R , reveals wave-like structures that propagate outwards in the inner disk but reverse direction and move inwards beyond R _g = 11.5 kpc. The most prominent feature is the G1 group, distinguished by its wider spread and negative median V _R at R _g ∼ 15 kpc, contrasting with the narrower G0 group that exhibits positive median V _R at R _g < 11.5 kpc. Furthermore, the [C/N] versus [Fe/H] relationship for the G1 group mirrors the opposite trend compared to the G0 group. Since [C/N] serves as a proxy for age, this contrasting behavior suggests an inverse age–metallicity relation for the G1 group. Comparison with the Sagittarius dwarf galaxy reveals that the G1 group possesses distinct [Mg/Fe] and [Al/Fe] ratios, yet its [C/N] versus [Fe/H] pattern is similar to that of the Sgr dwarf galaxy. Based on these observations, we proposed that the inward-moving and asymmetric velocity wave G1 might be linked to the minor merge of the Sagittarius galaxy.

Detailed analysis of kinematics of the Milky Way disk in the solar neighborhood based on the GAIA DR3 catalog reveals the existence of peculiarities in the stellar velocity distribution perpendicular to the galactic plane. We study the influence of resonances—the outer Lindblad resonance and the outer vertical Lindblad resonance—of a rotating bar with stellar oscillations perpendicular to the plane of the disk, and their role in shaping the spatial and the velocity distributions of stars. We find that the Z and V Z distributions of stars with respect to L Z are affected by the outer Lindblad resonance. The existence of bar resonance with stellar oscillations perpendicular to the plane of the disk is demonstrated for a long (large semi-axis 5 kpc) and fast rotating bar with Ω b = 60.0 km s − 1 kpc − 1 . We show also that, in the model with the long and fast rotating bar, some stars in the 2:1 OLR region deviate far from their original places, entering the bar region. A combination of resonance excitation of stellar motions at the 2:1 OLR region together with strong interaction of the stars with the bar potential leads to the formation of the group of 'escapees', i.e., stars that deviate in R and Z—directions at large distances from the resonance region. Simulations, however, do not demonstrate any noticeable effect on V Z -distribution of stars in the solar neighborhood. [ABSTRACT FROM AUTHOR]

We present new maps of the Milky Way disk showing the distribution of metallicity ([Fe/H]), α -element abundances ([Mg/Fe]), and stellar age, using a sample of 66,496 red giant stars from the final data release (DR17) of the Apache Point Observatory Galactic Evolution Experiment survey. We measure radial and vertical gradients, quantify the distribution functions for age and metallicity, and explore chemical clock relations across the Milky Way for the low- α disk, high- α disk, and total population independently. The low- α disk exhibits a negative radial metallicity gradient of −0.06 ± 0.001 dex kpc−1, which flattens with distance from the midplane. The high- α disk shows a flat radial gradient in metallicity and age across nearly all locations of the disk. The age and metallicity distribution functions shift from negatively skewed in the inner Galaxy to positively skewed at large radius. Significant bimodality in the [Mg/Fe]–[Fe/H] plane and in the [Mg/Fe]–age relation persist across the entire disk. The age estimates have typical uncertainties of ∼0.15 in log(age) and may be subject to additional systematic errors, which impose limitations on conclusions drawn from this sample. Nevertheless, these results act as critical constraints on galactic evolution models, constraining which physical processes played a dominant role in the formation of the Milky Way disk. We discuss how radial migration predicts many of the observed trends near the solar neighborhood and in the outer disk, but an additional more dramatic evolution history, such as the multi-infall model or a merger event, is needed to explain the chemical and age bimodality elsewhere in the Galaxy. [ABSTRACT FROM AUTHOR]

We develop a linear perturbative formalism to compute the response of an inhomogeneous stellar disk embedded in a nonresponsive dark matter (DM) halo to various perturbations like bars, spiral arms, and encounters with satellite galaxies. Without self-gravity to reinforce it, the response of a Fourier mode phase mixes away due to an intrinsic spread in the vertical (Ω z ), radial (Ω r ), and azimuthal (Ω ϕ ) frequencies, triggering local phase-space spirals. The detailed galactic potential dictates the shape of phase spirals: phase mixing occurs more slowly and thus phase spirals are more loosely wound in the outer disk and in the presence of an ambient DM halo. Collisional diffusion due to scattering of stars by structures like giant molecular clouds causes superexponential damping of the phase spiral amplitude. The z – v z phase spiral is one-armed (two-armed) for vertically antisymmetric (symmetric) bending (breathing) modes. Only transient perturbations with timescales (τ P) comparable to the vertical oscillation period (τ z ∼ 1/Ω z ) can trigger vertical phase spirals. Each (n, l, m) mode of the response to impulsive (τ P < τ = 1/(n Ω z + l Ω r + m Ω ϕ )) perturbations is power-law (∼ τ P/ τ) suppressed, but that to adiabatic (τ P > τ) perturbations is exponentially weak ( ∼ exp − τ P / τ α ) except for resonant (τ → ∞) modes. Slower (τ P > τ z ) perturbations, e.g., distant encounters with satellite galaxies, induce stronger bending modes. Sagittarius (Sgr) dominates the solar neighborhood response of the Milky Way (MW) disk to satellite encounters. Thus, if the Gaia phase spiral was triggered by a MW satellite, Sgr is the leading contender. However, the survival of the phase spiral against collisional damping necessitates an impact ∼0.6–0.7 Gyr ago. [ABSTRACT FROM AUTHOR]

Using a complete sample of about 330,000 dwarf stars, well measured by Gaia DR3, limited to the galactic north and south solid angles | b | < 75 ° and up to a vertical distance of 2 kpc, we analyze the vertical structure of the Milky Way stellar disks, based on projected tangential velocities. From selected subsamples dominated by their corresponding population, we obtain the thin and thick disk scale heights as h Z = 279.76 ± 12.49 pc and H Z = 797.23 ± 12.34 pc, respectively. Then from the simultaneous fitting of the sum of two populations over the whole sample, assuming these scale heights, we estimate the thick-to-thin disk number density ratio at the galactic plane to be ρ T / ρ t = 0.750 ± 0.049 , which is consistent with a previous result by the authors: in the galactic plane there is a significant number of thick disk stars, possibly as many as thin disk ones, which also points to the existence of more thick disk stars than generally thought. The overall fit does not closely follow the data for | Z | > 700 pc and points to the presence of more stars beyond the thin disk that cannot be accounted for by the two-disk model. [ABSTRACT FROM AUTHOR]

© 2020. The American Astronomical Society. All rights reserved. Microlensing surveys have discovered thousands of events, with almost all events discovered within the Galactic bulge or toward the Magellanic Clouds. The Zwicky Transient Facility (ZTF), while not designed to be a microlensing campaign, is an optical time-domain survey that observes the entire northern sky, including the Galactic plane, every few nights. The ZTF observes ∼109 stars in the g and r bands and can significantly contribute to the observed microlensing population. We predict that the ZTF will observe ∼1100 microlensing events in 3 yr of observing within 10° latitude of the Galactic plane, with ∼500 events in the outer Galaxy (ℓ ≥ 10°). This yield increases to ∼1400 (∼800) events by combining every three ZTF exposures, ∼1800 (∼900) events if the ZTF observes for a total of 5 yr, and ∼2400 (∼1300) events for a 5 yr survey with postprocessing image stacking. Using the microlensing modeling software PopSyCLE, we compare the microlensing populations in the Galactic bulge and the outer Galaxy. We also present an analysis of the microlensing event ZTF18abhxjmj to demonstrate how to leverage these population statistics in event modeling. The ZTF will constrain Galactic structure, stellar populations, and primordial black holes through photometric microlensing.

Age is the most difficult fundamental stellar parameter to infer for isolated stars. While isochrone-based ages are in general imprecise for both main-sequence dwarfs and red giants, precise isochrone-based ages can be obtained for stars on the subgiant branch transitioning from core to shell hydrogen burning. We synthesize Gaia Data Release 3–based distance inferences, multiwavelength photometry from the ultraviolet to the mid-infrared, and three-dimensional extinction maps to construct a sample of 289,759 solar-metallicity stars amenable to accurate, precise, and physically self-consistent age inferences. Using subgiants in the solar-metallicity open clusters NGC 2682 (i.e., M67) and NGC 188, we show that our approach yields accurate and physically self-consistent ages and metallicities with median statistical precisions of 8% and 0.06 dex, respectively. The inclusion of systematic uncertainties resulting from nonsingle or variable stars results in age and metallicity precisions of 9% and 0.12 dex, respectively. We supplement this solar-metallicity sample with an additional 112,062 metal-poor subgiants, including over 3000 stars with [Fe/H] ≲ −1.50, 7% age precisions, and apparent Gaia G -band magnitudes G < 14. We further demonstrate that our inferred metallicities agree with those produced by multiplexed spectroscopic surveys. As an example of the scientific potential of this catalog, we show that the solar neighborhood star formation history has three components at ([Fe/H], τ /Gyr) ≈ (+0.0, 4), (+0.2, 7), and a roughly linear sequence in age–metallicity space beginning at ([Fe/H], τ /Gyr) ≈ (+0.2, 7) and extending to (−0.5, 13). Our analyses indicate that the solar neighborhood includes stars on disk-like orbits even at the oldest ages and lowest metallicities accessible by our samples.

Decoding the key dynamical processes that shape the Galactic disk structure is crucial for reconstructing the Milky Way’s evolution history. The second Gaia data release unveils a novel wave pattern in the L _Z −〈 V _R 〉 space, but its formation mechanism remains elusive due to the intricate nature of involved perturbations and the challenges in disentangling their effects. Utilizing the latest Gaia DR3 data, we find that the L _Z −〈 V _R 〉 wave systematically shifts toward lower L _Z for dynamically hotter stars with larger J _Z values. The amplitude of this phase shift between stars of different dynamical hotness (Δ L _Z ) peaks at around 2100 km s ^−1 kpc. To differentiate the role of different perturbations, we perform three sets of test particle simulations, wherein a satellite galaxy, transient spiral arms, and a bar plus the transient spiral arms act as the sole perturber, respectively. Under the satellite impact, the phase shift amplitude Δ L _Z decreases toward higher L _Z , which we interpret through a toy model of radial phase mixing. While neither the transient spiral arms nor the bar generates an azimuthally universal phase shift variation pattern, combining the bar and spirals generates a characteristic Δ L _Z peak at the 2:1 outer Lindblad resonance (OLR) of the bar, qualitatively resembling the observed feature. Therefore, the L _Z −〈 V _R 〉 wave is more likely of internal origin. Furthermore, linking the Δ L _Z peak to the 2:1 OLR offers a novel approach to constraining the pattern speed of the Galactic bar, supporting the long/slow bar model.

The cold neutral medium (CNM) is where neutral atomic hydrogen (H i ) is converted into molecular clouds, so the structure and kinematics of the CNM are key drivers of galaxy evolution. Here we provide new constraints on the vertical distribution of the CNM using the recently developed kinematic _ scaleheight software package and a large catalog of sensitive H i absorption observations. We estimate the thickness of the CNM in the solar neighborhood to be σ _z ∼ 50–90 pc, assuming a Gaussian vertical distribution. This is a factor of ∼2 smaller than typically assumed, indicating that the thickness of the CNM in the solar neighborhood is similar to that found in the inner Galaxy, consistent with recent simulation results. If we consider only structures with H i optical depths τ > 0.1 or column densities N (H i ) > 10 ^19.5 cm ^−2 , which recent work suggests are thresholds for molecule formation, we find σ _z ∼ 50 pc. Meanwhile, for structures with τ < 0.1 or column densities N (H i ) < 10 ^19.5 cm ^−2 , we find σ _z ∼ 120 pc. These thicknesses are similar to those derived for the thin- and thick-disk molecular cloud populations traced by CO emission, possibly suggesting that cold H i and CO are well mixed. Approximately 20% of CNM structures are identified as outliers, with kinematics that are not well explained by Galactic rotation. We show that some of these CNM structures—perhaps representing intermediate-velocity clouds—are associated with the Local Bubble wall. We compare our results to recent observations and simulations, and we discuss their implications for the multiphase structure of the Milky Way’s interstellar medium.

Molecular clouds are regions of dense gas and dust in space where new stars and planets are born. There is a strong correlation between the distribution of dust and molecular gas in molecular clouds. The present work focuses on the three-dimensional morphological comparisons between dust and gas in 567 molecular clouds identified in a previously published catalog. We confirm a sample of 112 molecular clouds, where the cloud morphology based on CO observations and dust observations displays good overall consistency. There are up to 334 molecular clouds whose dust distribution might be related to the distribution of gas. We are unable to find gas structures that correlate with the shape of the dust distribution in 24 molecular clouds. For the 112 molecular clouds where the dust distribution correlates very well with the distribution of gas, we use CO observational data to measure the physical properties of these molecular clouds and compare them with the results derived from dust, exploring the correlation between gas and dust in the molecular clouds. We found that the gas and dust in the molecular clouds have a fairly good linear relationship, with a gas-to-dust ratio of $\mathrm{GDR}=({2.80}_{-0.34}^{+0.37})\times {10}^{21}\,{\mathrm{cm}}^{-2}\,{\mathrm{mag}}^{-1}$ . The ratio varies considerably among different molecular clouds. We measured the scale height of dust-CO clouds exhibiting strong correlations, finding ${h}_{Z}={43.3}_{-3.5}^{+4.0}$ pc.

Previous studies have revealed that the Galactic warp is a long-lived, nonsteady, and asymmetric structure. There is a need for a model that accounts for the warp’s long-term evolution. Given that this structure has persisted for over 5 Gyr, its timeline may coincide with the completion of the Gaia–Sausage–Enceladus (GSE) merger. Recent studies indicate that the GSE, the significant merger of our Galaxy, was likely a gas-rich merger and the large amount of gas introduced could have created a profound impact on the Galactic morphology. This study utilizes GIZMO simulation code to construct a gas-rich GSE merger. By reconstructing the observed characteristics of the GSE, we successfully reproduce the disk warp and capture nearly all of its documented features, which align closely with observational data from both stellar and gas disks. This simulation demonstrates the possibility that a single major merger could generate the Galactic warp amplitude and precession. Furthermore, the analysis of the warp’s long-term evolution may offer more clues into the formation history of the Milky Way.

We present a novel, deep-learning-based method—dubbed Galactic-Seismology Substructures and Streams Hunter, or GS ^3 Hunter for short—to search for substructures and streams in stellar kinematics data. GS ^3 Hunter relies on a combined application of Siamese neural networks to transform the phase space information and the K -means algorithm for the clustering. As a validation test, we apply GS ^3 Hunter to a subset of the Feedback in Realistic Environments (FIRE) cosmological simulations. The stellar streams and substructures thus identified are in good agreement with corresponding results reported earlier by the FIRE team. In the same vein, we apply our method to a subset of local halo stars from the Gaia Early Data Release 3 and GALAH DR3 data sets and recover several previously known dynamical groups, such as Thamnos 1+2, the hot thick disk, ED-1, L-RL3, Helmi 1+2, Gaia-Sausage-Enceladus, Sequoia, Virgo Radial Merger, Cronus, and Nereus. Finally, we apply our method without fine-tuning to a subset of K giant stars located in the inner halo region, obtained from the LAMOST Data Release 5 data set. We recover three previously known structures (Sagittarius, Hercules-Aquila Cloud, and the Virgo Overdensity), but we also discover a number of new substructures. We anticipate that GS ^3 Hunter will become a useful tool for the community dedicated to the search for stellar streams and structures in the Milky Way (MW) and the Local Group, thus helping advance our understanding of the stellar inner and outer halos and the assembly and tidal stripping history in and around the MW.

Type Ia supernovae (SNe Ia) produce most of the Fe-peak elements in the Universe and therefore are a crucial ingredient in galactic chemical evolution models. SNe Ia do not explode immediately after star formation, and the delay-time distribution (DTD) has not been definitively determined by supernova surveys or theoretical models. Because the DTD also affects the relationship among age, [Fe/H], and [ α /Fe] in chemical evolution models, comparison with observations of stars in the Milky Way is an important consistency check for any proposed DTD. We implement several popular forms of the DTD in combination with multiple star formation histories for the Milky Way in multizone chemical evolution models that include radial stellar migration. We compare our predicted interstellar medium abundance tracks, stellar abundance distributions, and stellar age distributions to the final data release of the Apache Point Observatory Galactic Evolution Experiment. We find that the DTD has the largest effect on the [ α /Fe] distribution: a DTD with more prompt SNe Ia produces a stellar abundance distribution that is skewed toward a lower [ α /Fe] ratio. While the DTD alone cannot explain the observed bimodality in the [ α /Fe] distribution, in combination with an appropriate star formation history it affects the goodness of fit between the predicted and observed high- α sequence. Our model results favor an extended DTD with fewer prompt SNe Ia than the fiducial t ^−1 power law.

Diffuse ionized gas pervades the disk of the Milky Way. We detect extremely faint emission from this Galactic warm ionized medium (WIM) using the Green Bank Telescope to make radio recombination line (RRL) observations toward two Milky Way sight lines: G20, ( ℓ , b ) = (20°, 0°), and G45, ( ℓ , b ) = (45°, 0°). We stack 18 consecutive Hn α transitions between 4.3 and 7.1 GHz to derive 〈Hn α 〉 spectra that are sensitive to RRL emission from plasmas with emission measures EM ≳ 10 cm ^−6 pc. Each sight line has two Gaussian-shaped spectral components with emission measures that range between ∼100 and ∼300 cm ^−6 pc. Because there is no detectable RRL emission at negative LSR velocities, the emitting plasma must be located interior to the solar orbit. The G20 and G45 emission measures imply rms densities of 0.15 and 0.18 cm ^−3 , respectively, if these sight lines are filled with homogeneous plasma. The observed 〈Hn β 〉/〈Hn α 〉 line ratios are consistent with LTE excitation for the strongest components. The high-velocity component of G20 has a narrow line width, 13.5 km s ^−1 , that sets an upper limit of ≲4000 K for the plasma electron temperature. This is inconsistent with the ansatz of a canonically pervasive, low-density, ∼10,000 K WIM plasma.

Stars in the Galactic disk are born on cold, nearly circular orbits with small vertical excursions. After their birth, their orbits evolve, driven by small- or large-scale perturbations in the Galactic disk’s gravitational potential. Here, we study the vertical motions of young stars over their first few orbital periods, using a sample of OBA stars from Gaia E/DR3, which includes radial velocities and ages τ from LAMOST. We construct a parametric model for the time evolution of the stellar orbits’ mean vertical actions J _z as a function of Galactocentric radius, R _GC . Accounting for data uncertainties, we use Markov Chain Monte Carlo analysis in annuli of Galactocentric radius ( R _GC ) to constrain the model parameters. Our best-fit model shows a remarkably linear increase of vertical actions with age across all Galactocentric radii examined. Orbital heating by random scattering could offer a straightforward interpretation for this trend. However, various other dynamical aspects of the Galactic disk, such as stars being born in a warped disk, might offer alternative explanations that could be tested in the future.

We illustrate the formation and evolution of the Milky Way over cosmic time, utilizing a sample of 10 million red giant stars with full chemodynamical information, including metallicities and α -abundances from low-resolution Gaia XP spectra. The evolution of angular momentum as a function of metallicity—a rough proxy for stellar age, particularly for high-[ α /Fe] stars—displays three distinct phases: the disordered and chaotic protogalaxy, the kinematically hot old disk, and the kinematically cold young disk. The old high- α disk starts at [Fe/H] ≈ −1.0, “spinning up” from the nascent protogalaxy, and then exhibiting a smooth “cooldown” toward more ordered and circular orbits at higher metallicities. The young low- α disk is kinematically cold throughout its metallicity range, with its observed properties modulated by a strong radial gradient. We interpret these trends using Milky Way analogs from the TNG50 cosmological simulation, identifying one that closely matches the kinematic evolution of our galaxy. This halo’s protogalaxy spins up into a relatively thin and misaligned high- α disk at early times, which is subsequently heated and torqued by a major gas-rich merger. The merger contributes a large amount of low-metallicity gas and angular momentum, from which the kinematically cold low- α stellar disk is subsequently born. This simulated history parallels several observed features of the Milky Way, particularly the decisive Gaia–Sausage–Enceladus merger that likely occurred at z ≈ 2. Our results provide an all-sky perspective on the emerging picture of our galaxy’s three-phase formation, impelled by the three physical mechanisms of spinup, merger, and cooldown.

Using a representative sample of Milky Way (MW)–like galaxies from the TNG50 cosmological simulation, we investigate physical processes driving the formation of galactic disks. A disk forms as a result of the interplay between inflow and outflow carrying angular momentum in and out of the galaxy. Interestingly, the inflow and outflow have remarkably similar distributions of angular momentum, suggesting an exchange of angular momentum and/or outflow recycling, leading to continuous feeding of prealigned material from the corotating circumgalactic medium. We show that the disk formation in TNG50 is correlated with stellar bulge formation, in qualitative agreement with a recent theoretical model of disk formation facilitated by steep gravitational potentials. Disk formation is also correlated with the formation of a hot circumgalactic halo with around half of the inflow occurring at subsonic and transonic velocities corresponding to Mach numbers of ≲2. In the context of recent theoretical works connecting disk settling and hot halo formation, our results imply that the subsonic part of the inflow may settle into a disk while the remaining supersonic inflow will perturb this disk via the chaotic cold accretion. We find that disks tend to form when the host halos become more massive than ∼(1–2) × 10 ^11 M _⊙ , consistent with previous theoretical findings and observational estimates of the predisk protogalaxy remnant in the MW. Our results do not prove that either corotating outflow recycling, gravitational potential steepening, or hot halo formation cause disk formation, but they show that all these processes occur concurrently and may play an important role in disk growth.

Photometric stellar surveys now cover a large fraction of the sky, probe to fainter magnitudes than large-scale spectroscopic surveys, and are relatively free from the target selection biases often associated with such studies. Photometric-metallicity estimates that include narrow/medium-band filters can achieve comparable accuracy and precision to existing low-resolution spectroscopic surveys such as Sloan Digital Sky Survey/SEGUE and LAMOST. Here we report on an effort to identify likely members of the Galactic disk system among the very metal-poor (VMP; [Fe/H] ≤ −2) and extremely metal-poor (EMP; [Fe/H] ≤ −3) stars. Our analysis is based on an initial sample of ∼11.5 million stars with full space motions selected from the SkyMapper Southern Survey (SMSS) and Stellar Abundance and Galactic Evolution Survey (SAGES). After applying a number of quality cuts to obtain the best available metallicity and dynamical estimates, we analyze a total of ∼5.86 million stars in the combined SMSS/SAGES sample. We employ two techniques that, depending on the method, identify between 876 and 1476 VMP stars (6.9%−11.7% of all VMP stars) and between 40 and 59 EMP stars (12.4%−18.3% of all EMP stars) that appear to be members of the Galactic disk system on highly prograde orbits ( v _ϕ > 150 km s ^−1 ). The total number of candidate VMP/EMP disklike stars is 1496, the majority of which have low orbital eccentricities, ecc ≤ 0.4; many have ecc ≤ 0.2. The large fractions of VMP/EMP stars associated with the Milky Way disk system strongly suggest the presence of an early-forming “primordial” disk.

In this study, we estimate the stellar ages of LAMOST DR8 red giant branch (RGB) stars based on the gradient boosting decision tree (GBDT) algorithm. We used 2643 RGB stars extracted from the APOKASC-2 asteroseismological catalog as the training data set. After selecting the parameters ([ α /Fe], [C/Fe], T _eff , [N/Fe], [C/H], log g ) highly correlated with age using GBDT, we apply the same GBDT method to the new catalog of more than 590,000 stars classified as RGB stars. The test data set shows that the median relative error is around 11.6% for the method. We also compare the predicted ages of RGB stars with other studies (e.g., based on APOGEE) and find some systematic differences. The final uncertainty is about 15%–30% compared to the ages of open clusters. Then, we present the spatial distribution of the RGB sample with an age determination, which could recreate the expected result, and discuss systematic biases. All these diagnostics show that one can apply the GBDT method to other stellar samples to estimate atmospheric parameters and age.

The vertical distribution of cold neutral hydrogen (H i ) clouds is a constraint on models of the structure, dynamics, and hydrostatic balance of the interstellar medium. In 1978, Crovisier pioneered a method to infer the vertical distribution of H i absorbing clouds in the solar neighborhood. Using data from the Nançay 21 cm absorption survey, Crovisier determined the mean vertical displacement of cold H i clouds, 〈∣ z ∣〉. We revisit that author’s analysis and explore the consequences of truncating the H i absorption sample in Galactic latitude. For any nonzero latitude limit, we find that the quantity inferred by Crovisier is not the mean vertical displacement but rather a ratio involving higher moments of the vertical distribution. The resultant distribution scale heights are thus ∼1.5 to ∼3 times smaller than previously determined. In light of this discovery, we develop a Bayesian Monte Carlo Markov Chain method to infer the vertical distribution of H i absorbing clouds. We fit our model to the original Nançay data and find a vertical distribution moment ratio 〈∣ z ∣ ^3 〉/〈∣ z ∣ ^2 〉 = 97 ± 15 pc, which corresponds to a Gaussian scale height σ _z = 61 ± 9 pc, an exponential scale height λ _z = 32 ± 5 pc, and a rectangular half-width W _z _,1/2 = 129 ± 20 pc. Consistent with recent simulations, the vertical scale height of cold H i clouds appears to remain constant between the inner Galaxy and the Galactocentric distance of the solar neighborhood. Local fluctuations might explain the large-scale height observed at the same Galactocentric distance on the far side of the Galaxy.

A detailed measurement is made of the metallicity distributions, kinematics, and dynamics of the thin and thick disks across a large disk volume (5.0 ≤ R ≤ 15.0 kpc and ∣ Z ∣ ≤ 3.0 kpc) by using the LAMOST–APOGEE red clump stars. The metallicity distribution results show that the radial metallicity gradient Δ[Fe/H]/Δ R of the thin disk weakens with ∣ Z ∣ from −0.06 dex kpc ^−1 at around ∣ Z ∣ < 0.25 kpc to −0.02 dex kpc ^−1 at around ∣ Z ∣ > 2.75 kpc, while the thick disk displays a global weak positive Δ[Fe/H]/Δ R that is generally weaker than 0.01 dex kpc ^−1 . The vertical metallicity gradient Δ[Fe/H]/Δ∣ Z ∣ steadily weakened from −0.36 dex kpc ^−1 at R ∼ 5.5 kpc to −0.05 dex kpc ^−1 at around R > 11.5 kpc for the thin disk, while the thick disk presents an almost constant value (nearly −0.06∼−0.08 dex kpc ^−1 ) for all the R bins. These results indicate the contribution of the radial migration to the disk evolution, and the obvious north–south asymmetry in [Fe/H] may be linked to disk warp and/or disk perturbation events. The oscillations in the corrected Δ[Fe/H]/Δ∣ Z ∣ with R likely arise from the resonances with the Galactic bar. Our detailed measurements of Δ V _ϕ /Δ[Fe/H] indicate an inside-out and upside-down star formation scenario for the thick disk. The results of eccentricity distributions and [ α /Fe]–velocity dispersion relations are likely to suggest that thick-disk stars require an obvious contribution from other heating mechanisms, such as mergers and accretion, or are born in the chaotic mergers of gas-rich systems and/or the turbulent interstellar medium.

Classical Cepheids (CCs) are excellent tracers for understanding the structure of the Milky Way disk. The latest Gaia Data Release 3 provides a large number of line-of-sight velocity information for Galactic CCs, offering an opportunity for studying the kinematics of the Milky Way. We determine the 3D velocities of 2057 CCs relative to the Galactic center. From the projections of the 3D velocities onto the X – Y plane of the Galactic disk, we find that the V _R and V _ϕ velocities of the northern and southern warps (directions with highest amplitude) are different. This phenomenon may be related to warp precession or asymmetry in warp structure. By investigating the kinematic warp model, we find that the vertical velocities of CCs are more suitable for constraining the warp precession rate than the line-of-node angle. Our results suggest that CCs at 12–14 kpc are the best sample for determining the Galactic warp precession rate. Based on the spatial structure parameters of Cepheid warp from Chen et al., we determine a warp precession rate of ω = 4.9 ± 1.6 km s ^−1 kpc ^−1 at 13 kpc, which supports a low precession rate in the warp model. In the future, more kinematic information on CCs will help to constrain the structure and evolution of the Milky Way better.

As the number of known Galactic structures mounts thanks to the Gaia Space Telescope, it is now pertinent to study methods for disentangling structures occupying the same regions of the Milky Way. Indeed, understanding the precise form of each individual structure and the interactions between structures may aid in understanding their origins and chronology. Moreover, accounting for known structures allows one to probe still finer Galactic structure. In order to demonstrate this, we have developed an odd low-pass filter (OLPF), which removes smaller, odd-parity structures like the vertical waves, and use the filtered data to examine the location of the Galaxy’s midplane. We find that the radial wave identified by Xu et al. continues inward to at least the Sun’s location, with an amplitude that decreases toward the inner, denser parts of the disk, consistent with a simple, qualitative simulation. Additionally, we employ the OLPF results to determine the solar offset, z _⊙ , with smaller structures filtered out. We find that z _⊙ = 34.2 ± 0.3 pc.

Recent results from spectroscopic and astrometric surveys of nearby stars suggest that the stellar disk of our Milky Way (MW) was formed quite early, within the first few billion years of its evolution. Chemokinematic signatures of disk formation in cosmological zoom-in simulations appear to be in tension with these data, implying that MW-like disk formation is delayed in simulations. We investigate the formation of galactic disks using a representative sample of MW-like galaxies from the cosmological volume simulation TNG50. We find that on average MW-mass disks indeed form later than the local data suggest. However, their formation time and metallicity exhibit a substantial scatter, such that ∼10% of MW-mass galaxies form disks early, similar to the MW. Thus, although the MW is unusual, it is consistent with the overall population of MW-mass disk galaxies. The direct MW analogs assemble most of their mass early, ≳10 Gyr ago, and are not affected by destructive mergers after that. In addition, these galaxies form their disks during the early enrichment stage when the interstellar medium metallicity increases rapidly, with only ∼25% of early-forming disks being as metal-poor as the MW was at the onset of disk formation, [Fe/H] ≈ −1.0. In contrast, most MW-mass galaxies either form disks from already enriched material or experience late destructive mergers that reset the signatures of galactic disk formation to later times and higher metallicities. Finally, we also show that earlier disk formation leads to more dominant rotationally supported stellar disks at redshift zero.

We present a comprehensive study of the kinematic properties of the different Galactic disk populations, as defined by the chemical abundance ratios and stellar ages, across a large disk volume (4.5 ≤ R ≤ 15.0 kpc and ∣ Z ∣ ≤ 3.0 kpc), by using the LAMOST-Gaia red clump sample stars. We determine the median velocities for various spatial and population bins, finding large-scale bulk motions; for example, the wave-like behavior in radial velocity, the north–south discrepancy in azimuthal velocity and the warp signal in vertical velocity, and the amplitudes and spatial dependences of these bulk motions show significant variations for different mono-age and mono-abundance populations. The global spatial behaviors of the velocity dispersions clearly show a signal of spiral arms and a signal of the disk perturbation event within 4 Gyr, as well as disk flaring in the outer region (i.e., R ≥ 12 kpc), mostly for young or alpha-poor stellar populations. Our detailed measurements of age/[ α /Fe]-velocity dispersion relations for different disk volumes indicate that young/ α -poor populations are likely to originate from dynamic heating by both giant molecular clouds and spiral arms, while old/ α -enhanced populations require an obvious contribution from other heating mechanisms, such as merger and accretion, or are born in the chaotic mergers of gas-rich systems and/or turbulent interstellar medium.

The orbit distribution of young stars in the Galactic disk is highly structured, from well-defined clusters to streams of stars that may be widely dispersed across the sky, but are compact in orbital action-angle space. The age distribution of such groups can constrain the timescales over which conatal groups of stars disperse into the “field.” Gaia data have proven powerful in identifying such groups in action-angle space, but the resulting member samples are often too small and have too narrow a color–magnitude diagram (CMD) coverage to allow robust age determinations. Here, we develop and illustrate a new approach that can estimate robust stellar population ages for such groups of stars. This first entails projecting the predetermined action-angle distribution into the 5D space of positions, parallaxes, and proper motions, where much larger samples of likely members can be identified over a much wider range of the CMD. It then entails isochrone fitting that accounts for: (a) widely varying distances and reddenings; (b) outliers and binaries; (c) sparsely populated main-sequence turnoffs, by incorporating the age information of the low-mass main sequence; and (d) the possible presence of an intrinsic age spread in the stellar population. When we apply this approach to 92 nearby stellar groups identified in 6D orbit space, we find that they are predominantly young (≲1 Gyr), mono-age populations. Many groups are established (known) localized clusters with possible tidal tails, while others tend to be widely dispersed and manifestly unbound. This new age-dating tool offers a stringent approach to understanding on which orbits stars form in the solar neighborhood and how quickly they disperse into the field.

Galactic gamma-ray diffuse emission (GDE) is emitted by cosmic rays (CRs), ultra-relativistic protons, and electrons, interacting with gas and electromagnetic radiation fields in the interstellar medium. Here we present the analysis of teraelectronvolt diffuse emission from a region of the Galactic plane over the range in longitude of l ∈ [43°, 73°], using data collected with the High Altitude Water Cherenkov (HAWC) detector. Spectral, longitudinal, and latitudinal distributions of the teraelectronvolt diffuse emission are shown. The radiation spectrum is compatible with the spectrum of the emission arising from a CR population with an index similar to that of the observed CRs. When comparing with the DRAGON base model , the HAWC GDE flux is higher by about a factor of 2. Unresolved sources such as pulsar wind nebulae and teraelectronvolt halos could explain the excess emission. Finally, deviations of the Galactic CR flux from the locally measured CR flux may additionally explain the difference between the predicted and measured diffuse fluxes.

In the epoch of deep photometric surveys, a large number of substructures—e.g., overdensities and streams—have been identified. With the help of astrometry and spectroscopy, the community has revealed a complex picture of our Milky Way (MW) after investigating their origins. The off-plane substructures the Anticenter Stream (ACS) and Monoceros Ring (MNC), once considered as dissolving dwarf galaxies, were later found to share similar kinematics and metallicity with the Galactic outer thin disk. In this work, we aim to chemically tag ACS and MNC with high-accuracy abundances from the APOGEE survey. By extrapolating chemical abundance trends in the outer thin-disk region (10 < R _GC < 18 kpc, 0 < ∣ Z _GC ∣ < 3 kpc), we found that ACS and MNC stars show consistent chemical abundances as the extrapolating values for 12 elements, including C, N, O, Mg, Al, Si, K, Ca, Cr, Mn, Co, and Ni. The similar chemical patterns indicate that ACS and MNC have a similar star formation history as the MW outer thin disk, while we also excluded their dwarf galaxy association, as they are distinctive in multiple chemical spaces. The ages of ACS and MNC stars are consistent with the time of the first Sagittarius dSph passage, indicating their possible connection.

The vertical phase-space spiral (snail) is a direct sign of disequilibrium of the Milky Way’s disk. Nevertheless, the wrapping of the phase snail contains information on the vertical potential. We propose a novel method to measure the vertical potential utilizing the intersections between the snail and z / V _z axes, for which we know the maximum vertical heights ( Z _max ) or velocities ( ${V}_{z,\max }$ ). Using a refined linear interpolation method, we directly obtain $({Z}_{\max },\tfrac{1}{2}{V}_{z,\max }^{2})$ for these snail intersections to constrain the vertical potential profile empirically. Our method is model-independent, since no assumptions about the snail shape or the vertical potential have been made. Although the snail binned by the guiding center radius ( R _g ) is more prominent, it traces a vertical potential shallower than that of the snail binned by the same galactocentric radius ( R ). We apply an empirical method to correct this difference. We measure the snail intersections in several R _g bins within 7.5 < R _g < 11.0 kpc for Gaia DR3 and apply the interpolation method to deduce the potential values at several vertical heights. The potential at the snail intersections, as well as the following mass modeling, is consistent with the popular Milky Way potentials in the literature. For the R _g -binned phase snail in the solar neighborhood, the mass modeling indicates a local dark matter density of ρ _dm = 0.0150 ± 0.0031 M _⊙ pc ^−3 , consistent with previous works. Our method could be applied to larger radial ranges in future works to provide independent and stronger constraints on the Milky Way’s potential.

Chemical abundances of Milky Way disk stars are empirical tracers of its enrichment history. However, they capture joint-information that is valuable to disentangle. In this work, we quantify how individual abundances evolve across the present-day Galactic radius, at fixed supernovae contribution ([Fe/H], [Mg/Fe]). We use 18,135 Apache Point Observatory Galactic Evolution Experiment Data Release 17 red clump stars and 7943 GALactic Archaeology with HERMES Data Release 3 main-sequence stars to compare the abundance distributions conditioned on ([Fe/H], [Mg/Fe]) across 3–13 kpc and 6.5–9.5 kpc, respectively. We examine 15 elements: C, N, Al, K (light), O, Si, S, Ca, (α), Mn, Ni, Cr, Cu, (iron-peak) Ce, Ba (s -process) and Eu (r -process). We find that the conditional neutron-capture and light elements most significantly trace variations in the disk's enrichment history, with absolute conditional radial gradients ≤0.03 dex kpc−1. The other elements studied have absolute conditional gradients ≲0.01 dex kpc−1. We uncover structured conditional abundance variations with [Fe/H] for the low- α, but not the high- α , sequence. The average scatter between the mean conditional abundances at different radii is σ intrinsic ≈ 0.02 dex (Ce, Eu, Ba σ intrinsic > 0.05 dex). These results serve as a measure of the magnitude via which different elements trace Galactic radial enrichment history once fiducial supernovae correlations are accounted for. Furthermore, we uncover subtle systematic variations in moments of the conditional abundance distributions and bimodal differences in [Al/Fe]. These suggest a nonuniform enrichment of each chemical cell, and will presumably constrain chemical evolution models of the Galaxy. [ABSTRACT FROM AUTHOR]

We present a Galactic disk vertical velocity analysis using OB type stars (OB stars), red clump (RC) stars, and main-sequence turnoff (MSTO) stars with different average age populations crossmatched with LAMOST DR5 and Gaia DR3. We show that the vertical velocities of the three populations clearly vary with the Galactocentric distance (R) and the younger stellar population has a stronger increasing trend in general. The bending and breathing modes indicated by the vertical motions are dependent on the populations and vary with spatial locations. These vertical motions may be due to the Galactic warp, or minor mergers, or nonequilibrium of the disk. Assuming the warp is the dominant component, we find that the amplitude of the warp (γ, Z ω ) of OB stars (younger population) is larger than that of RC stars (medium population) and the latter is also larger than that for MSTO stars (older population), which is in agreement with other independent analyses of stellar density distribution, and supports that the warp is a long-lived, nonsteady structure and is time evolving. This conclusion is robust whether the line of nodes ϕ w is fixed or is a free parameter (with ϕ w being around 3°−8.°5 as the best fit). Furthermore, we find that the warp is lopsided with asymmetries along the azimuthal angle (ϕ). [ABSTRACT FROM AUTHOR]

Detailed analysis of kinematics of the Milky Way disk in the solar neighborhood based on the GAIA DR3 catalog reveals the existence of peculiarities in the stellar velocity distribution perpendicular to the galactic plane. We study the influence of resonances—the outer Lindblad resonance and the outer vertical Lindblad resonance—of a rotating bar with stellar oscillations perpendicular to the plane of the disk, and their role in shaping the spatial and the velocity distributions of stars. We find that the Z and VZ distributions of stars with respect to LZ are affected by the outer Lindblad resonance. The existence of bar resonance with stellar oscillations perpendicular to the plane of the disk is demonstrated for a long (large semi-axis 5 kpc) and fast rotating bar with Ωb=60.0kms−1kpc−1. We show also that, in the model with the long and fast rotating bar, some stars in the 2:1 OLR region deviate far from their original places, entering the bar region. A combination of resonance excitation of stellar motions at the 2:1 OLR region together with strong interaction of the stars with the bar potential leads to the formation of the group of ‘escapees’, i.e., stars that deviate in R and Z—directions at large distances from the resonance region. Simulations, however, do not demonstrate any noticeable effect on VZ-distribution of stars in the solar neighborhood.

Using a complete sample of about 330,000 dwarf stars, well measured by Gaia DR3, limited to the galactic north and south solid angles |b|<75° and up to a vertical distance of 2 kpc, we analyze the vertical structure of the Milky Way stellar disks, based on projected tangential velocities. From selected subsamples dominated by their corresponding population, we obtain the thin and thick disk scale heights as hZ=279.76±12.49 pc and HZ=797.23±12.34 pc, respectively. Then from the simultaneous fitting of the sum of two populations over the whole sample, assuming these scale heights, we estimate the thick-to-thin disk number density ratio at the galactic plane to be ρT/ρt=0.750±0.049, which is consistent with a previous result by the authors: in the galactic plane there is a significant number of thick disk stars, possibly as many as thin disk ones, which also points to the existence of more thick disk stars than generally thought. The overall fit does not closely follow the data for |Z|>700 pc and points to the presence of more stars beyond the thin disk that cannot be accounted for by the two-disk model.

A relationship between chromospheric activity and age is calibrated for FGK dwarf stars using GALEX F U V magnitudes and Gaia (G B P − G) colors. Such a calibration between GALEX F U V magnitudes and stellar age has utility in population studies of dwarfs for further understanding of the chemical evolution of the Milky Way. As an illustration of one such application we have investigated a population of Sun-like, solar neighborhood stars for their metallicities and velocity dispersions; a cross-matched sample of FGK type dwarf stars from Casagrande et al. (2011) with the Gaia and GALEX catalogs. Using calibrated relationships between F U V magnitudes and age, we determined a chromospheric activity indicator, Q , and stellar age, τ , for each dwarf. We constructed age-velocity (AVR) and age-metallicity (AMR) relations with empirically-determined F U V ages. Power law fits to AVR plots are consistent with heating mechanism models within the literature. We further demonstrate that perigalactic distance and eccentricity versus F U V -age plots are consistent with an "inside out" formation history model. [ABSTRACT FROM AUTHOR]

Stellar ages are key for determining the formation history of the Milky Way, but are difficult to measure precisely. Furthermore, methods that use chemical abundances to infer ages may entangle the intrinsic evolution of stars with the chemodynamical evolution of the Galaxy. In this paper, we present a framework for making probabilistic predictions of stellar ages, and then quantify the contribution of both stellar evolution and Galactic chemical evolution to those predictions using SHapley Additive exPlanations. We apply this interpretable prediction framework to both a simulated Milky Way sample containing stars in a variety of evolutionary stages and an APOGEE-mocked sample of red clump stars. We find that in the former case, stellar evolution is the dominant driver for age estimates, while in the latter case, the more restricted evolutionary information causes the model to proxy ages through the chemical evolution model. We show that as a result of the use of nonintrinsic Galactic chemical information, trends estimated with the predicted ages, such as the age–metallicity relation, can deviate from the truth.

We explore the kinematic and chemical properties of the Monoceros stellar overdensity by combining data from the Two Micron All Sky Survey, Wide-field Infrared Survey Explorer, APOGEE, and Gaia. Monoceros is a structure located toward the Galactic anticenter and close to the disk. We have identified that its stars have azimuthal velocity in the range of 200 < v _ϕ (km s ^−1 ) < 250. Combining their kinematics and spatial distribution, we designed a new method to select stars from this overdensity. This method allows us to easily identify the structure in both hemispheres and estimate their distances. Our analysis was supported by comparison with simulated data from the entire sky generated by the Galaxia code. Furthermore, we characterized, for the first time, the Monoceros overdensity in several chemical abundance spaces. Our results confirm its similarity to stars found in the thin disk of the Galaxy and suggest an in situ formation. Furthermore, we demonstrate that the southern and northern regions of Monoceros exhibit indistinguishable chemical compositions.

With the unprecedented increase in the number of known star clusters, quick and modern tools are needed for their analysis. In this work, we develop an artificial neural network (ANN) trained on synthetic clusters to estimate the age, metallicity, extinction, and distance of Gaia open clusters. We implement a novel technique to extract features from the color–magnitude diagram of clusters by means of the QuadTree tool, and we adopt a multiband approach. We obtain reliable parameters for ∼5400 clusters. We demonstrate the effectiveness of our methodology in accurately determining crucial parameters of Gaia open clusters by performing a comprehensive scientific validation. In particular, with our analysis we have been able to reproduce the Galactic metallicity gradient as it is observed by high-resolution spectroscopic surveys. This demonstrates that our method reliably extracts information on metallicity from color–magnitude diagrams (CMDs) of stellar clusters. For the sample of clusters studied, we find an intriguing systematic older age compared to previous analyses present in the literature. This work introduces a novel approach to feature extraction using a QuadTree algorithm, effectively tracing sequences in CMDs despite photometric errors and outliers. The adoption of ANNs, rather than convolutional neural networks, maintains the full positional information and improves performance, while also demonstrating the potential for deriving cluster parameters from simultaneous analysis of multiple photometric bands, beneficial for upcoming telescopes like the Vera Rubin Observatory. The implementation of ANN tools with robust isochrone fit techniques could provide further improvements in the quest for open cluster parameters.

We use the largest catalog of open clusters in the post-Gaia era to provide an observational view of the Galactic disk. By compiling physical parameters such as age, distance, and kinematic information, we investigate the spatial distribution of open clusters and revisit the spiral arms and other asymmetries in the Galactic disk. Using young open clusters as a tracer of spiral arms, we map the spiral structure of the Galaxy and find that most of the clusters start migrating away from the spiral arms in about 10–20 Myr and fill the interarm regions as they age. Using the 3D kinematic information on 371 open star clusters, we derive different individual pattern speeds for spiral arms that closely follow the rotation curve of the Milky Way, hence favoring the transient nature of spiral arms in the Milky Way. The pattern rotation speeds of each spiral arm suggest that the spiral arms have not accelerated in the last 80 Myr. Based on the distribution of open clusters younger than 700 Myr above or below the Galactic plane, we found a solar offset of z _⊙ = 17.0 ± 0.9 pc north of the Galactic plane and estimated the scale height z _h = 91.7 ± 1.9 pc from the Galactic plane.

The recent discovery of a spiral pattern in the vertical kinematic structure in the solar neighborhood provides a prime opportunity to study nonequilibrium dynamics in the Milky Way from local stellar kinematics. Furthermore, results from simulations indicate that even in a limited volume, differences in stellar orbital histories allow us to trace variations in the initial perturbation across large regions of the disk. We present ESCARGOT , a novel algorithm for studying these variations in both simulated and observed data sets. ESCARGOT automatically extracts key quantities from the structure of a given phase spiral, including the time since perturbation and the perturbation mode. We test ESCARGOT on simulated data and show that it is capable of accurately recovering information about the time since the perturbation occurred as well as subtle differences in phase spiral morphology due to stellar locations in the disk at the time of perturbation. We apply ESCARGOT to kinematic data from data release 3 of the Gaia mission in bins of guiding radius. We show that similar structural differences in morphology occur in the Gaia phase spirals as a function of stellar orbital history. These results indicate that the phase spirals are the product of a complex dynamical response in the disk with large-scale coupling between different regions of phase space.

Varying oxygen abundance could impact modeling-inferred ages. This work aims to estimate the ages of dwarfs considering observed oxygen abundance. To characterize 67,503 LAMOST and 4006 GALAH FGK-type dwarf stars, we construct a grid of stellar models, which take into account oxygen abundance as an independent model input. Compared with ages determined with commonly used α -enhanced models, we find a difference of ∼9% on average when the observed oxygen abundance is considered. The age differences between the two types of models are correlated to [Fe/H] and [O/ α ], and they are relatively significant on stars with [Fe/H] ≲ −0.6 dex. Generally, varying 0.2 dex in [O/ α ] will alter the age estimates of metal-rich (−0.2 < [Fe/H] < 0.2) stars by ∼10% and relatively metal-poor (−1 < [Fe/H] < −0.2) stars by ∼15%. Of the low-O stars with [Fe/H] < 0.1 dex and [O/ α ] ∼−0.2 dex, many have fractional age differences of ≥ 10% and even reach up to 27%. The fractional age difference of high-O stars with [O/ α ] ∼0.4 dex reaches up to −33% to −42% at [Fe/H] ≲ −0.6 dex. We also analyze the chemical properties of these stars. We find a decreasing trend of [Fe/H] with ages from 7.5–9 Gyr to 5–6.5 Gyr for the stars from the LAMOST and GALAH. The [O/Fe] of these stars increases with decreasing age from 7.5–9 Gyr to 3–4 Gyr, indicating that the younger population is more O rich.

Metallicity gradients refer to the sloped radial profiles of the metallicities of gas and stars and are commonly seen in disk galaxies. A well-defined metallicity gradient of the Galactic disk is observed particularly well with classical Cepheids, which are good stellar tracers thanks to their period–luminosity relation, allowing precise distance estimation and other advantages. However, the measurement of the inner-disk gradient has been impeded by the incompleteness of previous samples of Cepheids and the limitations of optical spectroscopy in observing highly reddened objects. Here we report the metallicities of 16 Cepheids measured with high-resolution spectra in the near-infrared YJ bands. These Cepheids are located at 3–5.6 kpc in Galactocentric distance, R _GC , and reveal the metallicity gradient in this range for the first time. Their metallicities are mostly between 0.1 and 0.3 dex in [Fe/H] and more or less follow the extrapolation of the metallicity gradient found in the outer part, R _GC > 6.5 kpc. The gradient in the inner disk may be shallower or even flat, but the small sample does not allow the determination of the slope precisely. More extensive spectroscopic observations would also be necessary for studying minor populations, if any, with higher or lower metallicities that were reported in previous literature. In addition, the 3D velocities of our inner-disk Cepheids show a kinematic pattern that indicates noncircular orbits caused by the Galactic bar, which is consistent with the patterns reported in recent studies on high-mass star-forming regions and red giant branch stars.

We present new maps of the Milky Way disk showing the distribution of metallicity ([Fe/H]), α -element abundances ([Mg/Fe]), and stellar age, using a sample of 66,496 red giant stars from the final data release (DR17) of the Apache Point Observatory Galactic Evolution Experiment survey. We measure radial and vertical gradients, quantify the distribution functions for age and metallicity, and explore chemical clock relations across the Milky Way for the low- α disk, high- α disk, and total population independently. The low- α disk exhibits a negative radial metallicity gradient of −0.06 ± 0.001 dex kpc ^−1 , which flattens with distance from the midplane. The high- α disk shows a flat radial gradient in metallicity and age across nearly all locations of the disk. The age and metallicity distribution functions shift from negatively skewed in the inner Galaxy to positively skewed at large radius. Significant bimodality in the [Mg/Fe]–[Fe/H] plane and in the [Mg/Fe]–age relation persist across the entire disk. The age estimates have typical uncertainties of ∼0.15 in log(age) and may be subject to additional systematic errors, which impose limitations on conclusions drawn from this sample. Nevertheless, these results act as critical constraints on galactic evolution models, constraining which physical processes played a dominant role in the formation of the Milky Way disk. We discuss how radial migration predicts many of the observed trends near the solar neighborhood and in the outer disk, but an additional more dramatic evolution history, such as the multi-infall model or a merger event, is needed to explain the chemical and age bimodality elsewhere in the Galaxy.

We have for the first time identified the early stellar disk in the Milky Way by using a combination of elemental abundances and kinematics. Using data from APOGEE DR17 and Gaia we select stars in the Mg–Mn–Al–Fe plane with elemental abundances indicative of an accreted origin and find stars with both halo-like and disk-like kinematics. The stars with halo-like kinematics lie along a lower sequence in [Mg/Fe], while the stars with disk-like kinematics lie along a higher sequence. Combined with astroseismic observations, we determine that the stars with halo-like kinematics are old, 9–11 Gyr, and that the more evolved stellar disk is about 1–2 Gyr younger. We show that the in situ fraction of stars on deeply bound orbits is not small, in fact the inner Galaxy likely harbors a genuine in situ population together with an accreted one. In addition, we show that the selection of the Gaia–Sausage–Enceladus in the E _n − L _z plane is not very robust. In fact, radically different selection criteria give almost identical elemental abundance signatures for the accreted stars.

The tilt of the velocity ellipsoid is a helpful tracer of the gravitational potential of the Milky Way. In this paper, we use nearly 140,000 red clump (RC) stars selected from LAMOST and Gaia to make a detailed analysis of the tilt of the velocity ellipsoid for various populations, as defined by stellar ages and chemical information, within 4.5 ≤ R ≤ 15.0 kpc and ∣ Z ∣ ≤ 3.0 kpc, respectively. The tilt angle of the velocity ellipsoid of the RC sample stars is accurately described as α = α _0 $\arctan $ ( Z / R ) with α _0 = (0.68 ± 0.05). This indicates the alignment of the velocity ellipsoid is between cylindrical and spherical, implying that any deviation from the spherical alignment of the velocity ellipsoid may be caused by the gravitational potential of the baryonic disk. The results of various populations suggest that α _0 displays an age and population dependence, with values α _0 = (0.72 ± 0.08) and α _0 = (0.64 ± 0.07) for the thin and thick disks, respectively, and α _0 displays a decreasing trend as age (and [ α /Fe]) increases, meaning that the velocity ellipsoids of the kinematically relaxed stars are mainly dominated by the gravitational potential of the baryonic disk. We determine the α _0 –R relation for various populations, finding that α _0 displays oscillations with R for the different populations. The oscillations in α _0 appear in both kinematically hot and cold populations, indicating that resonances with the Galactic bar are the most likely origin for these oscillations.

We investigate the angular momentum of mono-abundance populations (MAPs) of the Milky Way thick disk by using a sample of 26,076 giant stars taken from APOGEE Data Release (DR) 17 and Gaia early DR3. The vertical and perpendicular angular momentum components, L _Z and L _P , of the MAPs in narrow bins have significant variations across the [ α /M]–[M/H] plane. L _Z and L _P systematically change with [M/H] and [ α /M] and can be alternatively quantified by the chemical gradients: d [M/H]/ dL _Z = 1.2 × 10 ^−3 dex kpc ^−1 km ^−1 s, d [M/H]/ dL _P = −5.0 × 10 ^−3 dec kpc ^−1 km ^−1 s, and d [ α /M]/ dL _Z = −3.0 × 10 ^−4 dex kpc ^−1 km ^−1 s, d [ α /M]/ dL _P = 1.2 × 10 ^−3 dec kpc ^−1 km ^−1 s. These correlations can also be explained as the chemical dependence of the spatial distribution shape of the MAPs. We also exhibit the corresponding age dependence of the angular momentum components. Under the assumption that the guiding radius ( R _g ) is proportional to L _Z , this provides direct observational evidence of the inside-out structure formation scenario of the thick disk, with dR _g / d Age = −1.9 kpc Gyr ^−1 . The progressive changes in the disk thickness can be explained by the upside-down formation or/and the consequent kinematical heating.