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 approximations that introduce ambiguities. This series of papers aims to establish a rigorous theoretical framework to guide the development of a QM/cDFT formulation that clarifies the approximations involved. In this first paper, we establish a comprehensive density functional theory (DFT) framework for mixed quantum (QM)-classical (MM) systems within the canonical ensemble. We derive the equilibrium density matrix for a system comprising N qm QM and N mm MM particles and, based on this derivation, propose a variational formulation of the Helmholtz free energy in terms of the full QM-MM density matrix. By leveraging permutational symmetry and thanks to constrained-search methods, we reformulate the variational problem using functionals of objects having reduced dimensionalities. Three different formulations are presented, each using quantities with progressively lower dimensionalities. This ultimately enables computation of the Helmholtz free energy using only the QM and MM one-body densities, generalizing both cDFT and electronic DFT (eDFT) for QM-MM systems. In a subsequent paper, this formalism will be extended to the grand canonical ensemble, allowing classical particle number fluctuations. This extension, in which cDFT is most naturally formulated, provides a framework for studying solvation problems.

Canonical and Grand Canonical Ensemble Expectation Values from Quantum Monte Carlo Simulations

Theory and Monte Carlo simulation of the ideal gas with shell particles in the canonical, isothermal-isobaric, grand canonical, and Gibbs ensembles

Canonical vs. Grand Canonical Ensemble for Bosonic Gases under Harmonic Confinement

Gaussian concentration bound and Ensemble equivalence in generic quantum many-body systems including long-range interaction

High-efficiency quantum Monte Carlo algorithm for extracting entanglement entropy in interacting fermion systems

Nonequilibrium grand-canonical ensemble built from a physical particle reservoir

Computing entanglement entropy of interacting fermions with quantum Monte Carlo: Why we failed and how to get it right

An integral algorithm of exponential observables for interacting fermions in quantum Monte Carlo simulation

Variational Microcanonical Estimator

Canonical Ensemble vs. Grand Canonical Ensemble in the Description of Multicomponent Bosonic Systems

Incremental SWAP Operator for Entanglement Entropy: Application for Exponential Observables in Quantum Monte Carlo Simulation

Grand Canonical Adaptive Resolution Simulation for Molecules with Electrons: A Theoretical Framework based on Physical Consistency

A quantum Monte Carlo algorithm for arbitrary spin-1/2 Hamiltonians

Non-equilibrium quantum Monte Carlo algorithm for stabilizer Renyi entropy in spin systems

A quantum expectation identity: Applications to statistical mechanics

Titration in Canonical and Grand-Canonical Ensembles

Canonical titration simulations

Extended Ensemble Monte Carlo

A Rigorous Theory of Quantum-Classical Systems at Finite Temperature, Part I: Exact Functional Formulation in the Canonical Ensemble