Characterizing quantum gases in time-controlled disorder realizations using cross-correlations of density distributions

The role of disorder on physical systems has been widely studied in the macroscopic and microscopic world. While static disorder is well understood in many cases, the impact of time-dependent disorder on quantum gases is still poorly investigated. In our experimental setup, we introduce and characterize a method capable of producing time-controlled optical-speckle disorder. Experimentally, coherent light illuminates a combination of a static and a rotating diffuser, thereby collecting a spatially varying phase due to the diffusers’ structure and a temporally variable phase due to the relative rotation. Controlling the rotation of the diffuser allows changing the speckle realization or, for future work, the characteristic time scale of the change of the speckle pattern, i.e. the correlation time, matching typical time scales of the quantum gases investigated. We characterize the speckle pattern ex-situ by measuring its intensity distribution cross-correlating different intensity patterns. In-situ, we observe its impact on a molecular Bose–Einstein condensate (BEC) and cross-correlate the density distributions of BECs probed in different speckle realizations. As one diffuser rotates relative to the other around the common optical axis, we trace the optical speckle’s intensity cross-correlations and the quantum gas’ density cross-correlations. Our results show comparable outcomes for both measurement methods. The setup allows us to tune the disorder potential adapted to the characteristics of the quantum gas. These studies pave the way for investigating nonequilibrium physics in interacting quantum gases using controlled dynamical-disorder potentials.


Introduction
Disorder is ubiquitous and can be found in a large variety of systems ranging from, e.g., fluid dynamics [1], diffusion processes [2] to solids [3,4], biology [5][6][7], spectroscopy [8,9], or towards simulations for quantum computing [10][11][12].The field of ultracold gases provides a versatile tool for combining quantum gases with disordered systems [13], offering a highly controllable model system for a broad range of applications.Soon after the first experimental realizations of disorder potentials for ultracold gases [14], observation of Anderson localization (AL) for quantum gases was also achieved experimentally [15].Initially predicted in the context of electron transport in crystals [16,17], AL describes the absence of wave propagation in disordered media due to destructive interference of possible paths originating from multiple scattering.While the first observations of AL with quantum gases were made in one dimension, localization can also occur in higher dimensions [18][19][20][21][22][23][24].Also, related disorder-induced observations, such as many-body localization in interacting, disordered lattice systems, have been reported [25][26][27][28][29].Moreover, disorder has proven valuable combined with topological transport studies, such as Thouless pumping [30,31].
An important recent research direction of ultracold quantum gases is the study of nonequilibrium dynamics in driven quantum systems.A key asset to studying such systems is tight control over all relevant experimental parameters.For dynamicallydisordered systems, this implies controlling the correlation time of the disordered potential in addition to engineering the spatial correlation length.For disordered systems, however, most experimental realizations considered a purely static disorder, and examples of dynamical disorder studies are scarce.For example, investigations of dynamical disorder have been reported using microwaves [32], where the time scales of the disorder dynamics are in the millisecond range.These time scales, however, are not suited for quantum gases, where relevant time scales are of the order of microseconds.The impact of time-dependent disorder on quantum systems is yet largely unexplored experimentally.Theoretical works predict many-body localization to resist turbulations created by time-dependent potentials [33][34][35] given that the driving frequency is sufficiently large.Theory also predicts universal transport in dynamic disorder in the absence of localization [36].At the same time, as shown for an optical setup in [37], AL seems to break down, and diffusive transport or hyper-broadening is expected to occur in time-dependent disorder.So far, it is unclear under which conditions localization could persist for slower dynamics.Furthermore, there are still open questions, e.g., on which timescales these processes take place and how the driving parameters influence the dynamics.Here, we show a simple scheme to create time-dependent optical disorder potentials where the spatial statistical properties of the potential are time-independent.To analyze the disorder properties of this scheme, we here experimentally measure the optical intensity cross-correlations, relating intensity distributions of two different disorder realizations, which decay as the disorder realization is experimentally changed in discrete steps.We address the extent to which these optical cross-correlations are reflected by the cross-correlations of quantum gas densities measured in different static disorder realizations.
While the former is a direct measurement of the disordered optical potential's cross-correlation, the latter is influenced by several effects.Theoretically, the density distribution of a BEC in disorder due to interaction effects has been studied; see Ref. [38].Experimentally, we address the influence of column integration, imaging noise, inhomogeneities, and finite temperature in images of quantum gases and whether the cross-correlations between density distributions taken in different realizations of the underlying disorder potential can be deduced.We find that, indeed, the atomic density cross-correlations between two absorption images of quantum gases trapped in different optical disorder realizations resemble the direct optical correlations.This result will allow the study of quantum gases in dynamical disorder when the optical speckle correlation is continuously reduced during the trapping of a quantum gas in future work.

Static and dynamic optical speckle
The prime example of random disorder phenomena is optical speckles.Speckles occur as interference fringes of random scattering at rough surfaces or within diffusive media [39].Applications include medical imaging [40], optical coherence tomography [41] or astronomy [42], but also forming disordered potentials for ultracold gases [43][44][45].For quantum gas applications, optical speckle patterns can conveniently be realized by transmitting a laser beam through a diffuser, e.g., a ground glass plate.This simple approach has evolved to an established method for optical disorder production [39,46,47].Microscopically, the surface structure acts as a random collection of scattering centers imprinting a random amplitude and phase distribution onto the light field as shown in Fig. 1.
When a collimated Gaussian-shaped laser beam with waist w d is transmitted through a diffusive plate, the light gets diffracted at every point of the intensity profile, and an angular distribution is formed with cone angle Θ d .An imaging lens with focal length f creates a disordered intensity distribution in the focal plane.This distribution in the imaging plane consists of individual speckle grains and has a waist of For the application as an optical potential, especially the local spatial and temporal properties of the speckle are important [44].The transversal correlation length σ s of the speckle measures the average size of the speckle grains in the image plane or, equivalently, the length scale on which the speckle preserves its coherence.In the focal plane of a lens, the grain size is minimized in the transversal direction.Here, the imaging resolution limits the achievable correlation length, which can be a few hundred nanometers.The speckle-grain extension in the longitudinal direction, i.e., along the propagation axis, is much larger and determined by the Rayleigh length of the focusing system.Experimentally, the speckle's correlation length can be determined from arbitrary light distributions by calculating each dimension's angle autocorrelation function of the speckle pattern.Knowing the magnitude of the correlation length is essential for using the speckle as a dipole potential for quantum gases.For a Bose-Einstein condensate, a relevant length scale is the healing length [48] ξ = 1/ √ 8πn 0 a of the gas in the trap center with peak density n 0 of the gas and s-wave scattering length a.The healing length represents the smallest length scale on which the condensate's wave function can react to perturbations.In our case, the correlation length σ s is larger than the healing length ξ by a factor of two, allowing the condensate to resolve all details of the speckle.Since the speckle imposes a potential onto the atomic cloud, it imprints a density modulation on the condensate.We can infer information about the speckle potential by recording the column densities of the condensate [38,49].Furthermore, for an ideal, fully-developed speckle, the distribution of intensities is the same for all realizations [39].It is given by the negative exponential intensity-probability-distribution [42] with intensity I = I(⃗ r) at the point ⃗ r.This leads to the fact that high-intensity speckle grains are scarce.The average intensity of a speckle realization is relatively low.A typical measured intensity-probability distribution and an intensity distribution for one column of a typical speckle image are shown in Fig. 2. The properties considered so far hold for speckle patterns in general.In our experiment, the speckle pattern can be applied as a static potential, or be changed in time.To this end, speckle patterns can, among others, be manipulated by modifications of the diffusive medium.In our case, the speckle is created with the help of an optical diffuser.The speckle pattern can change, e.g., by movements or rotation of the diffusive medium [50], which causes decorrelation effects.However, the entire speckle pattern is shifted or rotated for a simple translation or rotation.In our setup, we supplement the scheme shown in Fig. 1 by an identically-constructed second diffuser (Edmund Optics #47-991) in a distance of up to 10 cm to create a speckle pattern that can be changed via the relative angle between the diffusers (Fig. 3).Each diffuser plate imprints a random phase distribution onto the incoming laser beam.The first diffuser rotated around the principal optical axis, produces a rotating speckle pattern onto the second diffuser.This, in turn, imprints another random phase-and amplitude distribution onto the rotating distribution.This combination leads to a disordered interference pattern that can change depending on the relative rotation angle of the two diffusers around the principal optical axis.That changes the position of the speckle grains and their intensity.As a side effect, the second diffuser also widens the beam further, slightly reducing the transmitted light's intensity.The speckle pattern can change stepwise by rotating one of the diffusers for a specific angle or dynamically by rotating the diffuser with a certain angular velocity.In the present work, we use a stepwise change between different realizations of one experimental run for characterization purposes.The continuous, dynamic change enabled by this experimental scheme allows us to study quantum gases in dynamical disorder, which is beyond the scope of this work.
To characterize the speckle pattern produced independently from the quantum gas, a test setup was built that simulated the characteristics of our quantum gas setup.For details on the main setup, see Refs.[51,52].The test setup includes a laser beam of wavelength λ = 532 nm illuminating the abovementioned diffuser setup with waist w d , see Fig. 1.The speckle pattern results from imaging this randomized phase distribution using a focusing lens with a numerical aperture of 0.3.The image plane corresponds to the position of the molecules in the quantum gas setup.A microscope (Olympus RMS20X and Thorlabs TTL180A) with a magnification factor of 20 maps the pattern on a CCD camera (FLIR Blackfly BFS U3-63S4M-C), which is used to record and characterize the optical intensity distribution.A diffuser itself brings its own correlation length σ d , which is given by its surface structure.It can be calculated as [53] with wavelength λ of the incoming light.Including a second diffuser to the setup leads to a modified correlation length σ D of the whole setup compared to just using one diffuser in Eq. ( 3).The correlation length for configurations with two diffusers of different divergence angles can be assumed in a simple model as with the compound correlation length σ D , the single correlation lengths σ d,i corresponding to each plate, wavelength λ and the half divergence angles Θ d,i .

Ex-situ characterization of the correlation angle
As a first speckle characterization, we measure the speckle's transversal correlation lengths of the intensity distribution produced by the combined diffuser setup and characterize the grain size.We computationally analyze the autocorrelation of the speckle patterns, yielding a Gaussian decay of correlation in each direction.The transversal correlation length can be extracted by fitting a Gaussian function through the center of the two-dimensional angle autocorrelation function.The full width at half maximum of that function measures the speckle's grain size.These measured correlation lengths are not necessarily equal in the two directions of an image because  of the inhomogeneities of the camera sensor, the imaging system, and the speckle itself.
In our case, the correlation length was determined to be σ sx = 558.7 ± 3.6 nm in the x direction, and σ sy = 554.9± 3.6 nm in the y direction.We used the method described in [51] for the determination.
With our setup, we can, in particular, provide and control the decorrelation of optical speckles.Experimentally, this decorrelation is realized by the rotation of one of the diffusers by a given rotation angle, and it manifests itself in a change of the speckle pattern.However, specific details of the experimental setup will influence the particular decorrelation behavior.We determine different influencing factors and their impact on the speckle pattern by varying parameters in the experimental setup, including the beam waist illuminating the diffuser and the distance between the diffusers.For each of the following measurements in this section, the first diffuser is rotated in steps between 0.04 • and 0.10 • while the other is kept static.This method allows evaluating the most suitable parameters for later use in the quantum gas setup.The decorrelation can be quantified by numerically cross correlating pictures taken of the initial optical intensity pattern and after rotation of the diffuser by a given angle.This function is given by with rotation angle Θ, the speckle intensity distribution I Θ and position vector ⃗ r = (x, y).Eq. ( 5) therefore quantifies the change between two speckle patterns differing by the rotation angle Θ.In the image plane, the speckle pattern changes with rotation, and new peaks in the optical intensity appear while others move and eventually vanish.This is shown in Fig. 4 (b) for rotation angles of 0.00 • , 0.34 • , 0.60 • and 1.00 • .Plotting the maximum of the intensity cross-correlation versus the rotational angle results in a typical decorrelation curve [50] having the shape of a Gaussian function as exemplified in Fig. 4.This yields a characteristic angle at which the correlation has decayed to half its maximum (grey vertical line).That angle we define as the (rotational) correlation angle.Details of the illuminating beam and the diffusers determine this angle.To quantify this connection with a simple model, we use a geometric approach: In our setup, a rotation of the first diffuser by a small angle α causes a local displacement of If this shift is larger than the correlation length of the other diffuser, the amplitude and phase distributions and, therefore, the whole speckle pattern behind the two diffusers is significantly changed, hence decorrelated.For a phase distribution shift ∆s > σ d , this leads to a correlation angle In the remainder of this work, the correlation angle will be changed in discrete steps, and molecular BECs will experience a static disorder potential.The decorrelation can be precisely tuned by experimentally modifying the beam diameter of the laser beam passing through the diffusers as shown in Fig. 5.For example, the distance between the two diffusers and the diffusive angle will influence the correlation angle according to Eq. ( 1) and Eq.(3).Intuitively, the scattering of the first diffuser will lead to an additional divergence of the laser beam, modifying the beam diameter on the second diffuser depending on the distance.While this widening is hardly significant for small diffuser distances, it must be considered for large distances.The correlation angle as a function of beam diameter is shown in Fig. 5 (a).Additionally, the expectation from equation Eq. ( 7) is shown as solid orange line with π•σ d w d .The curves are in good agreement.However, the diffusers' divergence angles are another free parameter, see Eq. ( 4).The dependence of the correlation angle on σ d , i.e., the combined correlation length of the two diffusers, is presented in Fig. 5(b) for different combinations of diffusers with divergence angles of 0.5 • , 1.0 • , 5.0 • and 10.0 • (Edmund Optics #47-991).Here, the diffuser with the smaller divergence angle is placed first in the beam path for each combination.The resulting change in the correlation angle with combined correlation length is linear, indicated by the solid orange line in Fig. 5(b).Finally, the influence of the distance between the diffusers on the correlation angle is studied.As shown in Fig. 5 (c), the combination of different diffusers mainly influences the correlation length.The correlation angle stays the same within the errors for the distance in question.The combinations with smaller divergence angles, especially the combination of two 0.5 • divergence angles, lead to larger correlation angles.Moreover, the order of the two diffusers plays a role if diffusers with different divergence angles are used.If the plate with a smaller divergence angle is placed in front, this leads to a larger correlation angle.
In the following part for measurements with a molecular BEC, the distance between the two diffusers is approximately 5.5 cm.This is sufficient for optical access to the system.Two diffusive plates with divergence angles of 1.0 • (Edmund Optics #47-991) are chosen to compromise a large correlation angle and the desired beam width after the diffuser setup.With the beam waist given in the quantum gas setup approximately 4 mm, the test setup suggests correlation angles of roughly 0.5 • -0.6 • .

Preparation of a molecular Bose-Einstein condensate
We produce a molecular BEC comprising fermionic 6 Li atoms prepared in the two lowestlying hyperfine substates of the electronic ground-state 2 S 1/2 [51].The experimental sequence starts with laser cooling of atoms in a Zeeman slower and subsequently in a three-dimensional magneto-optical trap.After laser cooling, we evaporatively cool the sample at an external magnetic field of 763.6 G, close to a broad Feshbach-resonance centered at 832.2 G, which we use to tune the inter-particle scattering length.At the end of the evaporation, the final dipole trap power of 10 mW, corresponding to a trap depth of k B × 500 nK, is held constant for 100 ms, producing approximately N ≈ 4.3 × 10 5 bosonic molecules at temperatures of about 50 nK.Here, the s-wave scattering length lies at a = 2706 a 0 [54], resulting in a chemical potential of µ = 250 nK × k B and a healing length of ξ = 270 nm, which is smaller than the shortest correlation length of the speckle.Hence, the wave function of the molecular BEC resolves the microscopic details of the potential [38].The trapping frequencies of the combined optical and magnetic traps are (ω x , ω y , ω z ) = 2π × (195 Hz, 22.6 Hz, 129 Hz) with axial confinement in y direction by a magnetic saddle potential.The many-body relaxation time is h/µ ≈ 190 µs.
After preparing the molecular cloud, we apply the optical speckle pattern discussed above, serving as a random potential for the molecules.The speckle pattern here is created with a Gaussian laser beam at a wavelength of 532 nm as in the test setup, using a power of up to 10 W. We use the same ground-glass plates as tested before and insert one in a high-precision rotation stage (Owis DRTM65-D35-HiDS), providing rotation velocities up to 2100 • /s, as well as discrete steps with a minimal rotation angle of 0.036 • .The other diffusive plate is kept static.This combination provides a controllable dynamic disorder potential, which is focused on the position of the molecules.The static speckle potential is slowly introduced during 50 ms of linear ramping to minimize the creation of excitations in the system.Subsequently, it is held for 100 ms.The gas is trapped in the combined magnetic trap and speckle potential to allow the molecules to react to the disorder potential and equilibrate.Finally, we image the cloud in-situ to infer the column density distribution via resonant high-intensity absorption imaging.Before the next quantum gas is prepared, the diffuser is rotated by a given angle.Importantly, during every experimental realization, the molecular BEC is trapped in a static disorder realization, and a discrete angle changes the specific realization before the next quantum gas is prepared.Studying the evolution of quantum gases in continuously changing disorder, i.e., when the disorder realization and hence the disorder correlation changes while the molecular BEC is trapped, is beyond the scope of this work.

In-situ measurement of speckle waist
Inside the vacuum chamber, the magnetic field center determines the exact cloud's position.Ideally, this position should coincide with the center of the speckles' Gaussian envelope.To ensure both speckle and trap are adequately overlapped, we employ the loss of molecules upon quickly switching on the speckle, which depends on the local strength of the disorder potential.We characterize the loss of molecules depending on the speckle's exact positioning and the laser power creating the speckle.For the adjustment, a tiltable mirror between the diffusers and the objective is used to slightly alter the angle along both available directions at which the speckle beam enters the objective.With this procedure, we aim at the largest-possible molecule losses the speckle can induce when the speckle potential is suddenly switched on.Different from the usual measurement sequence, we switch the speckle rapidly on a time scale of less than 1 µs, and the freeevolution time in the speckle is 400 ms.For maximum losses, the disorder potential is centered on the molecular cloud.Once this optimal adjustment has been found, we map the fraction of molecules lost to the total power of the speckle pattern as shown in Fig. 6 (a).With this mapping, we can measure the waist of the disorder potential by shifting the relative position of speckle and cloud positions along two orthogonal directions via the objective; the position-dependent losses for constant speckle power as shown in Fig. 6  distribution.We extract the beam waist by fitting the data taken with a Gaussian curve.We find waists of 466 ± 25 µm in one direction and 414 ± 25 µm in the orthogonal direction.The speckle size is thus large enough for the molecular cloud with an extension of ∼ 200 µm on its long axis.

In-situ measurement of the correlation angle
In order to study the influence of the optical speckle potential on a quantum gas, we apply the speckle pattern and adjust it to the position of the molecular BEC as described above; for details, see [55].Focused onto the molecular cloud in the center of a vacuum system, the speckle pattern can not be characterized directly but only via its impact onto the ultracold cloud via the dipole force.However, detailed in-situ knowledge about the local speckle properties is crucial beyond the global properties, such as its waist, because small changes in the optical imaging system might lead to a significant change of, e.g., the correlation length.We explain the reasons for our assumption that optical and density cross-correlation and hence correlation angle are comparable in detail in the Appendix 7.
To compare the correlation properties inferred from molecular density distributions to the values obtained from optical measurements, we prepare molecular BECs described above and expose them to the disorder potential, ramped up linearly within 50 ms.After reaching the final disorder strength of 646.5 nK × k B , the cloud is held in the speckle for 100 ms, allowing the system to equilibrate.Finally, the sample's column density distribution n col (x, y) is measured using high-intensity-corrected absorption imaging [56].This procedure is repeated 380 times, where in-between repetitions, the angle φ of the second diffuser is increased in steps of 0.036 • up to a rotational angle of 13.640 • .For identical parameters, we repeat these measurements four times.That yields the column density distribution n col (φ, x, y) for each angle φ.One exemplary absorption image is shown in Fig. 7.We first fit the absorption images of the disturbed clouds with a two-dimensional Thomas-Fermi profile n fit (x, y).Fitting the corresponding molecular cloud ensures to avoid negative influences of shot-to-shot fluctuations compared to subtracting the density profile of a cloud without the speckle potential.The difference between this fit and the measured distribution yields the speckle-induced density fluctuations ∆n(x, y) = n col (x, y) − n fit (x, y) quantifying the speckle patterns' influence.The influence of the diffusers' relative angle is shown in Fig. 8 (a), showing the positionresolved fluctuations ∆n(x, y) for the three rotation angles of φ = 0.00 • , φ = 0.36 • , and φ = 0.72 • .Density maxima from the first picture move or vanish as the diffuser rotates to larger angles.In contrast, new maxima appear in other positions, as seen in the exsitu measurements.Although the basic structure of the original pattern can still be seen in the third image (φ = 0.72 • ), several differences have developed.With further rotation of the diffuser, these differences increase.Only random correlations to the initial pattern are found for rotation angles larger than 4.5 • .We quantify these differences between images for each rotation angle φ to the initial one via the autocorrelation within one image (for δφ = 0 • ) or the cross-correlation between two images (for δφ ̸ = 0 • ) in twodimensional space x and y and angle φ.Fig. 8 (b) shows the molecular cloud's positionresolved autocorrelation of the initial cloud (i.e., φ = 0 • ).The central (x = 0 = y) auto-correlation peak decreases height with increasing rotation angle.We normalize the peak height of the different images to this initial peak height of the autocorrelation function.The normalized autocorrelation function in the angle and space domain is given by AC ∆n (δφ, δx, δy) = (x i ,y j ) ∆n(φ, x i , y j )∆n(φ + δφ, x i + δx, y j + δy) where the summation runs over all pixel (x i , y j ) of the recorded image.We note that for δφ = 0, the autocorrelation of a single image is formed (Fig. 8 (b)).At the same time, for δφ ̸ = 0, the cross-correlation between two density distributions trapped in different speckle realizations is computed (Fig. 8 (c)).The power spectral density (PSD) in kspace can be calculated via the Fourier transformed density of the difference images with [57] PSD With this, we can calculate the autocorrelation function, which is proportional to the inverse Fourier transformation of the PSD AC ∆n (δφ, δx, δy) ∝ F −1 (PSD(k x , k y )) (10) From the decay of the maximum of each of these autocorrelation functions AC ∆n (δφ, 0, 0) as a function of the relative rotation angle δφ of the diffusers, the correlation angle can be determined, see Fig. 8 (c).
We again define the correlation angle as the angle where the cross-correlation has decayed to half the initial value at an angle of φ ≈ 1.14 • .This angle differs from the correlation angle determined ex-situ but is in the same order of magnitude.A combination of experimental limitations may explain the difference.The first constraint originates from limits in absorption imaging.We can only take two-dimensional images of the actual three-dimensional molecular cloud through absorption imaging.The third dimension of the cloud along the imaging direction is integrated over, which influences the density cross-correlations measured.Furthermore, the speckle grains are inhomogeneous.While they have their long axis along the direction of imaging, the extension of the molecular cloud is still larger than the length of the speckle grains as about two grains fit into the cloud along the imaging direction.Moreover, the positioning of the grains along the imaging direction is random.Thus, some lie only inside a small part of the cloud, while others may be entirely immersed in the cloud.Due to the two-dimensional image, it is unclear how many speckle grains are causing the measured density fluctuations.Some speckle grains are thus not perceived, systematically changing the correlation properties measured.In other cases, independent speckle grains may lie very close to each other and be perceived as a single grain in density-density cross-correlated images.This leads to a systematic error in which only a part of the speckle is perceived.Thus, the effects of the optical speckle on the molecular density distribution in the absorption image are reduced, causing an apparently larger decorrelation angle.

Conclusions
In conclusion, we have constructed a setup to generate time-dependent disorder with constant correlation length but time-controlled correlation angle.This method allows adjusting all relevant parameters of an optical disorder potential, including the potential depth and the time scales of the disorder dynamics, as required.We characterized a BEC in varying (but static) disorder realizations to extract the disorder correlation angle.Our results confirm that the density-density cross-correlations obtained from different quantum-gas images used as a probe for the disorder potential show the same behavior as the intensity distributions we recorded ex-situ.However, the absorptionimaging properties of the quantum gas limit the in-situ method, leading to apparently larger correlation angles.
An appealing extension of our scheme is a continuous evolution of the disorder, where the disorder realization changes while one and the same quantum gas is trapped.For a diffuser continuously rotating at angular frequency ω r , we can define the correlation time as the temporal counterpart of the correlation length via The correlation time determines when a strong decorrelation is achieved, and the disorder potential loses resemblance to the original intensity pattern.From Eq. ( 11), using typical experimental parameters, we can also quantify the correlation times accessible in our setup.For angular velocities of up to 2100 • s −1 , we find correlation times τ > 543 µs.This time scale is bigger than the corresponding minimum time scale of the quantum gas given by the inverse chemical potential h/µ = 190 µs.Hence, the gas can always adapt to the speckle potential's changes for this largest chemical potential.For a smaller chemical potential, however, we can realize an experimental situation where the many-body state of the gas cannot follow the fast changes of the speckle potential.
Because of the systematic errors in the measurement of correlation angles via density cross-correlations of the molecular absorption pictures, we assume the real correlation angle to be smaller than 1.14 • , which leads to shorter correlation times.However, we can not resolve this any finer with the current technical constraints.With the ex-situ determined correlation angle of 0.49 • , the minimal correlation time would be roughly 233 µs.This method of continuously transforming one disorder realization into another with an adjustable correlation time resembles naturally fluctuating disordered systems.It is, therefore, well suited for future simulation of quantum dynamics, including, for example, the motion of ions in a fluctuating polymer background [58], as found in Lithium-polymer batteries.
Complementary methods to model time-dependent disorder, e.g., ionic transport through biomembranes [59,60], use dynamic-disorder hopping models [61][62][63], where the disorder instantly changes after certain renewal times.While both models are expected to agree for short renewals and long observation times, this is not the case for the short-time behavior, i.e., when individual renewals and their time scales matter.These situations could be realized using digital mirror devices [27] or spatial light modulators [20], allowing for experiments with instantly changing potentials and different random patterns.
Our characterization measurements point towards realizing transport studies in time-dependent and controllable disorder in strongly interacting quantum gases, extending the initial work on superfluid transport in disorder [64], or in interacting thermal clouds.An interesting question concerns the role of localized, noninteracting quantum systems in disordered potentials [13,65,66] when the disorder is rendered time-dependent.While in the asymptotic limit, localization is expected to cease in time-dependent disorder [33][34][35], future experiments in our system might reveal how dynamics ensue in an originally localized quantum system, how the time scales and mobile fraction depend on the properties of the disorder potential, and if the microscopic time scales responsible for breaking the localization can be identified.where x 0 is the initial position and AC(•) denotes the auto-correlation function.The P SF is the same for all images because the basic optical configuration remains the same for all images and rotations.Therefore, its autocorrelation function also contributes a constant factor to the integral, and Section 7 leads to This means that the image correlations reflect the speckle pattern correlations except for a constant factor originating from the imaging system.With this, we can extract the correlation properties of the speckles from the in-situ density distributions of the molecular clouds.
The correlation angle is calculated as the angle where the intensity-intensity correlation function I Θ (x 0 ) has decayed to half the maximum value.One crucial point is that this only allows deducing and comparing the amount of correlations between two images.The opposite, i.e., inferring the complete information of the whole image from the correlation function, is not possible, especially since information on the molecular clouds can not be reconstructed.Reconstructions of the original images are possible but computationally intensive and usually require more detailed knowledge of the specific PSF [50,[72][73][74].To reconstruct the original speckle pattern (or density distribution) using this setup, one would need the exact point-spread function and not just a rough estimate.Nevertheless, our technique allows comparing the clouds for the molecular density.The inversion of the density distribution's minima to the maxima of the speckle does not change the differences in the correlation between two pictures evaluated the same way.

Figure 1 .
Figure 1.Schematic drawing of a simple diffuser setup consisting of a diffuser and a focusing lens.A collimated Gaussian laser beam with waist w d illuminates the diffuser, which imprints random amplitudes and phases onto the light.The beam is diffracted into a cone with a half-opening angle Θ d .Subsequently, an objective lens focuses on it with a focal length of f.The beam has a Gaussian-shaped intensity distribution envelope with the waist w s in the focal plane.

Figure 2 .
Figure 2. Defining speckle characteristics.(a) intensity-probability distribution for a recorded speckle image as an absolute count rate.The orange line depicts an exponential fit to the semi-logarithmic plotted measurement data.(b) One-dimensional intensity distribution for a speckle image with the mean value indicated as an orange dashed line.(c) 3d image of a typical speckle pattern of size (384 × 384) µm, where the speckle's correlation length is (575 ± 15) nm.

Figure 3 .
Figure 3. Sketch of the diffuser setup.Coherent light with a wavelength of λ = 532 nm is transmitted through a rotating diffuser, which widens the beam and imprints a random and rotating intensity distribution onto the beam.That illuminates a second diffuser, again imprinting a random phase and amplitude distribution.Combining two diffusers in a row leads to a dynamically changing disorder pattern.

Figure 4 .
Figure 4. Changes in the speckle patterns for diffuser rotation.(a) The maximum value of the cross correlation function of the speckle intensity distribution for different rotation angles with two diffusers of a 1 • divergence angle each.The blue points show the measurement data.Error bars mark the uncertainty of a fit to extract the maximum value from the ex-situ correlation.A solid orange line shows a Gaussian fit curve to the data points.The gray line marks the correlation angle at which the autocorrelation has dropped to the half maximum.bars denote the fitting error.(b) Contour plot of an extract of the resulting speckle pattern for the measurement points depicted in (a) with rotation angles of 0.00 • , 0.34 • , 0.60 • , and 1.00 • .

Figure 5 .
Figure 5. Influence of different setup changes on the correlation angle.Errors bars correspond to fitting errors.(a) Changes in the diameter of the incoming laser beam w d on the diffusers.The orange line depicts a theoretical course θD 2•w d , the blue dotted line a 1/x-fit.(b) Influence of the correlation length σ d .The orange line depicts a linear fit to the data points.(c) Influence of the distance between the two diffusers.Measurements for different combinations of diffusers are drawn.

Figure 6 .
Figure 6.Determination of speckle waist and alignment.(a) Losses in the molecular cloud versus speckle laser power.Error bars denote the standard deviation of five repetitions.The laser power measures the speckle intensity with the speckle potential centered at the molecular cloud.The solid line is a fifth-order polynomial function fitted to the data.(b) and (c) Disorder intensity envelope along two orthogonal directions in the focal plane of the objective.The speckle pattern is shifted with the same laser power, and the distribution is acquired via the mapping from (a).The error bars show the standard deviations for five repetitions.The resulting intensity distribution is fitted with a solid orange line as a Gaussian function, which gives waists of (b) 466 ± 25 µm and (c) 414 ± 25 µm.

Figure 7 .
Figure 7. Absorption image processing.Each recorded absorption image n col (x, y) (left) is fitted with a smooth, two-dimensional Thomas-Fermi profile n fit (x, y) (center), which drops to zero at the indicated Thomas-Fermi radii R x/y .Subtraction of the fitted profile from the measured column density distribution yields the speckle-induced density fluctuations ∆n(x, y) = n col (x, y) − n fit (x, y) (right), pixel values outside the Thomas-Fermi profile are set to zero in ∆n.

Figure 8 .
Figure 8. In-situ measurement of the correlation angle.(a) Speckle-induced density fluctuations for different rotational angles of the diffuser.Small rotations (≪ 1 • ) already lead to significant changes in the difference images.(b) Angle autocorrelation AC ∆n (δφ = 0, x, y) calculated from the column density distribution n col (φ, x, y) for the initial rotation angle with φ = 0 • .Without rotation, corresponding to the left image from (a), it shows one significant peak in the middle of the image.With increasing rotation angle, the angle autocorrelation function's peak height decreases.(c) Excerpt of angle autocorrelation function calculated from the first 192 column density distributions with φ = 0 • , . . ., 6.876 • .As the angle of rotation increases, the speckle pattern changes, and with it, the height of the autocorrelation maximum decreases in a Gaussian manner and approaches a finite value.This offset is due to random residual correlations between the images.The numbers indicate the corresponding values to the images from (a).The gray line marks the correlation angle of 1.14 • , at which the correlation function has dropped to half its initial value.