Mapping Hot Gas in the Universe using the Sunyaev-Zeldovich Effect

Transcription

1 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

2 Happy (belated; March 1) 75th birthday, Rashid!

3 Where is a galaxy cluster? Subaru image of RXJ (Medezinski et al. 2010) 3

4 Where is a galaxy cluster? Subaru image of RXJ (Medezinski et al. 2010) 4

5 Subaru Subaru image of RXJ (Medezinski et al. 2010) 5

6 Hubble Hubble image of RXJ (Bradac et al. 2008) 6

7 Chandra Chandra X-ray image of RXJ (Johnson et al. 2012) 7

8 1σ=17 μjy/beam =120 μkcmb ALMA! 5 resolution T. Kitayama Chandra X-ray image of RXJ (Johnson 2012) ALMA Band-3 Imageetofal. the Sunyaev-Zel dovich effect at 92 GHz (Kitayama et al. 2016) 8

9 A clear displacement between the X-ray and SZ images. What is going on? 9

10 10

11 Multi-wavelength Data IX = Z dl n2e (TX ) ISZ = T kb g me c2 Z dl ne Te Optical: X-ray: SZ [microwave]: galaxies hot gas (107 8 K) hot gas (107-8 K) velocity dispersion spectroscopic TX electron pressure gravitational lensing Intensity ~ ne2l Intensity ~ netel

12 A Story about RXJ Let me tell you a little story about this particular cluster, which highlights the unique power of the SZ data to study cluster astrophysics A massive cluster with Msun at z=0.45 The most X-ray luminous galaxy cluster found in the ROSAT All Sky Survey Very compact, cool core cluster 12

13 1997 ROSAT/HRI image [Schindler et al.] 5 resolution kev Looked pretty spherical Thought to be a typical, relaxed, cooling-flow cluster 13

14 2001 SZ w/ Nobeyama [Komatsu et al.] 12 resolution The highest angular resolution SZ mapping at that time (The record holder for a decade) Chandra X-ray image of RXJ (Johnson et al. 2012) A surprise!

15 2001 SZ w/ Nobeyama [Komatsu et al.] 12 resolution The highest angular resolution SZ mapping at that time (The record holder for a decade) Chandra X-ray image of RXJ (Johnson et al. 2012) A surprise!

16 2002 X-ray w/ Chandra [Allen et al.] kev An excess X-ray emission found at the location of the SZ excess A hot gas, missed by ROSAT due to the lack of sensitivity at high energies!

17 A lesson learned X-ray observations are band-limited They are not usually not sensitive to very hot gas with temperature >10(1+z) kev SZ observations are not band-limited They are in principle sensitive to arbitrarily high temperatures (more precisely, pressure) SZ data probe electron pressure: a good probe of shock-heated gas due to mergers RXJ was thought to be a relaxed cluster. Our Nobeyama data challenged it, and now it is accepted that this cluster is a merging system! 17

18 We have ALMA. Now what? What is a new science we can do with such high resolution, high sensitivity measurements? Finding shocks and hot clumps is fun, but can we do something new and more quantitative? One example: Pressure fluctuations 18

19 Let s subtract a smooth component X-ray SZ 19

20 Let s subtract a smooth component Ueda et al., in prep X-ray SZ 20

21 Let s subtract a smooth component Ueda et al., in prep X-ray SZ Gas density is stirred ( sloshed ), but no change in pressure! Not sound waves => Unique measurements of the effective equation of state of density fluctuations 21

22 Full-sky Thermal Pressure Map North Galactic Pole South Galactic Pole 22 Planck Collaboration

23 We can simulate this arxiv: [MNRAS, 463, 1797 (2016)] Klaus Dolag (MPA/LMU) Volume: (896 Mpc/h) 3 Cosmological hydro (P-GADGET3) with star formation and AGN feed back 2 x particles (mdm=7.5x10 8 Msun/h) 23

24 Dolag, EK, Sunyaev (2016) 24

25 The local universe simulation reproduces the observed structures pretty well 25

26 Dolag, EK, Sunyaev (2016) 1-point PDF fits!! 26

27 Dolag, EK, Sunyaev (2016) Power spectrum fits!! provided that we use: m = =

28 Simple Interpretation Cl [not l 2 Cl ] multipole Randomly-distributed point sources = Poisson spectrum = i(fluxi) 2 / 4π 28

29 Simple Interpretation Cl [not l 2 Cl ] multipole Extended sources = the power spectrum reflects intensity profiles 29

30 l(l+1)cl/2π [μk 2 ] >10 15 Msun >2x10 15 Msun >5x10 13 Msun >5x10 14 Msun Adding smaller clusters Multipole 30

31 Simple Formula C` = Z dz dv dz Z dm dn dm y`(m,z) 2 2d Fourier transform of pressure yl with small l just gives the total thermal pressure, MT ~ M 5/3 Heavily weighted by massive clusters The mass function, dn/dm, is sensitive to the amplitude of fluctuations, σ8 31

32 999 Degree-scale SZ power spectrum is less sensitive to astrophysics in cluster cores Komatsu & Kitayama (1999) 32

33 2014 confirmed by simulations with varying AGN feedback McCarthy et al. (2014) 33

34 999 It is very sensitive to the amplitude of fluctuations Komatsu & Kitayama (1999) Komatsu & Seljak (2002) 34

35 2014 tension? Planck13 parameters McCarthy et al. (2014) 35

36 2014 similar Planck13 to parameters planck15 McCarthy et al. (2014) 36

37 Dolag, EK, Sunyaev (2016) C` / 3 m 8 8 m = =0.829 vs m = =

38 Dolag, EK, Sunyaev (2016) C` / 3 m 8 8 m =0.315 ~20% too large 8 =0.829 vs m = =

39 Bolliet, Comis, EK, Macias-Perez (2017) Closer look at the B. Bolliet measurements The compton-y power spectrum of Planck contains various foreground sources What you saw as the data points were the raw data minus the bestfitting foreground components When fitting, the Planck team used Gaussian covariance ignoring the trispectrum term with trispecturm without How does this affect the results? 39

40 Bolliet, Comis, EK, Macias-Perez (2017) tsz power slightly lower without with trispecturm 40

41 Bolliet, Comis, EK, Macias-Perez (2017) Closer look at the parameter dependence Mass Bias Hubble σ8 Ωm w ns 41

42 Bolliet, Comis, EK, Macias-Perez (2017) model-independent 42 Closer look at the parameter dependence 2.6% measurement! Essentially cosmological

43 Bolliet, Comis, EK, Macias-Perez (2017) model-independent 43 Closer look at the parameter dependence 2.6% measurement! Essentially cosmological

44 <latexit sha1_base64="jg4q3clubbodwxz2vytq2xgjvfo=">aaacd3icbzc7sgnbfizn4y3g26qlzwaqlstuiqiwqoinjrdbmeaswuzkjbkye2hmrbcwfqqbx8xgqsxw1s63cxirnprdwmd/zuhm+b1ico2o82llzmbn5heyi7ml5zxvnxt941ahsejq4aemvc1jgqqioiicjdqibcz3jfs9/swwxr0dpuuy3oaggqbpuohocm7qwc17t4fctofetzjjx+eppf/g/ashfioqhjq9klxsvfnwrqj/wz1ankxubtkfjxbiyx8c5jjpxxedcjsjuyi4hdtxidvejpdzf+oga+adbiajg1k6y5w27ytkvadpyp05ktbf64hvmu6fyu9p14bmf7v6jj3tzikckeyi+hhrj5yuqzpmh7afao5yyibxjcxfke8xxtiadhmmbhf65l9qosycfdzro3yxnekjs7bintkjljkhrxjjyqrcolknj+szvfgp1pp1ar2nwzpwzgat/jl1/gugtpr4</latexit> <latexit sha1_base64="jg4q3clubbodwxz2vytq2xgjvfo=">aaacd3icbzc7sgnbfizn4y3g26qlzwaqlstuiqiwqoinjrdbmeaswuzkjbkye2hmrbcwfqqbx8xgqsxw1s63cxirnprdwmd/zuhm+b1ico2o82llzmbn5heyi7ml5zxvnxt941ahsejq4aemvc1jgqqioiicjdqibcz3jfs9/swwxr0dpuuy3oaggqbpuohocm7qwc17t4fctofetzjjx+eppf/g/ashfioqhjq9klxsvfnwrqj/wz1ankxubtkfjxbiyx8c5jjpxxedcjsjuyi4hdtxidvejpdzf+oga+adbiajg1k6y5w27ytkvadpyp05ktbf64hvmu6fyu9p14bmf7v6jj3tzikckeyi+hhrj5yuqzpmh7afao5yyibxjcxfke8xxtiadhmmbhf65l9qosycfdzro3yxnekjs7bintkjljkhrxjjyqrcolknj+szvfgp1pp1ar2nwzpwzgat/jl1/gugtpr4</latexit> <latexit sha1_base64="jg4q3clubbodwxz2vytq2xgjvfo=">aaacd3icbzc7sgnbfizn4y3g26qlzwaqlstuiqiwqoinjrdbmeaswuzkjbkye2hmrbcwfqqbx8xgqsxw1s63cxirnprdwmd/zuhm+b1ico2o82llzmbn5heyi7ml5zxvnxt941ahsejq4aemvc1jgqqioiicjdqibcz3jfs9/swwxr0dpuuy3oaggqbpuohocm7qwc17t4fctofetzjjx+eppf/g/ashfioqhjq9klxsvfnwrqj/wz1ankxubtkfjxbiyx8c5jjpxxedcjsjuyi4hdtxidvejpdzf+oga+adbiajg1k6y5w27ytkvadpyp05ktbf64hvmu6fyu9p14bmf7v6jj3tzikckeyi+hhrj5yuqzpmh7afao5yyibxjcxfke8xxtiadhmmbhf65l9qosycfdzro3yxnekjs7bintkjljkhrxjjyqrcolknj+szvfgp1pp1ar2nwzpwzgat/jl1/gugtpr4</latexit> <latexit sha1_base64="jg4q3clubbodwxz2vytq2xgjvfo=">aaacd3icbzc7sgnbfizn4y3g26qlzwaqlstuiqiwqoinjrdbmeaswuzkjbkye2hmrbcwfqqbx8xgqsxw1s63cxirnprdwmd/zuhm+b1ico2o82llzmbn5heyi7ml5zxvnxt941ahsejq4aemvc1jgqqioiicjdqibcz3jfs9/swwxr0dpuuy3oaggqbpuohocm7qwc17t4fctofetzjjx+eppf/g/ashfioqhjq9klxsvfnwrqj/wz1ankxubtkfjxbiyx8c5jjpxxedcjsjuyi4hdtxidvejpdzf+oga+adbiajg1k6y5w27ytkvadpyp05ktbf64hvmu6fyu9p14bmf7v6jj3tzikckeyi+hhrj5yuqzpmh7afao5yyibxjcxfke8xxtiadhmmbhf65l9qosycfdzro3yxnekjs7bintkjljkhrxjjyqrcolknj+szvfgp1pp1ar2nwzpwzgat/jl1/gugtpr4</latexit> C` = Planck Mass Bias Z dz dv dz Z dm dn dm y`(m,z) 2 The key ingredient of the power spectrum is a profile of thermal pressure of electrons M 500c = M 500c,true /B 44

45 Mass Bias in ΛCDM Constraining the ΛCDM parameters by the Planck (TT+lowP) chain, we find B = 1.71 ± 0.17 (68%CL; Bolliet et al.) or, 1-b = 1/B = 0.58 ± 0.06 Adding the CMB lensing, we find B = 1.59 ± 0.13 (68%CL; Makiya, Ando & EK, in prep) Cf: Simulation by Dolag, EK & Sunyaev: B ~ Manifestation that the new Compton-Y power spectrum is lower

46 Towards Tomography Cross-correlating the Compton-Y map with galaxies with known redshifts! 46

47 Huchra et al. (2012) 2MASS Redshift Survey ~40K galaxies with the median redshift of

48 Huchra et al. (2012) 2MASS Redshift Survey ~40K galaxies with the median redshift of

49 Ando, Benoit-Levy & EK (2018) 2MRS Auto Power Dominated by 1-halo term in most of the angular scales => Good for cross-correlation with Compton-Y 49

50 Ando, Benoit-Levy & EK (2018) 2MRS Auto Power 50

51 Makiya, Ando & EK, in prep Cross-power! R. Makiya 51

52 Makiya, Ando & EK, in prep Mass-bias Consistency [for Planck TT+lowP+lensing] We get consistent mass bias from Compton-Y and 2MRS cross. Neat. 52

53 Makiya, Ando & EK, in prep Mass Dependence 53

54 Makiya, Ando & EK, in prep Mass Dependence Cross is sensitive to less massive halos: We can use this to explore the mass bias as a function of mass! 54

55 <latexit sha1_base64="jg4q3clubbodwxz2vytq2xgjvfo=">aaacd3icbzc7sgnbfizn4y3g26qlzwaqlstuiqiwqoinjrdbmeaswuzkjbkye2hmrbcwfqqbx8xgqsxw1s63cxirnprdwmd/zuhm+b1ico2o82llzmbn5heyi7ml5zxvnxt941ahsejq4aemvc1jgqqioiicjdqibcz3jfs9/swwxr0dpuuy3oaggqbpuohocm7qwc17t4fctofetzjjx+eppf/g/ashfioqhjq9klxsvfnwrqj/wz1ankxubtkfjxbiyx8c5jjpxxedcjsjuyi4hdtxidvejpdzf+oga+adbiajg1k6y5w27ytkvadpyp05ktbf64hvmu6fyu9p14bmf7v6jj3tzikckeyi+hhrj5yuqzpmh7afao5yyibxjcxfke8xxtiadhmmbhf65l9qosycfdzro3yxnekjs7bintkjljkhrxjjyqrcolknj+szvfgp1pp1ar2nwzpwzgat/jl1/gugtpr4</latexit> <latexit sha1_base64="jg4q3clubbodwxz2vytq2xgjvfo=">aaacd3icbzc7sgnbfizn4y3g26qlzwaqlstuiqiwqoinjrdbmeaswuzkjbkye2hmrbcwfqqbx8xgqsxw1s63cxirnprdwmd/zuhm+b1ico2o82llzmbn5heyi7ml5zxvnxt941ahsejq4aemvc1jgqqioiicjdqibcz3jfs9/swwxr0dpuuy3oaggqbpuohocm7qwc17t4fctofetzjjx+eppf/g/ashfioqhjq9klxsvfnwrqj/wz1ankxubtkfjxbiyx8c5jjpxxedcjsjuyi4hdtxidvejpdzf+oga+adbiajg1k6y5w27ytkvadpyp05ktbf64hvmu6fyu9p14bmf7v6jj3tzikckeyi+hhrj5yuqzpmh7afao5yyibxjcxfke8xxtiadhmmbhf65l9qosycfdzro3yxnekjs7bintkjljkhrxjjyqrcolknj+szvfgp1pp1ar2nwzpwzgat/jl1/gugtpr4</latexit> <latexit sha1_base64="jg4q3clubbodwxz2vytq2xgjvfo=">aaacd3icbzc7sgnbfizn4y3g26qlzwaqlstuiqiwqoinjrdbmeaswuzkjbkye2hmrbcwfqqbx8xgqsxw1s63cxirnprdwmd/zuhm+b1ico2o82llzmbn5heyi7ml5zxvnxt941ahsejq4aemvc1jgqqioiicjdqibcz3jfs9/swwxr0dpuuy3oaggqbpuohocm7qwc17t4fctofetzjjx+eppf/g/ashfioqhjq9klxsvfnwrqj/wz1ankxubtkfjxbiyx8c5jjpxxedcjsjuyi4hdtxidvejpdzf+oga+adbiajg1k6y5w27ytkvadpyp05ktbf64hvmu6fyu9p14bmf7v6jj3tzikckeyi+hhrj5yuqzpmh7afao5yyibxjcxfke8xxtiadhmmbhf65l9qosycfdzro3yxnekjs7bintkjljkhrxjjyqrcolknj+szvfgp1pp1ar2nwzpwzgat/jl1/gugtpr4</latexit> <latexit sha1_base64="jg4q3clubbodwxz2vytq2xgjvfo=">aaacd3icbzc7sgnbfizn4y3g26qlzwaqlstuiqiwqoinjrdbmeaswuzkjbkye2hmrbcwfqqbx8xgqsxw1s63cxirnprdwmd/zuhm+b1ico2o82llzmbn5heyi7ml5zxvnxt941ahsejq4aemvc1jgqqioiicjdqibcz3jfs9/swwxr0dpuuy3oaggqbpuohocm7qwc17t4fctofetzjjx+eppf/g/ashfioqhjq9klxsvfnwrqj/wz1ankxubtkfjxbiyx8c5jjpxxedcjsjuyi4hdtxidvejpdzf+oga+adbiajg1k6y5w27ytkvadpyp05ktbf64hvmu6fyu9p14bmf7v6jj3tzikckeyi+hhrj5yuqzpmh7afao5yyibxjcxfke8xxtiadhmmbhf65l9qosycfdzro3yxnekjs7bintkjljkhrxjjyqrcolknj+szvfgp1pp1ar2nwzpwzgat/jl1/gugtpr4</latexit> C` = Planck Mass Bias Z dz dv dz Z dm dn dm y`(m,z) 2 The key ingredient of the power spectrum is a profile of thermal pressure of electrons αp M 500c = M 500c,true /B 55

56 Makiya, Ando & EK, in prep Mass Dependence Nailed 56

57 Makiya, Ando & EK, in prep Redshift Dependence 57

58 Makiya, Ando & EK, in prep Redshift Dependence High-ell data of Compton-Y auto is the key. But foreground contamination 58

59 Makiya, Ando & EK, in prep Z-dependence Poorly Constrained 59

60 Summary New results on the SZ effect, from small to large: T. Kitayama 1. The first SZ image by ALMA - opening up a new study of cluster astrophysics via pressure fluctuations 2. The SZ power spectrum at l<1000 has been determined finally! And we can simulate it K. Dolag B. Bolliet R. Makiya 3. Detailed look at mass bias from the SZ power spectrum and cross-correlation tomography B = 1.5 ± 0.1 (68%CL) for Planck TT+lowP+lensing. Expect B ~ 1.2 for hydrostatic mass bias in the simulation. Origin? 60

61 Compton Y Map of RXJ ALMA on-source integration times 5.6 hours with 7-m array 2.6 hours with 12-m array Thank you TAC!