Systematic Coupled Cluster Study of Noncovalent Interactions Involving Halogens, Chalcogens, and Pnicogens

Structures and stability of I2 and ICl complexes with pyridine: Ab initio and DFT study

Theoretical investigation of thermodynamic stability and bonding features of possible isomers of the molecular and ionic complexes of pyridine with molecular iodine and iodine monochloride IX (X = I,Cl) is presented. M06-2X DFT functional is found to provide bond distances and dissociation energies which are close to those obtained at high-level ab initio CCSD(T)/aug-cc-pvtz//CCSD/aug-cc-pvtz benchmark computations for the most stable isomers, formed via donation of a lone pair of nitrogen atom of pyridine to the iodine atom. These isomers are by 23-33 kJ mol-1 (in case of I2) and by 39-56 kJ mol-1 (in case of ICl) more stable than other molecular complexes. T-shaped π-σ* bonded isomers turn out to be energetically comparable with van der Waals bound compounds. Among the ionic isomers, structures featuring [IPy2]+ cation with I3 - or ICl2 - counterions are more stable. Oligomerization favors ionic isomers starting from the tetrameric clusters of the composition (IX)4Py4.

Mechanistic Insights into S-Depalmitolyse Activity of Cln5 Protein Linked to Neurodegeneration and Batten Disease: A QM/MM Study

Ceroid lipofuscinosis neuronal protein 5 (Cln5) is encoded by the CLN5 gene. The genetic variants of this gene are associated with the CLN5 form of Batten disease. Recently, the first crystal structure of Cln5 was reported. Cln5 shows cysteine palmitoyl thioesterase S-depalmitoylation activity, which was explored via fluorescent emission spectroscopy utilizing the fluorescent probe DDP-5. In this work, the mechanism of the reaction between Cln5 and DDP-5 was studied computationally by applying a QM/MM methodology at the ωB97X-D/6-31G(d,p):AMBER level. The results of our study clearly demonstrate the critical role of the catalytic triad Cys280-His166-Glu183 in S-depalmitoylation activity. This is evidenced through a comparison of the pathways catalyzed by the Cys280-His166-Glu183 triad and those with only Cys280 involved. The computed reaction barriers are in agreement with the catalytic efficiency. The calculated Gibb's free-energy profile suggests that S-depalmitoylation is a rate-limiting step compared to the preceding S-palmitoylation, with barriers of 26.1 and 25.3 kcal/mol, respectively. The energetics were complemented by monitoring the fluctuations in the electron density distribution through NBO charges and bond strength alterations via local mode stretching force constants during the catalytic pathways. This comprehensive protocol led to a more holistic picture of the reaction mechanism at the atomic level. It forms the foundation for future studies on the effects of gene mutations on both the S-palmitoylation and S-depalmitoylation steps, providing valuable data for the further development of enzyme replacement therapy, which is currently the only FDA-approved therapy for childhood neurodegenerative diseases, including Batten disease.

Electron delocalization in defect-containing graphene and its influence on tetrel bond formation

Using the advanced analyses of electron density and fermionic potential, we show how electron delocalization influences the ability of defect-containing graphene to form tetrel bonds. The Cg atoms of a vacancy defect can produce one nonpolar interaction, alongside a peculiar polar Cg⋯Cg bond. The latter stems from the presence of a localized electron pair on a vacancy defect Cg atom and the local depletion of electron localization on another Cg atom. This interaction is an example of intralayer tetrel bond. In the presence of an absorbed molecule of bisphenol A diglycidyl ether (DGEBA), graphene is able to form incipient tetrel Cg⋯O bonds with an ether group oxygen. In contrast to an epoxy group oxygen, the disposition of the ether oxygen often causes the orientation of electron-rich π-domains of graphene carbon on the weakly expressed electrophilic region of the oxygen. In the case of graphene with a point Si defect, the Si atom can form quite strong Si⋯C interactions with the DGEBA aryl carbons. In contrast to other noncovalent bonds, this interaction significantly alters the electron (de)localization on the Si atom and in the aryl ring. The reliability of the obtained results is enhanced by the use of multiple 2D periodic models with defects located at different positions along the DGEBA skeleton.

Empowering non-covalent hydrogen, halogen, and [S-N]2 bonds in synergistic molecular assemblies on Au(111)

Non-covalent bonds are fundamental for designing self-assembled organic structures with potentially high responsiveness to mechanical, light, and thermal stimuli. The weak intermolecular interaction allows triggering charge-transport, energy-conversion, enzymatic, and catalytic activity, to name a few. Here, we discuss the synergistic action that multiple highly-directional and purely electrostatic bonds have in assembling one molecular specie, namely 4,7-dibromobenzo[c]-1,2,5-thiadiazole (2Br-BTD), in two different patterns on the Au(111) surface. We find, using scanning tunneling microscopy (STM) and density functional theory (DFT), that multiple secondary-interactions strengthen the electrostatic attraction between the pnicogen and chalcogen atoms forming [S-N]₂ heterocycles, the building block of the two networks. Among these interactions, there are halogen-halogen bonds that form characteristic supra-molecular synthons of 3, 4, or 6 molecules. However, not all these nodal structures contribute to the cohesion of the system. In such cases, other secondary bonds involving hydrogen or nitrogen compensate for the eventual deficiency.

A Theoretical Study of the Halogen Bond between Heteronuclear Halogen and Benzene

Halogen bonds play an important role in many fields, such as biological systems, drug design and crystal engineering. In this work, the structural characteristics of the halogen bond between heteronuclear halogen XD (ClF, BrCl, IBr, ICl, BrF and IF) and benzene were studied using density functional theory. The structures of the complexes between heteronuclear halogen and benzene have Cs symmetry. The interaction energies of the complexes between heteronuclear halogen XD (ClF, BrCl, IBr, ICl, BrF and IF) and benzene range from -27.80 to -37.18 kJ/mol, increasing with the increases in the polarity between the atoms of X and D, and are proportional to the angles of a between the Z axis and the covalent bond of heteronuclear halogen. The electron density (ρ) and corresponding Laplacian (∇²ρ) values indicate that the interaction of the heteronuclear halogen and benzene is a typical long-range weak interaction similar to a hydrogen bond. Independent gradient model analysis suggests that the van der Waals is the main interaction between the complexes of heteronuclear halogen and benzene. Symmetry-adapted perturbation theory analysis suggests that the electrostatic interaction is the dominant part in the complexes of C₆H₆⋯ClF, C₆H₆⋯ICl, C₆H₆⋯BrF and C₆H₆⋯IF, and the dispersion interaction is the main part in the complexes of C₆H₆⋯BrCl, C₆H₆⋯IBr.

Papaya Leaf Extracts as Potential Dengue Treatment: An In-Silico Study

Dengue fever (DF), dengue hemorrhagic fever (DHF), and dengue shock syndrome (DSS) cause serious public health problems, with nearly 390 million people affected and 20,000 deaths per year in tropical and subtropical countries. Despite numerous attempts, no antiviral drug or vaccine is currently available to combat the manifestation. The challenge of discovering an efficient vaccine is enhanced by the surplus presence of efficient vectors and drug resistance from the virus. For centuries, papaya (Carica papaya) extracts have been traditionally used to treat DF, DHF, and DSS. In the present study, we systematically investigated seven compounds isolated from papaya leaf extract with regard to their potential as inhibitors for non-structural (NS) proteins, NS3 and NS5, which play a crucial role in viral RNA replication. The computational tools applied stretched across classical molecular docking, molecular dynamics (MD) simulations and SwissADME used to calculate binding affinities; binding free energies; Absorption, Distribution, Metabolism, and Excretion (ADME); and drug-likeness properties, thus, identifying Kaempferol, Chlorogenic acid, and Quercetin as potential candidates, with Kaempferol and Quercetin scoring best. Therefore, for the Kaempferol and Quercetin complexes, hybrid quantum mechanical/molecular mechanical (QM/MM) geometry and frequency calculations were performed, followed by the local mode analysis developed in our group to quantify Kaempferol-NS and Quercetin-NS hydrogen bonding. Given the non-toxic nature and the wide availability of the Kaempferol and Quercetin papaya extract in almost all of the susceptible regions, and our results showing high NS3 and NS5 binding affinities and energies, strong hydrogen bonding with both NS3 and NS5, and excellent ADME properties, we suggest Kaempferol and Quercetin as a strong NS3 and NS5 inhibitor to be further investigated in vitro.

Involvement of Arsenic Atom of AsF3 in Five Pnicogen Bonds: Differences between X-ray Structure and Theoretical Models

Bonding within the AsF₃ crystal is analyzed via quantum chemical methods so as to identify and quantify the pnicogen bonds that are present. The structure of a finite crystal segment containing nine molecules is compared with that of a fully optimized cluster of the same size. The geometries are qualitatively different, with a much larger binding energy within the optimized nonamer. Although the total interaction energy of a central unit with the remaining peripheral molecules is comparable for the two structures, the binding of the peripherals with one another is far larger in the optimized cluster. This distinction of much stronger total binding within the optimized cluster is not limited to the nonamer but repeats itself for smaller aggregates as well. The average binding energy of the cluster rises quickly with size, asymptotically approaching a value nearly triple that of the dimer.

Definition of the Pnictogen Bond: A Perspective

This article proposes a definition for the term “pnictogen bond” and lists its donors, acceptors, and characteristic features. These may be invoked to identify this specific subset of the inter- and intramolecular interactions formed by elements of Group 15 which possess an electrophilic site in a molecular entity.

Mechanistic Details of the Sharpless Epoxidation of Allylic Alcohols—A Combined URVA and Local Mode Study

In this work, we investigated the catalytic effects of a Sharpless dimeric titanium (IV)–tartrate–diester catalyst on the epoxidation of allylalcohol with methyl–hydroperoxide considering four different orientations of the reacting species coordinated at the titanium atom (reactions R1–R4) as well as a model for the non-catalyzed reaction (reaction R0). As major analysis tools, we applied the URVA (Unified Reaction Valley Approach) and LMA (Local Mode Analysis), both being based on vibrational spectroscopy and complemented by a QTAIM analysis of the electron density calculated at the DFT level of theory. The energetics of each reaction were recalculated at the DLPNO-CCSD(T) level of theory. The URVA curvature profiles identified the important chemical events of all five reactions as peroxide OO bond cleavage taking place before the TS (i.e., accounting for the energy barrier) and epoxide CO bond formation together with rehybridization of the carbon atoms of the targeted CC double bond after the TS. The energy decomposition into reaction phase contribution phases showed that the major effect of the catalyst is the weakening of the OO bond to be broken and replacement of OH bond breakage in the non-catalyzed reaction by an energetically more favorable TiO bond breakage. LMA performed at all stationary points rounded up the investigation (i) quantifying OO bond weakening of the oxidizing peroxide upon coordination at the metal atom, (ii) showing that a more synchronous formation of the new CO epoxide bonds correlates with smaller bond strength differences between these bonds, and (iii) elucidating the different roles of the three TiO bonds formed between catalyst and reactants and their interplay as orchestrated by the Sharpless catalyst. We hope that this article will inspire the computational community to use URVA complemented with LMA in the future as an efficient mechanistic tool for the optimization and fine-tuning of current Sharpless catalysts and for the design new of catalysts for epoxidation reactions.

Non-covalent interactions atlas benchmark data sets 4: σ-hole interactions

The SH250×10 dataset presented here extends the Non-Covalent Interactions Atlas database (https://www.nciatlas.org) to complexes bound by σ-hole interactions - halogen, chalcogen and pnictogen bonds. It comprises 250 complexes where Cl, Br, I, S, Se, P and As interact with diverse electron donors. An accurate CCSD(T)/CBS benchmark is provided for ten points along a dissociation curve of each complex. The SH250×10 set is used in testing a wide variety of DFT functionals and semiempirical quantum-mechanical methods. In DFT calculations, the new data set exposes large errors of some functionals related to exaggerated charge transfer. The size and diversity of the data set have also been exploited in the reparametrization of a halogen-bond correction for the PM6 semiempirical method.

Bonding in Nitrile Photo-dissociating Ruthenium Drug Candidates --A Local Vibrational Mode Study

In this work, we investigated bonding features 15 ruthenium complexes of the type [Ru(tpy)(L)-(CH3CN)]n+, containing the tridentate tpy ligand (tpy = 2,2':6',2'--terpyridine) and various bidentate ancillary ligands, 12 compounds originally synthesized by Loftus et al. (J. Phys. Chem. C 123, 10291-10299 (2019)) complemented with three additional complexes. The main focus of our work was to relate these local features to the experimental data of Loftus et al. which assess the efficiency of nitrile release in an indirect way via observed quantum yields for ruthenium water association after nitrile release. As a tool to quantitatively assess Ru-NC and Ru-L bonding we utilized the local vibrational mode analysis complemented by the topological analysis of the electron density and the natural bond orbital analysis. Interestingly, the stronger Ru-NC bonds have the greater observed quantum yields, leading to the conclusion that the observed quantum yields are a result of a complex interplay of several processes excluding a direct relationship between QY and Ru-NC or Ru-L bond strengths. We identified the ST splitting as one of the key players and not the Ru-NC bond strength, as one may have thought. In summary, this work has presented a modern computational tool set for the investigation of bonding features applied to nitrile photo-dissociating ruthenium drug candidates forming a valuable basis for future design and fine tuning of nitrile releasing ruthenium compounds, as well as for the understanding of how local properties affect overall experimental outcomes.

Uranium: The Nuclear Fuel Cycle and Beyond

This review summarizes the recent developments regarding the use of uranium as nuclear fuel, including recycling and health aspects, elucidated from a chemical point of view, i.e., emphasizing the rich uranium coordination chemistry, which has also raised interest in using uranium compounds in synthesis and catalysis. A number of novel uranium coordination features are addressed, such the emerging number of U(II) complexes and uranium nitride complexes as a promising class of materials for more efficient and safer nuclear fuels. The current discussion about uranium triple bonds is addressed by quantum chemical investigations using local vibrational mode force constants as quantitative bond strength descriptors based on vibrational spectroscopy. The local mode analysis of selected uranium nitrides, N≡U≡N, U≡N, N≡U=NH and N≡U=O, could confirm and quantify, for the first time, that these molecules exhibit a UN triple bond as hypothesized in the literature. We hope that this review will inspire the community interested in uranium chemistry and will serve as an incubator for fruitful collaborations between theory and experimentation in exploring the wealth of uranium chemistry.

BF3–Catalyzed Diels–Alder Reaction between Butadiene and Methyl Acrylate in Aqueous Solution—An URVA and Local Vibrational Mode Study

In this study we investigate the Diels–Alder reaction between methyl acrylate and butadiene, which is catalyzed by BF3 Lewis acid in explicit water solution, using URVA and Local Mode Analysis as major tools complemented with NBO, electron density and ring puckering analyses. We considered four different starting orientations of methyl acrylate and butadiene, which led to 16 DA reactions in total. In order to isolate the catalytic effects of the BF3 catalyst and those of the water environment and exploring how these effects are synchronized, we systematically compared the non-catalyzed reaction in gas phase and aqueous solution with the catalyzed reaction in gas phase and aqueous solution. Gas phase studies were performed at the B3LYP/6-311+G(2d,p) level of theory and studies in aqueous solution were performed utilizing a QM/MM approach at the B3LYP/6-311+G(2d,p)/AMBER level of theory. The URVA results revealed reaction path curvature profiles with an overall similar pattern for all 16 reactions showing the same sequence of CC single bond formation for all of them. In contrast to the parent DA reaction with symmetric substrates causing a synchronous bond formation process, here, first the new CC single bond on the CH2 side of methyl acrylate is formed followed by the CC bond at the ester side. As for the parent DA reaction, both bond formation events occur after the TS, i.e., they do not contribute to the energy barrier. What determines the barrier is the preparation process for CC bond formation, including the approach diene and dienophile, CC bond length changes and, in particular, rehybridization of the carbon atoms involved in the formation of the cyclohexene ring. This process is modified by both the BF3 catalyst and the water environment, where both work in a hand-in-hand fashion leading to the lowest energy barrier of 9.06 kcal/mol found for the catalyzed reaction R1 in aqueous solution compared to the highest energy barrier of 20.68 kcal/mol found for the non-catalyzed reaction R1 in the gas phase. The major effect of the BF3 catalyst is the increased mutual polarization and the increased charge transfer between methyl acrylate and butadiene, facilitating the approach of diene and dienophile and the pyramidalization of the CC atoms involved in the ring formation, which leads to a lowering of the activation energy. The catalytic effect of water solution is threefold. The polar environment leads also to increased polarization and charge transfer between the reacting species, similar as in the case of the BF3 catalyst, although to a smaller extend. More important is the formation of hydrogen bonds with the reaction complex, which are stronger for the TS than for the reactant, thus stabilizing the TS which leads to a further reduction of the activation energy. As shown by the ring puckering analysis, the third effect of water is space confinement of the reacting partners, conserving the boat form of the six-member ring from the entrance to the exit reaction channel. In summary, URVA combined with LMA has led to a clearer picture on how both BF3 catalyst and aqueous environment in a synchronized effort lower the reaction barrier. These new insights will serve to further fine-tune the DA reaction of methyl acrylate and butadiene and DA reactions in general.

Capturing Individual Hydrogen Bond Strengths in Ices via Periodic Local Vibrational Mode Theory: Beyond the Lattice Energy Picture

Local stretching force constants derived from periodic local vibrational modes at the vdW-DF2 density functional level have been employed to quantify the intrinsic hydrogen bond strength of 16 ice polymorphs, ices Ih, II, III, IV, V, VI, VII, VIII, IX, XI, XII, XIII, XIV, XV, XVII, and XIX, that are stable under ambient to elevated pressures. Based on this characterization on 1820 hydrogen bonds, relationships between local stretching force constants and structural parameters such as hydrogen bond length and angle were identified. Moreover, different bond strength distributions, from uniform to inhomogeneous, were observed for the 16 ices and could be explained in relation to different local structural elements within ices, that is, rings, that consist of different hydrogen bond types. In addition, criteria for the classification of hydrogen bonds as strong, intermediate, and weak were introduced. The latter was used to explore a different dimension of the water-ice phase diagram. These findings will provide important guidelines for assessing the credibility of new ice structures.

Triel bonds within anionanion complexes

The ability of two anions to interact with one another is tested in the context of pairs of TrX₄^- homodimers, where Tr represents any of the triel atoms B, Al, Ga, In, or Tl, and X refers to a halogen substituent F, Cl, or Br. None of these pairs engage in a stable complex in the gas phase, but the situation reverses in water where the two monomers are held together by Tr⋯X triel bonds, complemented by stabilizing interactions between X atoms. Some of these bonds are quite strong, notably those involving TrF₄^-, with interaction energies surpassing 30 kcal mol^-1. Others are very much weaker, with scarcely exothermic binding energies. The highly repulsive electrostatic interactions are counteracted by large polarization energies.

Hydrogen Bonding in Natural and Unnatural Base Pairs-A Local Vibrational Mode Study

In this work hydrogen bonding in a diverse set of 36 unnatural and the three natural Watson Crick base pairs adenine (A)-thymine (T), adenine (A)-uracil (U) and guanine (G)-cytosine (C) was assessed utilizing local vibrational force constants derived from the local mode analysis, originally introduced by Konkoli and Cremer as a unique bond strength measure based on vibrational spectroscopy. The local mode analysis was complemented by the topological analysis of the electronic density and the natural bond orbital analysis. The most interesting findings of our study are that (i) hydrogen bonding in Watson Crick base pairs is not exceptionally strong and (ii) the N-H⋯N is the most favorable hydrogen bond in both unnatural and natural base pairs while O-H⋯N/O bonds are the less favorable in unnatural base pairs and not found at all in natural base pairs. In addition, the important role of non-classical C-H⋯N/O bonds for the stabilization of base pairs was revealed, especially the role of C-H⋯O bonds in Watson Crick base pairs. Hydrogen bonding in Watson Crick base pairs modeled in the DNA via a QM/MM approach showed that the DNA environment increases the strength of the central N-H⋯N bond and the C-H⋯O bonds, and at the same time decreases the strength of the N-H⋯O bond. However, the general trends observed in the gas phase calculations remain unchanged. The new methodology presented and tested in this work provides the bioengineering community with an efficient design tool to assess and predict the type and strength of hydrogen bonding in artificial base pairs.

Halogen Bonding Involving I2 and d8 Transition-Metal Pincer Complexes

We systematically investigated iodine–metal and iodine–iodine bonding in van Koten’s pincer complex and 19 modifications changing substituents and/or the transition metal with a PBE0–D3(BJ)/aug–cc–pVTZ/PP(M,I) model chemistry. As a novel tool for the quantitative assessment of the iodine–metal and iodine–iodine bond strength in these complexes we used the local mode analysis, originally introduced by Konkoli and Cremer, complemented with NBO and Bader’s QTAIM analyses. Our study reveals the major electronic effects in the catalytic activity of the M–I–I non-classical three-center bond of the pincer complex, which is involved in the oxidative addition of molecular iodine I2 to the metal center. According to our investigations the charge transfer from the metal to the σ* antibonding orbital of the I–I bond changes the 3c–4e character of the M–I–I three-center bond, which leads to weakening of the iodine I–I bond and strengthening of the metal–iodine M–I bond, facilitating in this way the oxidative addition of I2 to the metal. The charge transfer can be systematically modified by substitution at different places of the pincer complex and by different transition metals, changing the strength of both the M–I and the I2 bonds. We also modeled for the original pincer complex how solvents with different polarity influence the 3c–4e character of the M–I–I bond. Our results provide new guidelines for the design of pincer complexes with specific iodine–metal bond strengths and introduce the local vibrational mode analysis as an efficient tool to assess the bond strength in complexes.

Relative Strengths of a Pnicogen and a Tetrel Bond and Their Mutual Effects upon One Another

The ability of the T and Z atoms of TR₃ZR₂ to engage in a noncovalent interaction with NH₃ is assessed by DFT calculations, where the T atom refers to C, Si, and Ge; Z = As, Sb, and P; and substituents R = H and F. In most instances, the tetrel bond (TB) is both stronger and shorter than the pnicogen bond (ZB). These two bond strengths can be equalized, or preference shifted to the ZB, if F substituents are placed on the Z and H on the T atoms. Employing C as the T atom results in a very weak TB, with the ZB clearly favored energetically. The simultaneous formation of both TB and ZB weakens both, particularly the latter, but both bonds survive intact. Geometric and spectroscopic perturbations of the subunits reflect the two types of noncovalent bonds.

Comment on "Exploring nature and predicting strength of hydrogen bonds: A correlation analysis between atoms-in-molecules descriptors, binding energies, and energy components of symmetry-adapted perturbation theory"

We evaluate the correlation between binding energy (BE) and electron density ρ(r) at the bond critical point for 28 neutral hydrogen bonds, recently reported by Emamian and co-workers (J. Comput. Chem., 2019, 40, 2868). As an efficient tool, we use local stretching force constant k HB a derived from the local vibrational mode theory of Konkoli and Cremer. We compare the physical nature of BE versus k HB a , and provide an important explanation for cases with significant deviation in the BE- k HB a relation as well as in the BE-ρ(r) correlation. We also show that care has to be taken when different hydrogen bond strength measures are compared. The BE is a cumulative hydrogen bond strength measure while k HB a is a local measure of hydrogen bond strength covering different aspects of bonding. A simplified and unified description of hydrogen bonding is not always possible and needs an in-depth understanding of the systems involved.

The HOX⋯SO2 (X=F, Cl, Br, I) Binary Complexes: Implications for Atmospheric Chemistry

Sulfur dioxide and hypohalous acids (HOX, X=F, Cl, Br, I) are ubiquitous molecules in the atmosphere that are central to important processes like seasonal ozone depletion, acid rain, and cloud nucleation. We present the first theoretical examination of the HOX⋯SO₂ binary complexes and the associated trends due to halogen substitution. Reliable geometries were optimized at the CCSD(T)/aug-cc-pV(T+d)Z level of theory for HOF and HOCl complexes. The HOBr and HOI complexes were optimized at the CCSD(T)/aug-cc-pV(D+d)Z level of theory with the exception of the Br and I atoms which were modeled with an aug-cc-pwCVDZ-PP pseudopotential. 27 HOX⋯SO₂ complexes were characterized and the focal point method was employed to produce CCSDT(Q)/CBS interaction energies. Natural Bond Orbital analysis and Symmetry Adapted Perturbation Theory were used to classify the nature of each principle interaction. The interaction energies of all HOX⋯SO₂ complexes in this study ranged from 1.35 to 3.81 kcal mol^-1 . The single-interaction hydrogen bonded complexes spanned a range of 2.62 to 3.07 kcal mol^-1 , while the single-interaction halogen bonded complexes were far more sensitive to halogen substitution ranging from 1.35 to 3.06 kcal mol^-1 , indicating that the two types of interactions are extremely competitive for heavier halogens. Our results provide insight into the interactions between HOX and SO₂ which may guide further research of related systems.

Modeling Hydrogen Release from Water with Borane and Alane Catalysts: A Unified Reaction Valley Approach

The unified reaction valley approach combined with the local vibrational mode and ring puckering analysis is applied to investigate the hydrogen evolution from water in the presence of small hydrides such as BH₃, metal hydrides as AlH₃, and their derivatives. We studied a series of reactions involving BH₃, AlH₃, B₂H₆, Al₂H₆, and AlH₃BH₃ with one- and two-water molecules, considering multiple reaction paths. In addition, the influence of the aqueous medium was examined. A general reaction mechanism was identified for most of the reactions. Those that deviate could be associated with unusually high reaction barriers with no hydrogen release. The charge transfer along the reaction path suggests that a viable hydrogen release is achieved when the catalyst adopts the role of a charge donor during the chemical processes. The puckering analysis showed that twistboat and boat forms are the predominant configurations in the case of an intermediate six-membered ring formation, which influences the activation barrier. The local mode analysis was used as a tool to detect the H-H bond formation as well as to probe catalyst regenerability. Based on the correlation between the activation energy and the change in the charge separation for cleaving O-H and B(Al)-H bonds, two promising subsets of reactions could be identified along with prescriptions for lowering the reaction barrier individually with electron-donating/withdrawing substituents.

Halogen-bonded haloamine trimers - modelling the X3 synthon

Halogen-halogen bonded haloamine trimers serve as model systems for the X3 synthon present in numerous crystal structures and in two-dimensional self-assembled nanoarchitectures. Halogen bonds forming the synthon are often considered to display cooperativity. Synergy effects were previously found for halogen-halogen bonded bromoamine and iodoamine tetramers. In the present study comparison between haloamine cyclic trimers and tetramers is made. The cooperativity occurring in bromoamine and iodoamine clusters is significantly weaker in the case of the trimers. The present study demonstrates that the bromoamine and iodoamine trimers display weaker cooperativity due to a smaller number of synergy components in comparison to the corresponding tetramers of stronger cooperativity. Moreover, the halogen-halogen interactions in bromoamine and iodoamine dimers with the geometries of the corresponding trimers and tetramers are examined using energy decomposition analysis methods (supermolecular, canonical EDA and SAPT) and the Kohn-Sham molecular orbital model. The results of the analysis indicate that although the interaction energy values for the dimers of the different spatial arrangement are very close to each other, their origin is substantially different. For pairs with the geometry of the trimers orbital interactions and electrostatic attraction are both weaker than for the corresponding dimers with the geometry of the tetramers. This is especially important because both donor-acceptor interactions and electrostatic attraction were previously proven to be responsible for cooperative effects occurring in the bromoamine and iodoamine tetramers.

Computational analysis of vibrational frequencies and rovibrational spectroscopic constants of hydrogen sulfide dimer using MP2 and CCSD(T)

Previous studies have shown that the weakly bonded H₂S dimer demands high level quantum chemical calculations to reproduce experimental values. We investigated the hydrogen bonding of H₂S dimer using MP2 and CCSD(T) levels of theory in combination with aug-cc-pV(D,T,Q)Z basis sets. More precisely, the binding energies, potential energy curves, rovibrational spectroscopic constants, decomposition lifetime, and normal vibrational frequencies were calculated. In addition, we introduced the local mode analysis of Konkoli-Cremer to quantify the hydrogen bonding in the H₂S dimer as well as providing for the first time the comprehensive decomposition of normal vibrational modes into local modes contributions, and a decomposition lifetime based on rate constant. The local mode force constant of the H₂S dimer hydrogen bond is smaller than that of the water dimer, in accordance with the weaker hydrogen bonding in the H₂S dimer.

Tweaking the Charge Transfer: Bonding Analysis of Bismuth(III) Complexes with a Flexidentate Phosphane Ligand

To account for the charge transfer and covalent character in bonding between P and Bi centers, the electronic structures of [P(C6H4-o-CH2SCH3)3BiCln](3-n)+ (n = 0-3) model species have been investigated computationally. On the basis of this survey a synthetic target compound with a dative P→Bi bond has been selected. Consecutively, the highly reactive bismuth cage [P(C6H4-o-CH2SCH3)3Bi]3+ has been accessed experimentally and characterized. Importantly, our experiments (single-crystal X-ray diffraction and solid-state NMR spectroscopy) and computations (NBO and AIM analysis) reveal that the P···Bi bonding in this trication can be described as a dative bond. Here we have shown that our accordion-like molecular framework allows for tuning of the interaction between P and Bi centers.

Mutual Influence of Pnicogen Bonds and Beryllium Bonds: Energies and Structures in the Spotlight

Pnicogen bonds, which are weak noncovalent interactions (NCIs), can be significantly modified by the presence of beryllium bonds, one of the strongest NCIs known. We demonstrate the importance of this influence by studying ternary complexes in which both NCIs are present, that is, the ternary complexes formed by a nitrogen base (NH₃, NHCH₂, and NCH), a phosphine (fluorophosphane, PH₂F) and a beryllium derivative (BeH₂, BeF₂, BeCl₂, BeCO₃, and BeSO₄). Energies, structures, and nature of the chemical bonding in these complexes are studied by means of ab initio computational methods. The pnicogen bond between the nitrogen base and the phosphine and the beryllium bond between the fluorine atom of fluorophosphane and the beryllium derivative show large cooperativity effects both on energies and geometries, with dissociation energies up to 296 kJ mol^-1 and cooperativity up to 104 kJ mol^-1 in the most strongly bound complex, CH₂HN:PH₂F:BeSO₄. In the complexes between the strongest nitrogen bases and the strongest beryllium donors, phosphorus-shared and phosphorus-transfer bonds are found.

Local Vibrational Mode Analysis of π–Hole Interactions between Aryl Donors and Small Molecule Acceptors

11 aryl–lone pair and three aryl–anion π –hole interactions are investigated, along with the argon–benzene dimer and water dimer as reference compounds, utilizing the local vibrational mode theory, originally introduced by Konkoli and Cremer, to quantify the strength of the π –hole interaction in terms of a new local vibrational mode stretching force constant between the two engaged monomers, which can be conveniently used to compare different π –hole systems. Several factors have emerged which influence strength of the π –hole interactions, including aryl substituent effects, the chemical nature of atoms composing the aryl rings/ π –hole acceptors, and secondary bonding interactions between donors/acceptors. Substituent effects indirectly affect the π –hole interaction strength, where electronegative aryl-substituents moderately increase π –hole interaction strength. N-aryl members significantly increase π –hole interaction strength, and anion acceptors bind more strongly with the π –hole compared to charge neutral acceptors (lone–pair donors). Secondary bonding interactions between the acceptor and the atoms in the aryl ring can increase π –hole interaction strength, while hydrogen bonding between the π –hole acceptor/donor can significantly increase or decrease strength of the π –hole interaction depending on the directionality of hydrogen bond donation. Work is in progress expanding this research on aryl π –hole interactions to a large number of systems, including halides, CO, and OCH3− as acceptors, in order to derive a general design protocol for new members of this interesting class of compounds.

Exploring the Mechanism of Catalysis with the Unified Reaction Valley Approach (URVA)—A Review

The unified reaction valley approach (URVA) differs from mainstream mechanistic studies, as it describes a chemical reaction via the reaction path and the surrounding reaction valley on the potential energy surface from the van der Waals region to the transition state and far out into the exit channel, where the products are located. The key feature of URVA is the focus on the curving of the reaction path. Moving along the reaction path, any electronic structure change of the reacting molecules is registered by a change in their normal vibrational modes and their coupling with the path, which recovers the curvature of the reaction path. This leads to a unique curvature profile for each chemical reaction with curvature minima reflecting minimal change and curvature maxima, the location of important chemical events such as bond breaking/forming, charge polarization and transfer, rehybridization, etc. A unique decomposition of the path curvature into internal coordinate components provides comprehensive insights into the origins of the chemical changes taking place. After presenting the theoretical background of URVA, we discuss its application to four diverse catalytic processes: (i) the Rh catalyzed methanol carbonylation—the Monsanto process; (ii) the Sharpless epoxidation of allylic alcohols—transition to heterogenous catalysis; (iii) Au(I) assisted [3,3]-sigmatropic rearrangement of allyl acetate; and (iv) the Bacillus subtilis chorismate mutase catalyzed Claisen rearrangement—and show how URVA leads to a new protocol for fine-tuning of existing catalysts and the design of new efficient and eco-friendly catalysts. At the end of this article the pURVA software is introduced. The overall goal of this article is to introduce to the chemical community a new protocol for fine-tuning existing catalytic reactions while aiding in the design of modern and environmentally friendly catalysts.

The challenge of non-covalent interactions: theory meets experiment for reconciling accuracy and interpretation

In the past decade, many gas-phase spectroscopic investigations have focused on the understanding of the nature of weak interactions in model systems. Despite the fact that non-covalent interactions play a key role in several biological and technological processes, their characterization and interpretation are still far from being satisfactory. In this connection, integrated experimental and computational investigations can play an invaluable role. Indeed, a number of different issues relevant to unraveling the properties of bulk or solvated systems can be addressed from experimental investigations on molecular complexes. Focusing on the interaction of biological model systems with solvent molecules (e.g., water), since the hydration of the biomolecules controls their structure and mechanism of action, the study of the molecular properties of hydrated systems containing a limited number of water molecules (microsolvation) is the basis for understanding the solvation process and how structure and reactivity vary from gas phase to solution. Although hydrogen bonding is probably the most widespread interaction in nature, other emerging classes, such as halogen, chalcogen and pnicogen interactions, have attracted much attention because of the role they play in different fields. Their understanding requires, first of all, the characterization of the directionality, strength, and nature of such interactions as well as a comprehensive analysis of their competition with other non-covalent bonds. In this review, it is shown how state-of-the-art quantum-chemical computations combined with rotational spectroscopy allow for fully characterizing intermolecular interactions taking place in molecular complexes from both structural and energetic points of view. The transition from bi-molecular complex to microsolvation and then to condensed phase is shortly addressed.

In Situ Assessment of Intrinsic Strength of X-I⋯OA-Type Halogen Bonds in Molecular Crystals with Periodic Local Vibrational Mode Theory

Periodic local vibrational modes were calculated with the rev-vdW-DF2 density functional to quantify the intrinsic strength of the X-I⋯OA-type halogen bonding (X = I or Cl; OA: carbonyl, ether and N-oxide groups) in 32 model systems originating from 20 molecular crystals. We found that the halogen bonding between the donor dihalogen X-I and the wide collection of acceptor molecules OA features considerable variations of the local stretching force constants (0.1–0.8 mdyn/Å) for I⋯O halogen bonds, demonstrating its powerful tunability in bond strength. Strong correlations between bond length and local stretching force constant were observed in crystals for both the donor X-I bonds and I⋯O halogen bonds, extending for the first time the generalized Badger’s rule to crystals. It is demonstrated that the halogen atom X controlling the electrostatic attraction between the σ-hole on atom I and the acceptor atom O dominates the intrinsic strength of I⋯O halogen bonds. Different oxygen-containing acceptor molecules OA and even subtle changes induced by substituents can tweak the n→σ∗(X-I) charge transfer character, which is the second important factor determining the I⋯O bond strength. In addition, the presence of the second halogen bond with atom X of the donor X-I bond in crystals can substantially weaken the target I⋯O halogen bond. In summary, this study performing the in situ measurement of halogen bonding strength in crystalline structures demonstrates the vast potential of the periodic local vibrational mode theory for characterizing and understanding non-covalent interactions in materials.

Modified Density Functional Dispersion Correction for Inorganic Layered MFX Compounds (M = Ca, Sr, Ba, Pb and X = Cl, Br, I)

MFX (M = Ca, Ba, Sr, Pb and X = Cl, Br, I) compounds have received considerable attention due to their technological application as X-ray detectors, pressure sensors, and optical data storage materials, when doped with rare-earth ions. MFX compounds belong to the class of layered materials with a tetragonal Matlockite crystal structure, characterized by weakly stacked double-halide layers along the crystallographic c-axis. These layers predominantly determine phase transitions, elastic, and mechanical properties. However, the correct description of the lattice parameter c is a challenge for most standard DFT functionals, which tend to overestimate the lattice parameter c. Because of the weak interactions between the halide layers, dispersion-corrected functionals seem to be a better choice. We investigated 11 different inorganic layered MFX compounds for which experimental data are available, with standard and dispersion-corrected functionals to assess their performance in reproducing the lattice parameter c, structural, and vibrational properties of the MFX compounds. Our results revealed that these functionals do not describe the weak interactions between the halide layers in a balanced way. Therefore, we modified Grimme's popular DFT-D2 dispersion correction scheme in two different ways by (i) replacing the dispersion coefficients and van der Waals radii with those of noble gas atoms or (ii) increasing the van der Waals radii of the MFX atoms up to 40%. Comparison with the available experimental data revealed that the latter approach applied to the PBE (Perdew-Burke-Ernzerhof)-D2 functional with 30% increased van der Waals radii, which we coined PBE-D2* (S_rvdw 1.30) is best suited to fine-tune the description of the weak interlayer interactions in MFX compounds, thus significantly improving the description of their structural, vibrational, and mechanical properties. Work is in progress applying this new, computationally inexpensive scheme to other inorganic layered compounds and periodic systems with weakly stacked layers.

Halogen⋯halogen interactions in decahalo-closo-carboranes: CSD analysis and theoretical study

We theoretically (PBE0-D3/def2TZVP) and experimentally (CSD analysis) demonstrate the importance of “like–like” halogen interactions for the stability of several decahalo-closo-carborane dimers.

Metal-Halogen Bonding Seen through the Eyes of Vibrational Spectroscopy

Incorporation of a metal center into halogen-bonded materials can efficiently fine-tune the strength of the halogen bonds and introduce new electronic functionalities. The metal atom can adopt two possible roles: serving as halogen acceptor or polarizing the halogen donor and acceptor groups. We investigated both scenarios for 23 metal–halogen dimers trans-M(Y₂)(NC₅H₄X-3)₂ with M = Pd(II), Pt(II); Y = F, Cl, Br; X = Cl, Br, I; and NC₅H₄X-3 = 3-halopyridine. As a new tool for the quantitative assessment of metal–halogen bonding, we introduced our local vibrational mode analysis, complemented by energy and electron density analyses and electrostatic potential studies at the density functional theory (DFT) and coupled-cluster single, double, and perturbative triple excitations (CCSD(T)) levels of theory. We could for the first time quantify the various attractive contacts and their contribution to the dimer stability and clarify the special role of halogen bonding in these systems. The largest contribution to the stability of the dimers is either due to halogen bonding or nonspecific interactions. Hydrogen bonding plays only a secondary role. The metal can only act as halogen acceptor when the monomer adopts a (quasi-)planar geometry. The best strategy to accomplish this is to substitute the halo-pyridine ring with a halo-diazole ring, which considerably strengthens halogen bonding. Our findings based on the local mode analysis provide a solid platform for fine-tuning of existing and for design of new metal–halogen-bonded materials.

Quantitative Assessment of B-B-B, B-Hb -B, and B-Ht Bonds: From BH3 to B12 H12 2

We report the thermodynamic stabilities and the intrinsic strengths of three-center-two-electron B-B-B and B-H_b -B bonds ( H b : bridging hydrogen), and two-center-two-electron B-H_t bonds ( H t : terminal hydrogen) which can be served as a new, effective tool to determine the decisive role of the intermediates of hydrogenation/dehydrogenation reactions of borohydride. The calculated heats of formation were obtained with the G4 composite method and the intrinsic strengths of B-B-B, B-H_b -B, and B-H_t bonds were derived from local stretching force constants obtained at the B3LYP-D2/cc-pVTZ level of theory for 21 boron-hydrogen compounds, including 19 intermediates. The Quantum Theory of Atoms in Molecules (QTAIM) was used to deepen the inside into the nature of B-B-B, B-H_b -B, and B-H_t bonds. We found that all of the experimentally identified intermediates hindering the reversibility of the decomposition reactions are thermodynamically stable and possess strong B-B-B, B-H_b -B, and B-H_t bonds. This proves that thermodynamic data and intrinsic B-B-B, B-H_b -B, and B-H_t bond strengths form a new, effective tool to characterize new (potential) intermediates and to predict their role for the reversibility of the hydrogenation/dehydrogenation reactions.

A Continuum from Halogen Bonds to Covalent Bonds: Where Do λ3 Iodanes Fit?

The intrinsic bonding nature of λ 3 -iodanes was investigated to determine where its hypervalent bonds fit along the spectrum between halogen bonding and covalent bonding. Density functional theory with an augmented Dunning valence triple zeta basis set ( ω B97X-D/aug-cc-pVTZ) coupled with vibrational spectroscopy was utilized to study a diverse set of 34 hypervalent iodine compounds. This level of theory was rationalized by comparing computational and experimental data for a small set of closely-related and well-studied iodine molecules and by a comparison with CCSD(T)/aug-cc-pVTZ results for a subset of the investigated iodine compounds. Axial bonds in λ 3 -iodanes fit between the three-center four-electron bond, as observed for the trihalide species IF 2 − and the covalent FI molecule. The equatorial bonds in λ 3 -iodanes are of a covalent nature. We explored how the equatorial ligand and axial substituents affect the chemical properties of λ 3 -iodanes by analyzing natural bond orbital charges, local vibrational modes, the covalent/electrostatic character, and the three-center four-electron bonding character. In summary, our results show for the first time that there is a smooth transition between halogen bonding → 3c–4e bonding in trihalides → 3c–4e bonding in hypervalent iodine compounds → covalent bonding, opening a manifold of new avenues for the design of hypervalent iodine compounds with specific properties.

Weak Pnictogen Bond with Bismuth: Experimental Evidence Based on Bi-P Through-Space Coupling

To study pnictogen bonding involving bismuth, flexible accordion-like molecular complexes of the composition [P(C6 H4 -o-CH2 SCH3 )3 BiX3 ], (X=Cl, Br, I) have been synthesised and characterised. The strength of the weak and mainly electrostatic interaction between the Bi and P centres strongly depends on the character of the halogen substituent on bismuth, which is confirmed by single-crystal X-ray diffraction analyses, DFT and ab initio computations. Significantly, 209 Bi-31 P through-space coupling (J=2560 Hz) is observed in solid-state 31 P NMR spectra, which is so far unprecedented in the literature, delivering direct information on the magnitude of this pnictogen interaction.

Halogen Bonding: A Halogen-Centered Noncovalent Interaction Yet to Be Understood

In addition to the underlying basic concepts and early recognition of halogen bonding, this paper reviews the conflicting views that consistently appear in the area of noncovalent interactions and the ability of covalently bonded halogen atoms in molecules to participate in noncovalent interactions that contribute to packing in the solid-state. It may be relatively straightforward to identify Type-II halogen bonding between atoms using the conceptual framework of σ-hole theory, especially when the interaction is linear and is formed between the axial positive region (σ-hole) on the halogen in one monomer and a negative site on a second interacting monomer. A σ-hole is an electron density deficient region on the halogen atom X opposite to the R–X covalent bond, where R is the remainder part of the molecule. However, it is not trivial to do so when secondary interactions are involved as the directionality of the interaction is significantly affected. We show, by providing some specific examples, that halogen bonds do not always follow the strict Type-II topology, and the occurrence of Type-I and -III halogen-centered contacts in crystals is very difficult to predict. In many instances, Type-I halogen-centered contacts appear simultaneously with Type-II halogen bonds. We employed the Independent Gradient Model, a recently proposed electron density approach for probing strong and weak interactions in molecular domains, to show that this is a very useful tool in unraveling the chemistry of halogen-assisted noncovalent interactions, especially in the weak bonding regime. Wherever possible, we have attempted to connect some of these results with those reported previously. Though useful for studying interactions of reasonable strength, IUPAC’s proposed “less than the sum of the van der Waals radii” criterion should not always be assumed as a necessary and sufficient feature to reveal weakly bound interactions, since in many crystals the attractive interaction happens to occur between the midpoint of a bond, or the junction region, and a positive or negative site.

The chalcogen bond: can it be formed by oxygen?

This study theoretically investigates the possibility of oxygen-centered chalcogen bonding in several complexes. Shown in the graph is such a bonding scenario formed between the electrophile on O in OF₂ and the nucleophile on O in H₂CO.