Skip to main content

MC-PDFT

Multiconfiguration pair-density functional theory (MC-PDFT) is a cost-effective post-SCF method which is able to recover the missing dynamic correlation from multiconfigurational wave functions generated through multireference methods such as CASCI or CASSCF. MC-PDFT computes a correct energy of a given reference wave function using an energy expression which contains a on-top energy functional term which depends on both the wave function’s density and on-top pair density. This approach adds minimal overhead cost to the reference multireference method while obtaining accuracy which rivals more traditional (and expensive) multireference perturbation theories such as CASPT2 or NEVPT2. Our group continues to develop MC-PDFT to tackle a wide variety of problems in photochemistry, photodynamics, catalysis, and biochemistry and integrate MC-PDFT with our automated active space selection scheme (DVS) and efficient multireference wave function ansatz (LASSCF and DMET).

The MC-PDFT program is available in several packages(Open-Molcas, Gamess, PySCF and SHARC-md) as described below.

1. Pyscf-forge

Pyscf-forge is an extension module for the PySCF electronic structure package. It can be installed using the following command

pip install pyscf-forge

To access the newest features of Pyscf-forge, you can install them from the git repository by running the command

pip install git+https://github.com/pyscf/pyscf-forge

Developers: The MC-PDFT module in PySCF-forge has been developed by:

  • Matthew R Hermes (University of Chicago)
  • Dayou Zhang (University of Minnesota)
  • Aleksandr Lykhin (University of Chicago)
  • Thais R Scott (University of Chicago)
  • Matthew R Hennefarth (University of Chicago)

Features:

  • MC-PDFT total electronic energy calculations for wave functions of various types.
    • CASCI
    • CASSCF
    • State-averaged CASSCF (including “mixed” solver with different spins and/or point groups)
  • Multi-state Extensions of MC-PDFT
  • On-the-fly generation of on-top density functionals from underlying KS-DFT ‘LDA’ or ‘GGA’ exchange-correlation functionals as defined in Libxc.
    • Translated functionals: JCTC 2014, 10, 3669
    • Fully-translated functionals: JCTC 2015, 11, 4077
    • Global hybrid functionals: JPCL 2020, 11, 10158 and JCTC 2020, 16, 2274
    • Notes:
      1. Translation of ‘meta’ KS-DFT functionals which depend on the kinetic energy density and/or Laplacian is not supported.
      2. Range-separated hybrid on-top functionals are not supported.
      3. Translation of functionals defined as global hybrids at the Libxc or PySCF level is not supported, except for ‘tPBEo’ and ‘ftPBE0’.
  • Additional properties
    • Decomposition of total electronic energy into core, Coulomb, on-top components
    • Analytical nuclear gradients (non-hybrid functionals only) for:
      1. Single-state CASSCF wave function: JCTC 2018, 14, 126
      2. State-averaged CASSCF wave functions: JCP 2020, 153, 014106
      3. CMS-PDFT: Mol Phys 2022, 120
      4. L-PDFT J Chem Theory Comput, 2024, 9, 20
    • Permanent electric dipole moment (non-hybrid functionals only) for:
      1. Single-state CASSCF wave function: JCTC 2021, 17, 7586
      2. State-averaged CASSCF wave functions
      3. CMS-PDFT
    • Transition electric dipole moment (non-hybrid functionals only) for:
      1. CMS-PDFT
  • Multi-configuration density-coherence functional theory (MC-DCFT) total energy: JCTC 2021, 17, 2775

License: PySCF and PySCF-forge are both licensed under the Apache License, Version 2.0 (the “License”); you may not use this file except in compliance with the License. You may obtain a copy of the License at http://www.apache.org/licenses/LICENSE-2.0 Unless required by applicable law or agreed to in writing, software distributed under the License is distributed on an “AS IS” BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. See the License for the specific language governing permissions and limitations under the License.

2. OpenMolcas

See: https://comp.chem.umn.edu/openmolcas

3. GAMESS

See: https://comp.chem.umn.edu/MC-PDFT_in_GAMESS

4. SHARC-md

The SHARC-md program suite is an ab-initio molecular dynamics software package developed to study the excited-state dynamics of molecules. The git repository of the latest version (and previous release versions) of SHARC-md is available at https://github.com/sharc-md/sharc. You may also download the latest version of the Gagliardi Group mirror with recent developments from those in the Gagliardi Group at https://github.com/GagliardiGroup/sharc

License: SHARC: Surface Hopping including Arbitrary Couplings Copyright (C) 2024 SHARC Developers
This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program. If not, see https://www.gnu.org/licenses.

Developers: The MC-PDFT module within SHARC-md has been developed by:

  • Matthew R Hennefarth (University of Chicago)
  • Paul B Calio (University of Chicago)

Features:

Interfaced with the PySCF/PySCF-forge electronic structure packages

  • Potential energy surfaces (and gradients) generated by the following electronic structure methods:
    • SA-CASSCF
    • MC-PDFT
    • L-PDFT
    • CMS-PDFT
  • Nonadiabatic coupling vectors available for the following methods:
    • SA-CASSCF
    • CMS-PDFT

Modified OpenMolcas interface to allow MC-PDFT methods for the on-the-fly electronic structure method.

  • Potential energy surfaces (and gradients) generated by the following electronic structure methods:
    • MC-PDFT
    • CMS-PDFT
    • XMS-PDFT (numerical gradients only)
  • Nonadiabatic coupling vectors available for the following methods:
    • CMS-PDFT

References

  1. JCTC 2020, 16, 7444: http://dx.doi.org/10.1021/acs.jctc.0c00908
  2. JCTC 2014, 10, 3669: http://dx.doi.org/10.1021/ct500483t
  3. JCTC 2015, 11, 4077: http://dx.doi.org/10.1021/acs.jctc.5b00609
  4. JPCL 2020, 11, 10158: http://dx.doi.org/10.1021/acs.jpclett.0c02956
  5. JCTC 2020, 16, 2274: http://dx.doi.org/10.1021/acs.jctc.9b01178
  6. JCTC 2018, 14, 126: http://dx.doi.org/10.1021/acs.jctc.7b00967
  7. JCP 2020, 153, 014106: http://dx.doi.org/10.1063/5.0007040
  8. JCTC 2021, 17, 7586: http://dx.doi.org/10.1021/acs.jctc.1c00915
  9. JCTC 2021, 17, 2775: http://dx.doi.org/10.1021/acs.jctc.0c01346
  10. Mol Phys 2022, 120: http://dx.doi.org/10.1080/00268976.2022.2110534
  11. Faraday Discuss 2020, 224, 348-372: http://dx.doi.org/10.1039/D0FD00037J
  12. JCTC 2023, 19, 3172: https://dx.doi.org/10.1021/acs.jctc.3c00207
  13. LASSCF and DMET: https://gagliardigroup.uchicago.edu/lasscf-and-cas-dmet-in-pyscf/
  14. DVS: https://dx.doi.org/10.1021/acs.jctc.3c00792
  15. libxc: https://libxc.gitlab.io
  16. PySCF-forge: https://github.com/pyscf/pyscf-forge
  17. PySCF: https://pyscf.org
  18. SHARC-md: https://sharc-md.org
  19. GAMESS: https://www.msg.chem.iastate.edu/gamess/
  20. OpenMolcas: https://molcas.gitlab.io/OpenMolcas/sphinx/