Full Parallel Implementation of an All-Electron Four-Component Dirac–Kohn–Sham Program

Cooperative small molecule activation by apolar and weakly polar bonds through the lens of a suitable computational protocol

Small molecule activation processes are central in chemical research and cooperativity is a valuable tool for the fine-tuning of the efficiency of these reactions. In this contribution, we discuss recent and remarkable examples in which activation processes are mediated by bimetallic compounds featuring apolar or weakly polar metal-metal bonds. Relevant experimental breakthroughs are thoroughly analyzed from a computational perspective. We highlight how the rational and non-trivial application of selected computational approaches not only allows rationalization of the observed reactivities but also inferring of general principles applicable to activation processes, such as the breakdown of the structure-reactivity relationship in carbon dioxide activation in a cooperative framework. We finally provide a simple yet unbiased computational protocol to study these reactions, which can support experimental advances aimed at expanding the range of applications of apolar and weakly polar bonds as catalysts for small molecule activation.

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

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 NH₃-H₂O system, in comparison with a reference nonrelativistic implementation. The computational performance has been evaluated on a series of gold clusters (Au_n, 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 C₆₀ cage to explore the confinement effect induced by C₆₀ 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.

Local Charge-Displacement Analysis: Targeting Local Charge-Flows in Complex Intermolecular Interactions

Charge-displacement (CD) analysis has recently proven to be a simple and powerful scheme for quantitatively analyzing the profile of the charge redistribution occurring upon intermolecular interactions along a given interaction axis. However, when two molecular fragments bind through complex interactions involving multiple concurrent charge flows, ordinary CD analysis is capable of providing only an averaged picture of the related charge-flow profiles and no detailed information on each of them. In this article, we combine CD analysis with a Hirshfeld partitioning of the molecular charge redistribution for a local analysis on focused portions of the molecule, allowing for a detailed characterization of one charge flow at a time. The resulting scheme - the local charge-displacement (LCD) analysis - is tested on the intriguing case of the dimethyl sulfide (DMS)-sulfur dioxide (SO2) complex, characterized by concurrent charge flows relating to a sulfur-sulfur homo-chalcogen interaction and a pair of hydrogen bonds. The LCD scheme is then applied to the analysis of multiple hydrogen bonding in the acetic acid dimer, of base-pairing interactions in DNA, and of ambifunctional hydrogen bonding in the ammonia-pyridine complex.

Charge-Flow Profiles along Curvilinear Paths: A Flexible Scheme for the Analysis of Charge Displacement upon Intermolecular Interactions

The Charge-Displacement (CD) analysis has proven to be a powerful tool for a quantitative characterization of the electron-density flow occurring upon chemical bonding along a suitably chosen interaction axis. In several classes of interesting intermolecular interactions, however, an interaction axis cannot be straightforwardly defined, and the CD analysis loses consistency and usefulness. In this article, we propose a general, flexible reformulation of the CD analysis capable of providing a quantitative view of the charge displacement along custom curvilinear paths. The new scheme naturally reduces to ordinary CD analysis if the path is chosen to be a straight line. An implementation based on a discrete sampling of the electron densities and a Voronoi space partitioning is described and shown in action on two test cases of a metal-carbonyl and a pyridine-ammonia complex.

Spin-orbit coupling is the key to unraveling intriguing features of the halogen bond involving astatine

The effect of spin-orbit coupling (SOC) on the halogen bond involving astatine has been investigated using state-of-the-art two- and four-component relativistic calculations. Adducts between Cl-X (X = Cl, Br, I and At) and ammonia have been selected to establish a trend on going down the periodic table. The SOC influence has been explored not only on the geometric and energetic features that can be used to characterize the halogen bond strength but also on the three main contributions to it that are the charge transfer, the "σ-hole" (i.e. the localized region with a net positive electrostatic potential at the halogen site) and the "polar flattening" (which is related to the effective shape of the halogen site). A surprisingly large increase of the Cl-At dipole moment, due to the inclusion of SOC, has been worked out using four-component CCSD(T) reference calculations, indicating that this bond is significantly more ionic than one may predict. Due to the SOC effect, which induces a peculiar charge accumulation on the At side in the Cl-At dimer, a weakening of the astatine-mediated halogen bond occurs arising from the (i) reduced amount of charge transfer, (ii) decrease of the polar flattening and (iii) lowering of the short-range Coulomb potential. The analysis of the electronic structure of the Cl-At moiety allows for a rationalization of the SOC effects on all the considered features of the halogen bond, including an unprecedented unsymmetrical charge back-donation from Cl-At to ammonia.

Chemical bonding in cuprous complexes with simple nitriles: octet rule and resonance concepts versus quantitative charge-redistribution analysis

Resonance structures for six cuprous complexes with simple nitriles are interpreted by means of a quantitative analysis of charge redistribution upon copper-nitrile bonding.

How π back-donation quantitatively controls the CO stretching response in classical and non-classical metal carbonyl complexes

The CO stretching response upon coordination to a metal M to form [(L) n M(CO)] m complexes (L is an auxiliary ligand) is investigated in relation to the σ donation and π back-donation components of the M-CO bond and to the electrostatic effect exerted by the ligand-metal fragment. Our analysis encompasses over 30 carbonyls, in which the relative importance of donation, back-donation and electrostatics are varied either through the ligand in a series of [(L)Au(CO)]0/+ gold(i) complexes, or through the metal in a series of anionic, neutral and cationic homoleptic carbonyls. Charge-displacement analysis is used to obtain well-defined, consistent measures of σ donation and π back-donation charges, as well as to quantify the σ and π components of CO polarization. It is found that all complexes feature a comparable charge flow of σ symmetry (both in the M-CO bonding region and in the CO fragment itself), which is therefore largely uncorrelated to CO response. By contrast, π back-donation is exceptionally variable and is found to correlate tightly with the change in CO bond distance, with the shift in CO stretching frequency, and with the extent and direction (C → O or C ← O) of the CO π polarization. As a result, we conclusively show that π back-donation can be an important bond component also in non-classical carbonyls and we provide the framework in which the spectroscopic data on coordinated CO can be used to extract quantitative information on the π donor properties of metal-ligand moieties.