Computational Laboratory Astrophysics for Star Formation: getting the microphysics right Phillip C. Stancil Department of Physics and Astronomy and the Center for Simulational Physics University of Georgia
Introduction: Astro-photos to Astro-physics Input data Model Intensity 6 10 14 Wavelength Output
Introduction: Astro-photos to Astro-physics Input Garbage data Model Intensity 6 10 14 Wavelength Output??
Outline Introduction Making Molecules: sticking on grains Exciting Molecules: molecular collisions Breaking Molecules: photodissociation Summary
Collaborators and Funding Junko Takahashi Benhui Yang Vijay Veeraghattam Samantha Fonseca dos Santos N. Balakrishnan Katie Manrodt Bob Forrey Joel Bowman Steve Lewis
Sticking Coefficients on Ice Mantles: H2 1 10 K amorphous ice Sticking coefficient 0.8 0.6 0.4 Ω 2 = 2 Ω 2 = 1 100 K 3 K Hollenbach & Salpeter (1970) Leitch-Devlin & Williams (1985) Current Work 10 K ice LAMMPS Molecular Dynamics code TIP4P Water-Water Potential H2-H2O potential of Zhang et al. (1992) 0.2 H2(v=0,J=0) 0 0 50 100 150 200 250 300 350 400 H 2 Kinetic Energy (K) S 1 Veeraghattam et al. 2013a, in prep.
Sticking Coefficients on Ice Mantles: benchmarking experiment Sticking coefficient 0.8 0.7 0.6 0.5 0.4 Current Work Matar et. al (2010) 10 K amorphous ice Experiment: H2 beam at 62 o from surface normal 0.3 0.2 0 50 100 150 200 250 300 350 400 H 2 Kinetic Energy (K) Veeraghattam et al. 2013a, in prep.
Sticking Coefficients on Ice Mantles: atomic H Early days Sticking coefficient 1 0.8 0.6 0.4 0.2 MD: Buch & Zhang 1991 QM: Hollenbach & Salpeter 1970 (Ω 2 =1) QM: Hollenbach & Salpeter 1970 (Ω 2 =2) 10 K amorphous ice Hollenbach & Salpeter used in most models MD = molecular dynamics Buch & Zhang considered small cluster 0 0 100 200 300 400 500 600 H atom kinetic energy (K)
Sticking Coefficients on Ice Mantles: atomic H 1990s-2000s Sticking coefficient 1 0.8 0.6 0.4 0.2 MD: Al-Halabi & van Dishoeck 2007 MD: Masuda et al. 1998 MD: Buch & Zhang 1991 QM: Hollenbach & Salpeter 1970 (Ω 2 =1) QM: Hollenbach & Salpeter 1970 (Ω 2 =2) Exp: Manico et al. 2001 10 K amorphous ice More MD simulations H2 recombination experiment 0 0 100 200 300 400 500 600 H atom kinetic energy (K)
Sticking Coefficients on Ice Mantles: atomic H Today Sticking coefficient 1 0.8 0.6 0.4 0.2 MD: Al-Halabi & van Dishoeck 2007 MD: Masuda et al. 1998 MD: Buch & Zhang 1991 QM: Hollenbach & Salpeter 1970 (Ω 2 =1) QM: Hollenbach & Salpeter 1970 (Ω 2 =2) MD: Current Work Exp: Manico et al. 2001 Exp: Watanabe et al. 2010 10 K amorphous ice TIP4P water-water potential. H-H2O potential of Zhang et al. (1991) MD code of Masuda et al. (1998) S=1, below 100 K 0 0 100 200 300 400 500 600 H atom kinetic energy (K) Veeraghattam et al. 2013b, to be submitted
Non-LTE Molecular Spectra: Role of Collisional Excitation H2(J1) + CO(J2) Critical density (cm -3 ) 10 7 10 6 10 5 10 4 10 3 10 2 CO(v=0,J2) J 2 =1 J 2 =40 10 100 1000 Temperature (K) Density is below critical density in most environments Upper levels not thermalized Level populations depend on collisional rates Yang et al. 2010, ApJ, 718, 1062
J =10 Pure Rotational De-excitation: H2(v1=0,J1=1) + CO(v2=0,J2=20) Rate coefficient (cm 3 /s) 10-10 10-11 Symbols: Flower 2001 17 16 18 J2 =10 J2 =19 Rigid-rotor approximation Close-coupling calculations with new surface (JS04) J2=1-40, paraand ortho-h2 10-12 1 10 100 1000 Temperature (K) 1-3000 K Yang et al. 2010, ApJ, 718, 1062
Pure Rotational De-excitation: sensitivity to potential energy surfaces H + CO(v=0,J=1) Rigid-rotor Rate coefficient (cm 3 s -1 ) 10-10 10-11 1->0 1975, Chu & Dalgarno Close-coupling calculations Semi-empirical potential surface 1 10 100 1000 Temperature (K)
Pure Rotational De-excitation: sensitivity to potential energy surfaces H + CO(v=0,J=1) Rigid-rotor Rate coefficient (cm 3 s -1 ) 10-10 10-11 1->0 1975, Chu & Dalgarno 1976, Green & Thaddeus Close-coupling calculations Alternate semiempirical potential surface 1 10 100 1000 Temperature (K)
Pure Rotational De-excitation: sensitivity to potential energy surfaces H + CO(v=0,J=1) 2002, Balakrishnan et al. Rigid-rotor Rate coefficient (cm 3 s -1 ) 10-10 10-11 1->0 1975, Chu & Dalgarno 1976, Green & Thaddeus Close-coupling calculations New ab initio potential surface (WKS) 1 10 100 1000 Temperature (K)
Pure Rotational Excitation: sensitivity to potential energy surfaces H + CO(v=0,J=0) Cross section (10-16 cm 2 ) 20 15 10 5 400 cm -1 PES: J=0 J MRCI CCSD(T) BBH WKS Homonuclear: forbidden odd ΔJ transitions Near-homonuclear: odd ΔJ transitions suppressed (Lee & Bowman 1987) New H-CO results violated this trend 0 0 1 2 3 4 5 6 7 8 9 10 11 J' Shepler et al. 2007, A&A, 475, L15
Pure Rotational De-excitation: sensitivity to potential energy surfaces H + CO(v=0,J=1) Rate coefficient (cm 3 s -1 ) 10-10 10-11 CCSD(T) MRCI 2002, Balakrishnan et al. 1975, Chu & Dalgarno 2013 1976, Green & Thaddeus 1->0 Rigid-rotor Close-coupling calculations J=1-5 1 10 100 1000 Temperature (K) T=1-3000 K Yang et al. 2013a, ApJ, submitted
Pure Rotational De-excitation: 10-10 He + HF(v=0,J=6) Full 3D Rate coefficient (cm 3 s -1 ) 10-11 10-12 10-13 10-14 J'=5 4 3 2 1 0 10-15 1 10 100 1000 Temperature (K) Close-coupling calculations v=0-1, range of J T=1-3000 K Yang et al. 2013b, in prep.
Pure Rotational De-excitation: effect of vibrational state 10-10 He + HF(v=0,1;J=6) Full 3D Rate coefficient (cm 3 s -1 ) 10-11 10-12 10-13 10-14 J'=5 4 3 2 1 0 v=0 v=1 10-15 1 10 100 1000 Temperature (K) Close-coupling calculations v=0-1, range of J T=1-3000 K Yang et al. 2013b, in prep.
Pure Rotational De-excitation: effect of vibrational state He + HF(v=0,1;J=6) Rate coefficient ratio (v=1 / v=0) 3.0 2.5 2.0 1.5 J'=5 4 2 3 1 0 Scaling via exponential energy gap law 1.0 1 10 100 1000 Temperature (K) Yang et al. 2013b, in prep.
Rovibrational Transitions for Diatom-Diatom: Collisions in 6D H2(v1J1) + H2(v2J2) Rate Coefficients (cm 3 s -1 ) 10-13 10-14 10,10 -> 20,00 10,01 -> 00,11 20,00 -> 10,10 6D scattering calculation Not rigidrotor 6D potential energy surface (Hinde 2008) 10-15 100 1000 Temperature (K) First for v=2 Fonseca dos Santos et al. 2013, J. Chem. Phys., 138, 104302
Summary Computational approaches can provide reliable data for astrophysical modeling Methods use state-of-the-art quantum techniques However, an on-going effort, with challenges... Many other processes need explicit investigation: Atom/ion fine-structure excitation by H, H2 Other grain-catalyzed reaction parameters (diffusion barriers, activation energies, etc.) Gas-phase neutral-neutral reactions, etc.
Pure Rotational Excitation: benchmarking experiment He + H2O(v=0,Jk-1,k+1) 10 Exp. (C. Yang et al. 2010) Cross section (10-16 cm 2 ) 1 0.1 SAPT-H SAPT-P VB SAPT-H (C. Yang et al. 2010) 1 10 2 12 2 21 3 03 3 12 3 21 4 14 Final state Yang et al. 2013, ApJ, 765, 77