![Page 1: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/1.jpg)
Progress of code development
2nd-order Einstein-Boltzmann solver for CMB anisotropy
Yukawa Institute for Theoretical Physics (YITP)Kyoto University
Takashi Hiramatsu
Collaboration with Ryo Saito (APC), Atsushi Naruko (TITech), Misao Sasaki (YITP)
Relativistic Cosmology, 8-19 Sep 2014 @ YITP
![Page 2: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/2.jpg)
2/25
Takashi HiramatsuCurrent status
1st-order perturbations 2nd-order perturbations
Implementing basic equations
History of electron density
line-of-sight integral+angular power spectrum
Qualitative check
Quantitative check using CAMB
Implementing basic equations 1
2nd-order line-of-sight integral
Observables
1st-order
Speed-up + Optimisation
2nd-order
Implementing basic equations 2
NOW
95%
15%
NOWR.Saito, Naruko, Hiramatsu, Sasaki, arXiv:1409.2464
![Page 3: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/3.jpg)
3/25
Takashi HiramatsuEvolution of fluctuations
INFLATION RADIATION
MATTERDARK-ENERGY
curvature perturbation
Photon temperature/polarisation fluc.
Neutrino temperature fluc.
CDM/baryon density fluc.
CDM/baryon velocity fluc.
Gravity field fluc.
spatial fluctuationof gravity potential
last scattering surface
Boltzmann equation
Euler/continuity equation
Einstein equation
reheating
initial fluctuations
CosmoLib : Huang, JCAP 1206 (2012) 012
CMBFAST : Seljak, Zaldarriaga, APJ469 (1996) 437CAMB : Lewis, Challinor, APJ538 (2000) 473CLASS II : Blas, Lesgourgues, Tram, JCAP 1107 (2011) 034
CAMB
CMBFASTCLASSIICosmoLib
![Page 4: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/4.jpg)
4/25
Takashi HiramatsuPhotons/Neutrinos
They are described by the Bose-Einstein/Fermi-Dirac distribution functions:
Directional cosine Fourier trans.
Temperature fluctuation has three independent variables
Multipole expansion
Multipole exp.
![Page 5: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/5.jpg)
5/25
Takashi HiramatsuPhotons/Neutrinos
Distribution function satisfies the Boltzmann equation : Liouville term Collision term
The efficiency of collision term is controlled by , time-derivative of optical depth
: Thomson cross-section
Deep in radiation dominant epoch
Last-scattering surface
Late-time Universe
Extremely large tight-coupling
suddenly decays
No collision, free-streaming
S. Dodelson, “Modern Cosmology”
Reionisation is revived, but not so significant
![Page 6: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/6.jpg)
6/25
Takashi HiramatsuCDM/baryons
CDMbaryon
EB equations for CDM/baryons + photon coupling.
Relativistic Fluid equationsfor CDM and baryons
with source terms.Fluid approx.
Fouriertransform
![Page 7: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/7.jpg)
7/25
Takashi HiramatsuGravity
Conformal Newton gauge :
NOTE : CAMB, CMBFAST use sychronous gauge :
(non-dynamical)
![Page 8: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/8.jpg)
8/25
Takashi Hiramatsu1st-order perturbation equations
CDM, baryon
Photon temperature
Gravity
Photon polarisation Massless neutrino temperature
![Page 9: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/9.jpg)
9/25
Takashi HiramatsuLine-of-sight integral
Truncating the Boltzmann hierarchy at to ensureaccurate calculation of
Large Boltzmann hierarchy, say , is required, but it is too hard to calculate...
For larger , it is useful to use the integral representation, line-of-sight integral.
visibility function
cf.
Seljak, Zaldarriaga, APJ 469 (1996) 437
optically thick
suppressing early-timecontributions
![Page 10: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/10.jpg)
10/25
Takashi HiramatsuLine-of-sight integral
Monopole
Dipole
Quadrupole
Integrated Sachs-Wolfe effect
Sourced by fluctuations at LSS ( including Sachs-Wolfe effect )
directly solving
LoS integral
![Page 11: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/11.jpg)
11/25
Takashi HiramatsuInitial conditions
Deep in radiation dominant epoch where all modes are larger than horizon scale.which crosses the horizon at
RD,Unchanged potential, Similarly fluctuated,Superhorizon,Tight-coupling,Negligible photon's quadrupole, Non-negligible neutrino's.
= almost flatduring inflation
after inflation
Promordial perturbations
= almost flat
![Page 12: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/12.jpg)
12/25
Takashi HiramatsuInitial conditions
All quantities can be written by the primordial curvature perturbation.Hence it is convenient to define the transfer functions such as
The governing equations are reinterpreted into those of the transfer functions.
Not used untilcalculating
@
![Page 13: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/13.jpg)
13/25
Takashi HiramatsuSupporting actors : Background geometry
Flat FLRW model
It can be easily extended to include non-flat case.
fiducial parameters
NOTE : we use as the time variable instead of in solving EB equations. Then we don't have to solve the Friedmann equation.
![Page 14: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/14.jpg)
14/25
Takashi HiramatsuSupporting actors : Recombination/Ionisation
Number density of free electrons- Strength of Photon-Baryon coupling
- Opacity of the Universe (affects ISW)
At recombination : instanteneous / polynomial fitting of / Weinberg's textbook/ Recfast output
controls ...
CAMB uses more realistic model by Seager, but the difference is small.
cf.
He2+He+He+He
H+H Seager,Sasselov,Scott, APJS 128 (2000) 407
Peebles, APJ 153 (1968) 1
Weinberg :
![Page 15: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/15.jpg)
15/25
Takashi HiramatsuSupporting actors : Recombination/Ionisation
Number density of free electrons- Strength of Photon-Baryon coupling
- Opacity of the Universe (affects ISW)
controls ...
At reionisation : 2 stages / instanteneous / none
cf.
H+ HHe+ He
He2+ He+
2-stage reionisation
During each stage, the fraction is interpolated by tanh() function.
Same scheme as CAMB's.....(maybe)
![Page 16: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/16.jpg)
16/25
Takashi HiramatsuSupporting actors : Spherical Bessel functions
Use three different methods to maintain an accuracy of for
Descending reccurence :
Debye's expansion :
Taylor expansion :
is given by the partitioned polynomial interpolation.
[Descending]
[Debye]
[Tay
lor]
exponentially small
oscillatory behaviour
with arbitrary argumentOnce storing for ,
[Debye]
[Descending]
[Taylor]
![Page 17: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/17.jpg)
17/25
Takashi Hiramatsu
Angular power spectrum
Angular power spectrum (1st-order)
We have some existing codes named CMBFAST, CAMB, CosmoLib, CLASS, etc..
CAMB
PRELIMINARY
CosmoLib : Huang, JCAP 1206 (2012) 012
CMBFAST : Seljak, Zaldarriaga, APJ469 (1996) 437 CAMB : Lewis, Challinor, APJ538 (2000) 473
CLASS II : Blas, Lesgourgues, Tram, JCAP 1107 (2011) 034
still has 14% undesiredoverestimation...
![Page 18: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/18.jpg)
18/25
Takashi HiramatsuAngular power spectrum (1st-order)
2014/07/30
2014/06/20
2014/06/28
2014/08/18
2014/07/02
![Page 19: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/19.jpg)
19/25
Takashi HiramatsuGo to 2nd-order
Einstein-Boltzmann equations for 2nd-order quantities sourced by [1st-order]2
Line-of-sight integral sourced by [1st-order]2
CDM+Baryon+Gravity have been implemented,but Baryon-Photon/Gravity-Photon couplings are not considered yet.
Existing 2nd-order Boltzmann solver
2nd-order contributions appear in...
CMBquickSONG CosmoLib2nd
R.Saito, Naruko, Hiramatsu, Sasaki, arXiv:1409.2464
CosmoLib2nd : Huang, Vernizzi, arXiv:1212.3573
CMBquick : Creminelli, Pitrou, Vernizzi, arXiv:1109.1822
SONG : Pettinari, Filder, Crittenden, Koyama, Wands; Pettinari, arXiv:1405.2280 (thesis)
Formulations have been completed by R.Saito.
![Page 20: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/20.jpg)
20/25
Takashi HiramatsuExpansion up to 2nd-order
Poisson gauge + neglecting 1st-order and
Expanding up to 2and-order
![Page 21: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/21.jpg)
21/25
Takashi Hiramatsu2nd-order perturbation equtaions
CDM + gravity (contribution from radiation is work in progress....)
is a fractional expression like
'Q'uadratic terms of 1st-order quantities
Pending until deciding how to treat 2nd-orderradiation (temperature ? brightness ?)
and
![Page 22: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/22.jpg)
22/25
Takashi Hiramatsu2nd-order perturbation equtaions
CDM + gravity (contribution from radiation is work in progress....)
![Page 23: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/23.jpg)
23/25
Takashi Hiramatsu2nd-order line(curve)-of-sight integral
SW x ISWSW x LensingSW x Time-delaySW x Deflection
ISW x ISW
ISW x Time-delayISW x Lensing
We found 7 combinations in this formula,
R.Saito, Naruko, Hiramatsu, Sasaki, arXiv:1409.2464
![Page 24: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/24.jpg)
24/25
Takashi Hiramatsu2nd-order line(curve)-of-sight integral
SW x Lensing
Bispectrum
R.Saito, Naruko, Hiramatsu, Sasaki, arXiv:1409.0000
![Page 25: 2nd-order Einstein-Boltzmann solver for CMB anisotropy · 2nd-order Einstein-Boltzmann solver for CMB anisotropy Yukawa Institute for Theoretical Physics (YITP) Kyoto University Takashi](https://reader035.vdocument.in/reader035/viewer/2022081521/5e164b4fad791449cb22a543/html5/thumbnails/25.jpg)
25/25
Takashi HiramatsuOverview of my code
- Full scratch development, completely independent of existing codes
- Time evolution : 1-stage 2nd-order implicit Runge-Kutta (Gauss-Legendre) method (implementing up to 4th-order schemes)
- Line-of-sight Integration : Trapezoidal rule
- Parallelised by OpenMP
- C++
- Ready for implementing a variety of recombination/reionisation simulators
- We are now suffered from a small mismatch between results of our code and CAMB at the 1st-order.- Implemented 2nd-order perturbations only for gravity and matter.- Moved to evaluate the bi-spectrum using the 2nd-order LOS integral.
- Interpolation scheme : Polynomial approximation (up to )