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.

A Diffusion Quantum Monte Carlo Approach to the Polaritonic Ground State

Path Integral Approach for Multi-Dimensional Polarons in A Symmetric Quantum Dot

Polaritonic Unitary Coupled Cluster for Quantum Computations

A comprehensive theory for relativistic polaritonic chemistry: a four components ab initio treatment of molecular systems coupled to quantum fields

Understanding the polaritonic ground state in cavity quantum electrodynamics

Ab initio methods for polariton chemistry

Polaritonic Chemistry using the Density Matrix Renormalization Group Method

Towards polaritonic molecular orbitals for large molecular systems

Ab Initio Approaches to Simulate Molecular Polaritons: Properties and Quantum Dynamics

A perspective on ab initio modeling of polaritonic chemistry: The role of non-equilibrium effects and quantum collectivity

Polaritonic response theory for exact and approximate wave functions

Cavity Quantum Electrodynamics Complete Active Space Configuration Interaction Theory

Simulating polaritonic ground states on noisy quantum devices

Toward Polaritonic Molecular Orbitals for Large Molecular Systems

Quantized Embedding Approaches for Collective Strong Coupling -- Connecting ab initio and macroscopic QED to Simple Models in Polaritonics

Polaritonic coupled-cluster theory

Quantum Dynamics Simulations of the Polariton Transport Process

Theoretical Advances in Polariton Chemistry and Molecular Cavity Quantum Electrodynamics

Understanding Polaritonic Chemistry from Ab Initio Quantum Electrodynamics

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