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.

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

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

Extending Density Functional Embedding Theory for Covalently Bonded Systems.

Density Fitted Fragment Embedding – Principles and Applications

A Reduced Scaling Approach to Obtaining Embedding Orbital Coefficients: Benchmarking DMET-CCSD(T) Relative Energies in Water Clusters

Subspace Density Matrix Functional Embedding Theory: Theory, Implementation, and Applications to Molecular Systems

Enhancing robustness and efficiency of density matrix embedding theory via semidefinite programming and local correlation potential fitting

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

A Practical Guide to Density Matrix Embedding Theory in Quantum Chemistry

Efficient Multi-Configurational Quantum Chemistry Approach to Single-Ion Magnets Based on Density Matrix Embedding Theory

Embedding theory contributions to average atom models for warm dense matter

Distributed Memory, GPU Accelerated Fock Construction for Hybrid, Gaussian Basis Density Functional Theory

Low-rank Approximations to Accelerate Hybrid Functional Enabled Real-Time Time-Dependent Density Functional Theory Within Plane Waves

Efficient Implementation of Ab Initio Quantum Embedding in Periodic Systems: Density Matrix Embedding Theory

Projected density matrix embedding theory with applications to the two-dimensional Hubbard model

Static Embedding with Pair Coupled Cluster Doubles Based Methods

Overlapped Embedded Fragment Stochastic Density Functional Theory for Covalently Bonded Materials

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

Embedded density functional theory for covalently bonded and strongly interacting subsystems

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