# The first detection of the imprint of filaments on CMB lensing

Galaxy redshift surveys, such as 2dF 2dF (), SDSS SDSS (), 6df 6df (), GAMA GAMA (), and VIPERS VIPERS (), have shown that the spatial distribution of matter forms a hierarchical structure consisting of clusters, filaments, sheets and voids. This hierarchical structure is known as the cosmic web bond_theo (). The majority of galaxy survey analyses measure the 2-point correlation, but ignoring the information beyond a small number of summary statistics. Since the matter density field becomes highly non-Gaussian as structures evolve, we expect other statistical descriptions of the field to provide us with additional information. One way to study the non-Gaussianity is to study filaments, which evolve non-linearly from the initial density fluctuation. Several previous works have studied the gravitational lensing of filaments to detect filaments and learn their mass profile fil_lensing2 (); fil_lensing_maturi (); fil_lensing_clampitt (). In our study, we provide the first detection of CMB (Cosmic Microwave Background) lensed by filaments and we measure how filaments trace the matter distribution on large scales. More specifically, we assume that, on large scales, filaments trace matter with a constant filament bias, defined as the ratio between the filament overdensity and the mass overdensity. We propose a phenomenological model for the cross power spectrum between filaments and the CMB lensing convergence field. By fitting the model to the data, we measure filament bias.

The cross-correlations of CMB lensing with tracers of large-scale structure have been widely studiedSmith:2007rg (); Hirata:2008cb (); Bleem:2012gm (); Sherwin:2012mr (); Ferraro:2014msa (); Allison:2015fac (); Allison:2015fac (); Giannantonio:2015ahz (); Eg (); Doux:2016xhg (); Singh:2016xey () to test General Relativity, probe the galaxy-halo
connection, etc. In our study, we detect the imprint of filaments on CMB lensing by cross-correlating filaments with the CMB lensing convergence () map. We use a filament intensity map, which is derived from the Cosmic Web Reconstruction filament catalogue fil_cat () ^{1}^{1}1Public in https://sites.google.com/site/yenchicr/catalogue from the Sloan Digital Sky Survey (SDSS) SDSS () Baryon Oscillations Spectroscopic Survey (BOSS) BOSS () Data Release 12 (DR 12) DR (). The filament finder (See Filament Finder in Method section) partitions the universe from 0.450 to 0.700 into slices with 0.005; filaments are found in each redshift bin as the density ridge of the smoothed galaxy density field yenchiSCMS ().
At redshift , the filament intensity is defined as

(1) |

where is the angular position, is the projected (angular) position of onto the nearest filament and is the uncertainty of filament at the projected position . The intensity function can be viewed as applying a Gaussian smoothing to the filament field; the smoothing scale here is chosen to be equal to the local uncertainty of filaments. Using the intensity map at each redshift bin, we construct the filament intensity overdensity map via

(2) |

Namely, the filament intensity overdensity map is obtained by stacking all the intensity maps together from redshift 0.45 to 0.7 and then scaling it by the average of the intensity .

In this work we use the CMB lensing convergence map ^{2}^{2}2missing from the Planck Planck_overview () satellite experiment. The Planck mission has reconstructed the lensing potential of the CMB from a foreground-cleaned map synthesized from the raw Planck 2015 full-mission frequency maps using the SMICA code Planck (). The lensing convergence is defined in terms of the lensing potential as

(3) |

We measure the cross angular power spectrum of CMB lensing convergence and filaments using standard techniques (See Estimator in Method section). We compute the error for each power spectrum by jackknife resampling 77 equally weighted regions that comprise the CMASS galaxy survey from where the filaments are detected.

We construct a phenomenological model to describe the cross-correlation of filaments and the CMB lensing convergence field. Instead of modeling the filament profile on small scalesfil_lensing_clampitt (); 2017MNRAS.468.2605E (); Higuchi:2014eia (), our model studies how filaments trace matter distribution on large scales through the use of the filament bias. We assume a CDM cosmology with Planck parameters from the 2015 release Planck_parameter (), where . In a spatially flat Friedmann-Robertson-Walker universe described by general relativity, the convergence field is

(4) |

where is the comoving radial distance, is the redshift observed at radial distance, is Newton’s gravitational constant, is the present-day mean density of the universe, and is the comoving distance to the CMB. On linear scale, we assume that filaments trace the matter as , where is defined as the large-scale filament bias.

The filament overdensity can be rewritten as

(5) |

where is the mean filament intensity redshift distribution defined as

(6) |

where is the filament intensity defined in eq. (1) and is the width of redshift slice, which is set to be 0.005. In cross-correlation, on scales smaller than the typical filament length, using filaments introduces additional smoothing compared to the true matter density. We model the smoothing as follows: the filaments have typical length and we lose all small-scale information about fluctuations along the filament; therefore, we take the corresponding filament power spectrum to be exponentially suppressed below the filament scale 1/(filament length) in Fourier space. Similarly, any matter in between filaments is either assigned to a filament or eliminated from the catalog (in underdense regions) . For this reason we also introduce a suppression in power in the direction perpendicular to the filaments, with suppression scale . We use two ways to model . The detailed models are shown later in the paper. Using the Limber approximation and the smoothing scale for small scales, the filament-convergence cross-correlation can be written as

(7) |

where is the CMB lensing kernel, and is modeled as

(8) |

where is the matter power spectrum. We use CAMB (Code for Anisotropies in the Microwave Background) ^{3}^{3}3http://www.camb.info/ to evaluate the theoretical matter power spectrum . The measurement of filament length is shown in Fig. 1. As shown in Fig. 1, we obtain very long filaments, and the mean (median) length of the filaments increases as a function of redshift due to the combination of two factors. Firstly, there is growth of the structure. While the number of galaxy clusters increases as the structure evolves, the length of filaments, acting as the mass bridges between galaxy clusters, will decrease. Secondly, the number of filaments detected also depends on the number density of galaxies, which, in the CMASS sample, is low and decreases as a function of redshift (See Supplement Fig.2 for galaxy redshift distribution and filament intensity redshift distribution). The large difference in the mean and median values of filament length indicates that the distribution of the filament length in each redshift bin is not Gaussian. We plot in the background the 2d histogram of filament length distribution as a function of redshift and filament length.

To check the validity of our model, we also compare the results to simulations. The excellent agreement that we find provides an important consistency check. The theoretical prediction for is shown in eq. (7). The matter-filament correlation is defined as

(9) |

By taking the parameters that are slowly varying compared to , we can get

(10) |

where . And for the filament catalogue, the effective redshift, which is defined as the weighted mean redshift of filament intensity, is 0.56. This approximation is not perfect, leading to a systematic bias in the prediction for . We propose an estimator for this systematic bias in Supplement Section 1 Eg (). As shown in Supplement Fig. 1, the systematic bias is less than 5%. Thus, the approximation only causes a negligible bias. In our analysis, we measure using 10 realizations of sky mocks of dark matter and corresponding filaments constructed from 10 N-body simulations (See sky mocks for dark matter and filament in Method section).

Fig. 2a shows the cross angular power spectrum of filaments and the CMB lensing convergence field. We bin our sample into 16 bins. The results from observation and simulation are shown as the blue crosses and green circles. Comparing simulation with data, we get with all 16 data points and without the first data point. The deviation of the first data point from the prediction is likely due to cosmic variance given the small sky area ( = 0.065) covered by the simulations. Thus, the first data can be very noisy. The result shows that the simulation is consistent with the data, indicating that the systematics in the data would not significantly affect our result.

We fit the model of eq. 7 to the data with filament bias () as the fitting parameter. We use two different smoothing methods to find . The first method is to define the perpendicular smoothing scale as the filament spacing, since any scale smaller than the filament spacing is smoothed out. The filament spacing is approximately the filament length. So, in this case, filament length is the overall smoothing scale for the effective power spectrum in eq. (8). The result is shown by the red line in Fig. 2. The best fit gives = 1.68 0.334. Since filaments also have width, filament spacing may be an overestimate of the smoothing perpendicular to filaments. In the second model, we fit for smoothing scale in the perpendicular direction as a free parameter, where we assume and we fit both for and constant using least fit. We get = 0.65 and correspondingly, . The result is shown as the orange line in Fig. 2.

We measure the significance of the cross-correlation detection by measuring the signal-to-noise ratio (SNR). Our SNR is defined as follows

(11) |

where

(12) | |||

(13) |

where is the cross angular signal in bin , is the best-fit theoretical prediction for the cross signal in bin , and is the covariance matrix estimated from jackknife resampling. The final result is shown in Table 1. The SNR values for both models show a significant detection of the cross-correlation. On large scales, we find that the filaments trace the matter with the filament bias around 1.5, which is smaller than galaxy bias, as expected, since the filaments also trace underdense regions where there are fewer or even no galaxies.

In order to validate the detection of the cross power spectrum between filaments and CMB lensing convergence field, we perform a null test as follows. We rotate the CMB lensing convergence map by , and , and then we cross correlate these rotated CMB convergence maps with the filament intensity map. Fig. 2b shows that the cross signal with the rotated maps fluctuates around 0. for the three cross angular power spectra is 0.79, 0.75 and 1.04. in all three cases is close to 1, which means the cross-correlation between rotated CMB maps and the filament intensity map is consistent with 0.

We define the cross-correlation coefficient between the filament and galaxy maps as , where is the cross angular power spectrum of filaments and galaxies, is the auto angular power spectrum of filaments and is the auto angular power spectrum of galaxies. The result is shown in the left panel of Fig. 3a. The signal is highly correlated on large scales, since both galaxies and filaments trace the large-scale structure of the matter. However, the correlation decreases on small scales. Fig. (3)b shows the cross-correlation of and , where is the angular cross power spectrum of the CMB lensing convergence map and the CMASS galaxy catalogue. These two figures show that the maps are not totally correlated with large deviations at small scales. Establishing the amount of extra cosmological information present in the filaments field would require a joint analysis with galaxy clustering and lensing measurements; this is left to future work.

## Method

Filament finder. We obtain filaments from the publicly available Cosmic Web Reconstruction filament catalogue fil_cat (). This filament catalogue is shown to have strong agreement with the redMaPPer Catalogue, since most clusters in the redMaPPer Catalogue lie at the intersection of the filaments in the Cosmic Web Reconstruction filament catalogue, which is predicted by theory. The filament catalogue also has good consistency with the Voronoi model yenchiSCMS (). Furthermore, filaments in this catalogue show desired effects on nearby galaxies gal_property_fil (), and filaments within this catalogue are featured with a local uncertainty measurement at the positions of filaments yenchiSCMS (); fil_cat ().

Here is a brief description of how filaments are obtained. The filament finder begins by partitioning the universe from 0.450 to 0.700 into slices with 0.005. Then for each slice, the galaxies are smoothed into a density field by the kernel density estimator (KDE)KDE ():

(14) |

where is the total number of galaxies in the slice, is the Gaussian function, is the angular position of the i-th galaxy at the slice at redshift , and is the smoothing parameter that is chosen depending on the density of galaxies in each redshift bin

(15) |

where is the minimal value for the standard deviation of galaxies’ location within this slice along each coordinate. The detailed explanation for the selection method is shown in APPENDIX A of yenchiSCMS ().
The KDE can be viewed as applying a Gaussian smoothing to galaxies to reconstruct the underlying probability density field.
We then remove galaxies in low-density regions (details can be found in yenchiSCMS ())
and apply the subspace constrained mean shift algorithm (SCMS) SCMS (); yenchiSCMS () to find filaments. The SCMS algorithm returns points representing the ridges of the probability density function . The filament’s uncertainty is obtained by the bootstrap approach.

Filament length measurement.
We get the filament intersections from Chen et al^{4}^{4}4Public in https://sites.google.com/site/yenchicr/catalogue fil_cat (). For each redshift bin, we use the hierarchical clustering method yenchiSCMS () to determine the number of branches at each intersection. The parameters in the hierarchical clustering are chosen to be the same as yenchiSCMS ():

(16) |

where is the smoothing
bandwidth as in eq. 14. At each intersection, we find the nearest point to the intersection point from each branch, and we group the nearest point as the filament point belonging to that filament (See Supplement Fig.3 for an example of filaments grouped.). Then we keep finding the nearest point to the newly grouped filament to find the next filament point belonging to that branch.
We stop the loop if the distance between filament points is less than and the distance between a filament point and the other intersection point is larger than .

Estimator.
We construct the filament map using the HEALPix pixelization with Nside=512. The CMB lensing convergence map is given directly by PLANCK using the HEALPix pixelization with Nside= 2048. We downgrade the lensing convergence map resolution to Nside=512 to cross-correlate with the filament map. The choice of resolution considers both the noise of the filament map and the smallest scale induced by the resolution. The resolution Nside = 512 gives us a pretty clean filament map; the scale of the pixel at the relevant redshift is approximately 3.5 Mpc, which is much smaller than the other length scales considered in the paper.

We measure the cross angular power spectrum for the filament catalogue and the CMB lensing convergence field using a pseudo- estimator:

(17) |

where is the sky fraction observed by both the filament catalogue and the CMB lensing convergence field, is the spherical harmonic transform of the lensing convergence field, and is the spherical harmonic transform of the filament intensity overdensity. The spherical harmonic transform and are computed using HEALPY.

Sky mock for filaments and dark matter.
We use -body simulation runs using the TreePM method (Bagla2002, ; White2002, ; Reid14, ). We use 10 realizations of this simulation based on the CDM model with and .
These simulations are in a periodic box of side length 1380Mpc and particles. A friend-of-friend halo catalogue is constructed at an effective redshift of . This is appropriate for our measurement since the galaxy sample used has effective redshift of 0.57. We use a Halo Occupation Distribution (HOD) (Peacock2000, ; Seljak2000, ; Benson2000, ; White2001, ; Berlind2002, ; Cooray2002, ) to relate the observed clustering of galaxies with halos measured in the -body simulation. We have used the HOD model proposed in Beutler13 () to populate the halo catalogue with galaxies.

(18) |

where is the average number of central galaxies for a given halo mass and is the average number of satellite galaxies. We use the HOD parameter set () from Beutler13 (). We have populated central galaxies at the center of our halo. The satellite galaxies are populated with radius (distance from central galaxy) distributed out to as per the NFW profile; the direction is chosen randomly with a uniform distribution.

The sky mocks of dark matter and galaxy are obtained from the simulation box using the method described in 2014MNRAS.437.2594W (). We use publicly available “MAKE SURVEY” ^{5}^{5}5https://github.com/mockFactory/make_survey code to transform a periodic box into the pattern of survey. The first step of this transformation involves a volume remapping of the periodic box to sky coordinates preserving the structure in the simulation. This is achieved by using the publicly available package called “BoxRemap”^{6}^{6}6http://mwhite.berkeley.edu/BoxRemap 2010ApJS..190..311C (). The BoxRemap defines an efficient volume-preserving, structure-preserving and one-to-one map to transform a periodic cubic box to non-cubical geometry. The non-cubical box is then translated and rotated to cover certain parts of the sky. We then convert the cartesian coordinate to the observed coordinate, which is right ascension, declination and redshift. We down-sample the galaxies with redshift to match the mock redshift with the redshift distribution observed in data. We request the reader refer to 2014MNRAS.437.2594W () for more details. We then apply the filament detection algorithm to these simulated mocks using the method described in Filament Finder.

model 1 | model 2 | |

1.680.334 | 1.470.28 | |

5.0 | 5.2 | |

25.77 | 24.39 | |

15 | 14 |

## References

- (1) Colless, M. et al. The 2dF Galaxy Redshift Survey: Spectra and redshifts. Mon. Not. Roy. Astron. Soc., 328:1039, 2001.
- (2) Eisenstein, D.J. et al. SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems. Astron. J., 142:72, 2011.
- (3) Jones, H.D. et al. The 6dF Galaxy Survey: Final Redshift Release (DR3) and Southern Large-Scale Structures. Mon. Not. Roy. Astron. Soc., 399:683, 2009.
- (4) Liske, J. et al. Galaxy and mass assembly (gama): end of survey report and data release 2. Monthly Notices of the Royal Astronomical Society, 452(2):2087, 2015.
- (5) Scodeggio, M. et al. The VIMOS Public Extragalactic Redshift Survey (VIPERS). Full spectroscopic data and auxiliary information release (PDR-2). ArXiv e-prints, November 2016.
- (6) Bond, R.J., Kofman, L., and Pogosyan, D. How filaments are woven into the cosmic web. Nature, 380:603–606, 1996.
- (7) Mead, J.M.G., King, L.J., and McCarthy, I.G. Probing the cosmic web: inter-cluster filament detection using gravitational lensing. Mon. Not. Roy. Astron. Soc., 401:2257, 2010.
- (8) Maturi, M. and Julian, M. A&A, 559(A112), 2013.
- (9) Clampitt, J., Jain, B., Takada, M., and Miyatake, H. Detection of Stacked Filament Lensing Between SDSS Luminous Red Galaxies. Mon. Not. Roy. Astron. Soc., 457(3):2391–2400, 2016.
- (10) K. M. Smith, O. Zahn, and O. Dore. Detection of Gravitational Lensing in the Cosmic Microwave Background. Phys. Rev., D76:043510, 2007.
- (11) C. M. Hirata, S. Ho, N. Padmanabhan, U. Seljak, and N. A. Bahcall. Correlation of CMB with large-scale structure: II. Weak lensing. Phys. Rev., D78:043520, 2008.
- (12) L. E. Bleem et al. A Measurement of the Correlation of Galaxy Surveys with CMB Lensing Convergence Maps from the South Pole Telescope. Astrophys. J., 753:L9, 2012.
- (13) B. D. Sherwin et al. The Atacama Cosmology Telescope: Cross-Correlation of CMB Lensing and Quasars. Phys. Rev., D86:083006, 2012.
- (14) S. Ferraro, B. D. Sherwin, and D. N. Spergel. WISE measurement of the integrated Sachs-Wolfe effect. Phys. Rev., D91(8):083533, 2015.
- (15) R. Allison et al. The Atacama Cosmology Telescope: measuring radio galaxy bias through cross-correlation with lensing. Mon. Not. Roy. Astron. Soc., 451(1):849–858, 2015.
- (16) T. Giannantonio et al. CMB lensing tomography with the DES Science Verification galaxies. Mon. Not. Roy. Astron. Soc., 456(3):3213–3244, 2016.
- (17) Pullen, A.R., Alam, S., He, S., and Ho, S. Constraining Gravity at the Largest Scales through CMB Lensing and Galaxy Velocities. Mon. Not. Roy. Astron. Soc., 460(4):4098–4108, 2016.
- (18) C. Doux, E. Schaan, E. Aubourg, K. Ganga, K. Lee, D. N. Spergel, and J. TrÃ©guer. First detection of cosmic microwave background lensing and Lyman- forest bispectrum. Phys. Rev., D94(10):103506, 2016.
- (19) S. Singh, R. Mandelbaum, and J. R. Brownstein. Cross-correlating Planck CMB lensing with SDSS: Lensing-lensing and galaxy-lensing cross-correlations. Mon. Not. Roy. Astron. Soc., 464(2):2120–2138, 2017.
- (20) Chen, Y., Ho, S., Brinkmann, J., Freeman, P.E., Genovese, C.R., Schneider, D.P., and Wasserman, L. Cosmic Web Reconstruction through Density Ridges: Catalogue. Mon. Not. Roy. Astron. Soc., 461(4):3896–3909, 2016.
- (21) Dawson, K.S. et al. The Baryon Oscillation Spectroscopic Survey of SDSS-III. Astron. J., 145:10, 2013.
- (22) Alam, S. et al. The Eleventh and Twelfth Data Releases of the Sloan Digital Sky Survey: Final Data from SDSS-III. Astrophys. J. Suppl., 219(1):12, 2015.
- (23) Chen, Y., Ho, S., Freeman, P.E., Genovese, C.R., and Wasserman, L. Cosmic Web Reconstruction through Density Ridges: Method and Algorithm. Mon. Not. Roy. Astron. Soc., 454(1):1140–1156, 2015.
- (24) Adam, R. et al. Planck 2015 results. I. Overview of products and scientific results. Astron. Astrophys., 594:A1, 2016.
- (25) Ade, P. A. R. et al. Planck 2015 results. XV. Gravitational lensing. Astron. Astrophys., 594:A15, 2016.
- (26) S. D. Epps and M. J. Hudson. The Weak Lensing Masses of Filaments between Luminous Red Galaxies. Mon. Not. Roy. Astron. Soc., 468:2605–2613, March 2017.
- (27) Y. Higuchi, M. Oguri, and M. Shirasaki. Statistical properties of filaments in weak gravitational lensing. Mon. Not. Roy. Astron. Soc., 441(1):745–756, 2014.
- (28) Ade, P.A.R. et al. Planck 2015 results. XIII. Cosmological parameters. Astron. Astrophys., 594:A13, 2016.
- (29) Chen, Y., Ho, S., Mandelbaum, R., Bahcall, N.A., Brownstein, J.R., Freeman, P.E., Genovese, C.R., Schneider, D.P., and Wasserman, L. Detecting Effects of Filaments on Galaxy Properties in the Sloan Digital Sky Survey III. Mon. Not. Roy. Astron. Soc., 466:1880, 2017.
- (30) Chen, Y. A Tutorial on Kernel Density Estimation and Recent Advances. ArXiv e-prints, April 2017.
- (31) Ozertem, U. and Erdogmus, D. Locally defined principal curves and surfaces. J. Mach. Learn. Res., 12:1249–1286, jul 2011.
- (32) Bagla, J.S. A TreePM code for cosmological N-body simulations. J. Astrophys. Astron., 23:185–196, 2002.
- (33) White, M.J., Hernquist, L., and Springel, V. Simulating the Sunyaev-Zel’dovich effect(s): Including radiative cooling and energy injection by galactic winds. Astrophys. J., 579:16, 2002.
- (34) Reid, B.A., Seo, H., Leauthaud, A., Tinker, J.L., and White, M. A 2.5 per cent measurement of the growth rate from small-scale redshift space clustering of SDSS-III CMASS galaxies. Mon. Not. Roy. Astron. Soc., 444(1):476–502, 2014.
- (35) Peacock, J.A. and Smith, R.E. Halo occupation numbers and galaxy bias. Mon. Not. Roy. Astron. Soc., 318:1144, 2000.
- (36) Seljak, U. Analytic model for galaxy and dark matter clustering. Mon. Not. Roy. Astron. Soc., 318:203, 2000.
- (37) Benson, A.J., Cole, S., Frenk, C.S., Baugh, C.M., and Lacey, C.G. The Nature of galaxy bias and clustering. Mon. Not. Roy. Astron. Soc., 311:793–808, 2000.
- (38) White, M.J., Hernquist, L., and Springel, V. The Halo model and numerical simulations. Astrophys. J., 550:L129–L132, 2001.
- (39) Berlind, A.A. and Weinberg, D.H. The Halo occupation distribution: Towards an empirical determination of the relation between galaxies and mass. Astrophys. J., 575:587–616, 2002.
- (40) Cooray, A. and Sheth, R.K. Halo models of large scale structure. Phys. Rept., 372:1–129, 2002.
- (41) Beutler, F. et al. The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Testing gravity with redshift-space distortions using the power spectrum multipoles. Mon. Not. Roy. Astron. Soc., 443(2):1065–1089, 2014.
- (42) White, M., Tinker, J.L., and McBride, C. K. Mock galaxy catalogues using the quick particle mesh method. Mon. Not. Roy. Astron. Soc., 437(3):2594–2606, 2014.
- (43) Carlson, J. and White, M. Embedding realistic surveys in simulations through volume remapping. Astrophys. J. Suppl., 190:311–314, 2010.
- (44) Gorski, K.M., Hivon, E., Banday, A.J., Wandelt, B.D., Hansen, F.K., Reinecke, M., and Bartelman, M. HEALPix - A Framework for high resolution discretization, and fast analysis of data distributed on the sphere. Astrophys. J., 622:759–771, 2005.

Acknowledgements We thank Anthony Pullen and Elena Giusarma for helpful discussion, Martin White for providing us the N-body simulations, Alex Krolewski, Benjamin Horowitz for comments on the draft. S.Ho is supported by NASA and DOE for this work. S.He is supported by NSF-AST1517593 for this work. S.A. is supported by the European Research Council through the COSFORM Research Grant (#670193). S.F. thanks the Miller Institute for Basic Research in Science at the University of California, Berkeley for support. Some of the results in this paper have been derived using the HEALPix package 2005ApJ…622..759G ().

Author Contributions S.He led the project and most of the manuscript writing. S.A. provided the sky mocks for galaxies and dark matter particles as well as wrote the text relative to sky mock for Filaments and Dark Matter in the Method section. S.F. helped with the theoretical modeling and the interpretation of the results, as well as writing part of the manuscript. Y.C. provided the filament intensity maps for data and simulations. S.Ho conceived the idea of cross-correlating filaments with CMB lensing. All authors contributed to the interpretation of the data and commented on the manuscript.

Competing interests Authors declare no competing financial interests.

## Appendix A Supplementary information

### a.1 Error calculation for from simulation

We derive the relation between and in eq. (10). By removing the appropriate functions from the integrands, which are slowly varying compared to , the correct expression is

However, the approximations required to produce this expression are not perfect, causing the estimation of from simulation to slightly deviate from the true value of . We estimate the deviation by relating the theoretical prediction for (eq. (7)) and (eq. (9)) by the following equation

(19) |

The result is shown in Supplement Fig. 1.

We see that the is less than 5% from unity, which is much smaller than , where is the error for . Thus, the approximation for converting to only causes a negligible bias.