Properties of complexes H2C═(X)P:PXH2, for X = F, Cl, OH, CN, NC, CCH, H, CH3, and BH2: P···P pnicogen bonding at σ-holes and π-holes

Types of noncovalent bonds within complexes of thiazole with CF4 and SiF4

The five-membered heteroaromatic thiazole molecule contains a number of electron-rich regions that could attract an electrophile, namely the N and S lone pairs that lie in the molecular plane, and π-system areas above the plane. The possibility of each of these sites engaging in a tetrel bond (TB) with CF4 and SiF4, as well as geometries that encompass a CH⋯F H-bond, was explored via DFT calculations. There are a number of minima that occur in the pairing of thiazole with CF4 that are very close in energy, but these complexes are weakly bound by less than 2 kcal mol-1 and the presence of a true TB is questionable. The inclusion of zero-point vibrational energies alters the energetic ordering, which is further modified when entropic effects are added. The preferred geometry would thus be sensitive to the temperature of an experiment. Replacement of CF4 by SiF4 leaves intact most of the configurations, and their tight energetic clustering, the ordering of which is again altered as the temperature rises. But there is one exception in that by far the most tightly bound complex involves a strong Si⋯N TB between SiF4 and the lone pair of the thiazole N, with an interaction energy of 30 kcal mol-1. Even accounting for its high deformation energy and entropic considerations, this structure remains as clearly the most stable at any temperature.

Maximal occupation by bases of π-hole bands surrounding linear molecules

Linear molecules such as CO₂ contain a positive π-hole ring that surrounds C on the molecule's equator. Quantum calculations examine the question as to how many bases can simultaneously bind to this ring. Linear molecules examined are TO₂ , where T = C, Si, Ge, Sn; bases are NCH and NH₃ . CO₂ engages in the weakest of the tetrel bonds, and can bind up to three NCH and two NH₃ . Unlike σ-hole tetrel bonds, Si forms the strongest tetrel bonds, with interaction energies as high as 43 kcal/mol with NH₃ . But like GeO₂ , SiO₂ can sustain only two bases in its equatorial ring. The π-hole ring of SnO₂ can engage in up to four tetrel bonds with either NCH or NH₃ , even though these bonds are weaker than those with GeO₂ or SiO₂ . As all of these complexes cast TO₂ in the role of multiple electron acceptor, the resulting negative cooperativity makes each successive bond weaker than its predecessor as bases are added, as well as reducing the magnitude of the central molecule's π-hole.

Anatomy of π-hole bonds: Linear systems

The list of σ-hole bonds is long and growing, encompassing both H-bonds and its closely related halogen, chalcogen, etc., sisters. These bonds rely on the asymmetric distribution of electron density, whose depletion along the extension of a covalent bond leaves a positive region of electrostatic potential from which these bonds derive their name. However, the density distributions of other molecules contain analogous positive regions that lie out of the molecular plane known as π-holes, which are likewise capable of engaging in noncovalent bonds. Quantum calculations are applied to study such π-hole bonds that involve linear molecules, whose positive region is a circular belt surrounding the molecule, rather than the more restricted area of a σ-hole. These bonds are examined in terms of their most fundamental elements arising from the spatial dispositions of their relevant molecular orbitals and the π-holes in both the total electron density and the electrostatic potential to which they lead. Systems examined comprise tetrel, chalcogen, aerogen, and triel bonds, as well as those involving group II elements, with atoms drawn from various rows of the Periodic Table. The π-hole bonds established by linear molecules tend to be weaker than those of comparable planar systems.

Dissection of the Origin of π-Holes and the Noncovalent Bonds in Which They Engage

Accompanying the rapidly growing list of σ-hole bonds has come the acknowledgment of parallel sorts of noncovalent bonds which owe their stability in large part to a deficiency of electron density in the area above the molecular plane, known as a π-hole. The origins of these π-holes are probed for a wide series of molecules, comprising halogen, chalcogen, pnicogen, tetrel, aerogen, and spodium bonds. Much like in the case of their σ-hole counterparts, formation of the internal covalent π-bond in the Lewis acid molecule pulls density toward the bond midpoint and away from its extremities. This depletion of density above the central atom is amplified by an electron-withdrawing substituent. At the same time, the amplitude of the π*-orbital is enhanced in the region of the density-depleted π-hole, facilitating a better overlap with the nucleophile's lone pair orbital and a stabilizing n → π* charge transfer. The presence of lone pairs on the central atom acts to attenuate the π-hole and shift its position somewhat, resulting in an overall weakening of the π-hole bond. There is a tendency for π-hole bonds to include a higher fraction of induction energy than σ-bonds with proportionately smaller electrostatic and dispersion components, but this distinction is less a product of the σ- or π-character and more a function of the overall bond strength.

Study of Beryllium, Magnesium, and Spodium Bonds to Carbenes and Carbodiphosphoranes

The aim of this article is to present results of theoretical study on the properties of C⋯M bonds, where C is either a carbene or carbodiphosphorane carbon atom and M is an acidic center of MX2 (M = Be, Mg, Zn). Due to the rarity of theoretical data regarding the C⋯Zn bond (i.e., the zinc bond), the main focus is placed on comparing the characteristics of this interaction with C⋯Be (beryllium bond) and C⋯Mg (magnesium bond). For this purpose, theoretical studies (ωB97X-D/6-311++G(2df,2p)) have been performed for a large group of dimers formed by MX2 (X = H, F, Cl, Br, Me) and either a carbene ((NH2)2C, imidazol-2-ylidene, imidazolidin-2-ylidene, tetrahydropyrymid-2-ylidene, cyclopropenylidene) or carbodiphosphorane ((PH3)2C, (NH3)2C) molecule. The investigated dimers are characterized by a very strong charge transfer effect from either the carbene or carbodiphosphorane molecule to the MX2 one. This may even be over six times as strong as in the water dimer. According to the QTAIM and NCI method, the zinc bond is not very different than the beryllium bond, with both featuring a significant covalent contribution. However, the zinc bond should be definitely stronger if delocalization index is considered.

Noncovalent Bonds through Sigma and Pi-Hole Located on the Same Molecule. Guiding Principles and Comparisons

Over the last years, scientific interest in noncovalent interactions based on the presence of electron-depleted regions called σ-holes or π-holes has markedly accelerated. Their high directionality and strength, comparable to hydrogen bonds, has been documented in many fields of modern chemistry. The current review gathers and digests recent results concerning these bonds, with a focus on those systems where both σ and π-holes are present on the same molecule. The underlying principles guiding the bonding in both sorts of interactions are discussed, and the trends that emerge from recent work offer a guide as to how one might design systems that allow multiple noncovalent bonds to occur simultaneously, or that prefer one bond type over another.

Nitrogen as a pnicogen?: evidence for π-hole driven novel pnicogen bonding interactions in nitromethane-ammonia aggregates using matrix isolation infrared spectroscopy and ab initio computations

The role of nitrogen, the first member of the pnicogen group, as an electron donor in hypervalent non-covalent interactions has been established long ago, while observation of its electron accepting capability is still elusive experimentally, and remains quite intriguing, conceptually. In the light of minimal computational exploration of this novel class of pnicogen bonding so far, the present work provides experimental proof with unprecedented clarity, for the existence of N(acceptor)N(donor) interaction using the model nitromethane (NM) molecule with ammonia (AM) as a Lewis base in NM-AM aggregates. The NM-AM dimer, in which the nitrogen atom of NM (as a unique pnicogen) accepts electrons from AM (the traditional electron donor), was synthesized at low temperatures under isolated conditions within inert gas matrixes and was characterized using infrared spectroscopy. The experimental generation of the NM-AM dimer stabilized via NN interaction has strong corroboration from ab initio calculations. Furthermore, confirmation regarding the directional prevalence of this NN interaction over C-HN and N-HO hydrogen bonding is elucidated quantitatively by quantum theory of atoms in molecules (QTAIM), electrostatic potential mapping (ESP), natural bond orbital (NBO), non-covalent interaction (NCI) and energy decomposition (ED) analyses. The present study also allows the extension of σ-hole/π-hole driven interactions to the atoms of the second period, in spite of their low polarizability.

Unusual blue to red shifting of C-H stretching frequency of CHCl3 in co-operatively P⋯Cl phosphorus bonded POCl3-CHCl3 heterodimers at low temperature inert matrixes

Heterodimers of POCl₃-CHCl₃ were generated in Ne, Ar, and Kr matrixes at low temperatures and were studied using infrared spectroscopy. The remarkable role of co-operative pentavalent phosphorus bonding in the stabilization of the structure dictated by hydrogen bonding is deciphered. The complete potential energy surface of the heterodimer was scanned by ab initio and density functional theory computational methodologies. The hydrogen bond between the phosphoryl oxygen of POCl₃ and C-H group of CHCl₃ in heterodimers induces a blue-shift in the C-H stretching frequency within the Ne matrix. However, in Ar and Kr matrixes, the C-H stretching frequency is exceptionally red-shifted in stark contrast with Ne. The plausibility of the Fermi resonance by the C-H stretching vibrational mode with higher order modes in the heterodimers has been eliminated as a possible cause within Ar and Kr matrixes by isotopic substitution (CDCl₃) experiments. To evaluate the influence of matrixes as a possible cause of red-shift, self-consistent Iso-density polarized continuum reaction field model was applied. This conveyed the important role of the dielectric matrixes in inducing the fascinating vibrational shift from blue (Ne) to red (Ar and Kr) due to the matrix specific transmutation of the POCl₃-CHCl₃ structure. The heterodimer produced in the Ne matrix possesses a cyclic structure stabilized by hydrogen bonding with co-operative phosphorus bonding, while in Ar and Kr the generation of an acyclic open structure stabilized solely by hydrogen bonding is promoted. Compelling justification regarding the dispersion force based influence of matrix environments in addition to the well-known dielectric influence is presented.

Dominance of unique Pπ phosphorus bonding with π donors: evidence using matrix isolation infrared spectroscopy and computational methodology

Albeit the first account of hypervalentπ interactions has been reported with halogenπ interactions, the feasibility of their extension to other hypervalent atoms as possible Lewis acids is still open. In this work, the role of phosphorus as an acceptor from the π electron cloud (Pπ pnicogen or phosphorus bonding) in PCl3-C2H2 and PCl3-C2H4 heterodimers is explored, by combining matrix isolation infrared spectroscopy with ab initio and DFT computational methodologies. The respective potential energy surfaces of the PCl3-C2H2 and PCl3-C2H4 heterodimers reveal unique minima stabilized by a concert of reasonably strong to weak interactions, of which Pπ phosphorus bonding was energetically dominant. Heterodimers, trimers and tetramers bound primarily by this unique phosphorus bond were generated at low temperatures. The dominance of phosphorus bonding in the PCl3-C2H2 and PCl3-C2H4 heterodimers over other interactions (such as Hπ, HCl, HP, Clπ and lone pair-π interactions) was confirmed and substantiated using extended quantum theory of atoms in molecules, natural bond orbital, electrostatic potential mapping and energy decomposition analyses. The following inferences in correlation with results from non-covalent-interaction analysis offer a complete understanding of the nature of the Pπ phosphorus bonding interactions. The significance of electrostatic forces kinetically favoring the formation of phosphorus bonded heterodimers, in addition to thermodynamic stabilization, is demonstrated experimentally.

Complexes H2 CO:PXH2 and HCO2 H : PXH2 for X=NC, F, Cl, CN, OH, CCH, CH3 , and H: Pnicogen Bonds and Hydrogen Bonds

Ab initio MP2/aug'-cc-pVTZ calculations have been carried out to investigate H₂ CO : PXH₂ pnicogen-bonded complexes and HCO₂ H : PXH₂ complexes that are stabilized by pnicogen bonds and hydrogen bonds, with X=NC, F, Cl, CN, OH, CCH, CH₃ , and H. The binding energies of these complexes exhibit a second-order dependence on the O-P distance. DFT-SAPT binding energies correlate linearly with MP2 binding energies. The HCO₂ H : PXH₂ complexes are stabilized by both a pnicogen bond and a hydrogen bond, resulting in greater binding energies for the HCO₂ H : PXH₂ complexes compared to H₂ CO : PXH₂ . Neither the O-P distance across the pnicogen bond nor the O-P distance across the hydrogen bond correlates with the binding energies of these complexes. The nonlinearity of the hydrogen bonds suggests that they are relatively weak bonds, except for complexes in which the substituent X is either CH₃ or H. The pnicogen bond is the more important stabilizing interaction in the HCO₂ H : PXH₂ complexes except when the substituent X is a more electropositive group. EOM-CCSD spin-spin coupling constants ^1p J(O-P) across pnicogen bonds in H₂ CO:PXH₂ and HCO₂ H : PXH₂ complexes increase as the O-P distance decreases, and exhibit a second order dependence on that distance. There is no correlation between ^2h J(O-P) and the O-P distance across the hydrogen bond in the HCO₂ H : PXH₂ complexes. ^2h J(O-P) coupling constants for complexes with X=CH₃ and H have much greater absolute values than anticipated from their O-P distances.

On the ability of pnicogen atoms to engage in both σ and π-hole complexes. Heterodimers of ZF2C6H5 (Z = P, As, Sb, Bi) and NH3

When bound to a pair of F atoms and a phenyl ring, a pyramidal pnicogen (Z) atom can form a pnicogen bond wherein an NH₃ base lies opposite one F atom. In addition to this σ-hole complex, the ZF₂C₆H₅ molecule can distort in such a way that the NH₃ approaches on the opposite side to the lone pair on Z, where there is a so-called π-hole. The interaction energies of these π-hole dimers are roughly 30 kcal/mol, much larger than the equivalent quantities for the σ-hole complexes, which are only 4-13 kcal/mol. On the other hand, this large interaction energy is countered by the considerable deformation energy required for the Lewis acid to adopt the geometry necessary to form the π-hole complex. The overall energetics of the complexation reaction are thus more exothermic for the σ-hole dimers than for the π-hole dimers.

Elusive hypervalent phosphorus⋯π interactions: evidence for paradigm transformation from hydrogen to phosphorus bonding at low temperatures

A paradigm transformation from hydrogen to phosphorus bonding is found to depend on the proton affinity of the interacting π-systems.

Implications of monomer deformation for tetrel and pnicogen bonds

Monomer rearrangement raises the interaction energy by up to 20 kcal mol⁻¹and intensifies its σ-hole by a factor of 1.5–2.9.

Pentavalent phosphorus as a unique phosphorus donor in POCl3 homodimer and POCl3–H2O heterodimer: matrix isolation infrared spectroscopic and computational studies

Phosphorus, an important element among the pnicogen group, opens up avenues for experimental and computational explorations of its interaction in a variety of compounds.

The effect of anions on noncovalent interactions in model clusters of chalcogen-containing (CH3)2X (X = O, S, Se) molecules

A computational study of F⁻⋯(CH₃)₂O⋯CH₃F with F⁻ bound to the protons of the two methyl groups, found significant enhancement of the O⋯C interaction relative to the neutral (CH₃)₂O⋯CH₃F dyad.

Theoretical insights into the π-hole interactions in the complexes containing triphosphorus hydride (P3H3) and its derivatives

The π-hole of triphosphorus hydride (P₃H₃) and its derivatives Z₃X₃ (Z = P, As; X = H, F, Cl, Br) was discovered and analyzed. MP2/aug-cc-pVDZ calculations were performed on the π-hole interactions in the HCN...Z₃X₃ complexes and the mutual influence between π-hole interactions and the hydrogen bond in the HCN...HCN...Z₃X₃ and HCN...Z₃X₃...HCN complexes studied. The π-hole interaction belongs to the typical closed-shell noncovalent interaction. The linear relationship was found between the most positive electrostatic potential of the π-hole (V_S,max) and the interaction energy. Moreover, the V_S,max of the π-hole was also found to be linearly correlated to the electrostatic energy term, indicating the important contribution of the electrostatic energy term to the π-hole interaction. There is positive cooperativity between the π-hole interaction and the hydrogen bond in the termolecular complexes. The π-hole interaction has a greater influence on the hydrogen bond than vice versa. The mutual enhancing effect between the π-hole interaction and the hydrogen bond in the HCN...HCN...Z₃X₃ complexes is greater than that in the HCN...Z₃X₃...HCN complexes.

Non-covalent C-Cl…π interaction in acetylene-carbon tetrachloride adducts: Matrix isolation infrared and ab initio computational studies

Non-covalent halogen-bonding interactions between π cloud of acetylene (C2H2) and chlorine atom of carbon tetrachloride (CCl4) have been investigated using matrix isolation infrared spectroscopy and quantum chemical computations. The structure and the energies of the 1:1 C2H2-CCl4 adducts were computed at the B3LYP, MP2 and M05-2X levels of theory using 6-311++G(d,p) basis set. The computations indicated two minima for the 1:1 C2H2-CCl4 adducts; with the C-Cl…π adduct being the global minimum, where π cloud of C2H2 is the electron donor. The second minimum corresponded to a C-H…Cl adduct, in which C2H2 is the proton donor. The interaction energies for the adducts A and B were found to be nearly identical. Experimentally, both C-Cl…π and C-H…Cl adducts were generated in Ar and N2 matrixes and characterized using infrared spectroscopy. This is the first report on halogen bonded adduct, stabilized through C-Cl…π interaction being identified at low temperatures using matrix isolation infrared spectroscopy. Atoms in Molecules (AIM) and Natural Bond Orbital (NBO) analyses were performed to support the experimental results. The structures of 2:1 ((C2H2)2-CCl4) and 1:2 (C2H2-(CCl4)2) multimers and their identification in the low temperature matrixes were also discussed.

Modulating intramolecular P⋯N pnictogen interactions

The strength of P⋯N intramolecular pnictogen interactions can be modulated, enhanced or diminished upon substitution of different electron withdrawing or donor groups.

PCl3–C6H6 heterodimers: evidence for P⋯π phosphorus bonding at low temperatures

A phosphorous trichloride (PCl₃)–benzene (C₆H₆) heterodimer was generated in a low temperature N₂ matrix and was characterized using infrared spectroscopy.

Comparison of π-hole tetrel bonding with σ-hole halogen bonds in complexes of XCN (X = F, Cl, Br, I) and NH3

In addition to the standard halogen bond formed when NH₃approaches XCN (X = F, Cl, Br, I) along its molecular axis, a perpendicular approach is also possible, toward a π-hole that is present above the X–C bond.

NMR Investigations of Noncovalent Carbon Tetrel Bonds. Computational Assessment and Initial Experimental Observation

Group IV tetrel elements may act as tetrel bond donors, whereby a region of positive electrostatic potential (σ-hole) interacts with a Lewis base. The results of calculations of NMR parameters are reported for a series of model compounds exhibiting tetrel bonding from a methyl carbon to the oxygen or nitrogen atoms in various functional groups. The (13)C chemical shift (δiso) and the (1c)J((13)C,Y) coupling (Y = (17)O, (15)N) across the tetrel bond are recorded as a function of geometry. The sensitivity of the NMR parameters to the noncovalent interaction is demonstrated via an increase in δiso and in |(1c)J((13)C,Y)| as the tetrel bond shortens. Gauge-including projector-augmented wave density functional theory (DFT) calculations of δiso are reported for crystals that exhibit tetrel bonding in the solid state. Experimental δiso values for solid sarcosine and its tetrel-bonded salts corroborate the computational findings. This work offers new insights into tetrel bonding and facilitates the incorporation of tetrel bonds as restraints in NMR crystallographic structure refinement.

An ab initio study on tunability of σ-hole interactions in XHS:PH2Y and XH2P:SHY complexes (X = F, Cl, Br; Y = H, OH, OCH3, CH3, C2H5, and NH2)

Quantum chemical calculations are performed to investigate the tunability of σ-hole interactions in chalcogen-bonded XHS:PH2Y and pnicogen-bonded XH2P:SHY complexes, where X  =  F, Cl, Br and Y  =  H, OH, OCH3, CH3, C2H5, NH2. The formation of these binary complexes can be understood in terms of molecular electrostatic potentials (MEPs), considering the P and S atoms as an electron acceptor or an electron donor in the chalcogen and pnicogen bonds. The strength of the XHS:PH2Y and XH2P:SHY complexes for a given Y increases as follows: X = Br < Cl < F. In addition, an acceptable linear relationship is found between the interaction energies and the magnitudes of the product of most positive and negative MEPs. This finding along with the electron density difference maps provides a clear picture of the electrostatic nature of the interactions in the XHS:PH2Y and XH2P:SHY complexes. The calculated spin-spin coupling constants across the chalcogen bond interactions in the XHS:PH2Y complexes display a quadratic dependence with the P···S binding distance.

Synthesis, X-ray characterization and DFT studies of N-benzimidazolyl-pyrimidine–M(ii) complexes (M = Cu, Co and Ni): the prominent role of π-hole and anion–π interactions

Five Cu(ii), Co(ii) and Ni(ii) complexes with 2-(N-benzimidazolyl)-pyrimidine and nitrate co-ligands have been synthesized and characterized by X-ray diffraction.

Experimental observation and theoretical investigation of a novel Cd(ii) complex with π-hole interactions involving nitro groups

A cadmium(ii) complex, which represents a novel example of supramolecular system bearing a network of π-hole interactions, has been synthesized.

The dual role of pnicogen as Lewis acid and base and the unexpected interplay between the pnicogen bond and coordination interaction in H3N⋯FH2X⋯MCN (X = P and As; M = Cu, Ag, and Au)

Ternary systems H₃N⋯FH₂X⋯MCN (X = P and As; M = Cu, Ag, and Au) as well as the corresponding pnicogen-bonded and coordination-bonded binary systems have been studied.