A variational formulation of the free energy of mixed quantum-classical systems: coupling classical and electronic density functional theories

Combining classical density functional theory (cDFT) with quantum mechanics (QM) methods offers a computationally efficient alternative to traditional QM/molecular mechanics (MM) approaches for modeling mixed quantum-classical systems at finite temperatures. However, both QM/MM and QM/cDFT rely on somewhat ambiguous approximations, the two major ones being: i) the definition of the QM and MM regions as well as the description of their coupling, and ii) the choice of the methods and levels of approximation made to describe each region. This paper addresses the second point and develop an exact theoretical framework that allows us to clarify the approximations involved in the QM/cDFT formulation. We establish a comprehensive density functional theory (DFT) framework for mixed quantum-classical systems within the canonical ensemble. We start by recalling the expression of the adiabatic equilibrium density matrix for a mixed system made of Nqm quantum and Nmm classical particles. Then, we propose a variational formulation of the Helmholtz free energy in terms of the full, non-equilibrium, QM/MM density matrix. Taking advantage of permutational symmetry and thanks to constrained-search methods, we reformulate the computation of the Helmholtz free energy using only the quantum and classical one-body densities.This paper generalizes both cDFT and electronic DFT (eDFT) to QM/MM systems. We then reformulate the functional to make the standard eDFT and cDFT Levy-Lieb functionals explicitly appear, together with a new universal correlation functional for QM/MM systems. A mean-field approximation is finally introduced in the context of solvation problems and we discuss its connection with several existing mixed cDFT-eDFT schemes. An extension to the semi-grand canonical ensemble, where the number of classical particles is allowed to fluctuate, is provided in the supplementary materials.

💡 Research Summary

The paper presents a rigorous variational framework that unifies classical density‑functional theory (cDFT) and electronic density‑functional theory (eDFT) for mixed quantum‑classical (QM/MM) systems within the canonical ensemble. Starting from the full equilibrium density matrix of a system containing N_qm quantum particles (electrons and possibly quantum nuclei) and N_mm classical particles (solvent molecules, ions, etc.), the authors employ a partial Wigner transformation to express the quantum operators in a phase‑space representation that is isomorphic to the classical (Q,P) space. This allows them to write the Helmholtz free energy as a functional of a general, possibly non‑equilibrium, density matrix ρ̂, namely
 F