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.
Star Formation in Crowds
YSOs in protoclusters
Postdoc: Allison Towner
PhD: Desmond Jeff, Theo Richardson, Alyssa Bulatek, Nazar Budaiev
REU: Justin Otter, Danielle Bovie, Josh Machado
Undergrad: Madeline Hall, Michael Fero, Derod Deal, Parker Ormonde
Star formation drives the evolution of the universe
After the big bang, there was a brief dark period.
After the dark ages, almost all of the light we see comes from star formation in galaxies.
Galaxies are made of stars, the observable universe is made of galaxies
Star Formation oversimplified
Ṁ
The star formation rate, i.e., how much gas turns to stars
L / M
The light per unit mass, i.e., how stars and stellar populations turn matter into light
Just what's on the slides
High-mass stars produce photons & heavy elements low-mass stars live practically forever
Point color shows effective temperature, point size shows luminosity (left) and mass (right)
describe the plots!
Points are colored by an approximation of how they appear on the sky
100 msun stars produce about a million times the light of a 1 msun star
1 and lower Msun stars live as long as the universe
URL links to script to make figures
The stellar initial mass function (IMF)
Stars are randomly sampled from this distribution
Log-differential probability vs log mass
Sample cluster of 1000 Msun
Points are a real sampling from the IMF, shows what the curve means in practice
URL links to script to make figures
The IMF is generally assumed to be universal in time and space
primarily because of the weakness of evidence of variation
Almost all of the light in star-forming galaxies is produced by high-mass stars
UV image of M31, Andromeda, our closest massive neighbor
Blue light is produced by high-mass stars
(high-mass stars dominate anyway, but even more in UV)
The stars form in and from gas
Herschel FIR dust image of M31
Dust is closely correlated with blue starlight = young stars
Dust traces gas, stars form where there is gas
Most of what we know of star formation in detail comes from small local clouds
Herschel FIR dust image of Taurus mol cld
distance ~140 pc, local neighborhood
Most of what we know of star formation in detail comes from small local clouds
NIR extinction map
distance ~140 pc, local neighborhood
Cartoon of low-mass star formation
The next several slides will walk through this
A molecular cloud fragments
We start with a relatively smooth molecular cloud
Turbulence seeds overdensities
Overdensities become gravitationally unstable and collapse: this is cloud fragmentation
The core forms a central protostar
The protostar heats its parent core and forms a disk
It drives an outflow and consumes or blows out its core
Eventually, you end with just a star-disk system
Cartoon of low-mass star formation
Most of what we know of star formation in detail comes from small local clouds
They contain only low-mass stars and do not represent star formation in general
Local clouds like Taurus lack high-mass stars
Plot shows a histogram of masses; stops at 1 Msun
Within 0.5 kpc, there is only 1 star with more than 10 Msun currently forming
Most stars form in denser regions
Otter+, submitted soon
FOV: 0.07 pc (16000 AU)
72 YSOs
One "hot core"
Lots of stars & outflows in a small environment
Still no massive (OB) stars
Most stars form in denser regions
NGC 1333, an embedded low-mass cluster
Lada & Lada 2003: >70% in embedded clusters
Lots of stars & outflows in a small environment
Still no massive (OB) stars
Most stars form in denser regions
NGC 3603 is a high-mass (104 M⊙ ) cluster
Lada & Lada 2003:
5-10% in bound clusters
in our Galaxy
Dozens of massive OB stars within < 1 pc
Dense, massive clusters like these are bound to their own
self-gravity, they last a long time, until they are destroyed by
dynamical interactions with molecular clouds
Star formation drives the evolution of the universe
Most stars in most galaxies formed long ago
Peak of cosmic star formation was z~2
Galaxies were smaller & denser back then
Y-axis is normalized size
Transition out: How did stars form at the peak of cosmic star formation, when galaxies were denser?
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
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 the CMZ
Γ 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
How do we learn about clustering? The IMF?
Count objects.
Cores are sometimes countable.
Protostars are countable.
YSO counting has taught us a lot in local clouds
The California molecular cloud with protostars
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
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
Our own 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
Our own Galaxy's center, the CMZ, has denser gas than the Galactic average
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
Our own Galaxy's center, the CMZ, has denser gas than the Galactic average
Each zoom reveals more sources.
Most of the stars reside in bound proto-clusters
A higher-resolution view of Sgr B2 Deep South
We catalog the protostellar cores, shown in orange and red
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
revisit the YSO catalog...
these are only high-mass YSOs
... look without the YSOs ... (next slide we go deeper)
The 1mm version, at 0.1" resolution and ~0.3 Msun sensitivity
The same 3mm cores shown on 1mm.
The 1mm reveal many more fainter cores we didn't see before
Desmond Jeff: Hot Cores and YSOs in Sgr B2 DS
Are there appropriate numbers of low-mass YSOs?
Is the IMF 'normal'?
Desmond is a 3rd-year PhD student
Focus on hot core temperature measurement & YSO identification
MYSO counts let us investigate thresholds
Local cloud studies support the idea of a gas density threshold for star formation
Thresholds are used in simulations to say "if gas reaches this density, turn it into stars"
Left: NYSOs vs total cloud mass
Right: NYSOs vs cloud mass above a fixed density
Correlation is tighter in the right plot - stars form in dense
gas
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
From YSO counts to the IMF?
How do we measure the CMF if the cores all have YSOs in them?
Some CMF measurements in clusters have given us tantalizing hints of CMF variation...
YSO modeling → luminosity functions
NSF 2008101: "How are stellar masses set?"
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
Top-heavier IMFs are seen in high-mass clusters, CMFs in protoclusters
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
ALMA-IMF is the next step in CMF measurement & YSO counting
Continuum data paper in prep (Ginsburg+, with big data reduction
team: Roberto Galvan-Madrid, Nichol Cunningham, Timea Csengeri,
Patricio Sanhueza, Fernando Olguin, Thomas Nony, Jordan Molet, Ana
Lopez, Yohan Pouteau, Andrez Guzman, Manuel Fernandez, Melisse
Bonfand)
Self-calibration (10-500% dynamic range improvement)
Mosaicing
Continuum selection
Method comparison
Survey overview paper in prep (Motte+)
Catalogue paper in prep (Louvet+)
Time for a brief mental break, just enjoy the photos
How do we measure masses?
On the top end, mass measurement is difficult:
cores are optically thick
cores are confused & blended
the measured luminosity can be the sum of whole (proto)clusters
Dynamical mass measurements are the gold standard.
Sometimes, we can't measure dynamical masses
Salts in Orion
Artists conception of a disk around a high-mass protostar
There is an outflow
Note the blue hazy inner ring
Orion Source I
a disk around a 15 M⊙ YSO
Salt: NaCl
Orion Source I lives at the center of the BN/KL explosion
We observed the disk known as source I
...in salts! NaCl and KCl
We used that disk to measure the mass of the central source
Observing the Keplerian rotation profile of a disk is the most direct way to measure a protostar's mass
(we can only see the disk, not the star itself)
Left is a model position-velocity diagram, right is the data
Red curve traces the outermost position of orbits of point particles
We can use salts to measure HMYSO masses
NaCl, KCl are only in the disk, not the outflow (water traces both)
NaCl is detected in at least two other HMYSOs
Salts are observable with ALMA, the JVLA, and the future ngVLA
Future projects will involve observing and modeling salt disks to measure HMYSO masses
W51 e2e: Too optically thick at 1mm to measure disk
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?
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?
Looking forward:
PASHION: Paschen Alpha Survey of Hydrogen Ions
JWST: Deep Paα, Brα, and broadband imaging
PASHION: H2RG with Lockheed electronics, three narrow-band filters, 2.5" resolution, 25' FOV
A 24 cm dedicated survey telescope will be the most sensitive Galactic plane survey of ionized gas
These are fiducial numbers for a 1-year mission performing a 100 square degree blind survey. An extended mission may be possible.
PASHION, and JWST, recombination line science
Accretion onto YSOs
HII regions
Assuming typical AV ~2 per kpc
Summary
Most stars form in regions unlike the solar neighborhood
Greater clustering, higher SF thresholds in denser clouds
ALMA-IMF will expand the sample to match or exceed local clouds
We have, and are building more, tools to measure MYSO masses
Salt is a new tool to probe disks around high-mass stars
CS masers may track YSO disks; they require rare conditions
Things I did not talk about today:
MUSTANG Galactic Plane Survey
HCHII regions are ~1/3 as abundant as UCHII
Spectral line survey (B34567) of The Brick [WIP]
Feedback, YSOs in W51
Massive (200-300 M⊙ ), hot (200-600K) cores suppress fragmentation
Multidirectional accretion flows
Accretion flows, outflows in Sgr B2
Turbulence in Sgr B2 [WIP: student Madeline Hall]
Gas temperatures in W51 [WIP: student Josh Machado]
Possible future uses for these lines?
Metallicity measurement in deeply embedded star-forming environments? (at least of Na, K, Cl)
Disk kinematics of high-mass stars, which are otherwise unobservable (τ>1 at mm wavelengths)
Disk kinematic measurements at early stages?
Probe dust destruction (and/or formation?) in outflows, disks?
Probe radiation environment around HMYSOs?
Why do we see salt?
Previously, NaCl & KCl only in AGB* atmospheres,
associated with dust formation
Most likely dust destruction here
Dust destruction happens immediately as the outflow is launched?
What about excitation? We see vibrationally excited lines, which are not seen in AGB*s
We do not have a viable model to explain these temperatures
A strong non-blackbody radiation field from 25-40 µm may explain them.
Forsterite (MgSiO4 ) has some emission bands in that range. Maybe?
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)
YSO disk counts in W51
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
Ammonia Masers
Large scales again: 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
Measuring Line Profiles
SCOUSE uses pyspeckit for manual fits. Gausspy+ is machine-learning trained.
We're exploring more automated approaches.