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.

19 THE MCSCF PROGRAM MULTI 139<br />

of the C 2v calculation as reference for all other calculations. In this case one would have to take<br />

out the second-last input card, which sets reforb=2141.2.<br />

***,H2S diabatic A" states<br />

basis=VDZ<br />

symmetry,x,planeyz<br />

orient,noorient<br />

geometry={s;h1,s,r1;h2,s,r2,h1,theta}<br />

gprint,orbitals,civector<br />

text,reference calculation for C2V<br />

theta=92.12,r1=2.3,r2=2.3<br />

{hf;occ,7,2;wf,18,1}<br />

!use cc-pVDZ basis set<br />

!use Cs symmetry & fix orientation of the molecule<br />

!dont allow automatic reorientation<br />

!Z-matrix geometry input<br />

!global print options<br />

!reference geometry<br />

!scf calculation for ground state<br />

{multi;occ,9,2;closed,4,1;<br />

!define active and inactive spaces<br />

wf,18,2;state,2;<br />

!two A" states (1B1 and 1A2 in C2v)<br />

orbital,2140.2} !save orbitals to 2140.2<br />

reforb=2140.2<br />

text,calculations at displaced geometries<br />

rd=[2.4,2.5,2.6]<br />

do i=1,#rd<br />

r2=rd(i)<br />

{multi;occ,9,2;closed,4,1;<br />

wf,18,2;state,2;<br />

orbital,2141.2<br />

diab,reforb}<br />

reforb=2141.2<br />

enddo<br />

!define a range of bond distances<br />

!loop over displaced geometries<br />

!set r2 to current distance<br />

!same wavefunction definition as at reference geom.<br />

!save new orbitals to record<br />

!compute diabatic orbitals using reference orbitals<br />

!stored on record reforb<br />

!set variable reforb to the new orbitals.<br />

http://www.molpro.net/info/current/examples/h2s_diab.com<br />

See section 38 for the automatic generation of diabatic energies.<br />

19.6 Selecting the optimization methods<br />

By default, MULTI uses the non-linear optimization method developed by Werner, Meyer, and<br />

Knowles. Other methods, such as the Newton-Raphson procedure or the Augmented Hessian<br />

procedure, are also implemented and can be selected using the ITERATIONS directive (for<br />

state-averaged calculations, only the non-linear optimization method can be used). For CASSCF<br />

calculations, the CI problem is solved in a basis of Slater determinants, unless a CONFIG card<br />

is given. Some procedures may be disabled using the DONT directive.<br />

19.6.1 Selecting the CI method<br />

CONFIG,key;<br />

key may be DET or CSF, and defaults to CSF. If no CONFIG or SELECT card is given, the<br />

default is determinants (CASSCF).

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

Saved successfully!

Ooh no, something went wrong!