MOLPRO

The emphasis in the program is on highly accurate computations, with extensive treatment of the electron correlation problem through the multireference configuration interaction, coupled cluster and associated methods. Integral-direct local electron correlation methods reduce the increase of the computational cost with molecular size. Accurate ab initio calculations can then be performed for larger molecules. With new explicitly correlated methods the basis set limit can be very closely approached.

History

Molpro was designed and maintained by Wilfried Meyer and Peter Pulay in the late 1960s. At that moment, Pulay developed the first analytical gradient code called Hartree-Fock (HF),^[2][3][4] and Meyer researched his PNO-CEPA (pseudo-natural orbital coupled-electron pair approximation) methods.^[5][6] In 1980, Werner and Meyer developed a new state-averaged, quadratically convergent (MC-SCF) method, which provided geometry optimization for multireference cases.^[7] By the same year, the first internally contracted multireference configuration interaction (IC-MRCI) program was developed by Werner and Reinsch.^[8] About four years later (1984), Werner and Knowles developed on a new generation program called CASSCF (complete active space SCF). This new CASSCF program combined fast orbital optimization algorithms^[7] with determinant-based full CI codes,^[9] and additional, more general, unitary group configuration interaction (CI) codes. This resulted in the quadratically convergent MCSCF/CASSCF code called MULTI,^[10][11] which allowed modals to be optimized a weighted energy average of several states, and is capable of treating both completely general configuration expansions. In fact, this method is still available today. In addition to these organizational developments, Knowles and Werner started to cooperate on a new, more efficient, IC-MRCI method.^[12][13] Extensions for accurate treatments of excited states became possible through a new IC-MRCI method.^[14] In brief, the present IC-MRCI will be described as MRCI. These recently developed MCSCF and MRCI methods resulted in the basis of the modern Molpro. In the following years, a number of new programs were added. Analytic energy gradients can be evaluated with coupled-cluster calculations, density functional theory (DFT), as well as many other programs. These structural changes make the code more modular and easier to use and maintain, and also reduces the probability of input error.^[15]

See also

References