Biphenyl: A stress tensor and vector‐based perspective explored within the quantum theory of atoms in molecules

We use quantum theory of atoms in molecules (QTAIM) and the stress tensor topological approaches to explain the effects of the torsion φ of the C‐C bond linking the two phenyl rings of the biphenyl molecule on a bond‐by‐bond basis using both a scalar and vector‐based analysis. Using the total local energy density H( r _b), we show the favorable conditions for the formation of the controversial H–H bonding interactions for a planar biphenyl geometry. This bond‐by‐bond QTAIM analysis is found to be agreement with an earlier alternative QTAIM atom‐by‐atom approach that indicated that the H–H bonding interaction provided a locally stabilizing effect that is overwhelmed by the destabilizing role of the C‐C bond. This leads to a global destabilization of the planar biphenyl conformation compared with the twisted global minimum. In addition, the H( r _b) analysis showed that only the central torsional C‐C bond indicated a minimum for a torsion φ value coinciding with that of the conventional global energy minimum. The H–H bonding interactions are found to be topologically unstable for any torsion of the central C‐C bond away from the planar biphenyl geometry. Conversely, we demonstrate that for 0.0° < φ < 39.95° there is a resultant increase in the topological stability of the C nuclei comprising the central torsional C‐C bond. Evidence is found of the effect of the H–H bonding interactions on the torsion φ of the central C‐C bond of the biphenyl molecule in the form of the QTAIM response β of the total electronic charge density ρ( r _b). Using a vector‐based treatment of QTAIM we confirm the presence of the sharing of chemical character between adjacent bonds. In addition, we present a QTAIM interpretation of hyperconjugation and conjugation effects, the former was quantified as larger in agreement with molecular orbital (MO) theory. The stress tensor and the QTAIM H atomic basin path set areas are independently found to be new tools relevant for the incommensurate gas to solid phase transition occurring in biphenyl for a value of the torsion reaction coordinate φ ≈ 5°. © 2015 Wiley Periodicals, Inc.     

Congested molecules. Where is the steric repulsion? An analysis of the electron density by the method of interacting quantum atoms

The computed electron density of several congested saturated hydrocarbons and halogenated derivatives has been analyzed by the method of interacting quantum atoms (IQA). For all the molecules studied, the calculations show the existence of a bond path between the congested atoms and which, according to the Quantum Theory of Atoms in Molecules, indicates that there is a stabilizing interaction between these atoms. The bond path is found to exist up to interatomic distances well‐beyond the sum of the van der Waals radii. The IQA results indicate that steric hindrance is not a repulsive force between the congested atoms but that is the result of an increase in the intra‐atomic or self‐energy of the congested atoms. This increase in self‐energy is caused by the deformation of the atomic basin of the congested atoms. © 2013 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.     

Hydrogen–hydrogen bonding in biphenyl revisited

The evidence for the stabilizing nature of the H–H bonding in planar biphenyl is succinctly reviewed. The stabilizing nature of the H–H bonding is revealed through a comparison of the atomic energy of every atom in planar biphenyl with the same atom in the twisted equilibrium structure. It is shown that the barrier to rotation via the planar transition state is the net resultant of a stabilisation of the four ortho-hydrogen atoms (by 8 kcal/mol each), a stabilisation of the two para-carbon atoms (by 3 kcal/mol each) and by the dominant destabilisation of the two carbon atoms joining the two rings—the two junction carbon atoms—(by 22 kcal/mol each). The energetic stabilisation of the four ortho-hydrogen atoms is further shown to be in large proportion due to the formation of the hydrogen–hydrogen interatomic surface. Furthermore, neither the “bond order” between the two junction carbon atoms nor the total electron delocalisation between the two rings exhibit a significant change in going from the planar to the twisted equilibrium geometry. These findings are in contrast with the classical view of a balance between “steric non-bonded repulsion” and better electron delocalisation as a function of the twist dihedral angle. Similar conclusions have been recently reached by Pacios and Gómez through a study of the electrostatic potential at the position of the hydrogen nuclei. We dedicate this article to Professor TM Krygowski on the occasion of his 70th birthday wishing him a long and productive life.     

Atomic Partitioning of the MPn (n = 2, 3, 4) Dynamic Electron Correlation Energy by the Interacting Quantum Atoms Method: A Fast and Accurate Electrostatic Potential Integral Approach

Recently, the quantum topological energy partitioning method called interacting quantum atoms (IQA) has been extended to MPn (n = 2, 3, 4) wave functions. This enables the extraction of chemical insight related to dynamic electron correlation. The large computational expense of the IQA-MPn approach is compensated by the advantages that IQA offers compared to older nontopological energy decomposition schemes. This expense is problematic in the construction of a machine learning training set to create kriging models for topological atoms. However, the algorithm presented here markedly accelerates the calculation of atomically partitioned electron correlation energies. Then again, the algorithm cannot calculate pairwise interatomic energies because it applies analytical integrals over whole space (rather than over atomic volumes). However, these pairwise energies are not needed in the quantum topological force field FFLUX, which only uses the energy of an atom interacting with all remaining atoms of the system that it is part of. Thus, it is now feasible to generate accurate and sizeable training sets at MPn level of theory. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.     

Hydrogen‐Bond Cooperative Effects in Small Cyclic Water Clusters as Revealed by the Interacting Quantum Atoms Approach

The cooperative effects of hydrogen bonding in small water clusters (H₂O)_n (n=3–6) have been studied by using the partition of the electronic energy in accordance with the interacting quantum atoms (IQA) approach. The IQA energy splitting is complemented by a topological analysis of the electron density (ρ( r )) compliant with the quantum theory of atoms‐in‐molecules (QTAIM) and the calculation of electrostatic interactions by using one‐ and two‐electron integrals, thereby avoiding convergence issues inherent to a multipolar expansion. The results show that the cooperative effects of hydrogen bonding in small water clusters arise from a compromise between: 1) the deformation energy (i.e., the energy necessary to modify the electron density and the configuration of the nuclei of the isolated water molecules to those within the water clusters), and 2) the interaction energy (E_int) of these contorted molecules in (H₂O)_n. Whereas the magnitude of both deformation and interaction energies is enhanced as water molecules are added to the system, the augmentation of the latter becomes dominant when the size of the cluster is increased. In addition, the electrostatic, classic, and exchange components of E_int for a pair of water molecules in the cluster (H₂O)_n?1 become more attractive when a new H₂O unit is incorporated to generate the system (H₂O)_n with the last‐mentioned contribution being consistently the most important part of E_int throughout the hydrogen bonds under consideration. This is opposed to the traditional view, which regards hydrogen bonding in water as an electrostatically driven interaction. Overall, the trends of the delocalization indices, δ(Ω,Ω′), the QTAIM atomic charges, the topology of ρ( r ), and the IQA results altogether show how polarization, charge transfer, electrostatics, and covalency contribute to the cooperative effects of hydrogen bonding in small water clusters. It is our hope that the analysis presented in this paper could offer insight into the different intra‐ and intermolecular interactions present in hydrogen‐bonded systems.     

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.     

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.     

Experimental and theoretical study of the intramolecular C–H···N and C–H···S hydrogen bonding effects in the 1H and 13C NMR spectra of the 2‐(alkylsulfanyl)‐5‐amino‐1‐vinylpyrroles: a particular state of amine nitrogen

In the ¹H NMR spectra of the 1‐vinylpyrroles with amino‐ and alkylsulfanyl groups in 5 and 2 positions, an extraordinarily large difference between resonance positions of the H_A and H_B terminal methylene protons of the vinyl group is discovered. Also, the one‐bond ¹J(C_β,H_B) coupling constant is surprisingly greater than the ¹J(C_β,H_A) coupling constant in pyrroles under investigation, while in all known cases, there was a reverse relationship between these coupling constants. These spectral anomalies are substantiated by quantum chemical calculations. The calculations show that the amine nitrogen lone pair is removed from the conjugation with the π‐system of the pyrrole ring so that it is directed toward the H_B hydrogen. These factors are favorable to the emergence of the intramolecular C–H_B???N hydrogen bonding in the s‐cis(N) conformation. On the other hand, the spatial proximity of the sulfur to the H_B hydrogen provides an opportunity of the intramolecular C–H_B???S hydrogen bonding in the s‐cis(S) conformation. Presence of the hydrogen bond critical points as well as ring critical point for corresponding chelate ring revealed by a quantum theory of atoms in molecules (QTAIM) approach confirms the existence of the weak intramolecular C–H???N and C–H???S hydrogen bonding. Therefore, an unusual high‐frequency shift of the H_B signal and the increase in the ¹J(C_β,H_B) coupling constant can be explained by the effects of hydrogen bonding. Copyright © 2013 John Wiley & Sons, Ltd.     

Reinvestigation of intramolecular hydrogen bond in malonaldehyde derivatives: An ab initio,AIM and NBO study

The RAHB systems in malonaldehyde and its derivatives at MP2/ 6‐311++G(d,p) level of theory were studied and their intramolecular hydrogen bond energies by using the related rotamers method was obtained. The topological properties of electron density distribution in O? H···O intramolecular hydrogen bond have been analyzed in term of quantum theory of atoms in molecules (QTAIM). Correlations between the H‐bond strength and topological parameters are probed. The results of QTAIM clearly showed that the linear correlation between the electron density distribution at HB critical point and RAHB ring critical point with the corresponding hydrogen bond energies was obtained. Moreover, it was found a linear correlation between the electronic potential energy density, V(r_cp), and hydrogen bond energy which can be used as a simple equation for evaluation of HB energy in complex RAHB systems. Finally, the similar linear treatment between the geometrical parameters, such as O···O or O? H distance, and Lp(O)→σ*_OH charge transfer energy with the intramolecular hydrogen bond energy is observed. © 2010 Wiley Periodicals, Inc., Int J Quantum Chem, 2011     

A Comparison of the Interacting Quantum Atoms (IQA) Analysis of the Two-Particle Density-Matrices of MP4SDQ and CCSD

Within the quantum topological energy partitioning method called Interacting Quantum Atoms (IQA) we transition from Møller-Plesset (MP4SDQ) to CCSD in calculating intra- and interatomic electron correlation energies for a set of hydrides, diatomics, a few simple molecules and non-covalently bonded complexes, using the uncontracted basis set 6-31++G(2d,2p). CCSD-IQA allows a more rigorous analysis of atomic electron correlation than that offered by Møller-Plesset, which returns IQA contributions that are identical to Hartree–Fock counterparts except for two-electron terms. The CCSD-IQA analysis returns bond and other interatomic correlation energies that are typically much larger in magnitude than the MP4SDQ values. Crisp patterns of energy transferability are detected in water clusters, both for intra-atomic and interatomic correlation energies. CCSD determines that the intra-atomic correlation energy of an oxygen drops by 15 kJ · mol^–1 for donating a hydrogen and by 25 kJ · mol^–1 for accepting a hydrogen.     

Assessment of intermolecular interactions at three sites of the arylalkyne in phenylacetylene‐containing lithium‐bonded complexes: Ab initio and QTAIM studies

The intermolecular interactions existing at three different sites between phenylacetylene and LiX (X = OH, NH₂, F, Cl, Br, CN, NC) have been investigated by means of second‐order Møller?Plesset perturbation theory (MP2) calculations and quantum theory of “atoms in molecules” (QTAIM) studies. At each site, the lithium‐bonding interactions with electron‐withdrawing groups (? F, ? Cl, ? Br, ? CN, ? NC) were found to be stronger than those with electron‐donating groups (? OH and ? NH₂). Molecular graphs of C₆H₅C?CH···LiF and πC₆H₅C?CH···LiF show the same connectional positions, and the electron densities at the lithium bond critical points (BCPs) of the πC₆H₅C?CH···LiF complexes are distinctly higher than those of the σC₆H₅C?CH···LiF complexes, indicating that the intermolecular interactions in the C₆H₅C?CH···LiX complexes can be mainly attributed to the π‐type interaction. QTAIM studies have shown that these lithium‐bond interactions display the characteristics of “closed‐shell” noncovalent interactions, and the molecular formation density difference indicates that electron transfer plays an important role in the formation of the lithium bond. For each site, linear relationships have been found between the topological properties at the BCP (the electron density ρ_b, its Laplacian ?²ρ_b, and the eigenvalue λ₃ of the Hessian matrix) and the lithium bond length d(Li‐bond). The shorter the lithium bond length d(Li‐bond), the larger ρ_b, and the stronger the π···Li bond. The shorter d(Li‐bond), the larger ?²ρ_b, and the greater the electrostatic character of the π···Li bond. © 2012 Wiley Periodicals, Inc.     

Multipolar electrostatics for proteins: Atom–atom electrostatic energies in crambin

Accurate electrostatics necessitates the use of multipole moments centered on nuclei or extra point charges centered away from the nuclei. Here, we follow the former alternative and investigate the convergence behavior of atom‐atom electrostatic interactions in the pilot protein crambin. Amino acids are cut out from a Protein Data Bank structure of crambin, as single amino acids, di, or tripeptides, and are then capped with a peptide bond at each side. The atoms in the amino acids are defined through Quantum Chemical Topology (QCT) as finite volume electron density fragments. Atom‐atom electrostatic energies are computed by means of a multipole expansion with regular spherical harmonics, up to a total interaction rank of L = ?_A+ ?_B + 1 = 10. The minimum internuclear distance in the convergent region of all the 15 possible types of atom‐atom interactions in crambin that were calculated based on single amino acids are close to the values calculated from di and tripeptides. Values obtained at B3LYP/aug‐cc‐pVTZ and MP2/aug‐cc‐pVTZ levels are only slightly larger than those calculated at HF/6‐31G(d,p) level. This convergence behavior is transferable to the well‐known amyloid beta polypeptide Aβ_1–42. Moreover, for a selected central atom, the influence of its neighbors on its multipole moments is investigated, and how far away this influence can be ignored is also determined. Finally, the convergence behavior of AMBER becomes closer to that of QCT with increasing internuclear distance. © 2013 Wiley Periodicals, Inc.     

Massenspektrometrische Untersuchung von Amiden X—Notiz zum Mechanismus der H-Abspaltung aus N,N-Dimethylbenzamid

The genesis of the [M? H]⁺ ion from N, N-dimethylbenzamides proceeds only via an interaction of the amido group with the hydrogens in the ortho-Position. The other hydrogen atoms are not involved in the fragmentation.     

A theoretical study of the intramolecular interaction between proximal atoms in planar conformations of biphenyl and related systems

The nature of the interaction between proximal hydrogens in planar biphenyl has been recently a matter of debate as arguments in favor of and against the existence of “H–H” bonding have been recently put forward. This issue is addressed here through the study of both the electron density ρ(r) and the electron localization function (ELF) η(r) obtained in quantum calculations on molecular systems with F atoms replacing hydrogens in the moiety that presents this interaction. The analysis of geometries and properties of ρ(r) and η(r) at both planar and twisted equilibrium conformations of those systems along with biphenyl, permits to get information on this intramolecular interaction that is compared with the use of the traditional chemical concepts (steric hindrance and π-resonance effects) involved. It is shown that although the ELF gives information compatible with these classical terms, this does not preclude the existence of bonds between proximal atoms with features rather similar to those of well-established intramolecular hydrogen bonds.     

Unusual structure of a biphenyl fragment: the important role of weak interactions

1,5‐Bis([1,1′‐biphenyl]‐4‐yl)‐5‐hydroxy‐3‐methylheptane‐1,6‐dione, C₃₂H₃₀O₃, was investigated by X‐ray analysis, quantum chemical calculations and AIM (atoms in molecules) analysis. It was shown that four intramolecular C—H…π forces are established between the two biphenyl fragments. C—H…π interactions lead to a bending of the biphenyl part. The energy of the intramolecular interactions was estimated and the exclusive role of the moment of force was established.     

Correlated one‐electron wave functions

Using four basis sets, 6‐311G(d,p), 6‐31+G(d,p), 6‐311++G(2d,2p), and 6‐311++G(3df,3pd), the optimized structures with all real frequencies were obtained at the MP2 level for dimers CH₂O? HF, CH₂O? H₂O, CH₂O? NH₃, and CH₂O? CH₄. The structures of CH₂O? HF, CH₂O? H₂O, and CH₂O? NH₃ are cycle‐shaped, which result from the larger bend of σ‐type hydrogen bonds. The bend of σ‐type H‐bond O…H? Y (Y?F, O, N) is illustrated and interpreted by an attractive interaction of a chemically intuitive π‐type hydrogen bond. The π‐type hydrogen bond is the interaction between one of the acidic H atoms of CH₂O and lone pair(s) on the F atom in HF, the O atom in H₂O, or the N atom in NH₃. By contrast with above the three dimers, for CH₂O? CH₄, because there is not a π‐type hydrogen‐bond to bend its linear hydrogen bond, the structure of CH₂O? CH₄ is a noncyclic shaped. The interaction energy of hydrogen bonds and the π‐type H‐bond are calculated and discussed at the CCSD(T)/6‐311++G(3df,3pd) level. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Influence of substitution on the strength and nature of CH···N hydrogen bond in XCCH···NH3 complexes

The effect of substitution on the strength and nature of CH···N hydrogen bond in XCCH···NH₃ (X = F, Cl, Br, OH, H, Me) and NCH···NH₃ complexes were investigated by quantum chemical calculations. Ab initio calculations were performed using MP2 method with a wide range of basis sets. With tacking into account the BSSE and ZPVE, the values of BEs decrease. Replacement of the nonparticipatory hydrogen atom of HCCH by the electronegative atoms (F, Cl, and Br), lead to the BEs increases. The BE corresponding to the replacement of the nonparticipatory hydrogen atom of HCCH by the OH and CH₃ groups decreases. A far greater enhancement of the interaction energy arises from replacement of HCCH by the more acidic HCN. The natural bond orbital analysis and the Bader's quantum theory of atoms in molecules were also used to elucidate the interaction characteristics of these complexes. The electrostatic nature of H‐bond interactions is predicted from QTAIM analysis. In addition, the relationship between the isotropic and anisotropic chemical shifts of the bridging hydrogen and binding energy of complexes as well as electron density at N···H BCPs were investigated. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011