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

The theory of analytic gradients is presented for the projector-based density functional theory (DFT) embedding approach utilizing the Huzinaga-equation. The advantages of the Huzinaga-equation-based formulation are demonstrated. In particular, it is shown that the projector employed does not appear in the Lagrangian, and the potential risk of numerical problems is avoided at the evaluation of the gradients. The efficient implementation of the analytic gradient theory is presented for the approaches where hybrid DFT, second-order Moller-Plesset perturbation theory (MP2), or double hybrid DFT is embedded in lower-level DFT environments. To demonstrate the applicability of the method and to gain insight into its accuracy, it is applied to equilibrium geometry optimizations, transition state searches, and potential energy surface scans. Our results show that bond lengths and angles converge rapidly with the size of the embedded system. While providing structural parameters close to high-level quality for the embedded atoms, the embedding approach has the potential to relax the coordinates of the environment as well. Our demonstrations on a 171-atom zeolite and a 570-atom protein system show that the Huzinaga-equation-based embedding can accelerate (double) hybrid gradient computations by an order of magnitude with sufficient active regions and enables affordable force evaluations or geometry optimizations for molecules of hundreds of atoms.