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

TL;DR

This paper develops an exact variational framework within density functional theory for mixed quantum-classical systems, clarifying approximations and generalizing cDFT and eDFT to improve modeling of such systems at finite temperatures.

Contribution

It introduces a comprehensive DFT-based theoretical framework for quantum-classical systems, explicitly incorporating Levy-Lieb functionals and a new correlation functional for QM/MM.

Findings

Reformulation of Helmholtz free energy using quantum and classical densities

Generalization of cDFT and eDFT to mixed QM/MM systems

Introduction of a mean-field approximation for solvation problems

Abstract

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.