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.

4 GENERAL PROGRAM STRUCTURE 15<br />

4.2 Files<br />

<strong>MOLPRO</strong> uses three sequential text files, namely the input file, the output file, and the punch<br />

file. The punch file is a short form of the output which contains the most important data and<br />

results, such as geometries, energies, dipole moments, etc. The punch file can be processed<br />

by the separate program READPUN, which selects specific results by keywords and is able to<br />

produce ordered tables in user supplied format. Furthermore, there are up to 9 binary <strong>MOLPRO</strong><br />

files available, each one known to the program simply by its number (1 to 9). By default, they<br />

are temporary files, usually allocated dynamically by the program, but they can be connected<br />

to permanent files with the FILE command. Each file is direct access, and word addressable<br />

(word=64 bit usually), but is organised in records of any length. The name, address and length<br />

of each record is held in a directory at the start of the file.<br />

File 1 is the main file, holding basis set, geometry, and the one and two electron integrals. By<br />

default, file 2 is the dump file and used to store the wavefunction information, i.e. orbitals, CI<br />

coefficients, and density matrices. File 3 is an auxiliary file which can be used in addition to file<br />

2 for restart purposes. Often files 1 and 2 (and 3) are declared as permanent files (see FILE) to<br />

enable restarts. Storing the wavefunction information on file 2 is useful, since the integral file is<br />

overwritten at each new geometry, while the orbitals and CI coefficients of one calculation can<br />

be used as a starting guess for the next calculation at a neighbouring geometry. Files 4 to 8 are<br />

used as scratch space, e.g., for sorting the integrals, storage of transformed integrals and of the<br />

CI vectors. These files should normally not be made permanent.<br />

Note that the file name appearing in molpro input is always converted to lower case on unix<br />

machines.<br />

4.3 Records<br />

Record names are positive integers, and are usually referred to in the format record.file, e.g.,<br />

2100.2 means the record called 2100 on file 2. Note that these names are quite arbitrary, and<br />

their numerical values have nothing to do with the order of the records in the file. Record names<br />

≤ 2000 are reserved for standard quantities (e.g. integrals, properties etc.) and you should never<br />

use these in an input, unless you know exactly what you are doing. Some important default<br />

records to remember are<br />

2100 RHF dump record (closed and open-shell)<br />

2200 UHF dump record<br />

2140 MCSCF dump record<br />

4100 CPHF restart information<br />

5000 MCSCF gradient information<br />

5100 CP-MCSCF gradient information<br />

5200 MP2 gradient information<br />

5300 Hessian restart information<br />

5400 Frequencies restart information<br />

6300 Domain restart information<br />

If an input contains several wavefunction calculations of the same type, e.g., several MCSCF<br />

calculations with different active spaces, the record number will be increased by 1 for each<br />

calculation of the same type. Thus, the results of the first SCF calculation in an input are stored

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

Saved successfully!

Ooh no, something went wrong!