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.
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 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.
Multilevel Hartree-Fock can be used for single-point calculations, and as a reference wavefunction for reduced space coupled cluster calculations.
Quantum electrodynamics Hartree-Fock supports the same features as restricted closed-shell Hartree-Fock (RHF) except geometry optimization.
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
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.
Ground state energies at the CCSD(T) level of theory are implemented.
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.
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.
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.
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.
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.
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.
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.
Ground state energies at the MP2 level of theory are implemented.
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 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.
Analytical gradients are available on the HF level of theory. Numerical gradients are available either using forward, or central difference.