Lecture 24 Stability of Molecular Clouds

Similar documents
Lecture 22 Stability of Molecular Clouds

General momentum equation

Section 11. Timescales Radiation transport in stars

Today in Astronomy 142: the Milky Way s disk

But for simplicity, we ll define significant as the time it takes a star to lose all memory of its original trajectory, i.e.,

A thermodynamic degree of freedom solution to the galaxy cluster problem of MOND. Abstract

Astrophysical Fluid Dynamics Solution Set 5 By Eric Bellm, Jeff Silverman, and Eugene Chiang

Physics 235 Chapter 5. Chapter 5 Gravitation

Diffusion and Transport. 10. Friction and the Langevin Equation. Langevin Equation. f d. f ext. f () t f () t. Then Newton s second law is ma f f f t.

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 4: Toroidal Equilibrium and Radial Pressure Balance

PROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.

Determining solar characteristics using planetary data

Review: Electrostatics and Magnetostatics

Stellar Structure and Evolution

Chem 453/544 Fall /08/03. Exam #1 Solutions

Takuya Ohtani (Osaka University Theoretical Astrophysics Group D2) Collaborator: Toru Tsuribe(Osaka Univ.)

Black Body Radiation and Radiometric Parameters:

From Gravitational Collapse to Black Holes

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s

= 4 3 π( m) 3 (5480 kg m 3 ) = kg.

PROBLEM SET #3A. A = Ω 2r 2 2 Ω 1r 2 1 r2 2 r2 1

Magnetic Dipoles Challenge Problem Solutions

11) A thin, uniform rod of mass M is supported by two vertical strings, as shown below.

OSCILLATIONS AND GRAVITATION

working pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50

On the Sun s Electric-Field

The Last 40 Years of Accretion Disk Theory...(abridged)

F-IF Logistic Growth Model, Abstract Version

Classical Mechanics Homework set 7, due Nov 8th: Solutions

Physics 221 Lecture 41 Nonlinear Absorption and Refraction

Chapter 7-8 Rotational Motion

Right-handed screw dislocation in an isotropic solid

Gauss Law. Physics 231 Lecture 2-1

is the instantaneous position vector of any grid point or fluid

Physics 505 Homework No. 9 Solutions S9-1

DOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS

EXAM NMR (8N090) November , am

In the previous section we considered problems where the

Velocimetry Techniques and Instrumentation

Gaia s Place in Space

arxiv: v1 [physics.pop-ph] 3 Jun 2013

PES 3950/PHYS 6950: Homework Assignment 6

EM Boundary Value Problems

DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS IB PHYSICS

Downloaded from

15 ρ R Total angular momentum of the Jupiter-Sun System. momentum plus the spin momenta of Jupiter and the Sun.

Physics 2212 GH Quiz #2 Solutions Spring 2016

Liquid gas interface under hydrostatic pressure

Between any two masses, there exists a mutual attractive force.

Homework 7 Solutions

INFLUENCE OF GROUND INHOMOGENEITY ON WIND INDUCED GROUND VIBRATIONS. Abstract

Quantum Mechanics and Stellar Spectroscopy

15 Solving the Laplace equation by Fourier method

Radial Inflow Experiment:GFD III

Problems with Mannheim s conformal gravity program

Pressure Calculation of a Constant Density Star in the Dynamic Theory of Gravity

B. Spherical Wave Propagation

Phases of Matter. Since liquids and gases are able to flow, they are called fluids. Compressible? Able to Flow? shape?

Hydroelastic Analysis of a 1900 TEU Container Ship Using Finite Element and Boundary Element Methods

Magnetohydrodynamics (MHD) I

Entropy and Free Energy: Predicting the direction of spontaneous change The approach to Chemical equilibrium

Three dimensional flow analysis in Axial Flow Compressors

Newton s Laws, Kepler s Laws, and Planetary Orbits

Chapter 13 Gravitation

MODULE 5 ADVANCED MECHANICS GRAVITATIONAL FIELD: MOTION OF PLANETS AND SATELLITES VISUAL PHYSICS ONLINE

TheWaveandHelmholtzEquations

AP-C WEP. h. Students should be able to recognize and solve problems that call for application both of conservation of energy and Newton s Laws.

, and the curve BC is symmetrical. Find also the horizontal force in x-direction on one side of the body. h C

The geometric construction of Ewald sphere and Bragg condition:

PHYSICS NOTES GRAVITATION

1 Fundamental Solutions to the Wave Equation

1) Consider an object of a parabolic shape with rotational symmetry z

ECE 3318 Applied Electricity and Magnetism. Spring Prof. David R. Jackson ECE Dept. Notes 13

Rigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018

Preliminary Exam: Quantum Physics 1/14/2011, 9:00-3:00

d 2 x 0a d d =0. Relative to an arbitrary (accelerating frame) specified by x a = x a (x 0b ), the latter becomes: d 2 x a d 2 + a dx b dx c

Numerical solution of diffusion mass transfer model in adsorption systems. Prof. Nina Paula Gonçalves Salau, D.Sc.

Physics 2A Chapter 10 - Moment of Inertia Fall 2018

Pulse Neutron Neutron (PNN) tool logging for porosity Some theoretical aspects

Fresnel Diffraction. monchromatic light source

Ch 30 - Sources of Magnetic Field! The Biot-Savart Law! = k m. r 2. Example 1! Example 2!

Substances that are liquids or solids under ordinary conditions may also exist as gases. These are often referred to as vapors.

3.1 Random variables

Lecture 8 - Gauss s Law

Nuclear and Particle Physics - Lecture 20 The shell model

Physics 107 TUTORIAL ASSIGNMENT #8

10. Groundwater in geotechnical problems

Introduction to Arrays

Algebra-based Physics II

Phys-272 Lecture 17. Motional Electromotive Force (emf) Induced Electric Fields Displacement Currents Maxwell s Equations

Introduction to Dielectric Properties and Magnetism

Astronomy 111, Fall October 2011

Chapter 3 Optical Systems with Annular Pupils

17.1 Electric Potential Energy. Equipotential Lines. PE = energy associated with an arrangement of objects that exert forces on each other

A New Approach to General Relativity

General Railgun Function

Title :THERMAL TRANSFER AND FLUID MECHANICS IN THE THEORY OF ETHER Author:Thierry DELORT Date:1 st May 2013

To Feel a Force Chapter 7 Static equilibrium - torque and friction

Hopefully Helpful Hints for Gauss s Law

EFFECTS OF FRINGING FIELDS ON SINGLE PARTICLE DYNAMICS. M. Bassetti and C. Biscari INFN-LNF, CP 13, Frascati (RM), Italy

Transcription:

Lectue 4 Stability of Molecula Clouds 1. Stability of Cloud Coes. Collapse and Fagmentation of Clouds 3. Applying the iial Theoem Refeences Oigins of Stas & Planetay Systems eds. Lada & Kylafis http://cfa-www.havad.edu/cete Myes, Physical Conditions in Molecula Clouds Lada, The Fomation of Low-Mass Stas Shu, Ch. 1 & Stahle, Appendix D fo iial Theoem

1. Stability of Molecula Cloud Coes Stat with the molecula cloud coe popeties fom Lec. 3 and focus on the items elating to viial equilibium. 1. associated with sta fomation. elongated (aspect atio ~ :1) 3. dynamics may be dominated by themal o tubulent motion 4. tempeatue: T ~ 10 0 K 5. size: R ~ 0.1 pc 6. ionization: x e ~ 10-7 7. size-line width elation R ~ σ p, p = 0.5 +/- 0. 8. appoximate viial equilibium 9. mass spectum: simila to GMCs.

Coe Coelations (Myes 1999) v /R non-themal vs. themal line widths consistency with viial equilibium Fo pue themal suppot, GM/R ~ kt/m, m =.3 m H ), and R 0.1 (M / M ) (10 K / T) pc. N

Coe Stability iial equilibium seems an appopiate state fom which coes poceed to make stas. Howeve, as a significant faction of coes ae obseved to have embedded potostas (obseved by IRAS), they can t be completely quiescent. We ae theefoe concened with the issue of coe stability, e.g., why they ae stable (if they ae) and how they become de-stabilized and collapse unde gavity to fom stas. Although Mye s data pesentation suggests that andom motion (themal and/o non-themal) may stabilize coes against gavitational collapse, we should fist conside othe possibilities, in paticula otation and magnetic fields.

Effects of Rotation Molecula cloud coes have modest velocity gadients, 0.4-3 km s -1 pc -1 and angula speeds Ω ~ 10-14 -10-13 ad s -1. β = E ot / E gav ~ 0.0 gad v β Coes do not appea to be suppoted by otation. J/M NB Not all of the obseved gadients coespond to the R oveall otation of the coe Goodman et al. ApJ 406 58 1993

Magnetic Pessue Magnetic fields ae difficult to measue in cloud coes OH absoption measuements ae had because of the small chance of a backgound adio souce Example of dak cloud B1 B cos θ -19 ± 4 µg OH pobes n H 10 3 cm -3 B /8π ~ 3 x 10 5 K cm -3 Compaable to the themal pessue of a 10 K coe with n ~ 10 4 cm -3 1667 MHz 1665 MHz Stokes Stokes I Cutche 1993 ApJ 407 175 Recall Lectue 13 and adio measuements of the Zeeman effect.

. Collapse and Fagmentation How do coes condense out of the lowe density egions of GMCs? The conventional wisdom is local gavitational instabilities in a globally stable GMC via the classical Jeans instability. Conside a unifom isothemal gas in hydostatic equilibium with gavity balanced by the pessue gadient. Now suppose a small spheical egion of size is petubed ρ Xρ whee X > 1 The ove-density Xρ geneates an outwad pessue foce pe unit mass:

foce pe unit mass Gavity wins out if Simple Stability Citeion The inceased density also leads > c Gρ 0 o F ~ p ~ GM c Gρ 0 / ρ ~ Xc / The ight hand side is essentially the Jeans length: λ = c J F G p > π Gρ 0 / to an inwad gavtitational ~ GXρ. 0 ρ 0 Xρ 0

we get a The Jeans Length & Mass Fom the Jeans length whose numeical value is J coesponding Jeans mass M J M λ = J = λ = 7.5M 3 J sun πc Gρ 0 πc ρ0 = G 3/ ρ 1/ 0 ( ) ( 4-3 T /10K 3/ n(h ) /10 cm ) c 3 ρ 0 1/ Length scales > λ J o masses > M J ae unstable against collapse.

Relevance of the Jeans Analysis Many assumptions have been made implicitly in the above analysis that violate the basic fact that molecula clouds ae not static but ae tubulent on lage scales. Tubulence povides an effective pessue suppot, and its supesonic motions must be dissipated by shocks befoe the collapse can poceed: Is it enough to eplace the sound velocity c in the Jeans fomulae by c = m kt +σ tub As the gas condenses the Jeans mass deceases, which suggests fagmentation (Hoyle 1953 ApJ 118 513), a collapse cascade into smalle and smalle masses.

Does Fagmentation Occu The dispesion elation fo a petubation δρ ~ ei(ωt - kx) in the Jeans poblem is ω = c k -4πGρ 0 o ω = c (k - k J ), k J = 4π Gρ 0 / c k - k J < 0 makes ω imaginay: Exponential gowth occus fo k < k J Gowth ate -iω inceases monotonically with deceasing k, which implies that Longest wavelength petubations (lagest mass) gow the fastest The fast collapse of the lagest scales suggests that fagmentation is unlikely (Lason 1985 MNRAS 14 379)

Swindled by Jeans? Conside an altenate model, a thin sheet of suface density. The dispesion elation is ω = c k -π G k o ω = c ( k - k c k ) Exponential gowth occus fo k < k c = π G / c with gowth ate -i ω =( π G k - c k ) 1/ which is maximum at k f = k c /= π G / c and gives a pefeed mass, M f ~ (π / k f ) = 4 c 4 / G This is a diffeent esult than the pevious Jeans analysis fo an isotopic 3-d medium. The maximum gowth ate now occus on an intemediate scale that is smalle than pedicted by Jeans, possibly favoing fagmentation.

Pefeed Length & Mass Apply the thin sheet model to the Tauus dak cloud. T = 10 K, c = 0.19 km/s λ c ~ 0.05 pc A() 5 mag. o 0.03 g cm - λ f ~ 0.1 pc, M f ~ M t f = λ f /(π c) ~ x 10 5 y Moe ealistic analyses tend to show that, when themal pessue povides the dominant suppot against gavity, thee is a minimum length & mass scale which can gow If the cloud is non-unifom thee is a pefeed scale Typically ~ few times the chaacteistic length of the backgound, e.g., the scale height, H =c /πg, fo a gaseous equilibium sheet (Lason 1985 MNRAS 14 379)

3.The iial Theoem and Stability A. Geneal Consideations The deivation of viial theoem stats with the equation of motion fo the macoscopic velocity of the fluid afte aveaging ove the andom themal velocities. The mathematics is given by Shu Ch. 4 and Stahle Appendix D. Dv t ρ = P ρ φ + T Dt Tij = Bi B j / 4π B δij / 8π t 1 T B ( 1 = + B )B 8π 4π The pespective of the iial Theoem is that these equations contain too much infomation. Hence, take moments, specifically the fist moment:

iial Theoem: LHS Analysis Take the dot poduct with and integate the equation of motion ove a finite volume, and analyze the ight and left hand sides sepaately. The LHS side yields two tems since: ( ) + = = Dt D d v v d Dt D d and ) ( ) ( 1 ρ ρ ρ K E Dt I D d v d Dt D 1 1 and finally = ρ ρ The 1 st tem vanishes fo a static cloud.

iial Theoem: RHS Analysis The RHS yields volume & suface tems 3P d + S P + B 8π B d 8π ds + 1 4π φ ( ) S B dm B d ( )( S ) The 3 d tem becomes the gavitational enegy. The 1 st and 4 th tems measue the diffeence between the extenal and mean intenal pessues. The othe tems ae tansfomed into magnetic pessue and tension.

The iial Theoem ( ) and then inset appoximate expessions in the RHS. 3 1 3 1 3 pessue we can solve fo the extenal Fo a static cloud, 8 1 ) ( 4 1 8 whee 3 K ext ext K M W E P P B ds B B ds B d M M W P P E Dt I D + + = + = + + + = π π π P ex t becomes a function of the global vaiables of the volume and the gas inside. If thee is no otation, E K = 0.

Simplest Application of the iial Theoem Conside a spheical, isothemal, un-magnetized & nontubulent cloud in equilibium of mass M and adius R 3 3 4πR Pext = 3Mc 5 GM R whee we have used P = ρc. If we now also ignoe gavity, we check that the extenal and intenal pessues ae the same. In the pesence of gavity, the second tem seves to lowe the extenal pessue needed to confine the gas Inside volume. Next, divide by 4πR 3 and plot the RHS vs. R

Simplest Application of the iial Theoem P ext P = 3c M 4π R c 1 R 3 3GM 0π Thee is a minimum adius below which this model cloud cannot suppot itself against gavity. States along the left segment, whee P inceases with R, ae unstable. 1 R 4 P ~ 1/R 3 R

Stable and Unstable iial Equlibia Fo a given extenal pessue, thee ae two equilibia, but only one is stable. P A B P ex t R max R Squeeze the cloud at B (decease R) Requies moe pessue to confine it: e-expands (stable) Squeeze the cloud at A (decease R) Requies less pessue to confine it: contacts (unstable)

The Citical Pessue Diffeentiating the extenal pessue fo fixed M and c P ext = 3c M 4π gives the citical adius, R c 1 R 3 GM c 3GM 0π 1 R The coesponding density and mass ae 6 3 c c ρc M 3 c M This last esult is essentially the same as gotten fom the Jeans analysis. ρ c 4

Fagmentation vs. Stability (o Jeans vs. iial) Both methods yield a simila chaacteistic Jeans mass The diffeence is in the initial conditions Cloud scenaio assumes a self-gavitating cloud coe in hydostatic equilibium The cloud may o may not be close to the citical condition fo collapse In the gavitational fagmentation pictue thee is no cloud It cannot be distinguished fom the backgound until it has stated to collapse Both pespectives may be applicable in diffeent pats of GMCs Thee ae pessue confined clumps in GMCs Not stongly self-gavitating Coes make stas Must be self-gavitating If fagmentation poduces coes these must be collapsing and it is too late to apply viial equilibium

Fagmentation vs. Stability A vey elevant question is how quickly coes fom (elative to the dynamical time scale) NH 3 coes do seem to be close to viial equilibium but they could be contacting (o expanding) slowly Statistics of coes with and without stas should eveal thee ages If coes ae stable, thee should be many moe coes without stas than with young stas Compaison of NH 3 coes and IRAS souces suggests ~ 1/4 of coes have young stella objects (Wood et al. ApJS 95 457 1994), but Jijina et al. give a lage faction The coesponding life time is < 1 My (Onishi et al. 1996 ApJ 465 815)