Toward deformation densities for intramolecular interactions without radical reference states using the fragment,atom, localized,delocalized, and interatomic (FALDI) charge density decomposition scheme

A novel approach for calculating deformation densities is presented, which enables to calculate the deformation density resulting from a change between two chemical states, typically conformers, without the need for radical fragments. The Fragment, Atom, Localized, Delocalized, and Interatomic (FALDI) charge density decomposition scheme is introduced, which is applicable to static electron densities (FALDI‐ED), conformational deformation densities (FALDI‐DD) as well as orthodox fragment‐based deformation densities. The formation of an intramolecular NH⋅⋅⋅N interaction in protonated ethylene diamine is used as a case study where the FALDI‐based conformational deformation densities (with atomic or fragment resolution) are compared with an orthodox EDA‐based approach. Atomic and fragment deformation densities revealed in real‐space details that (i) pointed at the origin of density changes associated with the intramolecular H‐bond formation and (ii) fully support the IUPAC H‐bond representation. The FALDI scheme is equally applicable to intra‐ and intermolecular interactions. © 2017 Wiley Periodicals, Inc.     

Forced Bonding and QTAIM Deficiencies: A Case Study of the Nature of Interactions in He@Adamantane and the Origin of the High Metastability

Calculations within the framework of the interacting quantum atoms (IQA) approach have shown that the interactions of the helium atom with both tertiary, tC, and secondary, sC, carbon atoms in the metastable He@adamantane (He@adam) endohedral complex are bonding in nature, whereas the earlier study performed within the framework of Bader’s quantum theory of atoms in molecules (QTAIM) revealed that only He???tC interactions are bonding. The He???tC and He???sC bonding interactions are shown to be forced by the high pressure that the helium and carbon atoms exert upon each other in He@adam. The occurrence of a bonding interaction between the helium and sC atoms, which are not linked by a bond path, clearly shows that the lack of a bond path between two atoms does not necessarily indicate the lack of a bonding interaction, as is asserted by QTAIM. IQA calculations showed that not only the destabilization of the adamantane cage, but also a huge internal destabilization of the helium atom, contribute to the metastability of He@adam, these contributions being roughly equal. This result disproves previous opinions based on QTAIM analysis that only the destabilization of the adamantane cage accounts for the endothermicity of He@adam. Also, it was found that there is no homeomorphism of the ρ( r ) and ‐v( r ) fields of He@adam. Comparison of the IQA and QTAIM results on the interactions in He@adam exposes other deficiencies of the QTAIM approach. The reasons for the deficiencies in the QTAIM approach are analyzed.     

Hydrogen–hydrogen interaction in planar biphenyl: A theoretical study based on the interacting quantum atoms and Hirshfeld atomic energy partitioning methods

The nature of H‐H interaction between ortho‐hydrogen atoms in planar biphenyl is investigated by two different atomic energy partitioning methods, namely fractional occupation iterative Hirshfeld (FOHI) and interacting quantum atoms (IQA), and compared with the traditional virial‐based approach of quantum theory of atoms in molecules (QTAIM). In agreement with Bader's hypothesis of H? H bonding, partitioning the atomic energy into intra‐atomic and interatomic terms reveals that there is a net attractive interaction between the ortho‐hydrogens in the planar biphenyl. This falsifies the classical view of steric repulsion between the hydrogens. In addition, in contrast to the traditional QTAIM energy analysis, both FOHI and IQA show that the total atomic energy of the ortho‐hydrogens remains almost constant when they participate in the H‐H interaction. Although, the interatomic part of atomic energy of the hydrogens plays a stabilizing role during the formation of the H? H bond, it is almost compensated by the destabilizing effects of the intra‐atomic parts and consequently, the total energy of the hydrogens remains constant. The trends in the changes of intra‐atomic and interatomic energy terms of ortho‐hydrogens during H? H bond formation are very similar to those observed for the H₂ molecule. © 2014 Wiley Periodicals, Inc.     

A Bond Path and an Attractive Ehrenfest Force Do Not Necessarily Indicate Bonding Interactions: Case Study on M2X2 (M=Li,Na, K; X=H,OH, F,Cl)

Interactions in dimers of model alkali metal derivatives M₂X₂ (M=Li or Na or K; X=H or F, Cl, OH) are studied in the frame of the quantum theory of atoms in molecules (QTAIM) using the interacting quantum atoms approach (IQA). Contrary to opinion prevalent in QTAIM studies, the interaction between two anions linked by a bond path is demonstrated to be strongly repulsive. One may therefore say that a bond path does not necessarily indicate bonding interactions. The interactions between two anions or two cations that are not linked by a bond path are also strongly repulsive. The repulsive anion–anion and cation–cation interactions are outweighed by much stronger attractive anion–cation interactions, and the model molecules are therefore in a stable state. The attractive Ehrenfest forces (calculated in the frame of the QTAIM) acting across interatomic surfaces shared by anions in the dimers do not reflect the repulsive interactions between anions. Probable reasons of this disagreement are discussed. The force exerted on the nucleus and the electrons of a particular atom by the nucleus and the electrons of any another atom in question is proposed. It is assumed that this force unambiguously exposes whether basins of two atoms are attracted or repelled by each other in a polyatomic molecule.     

Intramolecular Hydrogen Bonds in Low‐Molecular‐Weight Polyethylene Glycol

We used static DFT calculations to analyze, in detail, the intramolecular hydrogen bonds formed in low‐molecular‐weight polyethylene glycol (PEG) with two to five repeat subunits. Both red‐shifted O?H???O and blue‐shifting C?H???O hydrogen bonds, which control the structural flexibility of PEG, were detected. To estimate the strength of these hydrogen bonds, the quantum theory of atoms in molecules was used. Car–Parrinello molecular dynamics simulations were used to mimic the structural rearrangements and hydrogen‐bond breaking/formation in the PEG molecule at 300 K. The time evolution of the H???O bond length and valence angles of the formed hydrogen bonds were fully analyzed. The characteristic hydrogen‐bonding patterns of low‐molecular‐weight PEG were described with an estimation of their lifetime. The theoretical results obtained, in particular the presence of weak C?H???O hydrogen bonds, could serve as an explanation of the PEG structural stability in the experimental investigation.     

Anion-Anion Complexes Established between Aspartate Dimers

Stable dimers aspartate-aspartate have been studied in aqueous and gas phase through theoretical simulations. The polarizable continuum model (PCM) has been applied to simulate the effect of the hydration on monomers and complexes. The quantum theory of atoms in molecules (QTAIM) and the interacting quantum atoms (IQA) scheme has been used to inquire into if, in the aqueous phase, individual hydrogen bonds have attractive electrostatic components. In all cases a spontaneous formation of the complexes in the aqueous phase are observed, while in the gas phase a considerable energy barrier must be overcome (between 100.8 to 263.2 kJ mol⁻¹). The intermolecular distance at which this barrier is indicates when the hydrogen-bond interactions begin to take importance between the dimers and the corresponding molecular recognition among them. The IQA analysis shows that in aqueous phase, the hydrogen bonds N−H⋅⋅⋅O are mainly electrostatic in nature with a certain covalent character which increases linearly with the decrease of internuclear distances H⋅⋅⋅O. The H⋅⋅⋅H interactions observed are stabilizing and they are mainly quantum in nature.     

Steric repulsions, rotation barriers, and stereoelectronic effects: a real space perspective

Widely used chemical concepts like Pauli repulsion or hyperconjugation, and their role in determining rotation barriers or stereoelectronic effects, are analyzed from the real space perspective of the interacting quantum atoms approach (IQA). IQA emerges from the quantum theory of atoms in molecules (QTAIM), but is free from the equilibrium geometry constraint of the former. A framework with both electronically unrelaxed and relaxed wavefunctions is presented that leads to an approximate correspondence between the IQA concepts and those used in the EDA (energy decomposition analysis) or NBO (natural bond orbital) procedures. We show that no net force acts upon the electrons in an electronically relaxed system, so that any reasonable definition of Pauli repulsion must involve unrelaxed state functions. Using antisymmetrized fragments clarifies that Pauli repulsions are energetically connected to the IQA deformation energies, leaving footprints in the finally relaxed states. Similarly, EDA or NBO hyperconjugative stabilizations are found to be naturally related to the IQA electron delocalization patterns. Applications to the rotation barrier of ethane and other simple systems are presented, and the very often forgotten role of electrostatic contributions in determining preferred conformations is highlighted.     

A bond path and an attractive Ehrenfest force do not necessarily indicate bonding interactions: case study on M2X2 (M = Li, Na, K; X = H, OH, F, Cl)

Interactions in dimers of model alkali metal derivatives M(2)X(2) (M = Li or Na or K; X = H or F, Cl, OH) are studied in the frame of the quantum theory of atoms in molecules (QTAIM) using the interacting quantum atoms approach (IQA). Contrary to opinion prevalent in QTAIM studies, the interaction between two anions linked by a bond path is demonstrated to be strongly repulsive. One may therefore say that a bond path does not necessarily indicate bonding interactions. The interactions between two anions or two cations that are not linked by a bond path are also strongly repulsive. The repulsive anion-anion and cation-cation interactions are outweighed by much stronger attractive anion-cation interactions, and the model molecules are therefore in a stable state. The attractive Ehrenfest forces (calculated in the frame of the QTAIM) acting across interatomic surfaces shared by anions in the dimers do not reflect the repulsive interactions between anions. Probable reasons of this disagreement are discussed. The force exerted on the nucleus and the electrons of a particular atom by the nucleus and the electrons of any another atom in question is proposed. It is assumed that this force unambiguously exposes whether basins of two atoms are attracted or repelled by each other in a polyatomic molecule.     

Do one-step mechanisms always involve simultaneous evolution of electron density? QTAIM/IQA analysis of the Curtius rearrangement

The Curtius rearrangement reaction is studied by using quantum theory of atoms in molecules (QTAIM) analysis of the electron density and the interacting quantum atoms (IQA) formalism. Although the rearrangements take place in one stage, two phases are distinguished when the rearranged atom is H: the first one corresponds to the separation of N₂, and the second one to the N-H/C-H bond rearrangement. The transition state (TS) for the reaction does not represent an intermediate between reagent and product for the migration but for the isolation of the N₂ molecule. When the migration is undergone by a fluorine atom, no electronic phases can be distinguished and the process is really concerted. As the migration happens closer to the TS, the TS is more similar to the product. The IQA analysis reveals different electron density evolutions for H and F migrations, and the scarce relevance (in terms of energy) of the point where BCPs appear or disappear.     

Performance of the density matrix functional theory in the quantum theory of atoms in molecules

The generalization to arbitrary molecular geometries of the energetic partitioning provided by the atomic virial theorem of the quantum theory of atoms in molecules (QTAIM) leads to an exact and chemically intuitive energy partitioning scheme, the interacting quantum atoms (IQA) approach, that depends on the availability of second-order reduced density matrices (2-RDMs). This work explores the performance of this approach in particular and of the QTAIM in general with approximate 2-RDMs obtained from the density matrix functional theory (DMFT), which rests on the natural expansion (natural orbitals and their corresponding occupation numbers) of the first-order reduced density matrix (1-RDM). A number of these functionals have been implemented in the promolden code and used to perform QTAIM and IQA analyses on several representative molecules and model chemical reactions. Total energies, covalent intra- and interbasin exchange-correlation interactions, as well as localization and delocalization indices have been determined with these functionals from 1-RDMs obtained at different levels of theory. Results are compared to the values computed from the exact 2-RDMs, whenever possible.     

Does Fluorine Participate in Halogen Bonding?

When R is sufficiently electron withdrawing, the fluorine in the R?F molecules could interact with electron donors (e.g., ammonia) and form a noncovalent bond (F ??? N). Although these interactions are usually categorized as halogen bonding, our studies show that there are fundamental differences between these interactions and halogen bonds. Although the anisotropic distribution of electronic charge around a halogen is responsible for halogen bond formations, the electronic charge around the fluorine in these molecules is spherical. According to source function analysis, F is the sink of electron density at the F ??? N BCP, whereas other halogens are the source. In contrast to halogen bonds, the F ??? N interactions cannot be regarded as lump–hole interactions; there is no hole in the valence shell charge concentration (VSCC) of fluorine. Although the quadruple moment of Cl and Br is mainly responsible for the existence of σ‐holes, it is negligibly small in the fluorine. Here, the atomic dipole moment of F plays a stabilizing role in the formation of F ??? N bonds. Interacting quantum atoms (IQA) analysis indicates that the interaction between halogen and nitrogen in the halogen bonds is attractive, whereas it is repulsive in the F ??? N interactions. Virial‐based atomic energies show that the fluorine, in contrast to Cl and Br, stabilize upon complex formation. According to these differences, it seems that the F ??? N interactions should be referred to as “fluorine bond” instead of halogen bond.     

Hydrogen bonding interactions between <Emphasis Type="Italic">N</Emphasis>,<Emphasis Type="Italic">N</Emphasis>-dimethylformamide and cysteine: DFT studies of structures,properties, and topologies

The hydrogen bonding interactions between cysteine and N,N-dimethylformamide (DMF) were studied at the extended hybrid functional DFT-X3LYP/6-311++G(d,p) level regarding their geometries, energies, vibrational frequencies, and topological features of the electron density. The quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analyses were employed to elucidate the interaction characteristics in the complexes. The results show that two intermolecular hydrogen bonds (H-bonds) are formed in one complex except few complexes with one intermolecular H-bond. The H-bonds involving O atom of DMF as H-bond acceptor usually are red-shifting H-bonds, while the blue-shifting H-bond usually involve methyl of DMF or methenyl of cysteine moiety as H-bond donors. Both hydrogen bonding interaction and structural deformation play important roles in the relative stabilities of the complexes. Due to the π-bond cooperativity, the strongest H-bond is formed between hydroxyl of cysteine moiety and O atom of DMF, however, the serious deformation counteract the hydrogen bonding interaction to a great extent. The complex involves a stronger hydrogen bonding interaction as well as the smaller deformation is the most stable one. The electron density (ρ_b) as well as its Laplacian (∇²ρ_b) at the H-bond critical point predicted by QTAIM is strongly correlated with the H-bond structural parameter (δR _H···Y) and the second-perturbation energies E(2) in the NBO scheme.     

Influence of the Li⋅⋅⋅π Interaction on the H/X⋅⋅⋅π Interactions in HOLi⋅⋅⋅C6H6⋅⋅⋅HOX/XOH (X=F,Cl, Br,I) Complexes

The influences of the Li???π interaction of C₆H₆???LiOH on the H???π interaction of C₆H₆???HOX (X=F, Cl, Br, I) and the X???π interaction of C₆H₆???XOH (X=Cl, Br, I) are investigated by means of full electronic second‐order Møller–Plesset perturbation theory calculations and “quantum theory of atoms in molecules” (QTAIM) studies. The binding energies, binding distances, infrared vibrational frequencies, and electron densities at the bond critical points (BCPs) of the hydrogen bonds and halogen bonds prove that the addition of the Li???π interaction to benzene weakens the H???π and X???π interactions. The influences of the Li???π interaction on H???π interactions are greater than those on X???π interactions; the influences of the H???π interactions on the Li???π interaction are greater than X???π interactions on Li???π interaction. The greater the influence of Li???π interaction on H/X???π interactions, the greater the influences of H/X???π interactions on Li???π interaction. QTAIM studies show that the intermolecular interactions of C₆H₆???HOX and C₆H₆???XOH are mainly of the π type. The electron densities at the BCPs of hydrogen bonds and halogen bonds decrease on going from bimolecular complexes to termolecular complexes, and the π‐electron densities at the BCPs show the same pattern. Natural bond orbital analyses show that the Li???π interaction reduces electron transfer from C₆H₆ to HOX and XOH.     

Density functional theory and topological analysis on the hydrogen bonding interactions in N-protonated adrenaline–DMSO complexes

In this study, complexes formed via hydrogen bond interactions between N-protonated adrenaline (AdH⁺) and DMSO have been studied by density functional theory (DFT). The relevant geometries, energies, and IR characteristics of the hydrogen bonds (H-bonds) have been systematically investigated. The natural bond orbital (NBO) and the quantum theory of atoms in molecule (QTAIM) analysis have also been applied to understand the nature of the hydrogen bonding interactions in complexes. The H-bonds involving amino or hydroxyls as H-donor are dominant H-bonds in complexes and are attributed to strong H-bonds. The weak H-bonds, such as π H-bonds and H-bonds involving methyl (DMSO) or methenyls (C2H6 and C5H7 of AdH⁺) as H-acceptors, were found in complexes as well. The complexes in which the dominant H-bond involves amino of AdH⁺ as H-donor are more stable than those with the dominant H-bond involving hydroxyls as H-donor. Some relationships between various properties of QTAIM, NBO, geometry as well as frequency were also investigated.     

Experimental and Theoretical Evidence for an Agostic Interaction in a Gold(III) Complex

The first agostic interaction in a gold complex is described. The presence of a bonding C?H???Au interaction in a cationic “tricoordinate” gold(III) complex was suggested by DFT calculations and was subsequently confirmed by NMR spectroscopy at low temperature. The agostic interaction was analyzed computationally using NBO and QTAIM analyses (NBO=natural bond orbital; QTAIM=quantum theory of atoms in molecules).     

Exact and exclusive electron localization indices within QTAIM atomic basins

Novel measures of electron (de)localization within the Quantum Theory of Atoms in Molecules (QTAIM) atomic basins are presented which, unlike orthodox localization indices (LIs), are fully exclusive and can be easily visualized. This work shows that QTAIM‐defined LIs describe a portion of interatomic delocalized electrons; hence, the chemical/physical interpretation of orthodox LIs is misleading. Using the recently introduced Fragment, Atomic, Localized, Delocalized, and Interatomic (FALDI) density decomposition technique we derive two novel sets of LIs and delocalization indices (DIs), by accounting for the overlap between localized and delocalized density functions. The FALDI‐based LIs and DIs perfectly recover chemically expected core and bonded electron count. Usefulness of new (de)localization indices and their 3D representations were demonstrated on a number of examples, including formamide and benzene. We therefore expect that the scheme reported in this work will provide a valuable stepping stone between classical conceptual chemistry and quantum chemical topology. © 2018 Wiley Periodicals, Inc.     

A Combined Experimental and Theoretical Electron Density Study of Intra‐ and Intermolecular Interactions in Thiourea S,S‐Dioxide

The thiourea S,S‐dioxide molecule is recognized as a zwitterion with a high dipole moment and an unusually long C? S bond. The molecule has a most interesting set of intermolecular interactions in the crystalline state—a relatively strong O???H? N hydrogen bond and very weak intermolecular C???S and N???O interactions. The molecule has C_s symmetry, and each oxygen atom is hydrogen‐bonded to two hydrogen atoms with O???H? N distances of 2.837 and 2.826 Å and angles of 176.61 and 158.38°. The electron density distribution is obtained both from Xray diffraction data at 110 K and from a periodic density functional theory (DFT) calculation. Bond characterization is made in terms of the analysis of topological properties. The covalent characters of the C? N, N? H, C? S, and S? O bonds are apparent, and the agreement on the topological properties between experiment and theory is adequate. The features of the Laplacian distributions, bond paths, and atomic domains are comparable. In a systematic approach, DFT calculations are performed based on a monomer, a dimer, a heptamer, and a crystal to see the effect on the electron density distribution due to the intermolecular interactions. The dipole moment of the molecule is enhanced in the solid state. The typical values of ρ_b and H_b of the hydrogen bonds and weak intermolecular C???S and N???O interactions are given. All the interactions are verified by the location of the bond critical point and its associated topological properties. The isovalue surface of Laplacian charge density and the detailed atomic graph around each atomic site reveal the shape of the valence‐shell charge concentration and provide a reasonable interpretation of the bonding of each atom.