A comparative study of some lithium and hydrogen‐bonded complexes: Ab initio and QTAIM studies

The nature of the interactions of cyanide with lithium and hydrogen halides was investigated using ab initio calculations and topological analysis of electron density. The computed properties of the lithium‐bonded complexes RCN···LiX (R = H, F, Cl, Br, C?CH, CH?CH₂, CH₃, C₂H₅; X = Cl, Br) were compared with those of corresponding hydrogen‐bonded complexes RCN···HX. The results show that both types of intermolecular interactions are “closed‐shell” noncovalent interactions. The effect of substitution on the interaction energy and electron density at the bond critical points of the lithium and hydrogen bonding interactions is similar. In comparison, the interaction energies of lithium‐bonded complexes are more negative than those of hydrogen‐bonded counterparts. The electrostatic interaction plays a more important role in the lithium bond than in the hydrogen bond. On complex formation, the net charge and energy of the Li atom decrease and the atomic volume increases, while the net charge and energy of the H atom increase and the atomic volume decreases. © 2013 Wiley Periodicals, Inc.     

A polarizable dipole–dipole interaction model for evaluation of the interaction energies for NH···OC and CH···OC hydrogen‐bonded complexes

In this article, a polarizable dipole–dipole interaction model is established to estimate the equilibrium hydrogen bond distances and the interaction energies for hydrogen‐bonded complexes containing peptide amides and nucleic acid bases. We regard the chemical bonds N? H, C?O, and C? H as bond dipoles. The magnitude of the bond dipole moment varies according to its environment. We apply this polarizable dipole–dipole interaction model to a series of hydrogen‐bonded complexes containing the N? H···O?C and C? H···O?C hydrogen bonds, such as simple amide‐amide dimers, base‐base dimers, peptide‐base dimers, and β‐sheet models. We find that a simple two‐term function, only containing the permanent dipole–dipole interactions and the van der Waals interactions, can produce the equilibrium hydrogen bond distances compared favorably with those produced by the MP2/6‐31G(d) method, whereas the high‐quality counterpoise‐corrected (CP‐corrected) MP2/aug‐cc‐pVTZ interaction energies for the hydrogen‐bonded complexes can be well‐reproduced by a four‐term function which involves the permanent dipole–dipole interactions, the van der Waals interactions, the polarization contributions, and a corrected term. Based on the calculation results obtained from this polarizable dipole–dipole interaction model, the natures of the hydrogen bonding interactions in these hydrogen‐bonded complexes are further discussed. © 2013 Wiley Periodicals, Inc.     

The Double‐Faced Nature of Hydrogen Bonding in Hydroxy‐Functionalized Ionic Liquids Shown by Neutron Diffraction and Molecular Dynamics Simulations

We characterize the double‐faced nature of hydrogen bonding in hydroxy‐functionalized ionic liquids by means of neutron diffraction with isotopic substitution (NDIS), molecular dynamics (MD) simulations, and quantum chemical calculations. NDIS data are fit using the empirical potential structure refinement technique (EPSR) to elucidate the nearest neighbor H???O and O???O pair distribution functions for hydrogen bonds between ions of opposite charge and the same charge. Despite the presence of repulsive Coulomb forces, the cation–cation interaction is stronger than the cation–anion interaction. We compare the hydrogen‐bond geometries of both “doubly charged hydrogen bonds” with those reported for molecular liquids, such as water and alcohols. In combination, the NDIS measurements and MD simulations reveal the subtle balance between the two types of hydrogen bonds: The small transition enthalpy suggests that the elusive like‐charge attraction is almost competitive with conventional ion‐pair formation.     

Hydrogen Bonding Between Ions of Like Charge in Ionic Liquids Characterized by NMR Deuteron Quadrupole Coupling Constants—Comparison with Salt Bridges and Molecular Systems

We present deuteron quadrupole coupling constants (DQCC) for hydroxyl‐functionalized ionic liquids (ILs) in the crystalline or glassy states characterizing two types of hydrogen bonding: The regular Coulomb‐enhanced hydrogen bonds between cation and anion (c–a), and the unusual hydrogen bonds between cation and cation (c–c), which are present despite repulsive Coulomb forces. We measure these sensitive probes of hydrogen bonding by means of solid‐state NMR spectroscopy. The DQCCs of (c–a) ion pairs and (c–c) H‐bonds are compared to those of salt bridges in supramolecular complexes and those present in molecular liquids. At low temperatures, the (c–c) species successfully compete with the (c–a) ion pairs and dominate the cluster populations. Equilibrium constants obtained from molecular‐dynamics (MD) simulations show van't Hoff behavior with small transition enthalpies between the differently H‐bonded species. We show that cationic‐cluster formation prevents these ILs from crystallizing. With cooling, the (c–c) hydrogen bonds persist, resulting in supercooling and glass formation.     

Polyion Covalency: Exotic Species from the Unexplored World of Electrostatically Shielded Molecular Ion Chemistry

Standard quantum chemical methods have been employed to describe a variety of kinetically stable polyionic molecular species that are trapped in appreciable potential wells by chemical bonding forces, despite powerful electrostatic opposition that challenges conventional chemical detection and characterization. The studied species are covalent or dative analogs of “anti‐electrostatic” hydrogen‐bonded (AEHB) species, all illustrating how short‐range quantum covalency can overcome the powerful “shielding” opposition of long‐range electrostatic forces to form highly charged molecular species, analogous to known neutral or singly ionic counterparts. Computational predictions of representative structural, spectroscopic, and NBO‐based electronic signatures of multiply charged analogs of common neutral species (CH₃CH₃, CO₂, FeCO) are provided to suggest the unique material properties characteristic of this shielded domain of polyionic chemical phenomena.     

“Anti-electrostatic” Halogen Bonding between Ions of Like Charge

Halogen bonding occurs between molecules featuring Lewis acidic halogen substituents and Lewis bases. It is often rationalized as a predominantly electrostatic interaction and thus interactions between ions of like charge (e. g., of anionic halogen bond donors with halides) seem counter-intuitive. Herein, we provide an overview on such complexes. First, theoretical studies are described and their findings are compared. Next, experimental evidences are presented in the form of crystal structure database analyses, recent examples of strong “anti-electrostatic” halogen bonding in crystals, and the observation of such interactions also in solution. We then compare these complexes to select examples of “counter-intuitive” adducts formed by other interactions, like hydrogen bonding. Finally, we comment on key differences between charge-transfer and electrostatic polarization.     

Theoretical Study on Cooperativity Effects between Anion–π and Halogen‐Bonding Interactions

This article analyzes the interplay between lone pair–π (lp–π) or anion–π interactions and halogen‐bonding interactions. Interesting cooperativity effects are observed when lp/anion–π and halogen‐bonding interactions coexist in the same complex, and they are found even in systems in which the distance between the anion and halogen‐bond donor molecule is longer than 9 Å. These effects are studied theoretically in terms of energetic and geometric features of the complexes, which are computed by ab initio methods. Bader′s theory of “atoms in molecules” is used to characterize the interactions and to analyze their strengthening or weakening depending upon the variation of charge density at critical points. The physical nature of the interactions and cooperativity effects are studied by means of molecular interaction potential with polarization partition scheme. By taking advantage of all aforementioned computational methods, the present study examines how these interactions mutually influence each other. Additionally, experimental evidence for such interactions is obtained from the Cambridge Structural Database (CSD).     

Theoretical Studies on the Dihydrogen Bonding Between Shortchain Hydrocarbon and Magnesium Hydride

The C--H…H dihydrogen-bonded complexes of methane, ethylene, acetylene, and their derivatives with magnesium hydride were systematically investigated at MP2/aug-cc-PVTZ level. The results confirm that the strength of dihydrogen bonding increases in the following order of proton donors： C（sp3）-H〈C（sp2）-H〈C（sp）-H and chlorine substituents enhance the C-H…H interaction. In the majority of the complexes with a cyclic structure, the Mg-H proton-accepting bond is more sensitive to the surroundings than C-H proton-donating bond. The nature of the electrostatic interaction in these C-H…H dihydrogen bonds was also unveiled by means of the atoms in mo- lecules（AIM） analysis. The natural bond orbital（NBO） analysis suggests that the charge transfer in the cyclic com- plexes is characteristic of dual-channel. The direction of the net charge transfer in the cyclic complexes is contrary to that previously found in dihydrogen bonded systems.     

Density functional study of indole–pyrrole heterodimer potential energy hypersurface

A density‐functional study of indole–pyrrole heterodimer potential energy hypersurface (PES) was performed. Eight stationary points were located on the B3LYP/6‐31++G(d,p) PES, three of which correspond to real minima, all of them being characterized with an N? H … π type hydrogen bonding. In two of these minima (the local ones), pyrrole subunit acts as a hydrogen bond proton donor, while the global minimum corresponds to indole–H … π(‐pyrrole) arrangement. Besides the interaction and dissociation energies corrected for BSSE and the monomer relaxation energies and the relevant structural parameters, anharmonic N? H and N? H … π vibrational frequencies were calculated for various N? H oscillators involved in this interaction from the 1‐D DFT vibrational potentials. On the basis of anharmonic vibrational frequency analysis, it was concluded that the two types of N? H … π hydrogen bonded dimers (indole vs. pyrrole being a proton donor) should be distinguishable with spectroscopic methods. Various contributions to the overall anharmonic frequency shifts upon hydrogen bonding were calculated and discussed as well. The charge field perturbation (CFP) technique was employed to study the electrostatic + polarization influence of the proton accepting unit on the N? H(… π) vibrational potential. The second‐order perturbation theory analysis (SOPT) of the Fock matrix (i.e., its Kohn–Sham analog) within the natural bond orbital (NBO) basis, as well as various NBO deletion analyses revealed an essentially one‐directional charge transfer (CT) of a π(C? C) → σ*(N? H) type in the case of all three minima. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Atomic dipole approximation and energies of interactions between purine and pyrimidine bases. I. Electrostatic interactions of adenine with uracil,thymine, thiouracils,dihydrouracil, and 5-fluorouracil

The atomic dipole approximation has been employed to calculate the energies of electrostatic interaction between adenine and various pyrimidine derivatives. Minima of the interaction energy for various planar configurations were determined. Inclusion of the “monopole–dipole” and “dipole–dipole” terms in the multipole expansion improves considerably the agreement with experimental data. The effect of sulfur substitution has been investigated in detail. Formation of N? H…?S hydrogen bonds is less favorable than of N? H…?O bonds, due largely to the lower atomic dipole of the sulfur atom resulting from the shift of the π-electron charge toward the neighboring carbon. The results are relevant to the interactions of thiouracils in nucleic acids.     

Attraction between electrophilic caps: A counterintuitive case of noncovalent interactions

Intermolecular attractive interaction between electrophilic sites is a counterintuitive phenomenon, as the electrostatic interaction therein is repulsive and destabilizing. Here, we confirm this phenomenon in four representative complexes, using state-of-the-art quantum mechanical methods. By employing the block-localized wavefunction (BLW) method, which can turn off intermolecular charge transfer interactions, we profoundly demonstrated the significance of charge transfer interactions in these seemingly counterintuitive complexes. Indeed, after being “turned off” the intermolecular charge transfer interaction in, for example, the FNSi···BrF complex, the originally attractive intermolecular interaction turns to be repulsive. The energy decomposition approach based on the BLW method (BLW-ED) can partition the overall stability gained on the formation of intermolecular noncovalent interaction into several physically meaningful components. According to the BLW-ED analysis, the electrostatic repulsion in these counterintuitive cases is overwhelmed by the stabilizing polarization, dispersion interaction, and most importantly, the charge transfer interaction, resulting in the eventual counterintuitive overall attraction. The present study suggests that, predicting bonding sites of noncovalent interactions using only the “hole” concept may be not universally sufficient, because other significant stabilizing factors will contribute to the stability and sometimes, play even bigger roles than the electrostatic interaction and consequently govern the complex structures. © 2018 Wiley Periodicals, Inc.     

“Anti‐Electrostatic” Halogen Bonding

Halogen bonding is often described as being driven predominantly by electrostatics, and thus adducts between anionic halogen bond (XB) donors (halogen‐based Lewis acids) and anions seem counterintuitive. Such “anti‐electrostatic” XBs have been predicted theoretically but for organic XB donors, there are currently no experimental examples except for a few cases of self‐association. Reported herein is the synthesis of two negatively charged organoiodine derivatives that form anti‐electrostatic XBs with anions. Even though the electrostatic potential is universally negative across the surface of both compounds, DFT calculations indicate kinetic stabilization of their halide complexes in the gas phase and particularly in solution. Experimentally, self‐association of the anionic XB donors was observed in solid‐state structures, resulting in dimers, trimers, and infinite chains. In addition, co‐crystals with halides were obtained, representing the first cases of halogen bonding between an organic anionic XB donor and a different anion. The bond lengths of all observed interactions are 14–21 % shorter than the sum of the van der Waals radii.     

Structure and interaction between the [BMIM][Ala] alanine anion and the 1-butyl-3-methylimidazolium cation in ion pairs

The equilibrium geometries and vibrational frequencies of the ionic liquid 1-butyl-3-methylimidazolium cation and the alanine anion [BMIM][Ala] are studied using density functional theory (DFT) at the B3PW91/6-311+G(d,p) leve1. The most stable structures of the anion, the cation, and the ion pairs are obtained and characterized, and the geometry parameters of the ion pairs confirm the presence of a hydrogen bonding interaction between the anion and the cation. Natural bond orbital (NBO) analysis is also performed to analyze the atomic charge distribution and charge transfer in the [BMIM]⁺ cation and [BMIM][Ala] ionic liquids. The results show that there are the electrostatic interaction and multiple hydrogen bond interactions between the cation and the anion of the ionic liquids, and the stability of the ground state of the ion pairs mostly results from the hydrogen bonding between the lone pairs of O atoms in the anion and H in the imidazole cycle of the cation. There are some changes in microstructures and the charge distribution during the formation of the ion pairs.     

The Quadruple Bonding in C2 Reproduces the Properties of the Molecule

Ever since Lewis depicted the triple bond for acetylene, triple bonding has been considered as the highest limit of multiple bonding for main elements. Here we show that C₂ is bonded by a quadruple bond that can be distinctly characterized by valence‐bond (VB) calculations. We demonstrate that the quadruply‐bonded structure determines the key observables of the molecule, and accounts by itself for about 90 % of the molecule's bond dissociation energy, and for its bond lengths and its force constant. The quadruply‐bonded structure is made of two strong π bonds, one strong σ bond and a weaker fourth σ‐type bond, the bond strength of which is estimated as 17–21 kcal mol^?1. Alternative VB structures with double bonds; either two π bonds or one π bond and one σ bond lie at 129.5 and 106.1 kcal mol^?1, respectively, above the quadruply‐bonded structure, and they collapse to the latter structure given freedom to improve their double bonding by dative σ bonding. The usefulness of the quadruply‐bonded model is underscored by “predicting” the properties of the ³ state. C₂’s very high reactivity is rooted in its fourth weak bond. Thus, carbon and first‐row main elements are open to quadruple bonding!     

Electronic and topological properties of interactions between imidazolium-based ionic liquids and thiophenic compounds: a theoretical investigation

To deepen the understanding the interactions of thiophenic compounds in ionic liquids, we have performed a systemic study on the electronic structures, and topological properties of interactions between N-ethyl-N-ethylimidazolium diethyl phosphate ([EEIM][DEP]) ionic liquid and 3-methylthiophene (3-MT), benzothiophene (BT), or dibenzothiophene (DBT) using density functional theory. From NBO atomic charges and electrostatic potential analyses, most of the positive charge is located on C2–H2 in the [EEIM] cation, and the negative charge is focused on oxygen atoms in [DEP] anion, implying oxygen atoms in [DEP] should easily attack C2–H2 in [EEIM]. The electrostatic interaction between anion and cation may be dominant for the formation of the [EEIM]–[DEP] ion pair. The large stabilizing effect is due to the strong orbital interactions between the antibonding orbital of proton donor σ*(C2–H2) in [EEIM] cation and the lone pairs of proton acceptor LP(O) in [DEP] anion. A common feature of [EEIM][DEP], [EEIM][DEP]-3-MT/BT/DBT complexes is the presence of hydrogen bonds between [EEIM] cation and [DEP] anion. This work has also given the interacting mechanism of 3-MT, BT, and DBT adsorption on [EEIM][DEP] ionic liquid. Both [EEIM] cation and [DEP] anion are shown to play important roles in interactions between 3-MT, BT, DBT and [EEIM][DEP], which has been corroborated by NBO and AIM analyses. The π···π, π···C–H and hydrogen bonding interactions occur between [EEIM][DEP] and 3-MT, BT, DBT. The strength of sulfur involved interactions between 3-MT, BT, DBT and [EEIM][DEP] follows the order of 3-MT > BT > DBT. The order of interaction energies between [EEIM][DEP] and 3-MT, BT, DBT is 3-MT < BT < DBT, in agreement with the order of extractive selectivity from fuel oils (DBT > BT > 3-MT) in terms of sulfur partition coefficients.     

Crystallographic characterization of the first reported crystalline form of the potent hallucinogen (R)‐2‐amino‐1‐(8‐bromobenzo[1,2‐b;5,4‐b′]difuran‐4‐yl)propane or `bromodragonfly': the 1:1 anhydrous proton‐transfer compound with 3,5‐dinitrosalicylic acid

The 1:1 proton‐transfer compound of the potent substituted amphetamine hallucinogen (R)‐2‐amino‐1‐(8‐bromobenzo[1,2‐b;5,4‐b′]difuran‐4‐yl)propane (common trivial name `bromodragonfly') with 3,5‐dinitrosalicylic acid, namely 1‐(8‐bromobenzo[1,2‐b;5,4‐b′]difuran‐4‐yl)propan‐2‐aminium 2‐carboxy‐4,6‐dinitrophenolate, C₁₃H₁₃BrNO₂⁺·C₇H₃N₂O₇⁻, forms hydrogen‐bonded cation–anion chain substructures comprising undulating head‐to‐tail anion chains formed through C(8) carboxyl–nitro O—H...O associations and incorporating the aminium groups of the cations. The intrachain cation–anion hydrogen‐bonding associations feature proximal cyclic R₃³(8) interactions involving both an N⁺—H...O_phenolate and the carboxyl–nitro O—H...O associations and aromatic π–π ring interactions [minimum ring centroid separation = 3.566 (2) Å]. A lateral hydrogen‐bonding interaction between the third aminium H atom and a carboxyl O‐atom acceptor links the chain substructures, giving a two‐dimensional sheet structure. This determination represents the first of any form of this compound and is in the (R) absolute configuration. The atypical crystal stability is attributed both to the hydrogen‐bonded chain substructures provided by the anions, which accommodate the aminium proton‐donor groups of the cations and give crosslinking, and to the presence of the cation–anion aromatic ring π–π interactions.     

A new unconventional halogen bond C?X···H?M between HCCX (X = Cl and Br) and HMH (M = Be and Mg): An ab initio study

In this article, a new type of halogen‐bonded complex YCCX···HMY (X = Cl, Br; M = Be, Mg; Y = H, F, CH₃) has been predicted and characterized at the MP2/aug‐cc‐pVTZ level. We named it as halogen‐hydride halogen bonding. In each YCCX···HMY complex, a halogen bond is formed between the positively charged X atom and the negatively charged H atom. This new kind of halogen bond has similar characteristics to the conventional halogen bond, such as the elongation of the C? X bond and the red shift of the C? X stretch frequency upon complexation. The interaction strength of this type of halogen bond is in a range of 3.34–10.52 kJ/mol, which is smaller than that of dihydrogen bond and conventional halogen bond. The nature of the electrostatic interaction in this type of halogen bond has also been unveiled by means of the natural bond orbital, atoms in molecules, and energy decomposition analyses. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Tuning anion‐π interaction via halogen substituent effects in cyanuric acids and its derivatives

Several anion‐π complexes of isocyanuric acid, thioderivatives and their halogen substituted derivatives with chloride anion have been studied. The geometric and energetic features, charges transfer from chloride anion to the aromatic rings and “atoms‐in‐molecules” analysis are performed and discussed for these complexes. The results show that the strength of the anion‐π interaction between cyanuric derivatives and chloride anion can be tuned by halogen‐substituting. The localized molecular orbital energy decomposition analyses shows that, in the total interaction, exchange and electrostatic energies are the dominant stabilizing forces, and the polarization energies also make a favorable contribution. Finally, solvent effect significantly weakens the anion‐π interaction between the isocyanuric derivatives and chloride anion, especially in polar solvents. © 2015 Wiley Periodicals, Inc.     

Thermal dissociation of 1,2-dioxethane. III. Localized molecular-orbital study

The localized MO 's (LMO 's) of 1,2-dioxethane in its ground state and along the dissociation reaction path (to formaldehyde products) are generated using Boys' criteria for localization. The total charge density in each LMO is partitioned into atomic and overlap densities and the binding or antibinding character of each LMO is discussed in terms of the forces exerted on the nuclei by these densities. The driving force for the dissociation reaction is shown to arise essentially from the atomic dipole forces exerted on the oxygen nuclei by their lone-pair LMO 's. The characterization of a saddle point on the potential energy surface has been discussed in terms of the electrostatic equilibrium between forces exerted by the electron clouds “incomplete following” and “preceding.” The differences between the LMO 's obtained from the two Hartree–Fock solutions to which the SCF procedure converges have been discussed.