Vortex conveyor belt for matter-wave coherent splitting and interferometry.

We numerically study a matter wave interferometer realized by splitting a trapped Bose-Einstein condensate with phase imprinting. We show that a simple step-like imprinting pattern rapidly decays into a string of vortices that can generate opposite velocities on the two halves of the condensate. We first study in detail the splitting and launching effect of these vortex structures, whose functioning resembles the one of a conveyor belt, and we show that the initial exit velocity along the vortex conveyor belt can be controlled continuously by adjusting the vortex distance. We finally characterize the complete interferometric sequence, demonstrating how the phase of the resulting interference fringe can be used to measure an external acceleration. The proposed scheme has the potential to be developed into compact and high precision accelerometers.

Introduction

Atom interferometry is one of the pillars of the emerging quantum technologies with applications ranging from the measurement of fundamental physical constants[1,2] to tests of the equivalence principle[3,4] and the detection of gravitational waves[5-7]. Additionally, atom interferometers are promising candidates for navigation and geophysics surveys as they have shown excellent performance in the field of acceleration and rotation measurements[8,9]. Atom interferometers can be implemented with collimated atomic beams[10,11] or cold atoms[12], for which coherent splitting and recombination of the matter wave can be achieved by using microfabricated gratings[11], Bragg scattering[10] or Kapitza-Dirac scattering[10,13] from standing light waves, or stimulated Raman transitions[12]. Nowadays free-falling atom interferometers can reach record accuracy up to of 3.9 μGal[14], surpassing other state-of-the-art techniques. However, as they are reaching their ultimate performances, it’s quite challenging to make any further improvement. For example, to improve their sensitivity, it is necessary to increase the interrogation time, but this comes at the price of increasing sizes of the atomic clouds or of the atomic beams, making it harder to manipulate the atoms and extract the useful signal. More recently, advanced cooling techniques have allowed to achieve interrogation time of 2.3 s[15], however as the atoms are in free fall, the vacuum enclosure has to be on the order of 10 m length[15]. In many applications, like inertial navigation and geodesy, such large scale apparatuses are not practically feasible. In addition, vibrations[16], wave-front distortion[17], and detection noise[16] are other known limiting factors in the performances of these atom interferometers.

In the constant quest for more sensitive and portable devices, a new generation of interferometers based on trapped and guided atomic samples is emerging. As the atoms are held against gravity by the confining potentials, a trapped and guided interferometer promises to provide very long interrogation time within a small apparatus size. This is helpful to suppress phase shifts due to stray fields and mechanical vibrations[18] and increase the sensitivity. Besides, it allows to have precise control over the atomic wavefunction, which makes trapped interferometers ideal sensors for inertial navigation[19-21], where a large dynamical range is desirable, and for testing (sub-)gravitational and surface-generated forces[22,23], due to the high spatial resolution. The many advantages offered by trapped interferometers bear the promise to deliver a new generation of devices able to exceed the capabilities of the free-space counterparts[24]. Additionally, they allow the use of Bose-Einstein condensates (BECs), that are particularly appealing systems for the realization of precise interferometers due to their properties of macroscopic phase coherence[25]. For example, besides Bragg scattering[6,26], and Kapitza-Dirac scattering[27], the BECs can also be coherently split by deforming a single-well potential to a double-well potential, which can be realized by changing the optical trap[28], or by applying radio-frequency[19] or microwave fields[29]. Additionally, the possibility of controlling the interactions, that are often viewed as a major drawback, allows to engineer entangled states[30] or to use bright solitons[31,32], further enhancing the sensitivity of the device.

Phase imprinting methods have been proposed to manipulate the phase of the BEC with optical potentials[33]. The manipulation can be almost arbitrary in two dimensions using the widely employed spatial light modulators[34,35]. This powerful technique has been mainly used to study vortices[33] and solitons[34-36] in BECs but has also potential for applications in atom interferometry[37]. For example, it has been shown that interference can be produced in momentum space by imprinting part of a trapped quasi-one-dimensional BEC with a detuned laser[37], similarly to the generation of solitons[34,35]. In this work, we study an alternative and simple way to build a trapped matter-wave interferometer using phase imprinting. In our scheme, two opposite velocities are imprinted on two halves of a trapped BEC. We find that this imprinting generates a string of vortices along the central line of the BEC. Our numerical simulations show that the vortices effectively act as a conveyor belt which moves the atoms on its two sides in opposite directions. As a result, the BEC is coherently split into two clouds that, after half oscillation in the trapping potential, recombine in the centre, forming an interference pattern. We demonstrate that the phase of the interference fringe can be used to infer the magnitude of an external acceleration over a large dynamical range.

Results

Superfluidity is one of the most spectacular consequences of Bose-Einstein condensation[38]. At the mean field level, the trapped BECs can be well described by the Gross-Pitaevskii equation, which is a nonlinear Schrödinger equation:where Φ(r, t) is the wave function of the condensate (order parameter), m is the mass of the atom, V(r) is the trapping potential, and ξ = 4πℏ2a/m in which a is the s-wave scattering length. The complex order parameter Φ(r, t) can be written aswhere n(r, t) is the condensate density distribution, and S(r, t) is the phase distribution. The gradient of the phase S(r, t) fixes the velocity of the superfluid through[38]

From Eq. (3) it is easy to see that manipulating the phase of a BEC allows to spatially and temporally control its velocity field. Such manipulation can be achieved by exploiting the dipole potential , where η is the complex polarizability of the atom, c the speed of light, ε0 the vacuum permittivity and I(r) is the intensity distribution of the “imprinting” laser beam, which can be arbitrarily manipulated using e.g. a spatial light modulator[36]. This locally sets the phase of the BEC as S(r) = U(r)τ/ℏ, where τ is the illumination time. By suitably adjusting the length of the pulse τ and the intensity of the “imprinting” laser beam I(r), it is possible to imprint on the condensate arbitrarily shaped phase pattern.

Let us now consider a BEC trapped in a harmonic potential with frequencies (ω, ω, ω) where z is the direction of propagation of the imprinting beam. To optimize the imprinting, it is convenient to set , leading to a pancake-shaped cloud. Throughout this work, unless otherwise stated, we will express distances in units of the harmonic oscillator length , time in units of and energies in units of . To illustrate a specific realistic case we consider a BEC of 3 × 104 87Rb atoms with ω = ω = 0.06ω. We numerically simulate the dynamics of our trapped interferometer starting by imprinting, at t = 0, two opposite velocities on the two halves of the condensate using the following pattern:where the coefficient l fixes the phase gradient and therefore the induced velocity field. In Fig. 1(a) we show the phase pattern resulting from imprinting φ1 on the BEC while the corresponding velocity field is reported in Fig. 1(c). The periodicity of the phase pattern along the y-axis is l. We simulate the evolution following the phase imprinting by solving Eq. (1) using a standard split-step Fourier algorithm. The resulting dynamics for the column density distribution (which is often the observable in BEC experiments) is shown in Fig. 2. The BEC is first split into two clouds that separate and recombine due to the effect of the harmonic trap. When recombining, after approximately half period, a clear interference pattern is formed in the density profile of the BEC.

Figure 1

(a) Phase pattern resulting from the imprinting of Eq. (4) with l = 1. (b) Phase pattern generated by Eq. (5) with d = 0.25. (c) and (d) are the corresponding velocity fields.

Figure 2

Column density profiles of the BEC at different evolution times t after the phase imprinting. (a)–(f) shows the dynamics after the imprinting of Eq. (4), while (g)–(l) of Eq. (5). Note that for t = 0.24 and 7.7 the axis range is reduced to better show the initial formation of the vortex conveyor belt and the final interference patterns.

We first analyse in detail the dynamics immediately following the initial imprinting. We observe that the initially uniform density distribution rapidly decays into a regular line of vortices pinned at x = 0 and separated by l/2 (first panel of Fig. 2(a)). This behaviour can be understood by looking at the imprinting pattern in Fig. 1(a), that presents a periodicity of l along y and a jump across x = 0. The vortex line then behaves as a conveyor belt for the matter wave (each vortex corresponds to a pulley in a belt conveyor system) launching the two halves of the BEC in opposite directions. To better characterize the splitting and launching action of the conveyor belt and isolate its contribution to the overall dynamics, we have performed a series of numerical simulations where we directly imprint such vortex string on the BEC. A vortex is a topological feature of a superfluid characterized by the fact that, in a closed path around the vortex, the phase of the condensate wave function undergoes a 2π winding[39], which is due to the single-valuedness of the condensate wave function. Therefore, to directly generate the vortex conveyor belt we imprint on the BEC the following pattern:where N is a positive integer, and d is a positive constant. This phase pattern therefore creates singly quantized vortices along the x = 0 line with an initial distance of 2d between them. Figure 1(b) shows an example of the imprinted phase pattern and Fig. 1(d) shows the corresponding velocity field. Comparing Fig. 1(c) and Fig. 1(d), we note that when d = l/4 the imprinted velocity fields are essentially equivalent for much larger than the healing length. On the contrary, close to x = 0 the phase pattern resulting from φ2 is clearly distorted and intrinsically presents vorticity, that instead is absent in the phase pattern generated by φ1. By comparing the evolution following the two imprinting patterns with the condition d = l/4 shown in Fig. 2, it is possible to appreciate the similarity between the two cases, confirming that the vortex conveyor belt appearing after the imprinting of φ1 is responsible for the splitting of the condensate and launching of the two halves in opposite directions. The main difference is represented by the fact that excess vorticity generated by φ2 induces an additional small rotation to the two separating clouds that results in a more distorted phase pattern.

Clearly, due to the similarity of the phase patterns φ1 and φ2 when d = l/4, the initial exit velocity along the y-axis is inversely proportional to the distance between the vortices. To study this influence, we performed a series of simulations with different vortex distances. The initial exit velocity along the y-axis v was extracted from the column density profiles with two steps. First, the centre of mass positions of the two clouds were determined by fitting the column density profiles with two 2D Gaussian profiles. Due to the existence of the harmonic trap, the centre of mass motion of the clouds is almost harmonic. Therefore, in the second step, we fit the centre of mass motion with a cosine function to obtain the initial exit velocity. Figure 3 illustrates the result we obtained, confirming that the initial exit velocity along the y-axis is inversely proportional to the vortex distance 2d. For comparison, in Fig. 3 we also plot the curves corresponding to v = dφ1/dy with l = 4d, which correspond to the exit velocities expected imprinting φ1. The discrepancy between the velocities obtained with the two methods is again due to the small additional rotation generated by φ2. Figure 3 shows that the energy injected into the BEC decreases with an increasing initial vortex distance. In fact, the energy injected can also be compared by looking at the number and distribution of the vortices —the energy of a single vortex is higher if it’s closer to the centre of the BEC. If we keep increasing the vortex distance, the scheme eventually breaks down as the conveyor belt does not have enough ‘pulleys’ to efficiently and uniformly split and then launch the two halves of the BEC. We found that for , it is not possible to generate an interference pattern at the end of the interferometric sequence. In these cases, the vortices created sit stably within the condensate. On the opposite extreme, there is no lower boundary for the vortex distance, meaning that multiply charged vortices can be generated launching the two clouds with increasing exit velocities.

Figure 3

Initial exit velocity along the y-axis v vs. half initial vortex distance d after the imprinting of the pattern in Eq. (5). The blue circles and the red squares represent the values obtained by fitting the simulated density profiles with the procedure explained in the text. The solid lines are calculated with v = dφ1/dy with l = 4d.

Let us now characterize the performance of the full interferometric scheme and in particular its ability to detect inertial forces. To this end, we add a constant acceleration β towards the negative direction of the y-axis. To demonstrate the potential of our scheme for applications requiring large dynamical ranges, we vary β in a range spanning 5 orders of magnitude —from 10−2 to 10−7. Even in the presence of the external force, the total potential is still harmonic, but its centre is shifted towards lower y values. As a consequence, the initial position of the BEC changes with the magnitude of the acceleration, while the phase imprinting pattern stays constant.

To measure the acceleration, we fit the density profile of the recombined clouds at the end of the whole interferometer with the product of a Gaussian and a sine function. From this we extract the phase of the interference fringe α, which is sensitive to the external force. The values obtained for different imprinted velocities are displayed in Fig. 4, together with the corresponding linear fits. It can be seen that the relationship between the phase and the acceleration is perfectly linear, demonstrating the validity of our interferometer. As expected, the sensitivity of the accelerometer, given by the slope of the linear fits in Fig. 4, can be increased by decreasing l. For example, for the simulations shown in Fig. 4, the sensitivity increases from 102 to 166 when l is decreased from 1.38 to 0.92. Our simulations again indicate that when we imprint φ2, i.e. the sole vortex conveyor belt, at least for high accelerations the performance is comparable with the one obtained by imprinting φ1 with l = 4d. However, as shown in Fig. 4 the sensitivity is reduced by roughly 10% (for the curve in figure we obtain 89). To give a specific example, if we choose Hz, which are frequencies that can be easily obtained in standard experiments, the curves in Fig. 4 range between 1 and 105 μgal.

Figure 4

Phase of the interference fringe α vs. acceleration β for different values of the phase imprinting coefficient l or d. The filled symbols are the results of the numerical simulations after imprinting the pattern in Eq. (4). The solid lines are linear fits to the data. The open points are data extracted from the simulations resulting from the imprint of Eq. (5). The dashed line is the corresponding linear fit.

Discussion

We have shown that a trapped BEC can be coherently split with a simple phase pattern that generates a vortex conveyor belt. The conveyor belt further launches the two halves of the condensate in opposite directions initiating an interferometric sequence driven by the external trapping potential. We have characterized the initial exit velocity of the two wave packets as a function of the distance between the vortices in the conveyor belt, finding the limit for which the scheme breaks down. We have demonstrated that the interferometer generated is sensitive to external accelerations and that the sensitivity can be improved by increasing the corresponding phase gradient of the imprinted phase pattern. Due to the wide dynamical range and the possibility of being implemented in a compact setup, our interferometer has the potential to be developed into compact and high precision devices. Towards its practical implementation, there are many factors that can affect the acceleration sensitivity. For example, the column density profiles of the BEC are often detected optically. In this case, the optical resolution of the imaging system will limit the resolution of the column density profile and therefore the phase extracted from the interference pattern. Additionally, the optical resolution of the phase imprinting system would set a lower boundary for the l we can realize. Besides, any factors that can influence the evolution of the BEC, such as the fluctuations of the external potential, may also affect the acceleration sensitivity.

Methods

Numerical simulations

For every value of the external acceleration, we first find the ground state of the condensate by finding the solutions of the corresponding time-independent Gross-Pitaevskii equation. This is achieved by iterating in the imaginary time along the steepest descent of the energy. Once the ground state is obtained, we use a standard split-step Fourier algorithm to simulate the dynamics governed by Eq. (1). In our simulations, we employ a grid of (512 × 512 × 8) points, except when we consider the imprinting of Eq. (4) with or the imprinting of Eq. (5) with d < 0.28, in which cases a grid of (512 × 1024 × 8) points is employed to ensure that the atoms are always in the calculation area while the spatial resolution stays the same.

Interference fringes

We fit the column density profiles at the end of the interferometric sequence with the following equation:where A describes the amplitute and B is the offset. It mainly includes a 2D elliptical Gaussian function times a function that describes an ideal interference pattern. In the Gaussian function, x0 and y0 show the center of mass position of the BEC, while a1, b1 and c1 fix the shape and size of the Gaussian distribution. In the function that describes the interference, a2 is the fringe spatial frequency, θ1 shows the angle of the interference fringe, and α is the phase of the interference fringe. An example of our fitting procedure is shown in the Supplementary Fig. S1.