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

The Journal of Chemical Physics

Time-dependent vibrational coupled cluster theory: Theory and implementation at the two-mode coupling level

Mads Bøttger Hansen, Niels Kristian Madsen, Alberto Zoccante, Ove Christiansen A

Published Oct 21, 2019• 25 citations
Publisher

Abstract

Equations are derived for the time evolution of time-dependent vibrational coupled cluster (TDVCC) wave functions covering both the TDVCC ket state and the associated so-called Λ bra state. The equations are implemented in the special case of both the Hamiltonian and the cluster operator containing at most two-mode coupling terms. The nontrivial behavior of the evolution of norm, energy, and expectation values due to the nonunitary time-evolution of the nonvariational TDVCC theory is analyzed theoretically and confirmed in numerical experiments that also include time-dependent Hamiltonians. In the spirit of time-independent size-consistency analysis, the separability of both the coupled cluster and Λ states for noninteracting systems is studied. While the coupled cluster state clearly has the correct behavior, the behavior of the Λ state is more intricate, and the consequence for different properties is shown theoretically and numerically. Overall, the numerical experiments show that TDVCC in incomplete expansions gives higher accuracy than a standard linear variational wave function parameterization with the same number of independent parameters, while equivalent results are obtained for complete expansions. The efficiency of the methodology is illustrated in computations on polycyclic aromatic hydrocarbons with up to 156 modes.

DOI
10.1063/1.5117207

Related papers

Time-dependent vibrational coupled cluster with variationally optimized time-dependent basis sets

Niels Kristian Madsen, Mads Bøttger Hansen, Ove Christiansen A +1 more

We develop time-dependent vibrational coupled cluster with time-dependent modals (TDMVCC), where an active set of one-mode basis functions (modals) is evolved in time alongside coupled-cluster wave-function parameters. A biorthogonal second quantization formulation of many-mode dynamics is introduced, allowing separate biorthogonal bases for the bra and ket states, thus ensuring complex analyticity. We employ the time-dependent bivariational principle to derive equations of motion for both the one-mode basis functions and the parameters describing the cluster (T) and linear de-excitation (L) operators. The choice of constraint (or gauge) operators for the modal time evolution is discussed. In the case of untruncated cluster expansion, the result is independent of this choice, but restricting the excitation space removes this invariance; equations for the variational determination of the constraint operators are derived for the latter case. We show that all single-excitation parts of T and L are redundant and can be left out in the case of variationally determined constraint-operator evolution. Based on a pilot implementation, test computations on Henon-Heiles model systems, the water molecule, and a reduced-dimensionality bi-thiophene model are presented, showing highly encouraging results for TDMVCC. It is demonstrated how TDMVCC in the limit of a complete cluster expansion becomes equivalent to multiconfiguration time-dependent Hartree for the same active-space size. Similarly, it is discussed how TDMVCC generally gives better and more stable results than its time-independent-modals counterpart, while equivalent results are obtained for complete expansions and full one-mode basis sets.

The Journal of Chemical Physics•Nov 3, 2020•16 citations

Time‐Dependent Vibrational Coupled Cluster Theory With Static and Dynamic Basis Functions

Mads Greisen Højlund, Alberto Zoccante, Andreas Buchgraitz Jensen +1 more

ABSTRACT In recent decades, coupled cluster theory has proven valuable in accurately describing correlation in many‐body systems, particularly in time‐independent computations of molecular electronic structure and vibrations. This review describes recent advancements in using coupled cluster parameterizations for time‐dependent wave functions for the efficient computation of the quantum dynamics associated with the motion of nuclei. It covers time‐dependent vibrational coupled cluster (TDVCC) and time‐dependent modal vibrational coupled cluster (TDMVCC), which employ static and adaptive basis sets, respectively. We discuss the theoretical foundation, including many‐mode second quantization, bivariational principles, and various parameterizations of time‐dependent bases. Additionally, we highlight key features that make TDMVCC promising for future quantum dynamical simulations. These features include fast configuration‐space convergence, the use of a compact adaptive basis set, and the possibility of efficient implementations with a computational cost that scales only polynomially with system size.

Wiley Interdisciplinary Reviews Computational Molecular Science•Jan 1, 2025•2 citations

Exponential parameterization of wave functions for quantum dynamics: Time-dependent Hartree in second quantization

Niels Kristian Madsen, Mads Bøttger Hansen, Alberto Zoccante +3 more

We derive equations for describing the time evolution of variational wave functions in linear and exponential parameterization with a second-quantization (SQ) formulation. The SQ formalism covers time-dependent Hartree (TDH), while exact states and approximate vibrational configuration interaction wave functions are described using state-transfer operators. We present detailed expressions for efficient evaluation of TDH in linear (L-TDH) and exponential (X-TDH) parametrization and an efficient implementation supporting linear scaling with respect to the number of degrees of freedom M when the Hamiltonian operator contains a constant number of terms per mode independently of the size of the system. The computational cost of the X-TDH method is reduced significantly compared to the L-TDH method for systems with many operator terms per mode such as is typical for accurate molecular potential-energy surfaces. Numerical results for L-TDH and X-TDH are presented which confirm the theoretical reduction of the M scaling compared to standard first-quantization formulations. Calculations on Henon-Heiles potentials with more than 105 dimensions and polycyclic aromatic hydrocarbons with up to 264 modes have been performed. Thus, the SQ formulation and the X-TDH method pave the way for studying the time-resolved quantum dynamics of large molecules.

The Journal of Chemical Physics•Oct 4, 2018•13 citations

Adaptive time-dependent coupled cluster method for multiparticle wavepacket dynamics

Simen Kvaal

A dynamical formulation of coupled cluster theory is derived using a variational principle. By allowing time-dependent single-particle functions, a high degree of adaptivity is introduced, allowing complex systems to be simulated with high accuracy. Equations of motion are derived which are shown to be suitable for computer implementation. The method, called adaptive time-dependent coupled cluster, is a strict generalization of the formulation used in standard coupled cluster response theory, and also represents a systematic hierarchy of size-consistent approximations, including standard time-dependent Hartree--Fock as a trivial case.

arXiv:1008.3835•Aug 23, 2010•0 citations

Comparing real-time coupled-cluster methods through simulation of collective Rabi oscillations

Andreas S. Skeidsvoll, Henrik Koch A

The time-dependent equation-of-motion coupled-cluster singles and doubles (TD-EOM-CCSD) and time-dependent coupled-cluster singles and doubles (TDCCSD) methods are compared by simulating Rabi oscillations for distant atoms in a classical electromagnetic field. While the TD-EOM-CCSD simulations are numerically stable, the oscillating time-dependent energy expectation value displays an incorrect scaling with the number of atoms resonant with the field. On the other hand, the TDCCSD simulations exhibit the correct scaling in the initial stages of the Rabi cycle but break down when the multiatom system approaches complete population inversion. The numerical stability of the TDCCSD simulations is shown to be very sensitive to any interactions between the atoms and to left amplitudes that do not participate in the expression for time-dependent expectation values for noninteracting atoms. Moreover, we introduce a general theoretical framework for describing the two methods, where the cluster amplitude time derivatives serve as auxiliary conditions related to a shift of the time-dependent Hamiltonian matrix. In this framework, time-dependent coupled-cluster methods exhibit a shifted Hamiltonian matrix with a block upper triangular structure in terms of the number of excited noninteracting subsystems, explaining the correct scaling properties of the TDCCSD method.

arXiv:2301.05546•Physical Review A•Sep 28, 2023•6 citations

Coupled cluster response functions revisited

Thomas Bondo Pedersen, Henrik Koch A

We introduce an inherently real coupled cluster time-dependent expectation value of a Hermitian operator. Based on the expansion of this expectation value in orders of the generally time-dependent perturbation, we subsequently identify the coupled cluster time-independent expectation value, the linear response function, and the quadratic response function. The response functions and their residues behave physically correctly. Spectroscopic observables are identified as residues, whereas the identification of individual transition matrix elements is prohibited. Thus the unphysical behavior of previously published coupled cluster response functions may be viewed not as a consequence of the projection, but rather that identifications are made on the basis of an unphysical expectation value or quasienergy.

The Journal of Chemical Physics•May 15, 1997•148 citations
Equations are derived for the time evolution of time-dependent vibrational coupled cluster (TDVCC) wave functions covering both the TDVCC ket state and the associated so-called Λ bra state. The equations are implemented in the special case of both the Hamiltonian and the cluster operator containing at most two-mode coupling terms. The nontrivial behavior of the evolution of norm, energy, and expectation values due to the nonunitary time-evolution of the nonvariational TDVCC theory is analyzed theoretically and confirmed in numerical experiments that also include time-dependent Hamiltonians. In the spirit of time-independent size-consistency analysis, the separability of both the coupled cluster and Λ states for noninteracting systems is studied. While the coupled cluster state clearly has the correct behavior, the behavior of the Λ state is more intricate, and the consequence for different properties is shown theoretically and numerically. Overall, the numerical experiments show that TDVCC in incomplete expansions gives higher accuracy than a standard linear variational wave function parameterization with the same number of independent parameters, while equivalent results are obtained for complete expansions. The efficiency of the methodology is illustrated in computations on polycyclic aromatic hydrocarbons with up to 156 modes.
We develop time-dependent vibrational coupled cluster with time-dependent modals (TDMVCC), where an active set of one-mode basis functions (modals) is evolved in time alongside coupled-cluster wave-function parameters. A biorthogonal second quantization formulation of many-mode dynamics is introduced, allowing separate biorthogonal bases for the bra and ket states, thus ensuring complex analyticity. We employ the time-dependent bivariational principle to derive equations of motion for both the one-mode basis functions and the parameters describing the cluster (T) and linear de-excitation (L) operators. The choice of constraint (or gauge) operators for the modal time evolution is discussed. In the case of untruncated cluster expansion, the result is independent of this choice, but restricting the excitation space removes this invariance; equations for the variational determination of the constraint operators are derived for the latter case. We show that all single-excitation parts of T and L are redundant and can be left out in the case of variationally determined constraint-operator evolution. Based on a pilot implementation, test computations on Henon-Heiles model systems, the water molecule, and a reduced-dimensionality bi-thiophene model are presented, showing highly encouraging results for TDMVCC. It is demonstrated how TDMVCC in the limit of a complete cluster expansion becomes equivalent to multiconfiguration time-dependent Hartree for the same active-space size. Similarly, it is discussed how TDMVCC generally gives better and more stable results than its time-independent-modals counterpart, while equivalent results are obtained for complete expansions and full one-mode basis sets.
ABSTRACT In recent decades, coupled cluster theory has proven valuable in accurately describing correlation in many‐body systems, particularly in time‐independent computations of molecular electronic structure and vibrations. This review describes recent advancements in using coupled cluster parameterizations for time‐dependent wave functions for the efficient computation of the quantum dynamics associated with the motion of nuclei. It covers time‐dependent vibrational coupled cluster (TDVCC) and time‐dependent modal vibrational coupled cluster (TDMVCC), which employ static and adaptive basis sets, respectively. We discuss the theoretical foundation, including many‐mode second quantization, bivariational principles, and various parameterizations of time‐dependent bases. Additionally, we highlight key features that make TDMVCC promising for future quantum dynamical simulations. These features include fast configuration‐space convergence, the use of a compact adaptive basis set, and the possibility of efficient implementations with a computational cost that scales only polynomially with system size.
We derive equations for describing the time evolution of variational wave functions in linear and exponential parameterization with a second-quantization (SQ) formulation. The SQ formalism covers time-dependent Hartree (TDH), while exact states and approximate vibrational configuration interaction wave functions are described using state-transfer operators. We present detailed expressions for efficient evaluation of TDH in linear (L-TDH) and exponential (X-TDH) parametrization and an efficient implementation supporting linear scaling with respect to the number of degrees of freedom M when the Hamiltonian operator contains a constant number of terms per mode independently of the size of the system. The computational cost of the X-TDH method is reduced significantly compared to the L-TDH method for systems with many operator terms per mode such as is typical for accurate molecular potential-energy surfaces. Numerical results for L-TDH and X-TDH are presented which confirm the theoretical reduction of the M scaling compared to standard first-quantization formulations. Calculations on Henon-Heiles potentials with more than 105 dimensions and polycyclic aromatic hydrocarbons with up to 264 modes have been performed. Thus, the SQ formulation and the X-TDH method pave the way for studying the time-resolved quantum dynamics of large molecules.
A dynamical formulation of coupled cluster theory is derived using a variational principle. By allowing time-dependent single-particle functions, a high degree of adaptivity is introduced, allowing complex systems to be simulated with high accuracy. Equations of motion are derived which are shown to be suitable for computer implementation. The method, called adaptive time-dependent coupled cluster, is a strict generalization of the formulation used in standard coupled cluster response theory, and also represents a systematic hierarchy of size-consistent approximations, including standard time-dependent Hartree--Fock as a trivial case.
The time-dependent equation-of-motion coupled-cluster singles and doubles (TD-EOM-CCSD) and time-dependent coupled-cluster singles and doubles (TDCCSD) methods are compared by simulating Rabi oscillations for distant atoms in a classical electromagnetic field. While the TD-EOM-CCSD simulations are numerically stable, the oscillating time-dependent energy expectation value displays an incorrect scaling with the number of atoms resonant with the field. On the other hand, the TDCCSD simulations exhibit the correct scaling in the initial stages of the Rabi cycle but break down when the multiatom system approaches complete population inversion. The numerical stability of the TDCCSD simulations is shown to be very sensitive to any interactions between the atoms and to left amplitudes that do not participate in the expression for time-dependent expectation values for noninteracting atoms. Moreover, we introduce a general theoretical framework for describing the two methods, where the cluster amplitude time derivatives serve as auxiliary conditions related to a shift of the time-dependent Hamiltonian matrix. In this framework, time-dependent coupled-cluster methods exhibit a shifted Hamiltonian matrix with a block upper triangular structure in terms of the number of excited noninteracting subsystems, explaining the correct scaling properties of the TDCCSD method.
We introduce an inherently real coupled cluster time-dependent expectation value of a Hermitian operator. Based on the expansion of this expectation value in orders of the generally time-dependent perturbation, we subsequently identify the coupled cluster time-independent expectation value, the linear response function, and the quadratic response function. The response functions and their residues behave physically correctly. Spectroscopic observables are identified as residues, whereas the identification of individual transition matrix elements is prohibited. Thus the unphysical behavior of previously published coupled cluster response functions may be viewed not as a consequence of the projection, but rather that identifications are made on the basis of an unphysical expectation value or quasienergy.