Search NASA⌕ Search

SEARCH · Search NASA

Results for “Milky Way disk”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Slant, Fan, and Narrow: The Response of Stellar Streams to a Tilting Galactic Disk

Stellar streams are sensitive tracers of the gravitational potential, which is typically assumed to be static in the inner Galaxy. However, massive mergers like Gaia-Sausage-Enceladus can impart torques on the stellar disk of the Milky Way that result in the disk tilting at rates of up to 10°–20° Gyr –1 . Here, we demonstrate the effects of disk tilting on the morphology and kinematics of stellar streams. Through a series of numerical experiments, we find that streams with nearby apocenters (r apo ≲ 20 kpc) are sensitive to disk tilting, with the primary effect being changes to the stream's on-sky track and width. Interestingly, disk tilting can produce both more diffuse streams and more narrow streams, depending on the orbital inclination of the progenitor and the direction in which the disk is tilting. Our model of Pal 5's tidal tails for a tilting rate of 15° Gyr –1 is in excellent agreement with the observed stream's track and width, and reproduces the extreme narrowing of the trailing tail. We also find that failure to account for a tilting disk can bias constraints on shape parameters of the Milky Way's local dark matter distribution at the level of 5%–10%, with the direction of the bias changing for different streams. Disk tilting could therefore explain discrepancies in the Milky Way's dark matter halo shape inferred using different streams.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

The AGORA High-resolution Galaxy Simulations Comparison Project. VIII. Disk Formation and Evolution of Simulated Milky Way Mass Galaxy Progenitors at 1 < z < 5

We investigate how differences in the stellar feedback produce disks with different morphologies in Milky Way–like progenitors over 1 ≤ z ≤ 5, using eight state-of-the-art cosmological hydrodynamics simulation codes in the AGORA project. In three of the participating codes, a distinct, rotation-dominated inner core emerges with a formation timescale of ≲300 Myr, largely driven by a major merger event, while two other codes exhibit similar signs of wet compaction—gaseous shrinkage into a compact starburst phase—at earlier epochs. The remaining three codes show only weak evidence of wet compaction. Consequently, we divide the simulated galaxies into two groups: those with strong compaction signatures and those with weaker ones. Galaxies in these two groups differ in size, stellar age gradients, and disk-to-total mass ratios. Specifically, codes with strong wet compaction build their outer disks in an inside-out fashion, leading to negative age gradients, whereas codes with weaker compaction feature flat or positive age gradients caused primarily by outward stellar migration. Although the stellar half-mass radii of these two groups diverge at z ∼ 3, the inclusion of dust extinction brings their sizes and shapes in mock observations closer to each other and to observed galaxies. We attribute the observed morphological differences primarily to variations in the stellar feedback implementations—such as delayed cooling timescales, and feedback strengths—that regulate both the onset and duration of compaction. Overall, our results suggest that disk assembly at high redshifts is highly sensitive to the details of the stellar feedback prescriptions in simulations.

Jung, Minyong [Seoul National Univ. (Korea, Republ↗
DOE OSTI2023

Two-point Correlation Function Studies for the Milky Way: Discovery of Spatial Clustering from Disk Excitations and Substructure

We introduce a two-particle correlation function (2PCF) for the Milky Way, constructed to probe spatial correlations in the orthogonal directions of the stellar disk in the Galactic cylindrical coordinates of R, ϕ, and z. We use this new tool to probe the structure and dynamics of the Galaxy using the carefully selected set of solar neighborhood stars (d ≲ 3 kpc) from Gaia Data Release 2 that we previously employed for studies of axial symmetry breaking in stellar number counts. We make additional, extensive tests, comparing to reference numerical simulations, to ensure our control over possibly confounding systematic effects. Supposing either axial or north–south symmetry, we divide this data set into two nominally symmetric sectors and construct the 2PCF, in the manner of the Landy–Szalay estimator, from the Gaia data. In so doing, working well away from the midplane region in which the spiral arms appear, we have discovered distinct symmetry-breaking patterns in the 2PCF in its orthogonal directions, thus establishing the existence of correlations in stellar number counts alone at subkiloparsec length scales for the very first time. In particular, we observe extensive wavelike structures of amplitude greatly in excess of what we would estimate if the system were in a steady state. We study the variations in these patterns across the Galactic disk, and with increasing ∣z∣, and we show how our results complement other observations of non-steady-state effects near the Sun, such as vertical asymmetries in stellar number counts and the Gaia snail.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A Study of the Properties and Dynamics of the Disk of Satellites in a Milky-Way-like Galaxy System

The dynamics of the satellite systems of Milky-Way-like galaxies offer a useful means by which to study the galaxy formation process in the cosmological context. It has been suggested that the currently observed anisotropic distribution of the satellites in such galaxy systems is inconsistent with the concordance ΛCDM cosmology model on the galactic scale if the observed satellites are random samples of the dark matter (DM) sub-halos that are nearly isotropically distributed around the central galaxy. In this study, we present original high-resolution zoom-in studies of central galaxies and satellite systems based upon initial conditions for the DM distribution from the Aquarius simulations but with substantial high-resolution baryon physics added. We find that the galaxy most like the Milky Way in this study does indeed contain a disk of satellites (DOS). Although one galaxy DOS system does not answer the question of how common such disks are, it does allow the opportunity to explore the properties and dynamics of the DOS system. Our investigation centers on the spatial arrangement (distances, angles, etc.) of satellites in this Milky-Way-like galaxy system with a specific emphasis on identifying and analyzing the disk-like structure along with its dynamical and morphological properties. Among the conclusions from this study, we find that the satellites and DM sub-halos in the galaxy simulations are anisotropically distributed. The dynamical properties of the satellites, however, indicate that the direction of the angular momentum vector of the whole satellite system is different from the normal direction of the fitted DOS and from the normal direction of the velocity dispersion of the system. Hence, the fitted DOS appears to be comprised of infalling sub-halos and is not a rotationally supported system.

Zhao, Xinghai↗
DOE OSTI2026

Vertical Structure of Local Disk Galaxies Revealed by DESI Imaging Data

The vertical structure of galactic disks is an important probe of the disk assembly history. Here we investigate the vertical structure of a sample of 79 local disk galaxies within 50 Mpc using data from the Dark Energy Spectroscopic Instrument Legacy Imaging Surveys. Vertical luminosity profiles as a function of radius in the g, r, and z bands are extracted and fitted using a single-component sech2 model to determine scale height and its radial variation (i.e., flaring). Our measurements indicate local galactic disks are overall thin with negligible flaring. The median scale heights at 1 R e are 0.21, 0.22, and 0.22 kpc in the g, r, and z bands, respectively, while the median radial gradients of scale height are −0.006, 0.003, and 0.001 in these bands. These values are consistent with that of the geometric thin disk of the Milky Way represented by metal-rich, low-[α/Fe] populations, confirming the weak flaring of the geometric thin disk. A clear positive correlation of scale height with stellar mass is observed down to the low-mass end of 10 7 M ⊙ . These results provide a homogeneous benchmark for the vertical structure of nearby disk galaxies and establish the cosmological representativeness of the Milky Way thin disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

EDEN: Exploring Disks Embedded in N-body Simulations of Milky Way–mass Halos from Symphony

We investigate the impact of galactic disks on the tidal stripping of cold dark matter subhalos within Milky Way (MW)-mass halos (M vir ∼ 10 12 M ⊙ ) using a new simulation suite, EDEN. By re-simulating 45 MW-mass zoomed-in halos from the N-body Symphony compilation with embedded disk potentials, which evolve according to star formation histories predicted by the UniverseMachine model, we self-consistently tie disk growth to halo accretion rate and significantly expand the range of disk masses and formation histories studied. We use the particle-tracking-based subhalo finder Symfind to enhance the robustness of subhalo tracking. We find that disks near the median disk-to-halo mass ratio of our sample (M *,Disk /M vir,host = 2%) reduce subhalo peak mass functions within 100 kpc by about 10% for peak masses above 10 8 M ⊙ . Heavier, MW/M31-like disks (M *,Disk /M vir,host ≳ 5%) lead to a reduction of more than 40%. Subhalo abundance suppression is more pronounced near halo centers, with particularly enhanced stripping for subhalos accreted over 8 Gyr ago on orbits with pericenters <100 kpc. Suppression is further amplified when disk mass is increased within fixed halo and disk assembly histories. In all cases, the suppression we measure should be interpreted as stripping below the mass resolution limit rather than complete subhalo disruption. This study reshapes our understanding of the MW’s impact on its satellites, suggesting it strips subhalos more efficiently than typical MW-mass galaxies due to its larger disk-to-halo mass ratio and earlier disk formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Ground Vibrational State SiO Emission in the VLA BAaDE Survey

Using a subsample of the Bulge Asymmetries and Dynamical Evolution (BAaDE) survey of stellar SiO masers, we explore the prevalence and characteristics of {sup 28}SiO J = 1 − 0, v = 0 emission. We identify 90 detections of maser, thermal, or composite {sup 28}SiO J = 1 − 0, v = 0 emission out of approximately 13,000 candidate spectra from the National Science Foundation’s Karl G. Jansky Very Large Array (VLA). We find that the detected sources are likely asymptotic giant branch stars belonging to a bright, foreground Milky Way stellar disk population. For the 32 sources showing thermal components, we extract values for outflow velocity by fitting thermal line profiles. We find a range of circumstellar envelope expansion velocities, and compare to previously recorded OH and CO expansion velocities. This preliminary survey is already the largest study of stellar ground vibrational state SiO masers to date, and will be expanded to include the entire VLA BAaDE data set when data reduction for the 18,988 target sources is completed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

New determination of the millisecond pulsar gamma-ray luminosity function and implications for the Galactic Center gamma-ray excess

It has been suggested that the Galactic Center gamma-ray excess (GCE) could be produced by a large number of centrally located millisecond pulsars. The fact that no such pulsar population has been detected implies that these sources must be very faint and very numerous. Here, in this study, we use the contents of Fermi’s recently released Third Pulsar Catalog (3PC) to measure the luminosity function of the millisecond pulsars in the Milky Way’s disk. We find that this source population exhibits a luminosity function with a mean γ -ray luminosity of ⟨ L γ ⟩ ∼ 6 × 10 32 erg / s (integrated above 0.1 GeV). If the GCE were generated by millisecond pulsars with the same luminosity function, we find that ∼ 20 such sources from the inner Galaxy population should have already been detected by Fermi and included in the 3PC. Given the lack of such observed sources, we exclude the hypothesis that the GCE is generated by pulsars with the same luminosity function as those in the Galactic disk with a significance of 3.4 σ . We conclude that either less than 39% of the GCE is generated by pulsars, or that the millisecond pulsars in the inner Galaxy are at least 5 times less luminous on average than those found in the Galactic disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Signatures of a Tidally Induced Spiral Arm at the Anticenter of the Milky Way and a Kinematically Extended Anticenter Stream Using DESI Data Release 2

Using the Dark Energy Spectroscopic Instrument (DESI) Milky Way Survey, we examine the six-dimensional space of the anticenter region of the Milky Way stellar disk (150° < Galactic longitude < 220°) using 61,883 main-sequence turnoff stars. We focus on two well-known stellar overdensities in the anticenter: the Monoceros Ring (MRi) and Anticenter Stream (ACS). We find that the MRi overdensity has kinematic signatures consistent with a tidally induced spiral arm, a type of dynamic spiral arm created by an interaction with a satellite galaxy, most likely the Sagittarius dwarf spheroidal galaxy (Sgr). We use the kinematics of the MRi to calculate the two most recent passage times of Sgr, finding 0.25 ± 0.09 Gyr and 1.10 ± 0.23 Gyr from the present day. We validate that the ACS is kinematically decoupled from the MRi because they are moving in opposite radial and vertical directions. We find that the kinematics associated with the ACS extends beyond our defined overdensity. The features we see in the ACS region are likely part of a broader distribution of stars with the same kinematic signature as detected in other places, like the vertical wave in the outer disk and phase spiral.

Lambert, Mika [University of California, Santa Cru↗
DOE OSTI2020

Gravitational Microlensing Event Statistics for the Zwicky Transient Facility

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 ~10 9 stars in the g and r bands and can significantly contribute to the observed microlensing population. In this paper, 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.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Functional Data Analysis for Extracting the Intrinsic Dimensionality of Spectra: Application to Chemical Homogeneity in the Open Cluster M67

High-resolution spectroscopic surveys of the Milky Way have entered the Big Data regime and have opened avenues for solving outstanding questions in Galactic archeology. However, exploiting their full potential is limited by complex systematics, whose characterization has not received much attention in modern spectroscopic analyses. In this work, we present a novel method to disentangle the component of spectral data space intrinsic to the stars from that due to systematics. Using functional principal component analysis on a sample of 18,933 giant spectra from APOGEE, we find that the intrinsic structure above the level of observational uncertainties requires ≈10 functional principal components (FPCs). Our FPCs can reduce the dimensionality of spectra, remove systematics, and impute masked wavelengths, thereby enabling accurate studies of stellar populations. To demonstrate the applicability of our FPCs, we use them to infer stellar parameters and abundances of 28 giants in the open cluster M67. We employ Sequential Neural Likelihood, a simulation-based Bayesian inference method that learns likelihood functions using neural density estimators, to incorporate non-Gaussian effects in spectral likelihoods. By hierarchically combining the inferred abundances, we limit the spread of the following elements in M67: Fe ≲ 0.02 dex; C ≲ 0.03 dex; O, Mg, Si, Ni ≲ 0.04 dex; Ca ≲ 0.05 dex; N, Al ≲ 0.07 dex (at 68% confidence). Our constraints suggest a lack of self-pollution by core-collapse supernovae in M67, which has promising implications for the future of chemical tagging to understand the star formation history and dynamical evolution of the Milky Way.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Quark stars with 2.6 $M_\odot $ in a non-minimal geometry-matter coupling theory of gravity

This work analyses the hydrostatic equilibrium configurations of strange stars in a non-minimal geometry-matter coupling (GMC) theory of gravity. Those stars are made of strange quark matter, whose distribution is governed by the MIT equation of state. The non-minimal GMC theory is described by the following gravitational action: f(R,L) = R/2 + L + σRL, where R represents the curvature scalar, L is the matter Lagrangian density, and σ is the coupling parameter. When considering this theory, the strange stars become larger and more massive. In particular, when σ = 50 km 2 , the theory can achieve the 2.6 M ⊙ , which is suitable for describing the pulsars PSR J2215+5135 and PSR J1614-2230, and the mass of the secondary object in the GW190814 event. The 2.6 M ⊙ is a value hardly achievable in General Relativity, even considering fast rotation effects, and is also compatible with the mass of PSR J0952-0607 (M = 2.35 ± 0.17 M ⊙ ), the heaviest and fastest pulsar in the disk of the Milky Way, recently measured, supporting the possible existence of strange quark matter in its composition. The non-minimal GMC theory can also give feasible results to describe the macroscopical features of strange star candidates.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A Search for Neutrino Point-source Populations in 7 yr of IceCube Data with Neutrino-count Statistics

The presence of a population of point sources in a data set modifies the underlying neutrino-count statistics from the Poisson distribution. This deviation can be exactly quantified using the non-Poissonian template fitting technique, and in this work we present the first application of this approach to the IceCube high-energy neutrino data set. Using this method, we search in 7 yr of IceCube data for point-source populations correlated with the disk of the Milky Way, the Fermi bubbles, the Schlegel, Finkbeiner, and Davis dust map, or with the isotropic extragalactic sky. No evidence for such a population is found in the data using this technique, and in the absence of a signal, we establish constraints on population models with source-count distribution functions that can be described by a power law with a single break. Here, the derived limits can be interpreted in the context of many possible source classes. In order to enhance the flexibility of the results, we publish the full posterior from our analysis, which can be used to establish limits on specific population models that would contribute to the observed IceCube neutrino flux.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Gaia-ESO Survey: Oxygen Abundance in the Galactic Thin and Thick Disks

We analyze the oxygen abundances of a stellar sample representative of the two major Galactic populations: the thin and thick disks. The aim is to investigate the differences between members of the Galactic disks and contribute to the understanding of the origin of oxygen chemical enrichment in the Galaxy. The analysis is based on the [O i] = 6300.30 Å oxygen line in high-resolution spectra (R ∼ 52,500) obtained from the Gaia-ESO public spectroscopic Survey (GES). By comparing the observed spectra with a theoretical data set computed in LTE with the SPECTRUM synthesis and ATLAS12 codes, we derive the oxygen abundances of 516 FGK dwarfs for which we have previously measured carbon abundances. Based on kinematic, chemical, and dynamical considerations, we identify 20 thin and 365 thick disk members. We study the potential trends of both subsamples in terms of their chemistry ([O/H], [O/Fe], [O/Mg], and [C/O] versus [Fe/H] and [Mg/H]), age, and position in the Galaxy. The main results are that (a) [O/H] and [O/Fe] ratios versus [Fe/H] show systematic differences between thin and thick disk stars with an enhanced O abundance of thick disk stars with respect to thin disk members and a monotonic decrement of [O/Fe] with increasing metallicity, even at metal-rich regime; (b) there is a smooth correlation of [O/Mg] with age in both populations, suggesting that this abundance ratio can be a good proxy of stellar ages within the Milky Way; and (c) thin disk members with [Fe/H] ≃ 0 display a [C/O] ratio smaller than the solar value, suggesting a possibly outward migration of the Sun from lower Galactocentric radii.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

The Stellar Velocity Distribution Function in the Milky Way Galaxy

The stellar velocity distribution function in the solar vicinity is reexamined using data from the Sloan Digital Sky Survey Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey’s DR16 and Gaia DR2. By exploiting APOGEE’s ability to chemically discriminate with great reliability the thin-disk, thick-disk, and (accreted) halo populations, we can, for the first time, derive the three-dimensional velocity distribution functions (DFs) for these chemically separated populations. We employ this smaller but more data-rich APOGEE+Gaia sample to build a data-driven model of the local stellar population velocity DFs and use these as basis vectors for assessing the relative density proportions of these populations over the 5 < R < 12 kpc and −1.5 < z < 2.5 kpc range as derived from the larger, more complete (i.e., all-sky, magnitude-limited) Gaia database. We find that 81.9% ± 3.1% of the objects in the selected Gaia data set are thin-disk stars, 16.6% ± 3.2% are thick-disk stars, and 1.5% ± 0.1% belong to the Milky Way stellar halo. We also find the local thick-to-thin-disk density normalization to be ρ {sub T}(R {sub ⊙})/ρ {sub t}(R {sub ⊙}) = 2.1% ± 0.2%, a result consistent with, but determined in a completely different way from, typical star-count/density analyses. Using the same methodology, the local halo-to-disk-density normalization is found to be ρ {sub H}(R {sub ⊙})/(ρ {sub T}(R {sub ⊙}) + ρ {sub t}(R {sub ⊙})) = 1.2% ± 0.6%, a value that may be inflated due to the chemical overlap of halo and metal-weak thick-disk stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Supernova-driven Turbulent Metal Mixing in High-redshift Galactic Disks: Metallicity Fluctuations in the Interstellar Medium and its Imprints on Metal-poor Stars in the Milky Way

The extent to which turbulence mixes gas in the face of recurrent infusions of fresh metals by supernovae (SN) could help provide important constraints on the local star formation conditions. This includes predictions of the metallicity dispersion among metal-poor stars, which suggests that the interstellar medium was not very well mixed at these early times. The purpose of this Letter is to help isolate, via a series of numerical experiments, some of the key processes that regulate turbulent mixing of SN elements in galactic disks. We study the gas interactions in small simulated patches of a galaxy disk with the goal of resolving the small-scale mixing effects of metals at parsec scales, which enables us to measure the turbulent diffusion coefficient in various galaxy environments. By investigating the statistics of variations of α elements in these simulations, we are able to derive constraints not only on the allowed range of intrinsic yield variations in SN explosions but also on the star formation history of the Milky Way. We argue that the observed dispersion of [Mg/Fe] in metal-poor halo stars is compatible with the star-forming conditions expected in dwarf satellites or in an early low-star-forming Milky Way progenitor. In particular, metal variations in stars that have not been phase-mixed can be used to infer the star-forming conditions of disrupted dwarf satellites.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A Blueprint for the Milky Way’s Stellar Populations. IV. A String of Pearls—the Galactic Starburst Sequence

We continue our series of papers on phase-space distributions of stars in the Milky Way based on photometrically derived metallicities and Gaia astrometry, with a focus on the halo-disk interface in the local volume. To exploit various photometric databases, we develop a method of empirically calibrating synthetic stellar spectra based on a comparison with observations of stellar sequences and individual stars in the Sloan Digital Sky Survey, the SkyMapper Sky Survey, and the Pan-STARRS1 surveys, overcoming band-specific corrections employed in our previous work. In addition, photometric zero-point corrections are derived to provide an internally consistent photometric system with a spatially uniform metallicity zero-point. Using our phase-space diagrams, we find a remarkably narrow sequence in the rotational velocity (v $\phi$ ) versus metallicity ([Fe/H]) space for a sample of high proper-motion stars (>25 mas yr -1 ), which runs along Gaia Sausage/Enceladus (GSE) and the Splash substructures and is linked to the disk, spanning nearly 2 dex in [Fe/H]. Notably, a rapid increase of v $\phi$ from a nearly zero net rotation to ~180 km s -1 in a narrow metallicity interval (-0.6 ≲ [Fe/H] ≲ -0.4) suggests that some of these stars emerged quickly on a short gas-depletion timescale. Through measurements of a scale height and length, we argue that these stars are distinct from those heated dynamically by mergers. This chain of high proper-motion stars provides additional support for recent discoveries suggesting that a starburst took place when the young Milky Way encountered the gas-rich GSE progenitor, which eventually led to the settling of metal-enriched gas onto the disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Detection of open cluster rotation fields from Gaia EDR3 proper motions

Context: Most stars from in groups which with time disperse, building the field population of their host galaxy. In the Milky Way, open clusters have been continuously forming in the disk up to the present time, providing it with stars spanning a broad range of ages and masses. Observations of the details of cluster dissolution are, however, scarce. One of the main difficulties is obtaining a detailed characterisation of the internal cluster kinematics, which requires very high-quality proper motions. For open clusters, which are typically loose groups with tens to hundreds of members, there is the additional difficulty of inferring kinematic structures from sparse and irregular distributions of stars. Aims: Here, we aim to analyse internal stellar kinematics of open clusters, and identify rotation, expansion, or contraction patterns. Methods: We use Gaia Early Data Release 3 (EDR3) astrometry and integrated nested Laplace approximations to perform vector-field inference and create spatio-kinematic maps of 1237 open clusters. The sample is composed of clusters for which individual stellar memberships were already known, thus minimising contamination from field stars in the velocity maps. Projection effects were corrected using EDR3 data complemented with radial velocities from Gaia Data Release 2 and other surveys. Results: We report the detection of rotation patterns in eight open clusters. Nine additional clusters display possible rotation signs. We also observe 14 expanding clusters, with 15 other objects showing possible expansion patterns. Contraction is evident in two clusters, with one additional cluster presenting a more uncertain detection. In total, 53 clusters are found to display kinematic structures. Within these, elongated spatial distributions suggesting tidal tails are found in five clusters. These results indicate that the approach developed here can recover kinematic patterns from noisy vector fields, as those from astrometric measurements of open clusters or other stellar or galactic populations, thus offering a powerful probe for exploring the internal kinematics and dynamics of these types of objects.

79 ASTRONOMY AND ASTROPHYSICS↗