Abstract:We recently introduced a polarizable embedding scheme based on an integral-exact reformulation of the direct reaction field method (IEDRF) that accounts for differential solvation of ground and excited states in QM/MM simulations. The polarization and dispersion interactions between the quantum-mechanical (QM) and molecular-mechanical (MM) regions are described by the DRF Hamiltonian, while the Pauli repulsion between explicitly treated QM electrons and the implicit electron density around MM atoms is modeled with effective core potentials. A single Hamiltonian is used for all electronic states, so that Born-Oppenheimer states belonging to the same geometry are orthogonal and state crossings are well-defined. In this work, we describe the implementation of the method using graphical processing unit acceleration in TeraChem, where it is combined with multiple electronic structure methods, including Hartree-Fock, Time-dependent Density Functional Theory, and Complete Active Space Self-Consistent Field. In contrast with older implementations of the DRF method, integrals of the polarization operators are evaluated exactly. Expressions for ingredients needed to construct analytical gradients and non-adiabatic coupling vectors are derived and tested by optimizing a conical intersection between two excited states in the presence of a polarizable solvent shell. The method is applied to estimating the solvent shifts of absorption energies of a series of donor-acceptor dyes having low-lying charge-transfer states. Even for a non-polar solvent such as n-hexane, the inclusion of its static polarizability leads to non-negligible shifts that improve the agreement to essentially quantitative levels (0.03 eV) with full-system calculations. Good agreement with the positions of experimental absorption maxima measured in solution is also observed.

What problem does this paper attempt to address?

### Problems the paper attempts to solve This paper aims to solve the problem in quantum mechanics/molecular mechanics (QM/MM) simulations that different electronic states react differently to the polarized environment. Specifically, the paper proposes a multi - state, polarized embedding scheme based on the Direct Reaction Field (DRF) method to more accurately describe the influence of solvents on photochemical reactions, especially on the excited state. ### Background and motivation 1. **Environmental effects in photochemical reactions**: - Photochemical reactions usually focus on smaller molecular components (chromophores) which are responsible for absorbing light energy. - In the liquid or condensed phase, the environment has a significant impact on the photophysical properties of chromophores. - Solvent effects can be partly explained by modeling the solvent as an unstructured dielectric continuum or a viscous medium, but these simple models cannot capture complex solvent effects in some cases. 2. **Limitations of existing methods**: - **Linear response method**: It only considers the interaction between the ground state and the excited state and cannot fully capture the influence of the solvent on the long - range charge transfer (CT) state. - **State - specific method**: Each electronic state needs to be calculated separately, resulting in different electronic states no longer being eigenfunctions of the same Hamiltonian, making the transition matrix elements and state crossings undefined. ### Solutions 1. **Direct Reaction Field method (DRF)**: - The DRF method overcomes the deficiencies of the linear response and state - specific methods by directly adding the reaction field to the Hamiltonian. - This method uses a unified Hamiltonian to describe all electronic states, ensuring the orthogonality of different electronic states and the clear definition of state crossings. 2. **Technical improvements**: - **Exact integration**: The integrals of all polarization operators are calculated exactly, avoiding the errors introduced by approximate methods. - **GPU acceleration**: Implemented on the TeraChem platform, using GPU acceleration to improve computational efficiency. - **Analytical gradient**: The analytical gradient expression is derived and implemented for optimizing conical intersections and other non - adiabatic coupling vectors. ### Applications and verification 1. **Solvatochromic shift**: - The solvatochromic shifts of the absorption energies of a series of donor - acceptor dyes in different solvents are calculated. - Even for non - polar solvents (such as hexane), its static polarizability also causes a non - negligible shift, improving the consistency with the full - system calculation. 2. **Comparison with experimental data**: - It is compared with the experimentally measured absorption maximum positions in solution, and the results show good consistency. ### Conclusion This paper provides a more accurate and efficient QM/MM embedding scheme through the improved DRF method, which can better describe the influence of solvents on photochemical reactions, especially the solvent effects in the excited state. This provides a powerful tool for studying photochemical processes in complex systems.

Multistate, polarizable QM/MM embedding scheme based on the direct reaction field method: Solvatochromic shifts, analytical gradients and optimization of conical intersections in solution

Multistate, Polarizable QM/MM Embedding Scheme Based on the Direct Reaction Field Method: Solvatochromic Shifts, Analytical Gradients and Optimizations of Conical Intersections in Solution

Efficient polarizable QM/MM using the direct reaction field Hamiltonian with electrostatic potential fitted multipole operators

Tinktep: A Fully Self-Consistent, Mutually Polarizable Qm/Mm Approach Based On The Amoeba Force Field

Fully Polarizable Multiconfigurational Self-consistent Field/Fluctuating Charge Approach

A polarizable CASSCF/MM approach using the interface between OpenMMPol library and CFour

Fully Polarizable QM/Fluctuating Charge Approach to Two-Photon Absorption of Aqueous Solutions

DRF90: a polarizable force field

Analytic gradients for the electrostatic embedding QM/MM in periodic boundary conditions using particle-mesh Ewald sums and electrostatic potential fitted charge operators

Nonequilibrium Solvent Polarization Effects in Real-Time Electronic Dynamics of Solute Molecules Subject to Time-Dependent Electric Fields: A New Feature of the Polarizable Continuum Model

A fully-polarizable KS-DFT/AMOEBA embedding scheme for plane wave basis sets through the MiMiC framework

Seamless integration of GEM, a density based-force field, for QM/MM simulations via LICHEM, Psi4 and Tinker-HP

Hybrid Combination of Quantum Mechanics with Quantum-Based Polarizable Reactive Force Field for Large Scale Full Solvent Simulations of Electrocatalysis

Incorporating surface polarization effects into large-scale coarse-grained Molecular Dynamics simulation

Self-consistent Quantum Linear Response with a Polarizable Embedding environment

A Polarizable QM/MM Model That Combines the State-Averaged CASSCF and AMOEBA Force Field for Photoreactions in Proteins

Coupling Molecular Density Functional Theory with Converged Selected Configuration Interaction Methods to Study Excited states in Aqueous Solution

Assessment of embedding schemes in a hybrid machine learning/classical potentials (ML/MM) scheme

A Diffusion Quantum Monte Carlo Approach to the Polaritonic Ground State

Effect of the QM Size, Basis Set, and Polarization on QM/MM Interaction Energy Decomposition Analysis

Solvation effects on halides core spectra with Multilevel Real-Time quantum embedding