Adiabatic cycles in quantum systems with contact interactions
Atushi Tanaka, Taksu Cheon Department of Physics, Tokyo Metropolitan University, Hachioji, Tokyo 192-0397, Japan tanaka-atushi@tmu.ac.jp, taksu.cheon@kochi-tech.ac.jp
PACS 03.65.-w,03.65.Vf,02.40.Pc DOI 10.17586/2220-8054-2017-8-1-121-126
An adiabatic cycle involving the ¿-function potential is shown to excite Bose particles confined in a one-dimensional box. We consider a cycle, during which a wall described by a <5 potential is applied and its strength and position are slowly varied. When the system is initially prepared in the equilibrium ground state, the adiabatic cycle brings all bosons into the first excited one-particle state, leaving the system in a nonequilibrium state. The absorbed energy during the cycle is proportional to the number of bosons.
Keywords: Bose particles, adiabatic cycle, < -function potential, boson.
Received: 10 July 2016 Revised: 30 July 2016
1. Introduction
Since the discovery of the geometric phase [1], it has been widely recognized that adiabatic cyclic parameter variation does not necessarily bring a quantum eigenstate of a system back to itself. In general, an "exotic" quantum holonomy [2,3], in which one eigenstate turns into another eigenstate belonging to the same Hamiltonian, can be observed after the cyclic variation of the system parameters. In this paper, we examine an adiabatic cycle that excites a system consisting of Bose particles confined in a one-dimensional box. During the cycle, we vary an additional wall adiabatically, while the interparticle interaction is kept fixed. This is in contrast to the scheme described in [4,5], where the interaction strength between Bose particles is an effective adiabatic parameter. In this study, we suppose that the wall is described by a ¿-function shaped potential [6-8]. We show that the first excited one-particle state is occupied by all the bosons to achieve the population inversion completely, if the system is prepared to be in the ground state. Namely, the energy gained by the bosons during the adiabatic cycle is proportional to the number of bosons.
The adiabatic cycle described above has a immediate relevance to the population inversion in Bose gas system [9], particularly in the super-Tonks-Girardeau gas, which has attracted considerable recent attention in both experimental and theoretical studies [4,10-12]. In the super-Tonks-Girardeau gas, which may be described by the Lieb-Liniger model [13] with strongly attractive interaction, the population inversion is created through an "adiabatic" process, where the interaction strength is suddenly flipped from infinitely repulsive to infinitely attractive [4,14]. Such a population inversion can be induced even by an adiabatic cycle, which can be obtained with an extension of the adiabatic process that connects Tonks-Girardeau and super-Tonks-Girardeau gases both to weaker repulsive and weaker attractive regime. The repetitions of this adiabatic cycle transform the ground state of non-interacting bosons into their higher excited states and achieve the population inversion [5]. This is counterintuitive, since there is no external field to drive the final state of the bosons away from the initial state.
The paper is organized as follows: in Section 2, the basic setup of our model is presented in the simplest setting of one boson; Section 3 and 4 treat several interacting and non-interacting bosons, respectively; Section 5contains the summary and the discussion. .
2. A particle in a box with a ¿-wall
In order to examine N Bose particles in a one-dimensional box with an additional ¿-wall, we review the single particle case, i.e., N =1 [8], where the system is described by the Hamiltonian:
2
H(g,X) = ^ + V(x) + g¿(x - X), (1)
2m
where m is the particle mass, V(x) is the confinement potential, and g and X are the strength and position of ¿-wall. In particular, we assume that V(x) describes an infinite square well with the length L, i.e., V(x) = 0 for 0 < x < L and V(x) = to otherwise [6,7].
We introduce an adiabatic cycle C, which consists of three adiabatic processes CI, Cn and Cm, as shown in Fig. 1. We suppose that the ¿-wall is initially absent, i.e., g = 0 in 1, and that the system is in a stationary state initially. In the first part of C, which will be called as CI, an impenetrable wall is inserted at x0 adiabatically. In
terms of the ¿-wall, the strength g is slowly increased from 0 to œ, while its position X is fixed at x0 during C/. Subsequently, in the second part Cn, the position X of the impenetrable wall is adiabatically changed from x0 to xi. In the last part Cm, the ¿-wall at X = xi is adiabatically turned off. At the end of the cycle C, the ¿-wall has no effect, again.
Fig. 1. The adiabatic cycle C of a one-dimensional box, which contains Bose particles. The strength and the position of an additional ¿-wall is adiabatically varied during C
In Fig. 2, we depict the parametric dependence of eigenenergies of the single-particle Hamiltonian H(g, X) (1) along C. Throughout this manuscript, we indicate the eigenenergy E using a normalized wavenumber k:
k -V £
where e = (hn/L)2/(2m) is the ground eigenenergy of the particle in the infinite square well.
(2)
Fig. 2. Parametric evolution of eigenenergies with N = 1 along the cycle C, which consists of CI (left part), C// (middle part) and Cm (right part). The eigenenergies are depicted by their normalized wavenumber k (2). We set x0 = 0.4703L and x1 = L — x0
The adiabatic time evolution of the single-particle system along C depends on x0 and x1. In the following, 2 1 3
we explain the case -L < x0 < -L < x1 < -L, which may be explained from Fig. 2. A more rigorous argument
5 2 5
is found in [8].
First, let us consider the case that the initial state is the ground state |1(g = 0, X = x0)} of the particle in the infinite square well, where |n(g, X)} denotes the n-th adiabatic eigenstate of H(g, X) during processes CI and Cin. We will omit to indicate (g, X) in the following. After the completion of CI, the state vector arrives at |R1}, the ground state of the right well, since we choose the right well in C// is slightly larger than the left well. During C//, there occurs a spectral degeneracy between |R1} and |L1}, the ground state of the left well. This is because the size of the left (right) well is increasing (decreasing) during Cii, and these sizes coincide at X = L/2. At the end of C//, |R1} becomes the first excited state, which adiabatically continued to |2}, which is the second excited state of the particle in an infinite square well, through C//I. Hence, the "population inversion" in the single-particle system occurs if the system is prepared to be in the ground state initially.
Second, we examine the case that the initial state is the first excited state |2}, which offers the "inverse" of the population inversion. Through the adiabatic cycle C, the system arrives at |L1} after the completion of C1, and then arrives at |1} at the end of the cycle C. Namely, either |1} and |2} return to the initial states after the completion of the adiabatic cycle C twice.
Third, let us examine the cases that the initial states are |3} and |4}, which are the first and second excited state, respectively. Now, C induces an interchange of these two states, through the intermediate states |R2} and | L2 } , which are localized the right and left well during the process II.
A similar interchange of initial eigenstates occurs as a result of the adiabatic cycle C, as long as we choose x0 and x1 appropriately. In general, the level crossing of the one-particle Hamiltonian (1) during the process Cii
plays an important role to determine which pairs of eigenstates are interchanged by C, while there is no level crossing generically during the processes CI and Cm [8].
We will now comment on the stability of the present scheme for the one-body population inversion. A crucial point is the stability of the adiabatic time evolution across the level crossing during Cn. The level crossing may be lifted due to an imperfection of the impenetrable wall, i.e., the ¿-wall with an infinite strength. If the level splitting is small enough, we may employ the diabatic process around the avoided crossing to realize the one-body population inversion. It has been shown that an open diabatic process made of a time-dependent potential well produces a second excited state from the ground state of a Bose particle [17]. This diabatic scheme is applied to create collective excitations of interacting bosons [17-19].
3. Non-interacting bosons
We first examine the case that the number of the Bose particles is N, assuming the absence of interparticle interaction. It is straightforward to extend the above result for N = 1, once we restrict the case that N bosons initially occupies the one-particle state |n). Hence, the system is in an adiabatic state of the N Bosons:
|n0N ) = |nn...n), (3)
where the one-particle adiabatic state |n) is occupied by N bosons, during CI and Cm.
If there is no interparticle interactions, the parametric evolution of averaged wavenumber & for the adiabatic N-particle state agrees with the one for the single-particle system. This suggests that the adiabatic cycle C of the N -particle system with no interaction delivers the ground state |10N) to the excited state |20N), i.e., the complete population inversion, as is seen in Figure 2. The energy that the particles acquire during the cycle C is proportional to the number of the particles.
4. Interacting bosons
We now examine the adiabatic cycle C for N interacting Bose particles. We mainly examine the case that the system is initially in the ground state. In order to confirm that the N-particle population inversion really occurs, we need to examine the effect of the interparticle interaction.
We assume that the interparticle interaction V consists of two-body contact interactions. Namely, we suppose that V takes the following form:
V (xi,x2, ...,xn )= A ¿(xi — xj), (4)
(i,f)
where A is the interaction strength, and the summation is taken over pairs. We also assume that the interparticle interaction is weak enough so that the topology of the parametric dependence of eigenenergy remain unchanged, except around the level crossing points of the noninteracting bosons. Namely, when the gaps of the eigenenergies between neighboring levels in the noninteracting system are larger than a constant value, the interparticle interaction shifts the eigenenergy at most O(A), according to the standard perturbation theory. For small enough perturbative energy correction, the corresponding adiabatic time evolution of the stationary state of the interacting bosons closely follows the evolution of noninteracting bosons.
Accordingly, under the weak interparticle interaction condition, the eigenstates of the interacting bosons can be labeled by the quantum numbers of the noninteracting bosons. For example, the ground state of the initial and final points of the adiabatic cycle C may be denoted as |10N (A)), whose overlapping integral with the unperturbed state |10N) is large. Also, |10N(A)) can be constructed by the standard perturbation theory with a small parameter A.
On the other hand, even a weak interparticle interaction can strongly influence the parametric evolution of energy levels in the vicinity of level crossings by making avoided crossings. Hence, we need to closely examine the level crossing of the non-interacting Bose particles.
In the following, we argue that the adiabatic time evolution closely follows that of the noninteracting system examined above, if the number of the particles is large enough. The key is the selection rule for the matrix element of V in the adiabatic representation in the vicinity of the level crossings of non-interacting Bosons.
4.1. "Tunneling" and direct contributions of the interaction in N = 2
We show that the effect of the interparticle interaction is significantly different, depending on whether a level crossing locates either in C//, or in C/ U Cm, as for the two body case. In the former case, the relevant matrix elements may be small since it involves only tunneling processes through the impenetrable wall. On the other hand, in processes C/ and C///, the matrix element cannot be negligible. However, it seems that there is no corresponding level crossing that affects the population inversion whose initial state is the ground state.
The parametric evolutions of eigenenergies of the noninteracting two particle system are depicted in Fig. 3, in terms of the averaged wavenumber k (see, (2)). The parametric evolution of the eigenenergy that connects 111) and 122) has a level crossing with two eigenenergies during Cn. The initial states of these energy levels are 122) and 112), which are |LiLi) and |RiLi) during Cn, respectively.
Fig. 3. Parametric evolution of normalized wavenumbers for N = 2. Full lines indicate the levels whose initial states are |nn) (n = 1,..., 4) Other levels are depicted by dashed lines. Parameters are the same as in Fig. 2
We examine the matrix elements of the interparticle interaction term V between the adiabatic basis vectors |R1R1), |L1L1) and |R1L1). Note that |R1R1) corresponds to the initial state 111) of the adiabatic cycle:
L
{R1,L1\V |R1,R1) = y/2\j {Vr (x)Vl (x)}*{Vri (x)}2dx, (5)
0
for example. Since the single-particle adiabatic eigenfunctions (x) and (x) are completely localized in the left and right wells, respectively, the overlapping integral is zero, if the ¿-wall is completely impenetrable during Cn. The level crossing accordingly remains even in the presence of the interparticle interaction. Thus the adiabatic cycle C induces complete population inversion from 111) to 122), as in the non-interacting case.
Let us examine the case that the ¿-wall during Cn allows the tunneling leakage of particles due to some imperfections. Still, we may expect that the matrix elements due to the tunneling corrections are exponentially small. Since the resultant energy gap of the avoided crossing is also exponentially small, we may expect that the diabatic process easily almost recovers the complete population inversion.
Also, during the second process Cn, the left and right part of the well may be separated. This allows us to make the tunneling correction arbitrarily small. Accordingly, the adiabatic limit that follows the extremely small avoided crossing would be difficult to realize.
Conversely, if the level crossing appears during C/ or Cin, the interparticle interaction destroys the level crossing. In Fig. 3, such an example is seen between the levels whose initial states are 133) and 124), which is delivered to |R2R2) and |L1L2), respectively, in the absence of V.
The matrix element (33 |V |24) does not vanish in general, since the relevant single-particle adiabatic eigen-functions extend the whole box. Accordingly, the level crossings are destroyed to form avoided crossing. Thus the adiabatic process CI for example, delivers 133) and 124) at the initial point of CI, to |L1L2) and |R2R2), respectively. This breaks the population inversion whose initial state is a higher excited state, e.g., the adiabatic cycle C delivers |33) to | 44) in the absence of interparticle interaction.
4.2. Selection rule for N = 3
Here, we show that the interparticle interactions do not suppress the population inversion for N > 2 due to a selection rule of V.
We explain this with the case N = 3 (Fig. 4). Let us examine the level whose initial state is |103) along C. The corresponding final state is |203) in the absence of interparticle interaction.
First, the interparticle interaction has no, or exponentially small effect on the level crossing during Cn, as shown in the case of N = 2.
Second, we examine the level crossing in CIH, where the levels whose final states are |203) and 1113) exhibit crossing. We examine the matrix element of interparticle interaction (113 |V |203), which vanishes since V is a two-body interaction, and the set of quantum numbers (1,1, 3) and (2,2, 2) has no common quantum number.
222 113
122
112
111
1223 2222 1123 1222 1113 f1122
1112
1111
\v ^
- 1223 -2222 -1123 -1222 1113 -1122
1112
1111
II
III
II
III
Fig. 4. Parametric evolution of normalized wavenumbers for N = 3. Other parameters are the same as in Fig. 3
Fig. 5. Parametric evolution of normalized wavenumbers for N = 4. Other parameters are the same as in Fig. 3
Still, there may be a tiny avoided crossing whose magnitude can be explained by the standard second-order perturbation theory. We may expect that the diabatic process induces complete population inversion whose final state is |203}. Also, even if the interaction strength A is moderately large, where the topology of the level diagram remains unchanged except that the avoided crossing becomes noticeable, the final state should be 1113}, whose energy is far larger than the ground state. In this sense, an incomplete population inversion should be realized.
4.3. The population inversion for N > 2
We shall prove that the adiabatic cycle C delivers |10N} to |20N} for N > 2, even in the presence of two-body interparticle interaction. Here, we explain the selection rule for arbitrary N (> 2), and examine each part of the cycle C. For example, N = 4 case is shown in Fig. 5.
We explain the selection rule of the two-body interparticle interaction for N > 2. Namely, we examine the matrix element (n|n2... n'N |V|nm ... nN}. The matrix element vanishes when the two sets of quantum numbers (n'1; n2,..., n'N) and (ni,n2,...nN) have at least three different elements, i.e., the number of the common quantum numbers is equal to N — 3 or less. In other words, non-vanishing matrix element has the following (n1n2n3 ... nN|V|n1n2n3 ... nN} where (n3,..., nN) are the common quantum numbers.
We examine the first part C/ of C. We assume that the system is initially in the ground state |10N(A)}. According to the selection rule, it is sufficient to examine |10Nwhere and are single particle adiabatic states, e.g. |2}. Now, we examine whether the eigenenergies of these states are degenerate. This is equivalent to comparing the eigenenergies corresponding to 111} and of the two particle system. As is seen in Fig. 3, there is no level crossing in C/. In this sense, there is no effective level crossing with the level |10N(A)}, during C/.
As for C///, we conclude from a similar argument above, that the energy level corresponding to |20N(A)} has no effective level crossing.
Next, we examine C//, where the system is in |RfN}. According to the selection rule, it suffices to examine |RfNwith single particle adiabatic states and To clarify the level crossing, we compare |R1R1} with There are three cases. First, the levels corresponding to |R1R1} and |Rn, Rn>} ((n, n') = (1,1)) do not
occur. Second, the levels corresponding to |R1R1} and |Rn, Ln} exhibit a degeneracy only when n = 1 and n' = 1, where the corresponding matrix element involves a single-particle tunneling. Third, the levels corresponding to |R1R1} and |Ln, Ln>} exhibit a degeneracy only when n =1 and n' = 1, where the corresponding matrix element involves two-particle tunneling. Since the matrix elements involving tunneling contribution are exponentially small, the resultant gap should also be small. Hence, the diabatic process should occur even when the speed of the impenetrable wall is moderately slow.
5. Discussion and summary
Here, we argue that the experimental realization of the population inversion suggested in this paper is feasible with the current state of the art. For example, we may utilize the scheme [20] to realize ¿-wall with an approximate Gaussian wall.
Another possibility is to use a heavy particle as a wall, whose position may be manipulated by, say, an optical tweezer. The effective interaction between the wall particle and other particle may be tuned by external fields.
We note that the present scheme may offer a way to realize other exotic nonequilibrium states. Let us suppose, for example, the state of bosons is in |20N(A)), which can be generated from the adiabatic cycle C. After the interparticle interaction A is adiabatically increased to то, the system arrives the higher excited state of the Tonks-Girardeau system, which may be described by the Lieb-Linigher model with the infinite interparticle interaction strength [21,22]. Similarly, after A is adiabatically decreased to -то, the system now arrives at the higher excited state of the super-Tonks-Girardeau system [4]. This state is a much more highly-excited state compared to the super-Tonks-Girardeau state, because the initial state |20N(A)) is a higher excited state of noninteracting bosons.
We note that the present scheme offers a way to realize the super-Tonks-Girardeau gas in a one-dimensional box. For example, once we obtain |20N) from |10N) through the adiabatic cycle C, the system arrives at the super-Tonks-Girardeau gas as the strong attractive interparticle interaction is imposed. This scheme do not require neither the preparation of the Tonks-Girardeau gas [21,22] nor the sudden flip of the interparticle potential [4].
In summary, we have shown that the adiabatic cycle C induces the nearly complete population inversion of the multi-boson system, when the interparticle interaction is not too strong. As pointed out in [8], for a single particle case, the present scheme may be extended to the case of an arbitrary shape of the confinement potential V(x).
Acknowledgements
This research was supported by the Japan Ministry of Education, Culture, Sports, Science and Technology under the Grant number 15K05216.
References
[1] Berry M.V. Quantal Phase Factors Accompanying Adiabatic Changes. Proc. Roy. Soc. A, 1984, 392, P. 45-57.
[2] Tanaka A., Miyamoto M. Quasienergy Anholonomy and its Application to Adiabatic Quantum State Manipulation. Phys. Rev Lett., 2007, 98, P. 160407.
[3] Tanaka A., Cheon T. Bloch Vector, Disclination and Exotic Quantum Holonomy. Phys. Lett. A, 2015, 379, P. 1693-1698.
[4] Haller E., Gustavsson M., et al. Realization of an Excited, Strongly Correlated Quantum Gas Phase. Science, 2009, 325, P. 1224-1227.
[5] Yonezawa N., Tanaka A., Cheon T. Quantum Holonomy in the Lieb-Liniger Model. Phys. Rev. A, 2013, 87, P. 062113.
[6] Fliigge S. Practical Quantum Mechanics, vol. 1. Berlin, Springer-Verlag, 1971.
[7] Ushveridze A.G. Analytic Properties of Energy Levels in Models with Delta-function Potentials. J. Phys. A, 1988, 21, P. 955-970.
[8] Kasumie S., Miyamoto M., Tanaka A. Adiabatic Excitation of a Confined Particle in One Dimension with a Variable Infinitely Sharp Wall. Phys. Rev. A, 2016, 93, P. 042105.
[9] Grynberg G., Aspect A., Fabre C. Introduction to Quantum Optics. Cambridge, Cambridge University Press, 2010.
[10] Astrakharchik G.E., Boronat J., Casulleras J., Giorgini S. Beyond the Tonks-Girardeau Gas: Strongly Correlated Regime in Quasi-One-Dimensional Bose Gases. Phys. Rev. Lett., 2005, 95, P. 190407.
[11] Olshanii M. Atomic Scattering in the Presence of an External Confinement and a Gas of Impenetrable Bosons. Phys. Rev. Lett., 1998, 81, P. 938-941.
[12] Guan X. Critical Phenomena in One Dimension from a Bethe Ansatz Perspective. Int. J. Mod. Phys. B, 2014, 28, P. 1430015.
[13] Lieb E.H., Liniger W. Exact Analysis of an Interacting Bose Gas. I. The General Solution and the Ground State. Phys. Rev, 1963, 130, P. 1605-1616.
[14] Girardeau M. Relationship between Systems of Impenetrable Bosons and Fermions in One Dimension. J. Math. Phys., 1960, 1, P. 516-523.
[15] Vitanov N.V., Halfmann T., Shore B.W., Bergmann K. Laser-induced Population Transfer by Adiabatic Passage Techniques. Annu. Rev. Phys. Chem., 2001, 52, P. 763-809.
[16] Guerin S., Yatsenko L.P., Jauslin H.R. Dynamical Resonances and the Topology of the Multiphoton Adiabatic Passage. Phys. Rev A, 2001, 63, P. 031403(R).
[17] Karkuszewski Z.P., Sacha K., Zakrzewski J. Method for Collective Excitation of a Bose-Einstein Condensate. Phys. Rev A, 2001, 63, P. 061601.
[18] Damski B., Karkuszewski Z.P., Sacha K., Zakrzewski J. Simple Method for Excitation of a Bose-Einstein Condensate. Phys. Rev. A, 2002, 65, P. 013604.
[19] Karkuszewski Z.P., Sacha K., Smerzi A. Mean Field Loops versus Quantum Anti-crossing Nets in Trapped Bose-Einstein Condensates. Eur. Phys. J. D, 2002, 21, P. 251-254.
[20] Meyrath T.P., Schreck F., et al. Bose-Einstein Condensate in a Box. Phys. Rev. A, 2005, 71, P. 041604(R).
[21] Paredes B., Widera A., et al. Tonks-Girardeau Gas of Ultracold Atoms in an Optical Lattice. Nature, 2004, 499, P. 277-281.
[22] Kinoshita T., Wenger T., Weiss D.S. Observation of a One-Dimensional Tonks-Girardeau Gas. Science, 2004, 305, P. 1125-1128.