Computer simulations of many-body quantum dynamics of indistinguishable particles is a challenging task for computational physics. In this paper we demonstrate that the method of coupled coherent states (CCS) developed previously for multidimensional quantum dynamics of distinguishable particles can be used to study indistinguishable bosons in the second-quantization formalism. To prove its validity, the technique termed here coupled coherent states for indistinguishable bosons (CCSB) is tested on two model problems. The first is a system-bath problem consisting of a tunneling mode coupled to a harmonic bath, previously studied by CCS and other methods in distinguishable representation in 20 dimensions. The harmonic bath is comprised of identical oscillators, and may be second quantized for use with CCSB, so that this problem may be thought of as a bosonic bath with an impurity. The cross-correlation function for the dynamics of the system and Fourier transform spectrum compare extremely well with a benchmark calculation, which none of the prior methods of studying the problem achieved. The second model problem involves 100 bosons in a shifted harmonic trap. Breathing oscillations in the one-body density are calculated and shown to compare favorably to a multiconfigurational time-dependent Hartree for bosons calculation, demonstrating the applicability of the method as a new formally exact way to study the quantum dynamics of Bose-Einstein condensates.
Computer simulations of many-body quantum dynamics of indistinguishable particles is a challenging task for computational physics. In this paper we demonstrate that the method of coupled coherent states (CCS) developed previously for multidimensional quantum dynamics of distinguishable particles can be used to study indistinguishable bosons in the second-quantization formalism. To prove its validity, the technique termed here coupled coherent states for indistinguishable bosons (CCSB) is tested on two model problems. The first is a system-bath problem consisting of a tunneling mode coupled to a harmonic bath, previously studied by CCS and other methods in distinguishable representation in 20 dimensions. The harmonic bath is comprised of identical oscillators, and may be second quantized for use with CCSB, so that this problem may be thought of as a bosonic bath with an impurity. The cross-correlation function for the dynamics of the system and Fourier transform spectrum compare extremely well with a benchmark calculation, which none of the prior methods of studying the problem achieved. The second model problem involves 100 bosons in a shifted harmonic trap. Breathing oscillations in the one-body density are calculated and shown to compare favorably to a multiconfigurational time-dependent Hartree for bosons calculation, demonstrating the applicability of the method as a new formally exact way to study the quantum dynamics of Bose-Einstein condensates.
Computer simulations of many-body quantum dynamics of indistinguishable particles is a challenging task for computational physics. In this paper we demonstrate that the method of coupled coherent states (CCS) developed previously for multidimensional quantum dynamics can be used to study indistinguishable bosons in the second quantisation formalism. To prove its validity, the technique termed here coupled coherent states for indistinguishable bosons (CCSB) is tested on two model problems. The first is a system-bath problem consisting of a tunnelling mode coupled to a harmonic bath, previously studied by CCS and other methods in distinguishable representation in 20 dimensions. The harmonic bath is comprised of identical oscillators, and may be second quantised for use with CCSB. The cross-correlation function for the dynamics of the system and Fourier transform spectrum compare extremely well with a benchmark calculation, and are in much better agreement than prior methods of studying the problem. The second model problem involves 100 bosons in a shifted harmonic trap. Oscillations in the 1-body density are calculated and shown to compare favourably to a multiconfigurational time-dependent Hartree for bosons calculation, demonstrating the applicability of the method as a new formally exact way to study the quantum dynamics of Bose-Einstein condensates.
The coupled coherent states method has demonstrated itself as an accurate and efficient method of studying the quantum dynamics of various systems. In recent years, its applicability has been extended by incorporating a number of new numerical expansions and modifications to generate a closely related family of methods. In this thesis, two new augmentations are developed to further broaden the scope of problems that are able to be treated. The first of these is a 2-layer extension of coupled coherent states, capable of providing an increased mathematical description of a degree or degrees of freedom within a quantum mechanical system, as well as beneficial numerical and scalability properties. The newly developed method is tested on a model system-bath Hamiltonian consisting of a tunnelling mode governed by an asymmetric double well potential coupled to a harmonic bath. It is found to compare well to previous methods of studying the Hamiltonian, as well as a benchmark calculation on the system conducted in this thesis, and demonstrate the beneficial numerical and scalability properties expected. The second development is to extend coupled coherent states to treat systems of indistinguishable bosons in the second quantisation representation. The method is tested on the same Hamiltonian as the 2-layer coupled coherent states scheme, where the harmonic bath is second quantised as it is comprised of oscillators of the same frequency, so they may be thought of as indistinguishable. Exploiting this symmetry property is found to be extremely advantageous, with remarkable agreement to the benchmark calculation. The method is then tested on a model Hamiltonian consisting of 100 bosons in a shifted harmonic trap, with oscillations in the 1-body density calculated. The results are found to compare favourably with a multiconfigurational time-dependent Hartree for bosons calculation that is equivalent to the Gross-Pitaevskii equation, providing impetus for future studies on systems of Bose-Einstein condensates. The existing ab initio multiple cloning extension of coupled coherent states for nonadiabatic dynamics is also used to study the ultrafast photodissociation of 2-ethylpyrrole. The results are compared to experimental data, and a novel insight into the dissociation mechanism is obtained, with it shown to be composed of a two step process. Firstly, molecules that are able to dissociate immediately over the barrier along the N-H coordinate do so in < 50 fs, and this is followed by a second slower dissociation process from molecules that must sample the potential energy surface before finding a way around the barrier. This is not observed experimentally due to the temporal widths of the laser pulses obscuring the dynamics in the < 50 fs window.
A new theory is proposed to accurately simulate quantum dynamics in systems of identical particles. It is based on the second quantization formalism of many-body quantum theory, in which the Fock space is represented by occupation-number states. Within this representation the overall Fock space can be formally decomposed into smaller subspaces, and the wave function can be expressed as a multilayer multiconfiguration Hartree expansion involving subvectors in these subspaces. The theory unifies the multilayer multiconfiguration time-dependent Hartree theory for both distinguishable and indistinguishable particles. Specific formulations are given for systems of identical fermions, bosons, and combinations thereof. Practical implementations are discussed, especially for the case of fermions, to include the operator algebra that enforces the symmetry of identical particles. The theory is illustrated by a numerical example on vibrationally coupled electron transport.
The exactly solvable quantum many-particle model with harmonic one- and two-particle interaction terms is extended to include time dependency. We show that when the external trap potential and interparticle interaction have a time dependency, the numerically exact solutions of the corresponding time-dependent many-boson Schr\"odinger equation are still available. We use these exact solutions to benchmark the recently developed multiconfigurational time-dependent Hartree method for bosons (MCTDHB) [Phys. Rev. Lett. 99, 030402 (2007); Phys. Rev. A 77, 033613 (2008)]. In particular, we benchmark the MCTDHB method for (i) the ground state; (ii) the breathing many-body dynamics activated by a quench scenario where the interparticle interaction strength is suddenly turned on to a finite value; (iii) the nonequilibrium dynamic for driven scenarios where both the trap- and interparticle-interaction potentials are time-dependent. Excellent convergence of the ground state and dynamics is demonstrated. The great relevance of the self-consistency and time adaptivity, which are the intrinsic features of the MCTDHB method, is demonstrated by contrasting the MCTDHB predictions and those obtained within the standard full configuration interaction method spanning a Fock space of the same size, but utilizing as one-particle basis set the fixed-shape eigenstates of the one-particle potential. Connections of the model's results to ultracold Bose-Einstein condensed systems are addressed.
Following a ``bottom-up approach'' in understanding many-particle effects and dynamics we provide a systematic ab initio study of the dependence of the breathing dynamics of ultracold bosons in a one-dimensional (1D) harmonic trap on the number of bosons ranging from few to many. To this end, we employ the multilayer multiconfiguration time-dependent Hartree method for bosons (ML-MCTDHB) which has been developed very recently [Kr\"onke, Cao, Vendrell, and Schmelcher, New J. Phys. 15, 063018 (2013)]. The beating behavior for two bosons is found numerically and consequently explained by an analytical approach. Drawing on this, we show how to compute the complete breathing mode spectrum in this case. We examine how the two-mode breathing behavior of two bosons evolves to the single-frequency behavior of the many-particle limit when adding more particles. In the limit of many particles, we numerically study the dependence of the breathing mode frequency on both the interaction strength as well as on the particle number. We provide an estimate for the parameter region where the mean-field description provides a valid approximation.
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.