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.

Analytical Gradient Using Cluster-in-Molecule RI-MP2 Method for the Geometry Optimizations of Large Systems

Analytical Energy Gradients for the Cluster-in-Molecule MP2 Method and Its Application to Geometry Optimizations of Large Systems

High-performance multi-GPU analytic RI-MP2 energy gradients

Analytical Gradient Theory for Spin-Free State-Averaged Second-Order Driven Similarity Renormalization Group Perturbation Theory (SA-DSRG-MRPT2) and Its Applications for Conical Intersection Optimizations

Analytical Gradient Theory for Resolvent-Fitted Second-Order Extended Multiconfiguration Perturbation Theory (XMCQDPT2)

Localized Resolution of Identity Approach to the Analytical Gradients of Random-Phase Approximation Ground-State Energy: Algorithm and Benchmarks

Multi-GPU RI-HF Energies and Analytic Gradients─Toward High-Throughput Ab Initio Molecular Dynamics

Analytical Energy Gradient Evaluation in Relativistic and Nonrelativistic Density Functional Calculations

Analytic gradients for equation-of-motion coupled cluster with single, double, and perturbative triple excitations

Development of analytic gradients for the Huzinaga quantum embedding method and its applications to large-scale hybrid and double hybrid DFT forces

A refined cluster-in-molecule local correlation approach for predicting the relative energies of large systems

Massively parallel implementation of gradients within the random phase approximation: Application to the polymorphs of benzene

The Relative Energies of Polypeptide Conformers Predicted by Linear Scaling Second-Order Møller-Plesset Perturbation Theory

Improved Cluster-in-molecule Local Correlation Approach for Electron Correlation Calculation of Large Systems.

Massive-Parallel Implementation of the Resolution-of-Identity Coupled-Cluster Approaches in the Numeric Atom-Centered Orbital Framework for Molecular Systems

Multi-GPU RI-HF Energies and Analytic Gradients $-$ Towards High Throughput Ab Initio Molecular Dynamics

Fully Optimized Implementation of the Cluster-in-molecule Local Correlation Approach for Electron Correlation Calculations of Large 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

Analytical evaluation of ground state gradients in quantum electrodynamics coupled cluster theory