We present a systematic numerical iteration approach to study the evolution properties of the spin-boson systems, which works well in whole coupling regime. This approach involves the evaluation of a set of coefficients for the formal expansion of the time-dependent Schrödinger equation
∣t⟩=e−iH^t∣t=0⟩ by expanding the initial state
∣t=0⟩ in Fock space. The main advantage of this method is that this set of coefficients is unique for the Hamiltonian studied, which allows one to calculate the time evolution based on the different initial states. To complement our numerical calculations, the method is applied to the Buck-Sukumar model. Furthermore, we pointed out that, when the ground state energy of the model is unbounded and no ground state exists in a certain parameter space, the unstable time evolution of the physical quantities is the natural results. Furthermore, we test the performance of the numerical method to the Hamiltonian that use anti-Hermitian terms for modeling open quantum systems.
In this paper, we systematically review mathematical models, theories and numerical methods for ground states and dynamics of spinor Bose-Einstein condensates (BECs) based on the coupled Gross-Pitaevskii equations (GPEs). We start with a pseudo spin-1/2 BEC system with/without an internal atomic Josephson junction and spin-orbit coupling including (i) existence and uniqueness as well as non-existence of ground states under different parameter regimes, (ii) ground state structures under different limiting parameter regimes, (iii) dynamical properties, and (iv) efficient and accurate numerical methods for computing ground states and dynamics. Then we extend these results to spin-1 BEC and spin-2 BEC. Finally, extensions to dipolar spinor systems and/or general spin-F (F>=3) BEC are discussed.
The feasibility of using a basis set approach to the study of quantum dissipative dynamics is investigated for the spin-boson model, a system of two discrete states linearly coupled to a harmonic bath. The infinite Hamiltonian is discretized to a finite number of degrees of freedom. Traditional basis set approach, in a multiconfiguration time-dependent Hartree context, is used to solve the time-dependent Schrödinger equations by explicitly including all the degrees of freedom (“system”+“bath”). Quantities such as the reduced density matrix are then evaluated via a quadrature summation/Monte Carlo procedure over a certain number of time-dependent wave functions. Numerically exact results are obtained by systematically increasing the number of bath modes used to represent the condensed phase environment, as well as other variational parameters (number of basis functions, configurations, etc.). The potential of the current method is briefly discussed.
In this work we define, analyze, and compare different numerical schemes that can be used to study the ground state properties of Bose-Fermi systems, such as mixtures of different atomic species under external forces or self-bound quantum droplets. The bosonic atoms are assumed to be condensed and are described by the generalized Gross-Pitaevskii equation. The fermionic atoms, on the other hand, are treated individually, and each atom is associated with a wave function whose evolution follows the Hartree-Fock equation. We solve such a formulated set of equations using a variety of methods, including those based on adiabatic switching of interactions and the imaginary time propagation technique combined with the Gram-Schmidt orthonormalization or the diagonalization of the Hamiltonian matrix. We show how different algorithms compete at the numerical level by studying the mixture in the range of parameters covering the formation of self-bound quantum Bose-Fermi droplets. • Split operator method for the system of condensed bosons and degenerate fermions. • Ground state of the Bose-Fermi mixture by the Gram-Schmidt orthonormalization. • Ground state of the Bose-Fermi mixture by diagonalizing the Hamiltonian matrix. • Formation of Bose-Fermi droplets using the split operator method.
A simple model of a boson system is considered for which the exact Schr\"odinger equation reduces to a difference equation which can be solved numerically. Application of the Bogoliubov method and comparison of its results with the numerical results then yield some insight into the validity of the Bogoliubov approximations. The simple model considered has, in the absence of interactions, only a zero energy state and two states of unit energy available to each boson. Initially there is assumed to be no interaction between bosons in the degenerate excited states and all existing interactions are assumed to be repulsive. The ground-state energy is calculated in the Bogoliubov approximation, with and without the depletion effect considered, and in another numerical approximation. These results and the results of the numerical solution of the exact Schr\"odinger equation converge to the Bogoliubov result with depletion ignored when the number of particles in the system approaches infinity. It is surprising to note that, for the interaction strengths considered, the Bogoliubov result is within 3% of the numerical result for as few as 32 particles in the system, and within 10% for only 4 particles. A modification of the above system is considered in which there exists an additional two-body interaction between particles in the degenerate excited states which may be attractive or repulsive. It is shown that the ground-state energy, with this additional interaction present, deviates from the Bogoliubov value linearly with the strength of the added interaction, i.e., $W={W}_{\mathrm{Bog}}+\ensuremath{\alpha}F$ (with
F the strength parameter of the interaction). The factor $\ensuremath{\alpha}$ is found to be small and constant over a wide range of
F and in this range the deviation of
W from the Bogoliubov value is small. There is, however, a transition point
FT, beyond which $\ensuremath{\alpha}$ is constant and large and the usual Bogoliubov approximation is invalid for interactions more attractive than that characterized by
FT. A Bogoliubov-like approximation is shown to be quite accurate in this region. It is then shown that $F_{T}^{\ensuremath{-}1}$ is a linear function of the number of particles in the system.
The evolution of Bose-Einstein condensates is amply described by the time-dependent Gross-Pitaevskii mean-field theory which assumes all bosons to reside in a single time-dependent one-particle state throughout the propagation process. In this work, we go beyond mean field and develop an essentially exact many-body theory for the propagation of the time-dependent Schr\"odinger equation of
N interacting identical bosons. In our theory, the time-dependent many-boson wave function is written as a sum of permanents assembled from orthogonal one-particle functions, or orbitals, where both the expansion coefficients and the permanents (orbitals) themselves are time-dependent and fully determined according to a standard time-dependent variational principle. By employing either the usual Lagrangian formulation or the Dirac-Frenkel variational principle we arrive at two sets of coupled equations of motion, one for the orbitals and one for the expansion coefficients. The first set comprises of first-order differential equations in time and nonlinear integrodifferential equations in position space, whereas the second set consists of first-order differential equations with time-dependent coefficients. We call our theory multiconfigurational time-dependent Hartree for bosons, or
MCTDHB(M), where
M specifies the number of time-dependent orbitals used to construct the permanents. Numerical implementation of the theory is reported and illustrative numerical examples of many-body dynamics of trapped Bose-Einstein condensates are provided and discussed. The convergence of the method with a growing number
M of orbitals is demonstrated in a specific example of four interacting bosons in a double well.