A stable phase-locking-free single beam optical lattice with multiple configurations

Yirong Wang    \authormark1,2 Xiaoyu Dai    \authormark1,2 Xue Zhao    \authormark1,2 Guangren Sun    \authormark1,2 Kuiyi Gao    \authormark1,2,* and Wei Zhang\authormark1,2,3,\dagger \authormark1School of Physics and Beijing Key Laboratory of Opto-electronic Functional Materials & Micro-Nano Devices, Renmin University of China, Beijing 100872, China
\authormark2Key Laboratory of Quantum State Construction and Manipulation (Ministry of Education), Renmin University of China, Beijing 100872, China
\authormark3Beijing Academy of Quantum Information Sciences, Beijing 100093, China
\authormark*[email protected], \authormark\dagger[email protected]
journal: opticajournalarticletype: Research Article
{abstract*}

Optical lattices formed by interfering laser beams are widely used to trap and manipulate atoms for quantum simulation, metrology, and computation. To stabilize optical lattices in experiments, it is usually challenging to implement delicate phase-locking systems with complicated optics and electronics to reduce the relative phase fluctuation of multiple laser beams. Here we report a phase-locking-free scheme to implement optical lattices by passing a single laser beam through a prism with n𝑛nitalic_n-fold symmetric facets and large apex angles. The scheme ensures a stable optical lattice since the interference occurs among different deflected parts of a single laser beam without any moving component. Various lattice configurations, including a triangular lattice and a quasi-crystal lattice with ten-fold symmetry are demonstrated. In both cases, stability measurements show a change of lattice constant in less than 1.14%percent1.141.14\%1.14 %, and a drift of lattice position in less than 1.61%percent1.611.61\%1.61 %.

1 Introduction

Atoms trapped in different types of potentials are central subjects under investigation in the field of atomic physics and quantum optics. Due to the exceptional high purity and controllability[1], they have also received extensive attention in quantum simulation[2], quantum sensing[3] and quantum computation[4]. The pre-cooled cold atoms are usually trapped by conservative potentials such as magnetic and optical traps. To control the interaction of atoms for quantum simulation, or to isolate them for metrology applications, trapping potentials with different geometric configurations are implemented to prepare bulk gases in optical dipole traps or spatially ordered systems in optical lattices [5]or tweezer arrays. For example, optical lattices, known as spatially repeated optical potentials induced by position-dependent light shift of the inference of coherent laser beams, are widely used to control the spatial positions, the local density and even the interaction of the confined atoms, and trigger a burst of research in simulating important lattice models of various kinds in the last two decades [6, 7, 8].

The key mechanism of generating a periodic optical potential is using the interference pattern of coherent laser beams to exert periodic conservative forces to the atoms. The configuration of the resulting lattice is determined by the spatial arrangement, the polarization, and the relative phase of the laser beams. For instance, six counter-propagated beams from orthogonal directions can form a simple cubic lattice, in which pioneering experiments of superfluid to Mott insulator transition was first demonstrated [9]. Beyond that, more complex non-cubic lattice configurations, such as triangular lattice [10, 11], honeycomb/hexagonal lattice [12, 13, 14], Kagome lattice [15, 16], superlattice [17, 18, 19] and even quasi-crystalline lattice [20, 21] are also implemented in various experiments to explore the exotic quantum phases of matter.

Meanwhile, the splitting, controlling and overlapping multiple laser beams to form a complex optical lattice usually requires a complicated optical setup and faces some technical difficulties due to its low stability. For example, the geometry of a non-cubic lattice is very sensitive to the relative phases of the laser beams from non-orthogonal directions. A change of phase can deform the geometric configuration and even change the topology of the lattice. To overcome this issue, significant effort needs to be devoted to reduce the relative phase fluctuation by a phase-locking system [22] or even by an additional optical interferometry [23]. Instead of multiple-beam interference in a complicated optical system with a long optical path and delicate feedback control for their phases, multi-facet prism was proposed to create triangular and square lattices[24]. And a scheme of interfering multiple parallel beams focused by the same lens was reported to generate a dynamic cubic lattice [25]. These approaches naturally lead to more stable lattices since they are mostly composed by simple and fixed optical parts, which can significantly reduce relative phase fluctuation of the beams. So far, this scheme is yet to be implemented for more complex optical lattices.

In this work, we report a new scheme to generate complex non-cubic lattices, such as triangular lattice and quasi-crystalline lattice by shining a single Gaussian beam through a multi-facet azimuthal-symmetric prism. Every portion of the prism acts as a prism and deflects the corresponding part of the Gaussian beam with the same deflecting angle. All deflected parts of the beam meet in the Bessel region and interfere to form a lattice. The relative phases of different parts are determined by their corresponding optical paths before they meet. Since all the interfering beams are from the same Gaussian beam and almost share the same optical path, the consequent relative phases are small and stable. The interference pattern is measured in the Bessel region and also at the atom position by an imaging system with a high-resolution objective. We demonstrate a triangular lattice and a ten-fold quasi-crystalline lattice, which can be obtained and easily switched by just replacing the prism. The lattice patterns and their Fourier transform are analyzed to confirm their periodic characteristics. The time evolution of the interference pattern is recorded to show the stability of the beam, which also includes the pointing stability of the lattice beam itself. From the measurement, the root-mean-squared error of lattice spacing fluctuation is less than 1.14%percent1.141.14\%1.14 %, and a drift of lattice position in less than 1.61%percent1.611.61\%1.61 % for a measuring time up to 200 minutes. Moreover, the radial envelope of the lattice depth is no-longer harmonic but flat-top like, which helps to reduce effects of inhomogeneous density and phase separation of the trapped atoms. This phase-locking-free single beam optical lattice provides an efficient and compact way to access complex optical lattice for quantum simulation and metrology applications.

2 Concept and design

Refer to caption
Figure 1: (a) Illustration of setup to generate and picture optical lattices. The lattice potential is realized on plane X by shining a single beams through a prism, and zoomed to a desired size by a telescope ( F-75 lens and Mitutoyo objective). To further visualize the lattices, a second telescope (Olympus objective and F-250 lens) is used. (b) Schematic diagram of the multi-facet prism and light interference.

An overview of our phase-locking-free single beam optical lattice scheme is illustrated in Fig. 1(a). A collimated Gaussian beam is shined through a multi-facet prism, and is split into multiple beams. The multi-facet prism acquires an n𝑛nitalic_n-fold rotational symmetry, with symmetry axis along the optical axis. All beams are deflected toward the optical axis and overlap in the Bessel region, where they interference and form an optical lattice of desired configuration. As shown in Fig. 1(b), the Gaussian beam with electric field E0ei𝒌𝒓subscript𝐸0superscript𝑒𝑖𝒌𝒓{\vec{E}}_{0}e^{i\boldsymbol{k}\cdot\boldsymbol{r}}over→ start_ARG italic_E end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i bold_italic_k ⋅ bold_italic_r end_POSTSUPERSCRIPT incidents from the flat side of the prism, where the wavevector |k|=2π/λ𝑘2𝜋𝜆\left|k\right|=2\pi/\lambda| italic_k | = 2 italic_π / italic_λ and the beam waist is w0subscript𝑤0w_{0}italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. If the transmitted loss is negligible, it is then split into multiple beams with electric field components of E0iei𝒌i𝒓+δisubscript𝐸0𝑖superscript𝑒𝑖subscript𝒌𝑖𝒓subscript𝛿𝑖{\vec{E}}_{0i}e^{i\boldsymbol{k}_{i}\cdot\boldsymbol{r}+\delta_{i}}over→ start_ARG italic_E end_ARG start_POSTSUBSCRIPT 0 italic_i end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i bold_italic_k start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ⋅ bold_italic_r + italic_δ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUPERSCRIPT with i=1,2,,n𝑖12𝑛i=1,2,\dots,nitalic_i = 1 , 2 , … , italic_n. The wave vector and intensity of each scattered beam are assumed to be identical, i.e., |𝒌i||𝒌|subscript𝒌𝑖𝒌\left|\boldsymbol{k}_{i}\right|\approx\left|\boldsymbol{k}\right|| bold_italic_k start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | ≈ | bold_italic_k | and E0i=E0/nsubscript𝐸0𝑖subscript𝐸0𝑛{\vec{E}}_{0i}={\vec{E}}_{0}/\sqrt{n}over→ start_ARG italic_E end_ARG start_POSTSUBSCRIPT 0 italic_i end_POSTSUBSCRIPT = over→ start_ARG italic_E end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT / square-root start_ARG italic_n end_ARG. The n𝑛nitalic_n beams are deflected towards the axis with a small angle θ=(μ1)α𝜃𝜇1𝛼\theta=\left(\mu-1\right)\alphaitalic_θ = ( italic_μ - 1 ) italic_α, where μ𝜇\muitalic_μ is the refractive index of the fused silica prism. Assuming all the n𝑛nitalic_n beams have the same polarization, one can calculate the interference pattern after passing through the prism. Ignoring the radial variation of light intensity during light propagation, the distribution of interference light intensity in the radial plane is given by

I𝐼\displaystyle Iitalic_I =\displaystyle== |E0|2n[j=1nei(𝒌j𝒓+δj)]2superscriptsubscript𝐸02𝑛superscriptdelimited-[]superscriptsubscript𝑗1𝑛superscript𝑒𝑖subscript𝒌𝑗𝒓subscript𝛿𝑗2\displaystyle\frac{\left|E_{0}\right|^{2}}{n}\left[{\sum\limits_{j=1}^{n}e^{i{% (\boldsymbol{k}_{j}\cdot\boldsymbol{r}+\delta_{j}})}}\right]^{2}divide start_ARG | italic_E start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_n end_ARG [ ∑ start_POSTSUBSCRIPT italic_j = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_i ( bold_italic_k start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ⋅ bold_italic_r + italic_δ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (1)
=\displaystyle== |E0|2+2|E0|2nm>jnj=1mcos[(kρρcos(ϕϕj)cos(ϕϕm))+δjδm],superscriptsubscript𝐸022superscriptsubscript𝐸02𝑛superscriptsubscript𝑚𝑗𝑛superscriptsubscript𝑗1𝑚subscript𝑘𝜌𝜌italic-ϕsubscriptitalic-ϕ𝑗italic-ϕsubscriptitalic-ϕ𝑚subscript𝛿𝑗subscript𝛿𝑚\displaystyle\left|E_{0}\right|^{2}+\frac{2\left|E_{0}\right|^{2}}{n}{\sum% \limits_{m>j}^{n}{\sum\limits_{j=1}^{m}{\cos\left[\left(k_{\rho}\rho\cos\left(% \phi-\phi_{j}\right)-\cos\left(\phi-\phi_{m}\right)\right)+\delta_{j}-\delta_{% m}\right]}}},| italic_E start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG 2 | italic_E start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_n end_ARG ∑ start_POSTSUBSCRIPT italic_m > italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_j = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT roman_cos [ ( italic_k start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT italic_ρ roman_cos ( italic_ϕ - italic_ϕ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) - roman_cos ( italic_ϕ - italic_ϕ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT ) ) + italic_δ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT - italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT ] ,

where ρ𝜌\rhoitalic_ρ and kρsubscript𝑘𝜌k_{\rho}italic_k start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT are the spatial coordinate and projected wavevector in the radial plane, respectively, ϕitalic-ϕ\phiitalic_ϕ is the azimuthal angle, and ϕj=2π(j1)/nsubscriptitalic-ϕ𝑗2𝜋𝑗1𝑛\phi_{j}=2\pi(j-1)/nitalic_ϕ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT = 2 italic_π ( italic_j - 1 ) / italic_n. One can easily conclude that the interference pattern also acquires an n𝑛nitalic_n-fold rotational symmetry. The maximal overlapping area can be reached on plane X, located at a distance of w0/tanθsubscript𝑤0𝜃w_{0}/\tan\thetaitalic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT / roman_tan italic_θ from the vertex of the prism. For example, one can obtain a triangular lattice for n=3𝑛3n=3italic_n = 3, a cubic lattice for n=4𝑛4n=4italic_n = 4, and a ten-fold quasicrystalline lattice for n=5𝑛5n=5italic_n = 5. The cases for even higher value of n𝑛nitalic_n are not in the scope of this work.

3 Experiment

Next we focus on the triangular lattice and ten-fold quasicrystalline lattice, and demonstrate the implementation of our scheme and analyze the lattice stability. For that purpose, the prisms used are the 3- and 5-fold rotationally symmetric prisms. A 532nm laser beam (Coherent Verdi V18, coherent length  60m) is fiber coupled to the test system and collimated with a beam waist of w0=1.8mmsubscript𝑤01.8𝑚𝑚w_{0}=1.8mmitalic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1.8 italic_m italic_m. Fused silica prisms with a refractive index of μ=1.46𝜇1.46\mu=1.46italic_μ = 1.46 (for 532nm laser) are chosen for high power applications due to their lower lensing effect. The angle of the prisms is α=3𝛼superscript3\alpha=3^{\circ}italic_α = 3 start_POSTSUPERSCRIPT ∘ end_POSTSUPERSCRIPT, so the deflection angle is θ1.4𝜃superscript1.4\theta\approx 1.4^{\circ}italic_θ ≈ 1.4 start_POSTSUPERSCRIPT ∘ end_POSTSUPERSCRIPT.

A camera (IMAGINGSOURCE, DMK 23UX236) is used to record the light intensity distribution of the interference pattern on plane X. The pixel size of the camera is 2.8μm2.8𝜇𝑚2.8\mu m2.8 italic_μ italic_m. By utilizing 3- and 5-fold rotationally symmetric prisms, we obtain a triangular lattice and a ten-fold quasicrystalline lattice, as shown in Fig. 2. Additionally, a camera with smaller pixel size (Raspberry Pi, module V2) is also used to obtain higher spatial resolution of light intensity distribution for stability measurements. To obtain the lattice constant of the triangular lattice, two typical neighboring lattice sites are randomly chosen. As shown in Fig. 2(a), by Gaussian fitting the light intensity at each site and averaging distances crossing tens of nearest neighboring peaks, we obtain a lattice constant of 14.9μm14.9𝜇𝑚14.9\mu m14.9 italic_μ italic_m, which agrees well with the value of 14.7μm14.7𝜇𝑚14.7\mu m14.7 italic_μ italic_m in the calculated pattern, considering the possible machining precision of this prism. For the ten-fold quasicrystalline lattice demonstrated in Fig. 2(c), the light intensity distribution also consists with the calculated pattern. The distance between the two chosen neighboring lattice sites is 23.0μm23.0𝜇𝑚23.0\mu m23.0 italic_μ italic_m. Some defects in the ten-fold quasicrystalline lattice are attributed to the imperfection machining of the prism.

In order to further confirm the symmetry of the lattices, Fourier transformation of the light intensity distribution of both lattices are implemented. The six-fold rotational symmetric pattern of the Fourier transformed triangular lattice is clearly witnessed in Fig. 2(b). In the case of ten-fold quasi-crystalline lattice, the Fourier transformation shows a ten-fold rotational symmetry pattern in Fig. 2(d), which further proves the rotational symmetry of the quasi-crystalline lattice formed by the interference of five deflected beams after the 5-fold prism.

Refer to caption
Figure 2: The real-space light-intensity distribution (a) and its Fourier transform (b) of a triangular lattice. The lattice constant (red arrow) is 14.9μ𝜇\muitalic_μm. The real-space distribution (c) and its Fourier transform (d) of a ten-fold quasi-crystalline lattice. The distance between the two neighboring sites labeled by red arrow is 23.0μ𝜇\muitalic_μm.

Further, we find that the envelop of the intensity distribution of the optical lattices established on plane X is no longer a Gaussian shape, but a flat-top like function. This feature is desirable since it can partially overcome the long-standing problem of inhomogeneity, which is usually inevitable in traditional lattices formed by the interference of counter-propagating Gaussian beams along different directions. In our scheme, taking a Gaussian beam through a double-facet prism as an example for simplicity, the upper half of the cut Gaussian beam has lower intensity at the top and higher intensity at the bottom when it is deflected towards the optical axis. On the opposite, the lower half of the cut beam is stronger at the top and weaker at the bottom. In this way, the two beams compensate each other in their Gaussian intensity, and result in a more uniform overall light intensity distribution. The inhomogeneity can also be reduced by flatening the beam intensity distribution with a spatial light modulator (SLM), but with a much higher cost and technical complexity [26].

3.1 Lattice projection

To have an optical lattice being useful for quantum simulation in a cold atom experiment, the lattice constant needs to be significantly reduced to sub-micron scale so that atoms can hop between lattice sites and present many-body characteristics. To this aim, as shown in Fig. 1, a telescope with a lens of f1=75subscript𝑓175f_{1}=75italic_f start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 75mm and a commercial high NA microscope objective (Mitutoyo, G Plan Apo 50X, NA = 0.5, f2=4subscript𝑓24f_{2}=4italic_f start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = 4mm) are used to project the lattice potential to the plane where atom are supposed to be. The lattice constant in the atom plane is reduced by a factor of 18.75, leading to a lattice constant of 0.795μ0.795𝜇0.795\mu0.795 italic_μm for the triangular lattice, and a distance of 1.23μ1.23𝜇1.23\mu1.23 italic_μm between the two neighboring sites (labeled by red arrow in Fig. 2) for the ten-fold quasicrystalline lattice. We can also estimate the lattice depth from the power and beam waist of the original 532nm Gaussian beam. Taking an experiment with 6Li as an example, a Gaussian beam with power of 10W can implement a 100×100 triangular lattice with lattice constant of 0.795μ0.795𝜇0.795\mu0.795 italic_μm and a lattice depth up to 20Er20subscript𝐸𝑟20E_{r}20 italic_E start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT, which is high enough to observe many interesting quantum phenomena.

Refer to caption
Figure 3: Images of (a) the triangular lattice and (b) the ten-fold quasicrystalline lattice after magnification of the second telescope.

To further confirm the detailed structure of the projected optical lattice, we use a second telescope with another microscope objective (Olympus, MPlanFL N 50X, NA = 0.8, f3=3.6subscript𝑓33.6f_{3}=3.6italic_f start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT = 3.6mm) and a lens with a focal length of f4=250subscript𝑓4250f_{4}=250italic_f start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT = 250mm. The lattice images are magnified by a factor of 69.4, and are shown in Fig. 3. With the same Gaussian fitting of light intensity distribution as above, a period of 55.7μ55.7𝜇55.7\mu55.7 italic_μm is obtained from the magnified image of the triangular lattice, which gives a consistent result of 0.802μ0.802𝜇0.802\mu0.802 italic_μm for the lattice constant before magnification. The first microscope objective is infinity calibrated, so that the f1=75subscript𝑓175f_{1}=75italic_f start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 75mm lens can be replaced and the lattice constant can be freely chosen for different experimental requirements.

3.2 Lattice stability

Non-cubic lattices are usually sensitive to many factors, such as the pointing stability, the relative phase and the polarization of laser beams. In our setup, the total path length is about 200mm, which is about a factor of 10 shorter than that in typical experiments. Additionally, since all different parts of the beam from the multi-facet prism are deflected from the same Gaussian beam, the polarization are almost the same with a small-angled deflection and the relative phases barely change. Those reasons inherently lead to a super-stable lattice without any phase-locking system used in previous experiments [10, 12, 15].

Refer to caption
Figure 4: Stability of the lattice potential. (a) Time variation of the lattice constant of the triangular lattice, measured at interval Δt=Δ𝑡absent\Delta t=roman_Δ italic_t =40s, 5s, and 0.6s. (b) Time variation of distance between the two neighboring sites of the quasicrystalline lattice. The error bars are calculated from the 95%percent9595\%95 % confidence interval of Gaussian fits when we estimate the center of the interference peaks.

To evaluate other possible factors such as mechanical vibration and airflow disturbance that may bring in instability, we use a camera to monitor the variation of the optical lattice over a certain period of time. To eliminate possible problems of the adjustable mounts and translation stages of the objectives, a simple measurement of the lattice stability is taken by monitoring the lattice on plane X behind the prism. The measurement is repeated for every 0.6s, 5s and 40s, respectively. The variation of distances of peaks in the lattice pictures shown in Figs. 2(a) and 2(c) might be attributed to the noise of relative phases and polarizations. Fig.4(a) shows the variation of lattice constant of the triangular lattice. For short experiment cycles, the root mean square error (RMSE) of the lattice constant is 0.14μ0.14𝜇0.14\mu0.14 italic_μm, while it slightly grows to 0.17μ0.17𝜇0.17\mu0.17 italic_μm for long experiment cycles. For quasi-crystalline lattice, the RMSE of the characteristic lattice length for short experiment cycles is less than 0.13μ0.13𝜇0.13\mu0.13 italic_μm, and is 0.21μ0.21𝜇0.21\mu0.21 italic_μm for long experiment cycles, as shown in Fig. 4(b). The RMSE is less than 1.14%percent1.141.14\%1.14 % of the lattice constant for the triangular lattice, and is less than 0.91%percent0.910.91\%0.91 % of the distance between the two neighboring sites for the quasi-crystalline lattice. These results clearly demonstrate the stability of optical lattices obtained by our scheme. We also compare the relative position of the two kinds of lattice on plane X, and the drifts of lattice positions are within 0.24μm0.24𝜇𝑚0.24\mu m0.24 italic_μ italic_m for most of the time. The drift is about 1.61%percent1.611.61\%1.61 % of the triangular lattice constant, and 1.04%percent1.041.04\%1.04 % of the quasi-crystalline lattice site distance.

4 Summary and outlook

We experimentally demonstrate a scheme to realize optical lattices of different structures by shining a single beam to multi-facet prisms, and interfering beams diffracted from different facets. Benefits from the simple optical path without any mechanical moving components and phase-locking system, the scheme features ultra-high stability against fluctuations of polarization, optical intensity, and relative phase. For demonstrative purpose, we show a triangular lattice and a ten-fold symmetric quasi-crystalline lattice by using 3-fold and 5-fold symmetric prisms, respectively. In both case, the spacing between lattice sites can be demagnified down to a few hundreds of micrometers by an optical telescope, and the lattice depth can reach similar-to\sim20Ersubscript𝐸𝑟E_{r}italic_E start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT with a reasonable laser power of similar-to\sim10W. To quantify the stability of the optical lattices, we monitor the time variation of lattice spacing, and find a root-mean-squared error of lattice spacing less than 1.14%percent1.141.14\%1.14 % and a drift of lattice position less than 1.61%percent1.611.61\%1.61 % for a period of up to 200 minutes. The experimental scheme is very simple, low-cost and easy to be implemented for other configurations and applications. This single-beam lattice provides a super stable lattice potential for further exploration of quantum simulation of exotic phases of quantum matter in complicated lattice system, and might also find more applications in optical tweezers and optical lattice clocks.

\bmsection

Funding The National Key Research and Development Program of China (Grant No. 2022YFA1405301); The National Natural Science Foundation of China (Grants Nos.12274460, 12074428, and 92265208).

References

  • [1] I. Bloch, J. Dalibard, and S. Nascimbène, “Quantum simulations with ultracold quantum gases,” \JournalTitleNature Physics 8, 267–276 (2012).
  • [2] I. M. Georgescu, S. Ashhab, and F. Nori, “Quantum simulation,” \JournalTitleRev. Mod. Phys. 86, 153–185 (2014).
  • [3] C. L. Degen, F. Reinhard, and P. Cappellaro, “Quantum sensing,” \JournalTitleReviews of Modern Physics 89, 035002 (2017).
  • [4] D. S. Weiss and M. Saffman, “Quantum computing with neutral atoms,” \JournalTitlePhysics Today 70, 44–50 (2017).
  • [5] D. Jaksch, C. Bruder, J. I. Cirac, et al., “Cold bosonic atoms in optical lattices,” \JournalTitlePhysical Review Letters 81, 3108–3111 (1998).
  • [6] I. Bloch, J. Dalibard, and W. Zwerger, “Many-body physics with ultracold gases,” \JournalTitleReviews of Modern Physics 80, 885–964 (2008).
  • [7] F. Schäfer, T. Fukuhara, S. Sugawa, et al., “Tools for quantum simulation with ultracold atoms in optical lattices,” \JournalTitleNature Reviews Physics 2, 411–425 (2020).
  • [8] A. Browaeys and T. Lahaye, “Many-body physics with individually controlled rydberg atoms,” \JournalTitleNature Physics 16, 132–142 (2020).
  • [9] M. Greiner, O. Mandel, T. Esslinger, et al., “Quantum phase transition from a superfluid to a mott insulator in a gas of ultracold atoms,” \JournalTitleNature 415, 39–44 (2002).
  • [10] C. Becker, P. Soltan-Panahi, J. Kronjäger, et al., “Ultracold quantum gases in triangular optical lattices,” \JournalTitleNew Journal of Physics 12, 065025 (2010).
  • [11] J. Yang, L. Liu, J. Mongkolkiattichai, and P. Schauss, “Site-resolved imaging of ultracold fermions in a triangular-lattice quantum gas microscope,” \JournalTitlePRX Quantum 2, 020344 (2021).
  • [12] L. Tarruell, D. Greif, T. Uehlinger, et al., “Creating, moving and merging dirac points with a fermi gas in a tunable honeycomb lattice,” \JournalTitleNature 483, 302–305 (2012).
  • [13] T. Uehlinger, G. Jotzu, M. Messer, et al., “Artificial graphene with tunable interactions,” \JournalTitlePhysical Review Letters 111, 185307 (2013).
  • [14] G. Jotzu, M. Messer, R. Desbuquois, et al., “Experimental realization of the topological haldane model with ultracold fermions,” \JournalTitleNature 515, 237–240 (2014).
  • [15] G.-B. Jo, J. Guzman, C. K. Thomas, et al., “Ultracold atoms in a tunable optical kagome lattice,” \JournalTitlePhysical Review Letters 108, 045305 (2012).
  • [16] T.-H. Leung, M. N. Schwarz, S.-W. Chang, et al., “Interaction-enhanced group velocity of bosons in the flat band of an optical kagome lattice,” \JournalTitlePhysical Review Letters 125, 133001 (2020).
  • [17] S. Fölling, S. Trotzky, P. Cheinet, et al., “Direct observation of second-order atom tunnelling,” \JournalTitleNature 448, 1029–1032 (2007).
  • [18] M. Lohse, C. Schweizer, O. Zilberberg, et al., “A thouless quantum pump with ultracold bosonic atoms in an optical superlattice,” \JournalTitleNature Physics 12, 350–354 (2016).
  • [19] S. Nakajima, T. Tomita, S. Taie, et al., “Topological thouless pumping of ultracold fermions,” \JournalTitleNature Physics 12, 296–300 (2016).
  • [20] K. Viebahn, M. Sbroscia, E. Carter, et al., “Matter-wave diffraction from a quasicrystalline optical lattice,” \JournalTitlePhysical Review Letters 122, 110404 (2019).
  • [21] J.-C. Yu, S. Bhave, L. Reeve, et al., “Observing the two-dimensional bose glass in an optical quasicrystal,” \JournalTitleNature 633, 338–343 (2024).
  • [22] C. Becker, “Multi component bose-einstein condensates: from mean field physics to strong correlations,” Ph.D. thesis. (University of Hamburg, Hamburg, 2009).
  • [23] C. K. Thomas, “Quantum simulation of triangular, honeycomb and kagome crystal structures using ultracold atoms in lattices of laser light,” Ph.d.thesis. (University of California, Berkeley, 2017).
  • [24] M. Lei, B. Yao, and R. A. Rupp, “Structuring by multi-beam interference using symmetric pyramids,” \JournalTitleOptics Express 14, 5803–5811 (2006).
  • [25] R. A. Williams, J. D. Pillet, S. Al-Assam, et al., “Dynamic optical lattices: two-dimensional rotating and accordion lattices for ultracold atoms,” \JournalTitleOptics Express 16, 16977–16983 (2008).
  • [26] H.-J. Shao, Y.-X. Wang, D.-Z. Zhu, et al., “Antiferromagnetic phase transition in a 3d fermionic hubbard model,” \JournalTitleNature 632, 267–272 (2024).