Realistic Simulations of the Intergalactic Medium: The Search for - - PowerPoint PPT Presentation

realistic simulations of the intergalactic medium the
SMART_READER_LITE
LIVE PREVIEW

Realistic Simulations of the Intergalactic Medium: The Search for - - PowerPoint PPT Presentation

Realistic Simulations of the Intergalactic Medium: The Search for Missing Physics Michael Norman James Bordner San Diego Supercomputer Center University of California San Diego Collaborators: David Tytler, Pengfei Chen (UCSD) 6/6/2018 - M. L.


slide-1
SLIDE 1

Realistic Simulations of the Intergalactic Medium: The Search for Missing Physics

Michael Norman James Bordner San Diego Supercomputer Center University of California San Diego

Collaborators: David Tytler, Pengfei Chen (UCSD)

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 1

slide-2
SLIDE 2

Observing the intergalactic medium in quasar absorption line spectra

Source: M. Murphy

Lyman α forest

2 6/6/2018 - M. L. Norman Blue Waters Symposium 2018

slide-3
SLIDE 3

The Cosmic Web: Origin of the Ly α forest absorption

quasar Cosmic web formed by gravitational clustering of dark matter and baryons (H and He) Hundreds of millions of lightyears

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 3

slide-4
SLIDE 4

Simulation v. observation: good qualitative agreement

  • bserved

simulated Kirkman & Tytler (1997) Zhang et al. (1997)

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 4

slide-5
SLIDE 5

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 5

N=10243 L = 80 Mpc

Baryon Overdensity, z=3 Jena et al. (2005)

slide-6
SLIDE 6

A small but persistent discrepancy has emerged between sim. & obs. as higher precision has been achieved

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 6

Possible Explanation Finding Observational systematics NO Simulation not converged NO Simulation box size too small NO Missing IGM physics MAYBE Missing galaxy HI absorption MAYBE

This PRAC project

slide-7
SLIDE 7

Quantifying the discrepancy: Flux PDF

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 7

high transmission low transmission Too many Not enough

slide-8
SLIDE 8

An inverted T-ρ relation in the IGM improves agreement!

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 8

Bolton et al. (2008)

Inverted: lower density gas is hotter

slide-9
SLIDE 9

How to obtain an inverted T-ρ relation?

Treat quasars as time-dependent point sources rather than as a homogeneous UV background

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 9

Standard model

IGM ionized by homogeneous UVB

Quasar model

IGM ionized by quasar point sources UVB

slide-10
SLIDE 10
  • Enzo+MGFLD RT (Dan Reynolds)
  • QSOs randomly placed in halos

according to QLF every 45 Myr

  • 5 energy groups 54.4 – 400 eV

z=4 z=3.5 z=3 z=2

The Enzo Simulations

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 10

  • Range of spectral slopes F=(ν/ν0)-β within obs.

range

  • Reionization by galaxies handled with UVB
  • 10243 Cartesian grids; 80 & 126 Mpc boxes

z=4 z=3.5 z=3 z=2 temperature f(He III)

slide-11
SLIDE 11

Implicit Multi-Group Flux Limited Diffusion

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 11

Global linear system solves performed with the HYPRE library (LLNL)

slide-12
SLIDE 12

Why Multifrequency RT?

A: QSOs have hard UV spectrum

[54.4,65] eV photoheating [65,75] eV [75,125] eV [125,155] eV [155,400] eV

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 12

slide-13
SLIDE 13

Synthetic HI absorption spectra @ z=4

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 13

slide-14
SLIDE 14

Synthetic HI absorption spectra @ z=2

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 14

slide-15
SLIDE 15

Does it improve agreement with data?

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 15

slide-16
SLIDE 16

The flux power spectrum is still discrepant

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 16

slide-17
SLIDE 17

Summary and Conclusions

  • We have carried out simulations of the Lyman alpha

forest modeling QSOs as discrete point sources

  • Inhomogeneous He II reionization and heating have a

surprisingly small effect on the H I observables, compared to standard model which uses a UVB

  • Subtle and not-so-subtle discrepancies remain

between models and observations, suggesting there is still some Missing Physics

– Resolved galaxy halos? – New dark matter physics?

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 17

New PRAC

slide-18
SLIDE 18

Rollup

  • Key Challenges: high fidelity simulation of the IGM including

inhomogeneous quasar ionization and heating

  • Why it Matters: may explain discrepancy between theory and
  • bservation
  • Why Blue Waters: capability and capacity to develop and run

very large, computationally and memory intensive simulations

  • Accomplishments: high-res run completed and analyzed
  • Blue Waters team contributions: assistance with topology-

aware scheduling

  • Broader Impact: drove the development of MGFLD capability

for the Enzo community code

  • Shared Data: none yet
  • Products: MGFLD capability in Enzo code, papers in prep.

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 18

slide-19
SLIDE 19

RESERVE SLIDES

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 19

slide-20
SLIDE 20
  • Quasar model requires 3D time-dependent

multifrequency radiative transfer

– Very computationally intensive

  • Scale separation requires very large grids

– Must resolve Ly α forest absorbers (25 kpc) in a box large enough to contain hundreds of quasars (100 Mpc) 40003 grids (target)

  • Model development requires high throughput

for experimentation at scale

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 20

Why Blue Waters is Needed

Progress requires combination of capacity and capability that only Blue Waters can provide

slide-21
SLIDE 21

186 Mpc Projected He+ fraction

Helium reionization by time-dependent quasars: Enzo MGFLD simulation on Blue Waters

He++ He+

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 21

slide-22
SLIDE 22

Redshift Evolution of χ(He++) and T

He++ volume fraction Temperature

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 22

slide-23
SLIDE 23

Inhomogeneous heating of the IGM by quasars 10243 Enzo-MGFLD on Blue Waters

118 Mpc z=3.41 Z=3.41 Baryon density Baryon temperature

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 23

slide-24
SLIDE 24

118 Mpc He+ fraction Photoheating rate Z=3.41

Inhomogeneous heating of the IGM by quasars 10243 Enzo-MGFLD on Blue Waters

6/6/2018 - M. L. Norman Blue Waters Symposium 2018 24