Prediction of the aqueous redox properties of functionalized quinones using a new QM/MM variational formulation

We recently proposed a method coupling quantum mechanics (QM) methods and molecular density functional theory (MDFT) to describe mixed quantum-classical systems [J. Chem. Phys. 161, 014113 (2024)]. This approach is particularly appropriate to account for solvent effect into QM calculations. We introduce a new variational formulation for the grand potential of a mixed quantum-classical system. Within the Born-Oppenheimer approximation and neglecting electronic entropy, the quantum solute is described by a product of electronic and nuclear density matrices, both depending parametrically on coordinates of the classical solvent. It can then be shown that a functional of the total density matrix satisfies a variational principle for the grand potential. Using a mean-field approximation, we express the grand potential of the mixed quantum-classical system as a variational problem which depends only on the nuclear density matrix, which experiences an external field generated by the electronic and classical one-particle densities. In practice, the grand potential is computed by a series of coupled classical and quantum DFT calculations, together with geometry optimization.