Pnicogen-Bonded Cyclic Trimers (PH2X)3 with X = F, Cl, OH, NC, CN, CH3, H, and BH2

Ab initio investigation of the competition of pnicogen, halogen, and hydrogen bonds resulting from the interactions between cyanophosphine and hypohalous acids

CONTEXT

The complexes formed as a result of the interactions between cyanophosphine (CP, H2PCN) and hypohalous acid molecules (HOX, X = F, Cl, Br, and I) were studied by employing ab initio computations conducted at the MP2/aug-cc-pVTZ level. Three types of complexes were acquired (I, II, and III) as a result of the (O∙∙∙P) pnicogen bond, the (N∙∙∙H) hydrogen bond, and the (N∙∙∙X) halogen bond interaction, respectively. The results of harmonic vibrational frequency calculations with no imaginary frequencies confirmed the structures as minima. In addition, given the interaction energy of the complexes, hydrogen bond complexes of structure II have the highest stability compared to other structures. In all studied complexes, the strength of the interactions depended on the electronegativity of the halogen atoms. The characteristics and nature of the whole three types of complexes were examined and evaluated with natural bond orbital (NBO), atom in molecules (AIM), molecular electrostatic potential (MEP) maps, non-covalent interaction (NCI) index, and electron density difference (EDD) analyses.

METHOD

The optimization of all complexes and corresponding monomers was conducted through the ab initio method, employing the MP2 level along with the aug/cc-pVTZ basis set for all atoms, except for the iodine (I) atom, for which the aug-cc-pVTZ (PP) basis set was employed. Subsequent frequency calculations were executed to ascertain the minimum energy state of the complexes at the MP2 level and the aug/cc-pVTZ basis set, utilizing Gaussian09 software. The MEP maps of the monomers were generated using the analysis-surface suite (WFA-SAS) software package. To probe the orbital interactions within the studied complexes, NBO analysis was performed employing NBO software. The assessment of bond nature, topological features, and electron density values at critical points for the studied complexes was undertaken using AIMAll software. The NCI index was derived utilizing Multiwfn software, and its three-dimensional representation was rendered using VMD software.

Unique Dispersion-Induced Tetrel Bond with Co-operative σ-hole-Induced Pnicogen Bond in the POCl3-Acetone Heterodimer: Experimental Confirmation at Low Temperatures

Both tetrel and pnicogen bonds are known to be induced through σ-/π-holes. This work reports computational and experimental evidence of the carbonyl carbon of acetone hosting a tetrel bond by dispersion rather electrostatic forces, for the first time, while phosphorus of POCl₃ sustains pnicogen bonding via the σ-hole. Heterodimers of POCl₃ with acetone (CH₃COCH₃) have been isolated within inert gas matrixes of Ar and N₂ at 12 K. Characteristic vibrational bands at P═O stretching of POCl₃ and C═O stretching of CH₃COCH₃ have been obtained in support of the computations. The potential energy surface has been traced computationally using ab initio and density functional methods. CH₃COCH₃ harboring such a tetrel bond, in itself, is quite intriguing. The interplay of these interactions has been comprehended by the quantum theory of atoms in molecules, natural bond orbital, energy decomposition, electrostatic potential mapping, and noncovalent interaction analyses.

Properties and Stabilities of Cyclic and Open Chains of Halogen Bonds

Open and cyclic chains from two to eight units of ICl and IF are constructed and examined by density functional theory (DFT) calculations. These chains contain either I···I or I···X halogen bonds (XBs) where X refers to Cl or F. The closed rings are more stable than the open chains due to the presence of an additional XB and enhanced cooperativity. This pattern is true even for most trimers where there is sizable geometric distortion in the rings. I···F rings are generally more stable than the corresponding I···I cycles as the I···F bond is stronger than I···I even in the simple dimer. However, I···I rings are comparable in energy to I···Cl. It is possible to construct I···I rings of at least as large as eight units, which are held together exclusively by XBs. On the other hand, the maximum possible size of I···X rings is 6. Red shifts are observed in the I-X stretching frequency bands, which magnify as the chain, both cyclic and open, grows longer. The NMR chemical shielding of the I atoms increases for I···I chains but diminishes when I···Cl bonds are present.

Plausible Pnicogen Bonding of epi-Cinchonidine as a Chiral Scaffold in Catalysis

As a non-covalent interaction of a chiral scaffold in catalysis, pnicogen bonding of epi-cinchonidine (epi-CD), a cinchona alkaloid, was simulated to consider whether the interaction can have the potential controlling enantiotopic face like hydrogen bonding. Among five reactive functional groups in epi-CD, two stable complexes of the hydroxyl group (X-epi-CD1) at C₁₇ and of the quinoline ring (X-epi-CD2) at N₁₆ with pnictide family analytes [X = substituted phosphine (PX), i.e., F, Br, Cl, CF₃, CN, HO, NO₂, and CH₃, and pnictide family analytes, i.e., PBr₃, BiI₃, SbI₃, and AsI₃] were predicted with intermolecular interaction energies, charge transfer (Q_Mulliken and Q_NBO), and band gap energies of HOMO-LUMO (Eg) at the B3LYP/6-31G(d,p) level of density functional theory. It was found that the dominant site of pnicogen bonding in epi-CD is the quinoline ring (N₁₆ atom) rather than the hydroxyl group (O₃₆ atom). In addition, the UV-Vis spectra of the complex were calculated by time-dependent density functional theory (TD-DFT) at the B3LYP/6-31+G(d,p) level and compared with experimental measurements. Through these calculations, two intermolecular interactions (H-bond vs. pnicogen bond) of epi-CD were compared.

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.

Not Only Hydrogen Bonds: Other Noncovalent Interactions

In this review, we provide a consistent description of noncovalent interactions, covering most groups of the Periodic Table. Different types of bonds are discussed using their trivial names. Moreover, the new name “Spodium bonds” is proposed for group 12 since noncovalent interactions involving this group of elements as electron acceptors have not yet been named. Excluding hydrogen bonds, the following noncovalent interactions will be discussed: alkali, alkaline earth, regium, spodium, triel, tetrel, pnictogen, chalcogen, halogen, and aerogen, which almost covers the Periodic Table entirely. Other interactions, such as orthogonal interactions and π-π stacking, will also be considered. Research and applications of σ-hole and π-hole interactions involving the p-block element is growing exponentially. The important applications include supramolecular chemistry, crystal engineering, catalysis, enzymatic chemistry molecular machines, membrane ion transport, etc. Despite the fact that this review is not intended to be comprehensive, a number of representative works for each type of interaction is provided. The possibility of modeling the dissociation energies of the complexes using different models (HSAB, ECW, Alkorta-Legon) was analyzed. Finally, the extension of Cahn-Ingold-Prelog priority rules to noncovalent is proposed.

A quantum chemical perspective on the potency of electron donors and acceptors in pnicogen bonds (AS...N, P...N, N...N)

A quantum chemical perspective of 31 structures contains electron acceptors: ASCl₃ (arsenic trichloride), PCl₃ (phosphorous trichloride) and NCl₃ (nitrogen trichloride); forming non-covalent bond with various nitrogen-based electron donors that resulted in pnicogen bonds, AS...N, P...N and N...N were calculated at M062X/def2-QZVP level of theory. Besides the above method, MP2/def2-QZVP and CCSD(T)/def2-QZVP level of theories have also been analysed to have in depth knowledge about the bonds formed. The nature of the bonds was assumed from the electrostatic potential evaluated for all the monomers, where σ hole is positive for all the monomers. The strongest pnicogen bonds are ASCl₃-NF₂H, PCl₃-NCH₃CH₃CH₃ and NCl₃-NCH₃CH₃CH₃ having interaction energies as -4.15, -11.58 and -3.25 kcal/mol, respectively, at MP2/def2-QZVP level of theory. Further at CCSD(T)/def2-QZVP level of theory, ASCl₃-NF₂H and NCl₃-NCH₃CH₃CH₃ are found to be the most stable with interaction energies as -3.53 and -2.45 kcal/mol, respectively. The potential energy surface scan was performed for all the stable complexes in order to confirm the existences of energies are true minima. Moreover to confirm the halogen and pnicogen bonds, AIM analysis was carried out. The results from the above factors of pnicogen bond will help crystal growth, material science and engineering community to explore novel materials, which abide for modernization. Graphical abstract PCl₃-NCH₃CH₃CH₃ complex with 2.61 Å and pnicogen angle of 178.54° is strong, and interaction energy is -11.58 kcal/mol. Electron donors - ASCl₃, PCl₃ and NCl₃ and electron acceptors -NCH₃CH₃CH₃, NH₃C₂ and NHCO have strong electrostatic contribution. High and low values of (ρ) ∇²(ρ) reveal the strong and weak pnicogen bond. Schematic representation of acceptors surrounded by its donors and Electrostatic Potential map.

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.

Cooperativity between the hydrogen bonding and σ-hole interaction in linear NCX···(NCH)n=2–5 and O3Z···(NCH)n=2–5 complexes (X = Cl, Br; Z = Ar, Kr): a comparative study

Quantum chemical calculations are performed to investigate the cooperativity of hydrogen bonding with halogen or aerogen bonding interactions in linear NCX···(NCH)n=2–5 and O3Z···(NCH)n=2–5 clusters, where X = Cl, Br and Z = Ar, Kr. To understand the cooperativity mechanism in these systems, the corresponding binary NCX···NCH and O3Z···NCH complexes are also considered. The binding distances, interaction energies, and bonding properties of the NCX···(NCH)n=2–5 and O3Z···(NCH)n=2–5 clusters are analyzed in detail. It is found that the cooperative effects in the hydrogen bonding tend to strengthen X···N and Z···N interactions. For both NCX···(NCH)n and O3Z···(NCH)n clusters, a small bond shrinkage is observed from n = 4 to n = 5, which suggests that the cooperativity effects are almost saturated in the larger clusters (n > 5). As the size of the X or Z atom is increased, the magnitude of the cooperative energy in these systems is also increased, which is mainly ascribed to changes in electrostatic potentials and orbital interactions. Our results indicate that the cooperative effects lead to a substantial change in the 14N nuclear quadrupole coupling constants of the NCH molecule.

Modulating of the pnicogen-bonding by a H⋯π interaction: An ab initio study

An ab initio study of the cooperativity in XH₂P⋯NCH⋯Z and XH₂P⋯CNH⋯Z complexes (X=F, Cl, Br, CN, NC; Z=C₂H₂,C₆H₆) connected by pnicogen-bonding and H⋯π interactions is carried out by means of MP2 computational method. A detailed analysis of the structures, interaction energies and bonding properties is performed on these systems. For each set of the complexes considered, a favorable cooperativity is observed, especially in X=F and CN complexes. However, for a given X or Z, the amount of cooperativity effects in XH₂P⋯CNH⋯Z complexes are more important than XH₂P⋯NCH⋯Z counterparts. Besides, the influence of a H⋯π interaction on a P⋯N (C) bond is more pronounced than that of a P⋯N (C) bond on a H⋯π bond. The quantum theory of atoms in molecules shows that ternary complexes have increased electron densities at their bond critical points relative to the corresponding binary systems. The results also indicate that the strength of the P⋯N(C) and H⋯π interactions increases in the presence of the solvent.

The mutual influence of Y···N and H···H interactions in XHY···NCH···HM complexes (X = F, Cl, Br; Y = S, Se; M = Li, Na, BeH, MgH): tuning of the chalcogen bond by dihydrogen bond interaction

The equilibrium structures, interaction energies, and bonding properties of ternary XHY···NCH···HM complexes are studied by ab initio calculations, where X = F, Cl, Br, Y = S, Se, and M = Li, Na, BeH, MgH. The ab initio calculations are carried out at the MP2/aug-cc-pVTZ level. The results indicate that all optimized Y···N and H···H binding distances in the ternary complexes are smaller than the corresponding values in the binary systems. The calculated cooperative energies (Ecoop) are between −0.20 kcal/mol in BrHS···NCH···HBeH and −3.29 kcal/mol in FHSe···NCH···HNa. For a given Y and M, the estimated Ecoop values increase as X = F > Cl > Br. In addition, the selenium-bonded complexes exibit larger Ecoop values than those of the sulfur-bonded counterparts. The cooperativity between Y···N and H···H interactions is further analyzed by quantum theory of atoms in molecules and natural bond orbital methods. Cooperative effects make an increase in the J(Y–N) and J(H–H) spin–spin coupling constants of the ternary complexes with respect to the binary systems.

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.

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.