N-body Cosmology: The billion-body problem. Brian O Shea Michigan State University
|
|
- Diane Owens
- 5 years ago
- Views:
Transcription
1 N-body Cosmology: The billion-body problem Brian O Shea (oshea@msu.edu) Michigan State University TIARA Winter School, Jan. 2011
2 This lecture Goals of cosmological N-body simulations Overview of an N-body simulation Requirements and assumptions Equations solved Numerical techniques for N-body simulations Generating cosmological initial conditions Some example applications Useful references and codes
3 Goals of cosmological N-body sims Model the growth of structure in the universe Can be treated as an experiment to verify theories of the origin, evolution of the universe More specific applications: galaxy cluster mass function, evolution of dark energy EOS, substructure in DM halos, base for semi-analytic models of galaxy formation/evolution If other physics is included: galaxy formation, galaxy cluster formation
4 Overview of a cosmological N- body simulation Generate initial conditions at high redshift using matter power spectrum and perturbation theory Evolve simulation forward in time: calculate forces acting on particles calculate max. allowable timestep advance all particles Done iteratively! Compare to observations!
5 Requirements and assumptions We calculate the evolution of matter in an expanding universe! For this, we need and assume: Comoving coordinates (and accompanying cosmological model) Periodic boundary conditions Gravity is the only force considered Dark matter can be treated in a statistical way as a collection of mobile, discrete particles (so we can use an N-body method rather than solving the collisionless Boltzmann equation directly) Other physical processes (particularly gas, stars and related physics) are not important
6 An aside on cosmology WMAP Year 7 map of the microwave sky
7 SDSS galaxy map All galaxies from SDSS survey with declination of o
8 WMAP 7 matter power spectrum Larson et al. 2011, ApJS, 192, 16
9 combined P(k) as a function of wave number (c/o Max Tegmark)
10 Komatsu et al. 2011, ApJS, 192, 18
11 Equations
12 Equations of motion But, you must include comoving coordinates:
13 Where a(t) is calculated assuming a homogeneous, isotropic universe and a Robertson-Walker metric, and reduces to the Friedmann equations: Often, calculate p assuming a perfect fluid:
14 In comoving coordinates, equations of motion and Poisson equation become:
15 Numerical techniques
16 Orbital integration It is not necessarily desirable to calculate a very high-order position update (eighth order Runge Kutte, for example): Gravitational potential for N-body systems very noisy! A particle is not really intended to be a single object: pairwise interactions undesirable Keeping many intermediate steps is memory-intensive and computationally costly Essentially all cosmology N-body codes use a leapfrog integrator:
17 Time steps in cosmological simulations Gravity causes particles to clump. The relevant timescale is dynamical time! Choice of time step is critical: Too small: lots of wasted CPU time Too large: excessive errors in integration! Not all particles need to take the same time step, either! In general, time steps in cosmological simulations governed by local forces and resolution.
18 A brief aside: hydrodynamics Gas is important in astrophysics, and we ignore it at our peril! BUT, large scale structure formation is dominated by gravity: can ignore hydrodynamics for certain problems! Can solve equations of hydrodynamics using either Eulerian (grid-based) or Lagrangian (particle-based) formalism: the choice of formalism depends on many considerations. Different types of N-body solvers favor either Eulerian or Lagrangian hydro solutions: I will attempt to specify.
19 N-body techniques
20 Direct summation Discussed previously by Douglas Heggie Sum forces by direct particle-to-particle gravitational force:
21 Direct summation Advantage: Exact forces calculated on each particle! Disadvantage: Takes N 2 summations to calculate force in one timestep Difficult to apply periodic boundary conditions Doable for N ~ 10 6, impossible for N ~ 10 9
22 The first N-body simulation 37 particles in 2D: Holmberg 1941, ApJ, 94, 385
23 Integration method: analog (light bulbs, photocells, and graduate students) 37 particles in 2D: Holmberg 1941, ApJ, 94, 385
24 Results: tidal deformation of disk nebulae, also interacting galaxies tend to merge Figure 4b in Holmberg 1941, ApJ, 94, 385
25 Tree method Space divided recursively into a hierarchy of cells Particles see neighbors directly: distant groups of particles are aggregated into a single apparent object with mass equivalent to the sum of the particle masses Original reference: Barnes & Hut 1986, Nature, 324,
26 Tree method Advantages O(N log N) operations Force errors can be easily understood and reduced Relatively easy to implement Spatially adaptive - good force resolution Disadvantages Large amount of memory required Periodic boundary conditions challenging
27 Particle-Mesh Uses Cartesian grid to calculate forces on particles. For each timestep: 1. Compute mass density on grid 2. Solve Poisson s equation for gravitational potential using gridded density 3. Calculate force on particles using gridded potential See Computer Simulations Using Particles by Hockney and Eastwood
28 Particle-Mesh Advantages Very fast: O(Np) + O(Ng log Ng) operations using FFT Periodic boundary conditions trivial Very simple to write and parallelize Disadvantages Poor approximation of gravitational force for nearby particles: resolution ~ 2Δx Some force anisotropy due to grid
29 APM (Adaptive Particle-Mesh) Improves particle-mesh algorithm by placing subgrids in places where particles are clustered: better force resolution! Potential on subgrids calculated using BCs from coarse grid, density field from fine grid, and a multigrid relaxation method This is the basis for most adaptive mesh refinement code gravity solvers: Enzo, FLASH, RAMSES, MLAPM Jessop et al. 1994, J. Comp. Phys., 115, Kravtsov et al. 1997, ApJS, 111, 73-94
30 P 3 M Combines particle-particle and particle-mesh methods At large distances, use particle-mesh to calculate gravitational potentials Up close, force computed by direct summation Advantages: better approximation of gravitational force for close particles than PM faster for given Np than PP Disadvantage: becomes slow when clustering occurs (cost of direct-summation dominates) See Efstathiou et al. 1985, ApJS, 57,
31 AP 3 M Improves the P3M method by using additional subgrids at regions where clustering occurs (like APM, but with additional particle-particle summation) Advantages: Good accuracy of pair forces on almost all pairs Very fast: O(Np log Np) Disadvantages: limited dynamic range in subgrids some ambiguity in subgrid placement (spurious forces) still eventually run into particle-particle expense See Couchman et al. 1991, ApJL, 368, L23-26
32 Example AP 3 M simulation. From Couchman et al. 1995, ApJ, 452, 797
33 TreePM Synthesizes tree algorithm and particle-mesh algorithm: At large distances, use PM to calculate smooth contribution to gravitational potential At small distances (< 2dx) use tree to calculate local neighbors contributions Very fast and efficient, high dynamic range, excellent scaling: current state-of-the-art for pure N-body Example code: Gadget-2 (Springel 2005, MNRAS, 364, 1105; Original method papers: Xu 1995 (ApJS, 98, 355), Bode & Ostriker 2003 (ApJS, 145, 1)
34 Evolution of simulation (and computer) size with time
35 2009: (68.7 billion) particles Teyssier et al. 2010, A&A, 497, : multiple simulations (Habib, Heitmann et al., in prep) 2011: particles, Norman, O Shea et al. (in prep.), From Springel et al. 2005, Nature, 435, 629
36 From top500.org
37 From top500.org
38 Initial conditions
39 Background ICs for cosmological simulations are generated at high redshift (z ~ ), when the universe is approximately linear (density variations of ~5-10%) Inputs: cosmological parameters and matter power spectrum Grav. potential calculated on a lattice using a Gaussian random field Particles are typically initialized on a lattice (or in a glass force-free initial condition) and given spatial, velocity perturbations using the Zeldovich approximation (or something similar)
40 Gaussian random field Given a matter power spectrum, P(k), and a lattice with (Nx, Ny, Nz) cells and known volume. Sample the power spectrum discretely at each grid location, for a given k: At each position, the k-space fluctuation is given random phase and (complex) fluctuation amplitude: We then perform a Fourier transform to get the relative matter density fluctuation at each spatial grid point in the simulation volume
41 Particles start out distributed uniformly on a grid, or in a forcefree ( glass ) configuration. Then, they are displaced using the Zel dovich approximation (Zel dovich 1970, A&A, 5, 84-89): q = unperturbed position D(a) = growth factor H = Hubble parameter f = logarithmic growth rate of D(a) = irrotational displacement field Note: the Zel dovich approximation has incorrect second- and higher-order growth rates: it s worth considering higher-order approximation methods if you are doing something particularly sensitive to this (Crocce et al. 2006, MNRAS, 373, 369)
42 Example ICs (dark matter) Log dark matter density z=40 FOV: 32 Mpc/h comoving Depth: 8 Mpc/h comoving simulation has 256^3 (16.8 million) particles c/o Bode, Cen, Xu
43 Example ICs (dark matter) Log dark matter density z=0 FOV: 32 Mpc/h comoving Depth: 8 Mpc/h comoving simulation has 256^3 (16.8 million) particles c/o Bode, Cen, Xu
44 Some example applications
45 Formation of filaments in CDM Movie c/o Andrey Kravtsov Box 40 Mpc across
46 Formation of a galaxy group Zoom-in of previous sim. c/o Andrey Kravtsov
47 Next slide: fly-through of the Millennium Simulation (Springel et al. 2005, Nature, 435, )
48
49 Springel et al. 2005, Nature, 435,
50 Springel et al. 2005, Nature, 435,
51 The Aquarius Project (Virgo collaboration) Formation of MW-halo mass galaxy
52 Growth of six different MW-sized halos Boylan-Kolchin et al. 2010, MNRAS, 406, 896 Formation of metalpoor globular clusters as f(z) (blue: halos that form stars) Griffen et al. 2010, MNRAS, 405, 375
53 Formation of a Milky Way-type galaxy
54 Galaxy formation with a supermassive black hole Springel, Hopkins et al.
55 Galaxy merger with small (1:8) companion (retrograde) Younger et al. 2007, ApJ, 670, 269
56 Useful references and codes
57 References, I Computer Simulations Using Particles, by Hockney & Eastwood (Taylor & Francis, 1989) by Aarseth (Cambridge University Press, 2003) Gravitational N-body Simulations: Tools and Algorithms, Simulations of Structure Formation in the Universe, Bertschinger 1998, Annual Reviews of Astronomy and Astrophyscis, 36, Numerical Simulations in Cosmology, I: Methods, Klypin, astro-ph/ (also see papers II and III in this series) A survey of all known N-body methods:
58 References, II Institute for Advanced Studies Computational Astrophysics Summer School lectures (particularly those by V. Springel): University of California HIPACC Astro-computing summer school (part. lectures by Klypin): SummerSchool_archive.html Cosmological N-body simulation: Techniques, scope and status, Bagla 2005, Current Science, 88, Galactic Dynamics, by Binney & Tremaine (Princeton University Press, 1988)
59 Codes you can download and use Enzo: (PM/APM) Gadget: (Tree/TreePM) FLASH: (APM) RAMSES: TPM: (APM) (TreePM) Hydra: (AP 3 M)
Brad Gibson Centre for Astrophysics & Supercomputing Swinburne University
N-body Simulations: Tree and Mesh Approaches Brad Gibson Centre for Astrophysics & Supercomputing Swinburne University Demystifying the Jargon Trees Collisionless Boltzmann Equation Symplectic Time Integration
More informationComputational Astrophysics 1 Particle Mesh methods
Computational Astrophysics 1 Particle Mesh methods Oscar Agertz Outline - The Vlassov-Poisson equation - The Particle-Mesh method in a nut shell - Symplectic Integrators - Mass assignment and force interpolation
More informationHigh performance computing and numerical modeling
High performance computing and numerical modeling Volker Springel Plan for my lectures Lecture 1: Collisional and collisionless N-body dynamics Lecture 2: Gravitational force calculation Lecture 3: Basic
More informationNumerical 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 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 informationNumerical Cosmology & Galaxy Formation
Numerical Cosmology & Galaxy Formation Lecture 3: Generating initial conditions Benjamin Moster 1 Outline of the lecture course Lecture 1: Motivation & Historical Overview Lecture 2: Review of Cosmology
More informationN-body Simulations. Initial conditions: What kind of Dark Matter? How much Dark Matter? Initial density fluctuations P(k) GRAVITY
N-body Simulations N-body Simulations N-body Simulations Initial conditions: What kind of Dark Matter? How much Dark Matter? Initial density fluctuations P(k) GRAVITY Final distribution of dark matter.
More informationSimulations of structure-formation in theories beyond General Relativity
Simulations of structure-formation in theories beyond General Relativity Hans A. Winther Diving into the Dark, Cairns 21 July 2016 Hans Winther University of Oxford 2016 Motivation Why modified gravity?
More informationSimulations of Cosmic Structure Formation
Simulations of Cosmic Structure Formation Volker Springel The role of simulations in cosmology High-resolution N-body simulations Millennium XXL Hydrodynamic simulations and recent results for galaxy formation
More informationSwift: task-based hydrodynamics at Durham s IPCC. Bower
Swift: task-based hydrodynamics at Durham s IPCC Gonnet Schaller Chalk movie: Richard Bower (Durham) For the cosmological simulations of the formation of galaxies Bower Institute for Computational Cosmology
More informationPrinceton December 2009 The fine-scale structure of dark matter halos
Princeton December 2009 The fine-scale structure of dark matter halos Simon White Max Planck Institute for Astrophysics The dark matter structure of CDM halos A rich galaxy cluster halo Springel et al
More informationAstronomy 233 Winter 2009 Physical Cosmology Week 10 Simulations and Semi-Analytic Models Joel Primack University of California, Santa Cruz
Astronomy 233 Winter 2009 Physical Cosmology Week 10 Simulations and Semi-Analytic Models Joel Primack University of California, Santa Cruz Cosmological Simulation Methods Dissipationless Simulations Particle-Particle
More informationformation of the cosmic large-scale structure
formation of the cosmic large-scale structure Heraeus summer school on cosmology, Heidelberg 2013 Centre for Astronomy Fakultät für Physik und Astronomie, Universität Heidelberg August 23, 2013 outline
More informationSimulating the Universe
Simulating the Universe Christine Corbett Moran, Irshad Mohammed, Manuel Rabold, Davide Martizzi, Doug Potter, Aurel Schneider Oliver Hahn, Ben Moore, Joachim Stadel Outline - N-body codes: where do we
More informationStructure formation. Yvonne Y. Y. Wong Max-Planck-Institut für Physik, München
Structure formation Yvonne Y. Y. Wong Max-Planck-Institut für Physik, München Structure formation... Random density fluctuations, grow via gravitational instability galaxies, clusters, etc. Initial perturbations
More informationA Universe in Motion: Testing the Cosmological Paradigm with Galaxy Dynamics. John Dubinski, Toronto
A Universe in Motion: Testing the Cosmological Paradigm with Galaxy Dynamics John Dubinski, Toronto Outline Cosmology and galaxy dynamics Tools of the trade and computational challenges Case Studies: Triaxial
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 informationGalaxies 626. Lecture 3: From the CMBR to the first star
Galaxies 626 Lecture 3: From the CMBR to the first star Galaxies 626 Firstly, some very brief cosmology for background and notation: Summary: Foundations of Cosmology 1. Universe is homogenous and isotropic
More informationGalaxies & Introduction to Cosmology
Galaxies & Introduction to Cosmology Other Galaxies: How many are there? Hubble Deep Field Project 100 hour exposures over 10 days Covered an area of the sky about 1/100 the size of the full moon Probably
More informationGalaxy Formation! Lecture Seven: Galaxy Formation! Cosmic History. Big Bang! time! present! ...fluctuations to galaxies!
Galaxy Formation Lecture Seven: Why is the universe populated by galaxies, rather than a uniform sea of stars? Galaxy Formation...fluctuations to galaxies Why are most stars in galaxies with luminosities
More informationWriting a PM code. Andrey Kravtsov astro.uchicago.edu/ andrey/talks/pm/pm.pdf. February 20, 2002
Writing a PM code Andrey Kravtsov andrey@oddjob.uchicago.edu astro.uchicago.edu/ andrey/talks/pm/pm.pdf February 20, 2002 WRITING A PM CODE 1 Particle-Mesh (PM) technique: a brief history WRITING A PM
More informationMapping Hot Gas in the Universe using the Sunyaev-Zeldovich Effect
Mapping Hot Gas in the Universe using the Sunyaev-Zeldovich Effect Eiichiro Komatsu (Max-Planck-Institut für Astrophysik) Probing Fundamental Physics with CMB Spectral Distortions, CERN March 12, 2018
More informationTree Particle Mesh: an adaptive, efficient, and parallel code for collisionless cosmological simulation
to appear in ApJS v145, March 2003 Tree Particle Mesh: an adaptive, efficient, and parallel code for collisionless cosmological simulation Paul Bode and Jeremiah P. Ostriker Princeton University Observatory,
More informationA Comparison of AMR and SPH cosmology codes
Chapter 3 A Comparison of AMR and SPH cosmology codes 3.1 Summary We compare two cosmological hydrodynamic simulation codes in the context of hierarchical galaxy formation: the Lagrangian smoothed particle
More information21cmFAST A Fast, Semi-Numerical Simulation of the High-Redshift 21cm Signal
21cmFAST A Fast, Semi-Numerical Simulation of the High-Redshift 21cm Signal Mesinger, Furlanetto, & Cen (2010) Andrei Mesinger Princeton University Motivation We know next to nothing about high-z --> ENORMOUS
More informationGalaxy Bias from the Bispectrum: A New Method
Galaxy Bias from the Bispectrum: A New Method Jennifer E. Pollack University of Bonn Thesis Supervisor: Cristiano Porciani Collaborator: Robert E. Smith Motivation 1) Main Objective of state-of-the-art
More informationThe Millennium Simulation: cosmic evolution in a supercomputer. Simon White Max Planck Institute for Astrophysics
The Millennium Simulation: cosmic evolution in a supercomputer Simon White Max Planck Institute for Astrophysics The COBE satellite (1989-1993) Two instruments made maps of the whole sky in microwaves
More informationCOLA with scale dependent growth: applications to modified gravity and massive neutrinos
COLA with scale dependent growth: applications to modified gravity and massive neutrinos Kazuya Koyama Institute of Cosmology and Gravitation, University of Portsmouth Curvature Psaltis Living Rev. Relativity
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 informationCluster Formation in Lagrangian Perturbation Theory. Lavinia Heisenberg
Cluster Formation in Lagrangian Perturbation Theory Zentrum fu r Astronomie (ITA/ZAH) Universita t Heidelberg 18.12.2009, University of Canterbury, Christchurch outline 1 Introduction 2 LPT 3 ELL 4 PINOCCHIO
More informationThe galaxy population in cold and warm dark matter cosmologies
The galaxy population in cold and warm dark matter cosmologies Lan Wang National Astronomical Observatories, CAS Collaborators: Violeta Gonzalez-Perez, Lizhi Xie, Andrew Cooper, Carlos Frenk, Liang Gao,
More informationActive Galaxies and Galactic Structure Lecture 22 April 18th
Active Galaxies and Galactic Structure Lecture 22 April 18th FINAL Wednesday 5/9/2018 6-8 pm 100 questions, with ~20-30% based on material covered since test 3. Do not miss the final! Extra Credit: Thursday
More informationHunting for dark matter in the forest (astrophysical constraints on warm dark matter)
Hunting for dark matter in the forest (astrophysical constraints on warm dark matter) ICC, Durham! with the Eagle collaboration: J Schaye (Leiden), R Crain (Liverpool), R Bower, C Frenk, & M Schaller (ICC)
More informationarxiv: v1 [astro-ph.co] 1 Oct 2009
Mon. Not. R. Astron. Soc. 000, 000 000 (0000) Printed 10 November 2018 (MN LATEX style file v2.2) Second-order Lagrangian perturbation theory initial conditions for resimulations. Adrian Jenkins Institute
More informationTHE BOLSHOI COSMOLOGICAL SIMULATIONS AND THEIR IMPLICATIONS
GALAXY FORMATION - Durham -18 July 2011 THE BOLSHOI COSMOLOGICAL SIMULATIONS AND THEIR IMPLICATIONS JOEL PRIMACK, UCSC ΛCDM Cosmological Parameters for Bolshoi and BigBolshoi Halo Mass Function is 10x
More informationThe imprint of the initial conditions on large-scale structure
Stars, Planets and Galaxies 2018 Harnack House, Berlin The imprint of the initial conditions on large-scale structure Simon White Max Planck Institute for Astrophysics The Planck map of TCMB the initial
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 informationThe peculiar velocity and temperature profile of galaxy clusters
RAA 2014 Vol. 14 No. 6, 667 672 doi: 10.1088/1674 4527/14/6/005 http://www.raa-journal.org http://www.iop.org/journals/raa Research in Astronomy and Astrophysics The peculiar velocity and temperature profile
More informationA brief history of cosmological ideas
A brief history of cosmological ideas Cosmology: Science concerned with the origin and evolution of the universe, using the laws of physics. Cosmological principle: Our place in the universe is not special
More informationN-body Simulations and Dark energy
N-Body Simulations and models of Dark Energy Elise Jennings Supported by a Marie Curie Early Stage Training Fellowship N-body Simulations and Dark energy elise jennings Introduction N-Body simulations
More informationSome useful spherically symmetric profiles used to model galaxies
Some useful spherically symmetric profiles used to model galaxies de Vaucouleurs R 1/4 law for ellipticals and bulges A good fit to the light profile of many ellipticals and bulges: (constant such that
More informationCosmology and Astrophysics with Galaxy Clusters Recent Advances and Future Challenges
Cosmology and Astrophysics with Galaxy Clusters Recent Advances and Future Challenges Daisuke Nagai Yale University IPMU, July 15 th, 2010 Large-scale structure in the Universe SDSS (optical) Today δρ/ρ>>1
More informationAstronomy 114. Lecture35:TheBigBang. Martin D. Weinberg. UMass/Astronomy Department
Astronomy 114 Lecture35:TheBigBang Martin D. Weinberg weinberg@astro.umass.edu UMass/Astronomy Department A114: Lecture 35 09 May 2005 Read: Ch. 28,29 Astronomy 114 1/18 Announcements PS#8 due Monday!
More informationPhysical Cosmology 12/5/2017
Physical Cosmology 12/5/2017 Alessandro Melchiorri alessandro.melchiorri@roma1.infn.it slides can be found here: oberon.roma1.infn.it/alessandro/cosmo2017 Structure Formation Until now we have assumed
More informationLecture 37 Cosmology [not on exam] January 16b, 2014
1 Lecture 37 Cosmology [not on exam] January 16b, 2014 2 Structure of the Universe Does clustering of galaxies go on forever? Looked at very narrow regions of space to far distances. On large scales the
More informationDark Matter Detection Using Pulsar Timing
Dark Matter Detection Using Pulsar Timing ABSTRACT An observation program for detecting and studying dark matter subhalos in our galaxy is propsed. The gravitational field of a massive object distorts
More informationAdvanced Topics on Astrophysics: Lectures on dark matter
Advanced Topics on Astrophysics: Lectures on dark matter Jesús Zavala Franco e-mail: jzavalaf@uwaterloo.ca UW, Department of Physics and Astronomy, office: PHY 208C, ext. 38400 Perimeter Institute for
More informationMassive star clusters
Massive star clusters as a host of compact binaries Michiko Fujii ( 藤井通子 ) The University of Tokyo Outline Massive star clusters and compact binaries Dynamical evolution of star clusters Distribution of
More informationComments on the size of the simulation box in cosmological N-body simulations
Mon. Not. R. Astron. Soc. 358, 1076 1082 (2005) doi:10.1111/j.1365-2966.2005.08858.x Comments on the size of the simulation box in cosmological N-body simulations J. S. Bagla and Suryadeep Ray Harish-Chandra
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 informationExam 3 Astronomy 114
Exam 3 Astronomy 114 Select the answer that is the most appropriate among the choices given. 1. What is the Hubble law? (A) a relation between a galaxy s mass and radius. (B) a rule that describes the
More informationCurrent status of the ΛCDM structure formation model. Simon White Max Planck Institut für Astrophysik
Current status of the ΛCDM structure formation model Simon White Max Planck Institut für Astrophysik The idea that DM might be a neutral, weakly interacting particle took off around 1980, following a measurement
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 informationX-ray and Sunyaev-Zel dovich Effect cluster scaling relations: numerical simulations vs. observations
X-ray and Sunyaev-Zel dovich Effect cluster scaling relations: numerical simulations vs. observations Daisuke Nagai Theoretical Astrophysics, California Institute of Technology, Mail Code 130-33, Pasadena,
More informationVeilleux! see MBW ! 23! 24!
Veilleux! see MBW 10.4.3! 23! 24! MBW pg 488-491! 25! But simple closed-box model works well for bulge of Milky Way! Outflow and/or accretion is needed to explain!!!metallicity distribution of stars in
More information2. What are the largest objects that could have formed so far? 3. How do the cosmological parameters influence structure formation?
Einführung in die beobachtungsorientierte Kosmologie I / Introduction to observational Cosmology I LMU WS 2009/10 Rene Fassbender, MPE Tel: 30000-3319, rfassben@mpe.mpg.de 1. Cosmological Principles, Newtonian
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 informationPhysics of the Large Scale Structure. Pengjie Zhang. Department of Astronomy Shanghai Jiao Tong University
1 Physics of the Large Scale Structure Pengjie Zhang Department of Astronomy Shanghai Jiao Tong University The observed galaxy distribution of the nearby universe Observer 0.7 billion lys The observed
More informationLecture 05. Cosmology. Part I
Cosmology Part I What is Cosmology Cosmology is the study of the universe as a whole It asks the biggest questions in nature What is the content of the universe: Today? Long ago? In the far future? How
More informationStellar Dynamics and Structure of Galaxies
Stellar Dynamics and Structure of Galaxies Gerry Gilmore H47 email: gil@ast.cam.ac.uk Lectures: Monday 12:10-13:00 Wednesday 11:15-12:05 Friday 12:10-13:00 Books: Binney & Tremaine Galactic Dynamics Princeton
More informationAstronomical Computer Simulations. Aaron Smith
Astronomical Computer Simulations Aaron Smith 1 1. The Purpose and History of Astronomical Computer Simulations 2. Algorithms 3. Systems/Architectures 4. Simulation/Projects 2 The Purpose of Astronomical
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 informationarxiv:astro-ph/ v1 26 Nov 2004
Cosmological N-Body Simulations arxiv:astro-ph/0411730v1 26 Nov 2004 J.S.Bagla and T.Padmanabhan Inter-University Centre for Astronomy and Astrophysics, Post Bag 4, Ganeshkhind, Pune 411 007, INDIA Abstract
More informationBig Galaxies Are Rare! Cepheid Distance Measurement. Clusters of Galaxies. The Nature of Galaxies
Big Galaxies Are Rare! Potato Chip Rule: More small things than large things Big, bright spirals are easy to see, but least common Dwarf ellipticals & irregulars are most common Faint, hard to see Mostly
More informationPart two of a year-long introduction to astrophysics:
ASTR 3830 Astrophysics 2 - Galactic and Extragalactic Phil Armitage office: JILA tower A909 email: pja@jilau1.colorado.edu Spitzer Space telescope image of M81 Part two of a year-long introduction to astrophysics:
More informationLecture 22: The expanding Universe. Astronomy 111 Wednesday November 15, 2017
Lecture 22: The expanding Universe Astronomy 111 Wednesday November 15, 2017 Reminders Online homework #10 due Monday at 3pm Then one week off from homeworks Homework #11 is the last one The nature of
More informationLecture 14: Cosmological Principles
Lecture 14: Cosmological Principles The basic Cosmological Principles The geometry of the Universe The scale factor R and curvature constant k Comoving coordinates Einstein s initial solutions 3/28/11
More informationShedding Light on Dark Matter
Shedding Light on Dark Matter can t actually be done while these particles make up 82% of the matter in the Universe, they don t interact with light. So how do we know anything about these invisible particles?
More informationExam 3 Astronomy 100, Section 3. Some Equations You Might Need
Exam 3 Astronomy 100, Section 3 Some Equations You Might Need modified Kepler s law: M = [a(au)]3 [p(yr)] (a is radius of the orbit, p is the rotation period. You 2 should also remember that the period
More informationLicia Verde. ICREA & ICC-UB-IEEC CERN Theory Division.
Licia Verde ICREA & ICC-UB-IEEC CERN Theory Division http://icc.ub.edu/~liciaverde AIMS and GOALS Observational cosmology has been evolving very rapidly over the past few years Theoretical cosmology is
More informationMOND and the Galaxies
MOND and the Galaxies Françoise Combes Observatoire de Paris With Olivier Tiret Angus, Famaey, Gentile, Wu, Zhao Wednesday 1st July 2009 MOND =MOdified Newtonian Dynamics Modification at weak acceleration
More informationhalo formation in peaks halo bias if halos are formed without regard to the underlying density, then δn h n h halo bias in simulations
Physics 463, Spring 07 Bias, the Halo Model & Halo Occupation Statistics Lecture 8 Halo Bias the matter distribution is highly clustered. halos form at the peaks of this distribution 3 2 1 0 halo formation
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 Tools of Cosmology. Andrew Zentner The University of Pittsburgh
The Tools of Cosmology Andrew Zentner The University of Pittsburgh 1 Part Two: The Contemporary Universe 2 Contents Review of Part One The Pillars of Modern Cosmology Primordial Synthesis of Light Nuclei
More informationThe Cosmic Microwave Background
The Cosmic Microwave Background The Cosmic Microwave Background Key Concepts 1) The universe is filled with a Cosmic Microwave Background (CMB). 2) The microwave radiation that fills the universe is nearly
More informationAstronomy 422. Lecture 15: Expansion and Large Scale Structure of the Universe
Astronomy 422 Lecture 15: Expansion and Large Scale Structure of the Universe Key concepts: Hubble Flow Clusters and Large scale structure Gravitational Lensing Sunyaev-Zeldovich Effect Expansion and age
More informationIntroduction to Cosmology (in 5 lectures) Licia Verde
Introduction to Cosmology (in 5 lectures) Licia Verde http://icc.ub.edu/~liciaverde Program: Cosmology Introduction, Hubble law, Freedman- Robertson Walker metric Dark matter and large-scale cosmological
More informationASTR 610 Theory of Galaxy Formation Lecture 4: Newtonian Perturbation Theory I. Linearized Fluid Equations
ASTR 610 Theory of Galaxy Formation Lecture 4: Newtonian Perturbation Theory I. Linearized Fluid Equations Frank van den Bosch Yale University, spring 2017 Structure Formation: The Linear Regime Thus far
More informationAstro Assignment 1 on course web page (due 15 Feb) Instructors: Jim Cordes & Shami Chatterjee
Astro 2299 The Search for Life in the Universe Lecture 4 This time: Redshifts and the Hubble Law Hubble law and the expanding universe The cosmic microwave background (CMB) The elements and the periodic
More informationMichela Mapelli. LECTURES on COLLISIONAL DYNAMICS: 1. RELEVANT TIMESCALES, FORMATION OF STAR CLUSTERS, EQUILIBRIUM MODELS
Michela Mapelli LECTURES on COLLISIONAL DYNAMICS: 1. RELEVANT TIMESCALES, FORMATION OF STAR CLUSTERS, EQUILIBRIUM MODELS COLLISIONAL/COLLISIONLESS? Collisional systems are systems where interactions between
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 informationLarge-scale structure as a probe of dark energy. David Parkinson University of Sussex, UK
Large-scale structure as a probe of dark energy David Parkinson University of Sussex, UK Question Who was the greatest actor to portray James Bond in the 007 movies? a) Sean Connery b) George Lasenby c)
More informationThe Formation and Evolution of Galaxy Clusters
IAU Joint Discussion # 10 Sydney, July, 2003 The Formation and Evolution of Galaxy Clusters Simon D.M. White Max Planck Institute for Astrophysics The WMAP of the whole CMB sky Bennett et al 2003 > 105
More informationChapter 25 (and end of 24): Lecture Notes
Chapter 25 (and end of 24): Lecture Notes In order to understand the Hubble Law and its implications both for estimating distances and for interpreting the evolution of the universe, you have to be comfortable
More informationThe Next 2-3 Weeks. Important to read through Chapter 17 (Relativity) before I start lecturing on it.
The Next 2-3 Weeks [27.1] The Extragalactic Distance Scale. [27.2] The Expansion of the Universe. [29.1] Newtonian Cosmology [29.2] The Cosmic Microwave Background [17] General Relativity & Black Holes
More informationWhere are the missing baryons?
Where are the missing baryons? Qiang Wang (Brandon) B.Sc. Supervisor: Dr Scott Kay School of Physics and Astronomy, University of Manchester 27/29 June 2018 Image from Millennium Simulation 1 Content I.
More informationThe Los Cabos Lectures
January 2009 The Los Cabos Lectures Dark Matter Halos: 3 Simon White Max Planck Institute for Astrophysics Shapes of halo equidensity surfaces Group Jing & Suto 2002 Galaxy δ 100 2500 6250 Shapes of halo
More informationAstronomy 102: Stars and Galaxies Examination 3 Review Problems
Astronomy 102: Stars and Galaxies Examination 3 Review Problems Multiple Choice Questions: The first eight questions are multiple choice. Except where explicitly noted, only one answer is correct for each
More informationAnalyzing the CMB Brightness Fluctuations. Position of first peak measures curvature universe is flat
Analyzing the CMB Brightness Fluctuations (predicted) 1 st rarefaction Power = Average ( / ) 2 of clouds of given size scale 1 st compression 2 nd compression (deg) Fourier analyze WMAP image: Measures
More informationWriting a PM code. Andrey Kravtsov astro.uchicago.edu/ andrey/talks/pm/pmpotsdam.pdf. July 23, 2006
Writing a PM code Andrey Kravtsov andrey@oddjob.uchicago.edu astro.uchicago.edu/ andrey/talks/pm/pmpotsdam.pdf July 23, 2006 Cosmological simulations 1 Cosmological simulations: the problem The Universe
More informationDetailed energy (mass) budget of the Universe
Detailed energy (mass) budget of the Universe most baryons are in ionized gas outside galaxies Intergalactic baryons Quasar absorption lines zqso=3.18 The Lyα forest as a probe of intergalactic gas The
More informationThe Local Group Timing Argument
The Local Group Timing Argument Real galaxies Milky Way (MW) Andromeda (M31) Indranil Banik (ib45@st-andrews.ac.uk) Supervisor: Hongsheng Zhao University of Saint Andrews MNRAS, 467 (2), 2180 2198 Basic
More informationThe Milky Way in the cosmological context. Andrey Kravtsov The University of Chicago
The Milky Way in the cosmological context Andrey Kravtsov The University of Chicago Milky Way and Its Stars, KITP, 2 February 2015 Cosmological context: hierarchical structure formation from a Gaussian
More informationSPH simulations of galaxy formation
SPH simulations of galaxy formation Looking under the hood of GADGET-2 Volker Springel Large N-body runs Algorithmic aspects of a new simulation code: GADGET- II Non-standard physics in simulations of
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 informationModified Gravity (MOG) and Dark Matter: Can Dark Matter be Detected in the Present Universe?
Modified Gravity (MOG) and Dark Matter: Can Dark Matter be Detected in the Present Universe? John Moffat Perimeter Institute, Waterloo, Ontario, Canada Talk given at the Miami 2014 topical conference on
More informationGalaxy Hydrodynamic Simulations and Sunrise Visualizations
Galaxy Hydrodynamic Simulations and Sunrise Visualizations Joel Primack, UCSC Daniel Ceverino, HU Madrid Avishai Dekel, HU & UCSC Sandra Faber, UCSC Anatoly Klypin, NMSU Patrik Jonsson, Harvard CfA Chris
More informationarxiv: v1 [astro-ph.co] 24 Jun 2009
Cosmological simulations of galaxy clusters Stefano Borgani 1,2,3 & Andrey Kravtsov 4,5,6 arxiv:0906.4370v1 [astro-ph.co] 24 Jun 2009 1- Dipartimento di Astronomia dell Università di Trieste, via Tiepolo
More informationThe fine-scale structure of dark matter halos
COSMO11, Porto, August 2011 The fine-scale structure of dark matter halos Simon White Max-Planck-Institute for Astrophysics COSMO11, Porto, August 2011 Mark Vogelsberger The fine-scale structure of dark
More informationCo-Evolution of Central Black Holes and Nuclear Star Clusters
Co-Evolution of Central Black Holes and Nuclear Star Clusters Oleg Gnedin (University of Michigan) Globular clusters in the Galaxy median distance from the center is 5 kpc Resolved star cluster highest
More information