Abstract:Frozen density embedding (FDE) represents an embedding scheme in which environmental effects are included from first-principles calculations by considering the surrounding system explicitly by means of its electron density. In the present paper, we extend the full four-component relativistic Dirac-Kohn-Sham (DKS) method, as implemented in the BERTHA code, to include environmental and confinement effects with the FDE scheme (DKS-in-DFT FDE). The implementation, based on the auxiliary density fitting techniques, has been enormously facilitated by BERTHA's python API (PyBERTHA), which facilitates the interoperability with other FDE implementations available through the PyADF framework. The accuracy and numerical stability of this new implementation, also using different auxiliary fitting basis sets, has been demonstrated on the simple NH3-H2O system, in comparison with a reference nonrelativistic implementation. The computational performance has been evaluated on a series of gold clusters (Aun, with n = 2, 4, 8) embedded into an increasing number of water molecules (5, 10, 20, 40, and 80 water molecules). We found that the procedure scales approximately linearly both with the size of the frozen surrounding environment (consistent with the underpinnings of the FDE approach) and with the size of the active system (in line with the use of density fitting). Finally, we applied the code to a series of heavy (Rn) and super-heavy elements (Cn, Fl, Og) embedded in a C60 cage to explore the confinement effect induced by C60 on their electronic structure. We compare the results from our simulations, with respect to more-approximate models employed in the atomic physics literature. Our results indicate that the specific interactions described by FDE are able to improve upon the cruder approximations currently employed, and, thus, they provide a basis from which to generate more-realistic radial potentials for confined atoms.

Multi-State Formulation of the Frozen-Density Embedding Quasi-Diabatization Approach

Extending Density Functional Embedding Theory for Covalently Bonded Systems.

Frozen-Density Embedding for Including Environmental Effects in the Dirac-Kohn-Sham Theory: An Implementation Based on Density Fitting and Prototyping Techniques

Implementation of frozen density embedding in CP2K and OpenMolcas: CASSCF wavefunctions embedded in a Gaussian and plane wave DFT environment

Density Matrix Embedding Using Multiconfiguration Pair-Density Functional Theory

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

Density Fitted Fragment Embedding – Principles and Applications

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

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

Subsystem Density-Functional Theory (Update)

Fragment quantum embedding using the Householder transformation: a multi-state extension based on ensembles

Towards quantum-centric simulations of extended molecules: sample-based quantum diagonalization enhanced with density matrix embedding theory

Embedding approaches for bulk systems using projector-augmented waves

Can Density Matrix Embedding Theory with the Complete Activate Space Self-Consistent Field Solver Describe Single and Double Bond Breaking in Molecular Systems?

A Time-Dependent Random State Approach for Large-Scale Density Functional Calculations

A Practical Guide to Density Matrix Embedding Theory in Quantum Chemistry

Embedded density functional theory for covalently bonded and strongly interacting subsystems

Accelerating Density Matrix Embedding with Stochastic Density Fitting Theory – An Application to Hydrogen Bonded Clusters

Direct Diabatization of Electronic States by the Fourfold-Way: Including Dynamical Correlation by Multi-Configuration Quasidegenerate Perturbation Theory with Complete Active Space Self-Consistent-field Diabatic Molecular Orbitals

Embedding quantum statistical excitations in a classical force field

Understanding intermolecular interactions of large systems in ground state and excited state by using density functional based tight binding methods