Comparison of DAFH and FALDI-like approaches

Two complementary methodologies for extracting useful insights into electronic structure and bonding from contemporary wavefunctions are compared. The first of these, known as the analysis of domain-averaged Fermi holes (DAFH), mostly provides visually appealing descriptions of the role and the extent of electron sharing in chemical bonding. The second one, known as the fragment, atom, localized, delocalized and interatomic (FALDI) charge density decomposition scheme, uses the partitioning of certain localization and delocalization indices to focus on highly visual contributions associated with individual domains and with pairs of domains, respectively. Four variants of a FALDI-like approach are investigated here in some detail, mostly to establish which of them are the most reliable and the most informative. In addition to ‘full’ calculations that use the correlated pair density, the consequences for the DAFH and FALDI-like procedures of using instead a popular one-electron approximation are explored. Additionally, the geometry dependence of the degree of acceptability of the errors that this introduces for delocalization indices is assessed for different formal bond multiplicities. The familiar molecular test systems employed for these various linked investigations are the breaking of the bonds in H2 and in N2, as well as the nature of the bonding in B2H6, as a simple example of multicenter bonding. One of the key outcomes of this study is a clear understanding of how DAFH analysis and a particular variant of FALDI-like analysis could be most profitably deployed to extract complementary insights into more complex and/or controversial bonding situations.

Theoretical investigations of the chemical bonding in MM'O2 clusters (M, M' = Be, Mg, Ca)

Motivated by the known stability of the somewhat unusual Be₂O₂ rhombus, which features a short Be-Be distance but no direct metal-metal bonding, we investigate the nature of the bonding interactions in the analogous clusters MM'O₂ (M, M' = Be, Mg, Ca). CCSD/cc-pVTZ and CCSD(T)/cc-pVQZ calculations, amongst others, are used to determine optimized geometries and the dissociation energies for splitting the MM'O₂ clusters into metal oxide monomers. The primary tools used to investigate the chemical bonding are the analysis of domain-averaged Fermi holes, including the generation of localized natural orbitals, and the calculation of appropriate two- and three-center bond indices. Insights emerging from these various analyses concur with earlier studies of M₂O₂ rhombic clusters in that direct metal-metal bonding was not observed in the MM'O₂ rings whereas weak three-center (3c) bonding was detected in the MOM' moieties. In general terms, these mixed MM'O₂ clusters exhibit features that are intermediate between those of M₂O₂ and M'₂O₂, and the differences between the M and M' atoms appear to have little impact on the overall degree of 3c MOM' bonding. Graphical abstract Bonding situation in MM'O₂ clusters (M, M' = Be, Mg, Ca).

Domain overlap matrices from plane-wave-based methods of electronic structure calculation

Plane waves are one of the most popular and efficient basis sets for electronic structure calculations of solids; however, their delocalized nature makes it difficult to employ for them classical orbital-based methods of chemical bonding analysis. The quantum chemical topology approach, introducing chemical concepts via partitioning of real space into chemically meaningful domains, has no difficulties with plane-wave-based basis sets. Many popular tools employed within this approach, for instance delocalization indices, need overlap integrals over these domains-the elements of the so called domain overlap matrices. This article reports an efficient algorithm for evaluation of domain overlap matrix elements for plane-wave-based calculations as well as evaluation of its implementation for one of the most popular projector augmented wave (PAW) methods on the small set of simple and complex solids. The stability of the obtained results with respect to PAW calculation parameters has been investigated, and the comparison of the results with the results from other calculation methods has also been made.

The electron localization as the information content of the conditional pair density

In the present work, the information gained by an electron for "knowing" about the position of another electron with the same spin is calculated using the Kullback-Leibler divergence (DKL) between the same-spin conditional pair probability density and the marginal probability. DKL is proposed as an electron localization measurement, based on the observation that regions of the space with high information gain can be associated with strong correlated localized electrons. Taking into consideration the scaling of DKL with the number of σ-spin electrons of a system (N(σ)), the quantity χ = (N(σ) - 1) DKLfcut is introduced as a general descriptor that allows the quantification of the electron localization in the space. fcut is defined such that it goes smoothly to zero for negligible densities. χ is computed for a selection of atomic and molecular systems in order to test its capability to determine the region in space where electrons are localized. As a general conclusion, χ is able to explain the electron structure of molecules on the basis of chemical grounds with a high degree of success and to produce a clear differentiation of the localization of electrons that can be traced to the fluctuation in the average number of electrons in these regions.

Exploring the role of the 3-center–4-electron bond in hypervalent λ3-iodanes using the methodology of domain averaged Fermi holes

λ³-Iodanes are shown to express pronounced 3-center–4-electron bonds whose strength strongly correlates with the kind of ligands, thereby allowing to modulate the reactivity of these reagents.

Domain-averaged Fermi-hole analysis for solids

The domain-averaged Fermi hole (DAFH) orbitals provide highly visual representation of bonding in terms of orbital-like functions with attributed occupation numbers. It was successfully applied on many molecular systems including those with non-trivial bonding patterns. This article reports for the first time the extension of the DAFH analysis to the realm of extended periodic systems. Simple analytical model of DAFH orbital for single-band solids is introduced which allows to rationalize typical features that DAFH orbitals for extended systems may possess. In particular, a connection between Wannier and DAFH orbitals has been analyzed. The analysis of DAFH orbitals on the basis of DFT calculations is applied to hydrogen lattices of different dimensions as well as to the solids diamond, graphite, Na, Cu and NaCl. In case of hydrogen lattices, remarkable similarity is found between the DAFH orbitals evaluated with both the analytical approach and DFT. In case of the selected ionic and covalent solids the DAFH orbitals deliver bonding descriptions, which are compatible with classical orbital interpretation. For metals the DAFH analysis shows essential multicenter nature of bonding.

The rotational barrier in ethane: a molecular orbital study

The energy change on each Occupied Molecular Orbital as a function of rotation about the C-C bond in ethane was studied using the B3LYP, mPWB95 functional and MP2 methods with different basis sets. Also, the effect of the ZPE on rotational barrier was analyzed. We have found that σ and π energies contribution stabilize a staggered conformation. The σ(s) molecular orbital stabilizes the staggered conformation while the stabilizes the eclipsed conformation and destabilize the staggered conformation. The π(z) and molecular orbitals stabilize both the eclipsed and staggered conformations, which are destabilized by the π(v) and molecular orbitals. The results show that the method of calculation has the effect of changing the behavior of the energy change in each Occupied Molecular Orbital energy as a function of the angle of rotation about the C-C bond in ethane. Finally, we found that if the molecular orbital energy contribution is deleted from the rotational energy, an inversion in conformational preference occurs.

Exploring the origin of the internal rotational barrier for molecules with one rotatable dihedral angle

Continuing our recent endeavor, we systematically investigate in this work the origin of internal rotational barriers for small molecules using the new energy partition scheme proposed recently by one of the authors [S. B. Liu, J. Chem. Phys. 126, 244103 (2007)], where the total electronic energy is decomposed into three independent components, steric, electrostatic, and fermionic quantum. Specifically, we focus in this work on six carbon, nitrogen, and oxygen containing hydrides, CH(3)CH(3), CH(3)NH(2), CH(3)OH, NH(2)NH(2), NH(2)OH, and H(2)O(2), with only one rotatable dihedral angle [angle]H-X-Y-H (X,Y=C,N,O). The relative contributions of the different energy components to the total energy difference as a function of the internal dihedral rotation will be considered. Both optimized-geometry (adiabatic) and fixed-geometry (vertical) differences are examined, as are the results from the conventional energy partition and natural bond orbital analysis. A wealth of strong linear relationships among the total energy difference and energy component differences for different systems have been observed but no universal relationship applicable to all systems for both cases has been discovered, indicating that even for simple systems such as these, there exists no omnipresent, unique interpretation on the nature and origin of the internal rotation barrier. Different energy components can be employed for different systems in the rationalization of the barrier height. Confirming that the two differences, adiabatic and vertical, are disparate in nature, we find that for the vertical case there is a unique linear relationship applicable to all the six molecules between the total energy difference and the sum of the kinetic and electrostatic energy differences. For the adiabatic case, it is the total potential energy difference that has been found to correlate well with the total energy difference except for ethane whose rotation barrier is dominated by the quantum effect.