Issue 
A&A
Volume 633, January 2020



Article Number  A89  
Number of page(s)  12  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201936678  
Published online  15 January 2020 
A gravitational lensing detection of filamentary structures connecting luminous red galaxies
^{1}
Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK
email: qx211@roe.ac.uk
^{2}
Department of Astrophysics, University of Oxford, Keble Road, Oxford OX1 3RH, UK
email: naomi.robertson@physics.ox.ac.uk
^{3}
German Centre for Cosmological Lensing, Astronomisches Institut, RuhrUniversität Bochum, Universitätsstr, Bochum, Germany
^{4}
Kavli Institute for Particle Astrophysics & Cosmology, Stanford University, PO Box 2450, Stanford, CA 94305, USA
^{5}
ArgelanderInstitut für Astronomie, Universitẗ Bonn, Auf dem Hügel 71, 53121 Bonn, Germany
^{6}
Leiden Observatory, Leiden University, PO Box 9513, 2300 Leiden, The Netherlands
^{7}
Instituto de Física, Pontificia Universidad Católica de Valparaíso, Casilla, 4059 Valparaíso, Chile
Received:
12
September
2019
Accepted:
26
November
2019
We present a weak lensing detection of filamentary structures in the cosmic web, combining data from the KiloDegree Survey, the Red Cluster Sequence Lensing Survey, and the CanadaFranceHawaii Telescope Lensing Survey. The line connecting luminous red galaxies with a separation of 3 − 5 h^{−1} Mpc was chosen as a proxy for the location of filaments. We measured the average weak lensing shear around ∼11 000 candidate filaments selected in this way from the Sloan Digital Sky Survey. After nulling the shear induced by the dark matter haloes around each galaxy, we reported a 3.4σ detection of an anisotropic shear signal from the matter that connects them. Adopting a filament density profile, motivated from Nbody simulations, the average density at the centre of these filamentary structures was found to be 15 ± 4 times the critical density.
Key words: gravitational lensing: weak / largescale structure of Universe / dark matter / cosmology: observations
© ESO 2020
1. Introduction
Galaxy surveys, including the 2dF Galaxy Redshift Survey (Colless et al. 2001) and the Sloan Digital Sky Survey (SDSS; Zehavi et al. 2011), have shown that visible matter in our Universe is not uniformly distributed on intermediate scales ∼100 h^{−1} Mpc. Instead, a weblike structure is observed with clusters of galaxies identifying the densest regions. Nbody simulations predict the existence of these largescale structures (e.g. Bond et al. 1996; Springel et al. 2005), suggesting a hierarchical structure formation for the cosmic web. We can classify the web (e.g. Eardley et al. 2015) into regions of clusters, filaments, sheets and voids. In this cosmic web, large underdense regions (voids) are enclosed by anisotropically collapsed surface structures (sheets) and line structures (filaments) which intersect at the most overdense isotropic regions (clusters). The Zel’dovich approximation predicts that ∼42% of the mass of the Universe is in a filament environment (Zel’dovich 1970), and this has been confirmed by simulations (AragónCalvo et al. 2010). However, as filament environments do not display a very high density contrast, this makes direct observations challenging.
One way to observe filaments is from the Xray emission induced by the warm hot intergalactic medium (WHIM) with several intercluster filaments investigated in this way (Briel & Henry 1995; Kull & Böhringer 1999; Werner et al. 2008). There are also reported detections of filaments using overdensities of galaxies (Pimbblet & Drinkwater 2004; Ebeling et al. 2004). Recently, two independent studies (de Graaff et al. 2019; Tanimura et al. 2019) detected the Sunyaev–Zel’dovich (SZ) signal from the ionised gas in the cosmic web. They estimated the density of ionised gas to be ∼(28 ± 12)% of the total baryon density in the Universe, close to resolving the missing baryon problem (Bregman 2007).
In this paper, we investigate the use of weak gravitational lensing to detect filaments. Based on the distortion of light rays around massive objects, gravitational lensing probes the total mass traced by the largescale structures and is, therefore, highly complementary to the SZ detection of the gas. Although Dietrich et al. (2012) made a direct weak lensing detection of a filament connecting two massive Abell clusters, the direct detection of typical individual filaments is limited by the low signaltonoise measurement, and most studies stack large samples of candidate filaments and analyse the resulting average weak lensing signal instead (Mead et al. 2010; Clampitt et al. 2016; Epps & Hudson 2017; Kondo et al. 2019).
Clampitt et al. (2016; hereafter C16) determined the weak lensing signal at around 135 000 pairs of SDSS Luminous Red Galaxies (LRGs) with a projected separation 6 h^{−1} Mpc ≤ R_{sep} ≤ 14 h^{−1} Mpc and a redshift separation of Δz < 0.004. Using a “nulling” estimator that cancels the spherically symmetric contribution of the LRG haloes in the shear measurement, they reported a 4.5σ detection of the filament lensing signal. In another study, Epps & Hudson (2017) used ∼23 000 pairs of LRGs from the Baryon Oscillation Spectroscopic Survey (BOSS) LOWZ and CMASS samples as tracers of filaments. Using data from the CanadaFranceHawaii Telescope Lensing Survey (CFHTLenS), they performed a mass reconstruction of a set of stacked LRG pairs with a projected angular separation between 6 and 10 h^{−1} Mpc and a redshift separation Δz < 0.003. After subtracting the signal from a mass reconstruction of a set of stacked LRG pairs with the same separation on the sky but a greater redshift separation (0.033 < Δz < 0.04) such that haloes should not be physically connected, they reported a 5σ detection of a filament lensing signal. A more recent study (Kondo et al. 2019) used 70 210 pairs of LRGs from the CMASS sample with a projected separation between 6 and 14 h^{−1} Mpc and a lineofsight separation of less than 6 h^{−1} Mpc. Using the Subaru Hyper SuprimeCam (HSC) firstyear galaxy shape catalogue and adopting the C16 nulling approach they reported 3.9σ detection of a filament signal.
We note that the methodology taken in these previous studies can be understood as a threepoint galaxy–galaxyshear correlation function conditioned on specific intervals of separation between lens galaxies (Schneider & Watts 2005). The full suite of the galaxygalaxygalaxy lensing (GGGL) statistics have been applied to the RedSequence Cluster Survey (Simon et al. 2008) and CFHTLenS (Simon et al. 2019) to measure the excess mass around galaxy pairs separated by ≲300 h^{−1} kpc.
In this paper, we present the weak lensing signal measured between 11 706 LOWZ LRG pairs that have a separation of 3 − 5 h^{−1} Mpc, combining three public weak lensing surveys; the KiDS+VIKING450 survey (KV450; Hildebrandt et al. 2018; Wright et al. 2019), the Red Cluster Sequence Lensing Survey (RCSLenS; Hildebrandt et al. 2016), and the CFHTLenS (Heymans et al. 2012). We improve the nulling methodology described in C16 to deal with contamination from filament tracers and use a large suite of Nbody simulations to validate our pipeline and compare our results. A standard ΛCDM cosmology has been adopted throughout this study to calculate distances with a matter density Ω_{m} = 0.3, energy density Ω_{Λ} = 0.7, effective number of neutrino species N_{eff} = 3.04, baryon density Ω_{b} = 0.0 and current Hubble constant H_{0} = 100 h km s^{−1} Mpc^{−1} where h is the Hubble parameter h = 0.7.
This paper is structured as follows. In Sect. 2, we describe the survey data and simulations. Section 3 summarises the weak lensing formalism, the adopted filament model, and methodology. We show our results in Sect. 4 and conclude in Sect. 5. In Appendix A, we present a validation of our nulling technique. In Appendix B, we document the spherical rotation methodology that is required for high declination surveys.
2. Surveys and simulations
2.1. The lensing surveys
The properties of the three lensing surveys, KV450, RCSLenS, and CFHTLenS, are listed in Table 1. They share a similar data processing pipeline, where the shape measurement of galaxies was conducted using the lensfit model fitting code (Miller et al. 2013). This approach convolves the pixelised model PointSpreadFunction with an analytical surface brightness model consisting of bulge and disk components. It uses model fitting to estimate galaxy ellipticitices and with an associated inverse variance weight, w_{s}. The (reduced) shear (cf. Eq. (4)) is then given by the weighted average of ellipticities, . The observed shear is biased with respect to the true shear and is typically described by the linear bias model (Heymans et al. 2006) as
Properties of each of the three lensing surveys used in this analysis.
In all of these three surveys, the shear multiplicative bias terms were characterised as a function of the signaltonoise ratio and size of the galaxies, thereby allowing us to calculate the bias for an arbitrary selection of galaxies. The correction for this multiplicative bias was carried out as it was in Velander et al. (2014).
Photometric redshifts, z_{B}, were estimated using the Bayesian photometric redshift algorithm (BPZ; Benítez 2000) as detailed in Hildebrandt et al. (2012). Wright et al. (2019) and Hildebrandt et al. (2018) show how the photometric redshifts distributions for KV450 are then calibrated using the “weighted direct calibration” method, with weights estimated using a deep spectroscopic training sample in 9band ugriZYJHK_{s} magnitude space. No such calibration was performed for the 4band (RCSLenS) or 5band (CFHTLenS) surveys, and instead a probability distribution of true redshifts was estimated from the sum of the BPZ redshift probability distributions. This approach has been demonstrated to carry more systematic error (Choi et al. 2016; Hildebrandt et al. 2017). We discuss how we take this redshift uncertainty into account in our final analysis in Sect. 3.4.
2.2. The BOSS survey
We used the Baryon Oscillation Spectroscopic Survey (BOSS) galaxies from the 12th SDSS Data Release (Alam et al. 2015) to define candidate filaments. Among all three lensing surveys RCSLenS has the most SDSS overlap with almost double that of KV450 or CFHTLenS. Once a robust photometric redshift selection has been applied, however, RCSLenS has only 20%/30% the lensing source density in comparison to CFHTLenS/KV450, (see Table 1 for details).
Given the depth of the lensing surveys and the uncertainty in the high redshift tail of the redshift distribution for CFHTLenS and RCSLenS, we chose to limit our analysis to the LOWZ sample, selected based on colour and magnitude, using a redshift cut 0.15 < z < 0.43 (Ross et al. 2012). We did not consider the higher redshift CMASS sample. The typical virial halo mass of LOWZ galaxies is ∼5.2 × 10^{13} h^{−1} M_{⊙} (Parejko et al. 2013). These haloes have a typical virial radius ∼1 h^{−1} Mpc.
2.3. Simulations
We used the Scinet Light Cone Simulations (SLICS^{1}; HarnoisDéraps & van Waerbeke 2015; HarnoisDéraps et al. 2018) to test our methodology. This suite provided us with 819 independent light cones on a 10 × 10 deg^{2} patch of the sky. Each light cone was constructed from the full nonlinear evolution of 1536^{3} particles with m_{p} = 2.88 × 10^{9} h^{−1} M_{⊙}, within a 505^{3}(h^{−1} Mpc)^{3} cube. Particles were then projected onto mass sheets at 18 redshifts between 0 < z < 3, and subsequently inspected to identify dark matter haloes. For each simulation, 100 deg^{2} light cone mass sheets and haloes were extracted; the former were then raytraced into lensing shear maps, while the latter were used to generate mock LOWZ galaxies with a halo occupation distribution, optimised so that the clustering of mock galaxies would be consistent with the LOWZ data (see HarnoisDéraps et al. 2018 for details). We find the mock sample variance to be in good agreement with the Jackknife errors measured directly from the data. Source galaxy positions were drawn at random, with the shear assigned for a range of redshifts and high number density. We next randomly drew from the SLICS mock source galaxy sample so as to match the corresponding number density and redshift distribution of each of the three surveys, KV450, RCSLenS, and CFHTLenS. Intrinsic galaxy shapes were chosen to match the KiDS ellipticity dispersion (Hildebrandt et al. 2017), which is good description of the ellipticity dispersion that is also found in RCSLenS and CFHTLenS.
3. Summary of weak lensing formalism and methodology
In this section, we summarise the weak gravitational lensing theory following the more detailed derivations in Bartelmann & Schneider (2001). Assuming the thin lens approximation, a foreground object at a position θ has a 2D comoving surface mass density Σ(θ). The convergence is then defined as
where Σ_{crit} is the comoving critical surface density in a flat Universe given by
Here χ is the comoving distance and z_{l}, z_{s} are the redshifts of the lens and source, respectively. Since we are interested in the large scale comoving surface density around filaments, this is the appropriate choice for Σ_{crit} (see Dvornik et al. 2018, for a discussion on the different definitions of Σ_{crit}).
The deflection potential, ψ(θ), is connected to the convergence, κ, via Poisson’s equation ∇^{2}ψ = 2κ, and the complex shear is related to the second derivatives of the deflection potential via
where x_{1}, x_{2} are the horizontal and vertical displacements on the projected sky.
For a filament aligned with the x_{1} axis if we assume the deflection potential ψ and convergence κ are both invariant along the filament, then Eq. (4) immediately implies that partial derivatives with respect to the x_{1}axis will be equal to zero. This leads to the approximation that for filaments, we should expect to measure γ_{1} ≈ −κ and γ_{2} ≈ 0. Motivated by the simulation results of Colberg et al. (2005) and Mead et al. (2010) considered the power law density profile around filaments and suggested the model for the convergence at a distance r from filament centre, measured perpendicular to the major filament axis (x_{1}axis)
Here κ_{c} is the amplitude of the convergence at the filament centre (r = 0) and r_{c} is the halfmaximum radius of the density profile.
3.1. Filament candidiate
Colberg et al. (2005) showed that cluster pairs separated by < 5 h^{−1} Mpc are always connected by filamentary structures. We therefore select luminous red galaxy (LRG) pairs in the LOWZ catalogue with redshift separation δz < 0.002 and a projected separation 3 h^{−1} Mpc ≤ R_{sep} ≤ 5 h^{−1} Mpc as our candidate filaments^{2} which we will refer to as our physical pairs (PP). Nonphysical pairs (NP) are defined to have the same projected separation range but with large lineofsight separations with 0.033 < δz < 0.04 (corresponding to ∼100 h^{−1} Mpc). With such a large physical separation we would not expect to detect a filament signal. The NP therefore provide an important nulltest for our methodology.
Our candidates differ from the selection made by C16, Epps & Hudson (2017), de Graaff et al. (2019), Tanimura et al. (2019), and Kondo et al. (2019), who focused on separations of 6 − 10 h^{−1} Mpc. Our choice maximises signaltonoise, as shown in our analysis of numerical simulations in Sect. 4.1, but for completeness we also present an analysis of 6 − 10 h^{−1} Mpc filaments in Sect. 5.
3.2. Stacking method
For each lens filament candidate at redshift z_{l} = z_{f} = (z_{lens 1} + z_{lens 2})/2, we measure the Σ_{crit}weighted shear, E_{f}, on a grid (i, j) centred and oriented with the pair of LRGs, where
and the sum is taken over all sources^{3}, s, with z_{B} > z_{l} + 0.1 and
The complex ellipticity is the observed ellipticity of the source rotated into the reference frame where the filament lies along the x_{1} axis. The Σ_{crit} weight converts from a shear estimate to an estimate of the surface mass density Σ in Eq. (2). The grid (i, j) has an extent [ − 2R_{sep}, 2R_{sep}]×[−2R_{sep}, 2R_{sep}] and 129^{2} pixels, and the pair of lens galaxies that define the filament candidate are positioned to lie at the centre of the pixels at ( − 0.5R_{sep}, 0) and (0.5R_{sep}, 0).
We also construct a corresponding weight map for each filament candidate,
where the extra factor of provides optimal signaltonoise weighting (Velander et al. 2011).
When rotating each filament pair into a common reference frame, we note that at high declination, the tangent plane method used in Epps & Hudson (2017) and the direct cartesian approximation in C16 results in nonuniform grid cells. As the nulling approach requires a flat geometry on the grid, we found that these approximations lead to a biased result. To solve this problem for highdeclination patches, we use the spherical rotation method from de Graaff et al. (2019). This process is detailed in Appendix B, and illustrated in Fig. B.1, with the rotated shear map defined in Eq. (B.4).
As we have spectroscopic redshifts for the filaments but only photometric redshifts for the sources, the inverse critical surface mass density is calculated for each survey as
where is the probability distribution of the true redshift of the source galaxies that enter the measurement
For CFHTLenS and RCSLenS we use the persource p_{s}(z_{s}z_{B}) provided by each survey, even though this has been shown to introduce biases (Choi et al. 2016), which we account for in Sect. 3.4. For KV450 we use the weighted direct calibration method of Hildebrandt et al. (2017) to determine the source redshift distribution directly for an ensemble of sources. In practice we calculate in Eq. (10) at eight z_{l} values and interpolate to evaluate at each filament redshift. The effective n(z_{s}) for each survey, given by
is shown in Fig. 1, with CFHTLenS and KV450 providing a deeper source redshift than RCSLenS.
Fig. 1. The effective n(z) of all three lensing surveys, as defined in Eq. (12), when using LOWZ galaxies as lenses. Each curve is normalised such that ∫n(z) dz = 1. 
We correct the measured shear signal E_{f}, with the signal measured around “random” filaments. This is now a standard procedure in galaxy–galaxy lensing studies (e.g. Mandelbaum et al. 2005) which removes any hidden systematics and reduces sampling variance noise. We create random filament samples for each survey patch, listed in Table 2, by randomly shifting filaments within the same patch while preserving their redshifts, position angles and number density in the patch. As we do not expect any physical signal from the random catalogue, we subtract any measured “random” signal from the data as follows:
Effective area and number of filaments in each survey patch.
Here, for a patch with N_{p} filament candidates is given by
and
where E_{r} and W_{r} are the weighted shear (Eq. (6)) and weight (Eq. (8)) but measured around a random pair. w_{ran} is the normalisation weight where N_{R} is the number of realisations which ensures N_{R} × N_{p} exceeds 100 000 in each patch^{4}. This ensures that the random signal has low scatter so that we can take the mean as the random correction.
The randomcorrected shear map and weight map W_{f} are then optimally combined over all filament candidates to determine the total weighted shear signal, 𝒯, over the full sample,
where w_{f} = w_{lens 1} × w_{lens 2} is the product of the recommended SDSS completeness weights for the LOWZ galaxies to account for any extra galaxy pairs in the case of fibre collision. We have also applied the multiplicative calibration correction at this stage which is given by Velander et al. (2014) as
The total weighted shear 𝒯(i, j) is a combination of both the shear contribution from the haloes surrounding the LOWZ LRGs and the contribution from any filament that connects them.
In order to isolate the filament we apply the “nulling” procedure described in Appendix A to get a final measurement of the shear contribution from the filament only, ℱ(r), as a function of the distance, r, from the central filament axis,
Here 𝒩 is the nulling operator given in Eq. (A.11) and the summation over pixels from an i_{min} to i_{max} runs along the filament from −0.438R_{sep} to 0.438R_{sep}. This value was found to minimise any residual contribution from the haloes positioned at ±0.5R_{sep}, that remains after a nulling analysis of the SLICS simulation (see Sect. 4.1). Our nulling operator 𝒩 combines the shear values measured at 8 different positions (including 4 positions from a reflection about the filament axis) which alternatingly rotate around the two haloes. In this way the isotropic contribution from the parent haloes sum to zero (i.e. “null”) and any anisotropic contribution inbetween the two haloes can be recovered. In Appendix A, we provide a detailed proof and compare our nulling approach to that adopted in C16. Through tests on a fiducial model, we show that the C16 nulling approach produces a biased result on large scales.
In this derivation we have carried both components of the shear with ℱ = ℱ_{γ1} + iℱ_{γ2}. Given the Mead et al. (2010) filament model, where κ = −γ_{1}, we expect ℱ_{γ2} = 0. We will fit ℱ_{γ1} using the twoparameter model in Eq. (5) with the amplitude parameter replaced by ℱ_{c} which is equivalent to Σ_{crit}κ_{c} for a single lens–source pair.
3.3. Error estimation from SLICS
In order to estimate the error on the measured signal from observations, we use a large number of independent and representative lensing simulations from SLICS for each of the three surveys. SLICS allows us to correctly account for the sampling variance, which was found to be the dominant source of noise in Kondo et al. (2019). The source sample of galaxies differs for each filament pair owing to our source selection that z_{B} − z_{f} > 0.1. For KV450 we can apply this source selection criteria accurately as the SLICS simulations include mock KV450 photometric redshifts that reproduce the scatter, bias and catastrophic outlier populations found in the KV450 data (HarnoisDéraps et al. 2018). For CFHTLenS and RCSLenS, this information is not encoded. We therefore create mocks from SLICS, modelling source samples for four different filament bins with (z_{min}, z_{max}) = [0.1, 0.2],[0.2, 0.3],[0.3, 0.4], and [0.4, 0.5], respectively. For each filament bin, we calculate the source galaxy redshift distribution n(z), using Eq. (12) and the effective galaxy number density for sources with z_{B} − z_{max} > 0.1. We then populate 500 independent simulations using these distributions, and measure and combine the weighted shear and weight maps for each of the four filament bins. For KV450 we are able to verify that this binned approach is consistent to the unbinned methodology applied to the data using the KV450 SLICS simulations.
The covariance matrix of the signal is reweighted by the number of filament candidates for each survey, n_{fil, survey}, and estimated from the SLICS simulations as
where is the filament signal (Eq. (18)) averaged over all N_{sim} = 500 surveyspecific SLICS simulations, and is the average number of filament candidates in these simulations. The covariance is then used to calculate the χ^{2} when estimating parameters in the filament model as
where is our filament model defined in Eq. (5), calculated on a fine grid of parameters (ℱ_{c}, r_{c}). In analogy to Eq. (19), we also define a covariance for the ℱ_{γ2} component which serves as a systematic nulltest for our analysis.
We found that a simple bootstrap error analysis of the data, where the set of maps are resampled with repetitions before stacking, underestimates the true measurement error. This approach misses the sampling variance term which, like Kondo et al. (2019), we find is a significant component to the error for smallarea surveys such as KV450 and CFHTLenS.
3.4. Accounting for uncertainty in the redshift distributions
As discussed in Sect. 2.1, the probability distribution from Eq. (11) have not been calibrated for RCSLenS and CFHTLenS. A systematic uncertainty on the resulting n(z_{s}) is thus expected. In order to take this into account, we use a nuisance parameter δz_{s} = 0.1 for RCSLenS and δz_{s} = 0.04 for CFHTLenS that captures the p(z_{s}) uncertainty determined by Choi et al. (2016). For KV450 we use δz_{s} = 0.025 following Wright et al. (2019). We shift the by ±δz_{s} in Eq. (10) to yield two new functions , and repeat the full measurement and error analysis using both the and calibration. We then estimate the likelihood of the data D given the model (ℱ_{c}, r_{c}), assuming a prior on the uncertainty δz_{s} which consists of 3 δfunctions and marginalised over as
A full marginalisation where many samples are taken at different redshift offsets, spanning the δz range, is unfortunately unfeasible given the complexity of the measurement pipeline.
4. Results
4.1. Filaments in SLICS
We validated our pipeline using the SLICS simulations of mock LOWZ lens galaxies and mock KiDS sources. We analysed both a noisefree catalogue and a catalogue with shape noise and mock photometric redshifts. The result is shown in Fig. 2 where the upper panels show measurements of ℱ_{γ1}, and the lower panels show measurements of ℱ_{γ2}. The left and right columns correspond to results from the 3 − 5 h^{−1} Mpc and 6 − 10 h^{−1} Mpc filament candidates, respectively.
Fig. 2. Weighted shear measured between LRG pairs in the SLICS simulation of a random 400 deg^{2} degree survey. Left: results from 3 to 5 h^{−1} Mpc filament candidates. Upper panel: average surface mass density ℱ_{γ1}, and the lower panel shows the measured crossshear component. On both panels, the blue data shows the result for a KiDSlike survey depth and shape noise, and the orange data points show the measurement for a noisefree simulation with the errorbar given by the error on the mean of all 158 realisations. Right: set of results from 6 to 10 h^{−1} Mpc filament candidates. 
For the noisefree simulations, we find that ℱ_{γ2} is consistent with zero for both the 3 − 5 h^{−1} Mpc and 6 − 10 h^{−1} Mpc length filaments. This demonstrates that our nulling procedure correctly removes the contribution from the LRG haloes in the analysis. However, for the simulations with shape noise, both the ℱ_{γ1} and ℱ_{γ2} measurements are consistent with zero, which suggests that even though KV450 is deeper than the KiDS450 data simulated in SLICS, we should not expect a significant detection from KiDS alone. As the variance in the noisefree simulation reflects the level of sample variance, we also report that, by measuring the noise level from these two sets of simulations, the sample variance is comparable with the shape noise.
Constraining the parameters of Mead et al. (2010) model with the noisefree SLICS results we find and . The and demonstrate that the model is a good fit to the data (ν = 13 degrees of freedom). We find that the surface mass density of the filament is a factor of 2.5 smaller for the 6 − 10 h^{−1} Mpc filament and will therefore be more challenging to detect using gravitational lensing.
With the noisefree simulations we are able to analyse whether our signal depends on the redshift of the filament. We constrain the amplitude ℱ_{c} and the scale r_{c} parameters of the filament model for 4 redshift quantiles of the SLICS LOWZ filament samples using the same background sources, with the result shown in Fig. 3. The choice of parameterisation on the axes is motivated by the filament model equation as well as for visual simplicity. Here we see significant differences between the samples which could be caused by an evolution in the bias of the LOWZlike galaxy samples in the SLICS mocks, or the evolution of the filament density field. Given the different source redshift distributions of the three lensing surveys (see Fig. 1) which makes the effective redshift of the average lens differ, this result suggests that we should not necessarily expect the results of these surveys to agree perfectly.
Fig. 3. Redshift evolution of the filament signal in the noisefree SLICS simulation. Left panel: ℱ_{γ1} signal measured from filament candidates in SLICS for four redshift quantiles. Right panel: 68% and 95% confidence region of the model parameters (ℱ_{c}, r_{c}) from the corresponding signal in the left panel, with the result from all samples combined shown in grey. 
4.2. The detection of filaments with KV450, RCSLenS and CFHTLenS
Figure 4 presents our filament shear measurements and constraints on the two parameters of the filament model (Eq. (5)) for physical pairs, our filament candidates (upper panels), and nonphysical pairs, our control sample (lower panels). The left panel shows the nulled ℱ(r) shear signal as a function of the distance from the centre of the filament measured in units of h M_{⊙} pc^{−2}. The result is presented for each lensing survey. We also show the measurements from the threesurveys combined using inverse variance weighting. The blue data points are a measurement of ℱ_{γ1}, whereas the orange data points show the nulltest ℱ_{γ2}. The shaded region corresponds to the statistical noise from our fiducial measurements, and the capped errorbars correspond to the systematic uncertainty captured by the photometric redshift bias nuisance parameter δz (see Sect. 3.4). The right panel shows the 68% and 95% confidence region of parameters ℱ_{c} and r_{c} in Mead et al. (2010) filament model. These estimated parameters can be compared to the bestfit parameter from the noisefree SLICS analysis in Sect. 4.1 which is also represented by the cross in the right panel. We note that the noisefree SLICS bestfit is consistent with all three surveys. We also present joint constraints from the combined signal using a block covariance
Fig. 4. The detection of the cosmic web between neighbouring luminous red galaxies as detected through the weak lensing of background galaxies from different lensing surveys. Left: xaxis is the distance measured perpendicular to the filament axis scaled such that 1 is equivalent to the projected separation between the pair of LRGs. The yaxis is the nulled shear signal where ℱ_{γ1} (blue data points) measures the average surface mass density of the filament, and ℱ_{γ2} is expected to be consistent with zero and hence serves as a null test. Lower right small panel: measurements from the threesurveys combined using inverse variance weighting. We note that this additional panel is purely for illustration, however, as our joint surveyconstraints on the filament model are derived from a combination of the surveys on the likelihoodlevel. Right: estimated parameters in the filament density model Eq. (5) from the stacked signal for all surveys individually and their combination. Upper: results from 3 ∼ 5 h^{−1} Mpc Physical Pairs, our filament candidates; Lower: results from 3 ∼ 5 h^{−1} Mpc nonphysical Pairs, our control sample. 
and extended data vector ℱ_{ext} = (ℱ_{CF}, ℱ_{RC}, ℱ_{KV})^{T}, providing an estimate of the average filament profile from all the filament candidates across the three surveys. Equation (22) assumes that the surveys are uncorrelated, which is a good approximation to make given the lack of overlap between the different surveys.
To quantify the significance of our measurements, we use the likelihood ratio test between the null hypothesis H_{0} and the filament model H_{1}, where the likelihood ratio LR is
where B_{0} and B_{1} are the parameter space in each hypothesis, i.e., B_{0} has no free parameters and B_{1} = {ℱ_{c}, r_{c}}. By Wilks’ theorem (Wilks 1938; Williams 2001), the deviance defined as Dev = 2lnLR has an asymptotic chisquared distribution with dim(B_{1})−dim(B_{0}) = 2 degrees of freedom when H_{0} is true. Estimating the maximum likelihoods from and using Eq. (21) and computing the deviance, we report the significance level for each individual survey as well as the combined analysis in Table 3. The reported suggests our model is a reasonable fit to the data in all cases even for KV450, where p(χ^{2} < 7.40ν = 13) = 0.12. We find the bestfit model parameters for ℱ_{γ1} from all three surveys combined as ℱ_{c} = 10.5 ± 2.9 h M_{⊙} pc^{−2} and . We note that, the majority of the detection derives from CFHTLenS alone with a 3.1σ detection. Combining all three surveys we measure a 3.4σ detection of the filament weak lensing signal. CFHTLenS is the most constraining survey as it combines both survey depth with significant SDSS overlap. KV450, with roughly half the source density of CFHTLenS, and RCSLenS, with roughly 20% of the source density, will only start to add significant constraining power with the inclusion of additional overlapping SDSS area.
χ^{2} value and pvalue for all computed signals from each individual survey and their combination.
Our control sample of “nonphysical pairs” (NP) were selected to be pairs of lens galaxies with projected separations 3 − 5 h^{−1} Mpc, but distant in redshift space (0.033 < Δz < 0.04). These nonphysical pairs would not be connected by a filament, hence providing an important validation of our nulling approach to isolate the filament signal. We find that the measured signal for ℱ_{γ1, 2} is consistent with zero for all surveys and the combined survey as shown in Table 3 and the lower panel of Fig. 4.
For consistency with other analyses in the literature we also analysed 27 880 filament pairs in LOWZ that have a physical separation of 6 − 10 h^{−1} Mpc. In contrast to other studies, we do not detect a significant signal for these largerseparation filaments in any of the surveys individually. In combination we find a weak signal at 1.6 σ significance, with ℱ_{c} = 1.3 ± 0.6. The halfmaximum radius of the density profile, r_{c}, is, however, unconstrained.
5. Conclusions and discussion
In this paper, we present a 3.4σ detection of filamentary structure connecting luminous red galaxies separated by 3 − 5 h^{−1} Mpc. Through a series of null tests we have verified the robustness of this result. Our work extends the analysis presented in C16 by improving the methodology to null the weak lensing signal from the LRGs in order to isolate the weak lensing distortions induced by the filamentary structure alone. We note that this nulling method cannot distinguish between a pair of spherical haloes joined by a cylinder of matter, or two elliptical haloes which extend towards each other. Higuchi et al. (2014) shows that there is no hard line between a filament and its corresponding halo, with haloes typically extending along the filament. As we find a strong nulled signal on scales much larger than the typical (∼1 h^{−1} Mpc) virial radius of the LRGs (Parejko et al. 2013), we would argue that it is unlikely to originate from two perfectly aligned haloes (see Xia et al. 2017). But nevertheless we prefer to refer to our detection as that of filamentary structure, rather than that of a filament per se.
Previous studies have focused on LRGs separated by 6 − 10 h^{−1} Mpc. Using all three lensing surveys (KV450, CFHTLenS and RCSLenS), we do not detect a significant signal for these largerseparation filaments in either the surveys individually, or in combination. This is in contrast to the significant 5σ weak lensing detection of 6 − 10 h^{−1} Mpc separation filaments reported by Epps & Hudson (2017) using the same CFHTLenS dataset that has been analysed in this study. We report that we are unable to reproduce their result, even when adopting the same methodology. In comparison to C16, we find no comparable detection to their reported ∼4σ detection. As we have shown that sampling variance makes a significant contribution to the overall error budget, the increase by a factor of 5 in the number of filaments studied by C16, in contrast to this analysis, is key to their detection, even though the C16 lensing source galaxy density is significantly shallower than the source densities of KV450, CFHTLenS and RCSLenS. Neither Epps & Hudson (2017) nor C16 constrain the parameters of the Mead et al. (2010) filament model, but in the case of C16 we can compare the amplitudes of the measured signals. When adopting the C16 nulling approach, we find the two shear measurements to be fully consistent. The mean amplitude of our measurement is about four times larger than the amplitude reported in Kondo et al. (2019) (see discussion in Appendix A). Given our error budget, however, the results are consistent.
Authors de Graaff et al. (2019) calculate the average density between galaxy pairs separated by 6 − 14 h^{−1} Mpc using a CMB lensing convergence map, finding . This estimate assumes that the matter density follows a cylindrical filament model, with density
where ℓ defines the size of the filament in the lineofsight direction, r_{⊥} defines the distance perpendicular to the filament axis on the projected sky and σ is the intrinsic width. Integrating this density model over the lineofsight, we can relate this to the surface mass density at the centre of the filament ℱ_{c} as
For the value σ = 1.5 h^{−1} Mpc adopted by de Graaff et al. (2019), and which is our bestfit amplitude parameter ℱ_{c}, we find for the 3 − 5 h^{−1} Mpc filament sample at the average filament redshift z = 0.299. For the 6 − 10 h^{−1} Mpc filament sample we find , which is consistent with the de Graaff et al. (2019) result. Adapting the Mead et al. (2010) filament model, we can also integrate the model over the perpendicular distance and calculate the total mass enclosed between the two LRGs as
Taking the bestfit parameters ℱ_{c} and r_{c} for 3 − 5 h^{−1} Mpc measurement, we find M_{fil} = 4.9 ± 2.0 × 10^{13} h^{−1} M_{⊙}. It is worth noting that this estimate is based on the approximation that the deflection potential vanishes along the filament major axis. A more detailed analysis would attempt to obtain the excess mass map under the framework of galaxygalaxygalaxy lensing (Simon et al. 2008, 2019) but with a much larger separation of galaxy pairs.
Looking forward to upcoming deep weak lensing surveys, such as the European Space Agency’s Euclid mission^{5} and the Large Synoptic Survey Telescope (LSST^{6}), and deep spectroscopic surveys such as the Dark Energy Spectroscopic Instrument (DESI^{7}), the methodology that we present in this paper could be used to probe filamentary structure as a function of LRG mass and redshift. The combination of overlapping weak lensing surveys and spectroscopic surveys will provide the optimal datasets with which to fully explore the cosmic web.
We use a redshift cut z_{B} − z_{l} > 0.1 everywhere to ensure that the majority of our source galaxies are behind the foreground galaxies, and not associated with them. Fig. D3 in Amon et al. (2018) demonstrates that with this selection, contamination of the KiDS source sample is negligible for the scales we probe at > 0.1 h^{−1} Mpc.
Acknowledgments
We acknowledge support from the European Research Council under grant numbers 647112 (QX, CH, AA, MA, JHD, TT) and 770935 (HH, AW). CH also acknowledges the support from the Max Planck Society and the Alexander von Humboldt Foundation in the framework of the Max PlanckHumboldt Research Award endowed by the Federal Ministry of Education and Research. YC acknowledges the support of the Royal Society through the award of a University Research Fellowship and an Enhancement Award. We acknowledge support from the European Commission under a MarieSklodwoskaCurie European Fellowship under project numbers 656869 (JHD) and 797794 (TT). HH is supported by a Heisenberg grant of the Deutsche Forschungsgemeinschaft (Hi 1495/51). AK acknowledges support from Vici grant 639.043.512, financed by the Netherlands Organisation for Scientific Research (NWO). PS acknowledges support from the Deutsche Forschungsgemeinschaft in the framework of the TR33 “The Dark Universe”. Computations for the Nbody simulations were performed in part on the Orcinus supercomputer at the WestGrid HPC consortium (www.westgrid.ca), in part on the GPC supercomputer at the SciNet HPC Consortium. SciNet is funded by: the Canada Foundation for Innovation under the auspices of Compute Canada; the Government of Ontario; Ontario Research Fund – Research Excellence; and the University of Toronto. Funding for SDSSIII has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, and the U.S. Department of Energy Office of Science. The SDSSIII web site is http://www.sdss3.org/. We thank SDSSIII for making their data products so easily accessible. SDSSIII is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSSIII Collaboration including the University of Arizona, the Brazilian Participation Group, Brookhaven National Laboratory, Carnegie Mellon University, University of Florida, the French Participation Group, the German Participation Group, Harvard University, the Instituto de Astrofisica de Canarias, the Michigan State/Notre Dame/JINA Participation Group, Johns Hopkins University, Lawrence Berkeley National Laboratory, Max Planck Institute for Astrophysics, Max Planck Institute for Extraterrestrial Physics, New Mexico State University, New York University, Ohio State University, Pennsylvania State University, University of Portsmouth, Princeton University, the Spanish Participation Group, University of Tokyo, University of Utah, Vanderbilt University, University of Virginia, University of Washington, and Yale University. Author contributions: All authors contributed to the development and writing of this paper. The authorship list is given in two groups: the lead authors (QX, NR, CH), followed by an alphabetical group who contributed to either the scientific analysis or the data products.
References
 Alam, S., Albareti, F. D., Allende Prieto, C., et al. 2015, ApJS, 219, 12 [NASA ADS] [CrossRef] [Google Scholar]
 Amon, A., Heymans, C., Klaes, D., et al. 2018, MNRAS, 477, 4285 [NASA ADS] [CrossRef] [Google Scholar]
 AragónCalvo, M. A., van de Weygaert, R., & Jones, B. J. T. 2010, MNRAS, 408, 2163 [NASA ADS] [CrossRef] [Google Scholar]
 Bartelmann, M., & Schneider, P. 2001, Phys. Rep., 340, 291 [NASA ADS] [CrossRef] [Google Scholar]
 Benítez, N. 2000, ApJ, 536, 571 [NASA ADS] [CrossRef] [Google Scholar]
 Bond, J. R., Kofman, L., & Pogosyan, D. 1996, Nature, 380, 603 [NASA ADS] [CrossRef] [Google Scholar]
 Bregman, J. N. 2007, ARA&A, 45, 221 [NASA ADS] [CrossRef] [Google Scholar]
 Briel, U. G., & Henry, J. P. 1995, A&A, 302, L9 [NASA ADS] [Google Scholar]
 Cheng, H., & Gupta, K. C. 1989, ASME Trans. Ser. E: J. Appl. Mech., 56, 139 [NASA ADS] [CrossRef] [Google Scholar]
 Choi, A., Heymans, C., Blake, C., et al. 2016, MNRAS, 463, 3737 [NASA ADS] [CrossRef] [Google Scholar]
 Clampitt, J., Miyatake, H., Jain, B., & Takada, M. 2016, MNRAS, 457, 2391 [Google Scholar]
 Colberg, J. M., Krughoff, K. S., & Connolly, A. J. 2005, MNRAS, 359, 272 [NASA ADS] [CrossRef] [Google Scholar]
 Colless, M., Dalton, G., Maddox, S., et al. 2001, MNRAS, 328, 1039 [NASA ADS] [CrossRef] [Google Scholar]
 de Graaff, A., Cai, Y.C., Heymans, C., & Peacock, J. A. 2019, A&A, 624, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dietrich, J. P., Werner, N., Clowe, D., et al. 2012, Nature, 487, 202 [Google Scholar]
 Dolag, K., Meneghetti, M., Moscardini, L., Rasia, E., & Bonaldi, A. 2006, MNRAS, 370, 656 [NASA ADS] [Google Scholar]
 Dvornik, A., Hoekstra, H., Kuijken, K., et al. 2018, MNRAS, 479, 1240 [NASA ADS] [CrossRef] [Google Scholar]
 Eardley, E., Peacock, J. A., McNaughtRoberts, T., et al. 2015, MNRAS, 448, 3665 [NASA ADS] [CrossRef] [Google Scholar]
 Ebeling, H., Barrett, E., & Donovan, D. 2004, ApJ, 609, L49 [Google Scholar]
 Epps, S. D., & Hudson, M. J. 2017, MNRAS, 468, 2605 [NASA ADS] [CrossRef] [Google Scholar]
 HarnoisDéraps, J., & van Waerbeke, L. 2015, MNRAS, 450, 2857 [NASA ADS] [CrossRef] [Google Scholar]
 HarnoisDéraps, J., Amon, A., Choi, A., et al. 2018, MNRAS, 481, 1337 [NASA ADS] [CrossRef] [Google Scholar]
 Heymans, C., Van Waerbeke, L., Bacon, D., et al. 2006, MNRAS, 368, 1323 [NASA ADS] [CrossRef] [Google Scholar]
 Heymans, C., Van Waerbeke, L., Miller, L., et al. 2012, MNRAS, 427, 146 [NASA ADS] [CrossRef] [Google Scholar]
 Higuchi, Y., Oguri, M., & Shirasaki, M. 2014, MNRAS, 441, 745 [NASA ADS] [CrossRef] [Google Scholar]
 Hildebrandt, H., Erben, T., Kuijken, K., et al. 2012, MNRAS, 421, 2355 [NASA ADS] [CrossRef] [Google Scholar]
 Hildebrandt, H., Choi, A., Heymans, C., et al. 2016, MNRAS, 463, 635 [NASA ADS] [CrossRef] [Google Scholar]
 Hildebrandt, H., Viola, M., Heymans, C., et al. 2017, MNRAS, 465, 1454 [Google Scholar]
 Hildebrandt, H., Köhlinger, F., van den Busch, J. L., et al. 2018, A&A, accepted, https://doi.org/10.1051/00046361/201834878 [Google Scholar]
 Kondo, H., Miyatake, H., Shirasaki, M., Sugiyama, N., & Nishizawa, A. J. 2019, ArXiv eprints [arXiv:1905.08991] [Google Scholar]
 Kull, A., & Böhringer, H. 1999, A&A, 341, 23 [NASA ADS] [Google Scholar]
 Macciò, A. V., Dutton, A. A., van den Bosch, F. C., et al. 2007, MNRAS, 378, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Mandelbaum, R., Hirata, C. M., Seljak, U., et al. 2005, MNRAS, 361, 1287 [NASA ADS] [CrossRef] [Google Scholar]
 Mead, J. M. G., King, L. J., & McCarthy, I. G. 2010, MNRAS, 401, 2257 [NASA ADS] [CrossRef] [Google Scholar]
 Miller, L., Heymans, C., Kitching, T. D., et al. 2013, MNRAS, 429, 2858 [Google Scholar]
 Parejko, J. K., Sunayama, T., Padmanabhan, N., et al. 2013, MNRAS, 429, 98 [NASA ADS] [CrossRef] [Google Scholar]
 Pimbblet, K. A., & Drinkwater, M. J. 2004, MNRAS, 347, 137 [NASA ADS] [CrossRef] [Google Scholar]
 Ross, A. J., Percival, W. J., Sánchez, A. G., et al. 2012, MNRAS, 424, 564 [NASA ADS] [CrossRef] [Google Scholar]
 Schneider, P., & Watts, P. 2005, A&A, 432, 783 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Simon, P., Watts, P., Schneider, P., et al. 2008, A&A, 479, 655 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Simon, P., Saghiha, H., Hilbert, S., et al. 2019, A&A, 622, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Springel, V., White, S. D. M., Jenkins, A., et al. 2005, Nature, 435, 629 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Tanimura, H., Hinshaw, G., McCarthy, I. G., et al. 2019, MNRAS, 483, 223 [NASA ADS] [CrossRef] [Google Scholar]
 Velander, M., Kuijken, K., & Schrabback, T. 2011, MNRAS, 412, 2665 [NASA ADS] [CrossRef] [Google Scholar]
 Velander, M., van Uitert, E., Hoekstra, H., et al. 2014, MNRAS, 437, 2111 [NASA ADS] [CrossRef] [Google Scholar]
 Werner, N., Finoguenov, A., Kaastra, J. S., et al. 2008, A&A, 482, L29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Wilks, S. S. 1938, Ann. Math. Stat., 9, 60 [Google Scholar]
 Williams, D. 2001, Weighing the Odds: A Course in Probability and Statistics (Cambridge: Cambridge University Press) [CrossRef] [Google Scholar]
 Wright, A. H., Hildebrandt, H., Kuijken, K., et al. 2019, A&A, 632, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Wright, C. O., & Brainerd, T. G. 2000, ApJ, 534, 34 [NASA ADS] [CrossRef] [Google Scholar]
 Xia, Q., Kang, X., Wang, P., et al. 2017, ApJ, 848, 22 [NASA ADS] [CrossRef] [Google Scholar]
 Zehavi, I., Zheng, Z., Weinberg, D. H., et al. 2011, ApJ, 736, 59 [NASA ADS] [CrossRef] [Google Scholar]
 Zel’dovich, Y. B. 1970, A&A, 5, 84 [NASA ADS] [Google Scholar]
Appendix A: Remarks on the nulling technique
In this appendix, we review the C16 nulling technique and develop an improved methodology to isolate the weak lensing signal from filaments. To explain the motivation behind nulling, we start with a single circularly symmetric halo positioned at the origin (0,0), for which the complex shear is given by
where is the mean convergence inside r. We define its counterpart γ^{c} as the complex shear at the same radial position, with a clockwise rotation of 90°, such that
The counterpart is therefore able to “null” the shear, as γ^{c}(r)+γ(r) = 0.
For a twohalo system as shown in Fig. A.1, the shear at each position (r, θ) is composed with the shear from halo h_{1} and the shear from halo h_{2}. We can write this as
Fig. A.1. The bipolar configuration for two haloes centred symmetrically about origin O at h_{1} and h_{2}. 
where the coordinates (t, φ) are defined with halo h_{1} at the origin, and the coordinates (s, ϕ) are defined with halo h_{2} at the origin. Starting from position P_{0} = (r_{0}, θ_{0}), shear is given by
A clockwise rotation around halo h_{1} by 90° takes us to position P_{1} at (r_{1}, θ_{1}). The shear here is given by
At position P_{1}, we see the shear contribution from halo h_{1} is the counterpart to the shear contribution from halo h_{1} at position P_{0}. We next rotate around halo h_{2} by 90° to position P_{2} = (r_{2}, θ_{2}), where the shear is given by
Similarly another 90° rotation about halo h_{1} (see Fig. A.2) to position P_{3} gives
Fig. A.2. Illustration of the nulling technique. Every two line segments with the same colour represents an anticlockwise rotation with respect to one halo. Right: configuration described in C16. Left: our adopted “X19” configuration which starts below the horizontal axis. 
We note that, after another 90° rotation about halo h_{2}, we come to position P_{4} which is our starting point P_{0}. The sum of the shear from position P_{0}, P_{1}, P_{2} and P_{3} is given by
If we now add in a filament shear γ_{f} such that at each position γ = γ_{f} + γ_{h1} + γ_{h2}, then
As we expect the filament shear profile to be symmetric about the horizontal axis, we also sum over the shear from positions and P_{3}′ that are reflections of P_{0}, P_{1}, P_{2} and P_{3} about the horizontal axis respectively, in order to get the average shear value at any distance away from filament axis. We therefore define the nulling operator 𝒩 as
It is interesting to note that the above equations also apply when two halos are of different masses given their circular symmetry.
In Fig. A.2 we show two configurations for nulling. C16 chose to start P_{0} above the horizontal axis and our work starts P_{0} below the horizontal axis. We investigate the difference between these two approaches by first constructing a noiseless shear map from two NFW halo profiles using Eq. (14) in Wright & Brainerd (2000) with M_{vir} = 10^{13} h^{−1} M_{⊙} and assuming a mass–concentration relation from Macciò et al. (2007). We assume both haloes are located at z = 0.3 with background sources at redshift 0.7, close to the mean value of KiDS. The resultant γ_{1} map is shown in the topleft panel of Fig. A.3.
Fig. A.3. Left panel: shear maps from two NFW halos with or without fiducial filaments. Top row: γ_{1} map generated by two NFW haloes only, while the middle row shows the γ_{1} map with the addition of a fiducial filament profile. Lowest row: γ_{2} map generated by two NFW haloes. it Right panel: results of the nulling procedure corresponding to the shear map on the left. In the top row, we see under both the X19 configuration (solid line) and C16’s (dashed line), the resulting signal is consistent with zero. When adding a fiducial filament profile, in the middle row, we see our X19 configuration correctly recovers the input value whereas the C16 configuration is biased on large scales. In the lower row, we verify that both configurations null the γ_{2} signal from two NFW haloes. 
For each pixel (x_{1}, x_{2}) on the map, we calculate 𝒩[γ_{i}(x_{1}, x_{2})], and show in the topright panel in Fig. A.3, the average of the sum along the horizontal axis, i.e., . In the lowerright panel in Fig. A.3, we show the resulting 𝒩[γ_{2}(x_{1}, x_{2})]. As expected, the sum of nulling pixels are zeros (note the 10^{−16}/10^{−18} on the yaxis). In the middle panel we repeat the analysis with the inclusion of a fiducial filament modelled using the powerlaw profile model in Eq. (5). The filament contribution to the shear is included in the 2D shear map as a powerlaw, symmetric about the x_{1}axis (at x_{2} = 0) with k_{c} = 0.02 and r_{c} = 0.4 h^{−1} Mpc. These parameters were chosen to roughly replicate previous studies (Dolag et al. 2006), in order to model the true signal contrast in observations. If the nulling method is precise, we would expect the nulled filament signal measured from the map, to recover the input shear from the filament, , irrespective of the nulling method used. The middleright panel shows that this is indeed the case, with our nulling method. The nulling method in C16, however, recovers a lower expectation value for the shear induced by the filament. This difference results from the C16 configuration. As shown in Eq. (A.10), for the C16 configuration, since P_{0} and P_{1} are both in the filament region, the nulling operator effectively mixes scales. This has an effect in producing a significant signal on large scales which does not reflect the underlying filament density profile. As P_{1}, P_{2} and P_{3} in our adopted configuration lies outside the bridge between two haloes, γ_{f}(r_{1}, θ_{1})+γ_{f}(r_{2}, θ_{2})+γ_{f}(r_{3}, θ_{3}) is negligible. This enables us to recover the density profile accurately. We note that Kondo et al. (2019) adopted the C16 estimator but in their Eq. (9) (the equivalent of our Eq. (A.11)), they included an additional factor of 4 in the denominator, which, in our testcase in Fig. A.3 would result in the underestimation of the filament signal by a factor of 4.
Appendix B: Remarks on the spherical rotation
Here we detail the spherical rotation used in Sect. 3.2 to project all filaments onto the same reference frame. As illustrated in Fig. B.1, we rotate all galaxies about a given axis so that the filament (the shorter arc connected by solid pink diamonds) is transformed to lie horizontally on the equator (hollow pink diamonds). To do this, we first transfer the right ascension and declination onto a 3D vector on a unit sphere, such that their positions are g_{1} and g_{2}. The normal vector is defined as . Noting the rotation axis lies on the equator and is perpendicular to the normal vector, we can write down the rotation axis and angle β using components of , so that
Fig. B.1. Illustration of the spherical rotation for filament pairs (solid pink diamonds). Open symbols are the rotated positions. 
We note that by defining g_{1} to be always on the left of g_{2}, there is no ambiguity in the definition of the rotation axis, and the rotation angle will always lie between . Rodrigues’ rotation formula (Cheng & Gupta 1989) then allows us to rotate every point on the sky to the desired frame where the filament pair now lies on the equator, such that for each source galaxy at position g_{s}, the new position is at:
It is worth noting that because the shear was measured in each galaxy’s local (RA, Dec) coordinate frame, the angle of rotation is different for different source galaxies. The transformed shear map () for each galaxy is thus given by
where for each source galaxy, the rotation angle ϕ_{s} is defined by the angular change of its local coordinate frame, for example, the declination to the yaxis of the filament.
All Tables
χ^{2} value and pvalue for all computed signals from each individual survey and their combination.
All Figures
Fig. 1. The effective n(z) of all three lensing surveys, as defined in Eq. (12), when using LOWZ galaxies as lenses. Each curve is normalised such that ∫n(z) dz = 1. 

In the text 
Fig. 2. Weighted shear measured between LRG pairs in the SLICS simulation of a random 400 deg^{2} degree survey. Left: results from 3 to 5 h^{−1} Mpc filament candidates. Upper panel: average surface mass density ℱ_{γ1}, and the lower panel shows the measured crossshear component. On both panels, the blue data shows the result for a KiDSlike survey depth and shape noise, and the orange data points show the measurement for a noisefree simulation with the errorbar given by the error on the mean of all 158 realisations. Right: set of results from 6 to 10 h^{−1} Mpc filament candidates. 

In the text 
Fig. 3. Redshift evolution of the filament signal in the noisefree SLICS simulation. Left panel: ℱ_{γ1} signal measured from filament candidates in SLICS for four redshift quantiles. Right panel: 68% and 95% confidence region of the model parameters (ℱ_{c}, r_{c}) from the corresponding signal in the left panel, with the result from all samples combined shown in grey. 

In the text 
Fig. 4. The detection of the cosmic web between neighbouring luminous red galaxies as detected through the weak lensing of background galaxies from different lensing surveys. Left: xaxis is the distance measured perpendicular to the filament axis scaled such that 1 is equivalent to the projected separation between the pair of LRGs. The yaxis is the nulled shear signal where ℱ_{γ1} (blue data points) measures the average surface mass density of the filament, and ℱ_{γ2} is expected to be consistent with zero and hence serves as a null test. Lower right small panel: measurements from the threesurveys combined using inverse variance weighting. We note that this additional panel is purely for illustration, however, as our joint surveyconstraints on the filament model are derived from a combination of the surveys on the likelihoodlevel. Right: estimated parameters in the filament density model Eq. (5) from the stacked signal for all surveys individually and their combination. Upper: results from 3 ∼ 5 h^{−1} Mpc Physical Pairs, our filament candidates; Lower: results from 3 ∼ 5 h^{−1} Mpc nonphysical Pairs, our control sample. 

In the text 
Fig. A.1. The bipolar configuration for two haloes centred symmetrically about origin O at h_{1} and h_{2}. 

In the text 
Fig. A.2. Illustration of the nulling technique. Every two line segments with the same colour represents an anticlockwise rotation with respect to one halo. Right: configuration described in C16. Left: our adopted “X19” configuration which starts below the horizontal axis. 

In the text 
Fig. A.3. Left panel: shear maps from two NFW halos with or without fiducial filaments. Top row: γ_{1} map generated by two NFW haloes only, while the middle row shows the γ_{1} map with the addition of a fiducial filament profile. Lowest row: γ_{2} map generated by two NFW haloes. it Right panel: results of the nulling procedure corresponding to the shear map on the left. In the top row, we see under both the X19 configuration (solid line) and C16’s (dashed line), the resulting signal is consistent with zero. When adding a fiducial filament profile, in the middle row, we see our X19 configuration correctly recovers the input value whereas the C16 configuration is biased on large scales. In the lower row, we verify that both configurations null the γ_{2} signal from two NFW haloes. 

In the text 
Fig. B.1. Illustration of the spherical rotation for filament pairs (solid pink diamonds). Open symbols are the rotated positions. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.