Lecture 5 By Tom Wilson 1 Dust clouds
Lecture 5 By Tom Wilson 1
Dust clouds in Taurus Regions shown in gray are those where very few stars are seen (Becvar atlas) 2
The Discovery of Molecular Clouds -The discussion began in 1960’s when it was noticed that dark clouds had a deficiency of HI emission, based on an Av to HI relation obtained from studies of diffuse clouds. Some thought this could be due to high optical depth in HI, or self absorption, etc. Barrett et al. found intense OH lines 3
HI in Dark Clouds The dashed line is the boundary of the visible cloud 4
HI Self Absorption in Dark Clouds The broader HI emission is from a warm background 5
Molecular Clouds n n n The discovery of thermal OH emission from these clouds, and then NH 3, H 2 O and H 2 CO were found in warmer clouds. This made it likely that the gas in these clouds was mostly H 2 The final piece of evidence was the discovery of CO Still were problems with the abundance (even the very presence) of complex molecules. 6
The Life Cycle of H 2 -Make this molecule on grains, since must conserve energy and momentum -The exact production rate is known to a factor of 2 -Once made, H 2 enters the gas phase - H 2 is dissociated by 11 e. V photons - H 2 is destroyed by spectral line radiation, so is self shielded, and also shielded by dust - H 2 combines with H+ to form H 3+ which promotes ion-molecule chemistry 7
TMC-1: A special dust cloud in Taurus Location of TMC-1 8
Orion Hot Core A “chemistry factory” in the sky Trapezium stars, which ionize the HII region Bar, an Edge-on Photon Dominated Region Photo From HST (O’Dell) The molecular cloud in Orion is extended north-south, behind the HII region 9
Structures of a few interstellar molecules Methyl alcohol These are asymmetric top molecules which have Methyl formate many transitions in the millimeter range in warm sources Dimethyl ether 10
CLASSIFICATIONS OF PROTOSTARS m mm m 11
Sketch of a Photon Dominated Region From Hollenbach & Tielens 1999 12
Lecture 5 page 2 For molecules, the electronic states are given analogous assignments: S S P D 2 electron systems 3 S = “triplet S” parallel electron spins 1 S = “singlet S” auto-parallel electron spins 13
MOLECULES Vastly more complex than hydrogen atoms. The main assumption in a discussion of molecules is Born-Oppenheimer approximation: WEL (electronic energies): transitions with energies of a few e. V WVIB (vibrational energies): transitions have energies of 0. 1 to 0. 01 e. V WROT (rotational energies): transitions have energies of 0. 001 e. V To first order there no mixed terms in Hamiltonian, this means that the wavefunctions are product of: Y(electronic) * Y(vibration) * Y(rotation) Another way to say this, is that the electrons move much fasten than the rotation of the nuclei. Thus have a negative bowl of potential energy from electrons. In this, the nuclei rotate. What is the energy equivalent of 1 e. V in degrees Kelvin? 14
Morse Potential P(r) r 15
Approximations for P(r) Often can use 1 st order expansion Since for most radio lines, one is deep in potential well ROTATIONAL SPECTRA OF DIATOMIC MOLECULES The Hamiltonian is ╫ ╫ 16
Demand Assumed perfectly rigid ╫ ╫ ╫ Calculate the Be for 12 C 16 O and 13 C 16 O using the nuclear masses in Atomic mass units (AMU), rounded to integers. Compare the values and calculate the line frequencies using hn=EROT 17
Effects of distortion 18
So This is correct if molecule is perfectly rigid. If stretching, have an extra term Where D is > 0 and D / B =10 -5 EXAMPLE: For 12 C 16 O, Be = 57. 6360 GHz and De = 0. 185 MHz Use this relation given above to calculate the J = 1 – 0, 2 – 1, 3 – 2 line frequencies 19
For linear molecules, find From before And For total population of CO, under LTE, get Work out the Einstein A coefficients for CS, Si. O, HCO+, N 2 H+. Then calculate the critical density, n*, where A = n* < s v > = 10 -10 n*. Compare to n* for CO J = 1 → 0 20
SYMMETRIC AND ASYMMETRIC TOPS is constant, and projection on ax axis 1. - SYMMETRIC TOPS – projection on molecule axis of symmetry is also a constant. Definition: ╫ ╫ ╫ So ╫ ╫ ╫ 21
CH 3 CN: Prolate symmetric Top NH 3: Oblate symmetric Top Usually one defines ╫ ╫ Then for energy eigenvalues And the energy above the ground state is 22
Ammonia Quantum numbers are (J, K), where J is total angular momentum, K is component along molecule axis 23
Note that So 0 → 0 is forbidden, but An extra feature of NH 3 is inversion doubling. This is caused by the Q. M. tunneling of N through the plane of the 3 hydrogen atoms 24
The inversion doubling gives rise to many lines in the centimeter range from 20 GHz to 40 GHz. These are (J, K) → (J, K) lines Ammonia is an Oblate symmetric top. CH 3 CN and CH 3 C 2 H are Prolate symmetric tops. For prolate, “C” becomes “A”. Dipole moment along axis of symmetry of molecule. So no radiative transition involving a change in K quantum number 2. - ASYMMETRIC TOPS Most common type of polyatomic molecule: H 2 CO, H 2 O, … Now, only is a constant. Describe energy levels by the approximation that the asymmetric top is somewhere between a PROLATE and OBLATE symmetric top. That is between K values Ka and Kc. Levels are described by the quantum numbers Each of the K values is less than or equal to the J value 25
prolate oblate 26
Formaldehyde Almost a symmetric top molecule, with dipole moment along A axis. In an allowed radiative transition, no charge in Ka. There will be a change in J an Kc. So 101 → 000 (n = 72 GHz) Can also have 110 → 111 or 110 ← 111 (6 cm, or 4. 8 GHz) H 2 CO is nearly a prolate symmetric top H 2 O very asymmetric top, with dipole moment along B axis 27
Almost like a prolate symmetric top 28
EXAMPLE: H 2 CO A rotation about the A axis interchanges the 2 hydrogens. This must have an anti-symmetric wavefunction since these are FERMI particles. So Y(space) * Y(spin) = anti symmetric Definite relation of spin and (-1)Ka : if spins are parallel, ↑↑, then is ortho. If so, the parity of the space wavefunction must be anti-symmetric ↑↓ If spin is para, then the space wavefunction must be symmetric A rotation about B, A, and C, must leave the molecule in the initial state. So the product of parities of the 3 rotations must be symmetric. Then for ortho: the space part is anti symmetric in B since (-1)Ka is the symmetry. Then have Ka Kc = (even)(odd). Transitions for ortho: Since J must change in an allowed radiation transition, must connect in upper and lower states, Kc must change between states. Ka Kc=(e)(o)-(e)(e)-(e)(o) For para: the space symmetry of A is even, so wavefunctions must be odd. So para transitions for H 2 CO are (o)(o)-(o)(e) 29 (o)(e)-(o)(o)
Allowed radiative transitions are dipole. Thus parity of initial and final states must be different: For i and f states the parity must be different r: over all space r changes sign, but product be positive, so the partiy of the initial and final wavefunctions must be opposite. Symmetry of spatial part of wavefunctions is crucial: look at H 2 CO If proton spins parallel, have Ortho: Ortho the spin wavefunction is symmetric. Rotation of 180° about A axis interchanges protons. These are FERMI particles as total symmetry must be odd. Thus since the spin wavefunction is even so space wavefunction must be odd. Since the spin wavefunction varies as (-1)Ka. So when Ka = 0, the spin wavefunction is symmetric. When the spins ↑↓ this is a para state; the Ka = 1. For parallel spins, ↑↑, ortho states, Ka must be 0, 2, 4, 6, etc… 30
EXAMPLE: H 2 O A rotation about the B axis interchanges the 2 hydrogens. This must have an anti-symmetric wavefunction since these are FERMI particles. So Y(space) * Y(spin) = anti symmetric Definite relation of spin and (-1)Kb : ↑↑ ortho space: anti-symmetric ↑↓ para space: symmetric A rotation about B, A, and C, leave the molecule in the initial state. So the parities of the 3 rotations must be symmetric then ortho: ortho space anti symmetric in B since (-1)Ka is the symmetry have Ka Kc = (even)(odd) or (odd)(even) Ortho: Ortho Since J must change in an allowed radiation transition, must connect in upper and lower states (o)(e)-(e)(o)-(o)(e) Para: Para the symmetry of B is even, so A and C wavefunctions must be both even or both odd. So para transitions for H 2 O are (o)(o)-(e)(e)-(o)(o) 31
32
Methanol Not just more complex n n Has different notation The CH 3 and OH arm can rotate relative to each other, giving rise to ‘internal rotation’ 33
J, Total Angular momentum Angular Momentum Along Axis of Symmetry 34
OH Molecule Sketch is to show Not shown are The effect of L The Far IR Doubling Transitions connecting 84 micron lines rotationally excited states to ground 120 micron lines 35
How convert line intensity into a column density? Use t General Relation T is Tex Emission lines: TL=Tex(1 -e-t). If t<<1, TL=Tex t 36 36
CO: Ammonia: In General: Asymmetric Tops: T&S Table IV or Journal of Physical and Chemical Ref. Data for S(u, b) 37
LVG MODEL Assumes close coupling of molecule and radiation field. Also This makes the radiative transfer a local problem. Qualitative result: Instead of an escape from cloud, after a collision the photon is reabsorbed many times. Only happens if t > 1 38
39
The story of X The clumps are more dense. Each is assumes to be Virially stable, and then entire cloud is stable Virially. From LVG, Virial theorem: So 40
Need more data to get a secure total column density 1. Could measure lots of lines of a species, make “ Boltzmann Plot” 2. For CO, a number of models of excitation. From Mauersberger, one is based on large velocity gradient model Another is N(H 2): based on “clumpy cloud” model and a virialized cloud 41
THE END ! 42
- Slides: 42