Re-evaluation of the bond length-bond strength rule: The stronger bond is not always the shorter bond

Solvent Effect on Dative and Ionic Bond Strengths: A Unified Theory from Potential Analysis and Valence-Bond Computations

Solvent molecules interact with a solute through various intermolecular forces. Here we employed a potential energy surface (PES) analysis to interpret the solvent-induced variations in the strengths of dative (Me3NBH3) and ionic (LiCl) bonds, which possess both ionic and covalent (neutral) characteristics. The change of a bond is driven by the gradient (force) of the solvent-solute interaction energy with respect to the focused bond length. Positive force shortens the bond length and increases the bond force constant, leading to a blue-shift of the bond stretching vibrational frequency upon solvation. Conversely, negative force elongates the bond, resulting in a reduced bond force constant and red-shift of the stretching vibrational frequency. The different responses of Me3NBH3 and LiCl to solvation are studied with valence bond (VB) theory, as Me3NBH3 and LiCl are dominated by the neutral covalent VB structure and the ionic VB structure, respectively. The dipole moment of an ionic VB structure increases along the increasing bond distance, while the dipole moment of a neutral covalent VB structure increases with the decreasing bond distance. The roles of the dominating VB structures are further examined by the geometry optimizations and frequency calculations with the block-localized wavefunction (BLW) method.

Extraction of uranyl from spent nuclear fuel wastewater via complexation-a local vibrational mode study

CONTEXT

The efficient extraction of uranyl from spent nuclear fuel wastewater for subsequent reprocessing and reuse is an essential effort toward minimization of long-lived radioactive waste. N-substituted amides and Schiff base ligands are propitious candidates, where extraction occurs via complexation with the uranyl moiety. In this study, we extensively probed chemical bonding in various uranyl complexes, utilizing the local vibrational modes theory alongside QTAIM and NBO analyses. We focused on (i) the assessment of the equatorial O-U and N-U bonding, including the question of chelation, and (ii) how the strength of the axial U = O bonds of the uranyl moiety changes upon complexation. Our results reveal that the strength of the equatorial uranium-ligand interactions correlates with their covalent character and with charge donation from O and N lone pairs into the vacant uranium orbitals. We also found an inverse relationship between the covalent character of the equatorial ligand bonds and the strength of the axial uranium-oxygen bond. In summary, our study provides valuable data for a strategic modulation of N-substituted amide and Schiff base ligands towards the maximization of uranyl extraction.

METHOD

Quantum chemistry calculations were performed under the PBE0 level of theory, paired with the relativistic NESCau Hamiltonian, currently implemented in Cologne2020 (interfaced with Gaussian16). Wave functions were expanded in the cc-pwCVTZ-X2C basis set for uranium and Dunning's cc-pVTZ for the remaining atoms. For the bonding properties, we utilized the package LModeA in the local modes analyses, AIMALL in the QTAIM calculations, and NBO 7.0 for the NBO analyses.

A multi-descriptor analysis of substituent effects on the structure and aromaticity of benzene derivatives: π-Conjugation versus charge effects

This work provides a detailed multi-component analysis of aromaticity in monosubstituted (X = CH3, CH 2 - , CH 2 + , NH2, NH-, NH+, OH, O-, and O+) and para-homodisubstituted (X = CH3, CH2, NH2, NH, OH, and O) benzene derivatives. We investigate the effects of substituents using single-reference (B3LYP/DFT) and multireference (CASSCF/MRCI) methods, focusing on structural (HOMA), vibrational (AI(vib)), topological (ELFπ), electronic (MCI), magnetic (NICS), and stability (S0-T1 splitting) properties. The findings reveal that appropriate π-electron-donating and π-electron-accepting substituents with suitable size and symmetry can interact with the π-system of the ring, significantly influencing π-electron delocalization. While the charge factor has a minimal impact on π-electron delocalization, the presence of a pz orbital capable of interacting with the π-electron delocalization is the primary factor leading to a deviation from the typical aromaticity characteristics observed in benzene.

Iron-histidine bonding in bishistidyl hemoproteins-A local vibrational mode study

We investigated the intrinsic strength of distal and proximal FeN bonds for both ferric and ferrous oxidation states of bishistidyl hemoproteins from bacteria, animals, human, and plants, including two cytoglobins, ten hemoglobins, two myoglobins, six neuroglobins, and six phytoglobins. As a qualified measure of bond strength, we used local vibrational force constants ka (FeN) based on local mode theory developed in our group. All calculations were performed with a hybrid QM/MM ansatz. Starting geometries were taken from available x-ray structures. ka (FeN) values were correlated with FeN bond lengths and covalent bond character. We also investigated the stiffness of the axial NFeN bond angle. Our results highlight that protein effects are sensitively reflected in ka (FeN), allowing one to compare trends in diverse protein groups. Moreover, ka (NFeN) is a perfect tool to monitor changes in the axial heme framework caused by different protein environments as well as different Fe oxidation states.

Exploring Jahn-Teller distortions: a local vibrational mode perspective

The characterization of normal mode (CNM) procedure coupled with an adiabatic connection scheme (ACS) between local and normal vibrational modes, both being a part of the Local Vibrational Mode theory developed in our group, can identify spectral changes as structural fingerprints that monitor symmetry alterations, such as those caused by Jahn-Teller (JT) distortions. Employing the PBE0/Def2-TZVP level of theory, we investigated in this proof-of-concept study the hexaaquachromium cation case,[ Cr ( OH 2 ) 6 ] 3 + /[ Cr ( OH 2 ) 6 ] 2 + , as a commonly known example for a JT distortion, followed by the more difficult ferrous and ferric hexacyanide anion case,[ Fe ( CN ) 6 ] 4 - /[ Fe ( CN ) 6 ] 3 - . We found that in both cases CNM of the characteristic normal vibrational modes reflects delocalization consistent with high symmetry and ACS confirms symmetry breaking, as evidenced by the separation of axial and equatorial group frequencies. As underlined by the Cremer-Kraka criterion for covalent bonding, from[ Cr ( OH 2 ) 6 ] 3 + to[ Cr ( OH 2 ) 6 ] 2 + there is an increase in axial covalency whereas the equatorial bonds shift toward electrostatic character. From[ Fe ( CN ) 6 ] 4 - to[ Fe ( CN ) 6 ] 3 - we observed an increase in covalency without altering the bond nature. Distinct π back-donation disparity could be confirmed by comparison with the isolated CN- system. In summary, our study positions the CNM/ACS protocol as a robust tool for investigating less-explored JT distortions, paving the way for future applications.

Photocatalytic Sulfonyl Peroxidation of Alkenes via Deamination of N-Sulfonyl Ketimines

A photocatalytic three-component sulfonyl peroxidation of alkenes with N-sulfonyl ketimines and tert-butyl hydroperoxide is reported. The reaction takes place via the photoinduced EnT process, which allows the efficient synthesis of a variety of β-peroxyl sulfones under mild reaction conditions in the absence of a transition metal catalyst. The downstream derivatizations of the peroxides were also performed. Furthermore, the utility of this protocol was manifested by the synthesis of 11β-HSD1 inhibitor and the antiprostate cancer drug bicalutamide.

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.

Extreme NMR shielding in fluoro-nitrogen cations

The structure and NMR shielding of a set of N-F containing cations is reported to near-quantitative accuracy from extensive ab initio calculations. Currently, the shortest experimentally confirmed N-F bond is 1.2461(10) Å in NNF+, however CCSD(T)-F12b/cc-pVQZ-F12 optimised geometries suggest that even shorter N-F bonds are possible for both monocations (1.236 Å, HNF+) and dications (1.098 Å, NF2+). NMR shielding constants have been calculated in a composite manner with individual components from coupled-cluster expansions up to CCSDTQP and basis sets up to aug-cc-pCV8Z, together with vibrational and relativistic corrections. 15N and 19F NMR chemical shifts correlate well with available experimental data. Extreme 19F chemical shifts are predicted for HNF+ (1628.9 ppm) and NH2F2+ (1298.0 ppm), which are by far the largest 19F chemical shifts ever reported and well outside the known range of +865 ppm (F2O2) to -448 ppm (ClF). The 15N chemical shift of -1283.07 ppm in HNF+ is similarly extreme, being well outside the known range of 15N chemical shifts of -730 to 260 ppm (CH3NO2 reference). This work highlights the application of state-of-the-art theoretical techniques, and provides accurate NMR properties of both isolated and yet unknown N-F cations, which can serve to guide and supplement NMR experimentation.

Noncovalent π Interactions in Mutated Aquomet-Myoglobin Proteins: A QM/MM and Local Vibrational Mode Study

Protein dynamics and function is strongly connected to the energy flow taking place. Myoglobin (Mb) and its mutations are ideal systems to study the process of vibrational energy transfer (VET) at the molecular level. Anti-Stokes ultraviolet resonance Raman studies using a tryptophan (Trp) probe, introduced at different Mb positions by amino acid replacement, have suggested that the amount of VET depends on the position of the Trp probe relative to the heme group. Inspired by this experimental work, we explored the strength of noncovalent π interactions, as well as covalent interactions for both the axial and distal ligands bound to iron in aquomet-Mb with the local vibrational mode analysis (LMA), originally developed by Konkoli and Cremer. Two sets of noncovalent interactions were investigated: (1) the interaction between the water ligand and Trp rings and (2) the interaction between the Trp and the porphyrin rings of the heme group. We assessed the strength of these noncovalent interactions via a special local mode force constant. Various Trp-modified water-bound ferric Mb proteins in the ground state were studied (6 in total) using gas-phase and QM/MM calculations followed by LMA. Our results disclose that VET is indeed dependent on the position of the Trp probe relative to the heme group but also on the tautomeric nature of distal histidine. They provide new guidelines on how to assess noncovalent π interactions in proteins utilizing LMA and how to use these data to explore VET, and more generally protein dynamics and function.

EnT-Mediated N-S Bond Homolysis of a Bifunctional Reagent Leading to Aliphatic Sulfonyl Fluorides

Sulfur(VI) fluoride exchange (SuFEx) gives rise to a plethora of high-valent sulfur linkages; however, the availability of (aliphatic) sulfonyl fluoride manifolds lag behind, owing to the limited sources of introducing the SO₂F moiety via a classical two-electron approach. Recently, radical-based methodologies have emerged as a complementary strategy to increase the diversity of accessible click partners. In this work, synthesis of a bench-stable sulfamoyl fluoride reagent is presented, which may undergo sigma-bond homolysis upon visible-light-induced sensitization to form protected β-amino sulfonyl fluorides from alkene feedstocks. Notably, this offers an appealing strategy to access various building blocks for peptido sulfonyl fluorides, relevant in a medicinal chemistry context, as well as an intriguing entry to β-ammonium sulfonates and β-sultams, from alkenes. Densely functionalized 1,3-sultones were obtained by employing allyl alcohols as substrates. Surprisingly, allyl chloride-derived β-imino sulfonyl fluoride underwent S-O bond formation and ring closure to yield rigid cyclopropyl β-imino sulfonate ester under SuFEx conditions. Furthermore, by engaging a thiol-based hydrogen atom donor in the reaction, the reactivity of the same reagent can be tuned toward the direct synthesis of aliphatic sulfonyl fluorides. Mechanistic experiments indicate an energy transfer (EnT)-mediated process. The transient sulfonyl fluoride radical adds to the alkene and product formation occurs upon either radical-radical coupling or hydrogen atom transfer (HAT), respectively.

The effect of immediate environment on bond strength of different bond types-A valence bond study

The ability to design catalysis largely depends on our understanding of the electrostatic effect of the surrounding on the bonds participating in the reaction. Here, we used a simplistic model of point charges (PCs) to determine a set of rules guiding how to construct PC-bond arrangement that can strengthen or weaken different chemical bonds. Using valence bond theory to calculate the in situ bond energies, we show that the effect of the PC mainly depends on the bond's dipole moment irrespective of its type (being covalent or charge shift). That is, polar bonds are getting stronger or weaker depending on the sign and location of the PC, whereas non- or weakly polar bonds become stronger or weaker depending only on the location of the PC and to a smaller extent compared with polar bonds. We also show that for polar bonds, the maximal bond strengthening and weakening effect can be achieved when the PC is placed along the bond axis, as close as possible to the more and less polarizable atom/fragment, respectively. Finally, due to the stabilizing effects of polarizability, we show that, overall, it is easier to cause bond strengthening compared with bond weakening. Particularly, for polar bonds, bond strengthening is larger than bond weakening obtained by an oppositely signed PC. These rules should be useful in the future design of catalysis in, e.g., enzyme active sites.

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.

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.

A practical fluorosulfonylating platform via photocatalytic imidazolium-based SO2F radical reagent

Sulfonyl fluorides are key components in the fields of chemical biology, materials science and drug discovery. In this line, the highly active SO₂F radical has been employed for the construction of sulfonyl fluorides, but the utilization of gaseous ClSO₂F as radical precursor is limited due to the tedious and hazardous preparation. Meanwhile, the synthesis of sulfonyl fluorides from inert SO₂F₂ gas through a fluorosulfonyl radical (·SO₂F) process has met with inevitable difficulties due to the high homolytic bond dissociation energy of the S(VI)-F bond. Here we report a radical fluorosulfonylation strategy for the stereoselective synthesis of alkenyl sulfonyl fluorides and functional alkyl sulfonyl fluorides with an air-stable crystalline benzimidazolium fluorosulfonate cationic salt reagent. This bench-stable redox-active reagent offers a useful and operational protocol for the radical fluorosulfonylation of unsaturated hydrocarbons with good yield and high stereoselectivity, which can be further transformed into valuable functional SO₂F moieties.

Sulfonyl fluorides have potential application in chemical biology, materials science, and drug discovery, but their preparation remains challenging. Here, the authors report an air-stable fluorosulfonylating reagent that enables the radical fluorosulfonylation, hydrofluorosulfonylation and migratory SO₂F-difunctionalization of unsaturated hydrocarbons to construct a variety of sulfonyl fluoride compounds.

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.

Halogen-Exchange Reactions of Iminophosphonamido-Chlorosilylenes with Alkali Halides: Convenient Synthesis of Heavier Halosilylenes

Halogen-substituted silylenes are an important building block for synthesizing silicon-based low-valent and multiple-bond species. However, the number of reports on heavier halosilylenes that contain bromine and iodine is still limited. Here, we present a convenient synthesis for bromo- and iodosilylenes supported by an iminophosphonamide ligand. The heavier halosilylenes [Ph₂P(^tBuN)₂]SiX (2: X = Br, 3: X = I) were successfully synthesized via the halogen-exchange reaction of chlorosilylene 1 with alkali halides in THF. As a demonstration of the reactivity of 2 and 3, oxidative addition reactions of 2 and 3 with elemental selenium in C₆D₆ afforded the corresponding bromo- (5) or iodosilylene-selone (6) as colorless crystals. The molecular structures of 2, 3, 5, and 6 were fully characterized by spectroscopic means and single-crystal X-ray diffraction analysis. Furthermore, the effects of the halogen atom on the electronic state of halosilylenes 1-3 and halosilylene-selones 4-6 were investigated using density functional theory (DFT) calculations.

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.

A photoinduced transient activating strategy for late-stage chemoselective C(sp3)–H trifluoromethylation of azines

The direct functionalization of C(sp³)–H bonds is an ultimately ideal synthetic strategy with high atom economy and step efficiency. However, the direct trifluoromethylation of electron-deficient heteroaryl adjacent C(sp³)–H bonds remains a formidable challenge. We have described a transient activating strategy involving a Tf-shift process and π–π stacking interaction for catalyst-free direct benzylic C(sp³)–H trifluoromethylation of azines, such as pyridine, pyrimidine, quinoline, dihydropyridinone, tetrahydroisoquinoline and tetrahydroquinazoline, with an air-stable crystalline imidazolium sulfonate reagent IMDN-Tf. This bench-stable cationic reagent offers a scalable and practical protocol for the late-stage modification of drug molecules with high site selectivity, which avoids the prefunctionalization and the use of stoichiometric metals and strong oxidants. Furthermore, comprehensive mechanistic studies revealed the determining effect of π–π stacking for the activation of azinylic C(sp³)–H bonds.

Late-stage C(sp³)–H functionalization of unactivated azines: the traceless Tf switching process offers ample opportunities for site-selective derivatization of heteroaryls, allowing for the rapid increase of molecular complexity.

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.

The Chemical Bond: When Atom Size Instead of Electronegativity Difference Determines Trend in Bond Strength

We have quantum chemically analyzed element-element bonds of archetypal Hn X-YHn molecules (X, Y=C, N, O, F, Si, P, S, Cl, Br, I), using density functional theory. One purpose is to obtain a set of consistent homolytic bond dissociation energies (BDE) for establishing accurate trends across the periodic table. The main objective is to elucidate the underlying physical factors behind these chemical bonding trends. On one hand, we confirm that, along a period (e. g., from C-C to C-F), bonds strengthen because the electronegativity difference across the bond increases. But, down a period, our findings constitute a paradigm shift. From C-F to C-I, for example, bonds do become weaker, however, not because of the decreasing electronegativity difference. Instead, we show that the effective atom size (via steric Pauli repulsion) is the causal factor behind bond weakening in this series, and behind the weakening in orbital interactions at the equilibrium distance. We discuss the actual bonding mechanism and the importance of analyzing this mechanism as a function of the bond distance.

On the formation of CN bonds in Titan's atmosphere-a unified reaction valley approach study

In this work, we investigated the formation of protonated hydrogen cyanide HCNH⁺ and methylene amine cation CH[Formula: see text] (both identified in Titan's upper atmosphere) from three different pathways which stem from the interaction between CH₄ and N⁺(³P). As a mechanistic tool, we used the Unified Reaction Valley Approach (URVA) complemented with the Local Mode Analysis (LMA) assessing the strength of the CN bonds formed in these reactions. Our URVA studies could provide a comprehensive overview on bond formation/cleavage processes relevant to the specific mechanism of eight reactions R1- R8 that occur across the three pathways. In addition, we could explain the formation of CH[Formula: see text] and the appearance of HCNH⁺ and CHNH[Formula: see text] along these paths. Although only smaller molecules are involved in these reactions including isomerization, hydrogen atom abstraction, and hydrogen molecule capture, we found a number of interesting features, such as roaming in reaction R3 or the primary interaction of H₂ with the carbon atom in HCNH⁺ in reaction R8 followed by migration of one of the H₂ hydrogen atoms to the nitrogen which is more cost effective than breaking the HH bond first; a feature often found in catalysis. In all cases, charge transfer between carbon and nitrogen could be identified as a driving force for the CN bond formation. As revealed by LMA, the CN bonds formed in reactions R1-R8 cover a broad bond strength range from very weak to very strong, with the CN bond in protonated hydrogen cyanide HCNH⁺ identified as the strongest of all molecules investigated in this work. Our study demonstrates the large potential of both URVA and LMA to shed new light into these extraterrestrial reactions to help better understand prebiotic processes as well as develop guidelines for future investigations involving areas of complex interstellar chemistry. In particular, the formation of CN bonds as a precursor to the extraterrestrial formation of amino acids will be the focus of future investigations. Formation of CN bonds in Titan's atmosphere visualized via the reaction path curvature.

Relating Bond Strength and Nature to the Thermodynamic Stability of Hypervalent Togni-Type Iodine Compounds

The bond strength and nature of a set of 32 Togni-like reagents have been investigated at the M062X/def2-TZVP(D) level of theory in acetonitrile described with the SMD continuum solvent model, to rationalize the main factors responsible for their thermodynamic stability in different conformations, and trifluoromethylation capabilities. For the assessment of bond strength, we utilized local stretching force constants and associated bond strength orders, complemented with local features of the electron density to access the nature of the bonds. Bond dissociation energies varied from 31.6 to 79.9 kcal/mol depending on the polarizing power of the ligand trans to CF3 . Based on the analysis of the Laplacian of the density, we propose that the charge-shift bond character plays an important role in the stability of the molecules studied, especially for those containing I-O bonds. New insights on the trans influence and on possible ways to fine-tune the stability of these reagents are provided.

Highs and Lows of Bond Lengths: Is There Any Limit?

Two distinct points on the potential energy curve (PEC) of a pairwise interaction, the zero-energy crossing point and the point where the stretching force constant vanishes, allow us to anticipate the range of possible distances between two atoms in diatomic, molecular moieties and crystalline systems. We show that these bond-stability boundaries are unambiguously defined and correlate with topological descriptors of electron-density-based scalar fields, and can be calculated using generic PECs. Chemical databases and quantum-mechanical calculations are used to analyze a full set of diatomic bonds of atoms from the s-p main block. Emphasis is placed on the effect of substituents in C-C covalent bonds, concluding that distances shorter than 1.14 Å or longer than 2.0 Å are unlikely to be achieved, in agreement with ultra-high-pressure data and transition-state distances, respectively. Presumed exceptions are used to place our model in the correct framework and to formulate a conjecture for chained interactions, which offers an explanation for the multimodal histogram of O-H distances reported for hundreds of chemical systems.

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.

Redox-active benzimidazolium sulfonamides as cationic thiolating reagents for reductive cross-coupling of organic halides

Redox-active benzimidazolium sulfonamides as thiolating reagents have been developed for reductive C–S bond coupling. The IMDN-SO₂R reagent provides a bench-stable cationic precursor to generate a portfolio of highly active N–S intermediates, which can be successfully applied in cross-electrophilic coupling with various organic halides. The employment of an electrophilic sulfur source solved the problem of catalyst deactivation and avoided odorous thiols, featuring practical conditions, broad substrate scope, and excellent tolerance.

Redox-active benzimidazolium sulfonamides as thiolating reagents have been developed for reductive C–S bond coupling.

Intermolecular Interactions Involving Heavy Alkenes H2Si=TH2 (T = C, Si, Ge, Sn, Pb) with H2O and HCl: Tetrel Bond and Hydrogen Bond

The intermolecular interactions between the heavy alkenes H₂Si=TH₂ (T = C, Si, Ge, Sn, Pb) and H₂O or HCl have been explored at the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVDZ level. The various hydrogen bond (HB) and tetrel bond (TB) complexes can be located on the basis of molecular electrostatic potential maps of the isolated monomers. The competition between TB and HB interactions has been investigated through the relaxed potential energy surface scan. The results indicate that the HB complexes become more and more unstable relative to the TB complexes with the increase of the T atomic number, and cannot even retain as a minimum in some cases, for H₂Si=TH2···H₂O systems. In contrast, the HB complexes are generally more stable than TB complexes, and the TB complexes exhibit rather weak binding strength, for H₂Si=TH₂···HCl systems. The majority of the TB complexes formed between H₂Si=TH₂ and H₂O possesses very strong binding strength with covalent characteristics. The noncovalent TB complexes can be divided into two types on the basis of the orbital interactions: π-hole complexes, with binding angles ranging from 91 to 111°, and hybrid σ/π-hole complexes, with binding angles ranging from 130 to 165°. The interplay between different molecular interactions has been explored, and an interesting result is that the covalent TB interaction is significantly abated and becomes noncovalent because of the competitive effect.

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.

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.

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.

Integrated redox-active reagents for photoinduced regio- and stereoselective fluorocarboborylation

Vinylboronates and alkylboronates are key components in variegated transformations in all facets of chemical science. The synthesis of vinylboronates and alkylboronates suffers from step-tedious and poor stereoselective procedures. We have developed a regulated radical difunctionalization strategy for the construction of fluorine-containing vinylboronates and alkylboronates with an integrated redox-active reagent IMDN-SO2RF. This bench-stable imidazolium sulfonate cationic salt offers a scalable and operational protocol for the fluoroalkylation-borylation of unsaturated hydrocarbons in a high regio- and stereoselective manner. The products can be further transformed into valuable fluorinated building blocks.

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.

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.

Interplay of Ring Puckering and Hydrogen Bonding in Deoxyribonucleosides

The Cremer-Pople ring puckering analysis and the Konkoli-Cremer local mode analysis supported by the topological analysis of the electron density were applied for the first comprehensive analysis of the interplay between deoxyribose ring puckering and intramolecular H-bonding in 2'-deoxycytidine, 2'-deoxyadenosine, 2'-deoxythymidine, and 2'-deoxyguanosine. We mapped for each deoxyribonucleoside the complete conformational energy surface and the corresponding pseudorotation path. We found only incomplete pseudorotation cycles, caused by ring inversion, which we coined as pseudolibration paths. On each pseudolibration path a global and a local minimum separated by a transition state were identified. The investigation of H-bond free deoxyribonucleoside analogs revealed that removal of the H-bond does not restore the full conformational flexibility of the sugar ring. Our work showed that ring puckering predominantly determines the conformational energy; the larger the puckering amplitude, the lower the conformational energy. In contrast no direct correlation between conformational energy and H-bond strength was found. The longest and weakest H-bonds are located in the local minimum region, whereas the shortest and strongest H-bonds are located outside the global and local minimum regions at the turning points of the pseudolibration paths, i.e., H-bonding determines the shape and length of the pseudolibration paths. In addition to the H-bond strength, we evaluated the covalent/electrostatic character of the H-bonds applying the Cremer-Kraka criterion of covalent bonding. H-bonding in the puric bases has a more covalent character whereas in the pyrimidic bases the H-bond character is more electrostatic. We investigated how the mutual orientation of the CH₂OH group and the base influences H-bond formation via two geometrical parameters describing the rotation of the substituents perpendicular to the sugar ring and their tilting relative to the ring center. According to our results, rotation is more important for H-bond formation. In addition we assessed the influence of the H-bond acceptor, the lone pair (N, respectively O), via the delocalization energy. We found larger delocalization energies corresponding to stronger H-bonds for the puric bases. The global minimum conformation of 2'-deoxyguanosine has the strongest H-bond of all conformers investigated in this work with a bond strength of 0.436 which is even stronger than the H-bond in the water dimer (0.360). The application of our new analysis to DNA deoxyribonucleotides and to unnatural base pairs, which have recently drawn a lot of attention, is in progress.

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.