Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall...
Transcript of Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall...
![Page 1: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/1.jpg)
Dynamics of Nuclear Star Clusters David Merritt*
“Central Massive Objects”ESO Garching, 22-25 June 2010
*with:Fabio Antonini
Pat CoteRoberto Capuzzo-Dolcetta
Alessandra Mastrobuono BattistiSlawomir PiatekEugene Vasiliev
....Image: NASA/JPL-Caltech/S. Stolovy (Spitzer Science Center /Caltech)
Thursday, July 1, 2010
![Page 2: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/2.jpg)
I. Formation via cluster infall
II. Evolution over relaxation time scales
III. Triaxiality and its consequences
Thursday, July 1, 2010
![Page 3: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/3.jpg)
Rees 1988 Milosavljevic 2008
Black Holes Nuclear Star Clusters
I. Formation
Thursday, July 1, 2010
![Page 4: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/4.jpg)
4
Why Migration of Star Clusters?
1. Infall time scales for globular clusters are roughly correct
2. While they often contain young stars, the dominant populations of NSCs are old
3. Migration of 104-106 M⊙ objects to the center is common to both gas- and stellar-dynamical models
4. A massive “seed” is probably required for gas infall scenarios to work
5. Stars are easier than gas
Thursday, July 1, 2010
![Page 5: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/5.jpg)
5
Miocchi & Capuzzo-Dolcetta (2009)
Also:
Fellhauer & Kroupa(2002)
Bekki et al. (2004)
Thursday, July 1, 2010
![Page 6: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/6.jpg)
Ef = Ei + Eorb + Ecl
Mj+1 = (j + 1)M1
jEj+1 = (j + α)Ej + jE1
(j + 1)2R−1j+1 = j(j + α)R−1
j +R−11
Ei,f = initial, final energy of nucleus
Ecl = cluster internal energy = −Gm2/2r
Eorb = cluster orbital energy = −αGmMi/2Ri, α ≈ 1
(M1 = m)
Sequential Mergers of Clusters
where
Then
Energy conservation implies
Thursday, July 1, 2010
![Page 7: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/7.jpg)
7
α = 1.2
α = 0.8 Radius-mass relation
Sequential Mergers of Clusters
Thursday, July 1, 2010
![Page 8: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/8.jpg)
M1 = �mGC�
M1 = 10�mGC�
Predicted relation
Piatek et al. (unpublished)
Thursday, July 1, 2010
![Page 9: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/9.jpg)
r ≈�
M•MGC
�1/3
RGC
r
rinfl≈ 2
�M•
107M⊙
�−1/6 � MGC
106M⊙
�−1/3 �RGC
3pc
�i.e.
∴ The smallest NSCs should have sizes ~ a few x rinfl in galaxies with SMBHs
Complete disruption occurs when
E. g. M● = 106 M⊙, rmin ~ 3 pc ? M● = 107 M⊙, rmin ~ 10 pc ✓
Now Add a Supermassive Black Hole...
Thursday, July 1, 2010
![Page 10: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/10.jpg)
Cluster:
MGC = 4x106M⊙
(untruncated) = 1.1x106M⊙
(truncated) σ(0) = 35 km s-1
Galaxy:
ρ(1 pc) = 400 Msun pc-3
ρ ~ r -0.5
A. Battisti et al.
Thursday, July 1, 2010
![Page 11: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/11.jpg)
Cluster:
MGC = 4x106M⊙
(untruncated) = 1.1x106M⊙
(truncated) σ(0) = 35 km s-1
Galaxy:
ρ(1 pc) = 400 Msun pc-3
ρ ~ r -0.5
A. Battisti et al.
Thursday, July 1, 2010
![Page 12: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/12.jpg)
Cluster : Σ = Σ0
�1 +R2/r20
�−1
Galaxy : lnΣ = lnΣe − b(R/re)1/n + 1
11
Surface Density Profiles First Infallt Second Infallt Third Infallt
Battisti et al. (2010)
Thursday, July 1, 2010
![Page 13: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/13.jpg)
12
Mass-Radius Relation
Initial cluster radius BH influence radius
R ∝ M
Battisti et al. (2010)
Thursday, July 1, 2010
![Page 14: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/14.jpg)
13
First Infallt
Battisti et al. (2010)
Thursday, July 1, 2010
![Page 15: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/15.jpg)
14
Second Infallt
Battisti et al. (2010)
Thursday, July 1, 2010
![Page 16: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/16.jpg)
15
Third Infallt
Battisti et al. (2010)
Thursday, July 1, 2010
![Page 17: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/17.jpg)
16
Growth Timescales - dE Galaxies
Lotz et al. (2001):Examined globular cluster systems in 51 Virgo, Fornax dE galaxies.
Summed radial distributions Observed/predicted NSC magnitudes
Thursday, July 1, 2010
![Page 18: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/18.jpg)
ρ(r) = ρa(r/ra)−γ
∆t =C(γ)
lnΛ
r3aMGC
�ρaG
�1/2�r0ra
�(6−γ)/2
∆t1/2 ≈ 3× 1011yr
�Re
1kpc
�1.8 � MGC
105M⊙
�−1
≈ 9× 1010yr
�ra
1kpc
�3 � MGC
105M⊙
�−1 � ρa1M⊙pc−3
�1/2 � r0ra
�(6−γ)/2
17
Assume that the mass density of the galaxy follows
The time for a GC at initial radius r0 to spiral in to the center is
The time for GCs initially within Re to spiral to the center is
Growth Timescales - gE Galaxies
Thursday, July 1, 2010
![Page 19: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/19.jpg)
Cote et al. (2007)
Thursday, July 1, 2010
![Page 20: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/20.jpg)
Half-mass relaxation times*
• galaxy
✡ NSC
II. Relaxation
*assuming no SMBHs
Thursday, July 1, 2010
![Page 21: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/21.jpg)
τheat ≈�ρnucρgal
�1/2 �Vnuc
Vgal
�1/2
(tnuctgal)1/2
� (tnuctgal)1/2
Net effect of two-body relaxation depends on whether the galaxy is “hotter” or “colder” than the NSC.
If the galaxy is hotter, it transfers heat to the NSC on a timescale:
roughly the geometric mean of the galaxy and NSC relaxation times.
Dokuchaev & Ozernoi (1985)Kandrup (1990)Quinlan (1996)
Thursday, July 1, 2010
![Page 22: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/22.jpg)
Mnuc
Mgal� 104
�ξ−1
100
�2 �rnucrgal
�5
τheat � τcc ≡ ξ−1tnuc 10 � ξ−1 � 300
This heat transfer will reverse core collapse if
i.e.
Mnuc
Mgal= A
�rnucrgal
�B
More generally, one finds a critical size:
above which a NSC expands rather than contracts; A and B depend (weakly) on the galaxy density profile
Thursday, July 1, 2010
![Page 23: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/23.jpg)
Merritt (2009)
Core Collapse vs. Core Expansion
Thursday, July 1, 2010
![Page 24: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/24.jpg)
Mnuc
Mgal= A
�rnucrgal
�B
Evidence of NSC Evaporation?r n
uc/r g
al }expand
collapseO tevol ⪆ 10 Gyr
● tevol ⪅10 Gyr
Merritt (2009)
Thursday, July 1, 2010
![Page 25: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/25.jpg)
M•/Mgal = 0.03
Adding a Black Hole
Merritt (2009)
A large BH reverses the temperature gradient, slowing the transfer of heat from galaxy to NSC
A small BH inhibits core collapse, causing the NSC to expand more quickly
M•/Mgal = 10−4
Thursday, July 1, 2010
![Page 26: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/26.jpg)
25
van den Bergh (1986):
Nucleation fraction vs. galaxy magnitude
MB = -16 -14 -12 -10
Perhaps NSCs in fainter spheroids were destroyed by heating from the galaxy.
Thursday, July 1, 2010
![Page 27: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/27.jpg)
Orbits in Triaxial BH Nuclei(a) chaotic
(b-c) tubes
(g-i) pyramids
Poon & Merritt (2001)III. Triaxiality
Thursday, July 1, 2010
![Page 28: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/28.jpg)
γ = 1, T = 0.75 γ = 2, T = 0.75
ρ ∝ r−γ
Self-Consistent Triaxial NSCs
F = orbital fraction
T ≡ a2 − b2
a2 − c2
Pyramid
Tube
Chaotic Poon & Merritt (2004)
Thursday, July 1, 2010
![Page 29: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/29.jpg)
rGR � r � rinfl
Pyramid Orbits
Merritt & Vasiliev (2010)
• ~ Keplerian ellipses
• Librate about short axis
• Integrable (regular)*
• e ⇒ 1 at the corners! 1− e2
*for
Thursday, July 1, 2010
![Page 30: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/30.jpg)
2-body relaxation
Vector RR
Newtonian precession
Dynamica
l time
Scalar RR
Rel
ativ
istic
prec
essi
on
r, pc
T, yr
101101–
10–2
10–310
3
104
105
106
107
108
109
1010
Tdrain
=10–2
=10–3
=10–3
=10–2 precession
rbh
tpyr ≈ 3× 108yr
�T
10−2
��M•
106M⊙
�2/3
Pyramid Orbits
Merritt & Vasiliev (2010)
The time for pyramid orbits to “drain” is:
where T is the dimensionless coefficient of triaxiality.
The BH feeding rate can be much greater than that due to two-body relaxation.
Thursday, July 1, 2010
![Page 31: Dynamics of Nuclear Star Clusters · Dynamics of Nuclear Star Clusters David Merritt* ... Infall time scales for globular clusters are roughly correct 2. While they often contain](https://reader035.fdocumento.com/reader035/viewer/2022071420/611962f77935620968713f92/html5/thumbnails/31.jpg)
Summary
• Accretion of globular clusters appears to be a viable model for NSC formation, at least in bulge-dominated systems
• Disappearance of NSCs in spheroids with MB > -16 may be due to relaxation effects (heating from the galaxy)
• Stellar tidal disruption rates might be very high in triaxial NSCs (if SMBHs are present)
Thursday, July 1, 2010