Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Quantum Dynamics, Many-body methods and basis sets

Halvard Sutterud

Published Jan 1, 2020• 0 citations
OpenAlex

Abstract

Accurate descriptions of many-particle quantum systems subject to laser interactions can be found using real-time ab initio methods. Of these the arguably most popular and exact is the multi-configuration time-dependent Hartree-Fock (MCTDHF) method. However, MCTDHF suffers from computational limitations in that it quickly becomes too time consuming. The orbital-adaptive time-dependent coupled-cluster (OATDCC) method represents a hierarchy of approximations to MCTDHF that are less computationally expensive while retaining as much accuracy as possible. Building on an existing codebase we have in this thesis generalized the OATDCC method to include Q-space orbital equations. A novel ground state solver is implemented, employing adiabatic switching, since imaginary time propagation is not feasible "out of the box". Furthermore, we implement a sinc-discrete variable representation basis for one-dimensional model systems. We demonstrate that ionization and high-harmonic processes can be described using the OATDCCD method. Comparison with the more accurate, yet more expensive, multiconfigurational time-dependent Hartree-Fock method indicates that the OATDCCD method is an excellent approximation.

Related papers

AQUADUCT Ab Initio Quantum Dynamics Using Coupled Cluster in Time

Sebastian Gregorius Winther-Larsen

We advance the Coupled Cluster method’s in the time-dependent realm, by imple- menting a robust solver based on the orbital-adaptive time-dependent Coupled Cluster (OATDCC)[1] method. This involves implementing both a simplified static orbital time-dependent coupled cluster solver with single and double excitations (TDCCSD) and an orbital-adaptive scheme with double excitations (OATDCCD). To supplement the time-dependent methods we implement several ground state solvers based on the Lagrangian Coupled Cluster formulation, with single and double excitations, as well as a non-orthogonal orbital-optimised Coupled Cluster (NOCC) solver[2]. We construct several quantum dot basis sets with different potential functions in one- and two dimensions, including interactions with magnetic fields. What is more, we also implement an interface with popular quantum chemistry software modules PySCF[3] and Psi4[4] for extraction of additional basis sets for atoms and molecules. The quantum systems are allowed to vary with time by addition of a time-dependent term to the Hamiltonian, with which we simulate a laser field in the dipole approxima- tion. As a validation we reproduce results from the scientific literature, both for atoms, molecules and quantum dots. We show that our methods lead to convergence in the ever-increasing basis set size limit, for simple quantum dot systems. For the same quantum dot system, we show how sensitive a system is to changes in the frequency of a driving oscillating field. Frequencies closer to the resonant frequency lead to exctiations and increased energy. We are able to simulate systems that are fairly large - quantum dots in one- and two dimensions with up to twelve electrons. For systems that meander far from the reference state, we show that the orbital-adaptive method has far superior stability, compared with the method with static orbitals. For all quantum dot systems we find strong comformity with the harmonic potential theorem[5], yet we see a slight many-body effect for a two-dimensional double dot system. By subjecting the two-dimensional quantum dots to a homogoenous, static magnetic field in the form of an angular momentum operator, we see two frequencies in the dipole spectrum, instead of one frequency. This is also in accordance with the harmonic potential theorem. The difference between the two frequencies in the new spectrum is the same as the Larmor frequency of the magnetic field, within acceptable tolerance levels.

•Jan 1, 2019•1 citation

Colloquium: Multiconfigurational time-dependent Hartree approaches for indistinguishable particles

Axel U. J. Lode, Camille Lévêque, Lars Bojer Madsen +2 more

In this Colloquium, the wavefunction-based Multiconfigurational Time-Dependent Hartree approaches to the dynamics of indistinguishable particles (MCTDH-F for Fermions and MCTDH-B for Bosons) are reviewed. MCTDH-B and MCTDH-F or, together, MCTDH-X are methods for describing correlated quantum systems of identical particles by solving the time-dependent Schrödinger equation from first principles. MCTDH-X is used to accurately model the dynamics of real-world quantum many-body systems in atomic, molecular, and optical physics. The key feature of these approaches is the time-dependence and optimization of the single-particle states employed for the construction of a many-body basis set, which yields nonlinear working equations. We briefly describe the historical developments that have lead to the formulation of the MCTDH-X methods and motivate the necessity for wavefunction-based approaches. We sketch the derivation of the unified MCTDH-F and MCTDH-B equations of motion for complete and also specific restricted configuration spaces. The strengths and limitations of the MCTDH-X approach are assessed via benchmarks against an exactly solvable model and via convergence checks. We highlight some applications to instructive and experimentally-realized quantum many-body systems: the dynamics of atoms in Bose-Einstein condensates in magneto-optical and optical traps and of electrons in atoms and molecules. We discuss the current development and frontiers in the field of MCTDH-X: theories and numerical methods for indistinguishable particles, for mixtures of multiple species of indistinguishable particles, the inclusion of nuclear motion for the nonadia-batic dynamics of atomic and molecular systems, the so-called multilayer generalizations to the MCTDH-F and MCTDH-B methods, and the time-dependent orbital-adaptive coupled cluster theory are discussed.

•Nov 19, 2019•3 citations

Ab initio quantum dynamics using coupled-cluster

Simen Kvaal

The curse of dimensionality (COD) limits the current state-of-the-art ab initio propagation methods for non-relativistic quantum mechanics to relatively few particles. For stationary structure calculations, the coupled-cluster (CC) method overcomes the COD in the sense that the method scales polynomially with the number of particles while still being size-consistent and extensive. We generalize the CC method to the time domain while allowing the single-particle functions to vary in an adaptive fashion as well, thereby creating a highly flexible, polynomially scaling approximation to the time-dependent Schrödinger equation. The method inherits size-consistency and extensivity from the CC method. The method is dubbed orbital-adaptive time-dependent coupled-cluster, and is a hierarchy of approximations to the now standard multi-configurational time-dependent Hartree method for fermions. A numerical experiment is also given.

arXiv:1201.5548•The Journal of Chemical Physics•May 18, 2012•147 citations

Efficient computation of the MCTDHF approximation to the time-dependent Schrödinger equation

Othmar Koch

We discuss analytical and numerical properties of the multi-configuration time-dependent Hartree-Fock method for the approximate solution of the time-dependent multi-particle (electronic) Schrödinger equation which are relevant for an efficient implementation of this model reduction technique. Particularly, we focus on a discretization and low rank approximation in the evaluation of the meanfield terms occurring in the MCTDHF equations of motion, which is crucial for the computational tractability of the problem. We give error bounds for this approximation and demonstrate the achieved gain in performance.

Opuscula Mathematica•Jan 1, 2006•3 citations

One particle equations for many particle quantum systems: The MCTHDF method

Claude Bardos, Norbert J. Mauser

This contribution is devoted to the mathematical analysis of more or less sophisticated approximations of the time evolution of systems of $N$ quantum particles. New results for the Multiconfiguration Time Dependent Hartree-Fock (MCTDHF) method (which cover the material of the talk given by the first author at the “Nonlinear waves conference in honor of Walter Strauss”) are summarized and compared with the simpler Hartree and Hartree-Fock equations.

Quarterly of Applied Mathematics•Oct 19, 2009•1 citation

SOLVING THE SINGLE LAYER MULTICONFIGURATION TIME-DEPENDENT HARTREE–FOCK EQUATION IN SECOND QUANTIZATION REPRESENTATION BASED ON A DECOMPOSITION OF THE OVERALL FOCK SPACE

WENLIANG LI, Wenwu Xu, Keli Han

Very recently, the multilayer multiconfiguration time-dependent Hartree (MCTDH) method based on a decomposition of the overall Fock space in second quantization representation was proposed [Wang H, Thoss M, J Chem Phys131:024114, 2009], and they have presented demonstrative numerical example on vibrationally coupled electron transport. Followed by the thinking, the strategies of the implement of the theory to deal with multielectron dynamics with consideration coulomb interaction were discussed in detail. Some demonstrative calculations have been carried out by using the single layer MCTDHF method. The influences of the different choice way of initial condition, such as the species of basis function, the type of decomposition of the full Fock space, the choice of the spin orbitals, the number of the single particle function and so on, on the imaginary time propagation are mainly studied. The bridge between the generalized theory and the calculation application has been set up.

Journal of Theoretical and Computational Chemistry•Aug 26, 2012•4 citations
Accurate descriptions of many-particle quantum systems subject to laser interactions can be found using real-time ab initio methods. Of these the arguably most popular and exact is the multi-configuration time-dependent Hartree-Fock (MCTDHF) method. However, MCTDHF suffers from computational limitations in that it quickly becomes too time consuming. The orbital-adaptive time-dependent coupled-cluster (OATDCC) method represents a hierarchy of approximations to MCTDHF that are less computationally expensive while retaining as much accuracy as possible. Building on an existing codebase we have in this thesis generalized the OATDCC method to include Q-space orbital equations. A novel ground state solver is implemented, employing adiabatic switching, since imaginary time propagation is not feasible "out of the box". Furthermore, we implement a sinc-discrete variable representation basis for one-dimensional model systems. We demonstrate that ionization and high-harmonic processes can be described using the OATDCCD method. Comparison with the more accurate, yet more expensive, multiconfigurational time-dependent Hartree-Fock method indicates that the OATDCCD method is an excellent approximation.
We advance the Coupled Cluster method’s in the time-dependent realm, by imple- menting a robust solver based on the orbital-adaptive time-dependent Coupled Cluster (OATDCC)[1] method. This involves implementing both a simplified static orbital time-dependent coupled cluster solver with single and double excitations (TDCCSD) and an orbital-adaptive scheme with double excitations (OATDCCD). To supplement the time-dependent methods we implement several ground state solvers based on the Lagrangian Coupled Cluster formulation, with single and double excitations, as well as a non-orthogonal orbital-optimised Coupled Cluster (NOCC) solver[2]. We construct several quantum dot basis sets with different potential functions in one- and two dimensions, including interactions with magnetic fields. What is more, we also implement an interface with popular quantum chemistry software modules PySCF[3] and Psi4[4] for extraction of additional basis sets for atoms and molecules. The quantum systems are allowed to vary with time by addition of a time-dependent term to the Hamiltonian, with which we simulate a laser field in the dipole approxima- tion. As a validation we reproduce results from the scientific literature, both for atoms, molecules and quantum dots. We show that our methods lead to convergence in the ever-increasing basis set size limit, for simple quantum dot systems. For the same quantum dot system, we show how sensitive a system is to changes in the frequency of a driving oscillating field. Frequencies closer to the resonant frequency lead to exctiations and increased energy. We are able to simulate systems that are fairly large - quantum dots in one- and two dimensions with up to twelve electrons. For systems that meander far from the reference state, we show that the orbital-adaptive method has far superior stability, compared with the method with static orbitals. For all quantum dot systems we find strong comformity with the harmonic potential theorem[5], yet we see a slight many-body effect for a two-dimensional double dot system. By subjecting the two-dimensional quantum dots to a homogoenous, static magnetic field in the form of an angular momentum operator, we see two frequencies in the dipole spectrum, instead of one frequency. This is also in accordance with the harmonic potential theorem. The difference between the two frequencies in the new spectrum is the same as the Larmor frequency of the magnetic field, within acceptable tolerance levels.
In this Colloquium, the wavefunction-based Multiconfigurational Time-Dependent Hartree approaches to the dynamics of indistinguishable particles (MCTDH-F for Fermions and MCTDH-B for Bosons) are reviewed. MCTDH-B and MCTDH-F or, together, MCTDH-X are methods for describing correlated quantum systems of identical particles by solving the time-dependent Schrödinger equation from first principles. MCTDH-X is used to accurately model the dynamics of real-world quantum many-body systems in atomic, molecular, and optical physics. The key feature of these approaches is the time-dependence and optimization of the single-particle states employed for the construction of a many-body basis set, which yields nonlinear working equations. We briefly describe the historical developments that have lead to the formulation of the MCTDH-X methods and motivate the necessity for wavefunction-based approaches. We sketch the derivation of the unified MCTDH-F and MCTDH-B equations of motion for complete and also specific restricted configuration spaces. The strengths and limitations of the MCTDH-X approach are assessed via benchmarks against an exactly solvable model and via convergence checks. We highlight some applications to instructive and experimentally-realized quantum many-body systems: the dynamics of atoms in Bose-Einstein condensates in magneto-optical and optical traps and of electrons in atoms and molecules. We discuss the current development and frontiers in the field of MCTDH-X: theories and numerical methods for indistinguishable particles, for mixtures of multiple species of indistinguishable particles, the inclusion of nuclear motion for the nonadia-batic dynamics of atomic and molecular systems, the so-called multilayer generalizations to the MCTDH-F and MCTDH-B methods, and the time-dependent orbital-adaptive coupled cluster theory are discussed.
The curse of dimensionality (COD) limits the current state-of-the-art ab initio propagation methods for non-relativistic quantum mechanics to relatively few particles. For stationary structure calculations, the coupled-cluster (CC) method overcomes the COD in the sense that the method scales polynomially with the number of particles while still being size-consistent and extensive. We generalize the CC method to the time domain while allowing the single-particle functions to vary in an adaptive fashion as well, thereby creating a highly flexible, polynomially scaling approximation to the time-dependent Schrödinger equation. The method inherits size-consistency and extensivity from the CC method. The method is dubbed orbital-adaptive time-dependent coupled-cluster, and is a hierarchy of approximations to the now standard multi-configurational time-dependent Hartree method for fermions. A numerical experiment is also given.
We discuss analytical and numerical properties of the multi-configuration time-dependent Hartree-Fock method for the approximate solution of the time-dependent multi-particle (electronic) Schrödinger equation which are relevant for an efficient implementation of this model reduction technique. Particularly, we focus on a discretization and low rank approximation in the evaluation of the meanfield terms occurring in the MCTDHF equations of motion, which is crucial for the computational tractability of the problem. We give error bounds for this approximation and demonstrate the achieved gain in performance.
This contribution is devoted to the mathematical analysis of more or less sophisticated approximations of the time evolution of systems of NNN quantum particles. New results for the Multiconfiguration Time Dependent Hartree-Fock (MCTDHF) method (which cover the material of the talk given by the first author at the “Nonlinear waves conference in honor of Walter Strauss”) are summarized and compared with the simpler Hartree and Hartree-Fock equations.
Very recently, the multilayer multiconfiguration time-dependent Hartree (MCTDH) method based on a decomposition of the overall Fock space in second quantization representation was proposed [Wang H, Thoss M, J Chem Phys131:024114, 2009], and they have presented demonstrative numerical example on vibrationally coupled electron transport. Followed by the thinking, the strategies of the implement of the theory to deal with multielectron dynamics with consideration coulomb interaction were discussed in detail. Some demonstrative calculations have been carried out by using the single layer MCTDHF method. The influences of the different choice way of initial condition, such as the species of basis function, the type of decomposition of the full Fock space, the choice of the spin orbitals, the number of the single particle function and so on, on the imaginary time propagation are mainly studied. The bridge between the generalized theory and the calculation application has been set up.