manual   quickstart   instguide   update   basis
  Index  

Contents

• Introduction to MOLPRO
• MOLPRO on the WWW
• References
• 1 HOW TO READ THIS MANUAL
• 2 RUNNING MOLPRO
  • 2.1 Options
  • 2.2 Running MOLPRO on parallel computers
  
  
• 3 DEFINITION OF MOLPRO INPUT LANGUAGE
  • 3.1 Input format
  • 3.2 Commands
  • 3.3 Directives
  • 3.4 Global directives
  • 3.5 Options
  • 3.6 Data
  • 3.7 Expressions
  • 3.8 Intrinsic functions
  • 3.9 Variables
  • 3.10 Procedures
  
  
• 4 GENERAL PROGRAM STRUCTURE
  • 4.1 Input structure
  • 4.2 Files
  • 4.3 Records
  • 4.4 Restart
  • 4.5 Data set manipulation
  • 4.6 Memory allocation
  • 4.7 Multiple passes through the input
  • 4.8 Symmetry
  • 4.9 Defining the wavefunction
  • 4.10 Defining orbital subspaces
  • 4.11 Selecting orbitals and density matrices (ORBITAL, DENSITY)
  • 4.12 Summary of keywords known to the controlling program
  • 4.13 MOLPRO help
  
  
• 5 INTRODUCTORY EXAMPLES
  • 5.1 Using the molpro command
  • 5.2 Simple SCF calculations
  • 5.3 Geometry optimizations
  • 5.4 CCSD(T)
  • 5.5 CASSCF and MRCI
  • 5.6 Tables
  • 5.7 Procedures
  • 5.8 Do loops
  
  
• 6 PROGRAM CONTROL
  • 6.1 Starting a job (***)
  • 6.2 Ending a job (--)
  • 6.3 Restarting a job (RESTART)
  • 6.4 Including secondary input files (INCLUDE)
  • 6.5 Allocating dynamic memory (MEMORY)
  • 6.6 DO loops (DO/ENDDO)
  • 6.7 Branching (IF/ELSEIF/ENDIF)
  • 6.8 Procedures (PROC/ENDPROC)
  • 6.9 Text cards (TEXT)
  • 6.10 Checking the program status (STATUS)
  • 6.11 Global Thresholds (GTHRESH)
  • 6.12 Global Print Options (GPRINT/NOGPRINT)
  • 6.13 One-electron operators and expectation values (GEXPEC)
  
  
• 7 FILE HANDLING
  • 7.1 FILE
  • 7.2 DELETE
  • 7.3 ERASE
  • 7.4 DATA
  • 7.5 Assigning punch files (PUNCH)
  • 7.6 MOLPRO system parameters (GPARAM)
  
  
• 8 VARIABLES
  • 8.1 Setting variables
  • 8.2 Indexed variables
  • 8.3 String variables
  • 8.4 System variables
  • 8.5 Macro definitions using string variables
  • 8.6 Indexed Variables (Vectors)
  • 8.7 Vector operations
  • 8.8 Special variables
  • 8.9 Displaying variables
  • 8.10 Clearing variables
  • 8.11 Reading variables from an external file
  
  
• 9 TABLES AND PLOTTING
  • 9.1 Tables
  • 9.2 Plotting
  • 9.3 Diatomic potential curve analysis
  
  
• 10 MOLECULAR GEOMETRY
  • 10.1 Geometry specifications
  • 10.2 Symmetry specification
  • 10.3 Writing Gaussian, XMol or MOLDEN input (PUT)
  • 10.4 Geometry Files
  • 10.5 Lattice of point charges
  • 10.6 Redefining and printing atomic masses
  • 10.7 Dummy centres
  
  
• 11 BASIS INPUT
  • 11.1 Overview: sets and the basis library
  • 11.2 Cartesian and spherical harmonic basis functions
  • 11.3 The basis set library
  • 11.4 Default basis sets
  • 11.5 Default basis sets for individual atoms
  • 11.6 Basis blocks
  • 11.7 Auxiliary basis sets for density fitting or resolution of the identity
  • 11.8 Primitive set definition
  • 11.9 Contracted set definitions
  • 11.10 Examples
  
  
• 12 EFFECTIVE CORE POTENTIALS
  • 12.1 Input from ECP library
  • 12.2 Explicit input for ECPs
  • 12.3 Example for explicit ECP input
  • 12.4 Example for ECP input from library
  
  
• 13 CORE POLARIZATION POTENTIALS
  • 13.1 Input options
  • 13.2 Example for ECP/CPP
  
  
• 14 INTEGRATION
  • 14.1 Sorted integrals
  • 14.2 INTEGRAL-DIRECT CALCULATIONS (GDIRECT)
  
  
• 15 DENSITY FITTING
  • 15.1 Options for density fitting
  
  
• 16 THE SCF PROGRAM
  • 16.1 Options
  • 16.2 Defining the wavefunction
  • 16.3 Saving the final orbitals
  • 16.4 Starting orbitals
  • 16.5 Rotating pairs of orbitals
  • 16.6 Using additional point-group symmetry
  • 16.7 Expectation values
  • 16.8 Polarizabilities
  • 16.9 Miscellaneous directives
  
  
• 17 THE DENSITY FUNCTIONAL PROGRAM
  • 17.1 Options
  • 17.2 Directives
  • 17.3 Numerical integration grid control (GRID)
  • 17.4 Density Functionals
  • 17.5 Empirical damped dispersion correction
  • 17.6 Examples
  
  
• 18 ORBITAL LOCALIZATION
  • 18.1 Defining the input orbitals (ORBITAL)
  • 18.2 Saving the localized orbitals (SAVE)
  • 18.3 Choosing the localization method (METHOD)
  • 18.4 Delocalization of orbitals (DELOCAL)
  • 18.5 Localizing AOs(LOCAO)
  • 18.6 Selecting the orbital space
  • 18.7 Ordering of localized orbitals
  • 18.8 Localization thresholds (THRESH)
  • 18.9 Options for PM localization (PIPEK)
  • 18.10 Printing options (PRINT)
  
  
• 19 THE MCSCF PROGRAM MULTI
  • 19.1 Structure of the input
  • 19.2 Defining the orbital subspaces
  • 19.3 Defining the optimized states
  • 19.4 Defining the configuration space
  • 19.5 Restoring and saving the orbitals and CI vectors
  • 19.6 Selecting the optimization methods
  • 19.7 Calculating expectation values
  • 19.8 Miscellaneous options
  • 19.9 Coupled-perturbed MCSCF
  • 19.10 Optimizing valence bond wavefunctions
  • 19.11 Hints and strategies
  • 19.12 Examples
  
  
• 20 THE CI PROGRAM
  • 20.1 Introduction
  • 20.2 Specifying the wavefunction
  • 20.3 Options
  • 20.4 Miscellaneous thresholds
  • 20.5 Print options
  • 20.6 Examples
  • 20.7 Cluster corrections for multi-state MRCI
  
  
• 21 MULTIREFERENCE RAYLEIGH SCHRÖDINGER PERTURBATION THEORY
  • 21.1 Introduction
  • 21.2 Excited state calculations
  • 21.3 Multi-State CASPT2
  • 21.4 Modified Fock-operators in the zeroth-order Hamiltonian.
  • 21.5 Level shifts
  • 21.6 Integral direct calculations
  • 21.7 CASPT2 gradients
  • 21.8 Coupling MRCI and MRPT2: The CIPT2 method
  • 21.9 Further options for CASPT2 and CASPT3
  
  
• 22 NEVPT2 calculations
  • 22.1 General considerations
  • 22.2 Input description
  
  
• 23 MØLLER PLESSET PERTURBATION THEORY
  • 23.1 Expectation values for MP2
  • 23.2 Polarizabilities and second-order properties for MP2
  • 23.3 CPHF for gradients, expectation values and polarizabilities
  • 23.4 Density-fitting MP2 (DF-MP2, RI-MP2)
  • 23.5 Spin-component scaled MP2 (SCS-MP2)
  
  
• 24 THE CLOSED SHELL CCSD PROGRAM
  • 24.1 Coupled-cluster, CCSD
  • 24.2 Quadratic configuration interaction, QCI
  • 24.3 Brueckner coupled-cluster calculations, BCCD
  • 24.4 Singles-doubles configuration interaction, CISD
  • 24.5 The DIIS directive
  • 24.6 Examples
  • 24.7 Saving the density matrix
  • 24.8 Natural orbitals
  • 24.9 Dual basis set calculations
  
  
• 25 EXCITED STATES WITH EQUATION-OF-MOTION CCSD (EOM-CCSD)
  • 25.1 Options for EOM
  • 25.2 Options for EOMPAR card
  • 25.3 Options for EOMPRINT card
  • 25.4 Examples
  • 25.5 Excited states with CIS
  • 25.6 First- and second-order properties for CCSD
  
  
• 26 OPEN-SHELL COUPLED CLUSTER THEORIES
• 27 The MRCC program of M. Kallay (MRCC)
  • 27.1 Installing MRCC
  • 27.2 Running MRCC
  
  
• 28 SMILES
  • 28.1 INTERNAL BASIS SETS
  • 28.2 EXTERNAL BASIS SETS
  • 28.3 Example
  
  
• 29 LOCAL CORRELATION TREATMENTS
  • 29.1 Introduction
  • 29.2 Getting started
  • 29.3 Summary of options
  • 29.4 Summary of directives
  • 29.5 General Options
  • 29.6 Options for selection of domains
  • 29.7 Options for selection of pair classes
  • 29.8 Directives
  • 29.9 Doing it right
  • 29.10 Density-fitted LMP2 (DF-LMP2) and coupled cluster (DF-LCCSD(T0))
  
  
• 30 LOCAL METHODS FOR EXCITED STATES
  • 30.1 Local CC2 and ADC(2)
  • 30.2 Options for EOM
  • 30.3 Parameters on LAPLACE card
  • 30.4 Print options
  • 30.5 Examples
  
  
• 31 EXPLICITLY CORRELATED METHODS
  • 31.1 Reference functions
  • 31.2 Wave function Ansätze
  • 31.3 RI Approximations
  • 31.4 Basis sets
  • 31.5 Symmetry
  • 31.6 Options
  • 31.7 Choosing the ansatz and the level of approximation
  • 31.8 CABS Singles correction
  • 31.9 Pair specific geminal exponents
  • 31.10 CCSD(T)-F12
  • 31.11 DF-LMP2-F12 calculations with local approximations
  • 31.12 Variables set by the F12 programs
  
  
• 32 THE FULL CI PROGRAM
  • 32.1 Defining the orbitals
  • 32.2 Occupied orbitals
  • 32.3 Frozen-core orbitals
  • 32.4 Defining the state symmetry
  • 32.5 Printing options
  • 32.6 Interface to other programs
  • 32.7 Example
  
  
• 33 SYMMETRY-ADAPTED INTERMOLECULAR PERTURBATION THEORY
  • 33.1 Introduction
  • 33.2 First example
  • 33.3 DFT-SAPT
  • 33.4 High order terms
  • 33.5 Density fitting
  • 33.6 SAPT with ECP's
  • 33.7 Examples
  • 33.8 Options
  
  
• 34 PROPERTIES AND EXPECTATION VALUES
  • 34.1 The property program
  • 34.2 Distributed multipole analysis
  • 34.3 Mulliken population analysis
  • 34.4 Natural Bond Orbital Analysis
  • 34.5 Finite field calculations
  • 34.6 Relativistic corrections
  • 34.7 CUBE -- dump density or orbital values
  • 34.8 GOPENMOL -- calculate grids for visualization in gOpenMol
  
  
• 35 RELATIVISTIC CORRECTIONS
  • 35.1 Using the Douglas-Kroll-Hess Hamiltonian
  • 35.2 Example for computing relativistic corrections
  
  
• 36 DIABATIC ORBITALS
• 37 NON ADIABATIC COUPLING MATRIX ELEMENTS
  • 37.1 The DDR procedure
  
  
• 38 QUASI-DIABATIZATION
• 39 THE VB PROGRAM CASVB
  • 39.1 Structure of the input
  • 39.2 Defining the CASSCF wavefunction
  • 39.3 Other wavefunction directives
  • 39.4 Defining the valence bond wavefunction
  • 39.5 Recovering CASSCF CI vector and VB wavefunction
  • 39.6 Saving the VB wavefunction
  • 39.7 Specifying a guess
  • 39.8 Permuting orbitals
  • 39.9 Optimization control
  • 39.10 Point group symmetry and constraints
  • 39.11 Wavefunction analysis
  • 39.12 Controlling the amount of output
  • 39.13 Further facilities
  • 39.14 Service mode
  • 39.15 Examples
  
  
• 40 SPIN-ORBIT-COUPLING
  • 40.1 Introduction
  • 40.2 Calculation of SO integrals
  • 40.3 Calculation of individual SO matrix elements
  • 40.4 Approximations used in calculating spin-orbit integrals and matrix elements
  • 40.5 Calculation and diagonalization of the entire SO-matrix
  • 40.6 Modifying the unperturbed energies
  • 40.7 Examples
  
  
• 41 ENERGY GRADIENTS
  • 41.1 Analytical energy gradients
  • 41.2 Numerical gradients
  • 41.3 Saving the gradient in a variables
  
  
• 42 GEOMETRY OPTIMIZATION (OPTG)
  • 42.1 Options
  • 42.2 Directives for OPTG
  • 42.3 Using the SLAPAF program for geometry optimization
  • 42.4 Examples
  
  
• 43 VIBRATIONAL FREQUENCIES (FREQUENCIES)
  • 43.1 Options
  • 43.2 Printing options (PRINT)
  • 43.3 Saving the hessian and other information (SAVE)
  • 43.4 Restarting a hessian/Frequency calculation (START)
  • 43.5 Coordinates for numerical hessian calculations (COORD)
  • 43.6 Stepsizes for numerical hessian calculations (STEP)
  • 43.7 Numerical hessian using energy variables (VARIABLE)
  • 43.8 Thermodynamical properties (THERMO)
  • 43.9 Examples
  
  
• 44 MAGNETIC PROPERTIES OF MOLECULES
• 45 MINIMIZATION OF FUNCTIONS
  • 45.1 Examples
  
  
• 46 BASIS SET EXTRAPOLATION
  • 46.1 Options
  • 46.2 Extrapolation functionals
  • 46.3 Geometry optimization using extrapolated energies
  • 46.4 Harmonic vibrational frequencies using extrapolated energies
  
  
• 47 POTENTIAL ENERGY SURFACES (SURF)
  • 47.1 Options
  • 47.2 Multi-level calculations
  • 47.3 Restart capabilities
  • 47.4 Linear combinations of normal coordinates
  • 47.5 Scaling of individual coordinates
  • 47.6 Modeling of high-order $n$-body terms
  • 47.7 Quality Check
  • 47.8 Recommendations
  • 47.9 Standard Problems
  
  
• 48 POLYNOMIAL REPRESENTATIONS (POLY)
  • 48.1 Options
  
  
• 49 THE VSCF PROGRAM (VSCF)
  • 49.1 Options
  
  
• 50 THE VCI PROGRAM (VCI)
  • 50.1 Options
  • 50.2 Recommendations
  • 50.3 Examples
  
  
• 51 VIBRATIONAL MP2 (VMP2)
• 52 THE COSMO MODEL
  • 52.1 BASIC THEORY
  
  
• 53 QM/MM INTERFACES
  • 53.1 Chemshell
  
  
• 54 THE TDHF AND TDKS PROGRAMS
• 55 ORBITAL MERGING
  • 55.1 Defining the input orbitals (ORBITAL)
  • 55.2 Moving orbitals to the output set (MOVE)
  • 55.3 Adding orbitals to the output set (ADD)
  • 55.4 Defining extra symmetries (EXTRA)
  • 55.5 Defining offsets in the output set (OFFSET)
  • 55.6 Projecting orbitals (PROJECT)
  • 55.7 Symmetric orthonormalization (ORTH)
  • 55.8 Schmidt orthonormalization (SCHMIDT)
  • 55.9 Rotating orbitals (ROTATE)
  • 55.10 Initialization of a new output set (INIT)
  • 55.11 Saving the merged orbitals
  • 55.12 Printing options (PRINT)
  • 55.13 Examples
  
  
• 56 MATRIX OPERATIONS
  • 56.1 Calling the matrix facility (MATROP)
  • 56.2 Loading matrices (LOAD)
  • 56.3 Saving matrices (SAVE)
  • 56.4 Adding matrices (ADD)
  • 56.5 Trace of a matrix or the product of two matrices (TRACE)
  • 56.6 Setting variables (SET)
  • 56.7 Multiplying matrices (MULT)
  • 56.8 Transforming operators (TRAN)
  • 56.9 Transforming density matrices into the MO basis (DMO)
  • 56.10 Diagonalizing a matrix DIAG
  • 56.11 Generating natural orbitals (NATORB)
  • 56.12 Forming an outer product of two vectors (OPRD)
  • 56.13 Combining matrix columns (ADDVEC)
  • 56.14 Forming a closed-shell density matrix (DENS)
  • 56.15 Computing a fock matrix (FOCK)
  • 56.16 Computing a coulomb operator (COUL)
  • 56.17 Computing an exchange operator (EXCH)
  • 56.18 Printing matrices (PRINT)
  • 56.19 Printing diagonal elements of a matrix (PRID)
  • 56.20 Printing orbitals (PRIO)
  • 56.21 Assigning matrix elements to a variable (ELEM)
  • 56.22 Reading a matrix from the input file (READ)
  • 56.23 Writing a matrix to an ASCII file (WRITE)
  • 56.24 Examples
  • 56.25 Exercise: SCF program
  
  
• A. Installation Guide
  • A..1 Obtaining the distribution materials
  • A..2 Installation of pre-built binaries
  • A..3 Installation from source files
  
  
• B. Recent Changes
  • B..1 New features of MOLPRO2010.1
  • B..2 New features of MOLPRO2009.1
  • B..3 New features of MOLPRO2008.1
  • B..4 New features of MOLPRO2006.1
  • B..5 New features of MOLPRO2002.6
  • B..6 New features of MOLPRO2002
  • B..7 Features that were new in MOLPRO2000
  • B..8 Facilities that were new in MOLPRO98
  
  
• C. Density functional descriptions
  • C..1 B86: X $\alpha\beta\gamma$
  • C..2 B86MGC: X $\alpha\beta\gamma$ with Modified Gradient Correction
  • C..3 B86R: X $\alpha\beta\gamma$ Re-optimised
  • C..4 B88: Becke 1988 Exchange Functional
  • C..5 B88C: Becke 1988 Correlation Functional
  • C..6 B95: Becke 1995 Correlation Functional
  • C..7 B97DF: Density functional part of B97
  • C..8 B97RDF: Density functional part of B97 Re-parameterized by Hamprecht et al
  • C..9 BR: Becke-Roussel Exchange Functional
  • C..10 BRUEG: Becke-Roussel Exchange Functional -- Uniform Electron Gas Limit
  • C..11 BW: Becke-Wigner Exchange-Correlation Functional
  • C..12 CS1: Colle-Salvetti correlation functional
  • C..13 CS2: Colle-Salvetti correlation functional
  • C..14 DIRAC: Slater-Dirac Exchange Energy
  • C..15 ECERF: Short-range LDA correlation functional
  • C..16 EXERF: Short-range LDA correlation functional
  • C..17 G96: Gill's 1996 Gradient Corrected Exchange Functional
  • C..18 HCTH120: Handy least squares fitted functional
  • C..19 HCTH147: Handy least squares fitted functional
  • C..20 HCTH93: Handy least squares fitted functional
  • C..21 LTA: Local $\\ tau$ Approximation
  • C..22 LYP: Lee, Yang and Parr Correlation Functional
  • C..23 MK00: Exchange Functional for Accurate Virtual Orbital Energies
  • C..24 MK00B: Exchange Functional for Accurate Virtual Orbital Energies
  • C..25 P86: .
  • C..26 PBEC: PBE Correlation Functional
  • C..27 PBEX: PBE Exchange Functional
  • C..28 PBEXREV: Revised PBE Exchange Functional
  • C..29 PW86: .
  • C..30 PW91C: Perdew-Wang 1991 GGA Correlation Functional
  • C..31 PW91X: Perdew-Wang 1991 GGA Exchange Functional
  • C..32 PW92C: Perdew-Wang 1992 GGA Correlation Functional
  • C..33 STEST: Test for number of electrons
  • C..34 TH1: Tozer and Handy 1998
  • C..35 TH2: .
  • C..36 TH3: .
  • C..37 TH4: .
  • C..38 THGFC: .
  • C..39 THGFCFO: .
  • C..40 THGFCO: .
  • C..41 THGFL: .
  • C..42 VSXC: .
  • C..43 VWN3: Vosko-Wilk-Nusair (1980) III local correlation energy
  • C..44 VWN5: Vosko-Wilk-Nusair (1980) V local correlation energy
  
  
• D. License information
  • D..1 BLAS
  • D..2 LAPACK
  • D..3 Boost
  
  
• Index