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

arXiv:2407.04275 · The Journal of Chemical Physics

Time-dependent orbital-optimized coupled-cluster methods families for fermion-mixtures dynamics

Haifeng Lang, Takeshi Sato

Published Sep 18, 2024• 2 citations
arXivHTMLPDF

Abstract

Five time-dependent orbital optimized coupled-cluster methods, of which four can converge to the time-dependent complete active space self-consistent-field method, are presented for fermion-mixtures with arbitrary fermion kinds and numbers. Truncation schemes maintaining the intragroup orbital rotation invariance, as well as equations of motion of coupled-cluster (CC) amplitudes and orbitals, are derived. Present methods are compact CC-parameterization alternatives to the time-dependent multiconfiguration self-consistent-field method for systems consisting of arbitrarily different kinds and numbers of interacting fermions. Theoretical analysis of applications of present methods to various chemical systems is reported.

DOI
10.1063/5.0227236

Related papers

Non‐Bivariational Time‐Dependent Orbital‐Quasi‐Optimized Coupled‐Cluster Methods Family for Fermion‐Mixtures Dynamics

Haifeng Lang, Takeshi Sato

Recently, [JCP 161, 114114 (2024)] successfully derived numerically stable and converging (to the time-dependent complete active space self-consistent-field) time-dependent orthogonal coupled cluster (CC) methods, namely, time-dependent orbital optimized CC (TD-ooCC) methods for fermion-mixtures with arbitrary fermion kinds and numbers, via the time-dependent bivariational principle and setting constraints in the ansatz. In this work, we extend the TD-ooCC methods to eighteen non-bivariational time-dependent orbital-quasi-optimized coupled-cluster (TD-oqoCC-NB) methods. In TD-oqoCC-NB methods, the constraints in the ansatz are also used to avoid overparameterization in the complete active space limit and eliminate the numerical instability in the finite truncation case. The equations of motion are not determined by the time-dependent bivariational principle. Instead, all equations ensuring the convergence to the complete active space limit can be selected flexibly. The applications to electron dynamics and vibrational dynamics are also discussed. In particular, we find that all NB methods violate the Ehrenfest theorem even when they are gauge-invariant, which distinguishes them from TD-ooCC methods. The stationary condition of EOMs of each NB methods also provides novel CC methods for time-independent electronic structure calculations.

Journal of Computational Chemistry•Aug 5, 2026•0 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

An Alternative Algebraic Framework for the Simplification of Coupled Cluster Type Expectation Values

Daniel Pape

The Coupled-Cluster (CC) method is one of the most popular and efficient correlation methods in quantum chemistry. \nEspecially the CCSD(T) approximation, which includes single and double excitations by means of the application of the cluster operator \nto a reference determinant and triple excitations via a perturbative treatment, has become a standard tool \nin quantum chemical applications. However, the method is restricted to relatively small system sizes due to its unfavourable scaling (N^7, where N is the number of basis functions applied). \nFor the full treatment of triple excitations the scaling advances to N^8 and every further excitation level increases the exponent by two. \n \nThe goal of this work is to reduce the calculation time for closed shell systems at least by a factor growing with the excitation level for arbitrary truncation levels. \nThis is done by restricting the spatial parts of the spin orbitals and thus treat pairs of spin orbitals on the same footing. The restrictions can be easily constructed \nfor the CCSD model by spin integration. \nThe derivation of the restrictions arising in higher excited case will be done employing \nthe spatial orbital excitation operators E. \n \nIn the first part of this work an algorithm is derived that is capable of the derivation of the energy and amplitude equations for arbitrary excitation levels. In the second part \nan implementation of this algorithm is presented.

•Nov 12, 2013•0 citations

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

The extended coupled cluster method and the pairing problem

Chris Snape, Niels R. Walet

We study the application of various forms of the coupled cluster method to systems with paired fermions. The novel element of the analysis is the study of the breaking and eventual restoration of particle number in the CCM variants. We specifically include Arponen's extended coupled cluster method, which describes the normal Hartree-Fock-Bogoliubov mean field at lowest level of truncation. We show that all methods converge to the exact results as we increase the order of truncation, but that the breaking of particle number at an intermediate level means that this convergence occurs in a surprising way. We argue that the most straightforward form of the method seems to be the most stable approach to implement for realistic (large number of particles) pairing problems

arXiv:1011.5146•Nov 23, 2010•0 citations

Three Lagrangians for the complete-active space coupled-cluster method

Simen Kvaal

Three fully variational formulations of the complete-active space coupled-cluster method are derived. The formulations include the ability to approximate the model vectors by smooth manifolds, thereby opening up the possibility for overcoming the exponential wall of scaling for model spaces of complete-active space type. In particular, model vectors of matrix-product states are considered, and it is argued that the present variational formulation allows not only favorably scaling multireference coupled-cluster calculations but also systematic correction of tailored coupled-cluster calculations and of quantum chemical density-matrix renormalization group methods, which are fast and polynomial scaling but lack the ability to properly resolve dynamical correlation at chemical accuracy. The extension of the variational formulations to the time domain is also discussed, with derivations of abstract evolution equations.

arXiv:2205.08792•The Journal of Chemical Physics•Jun 26, 2023•3 citations
Five time-dependent orbital optimized coupled-cluster methods, of which four can converge to the time-dependent complete active space self-consistent-field method, are presented for fermion-mixtures with arbitrary fermion kinds and numbers. Truncation schemes maintaining the intragroup orbital rotation invariance, as well as equations of motion of coupled-cluster (CC) amplitudes and orbitals, are derived. Present methods are compact CC-parameterization alternatives to the time-dependent multiconfiguration self-consistent-field method for systems consisting of arbitrarily different kinds and numbers of interacting fermions. Theoretical analysis of applications of present methods to various chemical systems is reported.
Recently, [JCP 161, 114114 (2024)] successfully derived numerically stable and converging (to the time-dependent complete active space self-consistent-field) time-dependent orthogonal coupled cluster (CC) methods, namely, time-dependent orbital optimized CC (TD-ooCC) methods for fermion-mixtures with arbitrary fermion kinds and numbers, via the time-dependent bivariational principle and setting constraints in the ansatz. In this work, we extend the TD-ooCC methods to eighteen non-bivariational time-dependent orbital-quasi-optimized coupled-cluster (TD-oqoCC-NB) methods. In TD-oqoCC-NB methods, the constraints in the ansatz are also used to avoid overparameterization in the complete active space limit and eliminate the numerical instability in the finite truncation case. The equations of motion are not determined by the time-dependent bivariational principle. Instead, all equations ensuring the convergence to the complete active space limit can be selected flexibly. The applications to electron dynamics and vibrational dynamics are also discussed. In particular, we find that all NB methods violate the Ehrenfest theorem even when they are gauge-invariant, which distinguishes them from TD-ooCC methods. The stationary condition of EOMs of each NB methods also provides novel CC methods for time-independent electronic structure calculations.
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.
The Coupled-Cluster (CC) method is one of the most popular and efficient correlation methods in quantum chemistry. \nEspecially the CCSD(T) approximation, which includes single and double excitations by means of the application of the cluster operator \nto a reference determinant and triple excitations via a perturbative treatment, has become a standard tool \nin quantum chemical applications. However, the method is restricted to relatively small system sizes due to its unfavourable scaling (N^7, where N is the number of basis functions applied). \nFor the full treatment of triple excitations the scaling advances to N^8 and every further excitation level increases the exponent by two. \n \nThe goal of this work is to reduce the calculation time for closed shell systems at least by a factor growing with the excitation level for arbitrary truncation levels. \nThis is done by restricting the spatial parts of the spin orbitals and thus treat pairs of spin orbitals on the same footing. The restrictions can be easily constructed \nfor the CCSD model by spin integration. \nThe derivation of the restrictions arising in higher excited case will be done employing \nthe spatial orbital excitation operators E. \n \nIn the first part of this work an algorithm is derived that is capable of the derivation of the energy and amplitude equations for arbitrary excitation levels. In the second part \nan implementation of this algorithm is presented.
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.
We study the application of various forms of the coupled cluster method to systems with paired fermions. The novel element of the analysis is the study of the breaking and eventual restoration of particle number in the CCM variants. We specifically include Arponen's extended coupled cluster method, which describes the normal Hartree-Fock-Bogoliubov mean field at lowest level of truncation. We show that all methods converge to the exact results as we increase the order of truncation, but that the breaking of particle number at an intermediate level means that this convergence occurs in a surprising way. We argue that the most straightforward form of the method seems to be the most stable approach to implement for realistic (large number of particles) pairing problems
Three fully variational formulations of the complete-active space coupled-cluster method are derived. The formulations include the ability to approximate the model vectors by smooth manifolds, thereby opening up the possibility for overcoming the exponential wall of scaling for model spaces of complete-active space type. In particular, model vectors of matrix-product states are considered, and it is argued that the present variational formulation allows not only favorably scaling multireference coupled-cluster calculations but also systematic correction of tailored coupled-cluster calculations and of quantum chemical density-matrix renormalization group methods, which are fast and polynomial scaling but lack the ability to properly resolve dynamical correlation at chemical accuracy. The extension of the variational formulations to the time domain is also discussed, with derivations of abstract evolution equations.