PaperlayerView exact revision on arXiv

Source-linked AI summary

Evidence for Warped Disks of Young Stars in the Galactic Center

H. Bartko, F. Martins, T. K. Fritz, R. Genzel, Y. Levin, H. B. Perets, T. Paumard, S. Nayakshin, O. Gerhard, T. Alexander, K. Dodds-Eden, F. Eisenhauer, S. Gillessen, L. Mascetti, T. Ott, G. Perrin, O. Pfuhl, M. J. Reid, D. Rouan, A. Sternberg, S. Trippe

arXiv:0811.3903v3astro-ph

TL;DR

The study addresses how young massive stars near the Galactic Center are organized and formed despite the black hole’s strong tidal environment. Using improved NACO and SINFONI observations of 90 WR/O stars with statistical orbital analysis, it finds a warped clockwise system and a coherent counter-clockwise structure. These results support in situ formation in gaseous accretion disks, while eccentricity conclusions remain sensitive to systematic reconstruction effects.

  • Problem

    The origin and dynamical organization of young massive stars near the Galactic Center are uncertain because regular star formation is likely suppressed by the massive black hole’s tidal forces.

  • Method

    The study combines new NACO and SINFONI observations with Monte Carlo simulations and statistical analysis of orbital properties and orientations for 90 WR/O stars.

  • Results

    The clockwise stars form a thin but strongly warped disk-like system, while the counter-clockwise stars show a coherent non-isotropic structure at 98% confidence.

  • Takeaways & Limitations

    The observed warp and steep surface-density distribution favor in situ formation in gaseous accretion disks, with the counter-clockwise population possibly representing disrupted disks or streamers.

  • Takeaways & Limitations

    Firmer conclusions on the eccentricity distribution require more thorough analysis of systematic effects in eccentricity reconstruction.

Abstract

from arXiv · show

The central parsec around the super-massive black hole in the Galactic Center hosts more than 100 young and massive stars. Outside the central cusp (R~1") the majority of these O and Wolf-Rayet (WR) stars reside in a main clockwise system, plus a second, less prominent disk or streamer system at large angles with respect to the main system. Here we present the results from new observations of the Galactic Center with the AO-assisted near-infrared imager NACO and the integral field spectrograph SINFONI on the ESO/VLT. These include the detection of 27 new reliably measured WR/O stars in the central 12" and improved measurements of 63 previously detected stars, with proper motion uncertainties reduced by a factor of four compared to our earlier work. We develop a detailed statistical analysis of their orbital properties and orientations. Half of the WR/O stars are compatible with being members of a clockwise rotating system. The rotation axis of this system shows a strong transition as a function of the projected distance from SgrA*. The main clockwise system either is either a strongly warped single disk with a thickness of about 10 degrees, or consists of a series of streamers with significant radial variation in their orbital planes. 11 out of 61 clockwise moving stars have an angular separation of more than 30 degrees from the clockwise system. The mean eccentricity of the clockwise system is 0.36+/-0.06. The distribution of the counter-clockwise WR/O star is not isotropic at the 98% confidence level. It is compatible with a coherent structure such as stellar filaments, streams, small clusters or possibly a disk in a dissolving state. The observed disk warp and the steep surface density distribution favor in situ star formation in gaseous accretion disks as the origin of the young stars.

1. Introduction

The Galactic Center contains a young massive stellar population whose dynamical organization and formation history remain debated. New observations and statistical simulations analyze 90 WR/O stars to test disk-like structures and their orbital properties.

  • About 100 massive young stars occupy the central parsec, predominantly O-type supergiants and Wolf-Rayet stars with an estimated age of about 6 × 10^6 years.
  • Previous studies inferred two moderately thick counter-rotating disks, but later observations confirmed one disk and found limited evidence for the other.
  • The stars’ recent formation is surprising because tidal forces from the massive black hole likely suppress regular star formation near the Galactic Center.
  • New NACO and SINFONI observations added 27 reliably measured WR/O stars, improved earlier measurements, and reduced proper-motion uncertainties by a factor of four.
  • Using 90 well-measured WR/O stars, the study applies detailed orbital analysis and Monte Carlo simulations to test coherent structures against isotropic stellar distributions.

2. Data

The study assembles a quality-controlled sample of 90 early-type WR/O stars using new SINFONI spectroscopy and multi-epoch NACO astrometry. The resulting measurements support analysis of orbital orientations while retaining important completeness and line-of-sight constraints.

  • Observations: 27 new reliably measured stars supplement 63 previously measured stars in the central 12″.The observations combine SINFONI spectroscopy with proper motions derived from six NACO/CONICA epochs.
  • Orbital classification: 61 stars are classified as clockwise and 29 as counter-clockwise, based on the 90-star sample.The sample’s projected positions and normalized angular momenta are used to characterize the two orbital senses.
  • Sample selection: 90 WR/O stars remain after applying mK < 14 and radial-velocity uncertainty Δ(vz) ≤100 km/s.The sample excludes early B dwarfs and contains no reliably measured WR/O stars inside 0.8″ or beyond 12″.
  • Measurement quality: Proper-motion uncertainties decrease from a mean of 35 km/s to 5 km/s in the present analysis.The improvement is attributed to the larger data set, geometric-distortion correction, and reduced coordinate-system uncertainties.
  • Line-of-sight constraints: Five years of NACO data constrain stellar z-coordinates for only a handful of stars.Measured accelerations provide full orbital solutions for some stars, while acceleration upper limits provide lower limits on |z| for others.
  • Selection effects: The sample is incomplete, with 75% combined photometric and spectroscopic completeness for O stars with mK = 13−14.Completeness corrections are needed for reliable radial-density estimates but are considered unnecessary for the paper’s angular-momentum analysis.

3. MC Simulation of Signal and Background

The simulations generate bound stellar orbits under isotropic-cusp and thick-disk models, then propagate observational uncertainties to evaluate expected angular-momentum distributions. Their design isolates orbital-plane orientation while incorporating distributions of orbital elements and measurement errors.

  • Simulation goals: The simulations compare isotropic stellar distributions with thick-disk models to assess whether observed structures are compatible with isotropy.The analysis evaluates the probability of observed features under simulated background and disk populations.
  • Isotropic cusp: Isotropic models draw angular-momentum directions uniformly on a sphere and generate orbital elements including ω, τ, ϵ, and a.The semimajor-axis distribution follows dNstars/da ∝ a^−β+1 with β = 2 over 0.2″ ≤ a ≤ 40″.
  • Thick-disk model: Thick-disk models distribute angular-momentum directions around a disk axis with two-dimensional Gaussian sigma thickness σψ.σψ is defined as the angular thickness of the simulated disk in orientation space.
  • Model assumptions: The isotropic simulations use only the gravitational potential of Sgr A* and retain a strong prior from the chosen eccentricity and semimajor-axis distributions.The resulting angular-momentum expectations depend primarily on orbital-plane direction, represented by Ω and i.
  • Mock observations: Each simulated star is converted from orbital elements into positions and velocities before measurement errors are added.The injected errors follow the distributions measured for the observed data.

4. Analysis Method to Search for Features in the Star Distribution

The analysis reconstructs stellar angular-momentum directions from uncertain line-of-sight positions, maps their density, and tests coherent structures against isotropic Monte Carlo stars. Simulations show that disk-like angular-momentum concentrations are distinguishable from isotropy, while the adopted z-prior introduces smaller biases than alternatives.

  • Angular-momentum feature search: The disk search tests a common angular-momentum direction against isotropic stars using observed density maps, Monte Carlo mean and RMS maps, and significance maps.The analysis proceeds by comparing the observed angular-momentum density with expectations from simulated isotropic populations.
  • Angular-momentum feature search: Each star contributes 1000 Monte Carlo z-values, reconstructed angular-momentum directions, and densities within a fixed aperture.The line-of-sight coordinate is unknown, so the method propagates this uncertainty through angular-momentum reconstruction.
  • z-coordinate reconstruction: The adopted z-prior produces smaller reconstruction biases than uniform-z and uniform-acceleration priors, although line-of-sight distances remain important for possible biases.The analysis assumes β = 2 and restricts generated positions to bound orbits; tests with β = 1.5 and 2.5 gave similar results.
  • Angular-momentum feature search: 15° apertures maximize the signal-to-noise ratio for moderately thick disks with a 10° two-dimensional Gaussian thickness.The aperture uses flat weights, while a matched filter could provide higher signal-to-noise ratios.
  • Simulated distributions: For simulated isotropic stars, angular-momentum densities are relatively flat, whereas a disk produces a peak at its simulated angular-momentum direction.The isotropic and disk cases are compared using 4 × 10^5 and 4 × 10^4 simulated stars, respectively.
  • Simulated distributions: 17.2 stars per 15° aperture for a simulated thick disk contrasts with 1.6 stars for 90 isotropic stars, whose expected RMS is 0.55.This contrast supports distinguishing a thick disk from an isotropic stellar distribution.

5.1. The Clockwise System

The clockwise stellar system shows a statistically significant change in angular-momentum orientation with projected distance, consistent with a strongly warped disk or multiple radial streamers. Its candidate members have a nonzero mean eccentricity and a steep radial surface-density profile.

  • Significance of the clockwise excess: 12.2σ is the global maximum excess significance for the clockwise system, corresponding to 8.3 stars per 15° aperture versus 1.6 expected for 90 isotropic stars.The excess has a narrow core, and its HWHM is 18°.
  • Significance of the clockwise excess: 35 stars are estimated to contribute to the clockwise excess peak, whose statistical position uncertainty is approximately 3°.This estimate compares the observed peak with a simulated 90-star disk having a 10° two-dimensional sigma thickness.
  • Radial dependence of the excess position: The clockwise excess shifts systematically with radius: inner, middle, and outer intervals peak at distinct angular-momentum positions.The inner and outer peak positions are separated by (64 ± 6)°.
  • Radial dependence of the excess position: The change in excess position is significant at >10σ, and the local angular-momentum direction is compatible with a smooth function of projected distance.A quadratic-polynomial fit describes the radial change better than a great circle, with scatter attributed to Poisson noise and local disk thickness.
  • Interpretation of the clockwise structure: A warped single disk or at least two disks or planar streamers can explain the orientation change; the inner and outer regions are consistent with 10°-thick disks.The middle radial interval appears to be a transition region with a substantial counter-clockwise fraction.
  • Surface density: The best-fit surface-density profile of the clockwise system is Σ(r_disk) ∝ r^-1.95±0.25.
  • Orbital elements: The clockwise candidate stars have a mean eccentricity of 0.51, reduced to 0.42 ± 0.05 when five nearly radial, high-eccentricity candidates are excluded.The five stars may be chance-compatible rather than dynamically belonging to the system, but they cannot be excluded.

5.2. Hertzsprung-Russell Diagram

The Hertzsprung-Russell diagram is used to estimate the ages of the early-type stars by comparison with theoretical isochrones. The four stellar groups are consistent with a coeval population without an age trend with projected distance.

  • Diagram construction: The diagram includes O-type supergiants, giants and dwarfs, and B supergiants, while Wolf-Rayet stars are excluded because isochrones are ill-defined there.
  • Age estimate: 4–8 Myr is the common age range inferred for the four stellar groups in the Hertzsprung-Russell diagram.The estimate comes from comparing stellar positions with theoretical isochrones.
  • Age estimate: No variation of average stellar age with projected distance to Sgr A* is detected.

5.3. The Counter-Clockwise System

The counter-clockwise WR/O stars show a non-isotropic angular-momentum distribution, but the structure is less sharply defined than the clockwise system. The observed excess is compatible with a disk-like or otherwise coherent, possibly disrupted structure.

  • Evidence for a counter-clockwise structure: 29 of 90 WR/O stars, about 30%, are on counter-clockwise orbits, motivating a test for a coherent second structure.
  • Angular-momentum density: 3.8 stars per 25° aperture are found at the counter-clockwise density maximum, compared with 1.5 expected for 15 isotropic counter-clockwise stars.The maximum lies near the previously reported counter-clockwise system.
  • Statistical significance: 2% is the probability that an isotropic distribution produces the observed counter-clockwise density maximum when all aperture directions are searched.Including the full U-shaped excess makes the probability lower.
  • Angular-momentum density: The counter-clockwise excess has an extended U-like morphology with local maxima up to 3.7 stars per 25° aperture.
  • Caveats and alternatives: The counter-clockwise structure may reflect a thick or warped system, dissolution, separate streamers, projection effects, or a small non-isotropic stellar sample.The limited number of stars reduces the significance of radial trends.
  • Disk interpretation: A 10°-thick simulated disk can reproduce a similar U-like reconstructed-angular-momentum shape for a small sample at the observed sky positions.This demonstrates compatibility with a disk interpretation without requiring it.

5.4. Summary and Comparison to Previous Results

The expanded and more precise dataset confirms a warped clockwise system and a non-isotropic counter-clockwise population. The results support simple disk descriptions while allowing for more complex streamers or disrupted structures and leave eccentricity systematics unresolved.

  • Data improvements: 90 WR/O stars are analyzed after adding 27 reliably measured stars and reducing proper-motion uncertainties by a factor of four.
  • Clockwise system: The clockwise angular-momentum direction varies with projected distance, while its innermost edge remains compatible with the previously determined clockwise disk.
  • Clockwise system: The clockwise system has a disk thickness of (14 ± 4)° and an average orbit that is not circular.
  • Limitations: Firmer conclusions about the eccentricity distribution require further analysis of systematic effects in eccentricity reconstruction.
  • Counter-clockwise system: The counter-clockwise stars have a non-isotropic distribution at confidence beyond 98%, compatible with a second disrupted disk or streamer.
  • Interpretation: The stellar distribution may be a connected but complex system of separate streamers at different radii and orientations rather than simple flat disks or rings.
  • Comparison with previous work: The results of Lu et al. (2008) are not in contradiction with this analysis because the present dataset has more stars and smaller errors beyond 3.5″.

6. Discussion

The discussion evaluates disk-like structures around the Galactic Center through observable properties including number, thickness, orientation, warps, radial distribution, and eccentricity. The clockwise system is visualized as a warped and tilted disk, while the counter-clockwise population provides a comparison structure.

  • The clockwise stars form a relatively thin disk-like structure whose rotation axis varies with radius, indicating a strong warp or related streamers.
  • The counter-clockwise stars form a coherent structure that may consist of streamers or a more disrupted second disk at large angles to the primary system.
  • The models are evaluated against disk number, candidate fractions, thickness, orientation, warps, radial distribution, and eccentricity.
  • The three-dimensional visualization assumes locally flat circular orbits from a = 1”–10” and represents angular momentum as a smooth function of semi-major axis.

6.1. Observables

The observed young stars divide into a dominant warped clockwise system and a less massive counter-clockwise structure. Their orbital orientations vary with radius, while the clockwise stars show moderate non-circular eccentricities and a steep radial surface-density profile.

  • 6.1.1. Number of Stellar Structures and their Thickness: 55% of observed WR/O stars are candidate members of a warped clockwise system, while 20% belong to a counter-clockwise system inclined by about 100°.
  • 6.1.1. Number of Stellar Structures and their Thickness: 25% of WR/O stars have angular separations too large for candidate membership, leaving two main structures with different thicknesses and masses.
  • 6.1.2. Warps and Orientation: Δθ = 26° and Δφ = 41° separate the inner and intermediate clockwise regions, followed by Δθ = −34° and Δφ = 36° toward the outer region.
  • 6.1.2. Warps and Orientation: The clockwise angular-momentum direction deviates from a great-circle warp and is reasonably fitted by quadratic radial functions.
  • 6.1.3. Eccentricity: The clockwise disk candidates have mean eccentricity 0.36 ± 0.06 and exclude fully circular orbits beyond 5σ.
  • 6.1.4. Mass Function (MF) and Radial Distribution: The disk-star surface number density follows Σ(rdisk) ∝ r^-1.95±0.25, while current mass-function data remain insufficiently constrained.

6.2. Models for the origin of the young GC WR/O stars

The paper compares infalling-cluster and in-situ formation models for the young Galactic Center stars. The observed warp, steep radial profile, and inclined structures are more compatible with in-situ formation, although simulations leave important features unresolved.

  • Formation scenarios: The infalling-cluster model forms stars several parsecs away and strips them during inspiral, whereas the in-situ model fragments an infalling gas disk into stars.
  • In situ formation: Tidal forces require cloud densities above 6·10^9 M_SgrA*/(4·10^6 M⊙)(R/7”)^-3 cm^-3 for gravitational binding near Sgr A*.
  • In situ formation: A mean eccentricity of 0.36 ± 0.06 may require star formation on an orbital timescale after initial gas compression.
  • In situ formation: Simulations of cloud-cloud collisions can produce strongly warped disks and highly inclined filaments, streams, or small clusters.
  • Warping mechanisms: A counter-clockwise ring of about 5000 M⊙ at roughly 3” could induce precession periods of about 10^7 years in the clockwise system.
  • In situ formation: In-situ formation can produce a steep stellar surface-density profile, Σ(r) ∝ r^-2, compatible with the observed distribution.
  • Infalling cluster scenario: The observed warps strongly disfavor rapid infall of a cluster, which would preserve the initial angular-momentum orientation in a flat, non-warped distribution.
  • Infalling cluster scenario: The infalling-cluster model predicts Σ(r) ∝ r^-0.75, considerably shallower than the observed radial profile.

6.3. Summary

The observations support in-situ formation of WR/O stars in a clockwise disk plus a highly inclined counter-clockwise structure, while two-cluster infall does not readily produce the observed warp.

  • The observations are compatible with in-situ formation in a clockwise disk and another highly inclined counter-clockwise structure, possibly a disk.
  • Although two cluster infalls could explain two stellar systems, the stars would remain in disks too briefly to develop the observed warp.

7. Conclusions and Outlook

The 90-star sample reveals a warped, non-flat clockwise system and a coherent counter-clockwise feature, strengthening the case for in situ star formation in the Galactic Center.

  • 55% of WR/O stars are candidate members of a clockwise disk.
  • A (64 ± 6)° warp in the clockwise disk is detected at > 10 sigma, with angular momentum varying by projected distance.
  • The clockwise system is non-flat, with inclinations compatible with a two-dimensional Gaussian distribution having σ = 10°.
  • The clockwise system has mean eccentricity 0.36 ± 0.06, indicating non-circular stellar orbits.
  • 20% of WR/O stars form a coherent feature among counter-clockwise stars, with a 2% chance of arising from an isotropic cusp distribution.
  • The warped disk and coherent young-star structures strengthen the conclusion that these stars formed in situ from gas about 6×10^6 years ago.

A. Coordinate systems

The analysis uses Cartesian offsets from Sgr A* together with spherical coordinates, while orbital elements are defined relative to the chosen coordinate system.

  • Cartesian offsets use x eastward, y northward, and z along the line of sight away from the observer.
  • The coordinate system also uses spherical coordinates (x, y, z) ⇔ (r, θ, φ).
  • Figure 25 defines the classical orbital elements within this coordinate system.
Loading 0811.3903v3…