[Toward] Simulating Cosmological Reionization with Enzo
|
|
- Lynn Lawson
- 5 years ago
- Views:
Transcription
1 [Toward] Simulating Cosmological Reionization with Enzo Texas Cosmology Network Meeting October 2009 Daniel R. Reynolds M.L. Norman, P. Paschos, G. So [UCSD] J.C. Hayes [LLNL] Sponsor: NSF Astronomy & Astrophysics October 29, 2009 (1)
2 Personal Introduction Assistant professor in Mathematics at SMU. Computational/applied math background; expertise in nonlinear and linear solvers for multi-scale problems in computational physics. This work focuses on collaborations with UCSD s Laboratory for Computational Astrophysics (M. Norman et al.). Goal is to develop high-accuracy, scalable solvers for coupled simulations of radiation, hydrodynamics and ionization. Efforts focused on developing solver modules for Enzo, an open-source code for cosmological hydro., dark-matter, gravity and chemical ionization (advertised by Kim Tran earlier). Target applications are mainly in cosmological reionization, though goal is to develop numerical methods that span to other regimes as well. October 29, 2009 (2)
3 Enzo-RT Goals I. Extend Enzo to include radiation transport & ionization, to enable: studies of self-regulated star formation, predictions on the epoch of cosmic reionization, and predictions on the observed properties of early galaxies. II. Accurately model stiff cosmological RT and chemical ionization processes. III. Integrate new solver in Enzo code so that coupling respects hydrodynamics approach (shocks), but tightly couples physics. IV. Enable very-large-scale simulations (O(10 6 ) processors). October 29, 2009 (3)
4 Coupled Matter-Radiation System We consider the coupled cosmological PDE system, 2 φ = 4πG a (ρ b + ρ dm ρ 0 ), t ρ b + 1 a v b ρ b = 1 a ρ b v b, t v b + 1 a (v b ) v b = ȧ a v b 1 aρ b p 1 a φ, t e + 1 a v b e = 2ȧ a e 1 aρ b (pv b ) 1 a v b φ + G Λ, t n i + 1 a (n iv b )= 3 ȧ a n i n i Γ ph i + α rec i,j n en j, t E ν + 1 a (E νv b )= (D E ν )+ν ȧ a νe ν +4πη ν ck ν E ν. Here, G(E,n i ) and Λ(E,n i ) are the heating and cooling rates. The E ν equation approximates the radiative flux as a function of the energy density gradient, F ν = D E ν, where D(E ν, E ν ) IR 3 3 is the flux limiter. [Bryan et al., Comp. Phys. Comm., 1995; R. et al., J. Comput. Phys., 2009] October 29, 2009 (4)
5 Enzo Operator-Split Numerics We solve this coupled system in an operator-split framework, solving one component of the system at a time within a time step: (i) Project dark matter particles onto mesh to generate ρ dm. (ii) Solve for the gravitational potential φ via FFT or MG methods. (iii) Advect dark matter particles via Particle-Mesh method. (iv) Explicitly evolve (ρ b, v b,e) and advect (E ν, n i ) with a high-order PPM method. (v) Implicitly evolve a stiffreaction-di ffusion PDE system that updates (E ν, e, n i ). October 29, 2009 (5)
6 Enzo Operator-Split Numerics We solve this coupled system in an operator-split framework, solving one component of the system at a time within a time step: (i) Project dark matter particles onto mesh to generate ρ dm. (ii) Solve for the gravitational potential φ via FFT or MG methods. (iii) Advect dark matter particles via Particle-Mesh method. (iv) Explicitly evolve (ρ b, v b,e) and advect (E ν, n i ) with a high-order PPM method. (v) Implicitly evolve a stiffreaction-di ffusion PDE system that updates (E ν, e, n i ). October 29, 2009 (5)
7 Gray FLD Radiation Transfer Equation Prior to investigating the multi-frequency case, we begin using a single, frequency-integrated ( grey ) radiation energy, E(x, t) = Z ν 0 E ν (x, t,ν ) dν = Ẽ(x, t) Z ν 0 χ E (ν) dν, where χ E : IR IR is an assumed radiation energy density spectrum. With this approximation, the integrated radiation energy equation becomes t E + 1 a (Ev b) = (D E) ȧ E +4πη cke, a where the new coupling terms are integrated over frequency: η = k = Z ν 0 Z ν 0 η ν dν «. Z «k ν E ν dν E ν dν. ν 0 October 29, 2009 (6)
8 Rad-Hydro Verification Results Radiation Energy Density (erg cm -3 ) Distance (cm) Smearing of a radiation front in a vacuum due to FLD approximation (convergence wrt x). T Distance From Shock (cm) Radiating shock front convergence wrt x. [Lowrie & Edwards, ShWav, 2008]. Thermal equilibration: energy conserv. error, Z E tot dx«. Z «E tot (0)dx, L H and H L equilibration, under various nonlinear and linear tolerances (ε, δ). [Turner & Stone, ApJ, 2001] ε δ L H H L 1e-7 1e-9 7.0e e-12 1e-10 1e-9 3.1e e-14 1e-4 1e-9 3.6e e-11 1e-7 1e e e-12 1e-7 1e-3 7.0e e-12 October 29, 2009 (7)
9 Isothermal Radiative Ionization Test Tests of isothermal ionization of a static, initially-neutral hydrogen region: L =6.6 kpc. Monochromatic emission: Ṅ γ = photons s 1, at hν = 13.6 ev. Fixed gas temperature, T = 10 4 K. Case-B recombination rate, α B = cm 3 s 1. Constant number density n = 10 3 cm 3. Analytical solution given by r I = r Sˆ1 e tα B n H, where " # 1/3 3Ṅγ r S = 4πα B n 2. H [Iliev et al., MNRAS, 2006] October 29, 2009 (8)
10 Isothermal Ionization Results $!") $( & 6;-:< &% & 6;-:< (' & 6;-:< $%) & 6;-:< =3=>?*7.=> 1234-, /!8,23*692:7*723 6 NW: I-front position history.!"+!"&!"*!"% -./, , ,4,:,'!!,;80 ',!')!#!!#)!(!, / +, 0!"(!"'!"%! 6!!"# $ $"# % %"# *+*,-. & &"# ' '"# NE: Coarse radiation field (N x = 16).!")!"$!"(!"#!"'!,!!"#!"$!"%!"& '!()!$!!$)!)! /, 0!, 0 1 /+,2 3,,4,)56)0!7,46*)8495*546 # ()$!!& ) $: & );-9< '"# &% & );-9< ' :' & );-9< $%= & );-9< &"# & %"# % $"# $!"#! )!!"# $ $"# % *+*,-. %"# & &"# ' SW: r I convergence vs x. SE: Fine radiation field (N x = 128).!"+!"&!"*!"%!")!"$!"(!"#!"' -./, , ,4,:,'!!,;80 ',!#!!(!!$!!)!!%!!*!!&!!+!!,!!"#!"$!"%!"& ' October 29, 2009 (9)
11 Isothermal Cosmological Ionization Test Repeat of previous test, but in a cosmologically expanding universe. Four tests: q 0 z i L i [kpc] ρ b,i [g cm 3 ] H 0 Ωm Ω Λ Ω b e e e e Analytical solution given by r I (t) =r S,i Z 1/3 a(t) λe τ(t) e τ(b) [1 2q 0 +2q 0 (1 + z i )/b] db! 1/2, 1 where τ(a) =λ [F (a) F (1)] ˆ6q 0(1 2 + z i ) 2 1,»» 1+z i 1+z 1/2 i F (a) = 2 4q 0 2q 0 1 2q 0 +2q 0, a a λ = α B n H,i /H 0 /(1 + z i ). [Shapiro & Giroux, ApJ, 1987] October 29, 2009 (10)
12 $!"+!"*!") Cosmological Ionization Results : NW: I-front radii vs scaled z. ##,-#!!% #! + * >99/9-7?-9 7 2:569 = 2:5-@=-9<A=B7C:D-E!-F-!"'D-4!-F-& #( % -G<=B %$ % -G<=B (& % -G<=B #$* % -G<=B !"(!"#!"'!"&!"%!"$ NE: I-front error vs scaled z. 29 :9;<!9569 = ) ( ' & % $! :!!"# $ $"# % %"# &!,-./0$12340$ :;6<3=7>80?& 3 6 3!"( 8!"' SW: I-front radii for z i =4. # -!!"#!"$!"%!"&!"'!"(!")!./012#34562# Weak Scaling Results for Cosmological Ionization Test 3000!"& !!"#!!"#!"$!"%!"&!"'!"(!")!*+,-.#/012.#/ SE: Weak CPU scaling [Kraken] N src N CPU (O(N log N)) time (s) procs October 29, 2009 (11)
13 Hydrodynamic Radiative Ionization Test Re-do first test in a dynamic medium (L = 15 kpc); combines radiation, hydrodynamics and hydrogen ionization: Emission source has a T = 10 5 K blackbody spectrum. Initial temperature set to T = 10 2 K. Front transitions from R-type to D-type as it reaches Strömgren radius. Eventually stalls at r f = r S 2Ti T e 2/3, where Ti and T e are the temperatures behind and ahead of the i-front. R-type and D-type propgation have different analytical solutions, r R I = r S i» 1/3 h1 e tα B(T i )n H, r D I = r S 1+ 7c 4/7 st, 4r S with true solution somewhere between the two. [Whalen & Norman, ApJS, 2006; Iliev et al., MNRAS, 2009] October 29, 2009 (12)
14 Hydrodynamic Radiative Ionization Results %"# % $"# /012+*3+1,+4514-!6*01(47085(501 $9 & 4:+8; &% & 4:+8; 9' & 4:+8; $%< & 4:+8; 4 NW: I-front position history Temperature profile, t = 175 Myr * - )*. log(t) 3 $!"#! 4!!"# $ $"# % %"# ()( *+, & &"# ' '"# NE: T profile at 200 Myr r/l box -*.!*. / -)*0 1**2*3453.!6*25(37284(425!"!# 3 $9 & 3:+8;!"!'# &% & 3:+8;!"!' 9' & 3:+8;!"!&#!"!&!"!%#!"!%!"!$#!"!$!"!!#! 3!!"# $ $"# % %"# ()( *+, & &"# ' '"# SW: r I convergence wrt x. SE: T convergence wrt t /+ #!!# #!!$ #!!% #!!& #!!' #!!* #!!) #!!( 7 '1!& $1!& #1!& '1!' $1!' #1!' 819:1+34;+171++/+7:+/<521=747>7#)'7?@+!!"#!"$!"%!"&!"' +,-./0 7 October 29, 2009 (13)
15 The Need for Multi-Frequency However, in comparison with other RT codes, we get peculiar results: Right: n profiles at 250 Myr (cyan is Enzo) [1]. Below: Enzo n profiles results at 200, 500 Myr. Lower right: Zeus-MP results (solid is MF, dashed is monochromatic). '!!# 012,3)453/67894:)-;7<3 #!!4=9) >!!4=9) 4 / '!!( n (cm -3 ) 10-3 '!!$ 4!!"#!"$!"%!"& ' '"# )*+, r/l box October 29, 2009 (14)
16 Summary of Current Results We ve achieved a second-order accurate (space and time) coupled solver for cosmological radiation, hydrodynamics, self gravity and chemical ionization. Captures shock fronts due to trusted PPM hydrodynamics approach. Accurately solves couplings between radiation, ionization and gas energy, due to implicit formulation and coupled solvers. Highly scalable, with current tests using > 4000 processors, thanks to reliance on optimal multigrid methods. However, the grey radiation approximation has its shortcomings: Single radiation field allows full absorption by hydrogen, even though higher-frequency radiation should pass through. Leads to increased absorption in optically-thick regions, with no pre-heating ahead of I-fronts. October 29, 2009 (15)
17 Continuing Work Enhance stability of chemistry computations at large t Fastest dynamics occur in ionization, ostensibly allowing complete ionization of a volume in a single step [O(1) O(10 5 )]. Time inaccuracy can over-shoot to give negative densities. Extend radiation approximation to multi-frequency case, E E ν New physical couplings and interpolation of ν-space, based offof [Ricotti, Gnedin & Shull, ApJ, 2002] Large-scale solvers for coupled reaction-diffusion systems Extend implicit solver software to AMR grids Regular-grid geometric multigrid AMG or FAC methods October 29, 2009 (16)
Scalable Radiation-Induced Chemical Ionization Simulating Cosmological Reionization
Scalable Radiation-Induced Chemical Ionization Simulating Cosmological Reionization Daniel R. Reynolds reynolds@smu.edu Department of Mathematics Southern Methodist University KAUST Seminar, March 2011
More informationFlux Limited Diffusion in Enzo
Flux Limited Diffusion in Enzo Daniel R. Reynolds reynolds@smu.edu Department of Mathematics Southern Methodist University Enzo User Workshop, 29 June 2010 Reading The slides (Enzo2010 talk1 reynolds.pdf)
More informationBlack Hole Formation in the Early Universe
1 Black Hole Formation in the Early Universe Luo Yang (Osaka University) XC30 2 1. Introduction Observations of high-redshift quasars reveal that some supermassive black holes (SMBH) with masses exceeding
More informationCosmological Radiation Hydrodynamics with Enzo
Cosmological Radiation Hydrodynamics with Enzo Michael L. Norman,, Daniel R. Reynolds and Geoffrey C. So Physics Department, U.C. San Diego, La Jolla, CA 92093 Ctr. for Astrophysics and Space Sciences,
More informationRadiative Transfer in a Clumpy Universe: the UVB. Piero Madau UC Santa Cruz
Radiative Transfer in a Clumpy Universe: the UVB Piero Madau UC Santa Cruz The cosmic UVB originates from the integrated emission of starforming galaxies and QSOs. It determines the thermal and ionization
More informationNumerical Models of the high-z Universe
Texte Numerical Models of the high-z Universe Dominique AUBERT Observatoire Astronomique, Université de Strasbourg EOR Robertson et al. 2010 Epoch of Reionization ~200 Myrs - 1Gyr z~30-6! Challenge : Multiple
More informationX i t react. ~min i max i. R ij smallest. X j. Physical processes by characteristic timescale. largest. t diff ~ L2 D. t sound. ~ L a. t flow.
Physical processes by characteristic timescale Diffusive timescale t diff ~ L2 D largest Sound crossing timescale t sound ~ L a Flow timescale t flow ~ L u Free fall timescale Cooling timescale Reaction
More information(c) (a) 3kT/2. Cascade
1 AY30-HIITemp IV. Temperature of HII Regions A. Motivations B. History In star-forming galaxies, most of the heating + cooling occurs within HII regions Heating occurs via the UV photons from O and B
More information6. Cosmology. (same at all points) probably true on a sufficiently large scale. The present. ~ c. ~ h Mpc (6.1)
6. 6. Cosmology 6. Cosmological Principle Assume Universe is isotropic (same in all directions) and homogeneous (same at all points) probably true on a sufficiently large scale. The present Universe has
More informationFORMATION OF PRIMORDIAL STARS
Talk@INT, UW, July 5, 2006 FORMATION OF PRIMORDIAL STARS Naoki Yoshida Department of Physics Nagoya University Outline Thermal evolution of a primordial gas - Physics at high densities (cooling, chem.
More informationAdvanced Cosmological Simulations
Advanced Cosmological Simulations John Wise (Georgia Tech) Enzo Workshop 19 Oct 2013 1 Outline We will consider additional physics in Thursday s AMR (no nested grids) cosmology simulation. Refresher on
More informationII. HII Regions (Ionization State)
1 AY230-HIIReg II. HII Regions (Ionization State) A. Motivations Theoretical: HII regions are intamitely linked with past, current and future starforming regions in galaxies. To build theories of star-formation
More informationPhotoionized Gas Ionization Equilibrium
Photoionized Gas Ionization Equilibrium Ionization Recombination H nebulae - case A and B Strömgren spheres H + He nebulae Heavy elements, dielectronic recombination Ionization structure 1 Ionization Equilibrium
More informationLecture 27 The Intergalactic Medium
Lecture 27 The Intergalactic Medium 1. Cosmological Scenario 2. The Ly Forest 3. Ionization of the Forest 4. The Gunn-Peterson Effect 5. Comment on HeII Reionization References J Miralda-Escude, Science
More informationRadiative & Magnetohydrodynamic Shocks
Chapter 4 Radiative & Magnetohydrodynamic Shocks I have been dealing, so far, with non-radiative shocks. Since, as we have seen, a shock raises the density and temperature of the gas, it is quite likely,
More informationProtostars 1. Early growth and collapse. First core and main accretion phase
Protostars 1. First core and main accretion phase Stahler & Palla: Chapter 11.1 & 8.4.1 & Appendices F & G Early growth and collapse In a magnetized cloud undergoing contraction, the density gradually
More informationPreliminary Examination: Astronomy
Preliminary Examination: Astronomy Department of Physics and Astronomy University of New Mexico Spring 2017 Instructions: Answer 8 of the 10 questions (10 points each) Total time for the test is three
More informationComputational Astrophysics 7 Hydrodynamics with source terms
Computational Astrophysics 7 Hydrodynamics with source terms Oscar Agertz Outline - Optically thin radiative hydrodynamics - Relaxation towards the diffusion limit - Hydrodynamics with gravity source term
More informationRadiation-hydrodynamics from Mpc to sub-pc scales with RAMSES-RT
Radiation-hydrodynamics from Mpc to sub-pc scales with RAMSES-RT Joki Rosdahl Centre de Recherche Astrophysique de Lyon With Aubert, Blaizot, Bieri, Biernacki, Commercon, Costa, Dubois, Geen, Katz, Kimm,
More informationMichael Shull (University of Colorado)
Early Galaxies, Stars, Metals, and the Epoch of Reionization Michael Shull (University of Colorado) Far-IR Workshop (Pasadena, CA) May 29, 2008 Submillimeter Galaxies: only the brightest? How long? [dust
More informationOutline of lectures. Adding more physics. Doing Better: source terms in CTU unsplit integrator
Grid-based methods for hydrodynamics, MHD, and radiation hydrodynamics. Outline of lectures (Four lectures) Lecture 1. Introduction to physics and numerics Lecture 2. Operator split (ZEUS-like) methods
More informationH II Regions of the First Stars II: A Primer on I-front Instabilities. Dan Whalen UC San Diego
H II Regions of the First Stars II: A Primer on I-front Instabilities Dan Whalen UC San Diego Our Collaboration Daniel Whalen, T-6, LANL; UC San Diego; UIUC Brian O Shea, T-6, LANL Alex Heger, T-6, LANL
More informationCore Collapse Supernovae An Emerging Picture Stephen W. Bruenn
Core Collapse Supernovae An Emerging Picture Stephen W. Bruenn 19th Rencontres de Blois Matter and Energy in the Universe: from nucleosynthesis to cosmology Collaborators Anthony Mezzacappa John M. Blondin
More informationPhotoionization Modelling of H II Region for Oxygen Ions
Journal of Materials Science and Chemical Engineering, 2015, 3, 7-16 Published Online April 2015 in SciRes. http://www.scirp.org/journal/msce http://dx.doi.org/10.4236/msce.2015.34002 Photoionization Modelling
More informationThe improvement of START
The improvement of START Kenji Hasegawa (U. Tsukuba, CCS Kobe branch) Takashi Okamoto (U. Tsukuba, CCS Kobe branch) Cosmological Radiative Transfer Comparison Project Workshop IV @ Austin, Texas, Dec 11-14,
More informationGalaxies 626. Lecture 5
Galaxies 626 Lecture 5 Galaxies 626 The epoch of reionization After Reionization After reionization, star formation was never the same: the first massive stars produce dust, which catalyzes H2 formation
More informationThe Stellar Opacity. F ν = D U = 1 3 vl n = 1 3. and that, when integrated over all energies,
The Stellar Opacity The mean absorption coefficient, κ, is not a constant; it is dependent on frequency, and is therefore frequently written as κ ν. Inside a star, several different sources of opacity
More informationOutline. Walls, Filaments, Voids. Cosmic epochs. Jeans length I. Jeans length II. Cosmology AS7009, 2008 Lecture 10. λ =
Cosmology AS7009, 2008 Lecture 10 Outline Structure formation Jeans length, Jeans mass Structure formation with and without dark matter Cold versus hot dark matter Dissipation The matter power spectrum
More informationGasdynamical and radiative processes, gaseous halos
Gasdynamical and radiative processes, gaseous halos Houjun Mo March 19, 2004 Since luminous objects, such as galaxies, are believed to form through the collapse of baryonic gas, it is important to understand
More informationOpen-Source Astrophysics with the Enzo Community Code. Brian W. O Shea Michigan State University
Open-Source Astrophysics with the Enzo Community Code Brian W. O Shea Michigan State University What is Enzo? N-body + hydro (+MHD + RHD +...) blockstructured AMR code originally written at NCSA (by Greg
More informationNumerical Simulations. Duncan Christie
Numerical Simulations Duncan Christie Motivation There isn t enough time to derive the necessary methods to do numerical simulations, but there is enough time to survey what methods and codes are available
More informationThe first stars and primordial IMBHs
The first stars and primordial IMBHs Ab initio predictions of black hole merger rates by the time LISA flies? Tom Abel Penn State Initial Conditions Time Evolution z=100 z=24 z=20.4 10 comoving kpc Cosmological
More informationPractical in Numerical Astronomy, SS 2012 LECTURE 9
Practical in Numerical Astronomy, SS 01 Elliptic partial differential equations. Poisson solvers. LECTURE 9 1. Gravity force and the equations of hydrodynamics. Poisson equation versus Poisson integral.
More informationOpacity and Optical Depth
Opacity and Optical Depth Absorption dominated intensity change can be written as di λ = κ λ ρ I λ ds with κ λ the absorption coefficient, or opacity The initial intensity I λ 0 of a light beam will be
More informationThermal Equilibrium in Nebulae 1. For an ionized nebula under steady conditions, heating and cooling processes that in
Thermal Equilibrium in Nebulae 1 For an ionized nebula under steady conditions, heating and cooling processes that in isolation would change the thermal energy content of the gas are in balance, such that
More informationENZO+MORAY: radiation hydrodynamics adaptive mesh refinement simulations with adaptive ray tracing
Mon. Not. R. Astron. Soc. 414, 3458 3491 (2011) doi:10.1111/j.1365-2966.2011.18646.x ENZO+MORAY: radiation hydrodynamics adaptive mesh refinement simulations with adaptive ray tracing John H. Wise 1 and
More information80 2 Observational Cosmology L and the mean energy
80 2 Observational Cosmology fluctuations, short-wavelength modes have amplitudes that are suppressed because these modes oscillated as acoustic waves during the radiation epoch whereas the amplitude of
More informationReionization constraints post Planck-15
Reionization constraints post Planck-15 Tirthankar Roy Choudhury National Centre for Radio Astrophysics Tata Institute of Fundamental Research Pune CMB Spectral Distortions from Cosmic Baryon Evolution
More informationTheoretical Cosmology and Galaxy Formation at UMD
Theoretical Cosmology and Galaxy Formation at UMD Massimo Ricotti (Associate Professor, Dept. of Astronomy) Current group members: Owen Parry (Postdoc) Sam Leithner (CTC postdoc) Emil Polisensky (PhD student)
More informationBrief Introduction to Cosmology
Brief Introduction to Cosmology Matias Zaldarriaga Harvard University August 2006 Basic Questions in Cosmology: How does the Universe evolve? What is the universe made off? How is matter distributed? How
More informationP M 2 R 4. (3) To determine the luminosity, we now turn to the radiative diffusion equation,
Astronomy 715 Final Exam Solutions Question 1 (i). The equation of hydrostatic equilibrium is dp dr GM r r 2 ρ. (1) This corresponds to the scaling P M R ρ, (2) R2 where P and rho represent the central
More informationBaryon Acoustic Oscillations in the Lyman Alpha Forest
Baryon Acoustic Oscillations in the Lyman Alpha Forest Michael L. Norman, Pascal Paschos & Robert Harkness Laboratory for Computational Astrophysics, SERF 0424, University of California San Diego, 9500
More informationASTRONOMY QUALIFYING EXAM August Possibly Useful Quantities
L = 3.9 x 10 33 erg s 1 M = 2 x 10 33 g M bol = 4.74 R = 7 x 10 10 cm 1 A.U. = 1.5 x 10 13 cm 1 pc = 3.26 l.y. = 3.1 x 10 18 cm a = 7.56 x 10 15 erg cm 3 K 4 c= 3.0 x 10 10 cm s 1 σ = ac/4 = 5.7 x 10 5
More informationThe RAMSES code and related techniques 4. Source terms
The RAMSES code and related techniques 4. Source terms Outline - Optically thin radiative hydrodynamics - Relaxation towards the diffusion limit - Hydrodynamics with gravity source term - Relaxation towards
More informationHII regions. Massive (hot) stars produce large numbers of ionizing photons (energy above 13.6 ev) which ionize hydrogen in the vicinity.
HII regions Massive (hot) stars produce large numbers of ionizing photons (energy above 13.6 ev) which ionize hydrogen in the vicinity. Detailed nebular structure depends on density distribution of surrounding
More informationAtomic Physics 3 ASTR 2110 Sarazin
Atomic Physics 3 ASTR 2110 Sarazin Homework #5 Due Wednesday, October 4 due to fall break Test #1 Monday, October 9, 11-11:50 am Ruffner G006 (classroom) You may not consult the text, your notes, or any
More informationAST4320: LECTURE 10 M. DIJKSTRA
AST4320: LECTURE 10 M. DIJKSTRA 1. The Mass Power Spectrum P (k) 1.1. Introduction: the Power Spectrum & Transfer Function. The power spectrum P (k) emerged in several of our previous lectures: It fully
More informationCosmological simulations of X-ray heating during the Universe s Dark Ages
Cosmological simulations of X-ray heating during the Universe s Dark Ages Jordan Mirocha 1,5, Jack Burns 1,5, Eric Hallman 2,5, Steven Furlanetto 3,6, John Wise 4 1 University of Colorado at Boulder 2
More informationReally, what universe do we live in? White dwarfs Supernova type Ia Accelerating universe Cosmic shear Lyman α forest
Really, what universe do we live in? White dwarfs Supernova type Ia Accelerating universe Cosmic shear Lyman α forest White dwarf Core of solar mass star No energy from fusion or gravitational contraction
More informationa few more introductory subjects : equilib. vs non-equil. ISM sources and sinks : matter replenishment, and exhaustion Galactic Energetics
Today : a few more introductory subjects : equilib. vs non-equil. ISM sources and sinks : matter replenishment, and exhaustion Galactic Energetics photo-ionization of HII assoc. w/ OB stars ionization
More informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
More informationSome HI is in reasonably well defined clouds. Motions inside the cloud, and motion of the cloud will broaden and shift the observed lines!
Some HI is in reasonably well defined clouds. Motions inside the cloud, and motion of the cloud will broaden and shift the observed lines Idealized 21cm spectra Example observed 21cm spectra HI densities
More informationThe Deflagration Phase of Type Ia SNe
The Center for Astrophysical Thermonuclear Flashes The Deflagration Phase of Type Ia SNe Alan Calder ASC FLASH Center Type Ia Supernova Team Type Ia Supernovae and Cosmology August 5, 2004 An Advanced
More informationDiffuse Interstellar Medium
Diffuse Interstellar Medium Basics, velocity widths H I 21-cm radiation (emission) Interstellar absorption lines Radiative transfer Resolved Lines, column densities Unresolved lines, curve of growth Abundances,
More informationDetailed Study of a Turbulent multiphase multicomponent ISM
Detailed Study of a Turbulent multiphase multicomponent ISM Dieter Breitschwerdt Collaborators Miguel de Avillez (Evora, Portugal) Verena Baumgartner (Vienna, Austria) Jan Bolte (TU Berlin, Germany) Jenny
More informationThe Formation of Population III Protostars. Matthew Turk - Stanford/KIPAC with Tom Abel (Stanford/KIPAC)
The Formation of Population III Protostars Matthew Turk - Stanford/KIPAC with Tom Abel (Stanford/KIPAC) (preamble) Population III Stars Population III Population III Population III.1 Population III.2 Population
More informationCosmic Variance of Small-Scale Structure Formation: Large-Scale Density and CDM-Baryon Drift Velocity Environment
Cosmic Variance of Small-Scale Structure Formation: Large-Scale Density and CDM-Baryon Drift Velocity Environment Based on collaboration with: Paul Shapiro (Texas), Ilian Iliev (Sussex), Garrelt Mellema
More informationMaking Sense of the Universe with Supercomputers
Making Sense of the Universe with Supercomputers Tom Abel Kavli Institute for Particle Astrophysics and Cosmology, Stanford, SLAC Adaptive Mesh Refinement in Cosmology: First Stars Adaptive Ray-Tracing
More informationEnergy Balance in Clusters of Galaxies. Patrick M. Motl & Jack O. Burns Center for Astrophysics and Space Astronomy University of Colorado at Boulder
Energy Balance in Clusters of Galaxies Patrick M. Motl & Jack O. Burns Center for Astrophysics and Space Astronomy University of Colorado at Boulder X-ray and Radio Connections, February 6th, 2004 With
More informationNew Results from 3-D supernova models with spectral neutrino diffusion
New Results from 3-D supernova models with spectral neutrino diffusion Stuart C. Whitehouse and Matthias Liebendörfer Department of Physics, University of Basel, Switzerland Overview Introduction The Isotropic
More informationTHE ROLE OF RADIATION PRESSURE IN HIGH-Z DWARF GALAXIES
THE ROLE OF RADIATION PRESSURE IN HIGH-Z DWARF GALAXIES John Wise (Georgia Tech) Tom Abel (Stanford), Michael Norman (UC San Diego), Britton Smith (Michigan State), Matthew Turk (Columbia) 14 Dec 2012
More informationCosmological shock waves
Highlights of Spanish Astrophysics VII, Proceedings of the X Scientific Meeting of the Spanish Astronomical Society held on July 9-13, 2012, in Valencia, Spain. J. C. Guirado, L. M. Lara, V. Quilis, and
More informationElliptic Problems / Multigrid. PHY 604: Computational Methods for Physics and Astrophysics II
Elliptic Problems / Multigrid Summary of Hyperbolic PDEs We looked at a simple linear and a nonlinear scalar hyperbolic PDE There is a speed associated with the change of the solution Explicit methods
More informationEmitted Spectrum Summary of emission processes Emissivities for emission lines: - Collisionally excited lines - Recombination cascades Emissivities
Emitted Spectrum Summary of emission processes Emissivities for emission lines: - Collisionally excited lines - Recombination cascades Emissivities for continuum processes - recombination - brehmsstrahlung
More informationAstrochemistry. Lecture 10, Primordial chemistry. Jorma Harju. Department of Physics. Friday, April 5, 2013, 12:15-13:45, Lecture room D117
Astrochemistry Lecture 10, Primordial chemistry Jorma Harju Department of Physics Friday, April 5, 2013, 12:15-13:45, Lecture room D117 The first atoms (1) SBBN (Standard Big Bang Nucleosynthesis): elements
More informationAST242 LECTURE NOTES PART 7
AST242 LECTURE NOTES PART 7 Contents 1. HII regions and Ionization Fronts 1 1.1. The Strömgren Sphere 2 1.2. Early Evolution 3 1.3. Achieving Pressure equilibrium 3 1.4. Jump conditions on an ionization
More informationEmission lines from galaxies in the Epoch of Reioniza5on
Emission lines from galaxies in the Epoch of Reioniza5on Simona Gallerani Emanuele Sobacchi, Andrea Ferrara, Livia Vallini, Andrea PalloBni, CarloDa Gruppioni, Andrei Mesinger Monte Mario, Roma, SPICA
More informationM.Phys., M.Math.Phys., M.Sc. MTP Radiative Processes in Astrophysics and High-Energy Astrophysics
M.Phys., M.Math.Phys., M.Sc. MTP Radiative Processes in Astrophysics and High-Energy Astrophysics Professor Garret Cotter garret.cotter@physics.ox.ac.uk Office 756 in the DWB & Exeter College Radiative
More informationEffects of Massive Stars
Effects of Massive Stars Classical HII Regions Ultracompact HII Regions Stahler Palla: Sections 15.1, 15. HII Regions The salient characteristic of any massive star is its extreme energy output, much of
More informationTwo Phase Formation of Massive Galaxies
Two Phase Formation of Massive Galaxies Focus: High Resolution Cosmological Zoom Simulation of Massive Galaxies ApJ.L.,658,710 (2007) ApJ.,697, 38 (2009) ApJ.L.,699,L178 (2009) ApJ.,725,2312 (2010) ApJ.,744,63(2012)
More informationLec 22 Physical Properties of Molecular Clouds
Lec 22 Physical Properties of Molecular Clouds 1. Giant Molecular Clouds 2. Orion s Clouds 3. Correlations of Observed Properties 4. The X-Factor References Origins of Stars & Planetary Systems eds. Lada
More informationThe Interstellar Medium
http://www.strw.leidenuniv.nl/~pvdwerf/teaching/ The Interstellar Medium Lecturer: Dr. Paul van der Werf Fall 2014 Oortgebouw 565, ext 5883 pvdwerf@strw.leidenuniv.nl Assistant: Kirstin Doney Huygenslaboratorium
More informationJournal of Computational Physics
Journal of Computational Physics 228 (29) 6833 6854 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp Self-consistent solution of
More informationThe Probes and Sources of Cosmic Reionization Francesco Haardt University of Como INFN, Milano-Bicocca
1 The Probes and Sources of Cosmic Reionization Francesco Haardt University of Insubria@Lake Como INFN, Milano-Bicocca 2 TALK OUTLINE 1. Dark Ages and Reionization 2. Observations: QSO Absorption Lines
More information1. Why photons? 2. Photons in a vacuum
Photons and Other Messengers 1. Why photons? Ask class: most of our information about the universe comes from photons. What are the reasons for this? Let s compare them with other possible messengers,
More informationNUMERICAL METHODS IN ASTROPHYSICS An Introduction
-1 Series in Astronomy and Astrophysics NUMERICAL METHODS IN ASTROPHYSICS An Introduction Peter Bodenheimer University of California Santa Cruz, USA Gregory P. Laughlin University of California Santa Cruz,
More informationPoS(Cosmology2009)022
and 21cm Observations Max Planck Institute for Astrophysics E-mail: ciardi@mpa-garching.mpg.de With the advent in the near future of radio telescopes as LOFAR, a new window on the highredshift universe
More informationObservational Cosmology
(C. Porciani / K. Basu) Lecture 7 Cosmology with galaxy clusters (Mass function, clusters surveys) Course website: http://www.astro.uni-bonn.de/~kbasu/astro845.html Outline of the two lecture Galaxy clusters
More informationResearch Statement Daniel R. Reynolds
Research Statement Daniel R. Reynolds My research is in applied mathematics, specifically the sub-disciplines of scientific computation and numerical analysis. In the ever-evolving field of scientific
More informationPAPER 73 PHYSICAL COSMOLOGY
MATHEMATICAL TRIPOS Part III Wednesday 4 June 2008 1.30 to 4.30 PAPER 73 PHYSICAL COSMOLOGY Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal weight. STATIONERY
More informationAstronomy 102: Stars and Galaxies Sample Review Test for Examination 3
October 28, 2003 Name: Astronomy 102: Stars and Galaxies Sample Review Test for Examination 3 Do not open the test until instructed to begin. Instructions: Write your answers in the space provided. No
More informationarxiv: v2 [astro-ph.im] 15 Mar 2011
Mon. Not. R. Astron. Soc., 37 () Printed 3 September 8 (MN LATEX style file v.) Enzo+Moray: Radiation Hydrodynamics Adaptive Mesh Refinement Simulations with Adaptive Ray Tracing arxiv:.86v [astro-ph.im]
More informationThe Tempest Simulations Milky Way-type Galaxies, their Environments and Progenitors
The Tempest Simulations Milky Way-type Galaxies, their Environments and Progenitors John Wise Brian O Shea (PI, MSU) David Collins (FSU), Cameron Hummels (Caltech), Devin Silvia (MSU), Britton Smith (SDSC)
More informationGrowing and merging massive black holes
Growing and merging massive black holes Marta Volonteri Institut d Astrophysique de Paris S. Cielo (IAP) R. Bieri (MPA) Y. Dubois (IAP) M. Habouzit (Flatiron Institute) T. Hartwig (IAP) H. Pfister (IAP)
More information2. Basic Assumptions for Stellar Atmospheres
2. Basic Assumptions for Stellar Atmospheres 1. geometry, stationarity 2. conservation of momentum, mass 3. conservation of energy 4. Local Thermodynamic Equilibrium 1 1. Geometry Stars as gaseous spheres!
More informationCosmic ray feedback in hydrodynamical simulations. simulations of galaxy and structure formation
Cosmic ray feedback in hydrodynamical simulations of galaxy and structure formation Canadian Institute for Theoretical Astrophysics, Toronto April, 13 26 / Workshop Dark halos, UBC Vancouver Outline 1
More informationOverview spherical accretion
Spherical accretion - AGN generates energy by accretion, i.e., capture of ambient matter in gravitational potential of black hole -Potential energy can be released as radiation, and (some of) this can
More informationInterstellar Medium Physics
Physics of gas in galaxies. Two main parts: atomic processes & hydrodynamic processes. Atomic processes deal mainly with radiation Hydrodynamics is large scale dynamics of gas. Start small Radiative transfer
More informationDisk Galaxy Rotation Curves and Dark Matter Halos
Disk Galaxy Rotation Curves and Dark Matter Halos Thorben H. Mense 1 Galactic Astronomy SS 2018 1 Universität Bielefeld Outline Early Observations of Galaxy Rotation Curves Rotation of the Andromeda Nebula
More informationAstrophysical Quantities
Astr 8300 Resources Web page: http://www.astro.gsu.edu/~crenshaw/astr8300.html Electronic papers: http://adsabs.harvard.edu/abstract_service.html (ApJ, AJ, MNRAS, A&A, PASP, ARAA, etc.) General astronomy-type
More informationCosmic Dawn (CoDa): Radiation-hydrodynamics of galaxy formation during the EoR
Cosmic Dawn (CoDa): Radiation-hydrodynamics of galaxy formation during the EoR Ocvirk+2015: arxiv:1511.00011 P. Ocvirk Observatoire astronomique de Strasbourg Universite de Strasbourg D. Aubert, N. Gillet,
More informationSIMPLE RADIATIVE TRANSFER
ASTR 511/O Connell Lec 4 1 SIMPLE RADIATIVE TRANSFER The theory of radiative transfer provides the means for determining the emergent EM spectrum of a cosmic source and also for describing the effects
More information3D Simulations of Core-collapse Supernovae. Tomoya Takiwaki(NAOJ) Kei Kotake(Fukuoka U) Yudai Suwa(YITP) Tomohide Wada(vis) And many collaborators
2013/12/3 MMCOCOS@Fukuoka University 3D Simulations of Core-collapse Supernovae Tomoya Takiwaki(NAOJ) Kei Kotake(Fukuoka U) Yudai Suwa(YITP) Tomohide Wada(vis) And many collaborators Plan 1. Brief summary
More informationDiffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics
Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods
More informationWhat can we learn about reionization from the ksz
What can we learn about reionization from the ksz Andrei Mesinger Scuola Normale Superiore, Pisa IGM effect on CMB primary temperature anisotropies ionized IGM damps CMB temperature anisotropies through
More informationCOSMIC-RAY DRIVEN MAGNETIC FIELD DYNAMO IN GALAXIES
COSMIC-RAY DRIVEN MAGNETIC FIELD DYNAMO IN GALAXIES Michał Hanasz, Centre for Astronomy Nicolaus Copernicus University, Toruń MAGNETIC FIELDS IN SPIRAL GALAXIES - RADIO OBSERVATIONS M51 NGC891 A. Fletcher
More informationVery deep spectroscopy of planetary nebula NGC7009 The rich optical recombination spectrum and new effective recombination coefficients
Very deep spectroscopy of planetary nebula NGC7009 The rich optical recombination spectrum and new effective recombination coefficients Xuan Fang Department of Astronomy (DoA), Peking University Collaborators:
More informationThe Fragmentation of expanding shells. Kazunari Iwasaki (Nagoya Univ.) Collaborate with S. Inutsuka (Nagoya Univ.) T. Tsuribe (Osaka Univ.
.... The Fragmentation of expanding shells Kazunari Iwasaki (Nagoya Univ.) Collaborate with S. Inutsuka (Nagoya Univ.) T. Tsuribe (Osaka Univ.) 1 / 23 Outline Introduction Three Dimensional SPH Simulations
More informationAstrophysics with the Computer: Propagation of Ionization Fronts in Interstellar Gas
Astrophysics with the Computer: Propagation of Ionization Fronts in Interstellar Gas Joachim Köppen Heidelberg 1992 1 Astrophysics A hot star is born in an interstellar gas cloud. At first all the gas
More informationNotes on Photoionized Regions Wednesday, January 12, 2011
Notes on Photoionized Regions Wednesday, January 12, 2011 CONTENTS: 1. Introduction 2. Hydrogen Nebulae A. Ionization equations B. Recombination coefficients and cross sections C. Structure of the hydrogen
More information