Features

The eT program is first and foremost a coupled cluster program, with the CCS, CC2, CCSD, CCSD(T), and CC3 methods implemented. In addition to the standard coupled cluster methods, eT has multilevel CC2, multilevel CCSD, and coupled cluster time propagation. Besides coupled cluster methods, MP2, FCI and TDHF are implemented. Hartree-Fock calculations, both restricted closed shell, restricted open-shell, and unrestricted, are available. Furthermore, the multilevel Hartree-Fock method is implemented. QM/MM is available at both HF and CC level.

A detailed description of the features available for each method is given below.

Hartree-Fock

The available Hartree-Fock methods are

  • Restricted Hartree-Fock (RHF, ROHF)

  • Unrestricted Hartree-Fock (UHF, CUHF)

  • Multilevel Hartree-Fock (MLHF)

  • Quantum electrodynamics Hartree-Fock (QED-HF)

Restricted and unrestricted Hartree-Fock

Restricted closed-shell Hartree-Fock (RHF) can be used for single-point calculations, geometry optimizations, properties (dipole and quadrupole), and as a reference wavefunction for coupled cluster calculations.

Restricted open-shell Hartree-Fock (ROHF) can be used for single-point calculations.

Unrestricted Hartree-Fock (UHF) can be used for single-point calculations.

Setting up a Hartree-Fock calculation

Multilevel Hartree-Fock

Multilevel Hartree-Fock can be used for single-point calculations, and as a reference wavefunction for reduced space coupled cluster calculations.

Setting up a multilevel Hartree-Fock calculation

Quantum electrodynamics Hartree-Fock

Quantum electrodynamics Hartree-Fock supports the same features as restricted closed-shell Hartree-Fock (RHF) except geometry optimization.

Coupled cluster

The implemented coupled cluster methods in eT are

  • Standard methods: CCS and CCSD

  • Perturbative methods: CC2, CC3 and CCSD(T)

  • Multilevel CC2 and CCSD

Additionally, the code can perform time-propagation of some of the implemented coupled cluster wavefunctions.

The features implemented for the coupled cluster methods vary somewhat, and are detailed below.

Setting up a coupled cluster calculation

CCS, CC2, CCSD and CC3

For CCS, CC2, CCSD, and CC3, eT offers ground and excited state calculations in addition to dipole and quadrupole moments, and EOM oscillator strengths.

For CCS, CC2, and CCSD also EOM polarizabilities and linear response oscillator strengths as well as polarizabilities are implemented.

For CCS, CC2, and CCSD triplet excitation energies are implemented.

CCSD(T)

Ground state energies at the CCSD(T) level of theory are implemented.

Low-memory CC2

In eT, there are two versions of the CC2 code. The standard CC2 code has a memory requirement proportional to \(n_o^2 n_v^2\) , where \(n_o\) is the number of occupied orbitals and \(n_v\) is the number of virtual orbitals. However, the low-memory CC2 implementation has an \(N^2\) memory requirement, where \(N\) is the number of orbitals.

Currently, only ground and excited state energies are available with the low-memory CC2 code.

Multilevel coupled cluster

The multilevel CC2 (MLCC2) and multilevel CCSD (MLCCSD) methods are available with correlated natural transition orbitals, Cholesky orbitals (occupied and virtual), and projected atomic orbitals.

Currently, only ground and excited states are available at the MLCCSD level of theory.

Setting up a multilevel coupled cluster calculation

Time-dependent coupled cluster

Real-time propagation is available for CCS and CCSD methods. It solves the differential equations describing the time evolution of cluster amplitudes and Lagrange multipliers. The available integrators are

  • Euler

  • Gauss-Legendre (with order 2, 4 and 6)

  • Runge-Kutta (with order 4)

The time dependent quantities that are available as output are:

  • Amplitudes

  • Multipliers

  • Density matrix

  • Energy

  • Electric field

  • Dipole moment

In addition, visualizable time-dependent density and spectra given by the Fast Fourier Transform of the dipole moment and of the electric field can be requested as output.

Environments

The available QM/MM methods are

  • Electrostatic embedding (non-polarizable)

  • Polarizable QM/Fluctuating charges

  • Polarizable Continuum Model

Both the methods can be coupled both with HF or CC wavefunctions. In case of QM/FQ, only HF ground state Fock and MOs are affected by QM/FQ.

Cholesky decomposition of the ERIs

The coupled cluster code is based on Cholesky decomposed electron repulsion integrals (ERIs). For sufficiently small systems, T1-transformed ERIs are constructed from the Cholesky vectors and stored in memory. For systems where the ERIs cannot be placed in memory, the T1-transformed Cholesky vectors are stored in memory if possible and ERIs are constructed from these vectors on the fly. For larger systems, the Cholesky vectors are stored on disk.

Visualization

Ground state densities can be written to .plt or .cube files that are readable by Chimera for both Hartree-Fock and coupled cluster. For coupled cluster, it is also possible to plot transition densities.

FCI

Full CI calculations are available for both closed (RHF reference) and open (ROHF reference) shells. In addition, dipole and quadrupole moments can be calculated for the ground state.

MP2

Ground state energies at the MP2 level of theory are implemented.

Time-dependent Hartree-Fock

Excitation energies and oscillator strengths at the TDHF level of theory are available either within the random phase (RPA) or the Tamm-Dancoff approximation (TDA). In addition, static and frequency-dependent polarizabilities can be obtained.

Quantum electrodynamics time-dependent Hartree-Fock

Quantum electrodynamics time-dependent Hartree-Fock supports the same features as time-dependent Hartree-Fock except the polarizabilities. In addition, a measure of the photon character of each excitation is computed.

Geometry optimization

Ground state

Analytical gradients are available on the HF level of theory. Numerical gradients are available either using forward, or central difference.