skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on September 5, 2026

Title: Probing Turbulence, Gravity, Supernovae, and Magnetic Field Effects with the 6D Kinematics of Young Stars in Milky Way Star-forming Regions
Abstract The dynamics of star-forming gas can be affected by many physical processes, such as turbulence, gravity, supernova explosions, and magnetic fields. In this paper, we investigate several nearby star-forming regions (Orion, Upper Sco, Taurus, and Perseus) for kinematic imprints of these influences on the newly formed stars. Using Gaia DR3 astrometry and APOGEE DR17 radial velocities, we compute first-order velocity structure functions (VSFs) of young stars in galactic Cartesian coordinates in both 6D (3D positions and 3D velocities) and 4D (3D positions and each 1D velocity) to identify signatures of turbulence and anisotropic motion. We also construct 3D and 1D radial velocity profiles to identify coherent expansion trends, and compare stellar proper motions to plane-of-sky magnetic field orientations in Taurus and Perseus. We find that the VSFs are mildly anisotropic, with slightly different amplitudes, slopes, or features in different directions in several groups, but in general, they are all consistent with Larson’s Relation at intermediate length scales, especially in less compact groups. In several cases, the VSFs exhibit features suggestive of local energy injection from supernovae. Radial velocity profiles reveal clear anisotropic expansion in multiple groups, with the most extreme cases corresponding to those with the most anisotropic VSFs. In Perseus, we find that the motions of young stars are preferentially perpendicular to the local magnetic field. We find multiple, overlapping causes in each group for the observed kinematics. Our findings support that young stars remember more than just the turbulent state of their natal clouds.  more » « less
Award ID(s):
2219686
PAR ID:
10637883
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Publisher / Repository:
IOP
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
990
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
165
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The interstellar medium (ISM) is turbulent over vast scales and in various phases. In this paper, we study turbulence with different tracers in four nearby star-forming regions: Orion, Ophiuchus, Perseus, and Taurus. We combine the APOGEE-2 and Gaia surveys to obtain the full six-dimensional measurements of positions and velocities of young stars in these regions. The velocity structure functions (VSFs) of the stars show a universal scaling of turbulence. We also obtain Hαgas kinematics in these four regions from the Wisconsin H-Alpha Mapper. The VSFs of the Hαare more diverse compared to those of stars. In regions with recent supernova activities, they show characteristics of local energy injections and higher amplitudes compared to the VSFs of stars and of CO from the literature. Such difference in amplitude of the VSFs can be explained by the different energy and momentum transport from supernovae into different phases of the ISM, thus resulting in higher levels of turbulence in the warm ionized phase traced by Hα. In regions without recent supernova activities, the VSFs of young stars, Hα, and CO are generally consistent, indicating well-coupled turbulence between different phases. Within individual regions, the brighter parts of the Hαgas tend to have a higher level of turbulence than the low-emission parts. Our findings support a complex picture of the Milky Way ISM, where turbulence can be driven at different scales and inject energy unevenly into different phases. 
    more » « less
  2. For decades we have known that the Sun lies within the Local Bubble, a cavity of low-density, high-temperature plasma surrounded by a shell of cold, neutral gas and dust. However, the precise shape and extent of this shell, the impetus and timescale for its formation, and its relationship to nearby star formation have remained uncertain, largely due to low-resolution models of the local interstellar medium. Leveraging new spatial and dynamical constraints from the Gaia space mission, here we report an analysis of the 3D positions, shapes, and motions of dense gas and young stars within 200 pc of the Sun. We find that nearly all the star-forming complexes in the solar vicinity lie on the surface of the Local Bubble and that their young stars show outward expansion mainly perpendicular to the bubble's surface. Tracebacks of these young stars' motions support a scenario where the origin of the Local Bubble was a burst of stellar birth and then death (supernovae) taking place near the bubble's center beginning 14 Myr ago. The expansion of the Local Bubble created by the supernovae swept up the ambient interstellar medium into an extended shell that has now fragmented and collapsed into the most prominent nearby molecular clouds, in turn providing robust observational support for the theory of supernova-driven star formation. 
    more » « less
  3. ABSTRACT The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in (i) a 3D magnetohydrodynamic simulation, (ii) synthetic observations generated from the simulation at different viewing angles, and (iii) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc to core scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flows along the magnetic field towards dense cores. When comparing the observed cores identified from the Green Bank Ammonia Survey and Planck polarization-inferred magnetic field orientations, we find that the relative core–field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core–field orientation could be used to probe the relative significance of the magnetic field within the cloud. 
    more » « less
  4. Two main mechanisms have classically been proposed for the formation of runaway stars. In the binary supernova scenario (BSS), a massive star in a binary explodes as a supernova, ejecting its companion. In the dynamical ejection scenario, a star is ejected during a strong dynamical encounter between multiple stars. We propose a third mechanism for the formation of runaway stars: the subcluster ejection scenario (SCES), where a subset of stars from an infalling subcluster is ejected out of the cluster via a tidal interaction with the contracting gravitational potential of the assembling cluster. We demonstrate the SCES in a star-by-star simulation of the formation of a young massive cluster from a 106Mgas cloud using theTORCHframework. This star cluster forms hierarchically through a sequence of subcluster mergers determined by the initial turbulent, spherical conditions of the gas. We find that these mergers drive the formation of runaway stars in our model. Late-forming subclusters fall into the central potential, where they are tidally disrupted, forming tidal tails of runaway stars that are distributed highly anisotropically. Runaways formed in the same SCES have similar ages, velocities, and ejection directions. Surveying observations, we identify several SCES candidate groups with anisotropic ejection directions. The SCES is capable of producing runaway binaries: two wide dynamical binaries in infalling subclusters were tightened through ejection. This allows for another velocity kick via subsequent via a subsequent BSS ejection. An SCES-BSS ejection is a possible avenue for the creation of hypervelocity stars unbound to the Galaxy. The SCES occurs when subcluster formation is resolved. We expect nonspherical initial gas distributions to increase the number of calculated runaway stars, bringing it closer to observed values. The observation of groups of runaway stars formed via the SCES can thus reveal the assembly history of their natal clusters. 
    more » « less
  5. Precise Gaia measurements of positions, parallaxes, and proper motions provide an opportunity to calculate 3D positions and 2D velocities (i.e., 5D phase-space) of Milky Way stars. Where available, spectroscopic radial velocity (RV) measurements provide full 6D phase-space information, however there are now and will remain many stars without RV measurements. Without an RV it is not possible to directly calculate 3D stellar velocities; however, one can infer 3D stellar velocities by marginalizing over the missing RV dimension. In this paper, we infer the 3D velocities of stars in the Kepler field in Cartesian Galactocentric coordinates (vx, vy, vz). We directly calculate velocities for around a quarter of all Kepler targets, using RV measurements available from the Gaia, LAMOST, and APOGEE spectroscopic surveys. Using the velocity distributions of these stars as our prior, we infer velocities for the remaining three quarters of the sample by marginalizing over the RV dimension. The median uncertainties on our inferred vx, vy, and vz velocities are around 4, 18, and 4 km/s, respectively. We provide 3D velocities for a total of 148,590 stars in the Kepler field. These 3D velocities could enable kinematic age-dating, Milky Way stellar population studies, and other scientific studies using the benchmark sample of well-studied Kepler stars. Although the methodology used here is broadly applicable to targets across the sky, our prior is specifically constructed from and for the Kepler field. Care should be taken to use a suitable prior when extending this method to other parts of the Galaxy. 
    more » « less