Abstract:Chemical bonding is a ubiquitous concept in chemistry and it provides a common basis for experimental and theoretical chemists to explain and predict the structure, stability and reactivity of chemical species. Among others, the Energy Decomposition Analysis (EDA, also known as the Extended Transition State method) in combination with Natural Orbitals for Chemical Valence (EDA-NOCV) is a very powerful tool for the analysis of the chemical bonds based on a charge and energy decomposition scheme within a common theoretical framework. While the approach has been applied in a variety of chemical contexts, the current implementations of the EDA-NOCV scheme include relativistic effects only at scalar level, so simply neglecting the spin-orbit coupling effects and de facto limiting its applicability. In this work, we extend the EDA-NOCV method to the relativistic four-component Dirac-Kohn-Sham theory that variationally accounts for spin-orbit coupling. Its correctness and numerical stability have been demonstrated in the case of simple molecular systems, where the relativistic effects play a negligible role, by comparison with the implementation available in the ADF modelling suite (using the non-relativistic Hamiltonian and the scalar ZORA approximation). As an illustrative example we analyse the metal-ethylene coordination bond in the group 6-element series (CO)$_5$TM-C$_2$H$_4$, with TM =Cr, Mo, W, Sg, where relativistic effects are likely to play an increasingly important role as one moves down the group. The method provides a clear measure (also in combination with the CD analysis) of the donation and back-donation components in coordination bonds, even when relativistic effects, including spin-orbit coupling, are crucial for understanding the chemical bond involving heavy and superheavy atoms.

Revealing strong interactions with the reduced density gradient: a benchmark for covalent, ionic and charge-shift bonds

Formulating reduced density gradient approaches for noncovalent interactions

NCIPLOT and the analysis of noncovalent interactions using the reduced density gradient

Interaction Region Indicator (IRI): A Very Simple Real Space Function Clearly Revealing Both Chemical Bonds and Weak Interactions

Interpretation of the reduced density gradient

Noncovalent Interaction Analysis in Fluctuating Environments

A New Angle on Benchmarking Noncovalent Interactions

Chemical bond analysis for the entire periodic table: Energy Decomposition and Natural Orbitals for Chemical Valence in the Four-Component Relativistic Framework

Correlations between Density-Based Bond Orders and Orbital-Based Bond Energies for Chemical Bonding Analysis

Gradient Bundle Analysis: A Full Topological Approach to Chemical Bonding

Observing Noncovalent Interactions in Experimental Electron Density for Macromolecular Systems: A Novel Perspective for Protein–Ligand Interaction Research

Simultaneous identification of strong and weak interactions with Pauli energy, Pauli potential, Pauli force, and Pauli charge

Relative cooperativity in neutral and charged molecular clusters using QM/MM calculations

Visualization of Electron Density Changes along Chemical Reaction Pathways

Reaction Coordinates by Nonlinear Dimensionality Reduction

Multi-site reaction dynamics through multi-fragment density matrix embedding

Charge transfer interaction revisited by a Fermi-Dirac derived approach

Distinguishing between chemical bonding and physical binding using electron localization function (ELF)

Benchmarking DFT-based excited-state methods for intermolecular charge-transfer excitations

DeepNCI: DFT Noncovalent Interaction Correction with Transferable Multimodal Three-Dimensional Convolutional Neural Networks