Investigating Cluster Astrophysics and Cosmology with Cross-Correlation of Thermal Sunyaev-Zel dovich Effect and Weak Lensing 2017/7/14 13th Rencontres du Vietnam: Cosmology Ken Osato Dept. of Physics, the Univ. of Tokyo In collaboration with Daisuke Nagai (Yale), Samuel Flender (ANL), Masato Shirasaki (NAOJ), and Naoki Yoshida (UTokyo/Kavli IPMU) Based on ArXiv:1706.08972
Contents 1. Weak Lensing, thermal Sunyaev-Zel dovich Effect and cross-correlations 2. Construction of Model and Covariance 3. Results: Implications to Cluster Astrophysics and Cosmology 4. Summary
Weak Lensing: Cosmic Shear Cosmic Shear : Gravitational lensing by large-scale structure. unbiased tracer of dark matter. convergence galaxies position kernel density Convergence map from HSC (Oguri+, 2017)
Thermal Sunyaev-Zel dovich Effect CMB z ~ 1300 Hot gas (galaxy cluster) z~3-0 CMB photon Variation of T Compton-y Observer z=0 e- x T e +1 =y x x T ez 1 y= T kb me c2 h 4,x= kt ne Te dl Sunyaev and Zel dovich (1972, 1980)
tsz by Planck The all-sky thermal Sunyaev-Zel dovich effect has already been measured by Planck. The power spectrum of Compton-y can place a tight constraint on σ8. (c.f. Komatsu and Seljak, 2002). PDF power spectrum CMB+BAO+lensing multipole CMB+BAO σ8(ωm/0.28)3/8 Planck Col. (2015)
Cross-Correlation of tsz and WL Astrophysics and Cosmology with WL and tsz Power spectra of WL and tsz give us information of the large-scale structure in the Universe and are very useful to constrain cosmological models. However, the combination of them (cross-correlation) provides us with additional and independent information. Furthermore, tsz directly reflects gas distribution. We can obtain implications on cluster astrophysics as well. From the observational aspect Cross-correlation does not suffer from noise auto-correlation in the assumption that noises of different observables are uncorrelated. ha obs B obs i = habi, A obs = A + N A, B obs = B + N B
Measurements of Cross-Correlation Several groups have already reported the detection of the cross-correlation. CFHTLenS x Planck 2PCFs r eff p (at z=0.37) (Van Waerbeke+, 1 Mpc 10 Mpc 2014) 6 10-9 ξ yκ (r) 5 10-9 4 10-9 3 10-9 2 10-9 1 10-9 0 10 0-1 10-9 -2 10-9 4 10-9 3 10-9 Planck A 6σ significance Planck I RCSLenS x Planck (Hojjati+, 2016) 13σ significance ξ yκ (r) 2 10-9 1 10-9 0 10 0-1 10-9 0 10 20 30 40 50 60 r [arcmin] separation [arcmin] separation [arcmin]
Contents 1. Weak Lensing, thermal Sunyaev-Zel dovich Effect and cross-correlations 2. Construction of Model and Covariance 3. Results: Implications to Cluster Astrophysics and Cosmology 4. Summary
Models First, let us consider the model of the signal. Theoretical prediction of auto- and cross-power spectra Spectra can be decomposed into two terms based on halo model. C yapple ` C yapple(1h) ` = C yapple(2h) ` = = C yapple(1h) ` Z zdec 0 Z zdec 0 Z Mmax M min + C yapple(2h) dz d2 V dzd `, Z Mmax M min dz d2 V dzd P m(k = `/d A,z) dm dn dm b(m,z)apple`(m,z) well calibrated by N-body sim. Tinker et al. (2008, 2010) dm dn dm y`(m,z)apple`(m,z), Fourier transform of Compton-y of halo We need pressure profile of halo. Z Mmax M min dm dn dm b(m,z)y`(m,z).
Analytic Model of Gas Profiles We employ the analytic gas density/pressure profile model of individual halo, which is proposed by Shaw+ (2010) and improved by Flender+ (2017). The model contain six free parameters, each of which describes a physical process (e.g., SNe/AGN feedback, non-thermal pressure). Basic ideas s DM density follows NFW profile. DM = (r/r s )(1 + r/r s ) 2 Gas profile is determined from Euler eq. with polytope relation. dp tot = g (r) d (r) dr dr, P tot / g hydro. sim. suggests Γ~1.2 To determine the normalization, stellar and AGN feedback E p E g,f = E g,i + " DM E DM + " f M c 2 + Energy of gas dynamical friction between gas and DM work by gas expansion
Gas Profile Model Non-thermal pressure Turbulent motion also can support the self-gravity of the halo. This effect is parametrized as, P nth r (r) = (1 + z) P tot R 500 0.8 Note: only thermal pressure contributes as tsz. Free parameters are calibrated by gas density and gas fraction of X-ray clusters. We fix parameters other than alpha. gas density Flender+ (2017)
Covariance Estimation Let us move on covariance matrix estimation. Simulations In order to estimate covariance matrix, we employ N-body simulations. Box size: (1 Gpc/h) 3 # of particles: 2048 3 Cosmology: Planck 2015 The snapshots cover up to z=4.13. 0 1000 2000 3000 4000 5000 Radial comoving distance [Mpc/h] Mock observations WL: Multiple plane method (White & Hu, 2000) tsz: First, we find halos with Rockstar (Behroozi+, 2013). For each halo, we solve the analytic model, and assign thermal pressure to the member particles. Transverse comoving distance [Mpc/h] 400 200 0 200 400 Redshift 0.0 0.2 0.5 1.0 1.5 2.0 3.0 4.13
Covariance Estimation Let us move on covariance matrix estimation. WL tsz Simulations In order to estimate covariance matrix, we employ 10 degcosmological N-body simulations. Box size: (1 Gpc/h) 3 # of particles: 2048 3 Cosmology: Planck 2015 Mock observations WL: Multiple plane method (White & Hu, 2000) tsz: First, we find halos with Rockstar (Behroozi+, 2013). For each halo, we solve the analytic model, and assign thermal pressure to the member particles.
Covariance Estimation Let us move on covariance matrix estimation. Simulations In order to estimate covariance matrix, we employ 10 degcosmological N-body simulations. Box size: (1 Gpc/h) 3 # of particles: 2048 3 Cosmology: Planck 2015 x30 Mock observations WL: Multiple plane method (White & Hu, 2000) tsz: First, The we variance find halos of with maps Rockstar corresponds (Behroozi+, 2013). to For each halo, we solve the analytic model, and assign the covariance matrix thermal pressure to the member particles. x30
Contents 1. Weak Lensing, thermal Sunyaev-Zel dovich Effect and cross-correlations 2. Construction of Model and Covariance 3. Results: Implications to Cluster Astrophysics and Cosmology 4. Summary
Compton-y Spectra 10 12 `(` + 1)C yy (`)/2 10 0 10 1 10 2 simulation halo model (Planck col., 2013) halo model (Shaw et al., 2010) Planck col., 2015 ACT (Sievers et al., 2013) SPT (George et al., 2015) analytic model simulation Planck 2015 10 3 10 1 10 2 10 3 10 4 ` multipole KO+ (2017)
Cross-Correlation Function 3.5 3.0 2.5 simulation halo model (Planck col., 2013) halo model (Shaw et al., 2010) Hojjati et al. (2016) 10 9 yapple ( ) 2.0 1.5 1.0 0.5 0.0 RCSLenS x Planck 0 25 50 75 100 125 150 175 [arcmin] KO+ (2017)
Constraints on Non-thermal Pressure Non-thermal pressure Non-thermal contribution is hard to measure by conventional X-ray observations of clusters. However, the power spectrum and the cross-correlation are sensitive to it and can be used to estimate its contribution. Note: Non-thermal pressure and σ 8 are strongly degenerate. Power spectrum of Compton-y Cross-correlation function Small Large Varying the amplitude of non-thermal pressure
Constraints from Cross-Correlation We can constrain the amplitude of non-thermal pressure and σ 8 with power spectrum and cross-correlation. 8 1.0 0.9 0.8 0.7 0.6 0.5 P nth P tot (r) = (1 + z) Both 3 2 2 3 1 1 3 3 r R 500 0.8 Huge tension! Cross only 2 Power only 1 0.4 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0 KO+ (2017)
Removing Small Scales Correlations 1. Halo model assumption breaks down. 2. Incomplete separation from foreground contamination 8 1.0 0.9 0.8 0.7 C yy only C yy and yapple with small scale cuts yapple with small scale cuts 3 3 2 2 2 1 1 1 0.6 3 3 3 2 3 0.5 2 0.4 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0 KO+ (2017)
Source of Tension The analytic model is calibrated with galaxy clusters and low-z galaxy groups. However, high-z groups contribute to a substantial fraction of signal. At this range, analytic model inevitably contains uncertainty. Contribution from high-z groups (z > 0.2 and M 500c < 4 10 14 M sun /h) 10 12`(` + 1)C yy (`)/2 Ratio [%] 10 0 10 1 10 2 10 3 100 50 0 All z>0.2,m 500c < 4 10 14 M /h Around 50% of the signal comes 1.5 1.0 from unresolved halos! 0.5 10 1 10 2 10 3 10 4 ` 10 9 yk ( ) Ratio [%] 3.0 2.5 2.0 0.0 100 75 50 25 0 All z>0.2,m 500c < 4 10 14 M /h 0 20 40 60 80 100 [arcmin] KO+ (2017)
Possible Solution of Tension Let us consider an extreme case of enhanced star formation for group size halos. That corresponds to severe depletion of hot gas for such halos. 1.0 S M 500 M /M 500 = f 3 10 14 M 0.9 0.8 1 8 0.7 1 3 2 3 2 1 2 0.6 3 0.5 0.4 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0 KO+ (2017)
Contents 1. Weak Lensing, thermal Sunyaev-Zel dovich Effect and cross-correlations 2. Construction of Model and Covariance 3. Results: Implications to Cluster Astrophysics and Cosmology 4. Summary
Summary Weak lensing and thermal Sunyaev-Zel dovich effect are promising probes into the large-scale structure. Theoretical modeling and simulations can be used to estimate the signal and covariance matrix. The cross-correlation of them can provide us with additional information of cosmology and cluster astrophysics. We are currently working on the measurement of the cross-correlation with HSC and Planck data.