Advances in Optical TechnologiesVolume 2012 (2012), Article ID 647657, 9 pagesdoi:10.1155/2012/647657

Research Article

Generation of Atomic Optical Lattices by Dammann Gratings

Xianming Ji,¹ Shuwu Xu,1,2 and Songbo Gu¹

¹School of Science, Nantong University, Nantong 226007, China²State Key Laboratory of Precision Spectroscopy, Department of Physics, East China Normal University, Shanghai 200062, China

Received 10 October 2011; Accepted 5 February 2012

Academic Editor: Pierre Chavel

Copyright © 2012 Xianming Ji et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We first calculated the diffraction intensity distributions of the Dammann gratings illuminated by Gaussian light wave. The empirical equations were deduced by numerical calculations to calculate the parameters, such as the spatial period, the maximum intensity, and the maximum intensity gradient, of the optical trap array composed by a set of Dammann gratings and a focus lens. Thus, a novel type of optical trap array for trapping cold atoms (or molecules) was proposed and its features were discussed. The results showed the optical trap array with very short period could be generated. High optical dipole potential could be presented so as to have strong attractive force to the atoms to form atomic optical lattices of high lattice density. Compared with the optical lattices formed by standing wave interferences of CO₂ laser, there are many unique advantages of which are formed by Dammann gratings.

1. Introduction

In early 1970s, the concept of “Dammann grating” was first proposed by Dammann and Görtler [1] when they studied the copy of multiple imaging. The Dammann grating is a kind of phase-type grating with several groups of phase jump points in each period. The positions of the phase jump points can be computed by optimal designing. When the Dammann grating is illuminated by a coherent monochromatic light wave, many beams of diffraction light wave with equal intensity can be generated without the nonuniformity of beam intensity caused by sinc functions in common gratings. Many types [2–5] of Dammann gratings with high efficiency and large splitting ratio are developed in several decades to satisfy their applications. The orthogonal binary-phase gratings completing 64 × 64 splitting ratio have been successfully fabricated. Due to their advantages, such as the simplicity to be designed [6], the mature production process, and the easy calibration, they have been widely paid attentions to. In 1980s, Dammann grating was applied to the field of optical computing and optical interconnection. In 1990s, Tooley realized the image process of optical cells by using Dammann grating [7].

In recent years, a series of important developments and fruitful results have been obtained both in theory and in experiment on the cooling, trapping, and manipulating the neutral atoms (or molecules). As the deeply developing of the quantum information process and optical lattices [8–10], the surface microtrap array for trapping the cold atoms (or molecules) has become the research frontier in the field of atomic, molecular, and optical physics. Schemes of generating surface magnetic lattices, magneto-optical lattices, and Bose-Einstein condensation (BEC) array have been proposed [11–13]. Birkl and his coworkers realized all-optical array of trapping neutral cold atoms by using micro-trap array generated by a microlens array [14, 15]. The all-optical array has many advantages over the magnetotrap array. However, many disadvantages exist in the method of using the microlens array to generate the micro-trap array. Producing a large number of the sublens array requires high techniques. For example, the deficiency of the production of the micro-lens array and the nonuniformity of the illuminating light beam will directly affect the uniformity of the intensity distribution of the optical trap array, and additionally, the accompanying stray light will heat the atoms (or molecules) trapped. The difficulty can be effectively overcome by replacing the micro-lens array with a Dammann grating and a focusing Fourier lens to form the optical trap array. In Section 2, we calculated the intensity distributions of the optical trap array generated by Dammann grating. The disadvantages about using Dammann grating to form the optical trap array were analyzed, and thus, the scheme of generating the optical trap array with high density of the optical traps was proposed in Section 3. In Section 4, we discussed the potential applications of this optical trap array in atomic and molecular optics. And a summary was stated in the final section.

2. The Generation and the Intensity Distribution Computation of the Optical Trap Array with Dammann Grating

We only discuss the orthogonal phase-type Dammann grating with the phase value of 0 or in the present work. The splitting ratio of the grating is , and the period is . The coordinates of the phase jump points of 1D grating in each period are , . Here, is the number of the phase jump points. The phase distribution is shown in Figure 1. The amplitude of the corresponding transmittance function can be written as where rect(·) is the rectangular function [16]. To the phase-type Dammann gratings with the splitting ratio of 15 × 15 and 25 × 25, the values [17] of and are listed in Table 1 (the period of the grating is normalized to ). The transmittance function of the 2D Dammann grating can be written as

Table 1: Coordinates of the phase jump points of two Dammann gratings.

.

15
4
0, 0.14370; 0.26763, 0.61413; 0.70511, 0.78188; 0.87009, 0.91972.

25
7
0, 0.15387; 0.16452, 0.27298; 0.33211, 0.38120; 0.44909, 0.52349; 0.64396, 0.76823; 0.82082, 0.85980; 0.91639, 0.95832.

Figure 1: 1D phase distribution of Dammann grating.

Assuming the number of the period of the grating is , the diffraction light is focused by a lens with a focal length , the lens and the grating are superposed on (z = 0) plane, and the light axis of the lens is along axis, as is shown in Figure 2. The grating is perpendicularly illuminated by TEM₀₀ mode of a Gaussian light wave, and the function of the incident light is where P is the output power of the laser, k = 2π/, and is the wavelength of the laser. When , , where is the waist radius of the light, R(0)→∞, and ζ(0) = 0. According to the Fresnel diffraction theory, the complex amplitude of light field on the xoy plane, which is located in the back of the lens and has a distance to the lens, is given by (omitting the constant phase factors) where and are integers, , . and , the phase propagation factor of the Fresnel diffraction and the phase transmittance function of the lens, respectively, are as

Figure 2: Light path of the optical trap array employing a set of a Dammann grating and a confus lens.

Then, the intensity distribution on the back of lens is . On the focal plane (), the complex amplitude can be written as

It is shown that () is the Fourier transform of the product of the transmission function of the Dammann grating and the Gaussian wave function if the quadratic phase factor before the sum operators is omitted. Calculation shows that the intensity distribution on the focal plane forms mainly number of principal maximum bright spots, which we call as the optical trap array and the optical trap array centers at (, ), where and are integers with , . The period of the optical trap array is given by 　　　　　　　　　　　　

Figure 3 shows the density distribution of the intensity on the focal plane of a 15 × 15 optical trap array. The spatial intensity distributions of each optical trap in the array are basically identical. We numerically calculated the parameters that characterize the intensity distribution of the single optical trap. The results show that the intensity of the optical trap is related with several parameters of the optical system, such as the width of the grating (NT), the focal length of the lens , the laser power , and the waist radius of the Gaussian light wave. To a given laser power, when the waist radius is set to () , each optical trap of the array will reach their maximum intensity simultaneously. Additionally, the intensity of each optical trap has a good uniformity when the number of the grating periosd satisfies . To Dammann gratings with the splitting ratio of 15 and 25, taking , we numerically calculated many groups of the optical systems with different parameters and got an empirical formulae. The maximum intensity of the optical trap can be written as where is the effective relative aperture of the lens. Take 1/ of the maximum intensity as the boundary, the transverse and longitudinal widths of the optical trap are given by respectively. Regarding the optical trap as a rotating ellipsoid, the trapping volume can be written as 　

Figure 3: Intensity pattern of a 15 × 15 optical trap array.

The transverse and longitudinal gradients of the optical trap are given by respectively.

Equations (8)–(13) show that demanded optical trap arrays can be obtained by varying the parameters of the optical system, including relative aperture of the lens , the period of the grating , the number of the period , the splitting ratio , and the laser power . If the power and wavelength of the illuminating light are fixed, the intensity and its gradient of the optical trap increase when we properly increase the relative aperture, while the trapping volume increases when we decrease . Note that there is a certain degree of error in (9)–(13), which is relative to the number of the period and the relative aperture . When is large while is small, the error will be little. Equations (9)–(13) are valid when and . Furthermore, the maximum intensity of each optical trap in the array is not exactly identical. A factor Uni is defined to scale the maximum intensity uniformity of the optical traps:

Table 2 lists the array of optical trap formed by choosing different parameters of the optical system. Here, the intensity distribution parameters listed in rows of the 2nd to 5th are the results numerically calculated after (7). It is coincident with (8)–(13).

Table 2: Optical trap arrays generated by different optical system ().

Illuminating light parameters
(Wm⁻²)
(Wm⁻³)
(Wm⁻³)
(μm)
(μm)
(μm)
(Cm³)
Uni (%)
Trapped atoms

β
N
K
T (mm)
P (W)
λ (μm)

Atoms
(μK)
(s ⁻¹)

0.05
20
25
2.0
20
1.06
1.55 × 10⁸
2.12 × 10¹³
3.68 × 10¹¹
212
17.6
1020
1.65 × 10⁻⁷
0.207
¹³²Xe
45
0.313

0.1
10
25
8.0
20
1.06
6.19 × 10⁸
1.70 × 10¹⁴
5.89 × 10¹²
53
8.81
253
1.03 × 10⁻⁸
0.45
¹³³Cs
147
0.83

0.05
20
15
2.0
1
0.78
3.96 × 10⁷
7.40 × 10¹²
1.28 × 10¹¹
156
12.95
749
6.59 × 10⁻⁸
0.32
39 K
75.4
8.62

0.1
3
15
26.67
2
0.78
3.22 × 10⁸
11.84 × 10¹³
4.08 × 10¹²
39
6.47
187
4.11 × 10⁻⁹
2.2
⁷Li
51.1
0.434

0.2
4
15
40
10
0.53
1.34 × 10¹⁰
1.45 × 10¹⁶
1.02 × 10¹⁵
5.3
2.21
32.0
8.23 × 10⁻¹¹
5.6
⁸⁸Sr
340
1.65

0.2
3
15
53
4
0.53
5.4 × 10⁹
5.84 × 10¹⁵
4.12 × 10¹⁴
3.97
2.21
32.0
8.23 × 10⁻¹¹
8.7
⁴⁰Ca
503
8.96

Taking T = 2.0 mm, , f = 400 mm, = 20 mm, λ = 1.06 μm, and  W, which correspond to the 2nd row in Table 2, the intensity distribution of a single optical trap in the array is shown in Figure 4. Figure 4(a) gives the intensity distribution contours on the focal plane, and Figure 4(b) shows the corresponding 2D intensity distribution. The maximum intensity of the optical trap is  Wm⁻². The intensity distribution on the xoy plane is close to Gaussian distribution, and the intensity gradient along x (or y) direction is reached to a magnitude of 1013 Wm⁻³. Also, the intensity distribution along direction is close to Gaussian distribution. The corresponding volume of the optical trap is about 10⁻⁷ cm³. The spatial period of the optical array is about 0.2 mm.

(a)

(b)

Figure 4: Intensity distribution of a single optical trap in the array employing a set of a Dammann grating and a focus lens (, , , ,   μm). (a) Intensity isoline, (b) 2D intensity distribution.

Through analyzing the geometrical parameters of the optical trap , , , and , we can find that each optical trap in the array generated by Dammnn grating is quite similar with that generated by a single square lens. It implies that the combination of Dammann grating and a lens is equivalent to the lens array composed by number of microlens. However, generating the optical trap array by Dammnn grating has its unique advantages. The assembling is very simple; only one lens is needed to obtain the number of optical traps. Furthermore, the intensity distributions of the optical trap in the array are basically uniform when the intensity distribution of the illuminating light is inhomogeneous as the Gaussian light wave, while if adopting micro-lens array, the parameters of each micro-lens and the light wave illuminating each lens are required to be uniform.

3. Generation of High-Density Optical Trap Array by Dammann Grating

It is well known that Dammann grating is an optical element of far field diffraction. Calculating intensity distribution of diffraction by using (4) requires that the relative aperture, , of the lens should be small. It can be known from (8) that the distance between two optical traps in the array generated by Dammann grating is . Clearly, N, the number of the period of the grating, can not be small. So, when the wavelength of the illuminating light wave is not very short or the diffraction light is focused by a common spherical lens, the period of the optical trap evidently should be large. This will result in low optical trap density number in the array. In practical applications, however, we need high-density optical trap arrays. For example, when employing the array to form an optical lattice to trap cold atoms (molecules), we want the optical lattice constant, , to be small. It is well known that the wavelength of the illuminating light should be longer than the resonant wavelength of the cold atoms when the atoms are trapped at strong-field-seeking state. Consequently, large relative aperture of the focal lens is required to generate high-density optical trap array by Dammann grating. But the approximate condition of Fresnel diffraction is not valid any more in this case; that is, (5) or (7) or their deduced (8)–(13) are not valid any more. Sommerfeld diffraction theory is now needed to calculate the intensity distribution. Equation (5) should be corrected as and term in (4) now should also be replaced as . If the phase transmittance function of the lens is still as described in (6), there would be large error of the diffraction intensity calculated by Fresnel diffraction theory when , while identical intensity and identical spaced optical array cannot be obtained through calculations by Sommerfeld diffraction theory.

When the splitting ratio of the grating is not very large and the distance between two adjacent optical traps is small, the diffraction intensity mainly distributes near the focal point. According to aplanatic theory, (6) can be corrected as where and are two constants that can be determined by and . Choosing λ = 0.53 μm, which can be obtained by frequency doubling of YAG laser with wavelength of 1.06 μm, numerical calculation shows that, when , we can get that = 0.9990 and = 0.9299. Substituting (14) and (16) into (7) and calculating the diffraction intensity distribution focused by lens with large aperture show that relatively uniform optical trap array can still be obtained. Taking , two kinds of optical trap arrays were calculated and the results are listed at rows of the 6th and 7th in Table 2. In the 6th row, T = 40 mm, , the space between two optical traps is  μm; this is the same with the ones generated by CO₂ laser standing wave ( μm). Figure 5 gives the intensity distributions of the optical traps on the focal plane. Figures 5(a) and 5(b) show the intensity distributions of a single trap. Comparing with Figure 4, it can be found that the intensity distributions of the two kinds of optical traps are almost the same. Figure 5(d) gives the 1D intensity distribution passing through the center of the optical trap. It can be found that the intensity distributions of the optical trap array are still approximately the same as Fraunhofer diffraction if the system is focused by the modified lens, only the uniformity of the intensity distribution slightly decreases, which will affect little the formation of the optical lattices. The intensity  Wm⁻² in the table is the average value of the maximum light intensities of all optical traps in the array. In the 7th row where  mm and , the space between two optical traps is  μm; other parameters are almost the same as those in 6th row.

(a)

(b)

(c)

(d)

Figure 5: Intensity distribution of the array of optical trap with high density generated by Dammann grating (, , , ,   μm). (a) Intensity isoline of a single optical trap; (b) 2D intensity distribution of a single optical trap; (c) density distribution of the intensity for the array of optical trap (fraction); (d) 1D intensity distribution along direction that passes through the center of optical trap.

4. Potential Applications of the Optical Trap Array Generated by Dammann Grating

4.1. Optical Trapping of Cold Atoms

It is well known a two-level atom moved in an inhomogeneous light field will undergo dipole force owing to the ac Stark effect. When the light field is red detuned; that is, the laser frequency is less than the atom resonant frequency and the atom will experience gravitational effect of the light field and be attracted to the position of maximum intensity. When the absolute value of the detuning is large, the interactional dipole potential is given by [14]

On the other hand, the spontaneous emission will take place when a two-level atom moves in the light field. The scattering rate can be estimated as where is the spatial position, is the natural line width of the atoms, and is the spatial distribution of the intensity. The attractive force, , of the light field to the atoms is directly relative to the intensity gradient.

The dipole potential and spontaneous emission rate of the trapped cold atoms can be calculated after (17) and (18) with the values of and reported in references of atomic optics. Here, some atoms were chosen; and the results were calculated as listed in Table 2. It can be found that although the illuminating lights are the low-power laser in common use, the dipole potential is far higher than the temperature of the cold atoms in optical molasses (~20 μK), and the spontaneous emission rates are all low. Consequently, the optical trap can be used to trap cold atoms under this condition. Properly increasing the laser power, the array of optical trap can also be used to trap cold molecules.

4.2. Fabrication of Optical Lattices

Optical lattice, similar to the crystal structure in solid-state physics, is formed by ensemble cold atoms (or molecules) into the periodical optical potentials generated by optical trap arrays to form a periodical arrangement in spatial of cold atoms (or molecules) [18–20]. The optical lattice, on the one hand, can be used to study the dynamics process of cold atoms in the lattice, such as Bloch oscillations, quantum tunneling effect, coherent manipulation and control of atomic wave packet, and coherent transfer and control quantum states [21–23]. On the other hand, it can be used to experimentally study of quantum tunneling effect of BEC atoms and their macroscopic quantum interferences. In addition, it may have wide and important applications in the field of quantum calculation and quantum information processing [24–26].

Different kinds of optical lattices can be fabricated by using the optical trap arrays listed in Table 2. On the one hand, when the relative aperture of the optical system is relatively small and the wavelength of the illuminating light is relatively long, the lattice constant and the volume of the obtained optical trap are relatively large, as to the arrays listed in 2nd row. The lattice constant is about 0.2 mm, and the volume of a single optical trap can reach 1.65 × 10⁻⁷ cm³. To the BEC sodium atomic gas with a density of 10¹³~10¹⁴/cm3 [27], the trapped atoms of a single trap are more than 10⁶. On the other hand, when the relative aperture is relatively large and the wavelength of the illuminating light is relatively short, the lattice constant and the volume of the optical trap are very small. To the arrays listed in 6th and 7th row, the volume of a single optical trap is as small as 10⁻¹¹ cm³. For the atomic gas with thin density, each optical trap can ensemble one or two atoms. Atomic lattices with high density can be formed with small lattice constant (5.3 μm and 3.97 μm).

In recent years, researchers have proposed various schemes [28, 29] of generating optical lattices with standing wave interferences of CO₂ laser and studied them both theoretically and experimentally. Figure 6 shows the 2D array of optical trap generated by four CO₂ lasers, where two of them propagate reversely to each other. The lattice constants of the array of optical traps showed in Figures 5 and 6 are same. Comparing Figures 5(c) and 5(d) with Figures 6(a) and 6(b), one can easily find their differences. Firstly, the intensity surrounding each optical trap is zero in the optical trap array generated by Dammann grating. Though there have secondary maximum intensity between two neighbored optical trap, it is far lower than the maximum intensity. However, in the array of optical trap formed by standing wave interferences of CO₂ laser, the secondary maximum intensity is half of the maximum intensity to result in the effective depth of the optical trap which is half of the maximum dipole potential. The secondary maximum intensity will heat up cold atoms. Secondly, the intensity of the optical trap array formed by standing wave interferences of CO₂ laser changes by a rule of sine function. This will result in that, on the one hand, the intensity gradient is small and the attractive force used to trap cold atoms is correspondently small, on the other hand, the width of the optical trap is almost equal to the lattice constant resulting in thatthe condense ratio of the intensity distribution, which is the ratio between the total area of the array of optical trap and the area of the optical trap, is too small to condense the atoms together. However, in the optical trap array formed by Dammann grating, the transverse width of the optical trap is far smaller than the lattice constant. So, the condense ratio of the intensity distribution and the intensity gradient are large. Correspondingly, the attractive force of the light field to the cold atoms is large. It is easier to condense the cold atoms to a centralized space to form a lattice structure. Therefore, there are many advantages to generate array of optical trap by using Dammann grating over using standing wave interferences of CO₂ laser.

(a)

(b)

Figure 6: Array of optical trap generated by standing wave interferences of CO₂ laser ( μm). (a) Intensity distribution of 2D array of optical trap (fraction); (b) 1D intensity distribution along direction that passes through the center of optical trap.

5. Conclusion

In this paper, we calculated the diffraction intensity distributions of the Dammann gratings illuminated by Gaussian light wave. Also, we deduced the empirical equations for calculating the parameters of the optical trap array, such as the spatial period, the maximum intensity, the maximum intensity gradient, and the trap width, by numerical calculations. Different optical trap array can be obtained by adjusting the parameters of the optical system according to these equations. The optical trap array generated by Dammann grating is similar to that formed by micro-lens array however, the former has many advantages over the latter. Nevertheless, we analyzed the disadvantages of the scheme that generates the optical trap array by Dammann grating and proposed a novel experimental scheme to generate high-density optical trap array by using aplanatic lens to focus the diffraction light of Dammann grating. The result shows that the present array of optical trap with small period can be formed. Comparing with the optical lattices formed by standing wave interferences of CO₂ laser, there are many unique advantages. For instance, the optical lattices have high optical dipole potential, and the optical traps have strong attractive force to the cold atoms.

Acknowledgments

This work is supported by the National Natural Science Foundation of China under Grant no. 10674047, the Natural Science Foundation of Jiangsu Province under Grant no. BK2008183, and the State Key Laboratory for Precision Spectroscopy, East China Normal University, China.

References

H. Dammann and K. Görtler, “High-efficiency in-line multiple imaging by means of multiple phase holograms,” Optics Communications, vol. 3, no. 5, pp. 312–315, 1971.

A. Vasara, M. R. Taghizadeh, J. Turunen, et al., “Binary surface-relief gratings for array illumination in digital optics,” Applied Optics, vol. 31, no. 17, pp. 3320–3336, 1992.

S. Zhao and P. S. Chung, “Design of a circular Dammann grating,” Optics Letters, vol. 31, no. 16, pp. 2387–2389, 2006.

B. Luo, C. Wang, J. Du, C. Ma, Y. Guo, and J. Yao, “Design and analysis of phase gratings for laser beams coherent combination,” Microelectronic Engineering, vol. 83, no. 4–9, pp. 1368–1371, 2006.

S. Zhao, J. F. Wen, and P. S. Chung, “Simple focal-length measurement technique with a circular Dammann grating,” Applied Optics, vol. 46, no. 1, pp. 44–49, 2007.

P. Xi, C. H. Zhou, S. Zhao, H. S. Wang, and L. R. Liu, “Design and fabrication of 64×64× spot array Dammann grating,” Chinese Journal of Lasers, vol. 28, no. 4, pp. 369–371, 2001.

F. A. P. Tooley, S. Wakelin, and M. R. Taghizadeh, “Interconnects for a symmetric-self-electro-optic-effect-device cellular-logic image processor,” Applied Optics, vol. 33, no. 8, pp. 1398–1404, 1994.

X. He, P. Xu, J. Wang, and M. Zhan, “Rotating single atoms in a ring lattice generated by a spatial light modulator,” Optics Express, vol. 17, no. 23, pp. 21014–21021, 2009.

A. Ji, W. M. Liu, J. L. Song, and F. Zhou, “Dynamical creation of fractionalized vortices and vortex lattices,” Physical Review Letters, vol. 101, no. 1, Article ID 010402, 2008.

G.-P. Zheng, J.-Q. Liang, and W. M. Liu, “Phase diagram of two-species Bose-Einstein condensates in an optical lattice,” Physical Review A, vol. 71, no. 5, Article ID 053608, 2005.

M. Singh, M. Volk, A. Akulshin, A. Sidorov, R. McLean, and P. Hannaford, “One-dimensional lattice of permanent magnetic microtraps for ultracold atoms on an atom chip,” Journal of Physics B, vol. 41, no. 6, Article ID 065301, 2008.

G. N. Lewis, Z. Moktadir, C. Gollasch et al., “Fabrication of magnetooptical atom traps on a chip,” Journal of Microelectromechanical Systems, vol. 18, no. 2, pp. 347–353, 2009.

C. D. J. Sinclair, E. A. Curtis, I. L. Garcia et al., “Bose-Einstein condensation on a permanent-magnet atom chip,” Physical Review A, vol. 72, no. 3, Article ID 031603, 2005.

G. Birkl, F. B. J. Buchkremer, R. Dumke, and W. Ertmer, “Atom optics with microfabricated optical elements,” Optics Communications, vol. 191, no. 1-2, pp. 67–81, 2001.

R. Dumke, M. Volk, T. Müther, F. B. J. Buchkremer, G. Birkl, and W. Ertmer, “Micro-optical realization of arrays of selectively addressable dipole traps: a scalable configuration for quantum computation with atomic qubits,” Physical Review Letters, vol. 89, no. 9, Article ID 097903, 2002.

J. D. Gaskill, Linear Systems, Fourier Transforms, and Optics, John Wiley & Sons, New York, NY, USA, 1978.

C. Zhou and L. Li, “Numerical study of Dammann array illuminators,” Applied Optics, vol. 34, pp. 5961–5969, 1995.

T. Rom, T. Best, D. Van Oosten et al., “Free fermion antibunching in a degenerate atomic Fermi gas released from an optical lattice,” Nature, vol. 444, no. 7120, pp. 733–736, 2006.

R. Gerritsma, S. Whitlock, T. Fernholz et al., “Lattice of microtraps for ultracold atoms based on patterned magnetic films,” Physical Review A, vol. 76, no. 3, Article ID 033408, 2007.

J. K. Asbóth, H. Ritsch, and P. Domokos, “Optomechanical coupling in a one-dimensional optical lattice,” Physical Review A, vol. 77, no. 6, Article ID 063424, 2008.

K. D. Nelson, X. Li, and D. S. Weiss, “Imaging single atoms in a three-dimensional array,” Nature Physics, vol. 3, no. 8, pp. 556–560, 2007.

P. Kakashvili, S. G. Bhongale, H. Pu, and C. J. Bolech, “Realizing Luttinger liquids in trapped ultra-cold atomic Fermi gases using 2D optical lattices,” Physica B, vol. 404, no. 19, pp. 3320–3323, 2009.

C. Reichhardt and C. J. O. Reichhardt, “Pattern switching and polarizability for colloids in optical-trap arrays,” Physical Review E, vol. 80, Article ID 022401, 2009.

Y. Shin, C. Sanner, G. B. Jo et al., “Interference of Bose-Einstein condensates split with an atom chip,” Physical Review A, vol. 72, no. 2, Article ID 021604, pp. 1–4, 2005.

J. Mun, P. Medley, G. K. Campbell, L. G. Marcassa, D. E. Pritchard, and W. Ketterle, “Phase diagram for a bose-einstein condensate moving in an optical lattice,” Physical Review Letters, vol. 99, no. 15, Article ID 150604, 2007.

I. Bloch, “Quantum coherence and entanglement with ultracold atoms in optical lattices,” Nature, vol. 453, no. 7198, pp. 1016–1022, 2008.

K. B. Davis, M. O. Mewes, M. R. Andrews et al., “Bose-Einstein condensation in a gas of sodium atoms,” Physical Review Letters, vol. 75, no. 22, pp. 3969–3973, 1995.

S. Friebel, R. Scheunemann, J. Walz, T. W. Hänsch, and M. Weitz, “Laser cooling in a CO2-laser optical lattice,” Applied Physics B, vol. 67, no. 6, pp. 699–704, 1998.

S. Friebel, C. D'Andrea, J. Walz, M. Weitz, and T. W. Hänsch, “CO2-laser optical lattice with cold rubidium atoms,” Physical Review A, vol. 57, no. 1, pp. R20–R23, 1998.