Jérôme Gâteau

Ferdinand Claude

Gilles Tessier

Marc Guillon
email: marc.guillon@parisdescartes.fr

Jérôme Gateau

Topological transformations of speckles

come

Introduction

The propagation of coherent light through scattering media yields random wavefields with typical intensity structures called optical speckles. The control of light distribution inside and through complex media by wavefront modulation of the impinging beam is of critical importance for application ranging from bio-imaging [START_REF] Horstmeyer | Guidestarassisted wavefront-shaping methods for focusing light into biological tissue[END_REF] to telecommunications [START_REF] Rodenburg | Influence of atmospheric turbulence on states of light carrying orbital angular momentum[END_REF], for instance by multiplexing information with orbital angular momentum [START_REF] Gong | Optical orbital-angularmomentum-multiplexed data transmission under high scattering[END_REF]. Information transmission through diffusers is typically characterized in terms of field and intensity correlations [START_REF] Akkermans | Mesoscopic Physics of Electrons and Photons[END_REF]. For diffusers exhibiting so-called ''memory effect"correlations, important invariants were identified under specific spatial (tilt and shift) transformations [START_REF] Feng | Correlations and fluctuations of coherent wave transmission through disordered media[END_REF][START_REF] Freund | Memory effects in propagation of optical waves through disordered media[END_REF][START_REF] Judkewitz | Translation correlations in anisotropically scattering media[END_REF][START_REF] Osnabrugge | Generalized optical memory effect[END_REF]. Additionally, regardless of the wavefront of the impinging beam, critical points in random wavefields exhibit many topological correlations [START_REF] Freund | 1001 correlations in random wave fields[END_REF], which thus demand the development specific tools to be analyzed. Optical vortices are especially important critical points since they are centered on singular phase points coinciding with nodal points of the intensity. They spontaneously appear in random wavefields [START_REF] Nye | Dislocations in Wave Trains[END_REF], and thereby allow efficient super-resolution microscopy [START_REF] Pascucci | Superresolution Imaging of Optical Vortices in a Speckle Pattern[END_REF][START_REF] Pascucci | Compressive three-dimensional super-resolution microscopy with speckle-saturated fluorescence excitation[END_REF]. The present work aims at exploring the possibility to manipulate topological correlations between critical points in random wavefields under symmetry control and spiral wavefront modulation, in a Fourier plane of the impinging beam.

Critical points are characterized by their topological charge and their Poincaré number [START_REF] Dennis | Topical singularities in wave fields[END_REF]. They may typically be controlled by applying phase or amplitude masks in a Fourier plane. Any smooth and regular transform of the wavefield (either in phase or amplitude) induces changes preserving both the topological charge and the Poincaré number [START_REF] Freund | 1001 correlations in random wave fields[END_REF][START_REF] Nye | Phase saddles and dislocations in two-dimensional waves such as the tides[END_REF]. Noteworthy, these conservation rules account for the topological stability of isolated vortices of charge 1 in speckles since the creation or annihilation of vortices can only occur by pairs [START_REF] Nye | Phase saddles and dislocations in two-dimensional waves such as the tides[END_REF][START_REF] Pascucci | Superresolution Imaging of Optical Vortices in a Speckle Pattern[END_REF]. As opposed to smooth phase transforms, the addition of a spiral phase mask in a Fourier plane is a singular transform and results in a change of the total orbital angular momentum [START_REF] Larkin | Natural demodulation of two-dimensional fringe patterns. i. general background of the spiral phase quadrature transform[END_REF][START_REF] Dennis | Chapter 5 singular optics: Optical vortices and polarization singularities[END_REF]. Recently, considering correlations between the spatial distribution of critical points in a speckle under such spiral phase transforms [START_REF] Gateau | Complementary speckle patterns: Deterministic interchange of intrinsic vortices and maxima through scattering media[END_REF], we observed a strong inter-play between intensity maxima and optical vortices. More precisely, the obtained results suggested that the topological charge of these critical points were all incremented by applying a +1 spiral phase mask in the Fourier plane. The impossibility to spontaneously get +2-charged vortices (unstable and thus unlikely in random light structures [START_REF] Freund | Optical dislocation networks in highly random media[END_REF]) resulted in the observation of a partial cyclic permutation of the three populations of critical points (namely, maxima and ±1-charged vortices). Furthermore, as a third kind of possible transform, it was observed that the orbital angular momentum may be not conserved when using amplitude masks with 1 a high degree of symmetries [START_REF] Ferrando | Vortex transmutation[END_REF][START_REF] Gao | Experimental demonstration of free-space optical vortex transmutation with polygonal lenses[END_REF]. As a result, optical vortices can be created using simple amplitude masks [START_REF] Schoonover | Creating polarization singularities with an n-pinhole interferometer[END_REF][START_REF] Li | Generation of second-order vortex arrays with six-pinhole interferometers under plane wave illumination[END_REF]. This property proved to be of interest for imaging applications to reveal symmetries of an imaged object [START_REF] Yang | Digital spiral object identification using random light[END_REF][START_REF] Xie | Using a complex optical orbital-angular-momentum spectrum to measure object parameters[END_REF][START_REF] Zhang | Encoding and decoding of digital spiral imaging based on bidirectional transformation of light's spatial eigenmodes[END_REF] and for allowing topological charge measurements [START_REF] Hickmann | Unveiling a truncated optical lattice associated with a triangular aperture using light's orbital angular momentum[END_REF], especially in astronomy [START_REF] Berkhout | Method for probing the orbital angular momentum of optical vortices in electromagnetic waves from astronomical objects[END_REF].

Here, combining spiral phase transforms of order n with star-like amplitude masks having discrete point group symmetries D 3 and D 4 , we study experimentally the topological correlations between intensity maxima and optical vortices in speckles. A new co-localization criterion is proposed, inspired by statistical mechanics. Although random wavefields do not possess any symmetry, such a combination allows us to strengthen periodicity and even to control the period of the cyclic permutation. Noteworthy, for an amplitude mask of symmetry D 4 , a phase saddle point appears as a complementary critical point to complete a cycle of period 4. A transposition between vortices of charge -1 and vortices of charge +1 is also revealed when adding a 2-charged spiral phase mask.

Experimental procedure

The experimental procedure consisted in modulating a random phase pattern in a Fourier plane with an amplitude mask and a spiral phase mask. Here, spiral phase masks of order n, SP n (θ) = e i.n.θ (in polar coordinates), were applied for n ∈ -6; 6 . As amplitude masks, three Binary Amplitude (BA) masks were used: a disk and two periodic angular slits with a point group symmetry D 3 and D 4 . They are defined by the following angular transmission function (in polar coordinates):

BA N (θ) = 1, if |θ -(k -1 2 ). 2.π N | < π 32 with k ∈ 1; N . 0, otherwise.
(1) for N ∈ {3, 4}. By convention, BA ∞ defines the disk-shaped aperture (obtained for N 32). For N < 32, the orbital angular momentum content (or spiral spectrum [START_REF] Torner | Digital spiral imaging[END_REF]) of the BA N aperture exhibits discrete harmonics for the spiral modes n = p.N with p ∈ Z. Given the width of the angular slits and provided that N = 3 or 4, the aperture BA N can be considered as invariant by the addition of SP n when n = ±N or ±2.N. Such an invariance in the Fourier plane is thus necessarily associated with a periodic transform of the speckle pattern in the real space.

The experimental configuration is detailed in Fig. 1 [See Supplement 1, Section 1 for further details on the experimental methods]. A spatial light modulator (SLM) (LCOS, X10468, Hamamatsu) was illuminated with a collimated laser beam at 635 nm and Fourier conjugated to a camera (768x1024 pixels, pixel size: 4.65x4.65 µm 2 ) with a converging lens. The phase Φ n and amplitude A n of the modulated (SP n mask) random wave were measured at the camera plane by phase stepping interferometry [START_REF] Creath | Interferometry -phasemeasurement interferometry[END_REF]. To do so, the SLM (792x600 SLM pixels, pixel size: 20x20 µm 2 ) was split in two parts to generate both the modulated random wave (or signal wave) on one side and a reference wave on the other side (see Fig. 1a). The signal wave was generated by adding simultaneously the scattering random phase pattern, the spiral phase modulation SP n and the amplitude mask BA N . Adding a blazed grating achieved spatial separation of the imprinted signal wavefront from undiffracted light (the latter being sent to a beam-block). The signal speckle intensity I n could be measured directly by removing the contribution of the reference beam.

For phase-stepping interferometry, an additional Fresnel lens was added to the reference beam in order to cover the camera surface. The latter spherical contribution as well as the relative phase-tilt between the signal and the reference beams were removed in a numerical post-processing step. A stack of eight images was sequentially recorded by phase shifting the reference beam by 2π/8 phase-steps. All BA masks had the same radius of r = 170 pixels at the SLM, so yielding the same speckle grains size on the camera plane: λ/(2NA) = 70 µm -Full Width Half Maximum (FWHM), where λ is the wavelength and NA r/f 4.53 × 10 -3 the numerical aperture of illumination (with f = 750 mm the focal length of the lens L in Fig. 1a). The speckle grain size thus covered 15 camera pixels and ensured a fine sampling of the speckle patterns. Hereafter, all distances and spatial densities are expressed setting λ/(2NA) as the length unit.

In Fig. 1(b-d), an illustration of the speckle intensity and phase maps obtained for the three different geometries of BA masks is shown. In the following study, intensity maps I n and phase maps Φ n were measured for all three BA N masks (N ∈ {3, 4, ∞}) and for each SP n masks (n ∈ -6; 6 ). For comparison, the intensity-map I rand and the phase-map Φ rand obtained for a non-correlated scattering random pattern were acquired for all the BA N masks independently.

3 Statistical analysis of the topological correlation between critical points

Studied critical points

The field at the camera being linearly polarized, optical fields are here studied as scalar fields. The location of the main critical points of the experimental intensity and phase maps were measured [See Supplement 1, Section 1.C for details on the detection of critical points] and their statistical correlation distances were analyzed. Importantly, for a given BA mask, adding SP n masks preserves all statistical properties of the speckle patterns, such as the number-density of critical points. Phase saddle-points of Φ n are notated S p n and vortices of charge ±1: V ± n . Vortices of charge higher than 1 do not appear in Gaussian random wavefields [START_REF] Freund | Optical dislocation networks in highly random media[END_REF]. Maxima and saddle-points of I n are notated M n and S I n , respectively. Non-zero minima and phase extrema have not been considered here, since having significantly lower densities [START_REF] Freund | Amplitude topological singularities in random electromagnetic wavefields[END_REF]. All the notations are summarized in Table 1.

The measured average number-densities of the critical points are presented in Table 2. The density of the critical points of type X (X = V ± , M , S p or S I ) is notated ρ(X). As expected, ρ(V -) and ρ(V + ) are equal [START_REF] Freund | 1001 correlations in random wave fields[END_REF], and ρ(X) depends both on the type of critical point and the BA mask. To study statistical transformations of critical points quantitatively, new specific tools are presented. What we discuss as transformation of critical points by the addition of spiral phase masks refers to the mean nearest neighbor distances between populations of critical points and calls for a discrimination parameter. Two specific statistical tools are then described below: the radial density function (RDF) and the Weighted Median Normalized Distance (WMND).

WMD WMND

rand. a) b) c) d) = BA ∞ BA ∞ BA ∞ BA ∞
In our previous study [START_REF] Gateau | Complementary speckle patterns: Deterministic interchange of intrinsic vortices and maxima through scattering media[END_REF], correlations between critical points could be characterized by computing the radial probability density function (RPDF) of the nearest-neighbor distance. Fig. 2a presents RPDFs of the distance d V - 0 , X in the case of BA ∞ . We define d(Y, X) as the distance between a Y -point and the closest X-point. The RPDF(r) corresponds to the probability to find the closest X-point at the distance r from a Y -point, per unit area ( ∞ 0 RP DF (r)2πrdr = 1). One drawback associated with the use of the RPDF is that it may suggest paradoxes if improperly interpreted. Considering d V - 0 , V + 2 and d V - 0 , S I rand , it seems that V - 0 correlates both with V + 2 and S I rand since both RPDFs reach high values at zero distances.

While a correlation is expected in the former case (due to topological charge incrementation), no correlation is expected from the latter which involves two independent sets of random points. The reason why the amplitude of the RPDF of d V - 0 , S I rand is higher than the one of d V - 0 , V + 2 at zero distances, is just due to the ∼ 3-times higher spatial density of intensity saddle points S I rand as compared to vortices V + 2 (see Table 2): the probability to find a saddle point at close distance is thus larger. To quantitatively characterize topological correlations, we thus need to normalize RPDFs by the number densities ρ(X).

Our first statistical tool, the radial-distribution function (RDF) -well known in statistical mechanics [START_REF] Chandler | Introduction to Modern Statistical Mechanics[END_REF] -was extended here for nearest neighbor by normalizing the RPDF of d V - 0 , X by ρ(X), and the distances d V - 0 , X by the mean X-interpoint half-distance 2 ρ(X) -1 [START_REF] Tournus | Random nanoparticle deposition: Interparticle distances in 2D, 3D, and multilayer samples[END_REF]. Fig. 2c shows the RDF of the same data as in Fig. 2a. As a result, all the RDFs of d V - 0 , X rand .2. ρ(X) are superimposed for every X rand , and the spatial correlation betweeen V - 0 and V + 2 clearly appears.

To obtain a single binary parameter discriminating the spatial correlation between V - 0 and X, we further define the Weighted Median Normalized Distance (WMND) as a second statistical tool: the WMND V - 0 , X is the 50% weighted percentile of d V - 0 , X .2. ρ(X) with weights corresponding to the RDF values. A WMND V - 0 , X around 0.5 means that no spatial correlation exist between V - 0 and X, while WMND < 0.5 and WMND > 0.5 mean an attraction and a repulsion, respectively. A zero WMND value means perfect correlation while WMND=1 means perfect anti-correlation.

Fig 2d presents the WMND V - 0 , X for all the critical points considered in this study and for n ∈ 0; 6 . For comparison, the Weighted Median Distance (WMD) associated with the RPDFs -defined as the 50% weighted percentile of d V - 0 , X with weights corresponding to the RDPF values -is also computed and displayed in Fig. 2b. To validate this tool, taking BA ∞ as an illustrative example, we notice that WMND V - 0 , X rand is around 0.5 for all the X rand , as expected. By comparison, WMD V - 0 , S I rand = 0.35, which irrelevantly suggests correlations as discussed above. Moreover, for n > 3, the RDFs of d V - 0 , X n are observed to match the RDFs of d V - 0 , X rand : no noticeable spatial correlation is obtained for n > 3. As expected again, the WMND V - 0 , X rand are around 0.5 for n > 3. Conversely, we get WMD V - 0 , S I n < 0.38, which would falsely suggest correlations. All these observations validate the WMND as a parameter to assess the spatial correlation between pairs of critical points in a speckle pattern.

Topological correlations between critical points for the different amplitude masks

Fig. 3 presents the WMND(Y 0 , X n ) for all the critical points (Y 0 and X n ) screened in this study, for SP n masks with n ∈ -6; 6 and for the three considered BA apertures. The WMND was verified to be around 0.5 for all the amplitude masks and all the pairs (Y 0 ,X rand ), for which there is obviously no spatial correlation. For the sake of readability, in the following, we only discuss the interplay between critical points when adding positively charged SP masks but symmetrical behaviours are observed for negatively charged SP masks (Fig. 3). For the aperture BA ∞ , the WMND reveals several notice-able features (reported in Table 3). First, as expected from our previous study [START_REF] Gateau | Complementary speckle patterns: Deterministic interchange of intrinsic vortices and maxima through scattering media[END_REF], we notice some spatial correlations for the triplets (

(V - 0 ,X) (M 0 ,X) (V 0 + ,X) (S p 0 ,X) (S I 0 ,X) V - n M n V +
V - m-1 , M m , V + m+1 ). Because ρ(V - 0 ) < ρ(M 1
), a one-to-one transformation is impossible between vortices and maxima. Although V - 0 and V + 2 have the same number density, we also notice that WMND(V - 0 , V + 2 ) is significantly different from 0, indicating that the rate of the macroscopic transformation from V - 0 to V + 2 is below 1. In agreement with the 3-point cyclic permutation algebra observed in [START_REF] Gateau | Complementary speckle patterns: Deterministic interchange of intrinsic vortices and maxima through scattering media[END_REF], a weak attraction is found for the pair (V + 0 , V - +1 ), corresponding to the topological-charge equation 1+1 = -1. However, as an alternative transformation for V + 0 , a similar attraction is now also observed for (V + 0 , S p +1 ). V + 0 is thus subjet to a bifurcation between V - +1 and S p +1 , which implies two different mechanisms.

As a first possibility, some V + 0 transform into S p +1 . This transformation inspires the following interpretation: When adding a SP +1 phase mask to an isolated Laguerre-Gaussian beam with topological charge +1, a +2-charged vortex is obtained. Under weak perturbation, this +2 vortex splits into two +1 vortices, accompanied by the creation of both an intensity saddle point and a phase saddle point in between [START_REF] Nye | Phase saddles and dislocations in two-dimensional waves such as the tides[END_REF]. The creation of this pair of saddle points is governed by the Poincaré number conservation. In the frame of this model, a V + 0 vortex is expected to co-localize with both S I +1 and S p +1 and to anti-correlate with the two V + 1 that split away. Although no noticable spatial attraction was found for (V + 0 , S I +1 ), a co-localization is observed for (V + 0 , S p +1 ) and a weak repulsion is observed for the pair (V + 0 , V + +1 ), consistent with this interpretation. In speckles, where +2-charged vortices cannot be encountered since unstable [START_REF] Freund | Optical dislocation networks in highly random media[END_REF], the weak perturbation approximation cannot be fully valid, potentially accounting for the remaining discrepancy between experimental observations and the proposed model.

As a second possible transformation, the more surprising attraction of the pairs (V + 0 , V - +1 ) is observed, which calls for another mechanism. Since no such transformation can be imagined from isotropic V + 0 , it may only be interpreted by a mechanism dominated by strong perturbations. The statistically uniform mesh created by vortices and maxima in speckles [START_REF] Longuet-Higgins | Reflection and refraction at a random moving surface. i. pattern and paths of specular points[END_REF], together with strong correlations observed for pairs (M 0 , V + +1 ) and (V - 0 , M 1 ) seem to constrain V + 0 to colocalize with V - +1 . This transformation would deserve further analytical investigation but we anticipate that the creation mechanism of V - +1 from V + 0 can only be a many-body problem, involving the field structure (maxima, phase saddles and vortices) surrounding the initial V + 0 of interest. When adding a SP +2 mask for BA ∞ , V + 0 is not observed to significantly co-localize with any remarkable critical point (see Fig. 3 and Table 3), whereas two possible transformations might have been expected for V + 0 . On the one hand, from the 3-point cyclic permutation, we could expect that V + 0 would transform into M +2 . On the other hand, since in Table 3, maxima and phase saddle-points are noted to be simply exchanged (see pairs (M 0 , S p +2 ) and (S p 0 , M +2 )), a similar symmetrical echange between -1 and +1 vortices could be expected, yielding a transformation of V + 0 into V - 2 (as

V - 0 is transformed into V + 2 )
. However, no such correlation is observed either for the pair (V + 0 ,M 2 ) or for (V + 0 ,V - 2 ). Conversely, these correlations appear when applying amplitude masks BA 3 and BA 4 , respectively, as detailed in the following.

For |n| > 3, no significant spatial correlation with the addition of SP n is observed for BA ∞ . This aperture has a circular symmetry. Therefore, its spiral spectrum contains only the fondamental spiral mode n = 0, and is not invariant by the addition of any SP masks. All the described topological correlations associated with BA ∞ are summarized in Table 3. Table 3: Macroscopic transformations observed for the critical points Y 0 with the amplitude mask BA ∞ . The transformation rates are below 1.

Critical point Y 0 Adding SP +1 Adding SP +2 V - 0 → M 1 → V + 2 M 0 → V + 1 → S p 2 V + 0 → V - 1 + S p 1 → ∅ S p 0 → V - 1 → M 2 Adding SP -2 Adding SP -1 Critical point Y 0 ∅ ← V + -1 + S p -1 ← V - 0 S p -2 ← V - -1 ← M 0 V - -2 ← M -1 ← V + 0 M -2 ← V + -1 ← S p 0
As a possible solution to strengthen the 3-point cyclic permutation, we used the BA 3 amplitude mask, making the Fourier plane almost invariant with respect to the addition of SP ±3k (so long as 3k, with k integer, remains small enough).

Here, four main observations can be noted. First, as expected, the pairs (Y 0 , Y ±3 ) and (Y 0 , Y ±6 ) are observed to have a WMND very close to zero, indicating that the macroscopic transformation rate is close to 1 for all these pairs. Second, we observe that the cycle of period 3 reinforces the spatial correlations of the triplet (V - m-1 , M m , V + m+1 ) and even extends it to the 3 rd and 6 th spiral harmonics. Third, no noticeable correlation is observed between V + 0 and phase saddle-points S p +1 (although an anti-correlation is obtained between V + 0 and V + +1 ), contrary to the case of BA ∞ . The periodicity of 3 induces a strong correlation for the pairs (V + 0 , V - +1 ), and establishes a cyclic permutation of three populations of critical points V -, M and V + . Forth, as a consequence, the pairs (V + 0 , M 2 ) also exhibit strong spatial correlations, contrary to the case of BA ∞ . In these two latter permutations, it must be reminded that not all intensity maxima M may transform into vortices, because of the difference in spatial densities (Table 2) of these two populations of critical points [START_REF] Gateau | Complementary speckle patterns: Deterministic interchange of intrinsic vortices and maxima through scattering media[END_REF].

Next, we constrained the period to be equal to 4 by using the BA 4 mask. In this case, the WMND(Y 0 , Y ±4 ) are close to zero (transformation rate close to 1). As expected, this periodicity enhances the spatial correlation for the quadruplet (V - m-1 , M m , V + m+1 , S p m+2 ) and extends it to their 4 th spiral harmonics. Furthermore, in Fig. 3, strong correlations are observed for the pairs (V - 0 , M 1 ), (M 0 , V + +1 ) and (S p 0 , V - +1 ).

However, V + 0 is still observed to bifurcate between V - +1 and S p +1 with the same likelihood, similarly to the BA ∞ case. Therefore the cyclic permutation of four populations of critical points S p , V -, M and V + is not clearly established for BA 4 (conversely to the permutation obtained for BA 3 ).

When considering the addition of ±2-charged spiral masks for BA 4 , a transposition (2-cycle permuation) between V + and V -is clearly obtained (with a transposition rate below 1). A strong spatial correlation, reinforced as compared to BA ∞ , is also clearly observed between M and S P by addition of SP ±2k .

In summary, the results displayed in Fig. 3 reveals the fondamental topological transformations of critical points in a speckle with the addition of SP mask for a single-spiral mode aperture (BA ∞ ), and demonstrate the possible modification of topological transformation by the addition of BA masks with dihedral symmetry [See Supplement 1, Section 2.B for numerical confirmation of the experimental results and Section 2.C for details on transformations induced by BA masks with dihedral symmetries of orders higher than 4]. For the sake of simplicity, we chose here star-like amplitude masks with a dihedral symmetry which comprise spiral harmonics with equal amplitudes (at small enough spiral mode number) [See Supplement 1, Section 2.D for details on the influence of the width of the star branches].

Wavefield control in the vicinity of critial points

On a local scale, the transformations of critical points shown in Fig. 3 arise from the convolution in the imaging plane of the scattered field with the point spread function (PSF) associated with the combined amplitude and spiral phase masks.

Controlling the transformation of a critical point Y 0 of the complex wavefield A 0 , by adding a SP n mask (n ∈ Z * ) in a Fourier plane, requires that the PSF associated with the combination of the SP n and BA masks has a significant amplitude in the coherence area surrounding the critical point, or in other words: the area where the randomness of the speckle pattern has a limited influence as compared to the control by the incident wavefield. As a definition for the coherence area, we use the one proposed by Freund [START_REF] Freund | Optical vortices in gaussian random wave fields: statistical probability densities[END_REF] :

C area = (ρ(V + ) + ρ(V -)) -1 /2.
This definition avoids issues related to the shape of the aperture [START_REF] Freund | Optical vortices in gaussian random wave fields: statistical probability densities[END_REF] encountered when considering the area of the intensity autocorrelation peak [START_REF] Ochoa | Statistical properties of ray directions in a monochromatic speckle pattern[END_REF].

In our case, for all three BA masks, the coherence length was measured to be: C length = √ C area λ/(2.NA). Experimentally, the PSFs can be obtained (Fig. 4a) by computing the intensity cross-correlations of the measured speckle pattern I 0 and the measured speckle patterns I n associated with SP n masks (n ∈ 0; 6 ) and for the three BA masks. The mean values of I n were subtracted before computing the crosscorrelations. The intensity cross-correlations xcorr(I 0 ,I n ) are identical to the PSF of the combined BA and SP n masks. The centered spot of the autocorrelation xcorr(I 0 ,I 0 ) illustrates the spatial extent of the coherence area, and has the same dimension for all three BA masks since having the same radial aperture.

For BA ∞ , we observe that xcorr(I 0 ,I n ) has a circular symmetry with the highest values distributed on a ring whose radius (marked with a green line) increases with n (Fig. 4b), as observed for simple Laguerre-Gaussian beams [START_REF] Curtis | Structure of optical vortices[END_REF]. Interestingly, for n > 3, not only we observe that the ring radius is larger than twice C length but also that its amplitude is decreased to below 1/10 of the auto-correlation peak value (Fig. 4c). As a consequence, the transformation of critical points by applying SP n masks is inefficient (or "unlikely") and dominated by the surrounding random field. For this reason, no spatial correlation between pairs of critical points (Y 0 , X n ) could be found for n > 3. For BA 3 and BA 4 , the cross-correlation patterns xcorr(I 0 ,I n ) have dihedral symmetries D 3 and D 4 and a periodicity of N = 3 and 4, respectively. In both cases, the radial distance of the strongest peak remains below 1.4 * C length , and its amplitude always remains above 1/3 of the auto-correlation maximum value (Fig. 4c). The addition of the SP n mask thus allows controlling the field inside the coherence area surrounding critical points Y 0 , even for n > 3.

Conclusion

The critical points that naturally appear in a random wavefield can be transformed by the addition of a spiral phase mask in a Fourier plane. Here, we studied these transformations experimentally by imprinting spiral phase masks with a charge n ∈ -6; 6 to a laser beam impinging on a randomly scattering surface. In addition, these phase masks were combined with star-like amplitude masks with dihedral symmetries D 3 and D 4 in order to better control critical point transformations.

For a simple disk-shaped aperture carrying a single spiral mode n = 0, we experimentally demonstrated the topological correlation existing between the critical points of the initial wavefield A 0 , and the corresponding spiral transformed field A n . A partial transformation of vortices V - 0 into maxima M +1 was observed as well as a transformation of maxima M 0 into vortices V + +1 . Vortices V + 0 were observed to either correlate with phase saddle points S p +1 or with vortices of opposite sign V - +1 . For this statistical bifurcation, two transformation interpretations were suggested, calling for further future analytical studies. No simple topological correlation was found between the critical points of the wavefields A 0 and A n for |n| > 3. This result could be explained by the weak influence of spiral phase masks with a charge higher than 2 in the coherence area surrounding the critical points.

Furthermore, adding centered binary amplitude masks with dihedral symmetry D 3 or D 4 and Dirac-comb-like spiral spectra (of period 3 and 4), we demonstrated that it is possible to deeply modify the topological correlation between critical points. The observed changes arise from the introduction of a periodicity in the transformation between the critical points. We could thereby extend the correlation to spiral phase masks xcorr(I 0 ,I 6 )

BA 3 BA 4 BA ∞ c) BA ∞ BA 3 BA 4
Fig. 4: Cross-correlation (a) of the speckle patterns I 0 and I n with n ∈ 0; 6 . The cross correlations illustrate the spatial distribution of the intensity in the point spread function associated with the combined amplitude and spiral phase masks. The colorscale set for the autocorrelations xcorr(I 0 ,I 0 ) and was kept for the intercorrelations xcorr(I 0 ,I n ). A Nperiodicity is observed for the aperture BA N with N = 3 or 4.

Green circles mark the radial distance to the strongest peak. Radial distance (b) and normalized amplitude of the strongest peak (c) as a function of the charge of the spiral phase mask. The radial distance keeps increasing for the circular aperture BA ∞ while it remains close or below λ/(2.NA) (the coherence length of the speckle) even for SP n with n 3 for the aperture BA N and N = 3 or 4. The amplitude decreases for BA ∞ while it remains above 0.3 for BA 3 and BA 4 . with charges higher than 3, and reinforce some spatial correlation intrinsically present with a circular-aperture symmetry. For the amplitude mask with a D 3 -symmetry, a cyclic permutation between negatively charged vortices V -, maxima M , and positively charged vortices V + , is observed. For the amplitude mask with a D 4 -symmetry, phase saddle points participate as complementary points to complete the 4-periodic cycle. Considering the addition of 2-charged spiral masks, transpositions between V -and V + , and between M and S p were also revealed for D 4 -symmetry. The enhancement of the spatial correlation between the critical points of the wavefields A 0 and A n (compared to BA ∞ ) could be explained by the strong influence of the spiral phase mask in the coherence area surrounding each critical point, when the binary amplitude masks are added.

Here, cyclic permutations were controled using binary amplitude masks enforcing periodicity. Interestingly, our study may extend to other amplitude masks with a dihedral symmetry, such as polygonal [START_REF] Gao | Experimental demonstration of free-space optical vortex transmutation with polygonal lenses[END_REF] and triangular apertures [START_REF] Hickmann | Unveiling a truncated optical lattice associated with a triangular aperture using light's orbital angular momentum[END_REF], whose interactions with vortex beams were studied in free space. For N-gons, the N th spiral harmonics have a much lower amplitude than the fundamental spiral mode (n = 0). As a result, the spatial correlations between critical points for |n| > 3 is weaker than for BA N , and vanishes with the increasing charge of the SP mask [See Supplement 1, Section 2.E].

In a nutshell, we showed here that it is possible to manipulate the topological correlation between critical points and to control the transformation of critical points in random wavefields by combining amplitude masks and spiral phase transforms. Topological manipulation of critical points in random wavefields is of high importance to understand and control light propagation through scattering and complex media. The statistical study of correlations between permuted critical points provides a new tool to analyse seemingly informationless and random intensity patterns and thus to transmit information through complex media.

Supplementary information

The first section provides expended descriptions of the methods used to acquire and process the experimental data: details are given on the phase mask displayed on the spatial light modulator, the numerical treatment of the recorded speckle parterns is presented, and the algorithms for the detection of critical points are given. The second section contains numerical confirmations of the experimental results, as well as an extention (in numerical simulation) of the range of tested binary amplitude apertures: the Weighted Median Normalized Distances are shown for star-like amplitude masks with dihedral symmetries of orders higher than 4, for star-like amplitude masks with larger angular slits, and for polygonal apertures.

A Experimental data acquisition and processing

A.1 Phase mask on the SLM

The laser source was a laser diode emitting at λ = 635 mm. The laser beam was spatially filtered with a pinhole, expanded and collimated to cover the surface of the Spatial Light Modulator (SLM). The phase at the SLM was computed according to the scheme shown in Fig. 5. The scattering phase pattern at the SLM (or diffuser in Fig. 5) was engineered in order to both generate a fully developed speckle pattern at the camera plane and to minimize energy losses by controling the scattering angle. To achieve so, a complex matrix of dimension 396 × 396 pixels (half of the long-axis of the SLM) was generated numerically, with uniform amplitude and random phases evenly distributed between 0 and 2π. This random matrix was then multiplied by a disk-shaped amplitude mask of radius 60 SLM pixels (SLM pixel size : px 2 SLM = 20 × 20 µm 2 ) and Fourier transformed numerically to generate the random phase mask at the SLM plane. Provided the optical design of the setup and the focal distance of the lens f = 750 mm, this mask corresponds to a disk of diameter 60 400 × 2.λ.f px SLM 8.8 mm, slightly larger than the camera diagonal (7.62 mm), so ensuring uniform illumination of the camera chip (on average) and minimization of energy losses.

The BA amplitude mask were applied to the phase mask by setting the SLM phase to zero where the amplitude of the BA mask equals zero, and by adding a blazed grating where the amplitude of the BA mask equals one. The radius of the BA masks was chosen to have a speckle grain size of 15 camera pixels (camera pixel size: 4.65 × 4.65 µm 2 ). The blazed grating had a period of 4.34 SLM pixels (SLM pixel size : 20 × 20 µm 2 ), which yields a deflection angle of 7.31 mrad. With the focal distance of the lens f = 750 mm, this deflection angle provides a 5 mm shift of the zero order as compared to the speckle pattern. Because of the camera field of view (4.76 mm×5.95 mm), the undiffracted zero-order was centered at a distance 2.5 mm away from the camera sensor and could then efficiently be blocked with no spurious light observed on the camera.

For the reference beam, the same blazed grating was applied to a disk-shaped aperture separated from the scattering aperture. The radius of the disk was chosen so that the signal beam (scattered wavefield) and the reference beam have similar amplitudes at the camera plane. A defocus (Fresnel lens) was applied to the reference beam to cover the camera surface.

A.2 Numerical treatment of experimental data

The complex wavefield A n e iΦn at the camera plane was measured with the following procedure. The phase Φ n was obtained thanks to the intensity modulation induced by the phase-stepping interferometric measurements [START_REF] Robinson | Digital phase stepping speckle interferometry[END_REF][START_REF] Creath | Interferometry -phasemeasurement interferometry[END_REF]. The phase-stepping was performed by acquiring eight successive images with relative phase-shifts between the signal and reference waves: k × 2π/N , with N = 8 and k ∈ 0; 7 . The phase of the speckle pattern was obtained from phase-stepped intensity measurements:

I n,k = E R e 2ikπ N + A n e iΦn 2 (2)
= |E R | 2 + |A n | 2 + E R A n e -iΦn e 2ikπ N + E * R A n e iΦn e -2ikπ N
The phase can then be trivially retrieved by computing the argument of the following sum:

N -1 k=0 I n,k × e 2ikπ N = E * R A n e iΦn (3)
In our case, the reference beam did not exhibit a flat phase but a combination of a parabolic curvature as well as a relative phase tilt with respect to the speckled beam. Both of these profiles were numerically removed from the computed phase Φ n . The magnitude A n was not computed using this result to avoid amplitude uncertainties about the reference beam. Instead, A n was obtained from the signal speckle intensity I n without any reference beam: A n = √ I n . For all the apertures and all the speckle patterns, the Full Width Half Maximum (FWHM) of the point spread function (i.e. the speckle grain size) corresponded to 15 camera pixels, leading a fine spatial sampling of the speckle field. We set the length unit to 15 pixels and the corresponding spatial sampling frequency to 15 pixel -1 . For noise removal, the complex amplitude field was filtered in the Fourier domain by zeroing the values corresponding to spatial frequencies outside of the disk f 2

x + f 2 y 0.5 pixel -1 . Because of the spatial filtering by optical system, the frequencies outside of this disk contain only optical and electronic noise. We verified that the strongest spectral values of the complex field were within the disk. After this filtering process, experimental noise was not found to interfere with the detection of the critical points. A is applied to both beams in order to shift the speckle pattern of interest away from the unmodulated light. A defocus is added to the reference beam to cover the caerma surface, while the signal beams receives a random phase mask (or diffuser). Finally, the phase step and spiral wavefronts are added.

A.3 Detection of the critical points

The detection of critical points in the experimental speckle pattern is illustrated in Fig. 6. The pixel-precise locations of the critical points were determined in the spatial domain using the topology of each pixel neighborhood. For the intensity maxima M , the intensity saddles S I and the phase saddles S p , the eight neighbors of each pixel were used. We used the Matlab function developed by Tristan Ursell and made available in May 2013 (image extrema finder, https://fr.mathworks.com/matlabcentral/fileexchange/41955find-image-extrema). For the vortices V + and V -, we used the fact that the summed phase shifts on a closed loop around a vortex is greater than 2π. As a dicretised closed loop, we computed this sum over four neighbouring pixels, the phase differences between the adjacent corners of the square being computed and wrapped in the range -π, π . The summed phase shift was computed by adding the four phase differences, leading to a summed phase shift of ±2π around a vortex of charge ±1 and a summed phase shift close to zero elsewhere [START_REF] Brauer | Diffusers in digital holography[END_REF].

B Numerical simulations of scalar random wavefields

Numerical simulations of scalar random wavefields were performed to confirm our experimental results with larger statistics on critical points, independent random phase media and an independent method (section BB.2). Indeed, larger speckle patterns can be generated numerically directly from randomly distributed phase. Additionally, we extended our investigations to apertures with point group symmetries of orders higher than 4 (section BB.3), as well as to apertures with larger angular slits (section BB.4) and to polygon apertures (section BB.5).

B.1 Methods for the numerical simulation

The far-field of uniformly-illuminated apertures comprised of random phases (standard uniform distribution between 0 and 2π,) with spiral phase (SP) masks of order n ∈ -6; 6 , was computed to simulate the complex wavefield A n e iΦn . The speckle patterns were obtained by computing the twodimensional Fourier transform of the apertures addressed with random phases. As in the experimental procedure, an uncorrelated speckle pattern was also computed for each aperture. For all the apertures, the speckle grain size (FWHM) was set to λ/(2.N A) = 19 pixels where λ is the optical wavelength and NA is the numerical aperture, and a square grid of 17 mega pixels was computed. The pixel-precise location of the critical points were determined in the spatial domain using the topology of each pixel neighborhood, as for the experimental speckle patterns. These parameters lead to a count of critical points for the aperture BA ∞ of ∼ 10 4 vortices of each sign. The average number-densities of the critical points for the apertures BA ∞ , BA 3 and BA 4 are presented in Table 4. The slight differences in number-densities between experimental and simulation measurements can be attributed to the higher number of critical points in the simulation, which leads to better precision in the estimation of the critical point densities.

B.2 Validation of the experimental results

To support our experimental measurements, numerical simulations were performed for the apertures BA ∞ , BA 3 and BA 4 . The resulting Weighted Median Normalized Distances (WMND) are presented in Fig. 7. Despite the difference in the method to generate the speckle field, the size of the speckle grains, and the number of critical points involved in the statistical estimation of the WMNDs, our experimental results are remarkably consistent with numerical ones. Our experimental results are then validated in two regards. First, the results are robust to modifications in (i) the method to generate the speckle patterns, (ii) the scattering medium and (iii) the spatial sampling. Second, our experimental estimation of the WMND are statistically relevant since a higher average number of critical points did not lead to significantly different WMND.

B.3 Point group symetries of higher order

We investigated in numerical simulation the effect of an increase in the order of the dihedral point group symmetry. For this purpose, we used the apertures BA N with N ∈ {5, 6, 7, 8} ( point group symmetry D N ), as defined in the main article. Because of the strong anisotropy of the central peak of its point spread function, the aperture BA 2 could not be compared to the other apertures, and therefore was not computed.

The WMND are presented in Fig. 8. As expected for N= 5 or 6, the WMND(Y 0 , Y ±N ) are close to zero (transformation rate close to 1). We can also notice that, for BA 7 and n ∈ -3; 3 , as well as for BA 8 and n ∈ -4; 4 , the same WMND V + 0 , X ±n as for BA ∞ are found. The topological correlations for BA ∞ correspond to a single spiral mode, while the spiral spectrum of BA 8 contains harmonics. The first harmonics are the spiral modes m = ±8 and are for spiral modes high enough so that, for the pairs V + 0 , X ±n , the topological correlations of the fundamental mode (m = 0 as in BA ∞ ) dominates for n ∈ -4; 4 . For smaller N, the topological correlations between critical points are dictated by both the fundamental mode and the harmonics which leads to more complex WMND diagrams. This is this true in particular for BA 3 and BA 4 .

B.4 Influence of the width of the angular slits

The influence of the width of the angular slits was investigated in numerical simulation. Simulation were performed for periodic angular slits with a point group symmetry D 3 and three different widths of angular slits were tested : π 16 (i.e. BA 3 ), π 8 , and π 4 . The WMND are presented in Fig. 9. We can notice that the WMND(Y, Y ±6 ) increases with the increasing width of the slits, and even becomes equal to 0.5 for the largest slits. This results was expected considering the spiral spectrum of the apertures. All the considered apertures have harmonics at the spiral modes m = k.N with k ∈ Z. However, the amplitude of the harmonics is given by a sinc apodisation function which FWHM depends on the width of the angular slits. For the angular slit π 16 , the FWHM of the spectral apodization function corresponds to n= ±19. As a consequence, the aperture could be considered as spectrally invariant by the addition of SP n when n = ±3 or ±6. For angular slits of width : π 8 and π 4 , the FWHM of the spectral apodization function corresponds to n= ±10 and n= ±5, respectively. Consequently, the apertures are spectrally modified by the addition of SP ±6 .

In the experimental procedure, the width of the angular slits was chosen large enough so that the dark center induced by the spatial sampling of the SLM pixels has a radius below 2 SLM pixels. However, the width was chosen small enough to keep a spectral invariance in the range n ∈ -6; 6 . Both the dihedral symmetry and the width of the angular slit are important parmeters to observe the periodicity of the critical point transformation.

B.5 Polygonal apertures

Polygonal apertures with dihedral symmetries D 3 (equilateral triangle) and D 4 (square) were used to extend our study to other aperture shapes. The WMND computed from numerical simulations are shown in Fig. 10.

For the triangular aperture, topological correlations are observed for the pairs (M 0 , M ±3 ), but no cyclic permutation between V -, M and V + can be noticed. For the square, the WMND is very similar to the one of BA ∞ (disk). These results can be explained by the fact that, for N-gons, the N th spiral harmonics have a much lower amplitude than the fundamental spiral mode (m = 0). The spiral spectrum is therefore not invariant by the addition of a SP N mask. Moreover, for Fig. 8: WMND for all possible pairs of critical points screened here, for the addition of spiral phase mask with charges up to n =±6 and for the apertures BA N with N ∈ {5, 6, 7, 8}. The WMND were computed from numerical simulations.

Fig. 1 :

 1 Fig.1: Experimental setup (a) used to measure the intensity and the phase of speckle patterns corresponding to the different binary amplitude masks. A spatial light modulator (SLM) is illuminated with a collimated laser beam at 635nm. The phase is measured by phase stepping interferometry. Both the reference and the signal wavefronts are imprinted on the SLM. in addition to a blazed grating which allows sending undiffracted light to a beam block (BB) and the first order diffracted beam to a camera (Cam.). A stack of eight images was then sequentially recorded while phase shifting the reference beam. The measured intensity I 0 maps (top, green colorscale) and phase Φ 0 maps (bottom, gray colorscale) are presented for a circular aperture (BA ∞ ) (b), periodic angularslits with a point group symmetry D 3 (BA 3 ) (c), and periodic angular-slits with a point group symmetry D 4 (BA 4 ) (d). Miniatures of the BA masks are displayed for illustration.

Fig. 2 :

 2 Fig. 2: Statistical analysis of the separation distances between one set of critical points (here V - 0 ) and the closest point of another set (notated X). The distance between V - 0 and the closest X is notated d V - 0 , X . The Radial Probability Density Functions (RPDF) of the nearest neighbor (a) and the corresponding Weighted Median Distance (WMD) (b) are shown. Radial Distribution Functions (RDF) of the nearest neighbor (c) and the corresponding Weighted Median Normalized Distance (WMND) (d) provide a statistical toolbox to study the spatial correlation between pairs of critical points. The results were derived from experimental measurements of I n and Φ n obtained for the amplitude mask BA ∞ .

Fig. 3 :

 3 Fig. 3: Weighted Median Normalized Distance (WMND) for all possible pairs of critical points screened here and for the addition of spiral phase mask with charges up to n =±6. The WMND were computed from experimental measurements of I n and Φ n .

  0 ,I 0 ) xcorr(I 0 ,I 1 ) xcorr(I 0 ,I 2 ) xcorr(I 0 ,I 3 ) xcorr(I 0 ,I 4 ) xcorr(I 0 ,I 5 )

Fig. 5 :

 5 Fig.5: Computation of the phase displayed at the SLM. The SLM is virtually split into two parts of surfaces 396 × 396 pixels, each bearing pupils of different shapes, one for the reference beam (green), the other for the signal beam (blue). A is applied to both beams in order to shift the speckle pattern of interest away from the unmodulated light. A defocus is added to the reference beam to cover the caerma surface, while the signal beams receives a random phase mask (or diffuser). Finally, the phase step and spiral wavefronts are added.

Fig. 6 :

 6 Fig. 6: Detected critical points superimposed to the filtered experimental intensity (left) and phase (right) patterns obtained with the aperture BA ∞ . The displayed pattern correspond to an area of 150 × 150 pixels on the camera. Phase vortices (V + and V -) are also zeros of the intensity pattern. Notations for the critical points are identical to those introduced in the core of the main article.

Fig. 7 :

 7 Fig. 7: Weighted Median Normalized Distance (WMND) for all possible pairs of critical points screened here, for the addition of spiral phase mask with charges up to n =±6 and for the apertures BA N (N ∈ {3, 4, ∞}. The WMND were computed from numerical simulations.

Table 1 :

 1 Notations for the main critical points

	Phase	Maxima Saddle Vortices (charge ±1)
		-	S p	V -and V +
	Intensity Maxima Saddle	Zeros
		M	S I	V -and V +

Table 2 :

 2 Measured

		average number density of critical
	points (length unit: λ/(2.NA)). The average number
	of V -is 660.85 for the circular aperture (BA ∞ ).
	BA mask V -(or V + )	M	S p	S I
	BA ∞	0.19	0.32 0.36 0.65
	BA 3	0.20	0.39 0.36 0.79
	BA 4	0.19	0.33 0.40 0.70

Table 4 :

 4 Measured average number density of critical points from numerical simulations (length unit: λ/(2.NA)). The average number of V -is 9551.43 for the circular aperture (BA ∞ ).

	BA mask V -(or V + )	M	S p	S I
	BA ∞	0.20	0.32 0.39 0.67
	BA 3	0.21	0.38 0.39 0.75
	BA 4	0.20	0.34 0.43 0.70

Funding Information

This work was partially funded by the french Agence Nationale pour la Recherche (NEOCASTIP ANR-CE09-0015-01, SpeckleSTED ANR-18-CE42-0008-01).

WMND from numerical simulations

Fig. 9: Weighted Median Normalized Distance (WMND) for all possible pairs of critical points screened here, for the addition of spiral phase mask with charges up to n =±6 and for apertures with dihedral symmetry D 3 and angular slits of width : π 16 , π 8 , and π 4 , respectively (from top to bottom). The WMND were computed from numerical simulations. the triangular aperture, the point spread function (PSF) is an optical lattice which spatial extend expends as the charge of the SP mask increases [START_REF] Hickmann | Unveiling a truncated optical lattice associated with a triangular aperture using light's orbital angular momentum[END_REF]. Thereby, the maximum values of the PSF decreases, which strongly limits any possible control of the transformation of the critical points.

WMND from numerical simulations

Fig. 10: WMND for all possible pairs of critical points screened here, for the addition of spiral phase mask with charges up to n =±6 and for a disk aperture (top) and polygonal apertures : equilateral triangle (dihedral symmetry D 3 ) and square (dihedral symmetry D 4 ). The WMND were computed from numerical simulations.