External Trapping Potential Research Articles

Revised Enskog equation for hard rods

We point out that Percus’s collision integral for one-dimensional hard rods (Percus 1969 Phys. Fluids 12 1560–3) does not preserve the thermal equilibrium state in an external trapping potential. We derive a revised Enskog equation for hard rods and show that it preserves this thermal state exactly. In contrast to recent proposed kinetic equations for dynamics in integrability-breaking traps, both our kinetic equation and its thermal states are explicitly nonlocal in space. Our equation differs from earlier proposals at third order in spatial derivatives and we attribute this discrepancy to the choice of collision integral underlying our approach.

The application of the "inverse problem" method for constructing confining potentials that make N-soliton waveforms exact solutions in the Gross-Pitaevskii equation.

In this work, we discuss an application of the "inverse problem" method to find the external trapping potential, which has particular N trapped soliton-like solutions of the Gross-Pitaevskii equation (GPE) also known as the cubic nonlinear Schrödinger equation (NLSE). This inverse method assumes particular forms for the trapped soliton wave function, which then determines the (unique) external (confining) potential. The latter renders these assumed waveforms exact solutions of the GPE (NLSE) for both attractive (g<0) and repulsive (g>0) self-interactions. For both signs of g, we discuss the stability with respect to self-similar deformations and translations. For g<0, a critical mass Mc or equivalently the number of particles for instabilities to arise can often be found analytically. On the other hand, for the case with g>0 corresponding to repulsive self-interactions which is often discussed in the atomic physics realm of Bose-Einstein condensates, the bound solutions are found to be always stable. For g<0, we also determine the critical mass numerically by using linear stability or Bogoliubov-de Gennes analysis, and compare these results with our analytic estimates. Various analytic forms for the trapped N-soliton solutions in one, two, and three spatial dimensions are discussed, including sums of Gaussians or higher-order eigenfunctions of the harmonic oscillator Hamiltonian.

Synchronization Analysis of a Master-Slave BEC System via Active Control

This paper will focus on theoretical treatment of the dynamic of the Bose-Einstein Condensate (BEC) systems contained different external trapping potentials. We construct the phase space diagrams and Lyapunov Characteristic Exponents (LCEs) for master and slave systems depended on the system parameters and propose a nonlinear control for the synchronization of systems in their chaotic states. The synchronization is obtained in master-slave scheme for different initial values. Numerical results are also given to show the efficiency of the used control technique.

Baseband modulational instability and dynamics of rogue waves in coherently coupled Bose–Einstein condensates

We consider coupled dimensionless Gross-Pitaevskii (GP) equations that describes effective coherently coupled Bose–Einstein condensates (BECs) trapped in external complex potentials. Through a modified similarity transformation, the integrability condition is obtained and the system is reduced to one single integrable autonomous cubic-quintic nonlinear Schrödinger (NLS) equation with self-steepening and self-frequency shift. The baseband modulational instability (MI) for the obtained NLS equation is investigated in details. Through the exact first- and second-order rogue wave solutions with nonlinear chirp of the coupled GP equation, we investigate analytically the propagation dynamics of rogue waves on weakly perturbed continuous wave solutions in the parameter space where the baseband MI exist. Our results show that one can manipulate the motion of matter rogue waves by controlling the external trapping potential. Simulations provide supporting evidence that rogue waves can only emerge in regimes where baseband MI occurs. Our results show how multiplet first- and second-order rogue waves can be generated by varying some solution parameters such as the amplitude parameter.

Segregated solutions for some non-linear Schrödinger systems with critical growth

We find infinitely many non-radial solutions for a system of Schrödinger equations with critical growth in a fully attractive or repulsive regime in presence of an external radial trapping potential.

Staggered quantum phases of dipolar bosons at finite temperatures

The extended Bose-Hubbard model with correlated tunneling exhibits staggered superfluid and supersolid quantum phases. We study finite-temperature phase transitions of quantum phases of dipolar bosons in a two-dimensional optical lattice using Gutzwiller mean-field and quantum Monte Carlo approaches. When nearest-neighbor repulsion is comparable to the on-site interaction, we find that the two topologically distinct superfluids are separated by a normal-fluid phase, while at stronger off-site interactions, density-modulated insulating quantum phases appear. We estimate the critical temperature of the staggered-superfluid to normal-fluid transition and show that this transition is of the Kosterlitz-Thouless type. Finally, we elucidate the coexistence of staggered quantum phases in the presence of an external trapping potential. Our study paves the way to observe novel staggered quantum phases in dipolar optical lattice experiments.

Segregated solutions for nonlinear Schrödinger systems with weak interspecies forces

We find positive non-radial solutions for a system of Schrödinger equations in a weak fully attractive or repulsive regime in presence of an external radial trapping potential that exhibits a maximum or a minimum at infinity.

Nonequilibrium many-body dynamics in supersymmetric quenching

We study the dynamics induced by quenching an ultracold quantum many-body system between two supersymmetric Hamiltonians. Such a quench can be created by carefully changing the external trapping potential and leads to a situation where the eigenspectra before and after the quench are nearly identical. We show that the dynamics originating from this can be conveniently described using knowledge about the initial state only and apply this insight to the specific example of a fermionic gas that is initially trapped in an infinite box potential. When quenching to different, higher order supersymmetric partners potentials we observe the appearance of many-body revivals in the survival probability and show that some of these are robust at finite temperatures. This is in contrast to the well-known Talbot effect, which is the standard example for quenching into a system with a quadratic spectrum.

Quantum droplet of a two-component Bose gas in an optical lattice near the Mott insulator transition

We theoretically study dynamical formation of a quantum droplet in a two-component Bose-Hubbard system with an external trap potential. Specifically, the superfluid in the central region surrounded by the Mott insulator with double filling forms a quantum droplet, which is self-bound thanks to the discontinuous nature of the quantum phase transition between the two phases. We show how to induce the characteristic behavior of the droplet through the control of the trap potential by using the time-dependent Gutzwiller simulations in a two-dimensional system. The static and dynamical properties of the droplet can be described qualitatively by the effective Ginzburg-Landau field theory with cubic-quintic nonlinearities, where the attractive cubic nonlinearlity emerges although all the bare interparticle interactions are repulsive.

Vortex chains in rotating two-dimensional Bose-Einstein condensate in a harmonic plus optical lattices potential

&lt;sec&gt;Bose-Einstein condensate (BEC) is essentially a macroscopic quantum effect with quantum volatility, macroscopic quantum coherence and artificial controllability. Owing to its unique controllability, it becomes a new ideal platform for quantum simulations and studies of interacting quantum systems.&lt;/sec&gt;&lt;sec&gt;In this paper, the generation of vortices and the formation of vortex chains, as well as characteristics of vortex chains in rotating two-dimensional BEC in a potential composed of harmonic potential and optical lattice are studied numerically. Firstly, the generation of vortices, the formation and distribution of vortex chains and the effects of different physical parameters on the vortex chains in two-dimensional BEC are investigated by using the multigrid preconditioned conjugate gradient method. Secondly, the evolution of the vortex chains with time is studied by using the time-splitting spectral method. The results show that the generation of vortices in BEC trapped in the compound potential corresponds to the minimum value of the potential. When the depth of the optical lattice increases to a certain value, vortex chains are formed in the BEC. With the further increase of the depth of the optical lattice, the vortex depth in the vortex chain in the BEC decreases continuously, and finally the vortex chain disappears completely. When the interaction strength between atoms increases, the distribution range of the condensate expands, and the number of vortices and the number of vortex chains in the condensate also increase. When the interaction strength between atoms increases to a certain value, the symmetry of the vortex chains is broken. As the rotation frequency of the condensate increases, the distribution range of the condensate expands, and the number of vortices and the number of vortex chains in the condensate also increase. When the rotation frequency is close to the external trapping potential frequency, the linear alignment of the vortex chains is disrupted. It is also found that there are three stages in the evolution of the vortex chains in the BEC: in the first stage, vortex chains rotate together with the condensate, and the original chain distribution keeps unchanged; in the second stage, the phenomenon of vortex space extrusion appears, and the vortex chain is destroyed; in the third stage, the phenomenon of vortex space expansion occurs, and finally the vortex chains disappear. The results above show that the depth of the optical lattice, the interaction strength between atoms, and the rotation frequency of the condensate have important effects on the vortices and vortex chains in the condensate. By adjusting these physical quantities, the number of vortices and the shape of vortex chains in the BEC can be effectively manipulated. This may provide some theoretical reference and guidance for future experiments and applications.&lt;/sec&gt;

Quasiparticle kinetic theory for Calogero models

We show that the quasiparticle kinetic theory for quantum and classical Calogero models reduces to the free-streaming Boltzmann equation. We reconcile this simple emergent behaviour with the strongly interacting character of the model by developing a Bethe–Lax correspondence in the classical case. This demonstrates explicitly that the freely propagating degrees of freedom are not bare particles, but rather quasiparticles corresponding to eigenvectors of the Lax matrix. We apply the resulting kinetic theory to classical Calogero particles in external trapping potentials and find excellent agreement with numerical simulations in all cases, both for harmonic traps that preserve integrability and exhibit perfect revivals, and for anharmonic traps that break microscopic integrability. Our framework also yields a simple description of multi-soliton solutions in a harmonic trap, with solitons corresponding to sharp peaks in the quasiparticle density. Extensions to quantum systems of Calogero particles are discussed.

The mean-field limit of quantum Bose gases at positive temperature

We prove that the grand canonical Gibbs state of an interacting quantum Bose gas converges to the Gibbs measure of a nonlinear Schrödinger equation in the mean-field limit, where the density of the gas becomes large and the interaction strength is proportional to the inverse density. Our results hold in dimensions d ⩽ 3 d \leqslant 3 . For d &gt; 1 d &gt; 1 the Gibbs measure is supported on distributions of negative regularity and we have to renormalize the interaction. More precisely, we prove the convergence of the relative partition function and of the reduced density matrices in the L r L^r -norm with optimal exponent r r . Moreover, we prove the convergence in the L ∞ L^\infty -norm of Wick-ordered reduced density matrices, which allows us to control correlations of Wick-ordered particle densities as well as the asymptotic distribution of the particle number. Our proof is based on a functional integral representation of the grand canonical Gibbs state, in which convergence to the mean-field limit follows formally from an infinite-dimensional stationary phase argument for ill-defined non-Gaussian measures. We make this argument rigorous by introducing a white-noise-type auxiliary field, through which the functional integral is expressed in terms of propagators of heat equations driven by time-dependent periodic random potentials and can, in turn, be expressed as a gas of interacting Brownian loops and paths. When the gas is confined by an external trapping potential, we control the decay of the reduced density matrices using excursion probabilities of Brownian bridges.

Collective dipole oscillations in a bosonic ladder lattice with effective magnetic flux

In this work we study the model of flux lattice, specifically a two-leg Bose–Hubbard ladder lattice subject to effective magnetic flux, which can be implemented using spin–orbit coupled cold atoms trapped in an optical lattice. This model has attracted much recent research interest due to that it provides a platform to study the interplay between artificial gauge fields and interactions. We consider that collective dipole oscillations are activated in this model via quenching an external harmonic trapping potential and show that this system can exhibit unique dynamical behaviors in terms that oscillation amplitude and period critically dependent on the spin–orbit coupling. A dynamical slowing down in collective dipole oscillations will take place near the transition point between the Meissner superfuild and vortex superfluid phases, which agrees quantitatively with a variational analysis. We show that such a quantum phase transition can also be probed via the damping of the collective dipole oscillation. When the quenching displacement is large, the collective center-of-mass motion behaves as damped oscillation and the damping reaches its maximum near the phase transition point. The specific dynamical phenomena are explained with eigenspectra and level populations after the quenching. In the presence of interactions, the time-dependent density–matrix renormalization group (tDMRG) method is hired to study the collective dynamics. The case with weak interactions is similar to that without interactions, except that interactions tend to enhance damping. In addition to that, Bragg reflection can be observed when the quenching displacement exceeds a critical value. The case with strong interactions is also studied, in which case the phase transition points will be dramatically changed and the system can even move into the Mott insulator phase, which will greatly affect the corresponding collective dynamics. The phenomena predicted in this work can be readily observed in current available experiments on atom flux lattices.

Oscillating Quantum Droplets From the Free Expansion of Logarithmic One-dimensional Bose Gases

We analyse some issues related to the stability and free expansion of a one-dimensional logarithmic Bose–Einstein condensate, particularly its eventual relation to the formation of quantum droplet-type configurations. We prove that the corresponding properties, such as the energy of the associated N-body ground state, differ substantially with respect to its three-dimensional counterpart. Consequently, the free velocity expansion also shows remarkable differences with respect to the three-dimensional system when logarithmic interactions are taken into account. The one-dimensional logarithmic condensate tends to form quantum droplet-type configurations when the external trapping potential is turned off, i.e., the self-sustainability or self-confinement appears as in three-dimensions. However, we obtain that for some specific values of the self-interaction parameters and the number of particles under consideration, the cloud oscillates during the free expansion around to a specific equilibrium size. These results show that we are able to describe scenarios in which the one-dimensional cloud reaches stable configurations, i.e., oscillating quantum droplets.

Compensation process and generation of chirped femtosecond solitons and double-kink solitons in Bose–Einstein condensates with time-dependent atomic scattering length in a time-varying complex potential

We consider the one-dimensional (1D) cubic-quintic Gross–Pitaevskii (GP) equation, which governs the dynamics of Bose–Einstein condensate matter waves with time-varying scattering length and loss/gain of atoms in a harmonic trapping potential. We derive the integrability conditions and the compensation condition for the 1D GP equation and obtain, with the help of a cubic-quintic nonlinear Schrödinger equation with self-steepening and self-frequency shift, exact analytical solitonlike solutions with the corresponding frequency chirp which describe the dynamics of femtosecond solitons and double-kink solitons propagating on a vanishing background. Our investigation shows that under the compensation condition, the matter wave solitons maintain a constant amplitude, the amplitude of the frequency chirp depends on the scattering length, while the motion of both the matter wave solitons and the corresponding chirp depend on the external trapping potential. More interesting, the frequency chirps are localized and their feature depends on the sign of the self-steepening parameter. Our study also shows that our exact solutions can be used to describe the compression of matter wave solitons when the absolute value of the s-wave scattering length increases with time.

A chain of three dimensional scatterers under quasi-one dimensional quantum confinement: a Green’s function approach to energy band engineering

A mechanism to modify the energy band structure is proposed by considering a chain of periodic scatterers forming a linear lattice (akin to an optical lattice for ultracold atoms) around which an external cylindrical trapping potential is applied along the chain axis. It is analytically shown that, when this trapping (confining) potential is tight enough so as to reduce the dimensionality of the problem to an effective one dimension, it may modify the bound and scattering states of the local lattice potential, whose three-dimensional nature around each site is fully taken into account and going beyond the zero-range contact-potential approximation. Since these states contribute to the formation of the energy bands, the latter could thereby be continuously tuned by manipulating the confinement without the need to change the lattice potential. Such dimensionality reduction by quantum confinement is analyzed by using a Green’s functions method that can describe the scattering off the lattice sites in the presence of strong confinement and that can collect the contribuitions of several scattering partial waves.

Thermalization of a Trapped One-Dimensional Bose Gas via Diffusion.

For a decade the fate of a one-dimensional gas of interacting bosons in an external trapping potential remained mysterious. We here show that whenever the underlying integrability of the gas is broken by the presence of the external potential, the inevitable diffusive rearrangements between the quasiparticles, quantified by the diffusion constants of the gas, eventually lead the system to thermalize at late times. We show that the full thermalizing dynamics can be described by the generalized hydrodynamics with diffusion and force terms, and we compare these predictions to numerical simulations. Finally, we provide an explanation for the slow thermalization rates observed in numerical and experimental settings: the hydrodynamics of integrable models is characterized by a continuity of modes, which can have arbitrarily small diffusion coefficients. As a consequence, the approach to thermalization can display prethermal plateau and relaxation dynamics with long polynomial finite-time corrections.

Three-neutron bound and continuum states

The three-neutron ($3n$) system is studied by numerical calculations with the Faddeev three-body formalism for a realistic nucleon-nucleon (NN) potential. A response function for the transition from $^{3}\mathrm{H}$ to $3n$ continuum states by an isospin excitation operator is calculated, from which no evidence of $3n$ resonance state is found. Different methods to extrapolate the $3n$ energy from bound state energies with an extra attractive effect to the NN potential are examined. While extrapolations with attractive effects by enhanced NN potentials or three-body potentials result in the nonexistence of $3n$ resonance states, one by external trapping potentials leads to a positive $3n$ energy, which may be considered as a resonance state. It is found that this contradiction is due to a general defect of the trapping method.

Numerical implementation of a Mach-Zehnder interferometer for Bose-Einstein condensates

We numerically implement a Mach-Zehnder interferometer, where the coherence and oscillatory properties of Bose-Einstein condensates are explored and the system is modeled by the Gross-Pitaevskii equation. Several time-dependent external trapping potentials were engineered seeking the adiabatic regime which is quantified using fidelity measurements with respect to the actual ground-state of the trap. The dynamics of both conjugate variables, namely density and phase of the matter-wave function, are shown. Moreover, the density and fidelity profiles of the system are presented when the phase-shifter is switching-on and -off, being found in the presented profiles that the system exhibits three different regimes during the recombination stage among them even an orthogonal BEC to the original one is obtained. We achieve the numerical solution through an adequate implementation of the finite-difference method for the spatial discretization and a Runge-Kutta method for the time evolution.

Rogue waves of the dissipative Gross–Pitaevskii equation with distributed coefficients

By introducing a suitable ansatz and employing the similarity transformation technique, we construct the first- and second-order rational solutions for a quasi-one-dimensional (1D) dissipative Gross–Pitaevskii (GP) equation with a time-varying cubic nonlinearity and an external time-dependent potential. Then, by using these solutions, we engineer first- and second-order rogue waves in the Bose–Einstein condensate (BEC) contexts for the experimentally relevant systems when the gain/loss of atoms is taken into consideration. Our analysis shows that the control of the scattering length, the external harmonic, and the linear trapping potentials allows one to manage the motion and the background of dissipative rogue matter waves in BEC systems. We show that the wave amplitudes depend on the absolute value of s-wave scattering and the bias magnetic field, while its motion depends on the external trapping potentials. We show that unlike classical rogue waves, the nonzero continuous wave backgrounds of non-autonomous forced (damped) rogue matter waves in BECs with time-dependent complicated potential increases (decreases) during the wave motion. Our results also reveal that neither the gain nor the loss of the BEC atoms affects the amplitude of the rogue matter waves during their propagation. Our results may help to control and manage experimentally dissipative rogue waves in a BEC systems.