Spin unrestricted (RHF-UCCSD) and partially spin restricted (RHF-RCCSD) open-shell coupled cluster theories as described in J. Chem. Phys. 99 (1993) 5219 (see also erratum, J. Chem. Phys., 112 (2000) 3106) are available in MOLPRO. In both cases a high-spin RHF reference wavefunction is used. No coupled cluster methods based on UHF orbitals are implemented in MOLPRO (the only correlation method in MOLPRO which uses UHF orbitals is UMP2). In the description that follows, the acronyms RCCSD and UCCSD are used, but the theories should normally be referred to as RHF-RCCSD, RHF-UCCSD, in order to distinguish them from alternative ansätze based on spin-unrestricted orbitals. The program will accept either the full or abbreviated acronyms as input commands.
In the RCCSD theory certain restrictions among the amplitudes are introduced, such that the linear part of the wavefunction becomes a spin eigenfunction (this is not the case in the UCCSD method, even if an RHF reference function is used). At present, the implementation of RCCSD is only preliminary, and no CPU time is saved by as compared to UCCSD. However, improved algorithms, as described in the above publication, are currently being implemented, and will be available in the near future.
The input and the options are exactly the same as for closed-shell CCSD, except that RCCSD or UCCSD are used as keywords. By default, the open-shell orbitals are the same as used in the RHF reference function, but this can be modified using OCC, CLOSED, and WF cards.
Distinguishable cluster calculations can be performed using RDCSD or UDCSD commands.
Perturbative triples corrections are computed as follows:
In fact, all three contributions are always computed and printed. The following variables are used to store the results (here CCSD stands for either UCCSD or RCCSD):
It should be noted that in open-shell cases the triples energy slightly depends on the treatment of core orbitals. In MOLPRO pseudo-canonical alpha and beta spin orbitals (Chem. Phys. Letters 186 (1991) 130) are generated by block-diagonalizing the corresponding Fock matrices in the space of valence orbitals, leaving frozen core orbitals untouched. Some other programs include the frozen core orbitals in the canonicalization and transformation. Because of core-valence mixing this leads to slightly different energies. Neither of the two methods can be regarded as better or more justified -- it is just a matter of definition. However, the method in MOLPRO is more efficient since the subsequent integral transformation involves only valence orbitals and no core orbitals.