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.

47 POTENTIAL ENERGY SURFACES (SURF) 356<br />

1D: CCSD(T)/cc-pVTZ<br />

2D: MP4(SDQ)/cc-pVTZ<br />

3D: MP4(SDQ)/cc-pVTZ<br />

{SURF,Start1D=label1<br />

VMULT,Start2D=label2,Start3D=label2,Multi=3}<br />

MULTI=4: The 1D calculation provides all information for the 2D<br />

potential but does not so for the 3D part. Hence, the 1D contribution<br />

and the 2D contributions need to be computed twice. Examples:<br />

1D: CCSD(T)/cc-pVTZ<br />

2D: CCSD(T)/cc-pVTZ<br />

3D: MP4(SDQ)/cc-pVTZ<br />

{SURF,Start1D=label1<br />

VMULT,Start2D=label1,Start3D=label2,Multi=4}<br />

1D: CCSD(T)/cc-pVTZ<br />

2D: MP2/cc-pVTZ<br />

3D: MP2/cc-pVDZ<br />

{SURF,Start1D=label1<br />

VMULT,Start2D=label2,Start3D=label3<br />

VMULT,Var2D=EMP2,Multi=4}<br />

START2D=label<br />

VAR2D=variable<br />

In 2D and 4D calculations (i.e. NDIM=2,4) the VMULT command<br />

can be used as well. In 4D calculations the last level must always<br />

be identical to the 3D level. In 2D the meaning of MULTI=1 and<br />

MULTI=3 is the same. Likewise, MULTI=2 and MULTI=4 are the<br />

same in case of 2D calculations.<br />

START2D and START3D define labels in the input stream in order to<br />

compute the 2D and 3D terms at different levels of electronic structure<br />

theory than the 1D terms. The use of the START2D and START3D<br />

commands usually requests the use of GOTO commands in the input.<br />

The keywords VAR2D and VAR3D are defined in full analogy to the<br />

VAR1D option. They specify the internal variable (e.g. ENERGY,<br />

EMP2, CCSD, ...) to be read out for a given grid point.<br />

The following example shows a 1D:CCSD(T)/cc-pVTZ; 2D:MP4(SDQ)/cc-pVTZ and 3D:MP2/ccpVTZ<br />

multi-level calculation. As the MP2 energy is a byproduct of the CCSD(T) and MP4(SDQ)<br />

calculations only the 1D grid points will be computed twice (at the CCSD(T) and MP4(SDQ)<br />

levels). The 1D and 2D energies will be obtained from the internal variable ENERGY while the<br />

3D energies make use of the EMP2 variable.<br />

memory,50,m<br />

geomtyp=xyz<br />

orient,mass<br />

geometry={<br />

6<br />

Ethen<br />

C 0.0000000000 0.0000000000 -0.6685890718<br />

C 0.0000000000 0.0000000000 0.6685890718<br />

H 0.0000000000 -0.9240027061 -1.2338497710<br />

H 0.0000000000 0.9240027061 -1.2338497710

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

Saved successfully!

Ooh no, something went wrong!