# 4.2.32. MCPDFT¶

The MCPDFT program in Molcas performs multiconfiguration pair-density functional theory (MC-PDFT) calculations, as described in [91][92]. The MC-PDFT method involves two steps: (i) a CASSCF, RASSCF, or GASSCF wave function calculation to obtain the kinetic energy, classical Coulomb energy, total electron density, and on-top pair density; (ii) a post-SCF calculation of the remaining energy using an on-top density functional. In the current implementation, the on-top pair density functional is obtained by “translation” (t) of exchange-correlation functionals. Four translated functionals are currently available: tLSDA, tPBE, tBLYP, and trevPBE, in addition to the “fully-translated” (ft) variants [93]: ftLSDA, ftPBE, ftBLYP, and ftrevPBE. As multiconfigurational wave functions are used as input quantities, spin and space symmetry are correctly conserved.

The molecular orbitals and one- and two-body density matrices are read from the JOBIPH (or JOBOLD) file generated during a RASSCF run.

## 4.2.32.1. Dependencies¶

To start the MCPDFT module, the one-electron and two-electron integrals generated by SEWARD are required. For MC-PDFT calculations it is suggested to use a fine or ultrafine grid via the following input specifications (see the SEWARD section for further details):

&SEWARD
grid input
grid=ultrafine
end of grid input


Additionally, a JOBIPH file must be supplied to provide the molecular orbitals and one- and two-body density matrices.

## 4.2.32.2. Files¶

### 4.2.32.2.1. Input files¶

MCPDFT will use the following input files: ONEINT, ORDINT, RUNFILE, JOBOLD, JOBIPH.

If Cholesky options are selected, additional Cholesky-related files will also be used. Only JOBIPH or JOBOLD is needed; the code will first look for JOBOLD first and JOBIPH second.

### 4.2.32.2.2. Output files¶

JOBPDFT

This file is written in binary format and has the same structue of the JOBIPH file.

RUNFILE

The RUNFILE is updated with information from the MC-PDFT calculation.

MCDENS

This ASCII file is generated for MC-PDFT calculations. It contains spin densities, total density and on-top pair density values on grid (coordinates in a.u.).

## 4.2.32.3. Input¶

This section describes the input to the MCPDFT program in the Molcas program system. The input starts with the program name

&MCPDFT


The KSDFT is the only required keyword.

KSDFT

The functional choice follows. Currently available functionals are: tPBE, tBLYP, tLSDA, trevPBE, tOPBE, ftPBE, ftBLYP, ftLSDA, ftrevPBE and ftOPBE.

DFCF

Use this keyword to scale the exchange terms and/or correlation terms of the functional requested. This keyword should be followed by the scaling factor for the exchange terms and the scaling factor for the correlation terms, separated by a space. If the values are 1.0 (default), then the original functional is used. For an HLE-type functional, use 1.25 (for exchange) and 0.5 (for correlation). Example: DFCF=1.25 0.5

### 4.2.32.3.1. Input example¶

The following example shows the input to the RASSCF and MCPDFT programs for a calculation on the water molecule. The tPBE functional is used. The calculation is performed in $$C_{2v}$$ symmetry (symmetries: $$a_1$$, $$b_2$$, $$b_1$$, $$a_2$$, where the two last species are antisymmetric with respect to the molecular plane). Inactive orbitals are 1$$a_1$$ (oxygen 1s) 2$$a_1$$ (oxygen 2s) and 1$$b_1$$ (the $$\pi$$ lone-pair orbital). Two bonding and two anti-bonding $$\ce{OH}$$ orbitals are active, $$a_1$$ and $$b_2$$ symmetries. The calculation is performed for the $$^1A_1$$ ground state. Note that no information about basis set, geometry, etc. has to be given. Such information is supplied by the SEWARD integral program via the one-electron integral file ONEINT.

&RASSCF
Title= Water molecule. Active orbitals OH and OH* in both symmetries
Spin     = 1
Symmetry = 1
Inactive = 2 0 1 0
Ras2     = 2 2 0 0

&MCPDFT
KSDFT=TPBE


The first RASSCF run is a standard CASSCF calculation that leads to variationally optimized orbitals and CI coefficients. The MC-PDFT run will use the orbitals and density matrices optimized during the preceding RASSCF run.