19.01.2015 Views

MOLPRO

MOLPRO

MOLPRO

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

31 EXPLICITLY CORRELATED METHODS 231<br />

31.2.3 The fixed amplitude ansatz (FIX)<br />

The diagonal ansatz is used and the amplitudes of the explicitly correlated configurations are<br />

determined by the wavefunction cusp conditions, i.e.<br />

T i j,1<br />

i j = 1 2<br />

T i j,−1<br />

i j = 1 4<br />

This ansatz is orbital invariant, size consistent and free of GBSSE.<br />

31.3 RI Approximations<br />

Various approximations such as A, B, C, HY1, HY2 exist for the matrix elements of the firstorder<br />

Hamiltonian (see I). They differ in the way the RI approximations are made. In the limit<br />

of a complete RI basis, approximations B and C are identical and yield the exact result for a<br />

given wave function ansatz. We generally recommend approximation C, which is simpler and<br />

more efficient than approximation B. Normally, the union of the AO and RI basis sets is used<br />

to approximate the resolution of the identity (CABS approach). In the hybrid approximations<br />

(HY1, HY2, HX) only the AO basis is used in some less important terms. Together with the<br />

recommended approximation C, HY1 or HY2 can be used; HY2 is more accurate, HY1 more<br />

efficient. In most cases, approximation 3C(HY1) provides an excellent compromise between<br />

accuracy and efficiency. In approximation A, all terms involving exchange operators are neglected.<br />

This approximation is used along with local approximations in our low-order scaling<br />

LMP2-F12/3*A(loc) method that can be applied to large molecules (cf. section 31.11)<br />

If the extended Brillouin condition (EBC, see I) is assumed, the explicitly correlated and conventional<br />

amplitude equations decouple and can be solved independently. These approximations<br />

are denoted by a star, e.g. 3*C.<br />

31.4 Basis sets<br />

In <strong>MOLPRO</strong> the F12 integrals can only be computed using density fitting (DF) approximations.<br />

The many electron integrals are approximated by resolutions of the identity (RI) expansions.<br />

Thus, F12 calculations require three different basis sets: the orbital (AO) basis, the DF basis,<br />

and the RI basis.<br />

We recommend as AO basis sets the augmented correlation consistent basis sets (denoted AVnZ)<br />

or the specially optimized correlation consistent F12 basis sets (denoted VnZ-F12, cf. K.A.<br />

Peterson and H.-J. Werner, J. Chem. Phys. 128, 084102 (2008)). Normally, triples zeta basis<br />

sets (AVTZ or VTZ-F12) yield excellent results that are close to the basis set limit. Diffuse<br />

basis functions are rather essential both for the HF and MP2-F12 energies, and therefore the<br />

standard VTZ sets are not recommended. If the AVnZ or VnZ-F12 orbital basis sets are used,<br />

suitable density fitting (DF) basis and resolution of the identity (RI) basis sets are automatically<br />

chosen. For the AVnZ orbital basis sets, AVnZ/MP2FIT amd VnZ/JKFIT basis sets are used by<br />

default for the DF and RI, respectively. The associated optimized CABS basis set of Peterson<br />

et al. can be chosen by specifying RI BASIS=OPTRI. For the VnZ-F12 orbital basis, the<br />

associated CABS (OPTRI) basis sets are used by default. Other basis sets can be chosen using<br />

the DF BASIS, DF BASIS EXCH and RI BASIS options (cf. section 31.6). See section 15 for<br />

more details about density fitting.<br />

This is an example for using multiple basis sets for density fitting and resolution of the identity

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!