Cosmic Structure Formation on Supercomputers (and laptops)
|
|
- Donald Campbell
- 5 years ago
- Views:
Transcription
1 Cosmic Structure Formation on Supercomputers (and laptops) Lecture 4: Smoothed particle hydrodynamics and baryonic sub-grid models Benjamin Moster! Ewald Puchwein 1
2 Outline of the lecture course Lecture 1: Motivation & Initial conditions! Lecture 2: Gravity algorithms & parallelization! Lecture 3: Hydro schemes! Lecture 4: SPH, Radiative cooling & heating, Subresolution physics! Lecture 5: Halo and subhalo finders & Semi-analytic models! Lecture 6: Getting started with Gadget! Lecture 7: Example cosmological box! Lecture 8: Example galaxy collision 2 Benjamin Moster Numerical Structure Formation 4 IoA,
3 Outline of this lecture Smoothed particle hydrodynamics (SPH)! Idea & smoothing kernel! Equations of motion! Artificial viscosity! Advantages and problems (fluid-mixing, SPH blobs)! Modern SPH techniques (Pressure-entropy, artificial diffusion, etc.)! Baryonic sub-grid physics! Radiative Cooling and heating! Star formation! Stellar feedback (Multi-phase model, blastwave feedback)! AGN feedback 3 Benjamin Moster Numerical Structure Formation 4 IoA,
4 Smoothed particle hydrodynamics Idea: treat hydrodynamics in a completely mesh-free fashion Use: set of sampling particles to represent the fluid! Gas cloud of mass M g represented by N particles of mass m = Mg / N! Each particle i interacts with each other particle j through the symmetric SPH Kernel W(rij,h)! Depends on inter-particles distance and on smoothing-length h (defines Kernel width)! Most commonly the cubic spline Kernel is used: 4 W (r, h) = 8 h 3 r = ~x i ~x j 8 >< 1 6(r/h) 2 + 6(r/h) 3 0 apple (r/h) apple 1/2 2(1 r/h) >: 3 1/2 < (r/h) apple < (r/h) Benjamin Moster Numerical Structure Formation 4 IoA,
5 Smoothed particle hydrodynamics The density in SPH is given by smoothing the particle masses: (~x) = X j For any field F, we can define a smoothed interpolated version: F s (~x) = Z And the derivative is given by: ~rf s (~x) = X j mw(~x ~x j,h) d~x F(~x 0 ) W (~x ~x 0,h)= X j Usually smoothing length h is not constant:! Pressure is related to the density via: P i = A i i =( m j j F j ~ rw (~x ~xj,h) 1) i u i m j j F j W (~x ~x j,h) 4 3 h3 = N ngbm i i 5 Benjamin Moster Numerical Structure Formation 4 IoA,
6 Equations of motion Equations of motion derived from Lagrangian:! From we get where ~x i = ~ r j i d~v i m i dt = L = X i X j 1 2 m i~v i 2 m i u i m j P j 2 i and h is constant in the 1st term 6 Benjamin Moster Numerical Structure Formation 4 IoA,
7 Equations of motion Equations of motion derived from Lagrangian:! From we get where ~x i = ~ r j i Differentiation of j h 3 j = const : d~v i m i dt = L = X i X j 1 2 m i~v i 2 m i u i m j P j 2 i and h is constant in the 1st j = 1+ h 1 j ~ri i 3 j 7 Benjamin Moster Numerical Structure Formation 4 IoA,
8 Equations of motion Equations of motion derived from Lagrangian:! From we get where Differentiation of j h 3 j = const :! Using ~x i = ~ r j i ~r i j = m i ~ ri W ij (h j )+ ij X X j d~v we get: i dt = m j with f i = 1+ h i 3 i 1 f i P i 2 i d~v i m i dt = L = X i 1 2 m i~v i 2 m i u i, mass & energy automatically constant 8 Benjamin Moster Numerical Structure Formation 4 IoA, X j m j P j 2 i and h is constant in the 1st j = 1+ h 1 j ~ri i 3 j k m k ~ ri W ki (h i ) ~r i W ij (h i )+f j P j 2 j ~r j W ij (h j )!
9 Artificial viscosity Up to now: entropy is conserved along the flow (for each particle)! However, in shocks the entropy must increase! Artificial viscosity! Friction force:! d~v i dt vics = X j m j ij ~ ri Wij Since antisymmetric in i and j: (angular) momentum preserved To conserve total E: du i dt vics = 1 2 X m j ij ~v ij ri ~ Wij j 9 Benjamin Moster Numerical Structure Formation 4 IoA,
10 Artificial viscosity Up to now: entropy is conserved along the flow (for each particle)! However, in shocks the entropy must increase! Artificial viscosity! Friction force:! d~v i dt vics = X j m j ij ~ ri Wij Since antisymmetric in i and j: (angular) momentum preserved To conserve total E:! Usual choice: with and To prevent viscosity in shear flows (and conserve angular momentum): Balsara switch: du i dt vics = 1 2 with 10 Benjamin Moster Numerical Structure Formation 4 IoA, X j m j ij ~v ij ~ ri Wij ij =( c ij µ ij + µ 2 ij)/ ij µ ij = h ij~v ij ~x ij x ij 2 + h 2 ij ij =0 for! ~v ij ~x ij < 0 ij! f ij ij f AV i = ~ r ~v i ~ r ~v i + ~ r ~v i
11 Advantages and problems Advantages of SPH:! Numerically very robust! excellent conservation properties! Lagrangian nature trace flow & Galilean-invariant! Inherently adaptive (no ad-hoc refinement needed)! Couples properly to N-body gravity methods!! Problems?! SPH doesn t mix fluids (well)!,kaufmann-blobs form! Impacts accretion rates, SFRs disc sizes, gas fractions, etc Benjamin Moster Numerical Structure Formation 4 IoA,
12 Advantages and problems Advantages of SPH:! Numerically very robust! excellent conservation properties! Lagrangian nature trace flow & Galilean-invariant! Inherently adaptive (no ad-hoc refinement needed)! Couples properly to N-body gravity methods!! Problems?! SPH doesn t mix fluids (well)!,kaufmann-blobs form! Impacts accretion rates, SFRs disc sizes, gas fractions, etc Benjamin Moster Numerical Structure Formation 4 IoA,
13 Advantages and problems Advantages of SPH:! Numerically very robust! excellent conservation properties! Lagrangian nature trace flow & Galilean-invariant! Inherently adaptive (no ad-hoc refinement needed)! Couples properly to N-body gravity methods!! Problems?! SPH doesn t mix fluids (well)!,kaufmann-blobs form! Impacts accretion rates, SFRs disc sizes, gas fractions, etc Benjamin Moster Numerical Structure Formation 4 IoA,
14 Advantages and problems Advantages of SPH:! Numerically very robust! excellent conservation properties! Lagrangian nature trace flow & Galilean-invariant! Inherently adaptive (no ad-hoc refinement needed)! Couples properly to N-body gravity methods!! Problems?! SPH doesn t mix fluids (well)!,kaufmann-blobs form! Impacts accretion rates, SFRs disc sizes, gas fractions, etc Benjamin Moster Numerical Structure Formation 4 IoA,
15 The mixing problem P u ρ! Contact discontinuity in hydrostatic eq.: P = cst.! Left: u1 ρ1, Right: u2 ρ2 u1 = ρ2/ρ1 u2! SPH: smooth density gradient Energy not smoothed! Pressure blib introduced surface tension prevents mixing 12 Benjamin Moster Numerical Structure Formation 4 IoA,
16 Two gas clouds collide Density increases Energy decreases Density is smoothed Energy isn t smoothed Pressure blip Pressure minimum causes nearby particles to collapse blob The blob problem 13 Credit: A. Hobbs & J. Read Benjamin Moster Numerical Structure Formation 4 IoA,
17 Pressure-Entropy SPH Classic SPH: derive EoM from Lagrangian smoothing density:! h i Instead of using density, use pressure: P i = y i = h PN i j=1 m j A 1/ j W ij (h i )! # Automatically guarantees that pressure is single-valued removes pressure blip and surface tension! Equations will be well-behaved so long as pressure is smooth True by definition in contact discontinuities Not true at shocks, but neither is the density constant Do not lose any desirable behaviors of the density formulation 14 Benjamin Moster Numerical Structure Formation 4 IoA,
18 Physical Diffusion should be included in simulations (shouldn t include metal diffusion without energy diffusion! du/dt = P/ρ v! Diffusion removes the energy minimum pressure smooth! Surface tension vanishes Fluids can mix / no blobs form! However: Unclear how much diffusion to be added Artificial Diffusion 15 Benjamin Moster Numerical Structure Formation 4 IoA,
19 Kelvin Helmholtz Instability Clouds Saturn 16 Benjamin Moster Numerical Structure Formation 4 IoA,
20 Blob test Classic SPH Modern SPH 17 Benjamin Moster Numerical Structure Formation 4 IoA,
21 Isolated discs in hot halo Classic SPH Modern SPH 18 Benjamin Moster Numerical Structure Formation 4 IoA,
22 Outline of this lecture Smoothed particle hydrodynamics (SPH)! Idea & smoothing kernel! Equations of motion! Artificial viscosity! Advantages and problems (fluid-mixing, SPH blobs)! Modern SPH techniques (Pressure-entropy, artificial diffusion, etc.)! Baryonic sub-grid physics! Radiative Cooling and heating! Star formation! Stellar feedback (Multi-phase model, blastwave feedback)! AGN feedback 19 Benjamin Moster Numerical Structure Formation 4 IoA,
23 How to form a galaxy Up to now we have included Gravity (N-body), and Hydrodynamics! All other processed happen below the scale we can resolve! Include effective models to describe these physics:! Radiative cooling and heating! Star formation! Feedback from supernovae! Feedback from AGN! Feedback from other sources?! All processes are in general rather poorly understood 20 Benjamin Moster Numerical Structure Formation 4 IoA,
24 Radiative Cooling & Heating Stars can only form when gas is dense enough (i.e. cold enough)! Gas can cool via a number of mechanisms:! Inverse Compton scattering of CMB photons by electrons in the hot gas. Long cooling time (only important in the early Universe).! Excitation of rot. & vib. levels in molecular hydrogen through collisions, and subsequent decay (important in low mass haloes).! Photon emission after collisions between atoms and electrons, excited levels decay radiatively (important in int. mass haloes).! Bremsstrahlung: electrons are accelerated in an ionized plasma (important in massive clusters)! Gas is heated by a background of high energy UV photons from quasars and massive stars 21 Benjamin Moster Numerical Structure Formation 4 IoA,
25 Radiative Cooling & Heating Each particle (or cell) carries information on density and temperature! Need to compute abundances of the different ionic species of the gas Assume: gas optically thin, & Ionisation equilibrium! Compute cooling rate:!! = X i i (n i,t) Heating rate from photo- ionisation: H = X n i i i Katz Benjamin Moster Numerical Structure Formation 4 IoA,
26 Change in thermal energy With the cooling rate Λ and the heating rate H, we can compute the change in thermal energy:! du i dt = H i Up to now, only primordial cooling! In practice, cooling is computed on an element-by-element basis: where Λi is tabulated for nh, T, and z! Need to follow the creation of elements in SNe (usually 11 elements are tracked) i i = H,He + X i>he i 23 Benjamin Moster Numerical Structure Formation 4 IoA,
27 Star formation Stars should form in a gas cloud if: - the gas density is high enough (threshold density nmin) - the region is Jeans unstable! d dt = g t sf where ε is an efficiency and tsf is a characteristic timescale! The star formation rate is given by Often the dynamical timescale is used for t sf! Implemented in the code in a stochastic way: p = m g 1 e t/t sf m Draw random number r: if r < p create new star particle! No feedback gas can form stars on dynamical time overcooling 24 Benjamin Moster Numerical Structure Formation 4 IoA,
28 Stellar feedback: Multiphase model Springel & Hernquist (2003) developed a multi-phase model! Each gas particle is treated as a two-phase medium: cold and hot! Gas density is ρ = ρc + ρh with specific energy ε = εcρc + εhρh! d dt =(1 \epsilon_h where β is the fraction of short-lived stars and t* is a timescale! Star formation: cold clouds are converted to stars: ) c t SNe release energy and move gas from cold to hot phase. Cloud mass evaporated proportional to SN mass:! Cold clouds grow via radiative cooling d c dt RC = d h dt d c dt EV RC = A c t = net( h, h ) h c 25 Benjamin Moster Numerical Structure Formation 4 IoA,
29 Stellar feedback: Multiphase model The rates at which the hot and the cold phases evolve are:! d c dt =! d h! dt c t The star formation timescale t* is chosen to reproduce the observed Schmidt-Kennicutt relation: t ( ) =t 0 Free parameters are chosen such that - KS-relation reproduced A c + net( h, h ) t - threshold density ρth corresponds to observed one h c = c + A c net ( h, h ) t t th h 1/2 c 26 Benjamin Moster Numerical Structure Formation 4 IoA,
30 Galactic winds Feedback is still very inefficient (overcooling)! Introduce mechanical feedback given by galactic winds! The mass-loss rate is proportional to the SFR: η is called mass loading factor, observations suggest η = 2! Wind carries fixed fraction of SN energy: Set equal to the kinetic wind energy:! Particle added to the wind if random number is below p w =1 exp (1 ) c t t SN Ṁ Ṁ w = Ṁ v w = and! ~v 0 = ~v + v w ~n Numerical trick: decouple winds from hydro for short amount of time 27 Benjamin Moster Numerical Structure Formation 4 IoA, s 2 SN (1 )
31 Blastwave feedback Thermal feedback is usually quickly radiated away (in dense gas)! Stinson et al. (2006) implemented feedback based on a SN blast wave! In each SN explosion a fraction of the SN energy is added to the ISM! Cooling is temporarily switched of in neighboring particles within a radius of: for a time: R E = E 0.32 SN n 0.16 P 0.2 pc t E = E 0.31 SN n 0.27 P 0.64 yr based on work by Chevalier (1974) and McKee & Ostriker (1977)! Neighboring particles will become hot, reducing the SFR! High temperature leads to high pressure and drives a galactic wind 28 Benjamin Moster Numerical Structure Formation 4 IoA,
32 AGN feedback In massive haloes, SN feedback cannot reduce the SF efficiency Potential well too deep to drive massive winds! Energy released from accretion of gas onto BHs can be considerable! Relevant scales for the BH cannot be resolved, hence effective model:! Accretion rate given by Bondi-Hoyle rate: limited by Eddington rate.! Ṁ B = 4 G2 M 2 BH (c 2 s + v 2 ) 3/2 Some fraction εf of the radiated luminosity can couple thermally to surrounding gas: Ė FB = f L r = f r Ṁ BH c 2 (εf: M-E conversion)! Usually εf = % of the accreted rest mass energy is available 29 Benjamin Moster Numerical Structure Formation 4 IoA,
33 AGN feedback Implementation: BH particles are colissionless sink particles though compute gas density & temperature for each BH! For each gas particle around BH compute: p = wṁbh with w = kernel weight. If random number < p absorb gas particle in BH (including momentum)! t Add energy kernel-weighted to the thermal energy of the surrounding gas particles.! Merge BHs if they are within one smoothing length of each other and if the relative velocity is smaller than the sound speed 30 Benjamin Moster Numerical Structure Formation 4 IoA,
34 Results of feedback In low mass haloes SN feedback is very efficient! In massive haloes AGN feedback can reduce the conversion efficiency! To simulate realistic massive galaxies, a combination of both is needed 31 Benjamin Moster Numerical Structure Formation 4 IoA,
35 Text Books:! Final notes Cosmology: Galaxy Formation and Evolution (Mo, vdbosch, White)! Papers:! Hopkins (2013), MNRAS, 428, 2840! Hu et al. (2014), MNRAS, 443, 1173! Springel & Hernquist (2003), MNRAS, 339, 289! Stinson et al. (2006), MNRAS, 373, 1074! Springel et al. (2005), MNRAS, 62, 79! Gadget and N-GenIC website: 32 Benjamin Moster Numerical Structure Formation 4 IoA,
36 Up next Lecture 1: Motivation & Initial conditions! Lecture 2: Gravity algorithms & parallelization! Lecture 3: Hydro schemes! Lecture 4: Radiative cooling, photo heating & Subresolution physics! Lecture 5: Halo and subhalo finders & Semi-analytic models! Lecture 6: Getting started with Gadget! Lecture 7: Example cosmological box! Lecture 8: Example galaxy collision 33 Benjamin Moster Numerical Structure Formation 4 IoA,
37
38 Extra Slides
Numerical Cosmology & Galaxy Formation
Numerical Cosmology & Galaxy Formation Lecture 13: Example simulations Isolated galaxies, mergers & zooms Benjamin Moster 1 Outline of the lecture course Lecture 1: Motivation & Historical Overview Lecture
More informationSuperbubble Feedback in Galaxy Formation
Superbubble Feedback in Galaxy Formation Ben Keller (McMaster University) James Wadsley, Samantha Benincasa, Hugh Couchman Paper: astro-ph/1405.2625 (Accepted MNRAS) Keller, Wadsley, Benincasa & Couchman
More informationJoop Schaye (Leiden) (Yope Shea)
Overview of sub-grid models in cosmological simulations Joop Schaye (Leiden) (Yope Shea) Length Scales (cm) Subgrid models Cosmological simulations 8 0 2 11 18 20 22 24 28 interparticle distance in stars
More informationGalaxy Simulations Using the N-Body/SPH code GADGET
1 2010 HIPACC Astro-Computing Summer School Galaxy Simulations Using the N-Body/SPH code GADGET T.J. Cox (Carnegie Observatories) 2 Outline 1. Who am I and what am I doing here? My perspective, my science,
More informationMoving mesh cosmology: The hydrodynamics of galaxy formation
Moving mesh cosmology: The hydrodynamics of galaxy formation arxiv:1109.3468 Debora Sijacki, Hubble Fellow, ITC together with: Mark Vogelsberger, Dusan Keres, Paul Torrey Shy Genel, Dylan Nelson Volker
More informationGalaxy formation and evolution II. The physics of galaxy formation
Galaxy formation and evolution II. The physics of galaxy formation Gabriella De Lucia Astronomical Observatory of Trieste Outline: ü Observational properties of galaxies ü Galaxies and Cosmology ü Gas
More informationTheoretical ideas About Galaxy Wide Star Formation! Star Formation Efficiency!
Theoretical ideas About Galaxy Wide Star Formation Theoretical predictions are that galaxy formation is most efficient near a mass of 10 12 M based on analyses of supernova feedback and gas cooling times
More informationGalaxy Formation: Overview
Galaxy Formation: Overview Houjun Mo March 30, 2004 The basic picture Formation of dark matter halos. Gas cooling in dark matter halos Star formation in cold gas Evolution of the stellar populaion Metal
More informationCHAPTER 4. Basics of Fluid Dynamics
CHAPTER 4 Basics of Fluid Dynamics What is a fluid? A fluid is a substance that can flow, has no fixed shape, and offers little resistance to an external stress In a fluid the constituent particles (atoms,
More informationCosmological simulations of galaxy formation
Cosmological simulations of galaxy formation Modern galaxy formation simulations Mock gri SDSS composite image with dust absorption based on Draine opacity model. NGC4622 as seen from HST Outline - Impact
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, 11 26 / Colloquium University of Victoria Outline 1 Cosmic
More informationASTR 610 Theory of Galaxy Formation Lecture 17: Feedback
ASTR 610 Theory of Galaxy Formation Lecture 17: Feedback Frank van den Bosch Yale University, Fall 2018 Feedback In this lecture we discuss feedback, focussing mainly on supernove feedback. After describing
More informationDisk Formation and the Angular Momentum Problem. Presented by: Michael Solway
Disk Formation and the Angular Momentum Problem Presented by: Michael Solway Papers 1. Vitvitska, M. et al. 2002, The origin of angular momentum in dark matter halos, ApJ 581: 799-809 2. D Onghia, E. 2008,
More informationAstro-2: History of the Universe
Astro-2: History of the Universe Lecture 13; May 30 2013 Previously on astro-2 Energy and mass are equivalent through Einstein s equation and can be converted into each other (pair production and annihilations)
More informationAGN Feedback In an Isolated Elliptical Galaxy
AGN Feedback In an Isolated Elliptical Galaxy Feng Yuan Shanghai Astronomical Observatory, CAS Collaborators: Zhaoming Gan (SHAO) Jerry Ostriker (Princeton) Luca Ciotti (Bologna) Greg Novak (Paris) 2014.9.10;
More informationThe Initial Mass Function Elisa Chisari
The Initial Mass Function AST 541 Dec 4 2012 Outline The form of the IMF Summary of observations Ingredients of a complete model A Press Schechter model Numerical simulations Conclusions The form of the
More informationFeedback from Radiation Pressure during Galaxy Formation
Feedback from Radiation Pressure during Galaxy Formation Eliot Quataert (UC Berkeley) w/ Norm Murray, Jackson Debuhr, Phil Hopkins... Spitzer s view of Carina Outline Feedback: What is it good for? Absolutely
More informationBright Cluster Galaxy formation and the role of AGN feedback. Romain Teyssier
Bright Cluster Galaxy formation and the role of AGN feedback Romain Teyssier KITP 2011: Monster Inc. Romain Teyssier 1 Outline - Feedback and galaxy formation - The role of AGN feedback in Milky Way halos
More informationStar formation feedback in galaxy formation models. Yu Lu (KIPAC/Stanford)
Star formation feedback in galaxy formation models Yu Lu (KIPAC/Stanford) Overview Galaxies form in dark matter halos. CDM model predicts too many low-mass and high-mass halos. Feedback is needed to explain
More informationThe Illustris simulation: a new look at galaxy black hole co-evolution. Debora Sijacki IoA & KICC Cambridge
The Illustris simulation: a new look at galaxy black hole co-evolution Debora Sijacki IoA & KICC Cambridge LSS conference July 23 2015 Cosmological simulations of galaxy and structure formation Hierarchical
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 informationThe Iguaçu Lectures. Nonlinear Structure Formation: The growth of galaxies and larger scale structures
April 2006 The Iguaçu Lectures Nonlinear Structure Formation: The growth of galaxies and larger scale structures Simon White Max Planck Institute for Astrophysics z = 0 Dark Matter ROT EVOL Cluster structure
More informationKilling Dwarfs with Hot Pancakes. Frank C. van den Bosch (MPIA) with Houjun Mo, Xiaohu Yang & Neal Katz
Killing Dwarfs with Hot Pancakes Frank C. van den Bosch (MPIA) with Houjun Mo, Xiaohu Yang & Neal Katz The Paradigm... SN feedback AGN feedback The halo mass function is much steeper than luminosity function
More informationCoupling small and large scales: how massive and dying stars drive the formation of a galaxy
Coupling small and large scales: how massive and dying stars drive the formation of a galaxy P. Monaco, Trieste University and INAF-Osservatorio Astronomico di Trieste Micro-physics of star formation and
More informationMotivation Q: WHY IS STAR FORMATION SO INEFFICIENT? Ṁ M gas / dyn. Log SFR. Kennicutt Log. gas / dyn
Motivation Q: WHY IS STAR FORMATION SO INEFFICIENT? Ṁ 0.017 M gas / dyn Log SFR Kennicutt 1998 Log gas / dyn Motivation Q: WHY IS STAR FORMATION SO INEFFICIENT? Moster 2009 No Feedback 10% of baryons Log(
More informationAGN Feedback. Andrew King. Dept of Physics & Astronomy, University of Leicester Astronomical Institute, University of Amsterdam. Heidelberg, July 2014
AGN Feedback Andrew King Dept of Physics & Astronomy, University of Leicester Astronomical Institute, University of Amsterdam Heidelberg, July 2014 galaxy knows about central SBH mass velocity dispersion
More informationStellar evolution Part I of III Star formation
Stellar evolution Part I of III Star formation The interstellar medium (ISM) The space between the stars is not completely empty, but filled with very dilute gas and dust, producing some of the most beautiful
More informationChemodynamical Simulations Of the Universe & Elliptical Galaxies. Chiaki Kobayashi (Stromlo Fellow, RSAA, ANU)
Chemodynamical Simulations Of the Universe & Elliptical Galaxies Chiaki Kobayashi (Stromlo Fellow, RSAA, ANU) Chemodynamical Evolution AGN (negative & positive) Feedback? Gravity Hydrodynamics Star Formation?
More informationTowards Understanding Simulations of Galaxy Formation. Nigel Mitchell. On the Origin of Cores in Simulated Galaxy Clusters
Towards Understanding Simulations of Galaxy Formation Nigel Mitchell On the Origin of Cores in Simulated Galaxy Clusters Work published in the Monthly Notices of the Royal Astronomy Society Journal, 2009,
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 informationMulti-scale and multi-physics numerical models of galaxy formation
Multi-scale and multi-physics numerical models of galaxy formation M. Rieder and RT, 2016, MNRAS, 457, 1722 J. Rosdahl, J. Schaye, RT and O. Agertz, 2015, MNRAS, 451, 34 RAMSES: parallel Adaptive Mesh
More informationThe Impact of Quasar Feedback on the Formation & Evolution of Red Galaxies
The Impact of Quasar Feedback on the Formation & Evolution of Red Galaxies Philip Hopkins 07/17/06 Lars Hernquist, Volker Springel, Gordon Richards, T. J. Cox, Brant Robertson, Tiziana Di Matteo, Yuexing
More informationASTR 610 Theory of Galaxy Formation Lecture 14: Heating & Cooling
ASTR 610 Theory of Galaxy Formation Lecture 14: Heating & Cooling Frank van den Bosch Yale University, spring 2017 Heating & Cooling In this lecture we address heating and cooling of gas inside dark matter
More informationWhat Doesn t Quench Galaxy Formation?
What Doesn t Quench Galaxy Formation? Phil Hopkins Dusan Keres, Claude Faucher-Giguere, Jose Onorbe, Freeke van de Voort, Sasha Muratov, Xiangcheng Ma, Lena Murchikova, Norm Murray, Eliot Quataert, James
More informationFeedback and Galaxy Formation
Heating and Cooling in Galaxies and Clusters Garching August 2006 Feedback and Galaxy Formation Simon White Max Planck Institute for Astrophysics Cluster assembly in ΛCDM Gao et al 2004 'Concordance'
More informationFeedback, AGN and galaxy formation. Debora Sijacki
Feedback, AGN and galaxy formation Debora Sijacki Formation of black hole seeds: the big picture Planck data, 2013 (new results 2015) Formation of black hole seeds: the big picture CMB black body spectrum
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 informationMassive black hole formation in cosmological simulations
Institut d Astrophysique de Paris IAP - France Massive black hole formation in cosmological simulations Mélanie HABOUZIT Marta Volonteri In collaboration with Yohan Dubois Muhammed Latif Outline Project:
More information3 Observational Cosmology Evolution from the Big Bang Lecture 2
3 Observational Cosmology Evolution from the Big Bang Lecture 2 http://www.sr.bham.ac.uk/~smcgee/obscosmo/ Sean McGee smcgee@star.sr.bham.ac.uk http://www.star.sr.bham.ac.uk/~smcgee/obscosmo Nucleosynthesis
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 informationChapter 21 Galaxy Evolution. How do we observe the life histories of galaxies?
Chapter 21 Galaxy Evolution How do we observe the life histories of galaxies? Deep observations show us very distant galaxies as they were much earlier in time (old light from young galaxies). 1 Observing
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 informationModelling star formation in galaxy formation simulations
Modelling star formation in galaxy formation simulations Vadim Semenov (U.Chicago) Andrey Kravtsov University of Chicago Carving through the codes Davos, Switzerland 16 February, 2017 Nick Gnedin (Fermilab)
More informationBlack Hole Feedback. What is it? What does it do? Richard Bower Durham University
Black Hole Feedback What is it? What does it do? GALFORM: RGB + Benson, Lagos, Fanidakis, Frenk, Lacey, Baugh & Cole +++ EAGLE: Booth, Dalla Vecchia, Crain, Furlong, Rosas- Guevara, Schaye, RGB, Theuns,
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 informationSink particle accretion test
Sink particle accretion test David A. Hubber & Stefanie Walch 1 Objectives Simulate spherically-symmetric Bondi accretion onto a sink particle for an isothermal gas. Calculate the accretion rate onto a
More informationGas accretion in Galaxies
Massive Galaxies Over Cosmic Time 3, Tucson 11/2010 Gas accretion in Galaxies Dušan Kereš TAC, UC Berkeley Hubble Fellow Collaborators: Romeel Davé, Mark Fardal, C.-A. Faucher-Giguere, Lars Hernquist,
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 informationOrigin of Bi-modality
Origin of Bi-modality and Downsizing Avishai Dekel HU Jerusalem Galaxies and Structures Through Cosmic Times Venice, March 2006 Summary Q: z
More informationGravitational heating, clumps, overheating. Yuval Birnboim (Harvard Smithsonian Center for Astrophysics) Avishai Dekel (Hebrew University)
Gravitational heating, clumps, overheating Yuval Birnboim (Harvard Smithsonian Center for Astrophysics) Avishai Dekel (Hebrew University) Basic idea: Cooling flow Clusters need additional energy to reduce
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 informationCosmological simulations of our Universe Debora Sijacki IoA & KICC Cambridge
Cosmological simulations of our Universe Debora Sijacki IoA & KICC Cambridge Cosmo 17 Paris August 30th 2017 Cosmological simulations of galaxy and structure formation 2 Cosmological simulations of galaxy
More informationDynamics of Galaxies: Basics. Frontiers in Numerical Gravitational Astrophysics July 3, 2008
Dynamics of Galaxies: Basics Frontiers in Numerical Gravitational Astrophysics July 3, 2008 c = G = 1 h = 0 Outline (a) Galaxies as physical systems. (b) The Collisionless Boltzmann Equation. (c) N-body
More informationFluid Dynamics. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/14
Fluid Dynamics p.1/14 Fluid Dynamics Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/14 The equations of fluid dynamics are coupled PDEs that form an IVP (hyperbolic). Use
More informationPhysical Processes in Astrophysics
Physical Processes in Astrophysics Huirong Yan Uni Potsdam & Desy Email: hyan@mail.desy.de 1 Reference Books: Plasma Physics for Astrophysics, Russell M. Kulsrud (2005) The Physics of Astrophysics, Frank
More informationThe Birth Of Stars. How do stars form from the interstellar medium Where does star formation take place How do we induce star formation
Goals: The Birth Of Stars How do stars form from the interstellar medium Where does star formation take place How do we induce star formation Interstellar Medium Gas and dust between stars is the interstellar
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 informationThe importance of galactic fountain in galaxy evolution
The importance of galactic fountain in galaxy evolution Federico Marinacci in collaboration with: L. Armillotta, A. Marasco, F. Fraternali, J. Binney, L. Ciotti & C. Nipoti Sydney, 4 November 214 Galactic
More informationHydra A (x-ray) Perseus ` (optical)
The Biermann Lectures: Adventures in Theoretical Astrophysics I: The Physics of Galaxy Cluster Plasmas Eliot Quataert (UC Berkeley) w/ Mike McCourt, Ian Parrish, Prateek Sharma Perseus ` (optical) Hydra
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 informationObservational Evidence of AGN Feedback
10 de maio de 2012 Sumário Introduction AGN winds Galaxy outflows From the peak to the late evolution of AGN and quasars Mergers or secular evolution? The AGN feedback The interaction process between the
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 informationwhile the Planck mean opacity is defined by
PtII Astrophysics Lent, 2016 Physics of Astrophysics Example sheet 4 Radiation physics and feedback 1. Show that the recombination timescale for an ionised plasma of number density n is t rec 1/αn where
More informationOrianne ROOS CEA-Saclay Collaborators : F. Bournaud, J. Gabor, S. Juneau
Orianne ROOS CEA-Saclay Collaborators : F. Bournaud, J. Gabor, S. Juneau Bachelor of Physics, Master of Astrophysics Université de Strasbourg PhD, Université Paris-Diderot Observatoire de Strasbourg Les
More informationSubgrid Scale Physics in Galaxy Simulations
CRC 963 Astrophysical Turbulence and Flow Instabilities with thanks to Harald Braun, Jan Frederic Engels, Jens Niemeyer, IAG Ann Almgren and John Bell, LBNL Christoph Federrath, Monash University yt-project.org
More informationAccretion disks. AGN-7:HR-2007 p. 1. AGN-7:HR-2007 p. 2
Accretion disks AGN-7:HR-2007 p. 1 AGN-7:HR-2007 p. 2 1 Quantitative overview Gas orbits in nearly circular fashion Each gas element has a small inward motion due to viscous torques, resulting in an outward
More informationAST541 Lecture Notes: Galaxy Formation Dec, 2016
AST541 Lecture Notes: Galaxy Formation Dec, 2016 GalaxyFormation 1 The final topic is galaxy evolution. This is where galaxy meets cosmology. I will argue that while galaxy formation have to be understood
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 informationThe Effects of Radiative Transfer on Low-Mass Star Formation
The Effects of Radiative Transfer on Low-Mass Star Formation Stella Offner NSF Fellow, ITC Dense Cores in Dark Clouds Oct 23 2009 Collaborators: Chris McKee (UC Berkeley), Richard Klein (UC Berkeley; LLNL),
More informationSupernova Feedback in Low and High Mass Galaxies: Luke Hovey 10 December 2009
Supernova Feedback in Low and High Mass Galaxies: Luke Hovey 10 December 2009 Galactic Winds: Mathews, W. et al. 1971 Effects of Supernovae on the Early Evolution of Galaxies: Larson, R. 1974 The origin
More informationMoore et al Kenney et al. 2004
Moore et al. 1996 Kenney et al. 2004 (i) Interaction with other cluster members and/or with the cluster potential (ii) Interactions with the hot gas that permeates massive galaxy systems. The influence
More informationSunyaev-Zeldovich Contributions from Early Supernova Winds
Sunyaev-Zeldovich Contributions from Early Supernova Winds Sarah Benjamin Mentor: Steven Furlanetto Department of Physics and Astronomy, UCLA September 5, 2011 Abstract The purpose of this research was
More informationGravitational Radiation from Coalescing SMBH Binaries in a Hierarchical Galaxy Formation Model
Gravitational Radiation from Coalescing SMBH Binaries in a Hierarchical Galaxy Formation Model Motohiro ENOKI (National Astronomical Observatory of Japan) Kaiki Taro INOUE (Kinki University) Masahiro NAGASHIMA
More informationFLASH Code Tutorial. part III sink particles & feedback. Robi Banerjee Hamburger Sternwarte
FLASH Code Tutorial part III sink particles & feedback Robi Banerjee Hamburger Sternwarte banerjee@hs.uni-hamburg.de Motivation: modelling of dense regions in collapse simulations, e.g. star formation
More informationSchmidt-Kennicutt relations in SPH simulations of disc galaxies with effective SN thermal feedback
Schmidt-Kennicutt relations in SPH simulations of disc galaxies with effective SN thermal feedback Pierluigi Monaco Università di Trieste & INAF-Osservatorio Astronomico di Trieste In collaboration with:
More information8.1 Structure Formation: Introduction and the Growth of Density Perturbations
8.1 Structure Formation: Introduction and the Growth of Density Perturbations 1 Structure Formation and Evolution From this (Δρ/ρ ~ 10-6 ) to this (Δρ/ρ ~ 10 +2 ) to this (Δρ/ρ ~ 10 +6 ) 2 Origin of Structure
More informationChapter 16 Lecture. The Cosmic Perspective Seventh Edition. Star Birth Pearson Education, Inc.
Chapter 16 Lecture The Cosmic Perspective Seventh Edition Star Birth 2014 Pearson Education, Inc. Star Birth The dust and gas between the star in our galaxy is referred to as the Interstellar medium (ISM).
More informationGalaxy Evolution. Part 4. Jochen Liske Hamburger Sternwarte
Galaxy Evolution Part 4 Jochen Liske Hamburger Sternwarte jochen.liske@uni-hamburg.de Astronomical picture of the week NGC 1275 Central galaxy of the Perseus cluster Active Source: APOD Astronomical picture
More informationThe IGM in simulations
The IGM in simulations Institute for Computational Cosmology Ogden Centre for Fundamental Physics Durham University, UK and University of Antwerp Belgium 1 Menu: Simulated DLAs: column density and dynamics
More informationThe Black Hole in the Galactic Center. Eliot Quataert (UC Berkeley)
The Black Hole in the Galactic Center Eliot Quataert (UC Berkeley) Why focus on the Galactic Center? The Best Evidence for a BH: M 3.6 10 6 M (M = mass of sun) It s s close! only ~ 10 55 Planck Lengths
More informationHigh-Energy Astrophysics Lecture 6: Black holes in galaxies and the fundamentals of accretion. Overview
High-Energy Astrophysics Lecture 6: Black holes in galaxies and the fundamentals of accretion Robert Laing Overview Evidence for black holes in galaxies and techniques for estimating their mass Simple
More informationLocal photo-ionization radiation, Circum-galactic gas cooling and galaxy formation
Local photo-ionization radiation, Circum-galactic gas cooling and galaxy formation or A critical Star-Formation-Rate divides hot-mode from cold-mode accretion Sebastiano Cantalupo IMPS Fellow, UCSC Chandra
More informationWeek 8: Stellar winds So far, we have been discussing stars as though they have constant masses throughout their lifetimes. On the other hand, toward
Week 8: Stellar winds So far, we have been discussing stars as though they have constant masses throughout their lifetimes. On the other hand, toward the end of the discussion of what happens for post-main
More informationAGN Feedback. Andrew King. Dept of Physics and Astronomy, University of Leicester, UK Anton Pannekoek Institute, Amsterdam, NL Sterrewacht Leiden, NL
AGN Feedback Andrew King Dept of Physics and Astronomy, University of Leicester, UK Anton Pannekoek Institute, Amsterdam, NL Sterrewacht Leiden, NL Potsdam, 2018 SMBH affects galaxy bulge more than enough
More informationPATRICIA B. TISSERA. Institute for Astronomy and Space Physics Argentina
PATRICIA B. TISSERA Institute for Astronomy and Space Physics Argentina i Observational motivations Model Results GALAXY FORMATION In the last years, studies of chemical elements obtained in the Local
More informationCold Clouds in Cool Cores
Becky Canning NGC 1275, Perseus Cluster Cold Clouds in Cool Cores 1 Snowcluster 2015 Cold Clouds in Cool Cores Excitation mechanisms of the multi-phase cool/cold gas Why should we care? Heard that X-ray
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 informationChapter 11 The Formation and Structure of Stars
Chapter 11 The Formation and Structure of Stars Guidepost The last chapter introduced you to the gas and dust between the stars that are raw material for new stars. Here you will begin putting together
More informationFeedback from growth of supermassive black holes
Research Collection Other Conference Item Feedback from growth of supermassive black holes Author(s): Begelman, Mitchell C.; Ruszkowksi, Mateusz Publication Date: 2003 Permanent Link: https://doi.org/10.3929/ethz-a-004585094
More informationBlack Holes in the Early Universe Accretion and Feedback
1 1 Black Holes in the Early Universe Accretion and Feedback 1 1 Black Holes in the Early Universe Accretion and Feedback Geoff Bicknell & Alex Wagner Australian National University 1 1 High redshift radio
More informationBimodal regime in young massive clusters leading to formation of subsequent stellar generations
Bimodal regime in young massive clusters leading to formation of subsequent stellar generations Richard Wünsch J. Palouš, G. Tenorio-Tagle, C. Muñoz-Tuñón, S. Ehlerová Astronomical institute, Czech Academy
More informationImplementing sub-grid treatments of galactic outflows into cosmological simulations. Hugo Martel Université Laval
Implementing sub-grid treatments of galactic outflows into cosmological simulations Hugo Martel Université Laval Leiden, June 19, 2013 GALACTIC OUTFLOWS Optical image of galaxy (Hubble Space Telescope)
More informationThe Interstellar Medium
The Interstellar Medium Fall 2014 Lecturer: Dr. Paul van der Werf Oortgebouw 565, ext 5883 pvdwerf@strw.leidenuniv.nl Assistant: Kirstin Doney Huygenslaboratorium 528 doney@strw.leidenuniv.nl Class Schedule
More informationhigh density low density Rayleigh-Taylor Test: High density medium starts on top of low density medium and they mix (oil+vinegar) Springel (2010)
GAS MIXES high density Springel (2010) low density Rayleigh-Taylor Test: High density medium starts on top of low density medium and they mix (oil+vinegar) HOT HALO highest resolved density nth= 50x10
More informationOld pre-main-sequence stars
Old pre-main-sequence stars and a second chance for planet formation Peter Scicluna Uni Kiel & ESO, Garching 9th September 2014 Giovanni Rosotti (MPE, USM, IoA) Leonardo Testi (ESO, INAF) Jim Dale (Excellence
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 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 informationASTR 610 Theory of Galaxy Formation Lecture 16: Star Formation
ASTR 610 Theory of Galaxy Formation Lecture 16: Star Formation Frank van den Bosch Yale University, spring 2017 Star Formation In this lecture we discuss the formation of stars. After describing the structure
More informationThe Monster Roars: AGN Feedback & Co-Evolution with Galaxies
The Monster Roars: AGN Feedback & Co-Evolution with Galaxies Philip Hopkins Ø (Nearly?) Every massive galaxy hosts a supermassive black hole Ø Mass accreted in ~couple bright quasar phase(s) (Soltan, Salucci+,
More informationAstrofysikaliska Dynamiska Processer
Astrofysikaliska Dynamiska Processer VT 2008 Susanne Höfner hoefner@astro.uu.se Aims of this Course - understanding the role and nature of dynamical processes in astrophysical contexts and how to study
More information