Multistate, Polarizable QM/MM Embedding Scheme Based on the Direct Reaction Field Method: Solvatochromic Shifts,...

Area of Science:

• Computational Chemistry
• Theoretical Chemistry
• Quantum Chemistry

Background:

• Accurate modeling of solvation effects is crucial for understanding electronic states in condensed phases.
• Existing QM/MM methods often struggle to capture differential solvation of ground and excited states.
• Polarizable embedding schemes are essential for realistic simulations of molecular properties in solution.

Purpose of the Study:

• To implement and validate an integral-exact reformulation of the direct reaction field (IEDRF) method for QM/MM simulations.
• To enable accurate treatment of polarization and dispersion interactions between QM and MM regions.
• To investigate the solvent shifts of absorption energies for donor-acceptor dyes.

Main Methods:

• Developed an integral-exact polarizable embedding scheme (IEDRF) within the TeraChem software.
• Utilized graphical processing unit (GPU) acceleration for enhanced computational efficiency.
• Combined the IEDRF method with Hartree-Fock, time-dependent density functional theory, and complete active space self-consistent field electronic structure methods.
• Derived and tested analytical gradients and nonadiabatic coupling vectors.

Main Results:

• The IEDRF method accurately accounts for differential solvation of ground and excited states.
• GPU acceleration significantly speeds up QM/MM calculations.
• Optimized a conical intersection between excited states in a polarizable solvent.
• Achieved quantitative agreement (0.03 eV) with full-system calculations for solvent shifts of absorption energies.
• Demonstrated good agreement with experimental absorption maxima.

Conclusions:

• The implemented IEDRF method provides a robust and accurate approach for QM/MM simulations.
• The method significantly improves the prediction of solvatochromic shifts.
• This work paves the way for more reliable theoretical studies of electronic processes in solution.