Dynamics of Bose-Einstein Condensates Subject to the Pöschl-Teller Potential through Numerical and Variational Solutions of the Gross-Pitaevskii Equation.

We present for the first time an approach about Bose-Einstein condensates made up of atoms with attractive interatomic interactions confined to the Pöschl-Teller hyperbolic potential. In this paper, we consider a Bose-Einstein condensate confined in a cigar-shaped, and it was modeled by the mean field equation known as the Gross-Pitaevskii equation. An analytical (variational method) and numerical (two-step Crank-Nicolson) approach is proposed to study the proposed model of interatomic interaction. The solutions of the one-dimensional Gross-Pitaevskii equation obtained in this paper confirmed, from a theoretical point of view, the possibility of the Pöschl-Teller potential to confine Bose-Einstein condensates. The chemical potential as a function of the depth of the Pöschl-Teller potential showed a behavior very similar to the cases of Bose-Einstein condensates and superfluid Fermi gases in optical lattices and optical superlattices. The results presented in this paper can open the way for several applications in atomic and molecular physics, solid state physics, condensed matter physics, and material sciences.

1. Introduction

The advent of Bose–Einstein Condensation opened the way for a better understanding of ultracold atoms trapped [1,2]. The dynamic and static properties of Bose–Einstein condensates (BEC) in the ultracold temperature regime are well described by the mean field equation known as the Gross–Pitaevskii equation (GPE) [3]. The manipulation of GPE has generated several applications in many areas of physics such as condensed-matter physics [4], nonlinear optics [5], optical solitons in nonlinear media (optical fibers) [6], metamaterials [7,8], etc.

The first BECs were produced using magnetic traps which are described by parabolic potentials [9,10]. However, there are other types of traps that are used in the confinement of ultracold atoms. Among the traps used in BECs, we can mention the optical traps known as optical lattices: periodic potentials created through the interference between lasers [4]. Since the optical lattices configurations are similar to the geometric arrangements of atoms in crystalline solids, they can be used to study atomic behavior in a highly controlled environment [11,12,13,14]. Thus, the ultracold atoms present in the optical lattice become a perfect environment for the investigation of nonlinear phenomena and solid-state physics [15,16]. There is also the possibility that magnetic and/or optical traps can be combined experimentally together or with other potentials. An important example of such a combination is the double-well potential [17,18].

Although there are many papers on confinement of ultracold atoms through the trap potentials mentioned above, it is interesting, from a theoretical and experimental point of view, to explore other atomic confinement potentials. To date, there is little research that considers the Pöschl–Teller potential (PT) as a potential for trapping cold atoms.

The PT potential due to its interesting physical properties has been studied in several situations. Jia et al. modeled a diatomic molecule using the PT potential from the experimental Rydberg–Klein–Rees (RKR) potentials for six diatomic molecules and obtained accurate results in fitting experimental RKR potential curves over a large range of internuclear distances for six molecules examined [19]. In addition, Yahia et al. obtained thermodynamic properties such as vibrational mean energy, vibrational specific heat, vibrational mean free energy and vibrational entropy for the PT-type potential for the Schrödinger and Klein–Gordon equations using the parametric Nikiforov–Uvarov method [20]. According to Park [21], two-dimensional transport of quasiparticles in bilayer graphene were also studied through a PT potential. The PT potential is used to account for the physical properties of many systems which includes the excitons, quantum wires, and quantum dots [22]. A very interesting property of this potential is the absence of reflection [23,24]—from the perspective of materials science, it can have formidable applications such as waveguides, solitons propagation in nonlinear media, complex systems, etc.

In recent years, some studies have shown significant results about BECs subject to the PT potential. In Ref. [25], a BEC of noninteracting atoms in this trapping potential, vortices in low-density neutron matter, and cold Fermi gases in Ref. [26], theoretical studies on cigar-shaped BEC in Ref. [27], and also, within cosmology, black holes as a gravitons condensate in Ref. [28]. However, from a general point of view, we can conclude that there are still few studies that explore the dynamics of BECs trapped in this potential. Thus, from the PT potential, new phenomena can be studied and understood in the BECs.

Motivated by the properties mentioned above and other properties that will be mentioned in Section 2.3, we propose in this paper to obtain variational and numerical solutions of EGP for a one-dimensionally confined BEC in a PT potential with attractive interactions. For the variational approximation (VA), we will use a Gaussian ansatz to minimize the effective Lagrangian of the proposed model [29]. For the numerical approximation, we will consider algorithms involving real- and imaginary-time propagation based on a split-step Crank–Nicolson (SSCN) method [30].

The paper is structured as follows: In Section 2, we formulate the dynamical model for our system and present how to obtain a one-dimensional GPE through a transversal parabolic confinement (harmonic potential) and a longitudinal hyperbolic confinement (Pöschl–Teller potential). In Section 3, we present the variational and numerical results obtained for the 1D-GPE. In Section 4, we discuss the results obtained and compare them with analytical results. In addition, finally, in Section 5, we present the conclusions and final considerations about this paper.

2. The Dynamical Model

2.1. The Gross–Pitaevskii Equation

The static and dynamical properties of a pure BEC made of N dilute and ultracold atoms are described by the three-dimensional Gross–Pitaevskii equation (3D-GPE) as follows: where is the imaginary unit, is the reduced Planck constant, m is the mass of each atom. The operator is the Laplacian operator in three-dimensional cartesian coordinates given by , the term is the external trapping potential and the nonlinearity coefficient, is the strength of interatomic interaction, with the s-wave scattering length of atoms in the condensate; it is positive for repulsive interactions and negative for attractive interactions. In these equations, the quantity is the wave function of the BEC (also called the order parameter or macroscopic wave function) and, since is defined as the atomic density of the N condensed atoms [31], Equation (1) will be supplemented by the normalization condition

2.2. The External Potential

The external potential is responsible for BEC confinement. As mentioned in Section 1, there are several ways to confine ultracold atoms. In the case of the “standard” magnetic trap (MT), the external potential assumes the typical harmonic form where , , and are the trapping frequencies in the three different directions. The geometry of the trap and, hence, the shape of the condensate, may range from isotropic forms, to strongly anisotropic ones [32].

In this paper, we consider a BEC confined in a strongly elongated trap known as a “cigar-shaped” trap. In this configuration, which is highly anisotropic, the longitudinal and transverse confining frequencies (denoted by and , respectively) are such that , where and . This trap setting is interesting because the transverse confinement of the BEC is very narrow and causes the significant dynamics of cigar-shaped BEC to only occur in the longitudinal direction and thus the system effectively becomes one-dimensional. [31,32]. Thus, the external potential, , can be interpreted as a combination of a longitudinal, , and a transverse potential, , that is,

In this paper, our objective is to reduce the 3D-GPE to one-dimensional Gross–Pitaevskii equation (1D-GPE) form, assuming, as mentioned above, that the BEC is held in a relatively tight cigar-shaped that correspond to the following transverse potential:

Based on these assumptions, we may decompose the wave function of the cigar-shaped BEC, , in a longitudinal (along x) and a transverse (on the plane) component, and the system can be governed by 1D-GPE by assuming separable solutions on the form where the wave function is well described by a solution of the transverse quantum harmonic-oscillator problem, where is the transverse Laplacian operator, and is the energy associated with transverse confinement of the trap. Since the considered system is effectively 1D, it is natural to assume that the transverse component of the condensate wavefunction, , remains in the ground state [32]; Therefore, takes the form with being the ground state energy corresponding to transverse confinement.

Substituting Equation (8) into Equation (1) and integrating it over the transverse coordinate, we derive an effective 1D-GPE with being the strength of one-dimensional interatomic interaction, is the transversal harmonic oscillator ground state width and the coefficient characterizes the type of the s-wave interactions. In this paper, we propose a BEC with attractive interatomic interactions, i.e., . It is found in the literature that the term is sometimes omitted because it does not affect the dynamics of the condensate [33].

2.3. The Pöschl–Teller Potential

The Pöschl–Teller potential [34] (here, the longitudinal potential ) is a hyperbolic potential described mathematically by where and describe the amplitude/depth and width of this potential, respectively. The number of bound states of this potential is determined by a parameter l characterizing the depth of the potential energy [35]. This occurs when we rewrite the amplitude as and, as a consequence, by means of algebraic manipulations, we can reduce the one-dimensional Schrödinger equation subject to the PT potential to the Legendre differential equation [36,37]. Thus, the 1D-GPE for our model is given by

Finally, we can rewrite Equation (10) through the dimensionless variables:

Therefore, Equation (10) takes the following rescaled form: with the rescaled 1D wave function subject to normalization

3. Results

The main equation in our model (Equation (12)) that describes a BEC trapped by the potential PT in a regime with attractive interactions can be rewritten as where and .

Provided the external potential is time independent, , stationary solutions of the Equation (14) can be determined by expressing the condensate wave function as where is a real function and is the chemical potential which can be obtained from . Substitution of the Equation (15) into Equation (14) yields the following steady-state equation for :

3.1. Variational Results

From the point of view of the variational method, Equation (14), together with normalization condition (Equation (13)), can be derived from the following Lagrangian density

The variational solutions for Equation (14) were obtained by assuming an ansatz Gaussian form where the variational parameters of soliton are the norm M and width w. It is important to note that the norma of the wave function, , represents a conserved quantity.

The effective Lagrangian, L, is given by the integral where the chemical potential, , was introduced to ensure that parameter M maintains the correct normalization of the wave function [15,16]. Substituting Equation (29) in Equation (18), we get and the substitution of ansatz given by Equation (17) in above effective Lagrangian yields:

The evolution equations for variational parameters can be obtained from the Euler–Lagrange equation: where q are the generalized coordinates .

The first variational equation, , recovers the normalization adopted in Equation (13), that is, ; which is substituted in the final forms of other equations below, except for equation , where is substituted after the differentiation. The other equations, and , yield a set of coupled nonlinear integral equations which relates w and :

3.2. Numerical Results

The numerical solutions of the Equation (14) were obtained through the Split-Step Crank–Nicolson (SSCN) method, using imaginary time propagation in order to obtain the BEC stationary ground state wave function. This method is well explained in Refs. [30,38]. The dimensioneless 1D-GPE, given by Equation (14), can be rewritten as

The Hamiltonian, , can be expressed by , with the kinetic energy operator and the potential energy operator with and . The term represents the nonlinear EGP term which will be treated as a potential. The solution of Equation (24) is where . The operator can be approximated via the Baker–Campbell–Hausdorff formula [39]. Thus, Equation (25) becomes

To apply the SSCN method, first, Equation (24) is discretized in the domain using a spatial step and a temporal step . Then, we assume the initial condition, , and the boundary conditions and . For the stationary ground state, the wave function is essentially real and the method of propagating imaginary time (relaxation method) that deals with real variables seems to be convenient. The basic idea of the propagation of imaginary time is to replace the EGP temporal variable, t, with the imaginary temporal variable .

The temporal evolution of Equation (24) can be written in terms of the eigenfunctions with their respective eigenvalues

The propagation of Equation (24) in imaginary time alters the temporal evolution given by Equation (27) to

Thus, all eigenfunctions will decay exponentially over time. However, all states excited will decay exponentially faster compared to the ground state. Consequently, only the ground state survives. That is, this procedure will converge to the lowest energy ground state solution [40]:

The real-time propagation of the SSCN method preserves the normalization of the wave function, while the propagation of the imaginary time of the SSCN method does not preserve normalization. This problem can be solved by restoring the normalization of the wave function after each step is performed. Thus, the method of propagating imaginary time for stationary problems in the ground state produces very accurate results at low computational cost.

In this paper, we use normalized Gaussian as the initial condition and the Dirichlet boundary conditions to represent .

4. Discussion

In particular, in the case , for which results are presented below in Figure 1, Equation (16) with predicts the analytical solution whose corresponding chemical potential is . Both variational and numerical results proved to be accurate in obtaining the solution of Equations (14) and (16) when compared with the analytical solution. These results are important from the point of view of the precision analysis of the variational and numerical method because the analytical solution becomes a “calibration” parameter for both methods. Variational and numerical results were also obtained for other values of the PT potential depth. Figure 2 shows the wave functions obtained in the simulations.

Figure 1

Comparison of variational and numerical results with the analytical solution of Equation (16) for and . The analytical solution, represented by the continuous curve (green), is given by the expression . The variational solution, represented by the dashed curve (blue), was obtained using the Gaussian ansatz whose width, w, obtained by Equation (14) was . The numerical solution, represented by the dotted curve (magenta), was obtained through the propagation of imaginary time using the Gaussian function as an initial condition. Here, labels “anltc”, “var”, and “num” represent analytical, variational, and numerical results, respectively.

Figure 2

Variational and numerical results of the Equation (16) for different values of the Pöschl–Teller potential amplitude. It is noticed that, as the depth of the potential increases, the wave function width decreases and its amplitude increases, preserving the normalization. The nonlinearity coefficient remained fixed () in all simulations.

The results shown in Figure 2 indicate that the widths of the wave functions decrease and their amplitudes increase as the potential depth increases.

The results for the values of the chemical potentials obtained by varying the amplitude of the potential can be interpreted through Figure 3. The values adopted for the amplitude are within the range . We can see that, as the depth of the potential increases, the chemical potential decreases through linear behavior. In this adopted interval, the results obtained from the chemical potential versus amplitude showed similar results observed in studies of Bose–Einstein condensates and Superfluid Fermi gases in optical lattices and optical superlattices (doubly periodic and quasiperiodic) [41,42,43,44,45].

Figure 3

Variational and numerical results of chemical potential, , versus amplitude, , of PT potential . We consider in this simulation . The variational results are represented by the dashed curve (blue) while the numerical results are represented by the dotted curve (magenta). Similar behavior was observed in studies of stable gap solitons in SFGs trapped by optical lattices and optical superllatices.

5. Conclusions

In this paper, we consider the scenario where a BEC made up of diluted ultracold atoms with attractive interatomic interactions is transversally confined by a harmonic potential and longitudinally by the PT potential, assuming a cigar-shaped configuration. The results obtained demonstrated that the PT potential was able to trap BEC. The solutions showed stability and confirmed a behavior similar to the bright gap-solitons obtained in BECs and SFGs trapped by optical lattices. These results may open the way for a better understanding of ultracold quantum gases trapped by the PT potential and motivate new studies that may have applications in other areas of physics such as solid state physics, condensed matter physics, and materials sciences.