Slide 0 - Instructions :
To move slides, use the arrow keys or swipe on your mobile device
To see the speaker notes, press "s"
To go to full screen, press "f"
To print as PDF, go to this URL: ?print-pdf , then print.
To get a PDF with speaker notes, add ?print-pdf&showNotes=true to the URL.
ALMA-IMF vs Galactic Center 2026
Postdocs: Nazar Budaiev [PhD 2026] , Theo Richardson [PhD 2025] , Miriam Garcia Santa Maria (2024-2025), Allison Towner (2020-2023)
PhD: Desmond Jeff [PhD 2025] ,
Alyssa Bulatek [PhD 2026] ,
Savannah Gramze ,
Taehwa Yoo
Postbac: Aden Dawson
Undergrad: Antonio Daley, Ethan Bhula, Ani Vijarayaman, Avery Lacon, Laya Damaraju, Prashant Sikhdar, Giovanna Guida
Supported by NSF 2008101, 2206511, CAREER 2142300, STSCI 1905, 2221, 3523, 5365, 6151, Astropy
The Central Molecular Zone (CMZ) is one extreme of star-forming conditions in the Galaxy
The CMZ
$\sim10^8$ M$_\odot$ of gas in $\sim200$ pc, 10% of Galactic star formation
$>\frac{1}{3}$ of CMZ SF is in bound clusters (3-8$\times$ local)
$\sim$50% of CMZ SF occurs in the Sgr B2 cloud.
The CMZ: ACES (the ALMA CMZ Exploration Survey)
ACES: the first results
The survey overview, four data-release papers, and a wave of early-science results — the first ~10 ACES papers (2024–2026)
I. The survey: a contiguous ALMA view of the inner ~200 pc
~2″ (~0.07 pc) Band-3 mosaic of the whole CMZ; spatial dynamic range ~103.5 links 100-pc gas flows to protostellar cores.
II. CMZ-wide 3 mm continuum
A contiguous 3 mm continuum mosaic (dust + free-free) across the CMZ; catalog of compact sources forthcoming (Wallace+).
III. Gas kinematics at 0.2 km/s: a resolution leap
HNCO & HCO+ resolve CMZ-wide motions (~15″ Mopra → ~2″ ACES); ubiquitous parsec-scale linear HCO+ absorption features.
IV. Shock tracers & isotopologues at 1.7 km/s
SiO, SO, H13 CO+ , H13 CN, HN13 C, HC15 N maps.
Evaluation of image combination methods (feather is fine)
V. A rich spectrum: COMs and ionized gas
The broad windows (CS, SO, CH3 CHO, HC3 N, H40α) capture dozens of lines — complex molecules plus the H40α recombination line tracing ionized gas.
VI. The CMZ is highly filamentary
ACES resolves a population of small-scale filaments threading the entire CMZ — a new census of these structures.
The MUBLO: a millimeter ultra-broad-line object
A bright, very broad-line (~160 km/s), cold, dusty mm source with no IR/radio counterpart and no SiO — still unexplained.
The M0.8–0.2 ring: a hypernova remnant?
A ~106 M⊙ expanding shell (~20 km/s, >1051 erg) — most plausibly carved by a single high-energy hypernova.
An eccentric circum-nuclear gas flow around Sgr A*
CS(2–1) draws an elliptical track in the longitude–velocity diagram — an eccentric orbit of dense gas in the inner ~10 pc.
Magnetic fields support the large filaments
Comparing HNCO filaments to SOFIA/JCMT polarization: many filaments sit in magnetically dominated environments that resist collapse.
3 mm core luminosity functions
Distance-normalized 3 mm luminosities
3 mm luminosities, now including HII regions in the Sgr B2 sample
3 mm luminosities, Sgr B2 split into its sub-regions (N+M vs DS)
3 mm luminosities by sub-region, HII regions included
1 mm core luminosity functions
Distance-normalized 1 mm luminosities
1 mm luminosities, including HII regions in the Sgr B2 sample
1 mm luminosities, Sgr B2 split into its sub-regions (DS vs N+M)
1 mm luminosities by sub-region, HII regions included
The same diagnostic for the ALMA-IMF getsf cores, colored by protocluster (3 mm left, 1 mm right)
ALMA-IMF cores: spectral index (1 vs 3 mm) vs peak flux; vertical lines mark free-free (α=−0.1), optically-thick (α=2) and optically-thin (α=3.7) dust
Multi dataset CMD
Budaiev, Yoo, Ginsburg (in prep)
Most stars form in dense regions
Lada & Lada 2003:
5-10% in bound clusters
in our Galaxy
JWST life cycle of clusters: W51e → W51 IRS2 → Wd2 → Wd1
Massive star clusters are not isolated
Massive star clusters are not isolated
Massive star clusters are not isolated
forming cluster [ALMA] alongside the giant HII region
Cores fragment.
Massive stars cook their neighbors in
hot cores .
Hot Cores
Hot cores are chemically rich sites of high-mass star formation.
Bounded at ~100 K
(no YSO, 1 YSO, >1 YSO) = (25%, 30%, 45%) W51-E, (53%, 39%, 8%) IRS2
Taehwa Yoo+: fragmentation toward W51
Green ellipses are HII regions.
Arrows: cyan have alpha > 1.7, green < 1.7, yellow 1.3mm only, magenta 3 mm only
(no YSO, 1 YSO, >1 YSO) = (25%, 30%, 45%) W51-E, (53%, 39%, 8%) IRS2
Taehwa Yoo+: fragmentation toward W51
Green ellipses are HII regions.
Arrows: cyan have alpha > 1.7, green < 1.7, yellow 1.3mm only, magenta 3 mm only
(no YSO, 1 YSO, >1 YSO) = 17%, 37%, 45% Sgr B2
("cores" are bigger)
20% of YSOs in W51 lie within the $>100$ K contour
17+ PPOs in hot cores in IRS2 out of 93 =19%
13+ in e2
10+ in e1
out of 118 = 18%
Trend: More massive cores fragment more.
This is the same effect you'd see with sensitivity limits; e.g., all cores have the same fragmentation rate, but you can only see the fragments in the most massive cores.
Suggestive Trend: More massive cores make more massive single stars
Yoo+ resub.:
W51
Budaiev+ 2024 :
Sgr B2
The number of detectable fragments goes up, but the mass gets more
concentrated into the most massive object.
There is an overall asymmetry in the star formation seen both in JWST (HII regions) and ALMA (embedded YSOs)
What ALMA sees, JWST doesn't: 0/700 point source matches, though still investigating extended sources (HII regions, outflow cavities)
480
405
187
182
162
140
360
335
210
480
410
405
Interesting feature highlight reel
The extinction-producing layers seen at short wavelengths are not the ALMA star-forming gas:
this hints at edge-on filamentary (or planar) structures
The Galactic Center Dust ridge...
...has a foreground cloud in front
YSO model grid is testing accretion models
Sgr B2 is forming & fragmenting YSOs - asymmetrically
W51 fragments are dominated by most massive
JWST measures ice abundance across the galaxy
We can measure metallicity with ice
past here is probably stuff to exclude
W51 e2e: Massive YSO
No emission lines from a disk were detected
Continuum appears to be high optical depth in inner ~few hundred AU
CS v=0 J=1-0 and v=0 J=2-1 masers may trace the disk?
M = 24-10 +12 M⊙
if
the masers trace a disk
If CS traces a disk (which requires assuming very, very hard),
the implied mass is ~24 Msun, which is reasonable.
CS maser conditions
van der Walt+ 2020
Top: CS J=1-0, Bottom: CS J=2-1
Red: Consistent w/W51e2e observations
Masers do not coexist; require different specific CS column (N2-1 =1015.6 , N1-0 1016.1 cm-2 )
Require high abundance (XCS > 10-5 )
Hot (300-500 K), moderate-density (n~105 cm-3 ): Disk surface? Or outflow cavity wall?
ALMA-IMF: The CMF is shallow (top-heavy) in HMSFRs
ALMA-IMF: The CMF is shallow, and steepens with time?
Louvet+, subm
IRDC sample is from Sanhueza+ 2019 (ASHES)
ALMA-IMF Line Data: CH3 CN, CH3 CCH
Temperature measurements with per-pixel rotation diagrams
Hot cores in ALMA-IMF: From rare objects to a population
Cores with line forests
TD >50 K
TG ≳100K
Hot core overview:
9 HCs in W43-MM1 (Brouillet+ 2022 )
~60-70 HCs in ALMA-IMF sample from CH3 OCHO (Bonfand+ 2024 ; left)
~10% of continuum cores are within hot cores
CH3 CN temperature maps (Wyrowski+ in prep)
ALMA-IMF Key Results summary
Rich, science-ready data (Ginsburg+ 2022, Cunningham+ 2023)
CMF is top-heavy in HMSFRs (Pouteau+ 2022, Louvet+ 2024)
Core fragmentation is not 1-to-1 [WIP] (Budaiev+ (CMZ), Yoo+ (W51), Louvet+ (W43), ...)
10% of $M\gtrsim1\mathrm{~M}_\odot$ cores are hot cores (Brouillet+ 2022, Bonfand+ 2024, Wyrowski+)
Outflow feedback builds over time, sets initial conditions for many cores (Nony+ 2022, Towner+ 2023)
Briny chemistry:
Salts, SiS, SiO, and PN are seen together [but limited spectral coverage]
Clusters are sites of interactions & collisions
The BN/I/x interaction is the poster case of accretion ended by dynamical interaction.
...even though Orion is only a medium-mass proto-open-cluster
Classic model: a "core mass function" maps to the IMF
...but the naïve version doesn't work
ALMA-IMF "Core" Mass CDF
The CMF in protoclusters is shallower than the IMF
Louvet+ 2024
State of the art in local clouds is: Herschel Gould Belt Survey discovered 100s-1000s of 'cores', most of which are starless and all <10 Msun
Context: Recent CMF measurements
Many! Differ in: resolution, algorithm, selection.
Nearby surveys are only-starless; more distant are starless + protostellar
How does the CMF become the IMF?
Many cores are protostellar - so there's a mass correction
The path to better mass measurement:
YSO modeling → luminosity functions
Theo Richardson
Models are parametric
They do not necessarily represent reality
Their parameter space is limited
These models have a central temperature and luminosity, not a central mass: linking mass to luminosity requires prestellar evolutionary tracks
The path to better mass measurement:
YSO modeling → luminosity functions
NSF 2008101: "How are stellar masses set?"
Theo Richardson
Hot cores in the Galactic center: Distributed MYSOs
Desmond Jeff+ 2024
Ten hot cores in Sgr B2 DS
outside the massive clusters
TG ~ 200-500 K
M ~ 200 - 2900 M⊙
(proto-O-stars / clusters)
~5% of cores are hot cores
Desmond is a 6th-year PhD student, Reber
Focus on hot core temperature measurement & YSO identification
Sgr B2 DS: More massive cores than the Disk
Several different scenarios: mix of mechanisms
How is star formation in high-mass clusters different?
Feedback from one star affects many in clustered regions
IMF depends on density, feedback, global conditions
Total star formation efficiency is higher.
Collisions assemble the most massive stars?
Cartoon of high- and low-mass star formation
Main difference: massive stars affect their surroundings
This is the 'classical' picture of HMSF; I'll riff on it later
The key differences are the size of the feedback-affected region
and the type of feedback:
The classical picture is that high-mass stars end star formation
once they hit the main sequence / stop accreting (which are the same thing)
Classic HII region feedback: O-stars clear out their environment
Destructive feedback from expanding HII regions: the massive stars produce UV light,
ionizing and evacuating the gas
Right panel is by Anna McLeod using MUSE, done while she was a student at ESO
Accreting massive young stars affect their environment
The classic picture is incomplete: HMYSOs have a big effect
*while* they are accreting too
There is a HMYSO at the center (W51e2e)
The surrounding greyscale circular hot core is seen in
methanol emission
Data are from ALMA
Accreting massive young stars affect their environment
The outflow points to the central accreting star
The outflow implies the existence of a disk
Accreting massive young stars affect their environment
ALMA long-baseline data reveal this in more detail...
Radial profile of temperature around the source
Typical molecular clouds have temperatures 5-25 K, compared to the 100-500 K here
Methanol freezes onto grains, usually locked up in water ice:
it evaporates/desorbs at around 80-100 K
(transition) What effect does this high temperature have?
The characteristic fragmentation scale
The Jeans Mass MJ is the mass where gravity and thermal pressure are balanced.
MJ ∝ T3/2 ρ−1/2
Higher density means greater gravity, leading to collapse and fragmentation
Higher temperature means greater pressure, providing support
against gravitational collapse and suppressing it.
The characteristic fragmentation scale is larger
Jeans Mass
MJ ∝ T3/2 ρ−1/2
(describe plot)
In the warm neighborhood of accreting massive stars, the Jeans mass is higher
Despite the high density, the Jeans mass is 5-10 times greater than locally
Feedback affects dense gas
ALMA + VLA + GBT together give multiple temperature probes on multiple scales.
High-mass protoclusters are filled with gas warmed by feedback.
Floor temperature within 1 pc is ~30-40 K
Radial dependence indicates that the heating is internal
feedback-driven (as opposed to external, from surrounding OB
association)
The cartoon in the context of HMSF
These high mass cores suppress low-mass star formation (LMSF) in their vicinity.
They reduce or prevent LMSF in the cores of stellar clusters.
The 'hot cores' contain a substantial fraction of the molecular cloud mass
The densest regions are the most affected
More extreme: 'cooperative accretion'
With enough high-mass stars forming concurrently, massive stars may prevent fragmentation entirely.
If they still have enough gravity to bind the gas, the remaining gas is
forced onto the most massive gravitational sinks.
High-mass stars can themselves form in a clustered fashion
This is a speculative idea so far not backed up by simulation,
but it is a plausible way to obtain an overpopulation of very
massive stars at the cores of high-mass clusters
Is star formation different in other environments? In the Galactic Center?
Is star formation different in other environments? In the Galactic Center?
Is star formation different in other environments? In the Galactic Center?
Is star formation different in other environments? In the Galactic Center?
Is star formation different in other environments? In the Galactic Center?
The CMZ
$\sim10^8$ M$_\odot$ of gas in $\sim200$ pc, 10% of Galactic star formation
Zoom-in to the Sgr B2 cloud:
Most massive, actively star-forming cloud in the galaxy: 107 M⊙
ALMA 3mm image showing mix of dust and free-free emission
The Central Molecular Zone of the Galaxy represents one extreme of star forming conditions in the Galaxy
The proto-Super Star Cluster Sgr B2 dominates SF in the CMZ
Two $>10^4$ M$_\odot$ clusters within $<3$ pc; $10^6$ M$_\odot$ of gas in $\sim10$ pc
Sgr B2 accounts for most (50-80%?) of the CMZ's SF
Zoom-in to the Sgr B2 cloud:
Most massive, actively star-forming cloud in the galaxy: 107 M⊙
ALMA 3mm image showing mix of dust and free-free emission
In denser (parts of) galaxies, more stars form in clusters
Γ is the fraction of stars forming in bound clusters
Galaxy averages
Y-axis is fraction of stars forming in bound clusters, the "cluster formation efficiency" (CFE)
Milky Way point is the one w/large errors
The "Bound Cluster Fraction" is predicted higher in the CMZ
Γ is the fraction of stars forming in bound clusters
Galaxy averages
CMZ prediction
X-axis is gas surface density, can be approximately treated as same as other plot
Local model is blue dotted (depends on local mean density),
global model is red (depends on Galaxy scale height)
The "Bound Cluster Fraction" is higher in CMZs
Γ is the fraction of stars forming in bound clusters
Galaxy averages
CMZ prediction
Sgr B2 data
CFE is clearly higher in denser regions in our Galaxy
The "Bound Cluster Fraction" was higher in the past
At earlier times, more stars in clusters
Spikes occur during mergers, when density drives up
Observations of Sgr B2 help validate this model
The proto-Super Star Cluster Sgr B2 dominates SF in the CMZ
Two $>10^4$ M$_\odot$ clusters within $<3$ pc; $10^6$ M$_\odot$ of gas in $\sim10$ pc
Zoom-in to the Sgr B2 cloud:
Most massive, actively star-forming cloud in the galaxy: 107 M⊙
ALMA 3mm image showing mix of dust and free-free emission
The proto-Super Star Cluster Sgr B2 dominates SF in the CMZ
Two $>10^4$ M$_\odot$ clusters within $<3$ pc; $10^6$ M$_\odot$ of gas in $\sim10$ pc
Each zoom reveals more sources.
Most of the stars reside in bound proto-clusters
Comparison of Sgr B2 (one cloud, 200 pc2 ) with ALMA-IMF (15 SFRs, 53 pc2 )
3mm catalog → Simplistic mass inferences
At this sensitivity, all are M>8⊙ YSOs
Mtot = 8 M⊙
∫0 ∞ M N(M) dM
∫8 ∞ M N(M) dM
= 96 M ⊙
Compare to Herschel Orion Protostar Survey extrapolated to 3mm
Most lumionus YSO in HOPS is 2000 Lsun; as long as Lsun ~ L3mm, loose implication is that
all Sgr B2 3mm YSOs are L>104 Lsun, or > 8 Msun ish
The "Bound Cluster Fraction" is higher in CMZs
Γ is the fraction of stars forming in bound clusters
Galaxy averages
CMZ prediction
Sgr B2 data
CFE is clearly higher in denser regions in our Galaxy
The CMZ: ACES
The first complete survey of the CMZ with 2.4" resolution (previous best was ~15") between ~2 microns and 10 cm.
Results forthcoming! Continuum, line data papers, catalogs, kinematic analyses, filament identifications all in prep....
... but the first result is unrelated (?) to star formation:
Brinary disks
The SrcI disk has gas-phase salt (NaCl, KCl) and water (H2 O).
So it's brine.
(blame Adam Leroy for this term)
IRAS16547A/B (
Tanaka+ 2020 )
have (unresolved) salt water disks
Brine lines measure dynamical mass
NaCl v=1 J=18-17
Stack of v=[0,1] Ju =[18,17]
SrcI
15 M⊙
30 M⊙
40 M⊙
What governs the star formation rate?
Kennicutt-Schmidt relation ties gas surface density and
stellar surface density
Good correlation over several orders of magnitude, but the
scatter is substantial and important
The timescales are 100 Myr to 10 Gyr and show the "depletion
scales", the time for the gas to fully convert to stars
Turbulent ISM models
Turbulent models are the most popular and successful global
star formation models
Probability of a density vs density
Turbulent ISM models
Some of these authors describe a modified distribution
Turbulent ISM models
Gas above the threshold density forms stars
Continued turbulent driving repopulates that high density,
which gives the star formation rate
ALMA enables protostar counting in distant, massive clouds
Sgr B2: the most massive & star-forming cloud in the Galaxy
An example of one of the clouds from the previous plots
A column density map with the locations of protostars shown as red dots
Protostars are identified in nearby clouds with IR telescopes like Spitzer and Herschel
Column density measured with extinction or dust maps, resolution limited to 30" to a few arcminutes
Protostars are counted and compared to the density they reside in
Compare YSO counts in Sgr B2 and the CMC
Comparison of Sgr B2 to the California cloud on the same spatial scale
Contours show different levels of column density and K-band extinction
Lowest contour in Sgr B2 (green) is about an order of
magnitude above highest contour in CMC
Sgr B2 sources are high-mass (>8 Msun, approx)
They represent a population of sources
CMC stars are individuals
With ALMA, we can get down to the 0.5 Msun scale eventually;
future programs (i.e., ALMA-IMF) will do that
Sgr B2 sources reside mostly in the highest contours
Y-axis shows the column density in which the protostars reside, i.e., about half of the protostars are in regions with column
density > 1024 cm-2
Lada+ 2010 threshold is well below our lowest contour
A threshold separates Sgr B2 from The Brick
Dark curve shows the same cumulative background column density
distribution from the previous slide
"The Brick" is a CMZ cloud with little star formation but
still "very high" density compared to local clouds (a few protostars
seen, and they reside in that very high end)
Only a tiny fraction of The Brick's area has a column density
that overlaps with densities at which protostars form in Sgr B2
The difference between the clouds is consistent with there being
a threshold that The Brick has not met, but Sgr B2 has
Walker+ 2021
1mm cores in Sgr B2 DS
Compared to Dan Walker's
3mm Luminosity Function
3mm luminosity function - all sources are confirmed, high-confidence now
Some are HII regions, most are unclassified
3LF is not an IMF
What are the sources?
At this sensitivity, all are M>8⊙ YSOs
Compare to Herschel Orion Protostar Survey extrapolated to 3mm
Most lumionus YSO in HOPS is 2000 Lsun; as long as Lsun ~ L3mm, loose implication is that
all Sgr B2 3mm YSOs are L>104 Lsun, or > 8 Msun ish
Mass measurements: Optically thin, isothermal dust
T
D estimated with PPMAP fits to Herschel data
(~6" resolution)
Simple models assuming TD ~ f(M) don't change CMF much
We can do better with YSO models and TG measurements
From YSO counts to the IMF?
How do we measure the CMF if the cores all have YSOs in them?
Mapping accretion histories
(left: IS, right: TC)
onto the Robitaille 2017 model grid
Key addition: Envelope mass!
SPICY-ALMA-IMF:
Richardson-enhanced Robitaille+ 2017 model grid
fits including ALMA data
UG team:
Sydney Petz
Brice Tingle
Morgan Himes
Our Galaxy's center, the CMZ, has denser gas than the Galactic average
Cold Dust
Hot, ionized gas
Hot dust/PAHs
Cold dust is where new stars will form
Average molecular cloud densities about 1-2 orders of magnitude greater than in the Galactic disk
Clouds are warmer and more turbulent than the disk
Cloud conditions analogous to early universe [e.g., Kruijssen & Longmore 2013 ]
This is a "local"/"nearby" region where we can test theories that predict, e.g., higher
cluster formation efficiency at higher density
Final segment:
Astrochemistry:
Hot Cores
Ices
Line Surveys
F187N F182M F150W
F210M F187N F182M
F212N F210M F187W
F300M F212N F210M
F360M F300M F212N
F405N F360M F300M
F410M F405N F360M
F466N F410M F300M
F466N F410M F405N
F480M F410M F360M
F480M F466N F410M
F200W F182M F115W
F212N F200W F182M
F356W F212N F200W
F410M F356W F212N
F444W F410M F356W
F466N F444W F410M
The Brick isn't forming many stars
The Brick is icy
The Brick is icy
anticorrelation with extinction indicates ice absorption;
above a given threshold, no higher-extinction stars are seen (that pass
our criteria of being detected in all 6 filters)
Building new tools: Spectral Line Survey of The Brick
Alyssa Bulatek: Spectral Line Survey of The Brick
Dasar line at 107 GHz
CH3 OH only in absorption
Occurs at density $n<10^6 \mathrm{cm}^{-3}$
(me),
Desmond Jeff ,
Savannah Gramze ,
Theo Richardson ,
Nazar Budaiev ,
Alyssa Bulatek ,
Taehwa Yoo ,
Miriam Garcia Santa Maria
Higher CMF + top-heavy cluster IMF = top-heavy CMZ IMF
High-mass clusters (denser regions) have top-heavier IMFs
Observations of clusters support the idea that star formation
in the densest clusters is indeed different
Right panels show IMFs with high-end slopes of 2.3 and 1.75
respectively; they show 1000-Msun clusters with different IMFs
The lower-right plot extends out to 100 Msun stars, while the top
cut off around 30
Plots are produced using IMF figure code again
870 μm Dust
[APEX/ATLASGAL]
870 μm Dust
[APEX/ATLASGAL]
870 μm Dust
[APEX/ATLASGAL]
870 μm Dust
[APEX/ATLASGAL]
870 μm Dust
[APEX/ATLASGAL]
870 μm Dust
[APEX/ATLASGAL]
870 μm Dust
[APEX/ATLASGAL]
870 μm Dust
[APEX/ATLASGAL]
870 μm Dust
[APEX/ATLASGAL]
870 μm Dust
[APEX/ATLASGAL]
870 μm Dust
[APEX/ATLASGAL]
870 μm Dust
[APEX/ATLASGAL]
870 μm Dust
[APEX/ATLASGAL]
870 μm Dust
[APEX/ATLASGAL]
ALMA-IMF data highlights
Collapse is slower on larger scales, but fast enough to matter.
Time for a brief mental break, just enjoy the photos
ALMA-IMF: Continuum Data → core catalogs
Pouteau+ 2022 W43-MM2/3
ALMA-IMF: CMF measurement & YSO counting
(these papers are all now submitted or accepted, but this is a screenshot that includes everyone's photos...)
So how do high-mass stars form?
Is HMSF different from LMSF?
Yes! They have more neighbors.