Transcript Kuzmin

Kuzmin and Stellar Dynamics

Introduction

Dynamical models

G.G. Kuzmin’s pioneering work
– Mass models, orbits, distribution functions

Structure of triaxial galaxies

Conclusions
Galaxy Formation and Evolution

Galaxies form by hierarchical accretion/merging
–
–
–
–
–

Matter clumps through gravitation
Primordial gas starts forming first stars
Stars produce heavier elements (‘metals’)
Subsequent generations of stars contain more metals
Massive galaxies form from assembly of smaller units
Galaxy encounters still occur
– Deformation, stripping, merging
– Galaxies continue to evolve

Central black hole also influences evolution
Z=18
Z=0
Observational Approaches

Study very distant galaxies
– Observe evolution (far away = long ago)
– Objects faint and small: little information

Study nearby galaxies
–
–
–
–

Light not resolved in individual stars
Objects large & bright: structure accessible
Infer evolution through archaeology
Fossil record is cleanest in early-type galaxies
Study resolved stellar populations
– Ages, metallicities and motions of stars
– Archaeology of Milky Way and its neighbors
Dynamical Models

Aim: find phase-space distribution function f
–
–
–
–

Provides orbital structure
Mass-density distribution ρ = ∫∫∫ f d3v
Velocities v derive from gravitational potential V
Self-consistent model: 4πGρ= 2V
Approaches
– Assume f find ρ (but what to assume for f?)
– Assume ρ find f (solve integral equation)

Use Jeans theorem f = f(I) to make progress
– Provides f(E,L) for spheres, f(E,Lz) for axisymmetry
– f(E,I2,I3) for separable axisymmetric & triaxial models
Spheres

Hamilton-Jacobi equation separates in (r,θ,φ)
– Four integrals of motion: E, Lx, Ly, Lz
– All orbits regular: planar rosette’s

Mass model
– Defined by density profile ρ(r)

Gravitational potential by two single integrations

Selfconsistent models
– Isotropic models f=f(E) via Abel inversion (Eddington 1916)
– Circular orbit model: only orbits with zero radial action
– Many distribution functions: f=f(E), f=f(E+aL), f(E, L),
corresponding to different velocity anisotropies
– Constrain f further by measuring kinematics
Spheres

Large literature on construction of spherical models

Popular mass models include
– Hénon’s (1961) isochrone
– The -models (e.g., Dehnen 1993)



Already found by e.g., Franx in ~1988
Include the Jaffe (1982) and Hernquist (1990) models
Many of these were studied much earlier by Kuzmin
and collaborators
– In particular Veltmann (and later Tenjes)
– Density profiles and distribution functions
– Results not well known in Western literature, but
summarized in IAU 153, 363-366 (1993)
The Milky Way

Stellar motions near the Sun
– If Galaxy oblate and f=f(E, Lz) then vR2= vz2 and vRvz=0
– Observed: vR2 v2 vz2 and vRvz0
– Galactic potential must support a third integral of motion I3

Separable potentials known to have three exact
integrals of motion, E, I2 and I3, quadratic in velocities
– Stäckel (1890), Eddington (1915), Clark (1936)

Chandrasekhar assumed f=f(E+aI2+bI3) to find 
– This is the ‘Ellipsoidal Hypothesis’
– Model self-consistent only if  spherical: limited applicability

Little interest in opposite route: from  to f
– G.B.van Albada (1953): oblate separable potentials not
associated with sensible mass distributions ()
Kuzmin’s Contribution

Set of seminal papers based on his 1952 PhD thesis*

Considers mass models with potential
V 

F (  )  F ( )
 
in spheroidal coordinates (, , )
and F() a smooth function ( = , )

These potentials have
– Three exact integrals of motion E, Lz and I3
– Useful associated densities, given by simple formula
– (R, z)  0 if and only if (0, z)  0 (Kuzmin’s Theorem)
*Translated by Tenjes in 1996, including additions from 1969

Kuzmin’s Contribution

Assumption:

n=3
d
d
(   )
1/ 2
F ' ( )  
n / 2
– Fair approximation to Milky Way potential (no dark halo)
– Flattened generalisation of Hénon’s isochrone (1961)

n=4
– Exactly spheroidal model

 
0
(1  m )
2
2
with
In limit of extreme flattening
– Models  Kuzmin disk; surface density
– Rediscovered by Toomre (1963)

m  R  z /q

2
2
2
2
0
(1  R )
2
3/2
Model n=n0 is weighted sum of models with n>n0
– This built on his pioneering 1943 work on construction of
models by superposition of inhomogeneous spheroids
Kuzmin’s Contribution

Orbits in oblate separable models
– All short-axis tubes (bounded by coordinate surfaces)
– Similar to orbits in Milky Way found numerically by
Ollongren (1962) using Schmidt’s (1956) mass model

Distribution function f is function of single-valued
integrals of motion only
– Rediscovered by Lynden-Bell (1962)

f(E, Lz) for model n=3 (with Kutuzov, 1960)
– (R, z) can be written explicitly as (R, V) without any
reference to spheroidal coordinates
– Allows computing f(E, Lz) via series expansion à la Fricke
– f(E, Lz, I3) found by Dejonghe & de Zeeuw (1988) making
full use of the elegant properties of the model
Kuzmin 1972

Generalization of earlier work to triaxial shapes
– Very concise summary in Alma Ata conference 1972
– English translation in IAU 127, 553-556 (1987)

Potentials separable in ellipsoidal coordinates (,,)
–
–
–
–
–

Three exact integrals of motion E, I2 and I3
(x, y, z)  0 if and only if (0, 0, z)  0
Elegant formula for density
Includes ellipsoidal model:   (1 m ) with m  x
Four major orbit families
2
0
2
2
Rediscovered in 1982-1985 (de Zeeuw)
– Via completely independent route
2
 y / p  z /q
2
2
2
2
Separable Triaxial Models

Four orbit families
z
x
y
1. Box orbit

2. Inner longaxis tube orbit
3. Outer longaxis tube orbit
4. Short-axis
tube orbit
Same four orbit families found in Schwarschild’s
(1979) numerical model for stationary triaxial galaxy
Separable Triaxial Models

Mass models
–
–
–
–
Defined by short-axis density profile & central axis ratios
Stationary triaxial shape, with central core
Gravitational potential by two single integrations
Each model is weighted integral of constituent ellipsoids



Weight function follows via Stieltjes transform
Projection is same weighted integral of constituent elliptic disks:
new method for finding potential of disks
These properties shared by larger set of models
– Each ellipsoid
 
0
(1  m )
2
p
(p=n or n/2) generates similar family
de Zeeuw & Pfenniger (1988); Evans & de Zeeuw (1992)
Separable Triaxial Models

Jeans equations: obtain vi2 directly to ρ and V
– Three partial differential equations for three unknowns
– Equations written down by Lynden-Bell (1960), and solved
by van de Ven et al. (2003). No guarantee that f  0

Analytic selfconsistent models
– Thin-tube orbit models (only tubes with zero radial action)
– Existence of more than one major orbit family: f(E, I2, I3)
not uniquely defined by ρ(x, y, z)
– Abel models f = Σ fi(E+aiI2+biI3) Dejonghe; van de Ven et al. 2008

Through Kuzmin’s work and subsequent follow-up
the theory of stationary triaxial dynamical models is
now as comprehensive as that for spheres
Early-type Galaxies

Structure
– Mildly triaxial shape
– Central cusp in density profile
– Super-massive central black hole

Implications for orbital structure
– No global extra integrals I2 and I3
– Three tube orbit families
– Box orbits replaced by mix of boxlets (higher-order
resonant orbits) and chaotic orbits: slow evolution

Dynamical models
– Construct by numerical orbit superposition
– Use separable models for testing and insight
– Use kinematic data to constrain f
Stellar Orbits in Galaxies
T=1
T=10
T=50
T=200
Image of orbit on sky

Galaxies are made of stars
 Stars move on orbits (with integrals of motion)
 Galaxies are collections of orbits
Schwarzschild’s Approach
Images of model orbits
Observed
galaxy image

Many different orbits possible in a given galaxy

Find combination of orbits that are occupied by
stars in the galaxy  dynamical model (i.e. f)
Schwarzschild 1979; Vandervoort 1984
Numerical Orbit Superposition

No restriction on form of potential
– Arbitrary geometry
– Multiple components (BH, stars, dark halo)

No restriction on distribution function
– No need to know analytic integrals of motion
– Full range of velocity anisotropy

Include all kinematic observables
– Fit on sky plane
– Codes exist to do this for spherical, axisymmetric and
non-tumbling triaxial geometry
Leiden group: Cretton, Cappellari, van den Bosch; Gebhardt & Richstone; Valluri
The E3 Galaxy NGC 4365

Kinematically Decoupled Core
– Long-axis rotator, core rotates
around short axis (Surma & Bender 1995)

SAURON kinematics:
– Rotation axes of main body and
core misaligned by 82o
– Consistent with triaxial shape, both
long-axis & short-axis tubes occupied

Customary interpretation:
– Core is distinct, and remnant of
last major accretion ~12 Gyr ago
Triaxial Dynamical Model

Parameters
– Two axis ratios, two viewing
angles, M/L, MBH

Best-fit model
– Fairly oblate (0.7:0.95:1)
– Short axis tubes dominate,
but ~50% counter rotate,
except in core; cf NGC4550
– Net rotation caused by
long-axis tubes, except in core
– KDC not a physical subunit,
but appears so because of
embedded counter-rotating structure
van den Bosch et al. 2008
Dynamics of Slow Rotators

11 slow rotators in representative SAURON sample
– Range of triaxiality: 0.2  T  0.7  no prolate objects
– Mildly radially anisotropic
– Most have ‘KDC’

Dynamical structure
– Short axis tubes dominate
– Smooth variation with radius
– ~similar to dry merger simulations
Jesseit et al. 2005; Hoffman et al. 2010
– No sudden transition at RKDC

KDC not distinct from main body
– In harmony with smooth Mgb and Fe gradients
van den Bosch et al. 2011, in prep.
Conclusions

Kuzmin was a very gifted dynamicist

Much of this work was unknown in West
– Few read Russian; translations came later, but even today
most papers are not in, e.g., ADS
– Kuzmin sent short English synopses to key dynamicists,
but these were not widely distributed
– Perek’s (1962) review did help advertize the results, but
even so, much of his work was independently rediscovered

Kuzmin’s work has substantially increased our
understanding of galaxy dynamics

And increased the luminosity of Tartu Observatory