Wavefunction-Based Electrostatic-Embedding QM/MM Using CFOUR through MiMiC

Till Kirsch*, Jógvan Magnus Haugaard Olsen, Viacheslav Bolnykh, Simone Meloni, Emiliano Ippoliti, Ursula Rothlisberger, Michele Cascella, Jürgen Gauss

*Corresponding author for this work

Research output: Contribution to journalJournal articleResearchpeer-review

42 Downloads (Pure)


We present an interface of the wavefunction-based quantum chemical software CFOUR to the multiscale modeling framework MiMiC. Electrostatic embedding of the quantum mechanical (QM) part is achieved by analytic evaluation of one-electron integrals in CFOUR, while the rest of the QM/molecular mechanical (MM) operations are treated according to the previous MiMiC-based QM/MM implementation. Long-range electrostatic interactions are treated by a multipole expansion of the potential from the QM electron density to reduce the computational cost without loss of accuracy. Testing on model water/water systems, we verified that the CFOUR interface to MiMiC is robust, guaranteeing fast convergence of the self-consistent field cycles and optimal conservation of the energy during the integration of the equations of motion. Finally, we verified that the CFOUR interface to MiMiC is compatible with the use of a QM/QM multiple time-step algorithm, which effectively reduces the cost of ab initio MD (AIMD) or QM/MM-MD simulations using higher level wavefunction-based approaches compared to cheaper density functional theory-based ones. The new wavefunction-based AIMD and QM/MM-MD implementations were tested and validated for a large number of wavefunction approaches, including Hartree-Fock and post-Hartree-Fock methods like Møller-Plesset, coupled-cluster, and complete active space self-consistent field.
Original languageEnglish
JournalJournal of Chemical Theory and Computation
Issue number1
Pages (from-to)13–24
Number of pages12
Publication statusPublished - 2022


Dive into the research topics of 'Wavefunction-Based Electrostatic-Embedding QM/MM Using CFOUR through MiMiC'. Together they form a unique fingerprint.

Cite this