"Hydridic Hydrogen-Bond Donors" Are Not Hydrogen-Bond Donors.

Herein, we dismiss a recent proposal by Civis, Hobza, and co-workers to modify the IUPAC definition of hydrogen bonds in order to expand the scope from protonic Y-Hδ+ to hydridic Y-Hδ- hydrogen-bond donor fragments [J. Am. Chem. Soc. 2023, 145, 8550]. Based on accurate Kohn-Sham molecular orbital (KS-MO) analyses, we falsify the conclusion that interactions involving protonic and hydridic hydrogens are both hydrogen bonds; they are not. Instead, our quantitative KS-MO, energy decomposition, and Voronoi deformation density analyses reveal two fundamentally different bonding mechanisms for protonic Y-Hδ+ and hydridic Y-Hδ- fragments which go with charge transfer in opposite directions. On one hand, we confirm the IUPAC definition for regular hydrogen bonds in the case of protonic Y-Hδ+ fragments. On the other hand, complexes involving Y-Hδ- fragments are, in fact, acceptors in other well-known families of Lewis-acid/base interactions, such as halogen bonds, chalcogen bonds, and pnictogen bonds. These mechanisms lead to the same spectroscopic phenomenon in both the Y-Hδ+ and Y-Hδ- fragments, that is, the redshift in the Y-H stretching frequency, which is, thus, not an exclusive indicator for hydrogen bonding.

Chalcogen Atom Size: A Key Parameter in Modulating Carbonyl Compound Properties.

Exchanging oxygen in the functional group C=O (i. e., carbonyl) for the less electronegative Group 16 elements, sulfur or selenium, unexpectedly enhances the electronegativity of the C=X group in π-conjugated molecules and reduces the molecular π HOMO-LUMO energy gap. Quantum-chemical analyses revealed that the steric size of the chalcogen atom X is at the origin of this seemingly counterintuitive behavior. This tuning of the chemical properties of carbonyl compounds by varying the chalcogen atom size in the C=X bond can be applied in many fields of chemistry. This concept article delineates several useful applications in the fields of organocatalysis, supramolecular chemistry, and photo(electro)chemistry.

More Electropositive is More Electronegative: Atom Size Determines C=X Group Electronegativity.

Opposite to what one might expect, we find that the C=X group can become effectively more, not less, electronegative when the Pauling electronegativity of atom X decreases down Groups 16, 15, and 14 of the Periodic Table. Our quantum-chemical analyses, show that, and why, this phenomenon is a direct consequence of the increasing size of atom X down a group. These findings can be applied to tuning and improving the hydrogen-bond donor strength of amides H2 NC(=X)R by increasingly withdrawing density from the NH2 group. A striking example is that H2 NC(=SiR2 )R is a stronger hydrogen-bond donor than H2 NC(=CR2 )R.

Blueshift in Trifurcated Hydrogen Bonds: A Tradeoff between Tetrel Bonding and Steric Repulsion.

We have quantum chemically investigated the origin of the atypical blueshift of the H-C bond stretching frequency in the hydrogen-bonded complex X- •••H3 C-Y (X, Y=F, Cl, Br, I), as compared to the corresponding redshift occurring in Cl- •••H3 N and Cl- •••H3 C-H, using relativistic density functional theory (DFT) at ZORA-BLYP-D3(BJ)/QZ4P. Previously, this blueshift was attributed, among others, to the contraction of the H-C bonds as the H3 C moiety becomes less pyramidal. Herein, we provide quantitative evidence that, instead, the blueshift arises from a direct and strong X- •••C interaction of the HOMO of A- with the backside lobe on carbon of the low-lying C-Y antibonding σ* LUMO of the H3 C-Y fragment. This X- •••C bond, in essence a tetrel bond, pushes the H atoms towards a shorter H-C distance and makes the H3 C moiety more planar. The blueshift may, therefore, serve as a diagnostic for tetrel bonding.

Blueshift in Trifurcated Hydrogen Bonds: A Tradeoff between Tetrel Bonding and Steric Repulsion.

The front cover artwork is provided by the TheoCheM group at the Vrije Universiteit Amsterdam. The image shows how, in X- •••H3 C-Y complexes, the Lewis base X- tetrel-binds to the central C while sterically pushing the H atoms towards C; hence, the compression and blueshift of the H-C bonds. Read the full text of the Research Article at 10.1002/cphc.202300480.

Multiple hydrogen-bonded dimers: are only the frontier atoms relevant?

Non-frontier atom exchanges in hydrogen-bonded aromatic dimers can induce significant interaction energy changes (up to 6.5 kcal mol-1). Our quantum-chemical analyses reveal that the relative hydrogen-bond strengths of N-edited guanine-cytosine base pair isosteres, which cannot be explained from the frontier atoms, follow from the charge accumulation in the monomers.

The nature of metallophilic interactions in closed-shell d8-d8 metal complexes.

We have quantum chemically analyzed the closed-shell d8-d8 metallophilic interaction in dimers of square planar [M(CO)2X2] complexes (M = Ni, Pd, Pt; X = Cl, Br, I) using dispersion-corrected density functional theory at ZORA-BLYP-D3(BJ)/TZ2P level of theory. Our purpose is to reveal the nature of the [X2(CO)2M]⋯[M(CO)2X2] bonding mechanism by analyzing trends upon variations in M and X. Our analyses reveal that the formation of the [M(CO)2X2]2 dimers is favored by an increasingly stabilizing electrostatic interaction when the M increases in size and by more stabilizing dispersion interactions promoted by the larger X. In addition, there is an overlooked covalent component stemming from metal-metal and ligand-ligand donor-acceptor interactions. Thus, at variance with the currently accepted picture, the d8-d8 metallophilicity is attractive, and the formation of [M(CO)2X2]2 dimers is not a purely dispersion-driven phenomenon.

Understanding chemistry with the symmetry-decomposed Voronoi deformation density charge analysis.

The symmetry-decomposed Voronoi deformation density (VDD) charge analysis is an insightful and robust computational tool to aid the understanding of chemical bonding throughout all fields of chemistry. This method quantifies the atomic charge flow associated with chemical-bond formation and enables decomposition of this charge flow into contributions of (1) orbital interaction types, that is, Pauli repulsive or bonding orbital interactions; (2) per irreducible representation (irrep) of any point-group symmetry of interacting closed-shell molecular fragments; and now also (3) interacting open-shell (i.e., radical) molecular fragments. The symmetry-decomposed VDD charge analysis augments the symmetry-decomposed energy decomposition analysis (EDA) so that the charge flow associated with Pauli repulsion and orbital interactions can be quantified both per atom and per irrep, for example, for σ, π, and δ electrons. This provides detailed insights into fundamental aspects of chemical bonding that are not accessible from EDA.

Chalcogen Atom Size Dictates Stability of Benzene-1,3,5-triamide Polymers: Overlooked Role of Geometrical Fit for Enhanced Hydrogen Bonding.

Our quantum chemical analyses elucidated how the replacement of O in the amide bonds of benzene-1,3-5-tricarboxamides (OBTAs) with the larger chalcogens S and Se enhances the intermolecular interactions and thereby the stability of the obtained hydrogen-bonded supramolecular polymers due to two unexpected reasons: i) the SBTA and SeBTA monomers have a better geometry for self-assembly and ii) induce stronger covalent (hydrogen-bond) interactions besides enhanced dispersion interactions. In addition, it is shown that the cooperativity in benzene-1,3,5-triamide (BTA) self-assembly is caused by charge separation in the σ-electronic system following the covalency in the hydrogen bonds.

Strengthened cooperativity of DNA-based cyclic hydrogen-bonded rosettes by subtle functionalization.

Cooperative effects cause extra stabilization of hydrogen-bonded supramolecular systems. In this work we have designed hydrogen-bonded rosettes derived from a guanine-cytosine Janus-type motif with the aim of finding a monomer that enhances the synergy of supramolecular systems. For this, relativistic dispersion-corrected density functional theory computations have been performed. Our proposal involves a monomer with three hydrogen-bonds pointing in the same direction, which translates into shorter bonds, stronger donor-acceptor interactions, and more attractive electrostatic interactions, thus giving rise to rosettes with strengthened cooperativity. This newly designed rosette has triple the cooperativity found for the naturally occurring guanine quadruplex.

How the Chalcogen Atom Size Dictates the Hydrogen-Bond Donor Capability of Carboxamides, Thioamides, and Selenoamides.

Invited for the cover of this issue are Celine Nieuwland and Célia Fonseca Guerra of the Vrije Universiteit Amsterdam. The image depicts how the increasing atom size of the chalcogen from O to S to Se elongates the carbon-chalcogen bond in amides due to the increase in the steric Pauli repulsion and thereby enhances the amide hydrogen-bond donor strength. Read the full text of the article at 10.1002/chem.202200755.

How the Chalcogen Atom Size Dictates the Hydrogen-Bond Donor Capability of Carboxamides, Thioamides, and Selenoamides.

The amino groups of thio- and selenoamides can act as stronger hydrogen-bond donors than of carboxamides, despite the lower electronegativity of S and Se. This phenomenon has been experimentally explored, particularly in organocatalysis, but a sound electronic explanation is lacking. Our quantum chemical investigations show that the NH2 groups in thio- and selenoamides are more positively charged than in carboxamides. This originates from the larger electronic density flow from the nitrogen lone pair of the NH2 group towards the lower-lying π*C=S and π*C=Se orbitals than to the high-lying π*C=O orbital. The relative energies of the π* orbitals result from the overlap between the chalcogen np and carbon 2p atomic orbitals, which is set by the carbon-chalcogen equilibrium distance, a consequence of the Pauli repulsion between the two bonded atoms. Thus, neither the electronegativity nor the often-suggested polarizability but the steric size of the chalcogen atom determines the amide's hydrogen-bond donor capability.

Iodine Gauche Effect Induced by an Intramolecular Hydrogen Bond.

The gauche conformer in 1-X,2-Y-disubstituted ethanes, that is, the staggered orientation in which X and Y are in closer contact, is only favored for relatively small substituents that do not give rise to large X···Y steric repulsion. For more diffuse substituents, weakly attractive orbital interactions between antiperiplanar bonds (i.e., hyperconjugation) cannot overrule the repulsive forces between X and Y. Our quantum chemical analyses of the rotational isomerism of XCH2CH2Y (X = F, OH; Y = I) at ZORA-BP86-D3(BJ)/QZ4P reveal that indeed the anti conformer is generally favored due to a less destabilizing I···F and I···O-H steric repulsion. The only case when the gauche conformer is preferred is when the hydroxyl hydrogen is oriented toward the iodine atom in the 2-iodoethanol. This is because of the significantly stabilizing covalent component of the I···H-O intramolecular hydrogen bond. Therefore, we show that strong intramolecular interactions can overcome the steric repulsion between bulky substituents in 1,2-disubstituted ethanes and cause the gauche effect. Our quantum chemical computations have guided nuclear magnetic resonance experiments that confirm the increase in the gauche population as X goes from F to OH.

Clarifying notes on the bonding analysis adopted by the energy decomposition analysis.

We discuss the fundamental aspects of the EDA-NOCV method and address some critical comments that have been made recently. The EDA-NOCV method unlike most other methods focuses on the process of bond formation between the interacting species and not just only on the analysis of the finally formed bond. This is demonstrated using LiF as an example. There is a difference between the interactions between the initial species which form the bond and are also the final product of bond cleavage, and the interactions between the fragments in the eventually formed molecule. The flexibility of the method allows the choice of the interacting fragments which helps to identify the charge and electron configuration of the fragments which describe the bond. This is very helpful in cases where the bond may be described with several Lewis structures. We reject the idea that it would be a disadvantage to have "bond path functions" as the energy components in the EDA, which actually indicate the variability of the method. The bonding analysis in a different sequence of the bond formation gives important results for the various questions that can be asked. This is demonstrated by using CH2, CO2 and the formation of a guanine quartet as examples. The fact that a bond is always defined by the bound molecule, the fragments, and their states is universal and deeply physical, as we show here again for various examples. The results of the EDA-NOCV method are in full accordance with the physical mechanism of the chemical bond as revealed by Ruedenberg.

Computational Prediction of Tc-99 NMR Chemical Shifts in Technetium Complexes with Radiopharmaceutical Applications.

The Tc-99m nucleus is the most used nuclide in radiopharmaceuticals designed for imaging diagnosis. The metal can exist in nine distinct oxidation states and forms distinct coordination complexes with a variety of chelating agents and geometries. These complexes are usually characterized through Tc-99 NMR that is very sensitive to the Tc coordination sphere. Therefore, predicting Tc-99 NMR might be useful to assist experimentalists in structural characterization. In the present study, we propose three computational protocols for predicting Tc-99 NMR chemical shifts based on density functional theory calculations using relativistic and nonrelativistic Hamiltonians: the relativistic Model 1, the nonrelativistic Model 2, and the empirical nonrelativistic Model 3. In Models 2 and 3, the NMR-DKH basis set was used for all atoms, including the Tc, for which it was developed here. All models were applied for a set of 41 Tc-complexes with metal oxidation states 0, I, and V, for which the Tc-99 chemical shift was available experimentally. The mean absolute deviation and the mean relative deviation were 67 ppm and 4.8% (Model 1), 92 ppm and 6.2% (Model 2), and 65 ppm and 4.9% (Model 3), respectively. Last, the effect of the explicit solvent was evaluated for the [TcO2(en)2]+─Tc(V) complex. The calculated results for the Tc-99 NMR chemical shift at SO-ZORA-SSB-D/TZ2P-ZORA/COSMO//TPSS/def2-SVP/IEF-PCM(UFF) show that the inclusion of 14 water molecules (first solvation shell) together with the implicit solvation model leads to an absolute deviation of only 7 ppm (0.3%) from the experimental value, indicating that the solvent effects play a key role in predicting Tc-99 NMR.

Not Carbon s-p Hybridization, but Coordination Number Determines C-H and C-C Bond Length.

A fundamental and ubiquitous phenomenon in chemistry is the contraction of both C-H and C-C bonds as the carbon atoms involved vary, in s-p hybridization, along sp3 to sp2 to sp. Our quantum chemical bonding analyses based on Kohn-Sham molecular orbital theory show that the generally accepted rationale behind this trend is incorrect. Inspection of the molecular orbitals and their corresponding orbital overlaps reveals that the above-mentioned shortening in C-H and C-C bonds is not determined by an increasing amount of s-character at the carbon atom in these bonds. Instead, we establish that this structural trend is caused by a diminishing steric (Pauli) repulsion between substituents around the pertinent carbon atom, as the coordination number decreases along sp3 to sp2 to sp.

How Divalent Cations Interact with the Internal Channel Site of Guanine Quadruplexes.

The front cover artwork is provided by the TheoCheM group of the Vrije Universiteit Amsterdam. The image shows that guanine quadruplexes have a stronger binding affinity for divalent cations than monovalent cations. Read the full text of the Article at 10.1002/cphc.202100529.

Boron Tunneling in the "Weak" Bond-Stretch Isomerization of N-B Lewis Adducts.

Some nitrile-boron halide adducts exhibit a double-well potential energy surface with two distinct minima: a "long bond" geometry (LB, a van der Waals interaction mostly based on electrostatics, but including a residual charge transfer component) and a "short bond" structure (SB, a covalent dative bond). This behavior can be considered as a "weak" form of bond stretch isomerism. Our computations reveal that complexes RCN-BX3 (R=CH3 , FCH2 , BrCH2 , and X=Cl, Br) exhibit a fast interconversion from LB to SB geometries even close to the absolute zero thanks to a boron atom tunneling mechanism. The computed half-lives of the meta-stable LB compounds vary between minutes to nanoseconds at cryogenic conditions. Accordingly, we predict that the long bond structures are practically impossible to isolate or characterize, which agrees with previous matrix-isolation experiments.

How Divalent Cations Interact with the Internal Channel Site of Guanine Quadruplexes.

The formation of guanine quadruplexes (GQ) in DNA is crucial in telomere homeostasis and regulation of gene expression. Pollution metals can interfere with these DNA superstructures upon coordination. In this work, we study the affinity of the internal GQ channel site towards alkaline earth metal (Mg2+ , Ca2+ , Sr2+ , and Ba2+ ), and (post-)transition metal (Zn2+ , Cd2+ , Hg2+ , and Pb2+ ) cations using density functional theory computations. We find that divalent cations generally bind to the GQ cavity with a higher affinity than conventional monovalent cations (e. g. K+ ). Importantly, we establish the nature of the cation-GQ interaction and highlight the relationship between ionic and nuclear charge, and the electrostatic and covalent interactions. The covalent interaction strength plays an important role in the cation affinity and can be traced back to the relative stabilization of cations' unoccupied atomic orbitals. Overall, our findings contribute to a deeper understanding of how pollution metals could induce genomic instability.

The Gauche Effect in XCH2 CH2 X Revisited.

We have quantum chemically investigated the rotational isomerism of 1,2-dihaloethanes XCH2 CH2 X (X = F, Cl, Br, I) at ZORA-BP86-D3(BJ)/QZ4P. Our Kohn-Sham molecular orbital (KS-MO) analyses reveal that hyperconjugative orbital interactions favor the gauche conformation in all cases (X = F-I), not only for X = F as in the current model of this so-called gauche effect. We show that, instead, it is the interplay of hyperconjugation with Pauli repulsion between lone-pair-type orbitals on the halogen substituents that constitutes the causal mechanism for the gauche effect. Thus, only in the case of the relatively small fluorine atoms, steric Pauli repulsion is too weak to overrule the gauche preference of the hyperconjugative orbital interactions. For the larger halogens, X⋅⋅⋅X steric Pauli repulsion becomes sufficiently destabilizing to shift the energetic preference from gauche to anti, despite the opposite preference of hyperconjugation.