Explicitly correlated Fock-space coupled-cluster singles and doubles method for (1,1), (0,2), and (2,0) sectors.

A linearly approximated explicitly correlated coupled-cluster singles and doubles model for the Fock-space coupled-cluster method has been formulated and implemented. An extension of the Fock-space wave operators is introduced in order to treat the short-range correlation effects for excited and doubly electron-attached states. We show that an effective reduction in the number of active virtuals can also be obtained by improving how the short-range correlation is treated. Numerical results to gauge the performance for valence and Rydberg excitation energies, double ionization potentials, and double electron attachment energies of several molecules are obtained. Statistical measures of the errors in excitation energies show that the explicitly correlated results are within 0.1 eV from the complete basis set limit already at the double-ζ level unless the excitation energies are too close to the ionization thresholds. Similar accuracy is seen for the double ionization potentials and double electron attachment energies.

Similarity-transformed equation-of-motion coupled-cluster singles and doubles method with spin-orbit effects for excited states.

The similarity transformed equation-of-motion coupled-cluster method (STEOM-CCSD) for excited states is extended to treat spin-orbit coupling interactions (SOIs) using the spin-orbit mean field approximation of the Breit-Pauli Hamiltonian. Two possible schemes to include the spin-orbit splittings of excited states within the STEOM-CCSD model are formulated. They are identified as "diagonalize-then-perturb" and "perturb-then-diagonalize" approaches. The second approach is more suited for cases where SOI is larger, and the first approach breaks down. With the aid of the standard many-body diagrammatic techniques, expressions for all the necessary matrix elements can be derived. These new formulations are implemented in the ACES III suite of parallel coupled cluster programs, and benchmark studies are performed. Numerical tests for several atoms and molecules show a good agreement of calculated spin-orbit splittings to experiment, while also documenting the numerical differences between the two approaches.

Long-range dispersion C6 coefficient for SF6 dimer: Experimental and theoretical study.

The long-range dispersion C6 coefficient for the SF6 dimer is experimentally measured using a technique that uses the expansion of a supersonic pulse jet into a vacuum. A dynamic model of the jet enables us to correlate the position of the maximal peak in the time-of-flight spectrum with the initial conditions of the experiment and the parameters of the intermolecular interaction potential. Due to the low temperature of the jet target, the C6 coefficient can be extracted directly from the experimental results. Theoretical calculation of the C6 dispersion coefficient is also performed by using linearly approximated explicitly correlated coupled-cluster singles and doubles (CCSD(F12)) method with the subsequent utilization of the Casimir-Polder formula. Good agreement of experimental and theoretical results confirms the reliability of results.

Excitation energies with spin-orbit couplings using equation-of-motion coupled-cluster singles and doubles eigenvectors.

A method of calculation of excited states with spin-orbit couplings, which utilizes left and right eigenvectors of equation-of-motion coupled-cluster singles and doubles model has been formulated and implemented. The spin-orbit interactions are introduced by using the spin-orbit mean field approximation of the Briet-Pauli Hamiltonian. In order to evaluate all the necessary matrix elements, a scheme based on the diagrammatic representation of the second-quantized form of the spin-orbit interaction operator and the standard rules of second-quantized algebra is presented. We posit that this scheme is general and much simpler to use than the often used rules derived for the configuration state functions by using the Wigner-Eckart theorem. We show that the spin-orbit coupled states (i.e., target relativistic states) must satisfy specific conditions in order to classify them according to the double group symmetry. This interrelation between the structure of the target relativistic states and its double group symmetry is discussed in detail. An algorithm to classify the target states according to the irreducible representation of the double group symmetry is offered and implemented. Numerical tests for several atoms and molecules show good agreement of predicted and experimental spin-orbit splittings of the target excited states.

Explicitly correlated coupled-cluster theory for static polarizabilities.

A method of calculation of static polarizabilities with wavefunctions, corresponding to linearly approximated explicitly correlated coupled-cluster singles and doubles [CCSD(F12)] model, has been formulated and implemented. For the proper description of the response of system on applied electric field, modified ansatz is introduced for geminal part of cluster operators. Such extension of CCSD(F12) model provides balanced description of both perturbed and unperturbed wave functions, what leads to the increase of the accuracy of target polarizabilities. As a part of algorithm, explicitly correlated version of coupled-perturbed CCSD equations has also been derived and implemented. Numerical tests conducted for the set of eight molecules show good agreement between static polarizabilities, calculated with developed explicitly correlated approach and corresponding complete basis set results in regular CCSD already at triple-ζ level.

Electric multipole moments calculation with explicitly correlated coupled-cluster wavefunctions.

A method of calculation of expectation values of dipole and quadrupole moments with wavefunctions, corresponding to linearly approximated explicitly correlated coupled-cluster singles and doubles [CCSD(F12)] model has been formulated and implemented. As a part of algorithm, explicitly correlated version of Λ equations has also been derived and implemented. Numerical tests, conducted for sets of molecules, show that explicitly correlated results for expectation values of dipole moment are accurate up to 0.01 a.u. already at a double-ζ level compared to those in the complete basis set limit. The corresponding results for quadrupole moments at double-ζ level are accurate up to 0.1 a.u., while for the triple-ζ bases errors do not exceed 0.01 a.u.

Explicitly correlated similarity transformed equation-of-motion coupled-cluster method.

Similarity transformed equation-of-motion method, based on linearly approximated explicitly correlated coupled-cluster singles and doubles [CCSD(F12)] model, has been formulated and implemented. An extension of similarity transformation operator is introduced in order to treat short-range correlation effects for excited states. Additionally, effective reduction of the number of active virtuals can be obtained by such modification. Numerical tests for sets of valence and Rydberg excited states of several molecules are conducted. Statistical measures of errors in excitation energies show that explicitly correlated results are accurate up to 0.1 e.V already at a double-ζ level compared to those in the complete basis set limit, if the excitation energy is not too close to an ionization threshold. An example of long-range charge transfer excitation is also considered and highly accurate results are obtained.

Explicitly correlated second-order Møller-Plesset perturbation theory employing pseudospectral numerical quadratures.

We implemented explicitly correlated second-order Møller-Plesset perturbation theory with numerical quadratures using pseudospectral construction of grids. Introduction of pseudospectral approach for the calculation of many-electron integrals gives a possibility to use coarse grids without significant loss of precision in correlation energies, while the number of points in the grid is reduced about nine times. The use of complementary auxiliary basis sets as the sets of dealiasing functions is justified at both theoretical and computational levels. Benchmark calculations for a set of 16 molecules have shown the possibility to keep an error of second-order correlation energies within 1 milihartree (mH) with respect to MP2-F12 method with dense grids. Numerical tests for a set of 13 isogyric reactions are also performed.