The Millenium Simulation: Our Universe in a Supecomputer

Total Page:16

File Type:pdf, Size:1020Kb

The Millenium Simulation: Our Universe in a Supecomputer The Millenium Simulation: our Universe in a supecomputer Simon White Max Planck Institute for Astrophysics Moore's Law for Cosmological N-body Simulations Springel et al 2005 Millennium Run ● Computers double their speed every 18 months ● A naive N-body force 2 calculation needs N op's ● Simulations double their size every 16.5 months ● Progress has been roughly equally due to hardware and to improved algorithms Millennium Run Statistics Springel et al 2005: The Virgo Consortium 3 10 ● DM particle number: N = 2160 = 10,077,696,000 ≈ 10 5 ● Box size: L = 500 Mpc/h, Softening: = 5 kpc/h L/ = 10 ● Initial redshift: z = 127 init ● Cosmology: =1, =0.25, =0.045, h =0.73, n=1, =0.9 tot m b 8 ● 343,000 processor-hours on 512 nodes of an IBM Regatta (28 machine days @ 0.2 Tflops using 1 Tbyte RAM) ● Full raw and reduced data stored at 64 redshifts 27 Tbytes of stored data A testbed for simulating the formation of ~107 galaxies Goals for simulations of galaxy/AGN evolution ● Explore the physics of galaxy formation ● Understand the links between galaxy and SMBH formation ● Clarify why galaxy properties are related to clustering ● Determine how environment stimulates galaxy activity ● Interpret new multi-wavelength surveys of galaxies ● Check if such surveys can provide precision tests of and parameter estimates for the standard CDM paradigm Simulating galaxies /AGN in the Millennium Run Springel et al 2005; Croton et al 2005, De Lucia et al 2005 ● Build and store merger trees which encode the detailed assembly history of every z=0 halo and of the substructure within it ● Implement models for the formation/evolution of galaxies to follow -- accretion, shock-heating and cooling of diffuse gas into disks -- star formation from the ISM in disks -- stellar evolution [After Springel et al (2001) -- SN feedback and stellar winds and De Lucia et al (2004)] -- chemical enrichment/dust formation -- galaxy merging/morphological transformation ● Implement models for the growth of central black holes to follow -- formation and growth from ISM gas during mergers [ A f t e r K a u f f m an n -- black hole mergers following galaxy mergers & Haehnelt (2000)] ● Include “radio mode” feedback from BH at cooling flow centres -- energy feedback from BH in the central dominant galaxy depends on 1.5 BH mass, gas temperature and gas mass fraction f m T ∝ gas BH z = 0 Dark Matter Halo Mass Functions in the MS Springel et al 2005 Solid curves are the empirical fitting formula from Jenkins et al 2001 with no parameters adjusted At z = 0 half of all mass is in lumps of M > 2 × 1010M ⊙ Does halo clustering depend on formation history? Gao, Springel & White 2005 The 20% of halos with the lowest formation redshifts in a 30 Mpc/h thick slice M ~ 1011M halo ⊙ Does halo clustering depend on formation history? Gao, Springel & White 2005 The 20% of halos with the highest formation redshifts in a 30 Mpc/h thick slice M ~ 1011M halo ⊙ Does formation history depend on environment? Gao, Springel & White 2005 An equal number of randomly chosen DM particles Halo bias as a function of M = 1011M /h halo ⊙ mass and formation time Gao, Springel & White 2005 ● Bias increases smoothly with formation redshift ● The dependence on formation redshift is strongest at low mass ● This dependence is consistent neither with excursion set M = 6×1012M /h * ⊙ theory nor with HOD models z = 0 Dark Matter z = 0 Galaxy Light Springel, Frenk & White 2006 Galaxy autocorrelation function Springel et al 2005 mass For such a large simulation the purely statistical error bars are o][ 2dFGRS negligible on even for the galaxies Precise estimates of autocorrelation functions Luminous red galaxies in the SDSS Masjedi et al 2005 Precise estimates of autocorrelation functions Luminous red galaxies in the M.S. From public Millennium Simulation data archive M > 1011M * ⊙ g – r > 0.8 A test of the HOD model Galaxy autocorrelations are compared between two models (a) a full simulation of galaxy formation (b) the same simulation but with galaxy popula- tions shuffled among halos of fixed mass The bias plotted is the square root of the ratio of the correlation func's All galaxies A test of the HOD model Galaxy autocorrelations are compared between two models (a) a full simulation of galaxy formation (b) the same simulation but with galaxy popula- tions shuffled among halos of fixed mass The bias plotted is the square root of the ratio of the correlation func's Blue galaxies A test of the HOD model Galaxy autocorrelations are compared between two models (a) a full simulation of galaxy formation (b) the same simulation but with galaxy popula- tions shuffled among halos of fixed mass The bias plotted is the square root of the ratio of the correlation func's Red galaxies Baryon galaxies wiggles in mass the galaxy distribution linear Springel et al 2005 Power spectra from the Millennium run divided by a baryon- free CDM spectrum Galaxy samples are matched to plausible large observational surveys at given z The effects of “radio mode” feedback on z=0 galaxies Croton et al 2005 ● In the absence of a “cure” for the cooling flow problem, the most massive galaxies are: too bright too blue disk-dominated ● With cooling flows suppressed by “radio AGN” these galaxies are less massive red elliptical Effect of feedback on the Luminosity Function Full model with reionisation, AGN and SN feedback Croton et al 2005 Effect of feedback on the Luminosity Function Full model with reionisation, AGN and SN feedback Croton et al 2005 Effect of feedback on the Luminosity Function Full model with reionisation, AGN and SN feedback Croton et al 2005 Effect of feedback on the Luminosity Function Full model with reionisation, AGN and SN feedback Croton et al 2005 Effect of feedback on the Luminosity Function Full model with reionisation, AGN and SN feedback Croton et al 2005 Croton et al Other galaxy properties at z=0 2005 Evolution of feedback sources Croton et al 2005 B-band luminosity function evolution Mass function evolution Kitzbichler et al 2006 Formation histories of ellipticals De Lucia et al 2005 M > 4 109 * Distribution of star formation times M > 1011 * Distribution of assembly times Evolution of mass and galaxy correlations Springel, Frenk & White 2006 M – 5 log h < -20 I I < 24 “COSMOS” 1.4º x 1.4º Kitzbichler et al 2006 A bright quasar and its surroundings at 1 billion years One of the most massive dark matter clumps, containing one of the most massive galaxies and most massive black holes. M = 5 × 1012M M = 1011M h ⊙ * ⊙ SFR = 235 M /yr M = 108M ⊙ BH ⊙ The quasar's descendant and its surroundings today, at t = 13.7 billion years One of the most massive galaxy clusters. The quasar's descendant is part of the central massive galaxy of the cluster. M = 2 × 1015M h ⊙ http://www.g-vo.org Local Group “Timing Argument” mass estimates with a ΛCDM prior Li & White 2006 M =M (MW)+M (M31) Find all isolated pairs of disk- true 200 200 dominated galaxies with: ● 150 km/s < V < 300 km/s max,halo ● - 80 km/s > V > -180 km/s rad ● 500 kpc < D < 1 Mpc Find ratio of the true mass to the “Timing Argument” estimate for each pair Multiply by the TA estimate for the Local Group.
Recommended publications
  • Curriculum Vitæ 0000-0002-2769-9507
    Institut de Physique des 2 Infinis de Lyon 4 rue Enrico Fermi, Université Claude Bernard Lyon 1 Lyon, 69622 Villeurbanne B [email protected] Oliver Newton MusicalNeutron Musical-Neutron Curriculum Vitæ 0000-0002-2769-9507 Employment & voluntary work Vocational 2019– Postdoctoral Researcher, Institut de Physique des Deux Infinis, Lyon, France. 2014–2015 Graduate Trading Platform Engineer, Fidessa, Woking, UK. Volunteering 2014–2020 Founding trustee & Treasurer, UniBrass Foundation (charity number: 1159359), UK. Education 2015–2019 PhD, Institute for Computational Cosmology, Durham, UK. Probing the nature of dark matter with small-scale cosmology Supervisors: Prof. Adrian Jenkins and Prof. Carlos Frenk 2010–2014 BSc MPhys, University of Warwick, Coventry, UK, 1st Class (Hons). Masters project: Determining the fundamental properties of Higgs candidates at the LHC Awards and scholarships 2019 ICC Research Scholarship, Institute for Computational Cosmology, Durham, UK. 2015–2019 STFC Postgraduate Studentship, Institute for Computational Cosmology, Durham, UK. Conference contributions Contributed talks July 2020 EAS (Online), Leiden, Netherlands. Constraining the properties of WDM using the satellite galaxies of the Milky Way Jan 2020 VIRGO Consortium meeting, Durham, UK. Constraining the properties of WDM using the satellite galaxies of the Milky Way Sep 2019 CLUES Collaboration meeting, IN2P3, Lyon, France. Constraints on thermal relic WDM from satellites of the LG July 2019 Small Galaxies, Cosmic Questions, Durham, UK. Constraints on the mass of the thermal relic warm dark matter particle Jan 2019 DEX XV, Edinburgh, UK. Constraints on the mass of the thermal relic warm dark matter particle Dec 2018 VIRGO Consortium meeting, Leiden University, Netherlands. Constraints on the mass of the thermal relic warm dark matter particle Aug 2018 XXX IAU General Assembly, Vienna, Austria.
    [Show full text]
  • Gravitational Waves from Resolvable Massive Black Hole Binary Systems
    Mon. Not. R. Astron. Soc. 000, 000–000 (0000) Printed 22 October 2018 (MN LATEX style file v2.2) Gravitational waves from resolvable massive black hole binary systems and observations with Pulsar Timing Arrays A. Sesana1, A. Vecchio2 and M. Volonteri3 1Center for Gravitational Wave Physics, The Pennsylvania State University, University Park, PA 16802, USA 2School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK 3 Dept. of Astronomy, University of Michigan, Ann Arbor, MI 48109, USA Received — ABSTRACT Massive black holes are key components of the assembly and evolution of cosmic structures and a number of surveys are currently on-going or planned to probe the demographics of these objects and to gain insight into the relevant physical processes. Pulsar Timing Arrays (PTAs) currently provide the only means to observe gravita- > 7 tional radiation from massive black hole binary systems with masses ∼ 10 M⊙. The whole cosmic population produces a stochastic background that could be detectable with upcoming Pulsar Timing Arrays. Sources sufficiently close and/or massive gener- ate gravitational radiation that significantly exceeds the level of the background and could be individually resolved. We consider a wide range of massive black hole binary assembly scenarios, we investigate the distribution of the main physical parameters of the sources, such as masses and redshift, and explore the consequences for Pulsar Timing Arrays observations. Depending on the specific massive black hole population model, we estimate that on average at least one resolvable source produces timing residuals in the range ∼ 5 − 50 ns. Pulsar Timing Arrays, and in particular the future Square Kilometre Array (SKA), can plausibly detect these unique systems, although the events are likely to be rare.
    [Show full text]
  • The Illustristng Simulations: Public Data Release
    Nelson et al. Computational Astrophysics and Cosmology (2019)6:2 https://doi.org/10.1186/s40668-019-0028-x R E S E A R C H Open Access The IllustrisTNG simulations: public data release Dylan Nelson1*,VolkerSpringel1,6,7, Annalisa Pillepich2, Vicente Rodriguez-Gomez8,PaulTorrey9,5, Shy Genel4, Mark Vogelsberger5, Ruediger Pakmor1,6, Federico Marinacci10,5, Rainer Weinberger3, Luke Kelley11,MarkLovell12,13, Benedikt Diemer3 and Lars Hernquist3 Abstract We present the full public release of all data from the TNG100 and TNG300 simulations of the IllustrisTNG project. IllustrisTNG is a suite of large volume, cosmological, gravo-magnetohydrodynamical simulations run with the moving-mesh code AREPO. TNG includes a comprehensive model for galaxy formation physics, and each TNG simulation self-consistently solves for the coupled evolution of dark matter, cosmic gas, luminous stars, and supermassive black holes from early time to the present day, z = 0. Each of the flagship runs—TNG50, TNG100, and TNG300—are accompanied by halo/subhalo catalogs, merger trees, lower-resolution and dark-matter only counterparts, all available with 100 snapshots. We discuss scientific and numerical cautions and caveats relevant when using TNG. The data volume now directly accessible online is ∼750 TB, including 1200 full volume snapshots and ∼80,000 high time-resolution subbox snapshots. This will increase to ∼1.1 PB with the future release of TNG50. Data access and analysis examples are available in IDL, Python, and Matlab. We describe improvements and new functionality in the web-based API, including on-demand visualization and analysis of galaxies and halos, exploratory plotting of scaling relations and other relationships between galactic and halo properties, and a new JupyterLab interface.
    [Show full text]
  • Machine Learning and Cosmological Simulations I: Semi-Analytical Models
    MNRAS 000,1{18 (2015) Preprint 20 August 2018 Compiled using MNRAS LATEX style file v3.0 Machine Learning and Cosmological Simulations I: Semi-Analytical Models Harshil M. Kamdar1;2?, Matthew J. Turk2;4 and Robert J. Brunner2;3;4;5 1Department of Physics, University of Illinois, Urbana, IL 61801 USA 2Department of Astronomy, University of Illinois, Urbana, IL 61801 USA 3Department of Statistics, University of Illinois, Champaign, IL 61820 USA 4National Center for Supercomputing Applications, Urbana, IL 61801 USA 5Beckman Institute For Advanced Science and Technology, University of Illinois, Urbana, IL, 61801 USA Accepted 2015 October 1. Received 2015 September 30; in original form 2015 July 2 ABSTRACT We present a new exploratory framework to model galaxy formation and evolution in a hierarchical universe by using machine learning (ML). Our motivations are two- fold: (1) presenting a new, promising technique to study galaxy formation, and (2) quantitatively analyzing the extent of the influence of dark matter halo properties on galaxies in the backdrop of semi-analytical models (SAMs). We use the influential Millennium Simulation and the corresponding Munich SAM to train and test vari- ous sophisticated machine learning algorithms (k-Nearest Neighbors, decision trees, random forests and extremely randomized trees). By using only essential dark matter 12 halo physical properties for haloes of M > 10 M and a partial merger tree, our model predicts the hot gas mass, cold gas mass, bulge mass, total stellar mass, black hole mass and cooling radius at z = 0 for each central galaxy in a dark matter halo for the Millennium run. Our results provide a unique and powerful phenomenolog- ical framework to explore the galaxy-halo connection that is built upon SAMs and demonstrably place ML as a promising and a computationally efficient tool to study small-scale structure formation.
    [Show full text]
  • Indra: a Public Computationally Accessible Suite of Cosmological 푁-Body Simulations
    MNRAS 000,1–12 (2021) Preprint 5 August 2021 Compiled using MNRAS LATEX style file v3.0 Indra: a Public Computationally Accessible Suite of Cosmological #-body Simulations Bridget Falck,1¢ Jie Wang,2,3 Adrian Jenkins,4 Gerard Lemson,1 Dmitry Medvedev,1 Mark C. Neyrinck,5,6 Alex S. Szalay1 1Department of Physics and Astronomy, Johns Hopkins University, 3400 N Charles St, Baltimore, MD 21218, USA 2National Astronomical Observatories, Chinese Academy of Sciences, Beijing, 100012, China 3School of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing 100049, China 4Institute for Computational Cosmology, Department of Physics, Durham University, Durham DH1 3LE, UK 5Ikerbasque, the Basque Foundation for Science, 48009, Bilbao, Spain 6Department of Physics, University of the Basque Country UPV/EHU, 48080, Bilbao, Spain Accepted XXX. Received YYY; in original form ZZZ; Draft version 5 August 2021 ABSTRACT Indra is a suite of large-volume cosmological #-body simulations with the goal of providing excellent statistics of the large-scale features of the distribution of dark matter. Each of the 384 simulations is computed with the same cosmological parameters and different initial phases, with 10243 dark matter particles in a box of length 1 ℎ−1 Gpc, 64 snapshots of particle data and halo catalogs, and 505 time steps of the Fourier modes of the density field, amounting to almost a petabyte of data. All of the Indra data are immediately available for analysis via the SciServer science platform, which provides interactive and batch computing modes, personal data storage, and other hosted data sets such as the Millennium simulations and many astronomical surveys.
    [Show full text]
  • Virtu - the Virtual Universe
    VirtU - The Virtual Universe comprising the Theoretical Virtual Observatory & the Theory/Observations Interface A collaborative e-science proposal between Carlos Frenk University of Durham Ofer Lahav University College London Nicholas Walton Cambridge University Carlton Baugh University of Durham Greg Bryan University of Oxford Shaun Cole University of Durham Jean-Christophe Desplat Edinburgh Parallel Computing Centre Gavin Dalton University of Oxford George Efstathiou Cambridge University Nick Holliman University of Durham Adrian Jenkins University of Durham Andrew King University of Leicester Simon Morris University of Durham Jerry Ostriker Cambridge University John Peacock University of Edinburgh Frazer Pearce University of Nottingham Jim Pringle Cambridge University Chris Rudge University of Leicester Joe Silk University of Oxford Chris Simpson University of Durham Linda Smith University College London Iain Stewart University of Durham Tom Theuns University of Durham Peter Thomas University of Sussex Gavin Pringle Edinburgh Parallel Computing Centre Jie Xu University of Durham Sukyoung Yi University of Oxford Overseas collaborators: Joerg Colberg University of Pittsburgh Hugh Couchman McMaster University August Evrard University of Michigan Guinevere Kauffmann Max-Planck Institut fur Astrophysik Volker Springel Max-Planck Institut fur Astrophysik Rob Thacker McMaster University Simon White Max-Planck Institut fur Astrophysik 1 Executive Summary 1. We propose to construct a Virtual Universe (VirtU) consisting of the \Theoretical Virtual
    [Show full text]
  • The 2Df Galaxy Redshift Survey and Cosmological Simulations Carlos S
    University of Durham The 2dF galaxy redshift survey and cosmological simulations Carlos S. Frenk Institute for Computational Cosmology, Durham Institute for Computational Cosmology The 2dF Galaxy Redshift University of Durham Survey 1997- 2002 250 nights at 4m AAT 221,000 redshifts to bj<19.45 median z = 0.11 First 100k z's released June/01 Full catalogue released July/03 Institute for Computational Cosmology 2dF Galaxy Redshift Survey: University of Durham Team Members Ivan K. Baldry10 Carlton M. Baugh2 Joss Bland-Hawthorn1 Terry Bridges1 Russell Cannon1 Shaun Cole2 Matthew Colless3 Chris Collins13 Warrick Couch5 Nicholas Cross6 Gavin Dalton9 Kathryn Deely5 Roberto De Propris5 Simon P. Driver6 George Efstathiou8 Richard S. Ellis7 Carlos S. Frenk2 Karl Glazebrook10 Edward Hawkins12 Carole Jackson3 Ofer Lahav8 Ian Lewis9 Stuart Lumsden11 Steve Maddox12 Darren Madgwick8 Stephen Moody8 Peder Norberg2 John A. Peacock4 Will Precival4 Bruce A. Peterson3 Mark Seaborne9 Will Sutherland4 Keith Taylor7 Institutions 1Anglo-Australian Observatory 2University of Durham 3The Australian National University 4University of Edinburgh 5University of New South Wales 6University of St Andrews 7California Institute of Technology 8University of Cambridge 9University of Oxford 10Johns Hopkins University 33 people at 11University of Leeds 12University of Nottingham 13 Liverpool John Moores University 12 Institute for Computational Cosmology institutions The 2dF galaxy redshift survey University of Durham QuickTimeã and a YUV420 codec decompressor are needed to
    [Show full text]
  • Public Release of N-Body Simulation and Related Data by the Virgo Consortium
    Mon. Not. R. Astron. Soc. 000, 1–4 (2000) Printed 25 October 2018 (MN LATEX style file v1.4) Public Release of N-body simulation and related data by the Virgo consortium. C. S. Frenk,1 J. M. Colberg,2 H. M. P. Couchman,3 G. Efstathiou,4 A. E. Evrard,5 A. Jenkins,1 T. J. MacFarland,6 B. Moore,1 J. A. Peacock,7 F. R. Pearce,1 P. A. Thomas,8 S. D. M. White2 and N. Yoshida.2 (The Virgo Consortium)10 1 Dept Physics, South Road, Durham, DH1 3LE. 2 Max-Planck Inst. for Astrophysics, Garching, Munich, D-85740, Germany. 3 Department of Physics and Astronomy, McMaster University, Hamilton, Ontario, L8S 4M1, Canada 4 Institute of Astronomy, Madingley Road, Cambridge CB3 OHA 6 Now at 105 Lexington Avenue, Apt. 6F,New York, NY 10016 5 Dept Physics, University of Michigan, Ann Arbor, MI-48109-1120. 7 Royal Observatory, Institute of Astronomy, Edinburgh, EH9 3HJ 8 Astronomy Centre,CPES, University of Sussex, Falmer, Brighton, BN1 9QJ 10 http://star-www.dur.ac.uk/∼frazerp/virgo/ 25 October 2018 ABSTRACT We are making available on the WWW a selection of the archived data from N-body simulations carried out by the Virgo consortium and related groups. This currently includes: (i) time-slice, lightcone and cluster data from the two 109-particle Hubble volume simulations described by Evrard 1998; (ii) time-slice data from simula- tions of 4 different cold dark matter cosmological models with 2563 particles analysed by Jenkins et al 1998; (iii) Dark halo catalogs, merger trees and galaxy catalogs from arXiv:astro-ph/0007362v1 25 Jul 2000 the GIF project described by Kauffmann et al 1999.
    [Show full text]
  • Mocking the Universe: the Jubilee Cosmological Simulations
    Mocking the Universe: The Jubilee cosmological simulations Deep galaxy surveys, large-scale structure and dark energy. Spanish participation in future projects, Valencia, 29-30 March 2012 JUBILEE Juropa Hubble Volume Simulation Project Participants: P.I: S. Gottlober (AIP) I. I. Iliev (Sussex) II. G. Yepes (UAM) III. J.M. Diego (IFCA) IV. E. Martínez González (IFCA) Why do we need large particle simulations? Large Volume Galaxy Surveys DES, KIDS, BOSS, JPAS, LSST, BigBOSS, Euclid They will probe 10-100 Gpc^3 volumes • Need to resolve halos hosting the faintest galaxies of these surveys to produce realistic mock catalogues. Higher z surveys imply smaller galaxies and smaller halos->more mass resolution. • Fundamental tool to compare clustering properties Zero crossing of CF. 130Mpc of galaxies with theoretical predictions from Prada et al 2012 cosmological models at few % level. Not possible only with linear theory. Must do the full non-linear evolution for scales 100+ Mpc (BAO, zero crossing) Galaxy Biases: Large mass resolution is needed if BAO in BOSS, internal sub-structure of dm halos has to be properly 100Mpc scale resolved to map halos to galaxies. e.g. Using the Halo Abundance Matching technique (e.g. Trujillo-Gómez et al 2011, Nuza et al 2012) We will need trillion+ particle simulations Large Volume Galaxy Surveys A typical example: BOSS (z=0.1..0.7) Box size to host BOSS survey : 3.5 h^-1 Gpc BOSS completed down to galx with Vcir >350 km/s - > Mvir ~5x10^12 Msun. To properly resolve the peak of the Vrot in a dark matter halo we need a minimum of 500-1000 particles.
    [Show full text]
  • Numerical Cosmology: Recreating the Universe in a Supercomputer The
    March 2006 The Alexandria Lectures Numerical Cosmology: Recreating the Universe in a Supercomputer Simon White Max Planck Institute for Astrophysics The Three-fold Way to Astrophysical Truth OBSERVATION The Three-fold Way to Astrophysical Truth OBSERVATION THEORY The Three-fold Way to Astrophysical Truth OBSERVATION SIMULATION THEORY INGREDIENTS FOR SIMULATIONS ● The physical contents of the Universe Ordinary (baryonic) matter – protons, neutrons, electrons Radiation (photons, neutrinos...) Dark Matter Dark Energy ● The Laws of Physics General Relativity Electromagnetism Standard model of particle physics Thermodynamics ● Initial and boundary conditions Global cosmological context Creation of “initial” structure ● Astrophysical phenomenology (“subgrid” physics) Star formation and evoluti on The COBE satellite (1989 - 1993) ● Three instruments Far Infrared Absolute Spectroph. Differential Microwave Radiom. Diffuse InfraRed Background Exp Spectrum of the Cosmic Microwave Background Data from COBE/FIRAS a near-perfect black body! What do we learn from the COBE spectrum? ● The microwave background radiation looks like thermal radiation from a `Planckian black-body'. This determines its temperature T = 2.73K ● In the past the Universe was hot and almost without structure -- Void and without form -- At that time it was nearly in thermal equilibrium ● There has been no substantial heating of the Universe since a few months after the Big Bang itself. COBE's temperature map of the entire sky T = 2.728 K T = 0.1 K COBE's temperature map of the
    [Show full text]
  • The Millennium Simulation: Cosmic Evolution in a Supercomputer Simon White Max Planck Institute for Astrophysics the COBE Satellite (1989 - 1993)
    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 and in infrared radiation ● One instrument took a precise spectrum of the sky in microwaves COBE's temperature map of the entire sky T = 2.728 K T = 0.1 K COBE's temperature map of the entire sky T = 2.728 K T = 0.0034 K COBE's temperature map of the entire sky T = 2.728 K T = 0.00002 K Structure in the COBE map ● One side of the sky is `hot', the other is `cold' the Earth's motion through the Cosmos V = 600 km/s Milky Way ● Radiation from hot gas and dust in our own Milky Way ● Structure in the Microwave Background itself Structure in the Microwave Background Where is the structure? In the cosmic “clouds”, 40 billion light years away What are we seeing? Weak gravito-sound waves in the clouds When do we see these clouds? When the Universe was 400,000 years old, and was 1,000 times smaller and 1,000 times hotter than today How big are the structures? At least a billion light-years across (in the COBE maps) When were they made? During inflation, perhaps 10-30sec. after the Big Bang What did they turn into? Everything we see in the present Universe The WMAP Satellite at Lagrange-Point L2 The WMAP of the whole CMB sky Bennett et al 2003 What can we learn from these structures? The pattern of the structures is influenced by several things: --the Geometry of the Universe finite or infinite eternal or doomed to end --the Content of the Universe:
    [Show full text]
  • Optimising FLASH for Cosmological Simulations
    I V N E R U S E I T H Y T O H F G E R D I N B U Optimising FLASH for Cosmological Simulations Chris Daley August 24, 2007 M.Sc. in High Performance Computing The University of Edinburgh Year of Presentation: 2007 Abstract Runtime performance is a limiting factor in large cosmological simulations. A paral- lel application code called FLASH, which can be used for cosmological simulations is investigated in this report. FLASH is an Adaptive Mesh Refinement (AMR) code, parallelised using MPI, and incorporating a Particle-Mesh (PM) Poisson solver. Pro- files of a cosmological simulation based on FLASH, indicate that the procedure which accumulates particle mass onto the mesh is the major bottleneck. A parallel algorithm devised by the FLASH center may help to overcome this bottleneck. In this project, the algorithm is implemented for the latest beta release of FLASH, which cannot currently perform cosmological simulations. There is interest in FLASH because the adaptive mesh can be refined to resolve shocks in cosmological simulations. In addition, the new version of FLASH is designed with an emphasis on flexible and easy to extend code units. The communication strategy in the implemented parallel algorithm does not involve guard cell exchanges (halo swaps). As such, it uses less memory than the previous im- plementation in the previous version of FLASH. Due to time restrictions, the delivered implementation is only compatible with a uniform grid in both uniform and adaptive grid modes. In addition, results indicate it is approximately 3-5 times slower than the previous implementation.
    [Show full text]