Anionic CH⋅⋅⋅X− Hydrogen Bonds: Origin of Their Strength,Geometry, and Other Properties

CF₃H as a proton donor was paired with a variety of anions, and its properties were assessed by MP2/aug‐cc‐pVDZ calculations. The binding energy of monoanions halide, NO₃^?, formate, acetate, HSO₄^?, and H₂PO₄^? lie in the 12–17 kcal mol^?1 range, although F^? is more strongly bound, by 26 kcal mol^?1. Dianions SO₄^2? and HPO₄^2? are bound by 27 kcal mol^?1, and trianion PO₄^3? by 45 kcal mol^?1. When two O atoms are available on the anion, the CH???O^? H‐bond (HB) is usually bifurcated, although asymmetrically. The CH bond is elongated and its stretching frequency redshifted in these ionic HBs, but the shift is reduced in the bifurcated structures. Slightly more than half of the binding energy is attributed to Coulombic attraction, with smaller contributions from induction and dispersion. The amount of charge transfer from the anions to the σ*(CH) orbital correlates with many of the other indicators of bond strength, such as binding energy, CH bond stretch, CH redshift, downfield NMR spectroscopic chemical shift of the bridging proton, and density at bond critical points.     

Combined experimental and computational modeling studies on 4‐[(2‐hydroxy‐3‐methylbenzylidene) amino]‐1,5‐dimethyl‐2‐phenyl‐1,2‐dihydro‐3H‐pyrazol‐3‐one

The Schiff base compound, 4‐[(2‐hydroxy‐3‐methylbenzylidene)amino]‐1,5‐dimethyl‐2‐phenyl‐1,2‐dihydro‐3H‐pyrazol‐3‐one, has been synthesized and characterized by IR, UV–vis, and X‐ray single‐crystal determination. Molecular geometry from X‐ray experiment of the title compound in the ground state have been compared using the density functional method (B3LYP) with 6‐31G(d,p) basis set. Calculated results show that density functional theory (DFT) can well reproduce the structure of the title compound. The energetic behavior of the title compound in solvent media has been examined using B3LYP method with the 6‐31G(d,p) basis set by applying the Onsager and the polarizable continuum model (PCM). The results obtained with these methods reveal that the PCM method provided more stable structure than Onsager's method. By using TD‐DFT method, electronic absorption spectra of the title compound have been predicted and a good agreement with the TD‐DFT method and the experimental one is determined. The predicted nonlinear optical properties of the title compound are much greater than ones of urea. In addition, DFT calculations of the title compound, molecular electrostatic potential and NBO analysis were performed at B3LYP/6‐31G(d,p) level of theory. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

H‐bonding and stacking interactions between chloroquine and temozolomide

The interactions between temozolomide and chloroquine were examined via Dispersion‐Corrected Density Functional Theory and MP2 methods. Chloroquine was considered in both its lowest energy structure and in a local minimum where its aromatic system and secondary amine group are free to interact directly with temozolomide. The accessibility of these two components to intermolecular interaction makes the lowest energy dimer of this local monomer minimum competitive in total energy with that involving chloroquine's most stable monomer geometry. In either case, the most stable heterodimer places the aromatic ring systems of the two molecules parallel and directly above one another in a stacked geometry. Most of the local minima are also characterized by a stacked geometry as well. Comparison between B3LYP and B3LYP‐D binding energies confirms dispersion is a primary factor in stabilizing these structures. © 2016 Wiley Periodicals, Inc.     

Comparison between Hydrogen and Halogen Bonds in Complexes of 6-OX-Fulvene with Pnicogen and Chalcogen Electron Donors

Quantum chemical calculations are applied to complexes of 6-OX-fulvene (X=H, Cl, Br, I) with ZH₃/H₂Y (Z=N, P, As, Sb; Y=O, S, Se, Te) to study the competition between the hydrogen bond and the halogen bond. The H-bond weakens as the base atom grows in size and the associated negative electrostatic potential on the Lewis base atom diminishes. The pattern for the halogen bonds is more complicated. In most cases, the halogen bond is stronger for the heavier halogen atom, and pnicogen electron donors are more strongly bound than chalcogen. Halogen bonds to chalcogen atoms strengthen in the order O<S<Se<Te, whereas the pattern is murkier for the pnicogen donors. In terms of competition, most halogen bonds to pnicogen donors are stronger than their H-bond analogues, but there is no clear pattern with respect to chalcogen donors. O prefers a H-bond, while halogen bonds are favored by Te. For S and Se, I-bonds are strongest, followed Br, H, and Cl-bonds in that order.     

Predicting reduction potentials of 1,3,6‐triphenyl fulvenes using molecular electrostatic potential analysis of substituent effects

The influence of mono‐ and multiple substituent effect on the reduction potential (E⁰) of 1,3,6‐triphenyl fulvenes is investigated using B3LYP‐SMD/6‐311+G(d,p) level density functional theory. The molecular electrostatic potential (MESP) minimum at the fulvene π‐system (V_min) and the change in MESP at any of the fulvene carbon atoms (ΔV_C) for both neutral and reduced forms are used as excellent measures of substituent effect from the para and meta positions of the 1,3 and 6‐phenyl moieties. Substitution at 6‐phenyl para position has led to significant change in E⁰ than any other positions. By applying the additivity rule of substituent effects, an equation in ΔV_C is derived to predict E⁰ for multiply substituted fulvenes. Further, E⁰ is predicted for a set of 2000 hexa‐substituted fulvene derivatives where the substituents and their positions in the system are chosen in a random way. The calculated E⁰ agreed very well with the experimental E⁰ reported by Godman et al. Predicting E⁰ solely by substituent effect offers a simple and powerful way to select suitable combinations of substituents on fulvene system for light harvesting applications. © 2018 Wiley Periodicals, Inc.     

Determination of the electron density in methyl (±)‐(1S,2S,3R)‐2‐methyl‐1,3‐diphenylcyclopropanecarboxylate using refinements with X‐ray scattering factors from wavefunction calculations of the whole molecule

The molecule of the title compound, C₁₈H₁₈O₂, is a substituted cyclopropane ring. The electron density in this molecule has been determined by refining single‐crystal X‐ray data using scattering factors derived from quantum mechanical calculations. Topological analysis of the electron densities in the three cyclopropane C—C bonds was carried out. The results show the effects of this substitution on these C—C bonds.     

Nonexistence of local maxima in molecular electrostatic potential maps

It has been rigorously established by means of classical electrostatic arguments, that molecular electrostatic potential maps are devoid of local maxima. This forms a generalization of the earlier works of Politzer and co-workers which were restricted to the case of atoms.     

Quantum chemistry aspects of the solvent effects on the ene reaction of 1‐Phenyl‐1,3,4‐triazolin‐2,5‐dione and 2‐methyl‐2‐butene

A density functional theory study was used to investigate the quantum aspects of the solvent effects on the kinetic and mechanism of the ene reaction of 1‐phenyl‐1,3,4‐triazolin‐2,5‐dione and 2‐methyl‐2‐butene. Using the B3LYP/6–311++ G(d,p) level of the theory, reaction rates have been calculated in the various solvents and good agreement with the experimental data has been obtained. Natural bond orbital analysis has been applied to calculate the stabilization energy of N18? H19 bond during the reaction. Topological analysis of quantum theory of atom in molecule (QTAIM) studies for the electron charge density in the bond critical point (BCP) of N18? H19 bond of the transition states (TSs) in different solvents shows a linear correlation with the interaction energy. It is also seen form the QTAIM analysis that increase in the electron density in the BCP of N18? H19, raises the corresponding vibrational frequency. Average calculated ratio of 0.37 for kinetic energy density to local potential energy density at the BCPs as functions of N18? H19 bond length in different media confirmed covalent nature of this bond. Using the concepts of the global electrophilicity index, chemical hardness and electronic chemical potentials, some correlations with the rate constants and interaction energy have been established. Mechanism and kinetic studies on 1‐phenyl‐1,3,4‐triazolin‐2,5‐dione and 2‐methyl‐2‐butene ene reaction suggests that the reaction rate will boost with interaction energy enhancement. Interaction energy of the TS depends on the solvent nature and is directly related to electron density of the bonds involved in the reaction proceeding, global electrophilicity index and electronic chemical potential. However, the chemical hardness relationship is reversed. Finally, an interesting and direct correlation between the imaginary vibrational frequency of the N18? H19 critical bond and its electron density at the TS has been obtained. © 2014 Wiley Periodicals, Inc.     

Interactions between chlorinated dioxins and a positively charged molecular probe: New molecular interaction potential

Interaction with the ligand binding domain of receptors for natural chemicals present one potential mechanism for the biological effects of environmental chemicals. Evidence suggests that the electrostatic interaction between the ligand and the receptor is an important component for binding to some of the relevant receptors. The presence of charged residues near the binding site suggests that the charge distribution of the free ligand may be different from the charge distribution of the ligand as it approaches the binding domain of the protein. In this study a new type of potential is computed for a series of dibenzo-p-dioxin (dioxin) ligands. This quantum mechanically computed potential results from interaction between the ligand and a trimethyl ammonium probe at a set of grid points. This interaction potential is compared with the molecular electrostatic potential computed from the wave function of the isolated ligands. Three types of local minima are found: (1) above the oxygen; (2) above the conjugated ring; and (3) above the chlorine(s). The molecular electrostatic potential emphasizes the minima associated with the chlorine atoms and, in that potential, the minima associated with the oxygen atoms disappear with chlorination. In the new potential, the minima over the oxygen atoms are maintained even in tetrachlorodioxin. As chlorination is increased the differences between the two potentials increases. The new potential shows the influence of the π-cation interaction, which is largest when there is little substitution on the ring. The presence of the probe induces a dipole component of 1 debye perpendicular to the plane of the ligand. Local minima in the interaction potential are then used as starting structures for the determination of the most stable ligand–probe complexes. The most stable structures are obtained from the minima associated with the oxygen atoms. These structures are stabilized by a hydrogen bond formation between the probe and the oxygen and the molecule is bent by 30° about the O(SINGLE BOND)O axis. For this series of molecules, the new potential retains some of the features that determine the hydrogen bond whereas the molecular electrostatic potential does not. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 673–684, 1998     

Density functional computational studies on 2‐[(2,4‐Dimethylphenyl)iminomethyl]‐3,5‐dimethoxyphenol

Density functional calculations of the structure, molecular electrostatic potential, and thermodynamic functions have been performed at B3LYP/6‐31G(d) level of theory for the title compound of 2‐[(2,4‐dimethylphenyl)iminomethyl]‐3,5‐dimethoxyphenol ( I ). To investigate the tautomeric stability, optimization calculations at B3LYP/6‐31G(d) level were performed for the enol and keto forms of I . Calculated results reveal that the enol form of I is more stable than its keto form. The predicted nonlinear optical properties of I are much greater than ones of urea. The changes of thermodynamic properties for the formation of the title compound with the temperature ranging from 200 to 500 K have been obtained using the statistical thermodynamic method. At 298.15 K, the change of Gibbs free energy for the formation reaction of I is 32.973 kJ/mol. The title compound can not be spontaneously produced from the isolated monomers at room temperature. The tautomeric equilibrium constant is computed as 0.868 at 298.15 K for enol‐imine?keto‐amine tautomerization of I . In addition, natural bond orbital analysis of I was performed using the B3LYP/6‐31G(d) method. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Experimental and theoretical investigation of the molecular and electronic structure of N′‐benzylidene‐N‐[4‐(3‐methyl‐3‐phenyl‐cyclobutyl)‐thiazol‐2‐yl]‐chloro‐acetic acid hydrazide

The title compound, N′‐benzylidene‐N‐[4‐(3‐methyl‐3‐phenyl‐cyclobutyl)‐thiazol‐2‐yl]‐chloro‐acetic acid hydrazide, has been synthesized and characterized by elemental analysis, IR, ¹H and ¹³C NMR, and X‐ray single crystal diffraction. The compound crystallizes in the orthorhombic space group P 21 21 21 with a = 5.8671 (3) Å, b = 17.7182 (9) Å, and c = 20.6373 (8) Å. Moreover, the molecular geometry from X‐ray experiment, the molecular geometry, vibrational frequencies, and gauge‐including atomic orbital ¹H and ¹³C chemical shift values of the title compound in the ground state have been calculated by using the Hartree–Fock and density functional methods (B3LYP) with 6‐31G(d) and 6‐31G(d,p) basis sets. The results of the optimized molecular structure are exhibited and compared with the experimental X‐ray diffraction. Besides, molecular electrostatic potential, Frontier molecular orbitals, and thermodynamic properties of the title compound were determined at B3LYP/6‐31G(d) levels of theory. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Similarities and Differences between Crystal and Enzyme Environmental Effects on the Electron Density of Drug Molecules

The crystal interaction density is generally assumed to be a suitable measure of the polarization of a low-molecular weight ligand inside an enzyme, but this approximation has seldomly been tested and has never been quantified before. In this study, we compare the crystal interaction density and the interaction electrostatic potential for a model compound of loxistatin acid (E64c) with those inside cathepsin B, in solution, and in vacuum. We apply QM/MM calculations and experimental quantum crystallography to show that the crystal interaction density is indeed very similar to the enzyme interaction density. Less than 0.1 e are shifted between these two environments in total. However, this difference has non-negligible consequences for derived properties.     

Experimental and quantum chemical calculational studies on N‐[4‐(3‐Methyl‐3‐phenyl‐cyclobutyl)‐thiazol‐2‐yl]‐N′‐pyridin‐3ylmethylene‐hydrazine

The title molecule, N‐[4‐(3‐Methyl‐3‐phenyl‐cyclobutyl)‐thiazol‐2‐yl]‐N′‐pyridin‐3ylmethylene‐ hydrazine (C₂₀ H₂₀ N₄ S₁), was characterized by ¹H‐NMR, ¹³C‐NMR, IR, UV‐visible, and X‐ray determination. In addition to the molecular geometry from X‐ray experiment, the molecular geometry, vibrational frequencies and gauge including atomic orbital ¹H‐ and ¹³C‐NMR chemical shift values of the title compound in the ground state have been calculated using the Hartree‐Fock and density functional method (B3LYP) with 6‐31G(d, p) basis set. The calculated results show that optimized geometries can well reproduce the crystal structural parameters. By using time‐dependent density functional theory method, electronic absorption spectrum of the title compound has been predicted. © 2011 Wiley Periodicals, Inc.     

Off‐atomic site charges in some anions,metal cations,and complexes: Significance for electrostatic properties and binding

Electronic structures and properties of several anions, metal cations, and their complexes with neutral molecules were investigated at the HF/6‐31G** and B3LYP/6‐31G** levels of theory. Charges shifted from atomic sites due to atomic orbital hybridization called hybridization displacement charges (HDC) were investigated in detail. It has been found that many components of HDC are associated with each atom of ion that are shifted from the atomic sites, those associated with metal cations being shifted by large distances as found previously in electrically neutral systems. It is shown that atomic orbitals are appreciably rehybridized in going from neutral molecules to anions and cations. Molecular dipole moments and surface molecular electrostatic potentials (MEP) are obtained satisfactorily using HDC for the various types of species mentioned above. In the OH^?? H₂O complex, reversal of direction of shift of an HDC component associated with the hydrogen atom of H₂O involved in hydrogen bonding, indicates that the hydrogen bond between OH^? and H₂O would have some covalent character. Other atomic site‐based point charge models cannot provide such information about the nature of bonding. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem 2007     

Extreme Carbon Dioxide Sorption Hysteresis in Open‐Channel Rigid Metal–Organic Frameworks

A systematic study is presented of three closely related microporous metal‐organic frameworks the pore dimensions of which vary according to the choice of 4,4′‐bipyridyl linker. The tunable linker allows exploration of the effect of increasing pore dimensions on the sorption behavior of the frameworks. The MOFs described capture CO₂ under supercritical conditions and continue to sequester the gas under ambient conditions. Gas sorption isotherms for CO₂ are compared with thermogravimetric data, and the CO₂ molecules in the channels of the frameworks could be modeled using single‐crystal X‐ray diffraction analysis. Crystallographic data were used to construct a theoretical model based on DFT methods to calculate framework electrostatic potential maps with a view to understanding the nature of the sorbate–sorbent interactions.     

A comparative study on the nature and strength of O–O, S–S, and Se–Se bond

A computational study on dichalcogenide molecules (R₂X₂; X = O, S, Se; R = H, CH₃, NH₂) has been carried out employing B3LYP and MP2 levels using 6-31+G*, 6-311+G*, 6-311++G**, and PVDZ basis sets. The relative energies have been evaluated at G2MP2 also. The rotational barriers and bond dissociation energies indicate that S–S bond is stronger than Se–Se and O–O bond. NBO analysis at MP2/6-31+G* suggest the presence of partial π character between X–X bond that decreases in the order S–S > Se–Se > O–O. Fuki functions for nucleophilic and electrophilic attack fail to distinguish the reactivity of S and Se. The proton affinities of the O₂H₂, S₂H₂, Se₂H₂ decrease in the order Se > S > O.     

An experimental and theoretical approach to the molecular structure of 3‐{[4‐(3‐Methyl‐3‐phenyl‐cyclobutyl)‐thiazol‐2‐yl]‐hydrazono}‐1,3‐dihydro‐ indol‐2‐one

The title molecule, 3‐{[4‐(3‐methyl‐3‐phenyl‐cyclobutyl)‐thiazol‐2‐yl]‐hydrazono}‐1,3‐dihydro‐indol‐2‐one (C₂₂H₂₀N₄O₁S₁), was prepared and characterized by ¹H NMR, ¹³C NMR, IR, UV–visible, and single‐crystal X‐ray diffraction. The compound crystallizes in the monoclinic space group P21 with a = 8.3401(5), b = 5.6976(3), c = 20.8155(14) Å, and β = 95.144(5)°. Molecular geometry from X‐ray experiment and vibrational frequencies of the title compound in the ground state has been calculated using the Hartree–Fock with 6‐31G(d, p) and density functional method (B3LYP) with 6‐31G(d, p) and 6‐311G(d, p) basis sets, and compared with the experimental data. The calculated results show that optimized geometries can well reproduce the crystal structural parameters, and the theoretical vibrational frequencies values show good agreement with experimental data. Density functional theory calculations of the title compound and thermodynamic properties were performed at B3LYP/6‐31G(d, p) level of theory. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

A computational study on the ring stretching modes of halogen‐substituted pyridine involved in H‐bonding

Density functional method B3LYP plus the AUG‐cc‐pVDZ and AUG‐cc‐pVTZ basis sets is used to investigate ring normal modes of halogen‐substituted pyridines involved in the N ··· H? X H‐bonds with HX (X = F, Cl). The results demonstrated that the formation of hydrogen bond leads to an increase in the frequencies of the ring breathing mode v₁, the N‐para‐C stretching mode v_6a and the meta‐CC stretching mode v_8a, whereas there is no change in the triangle mode v₁₂ for free pyridine and a smaller blue shift for substituted pyridines. There is a strong coupling between the C? Y stretching vibration and the triangle mode (ortho‐ and para‐substituted) or the breathing mode (meta‐substituted) in substituted pyridines, which leads to the frequency decrease in the triangle or breathing modes. The natural bond orbital analysis suggests that electrostatic interaction and charge transfer caused by the intermolecular and intramolecular hyperconjugations are the origin of the frequency blue shift in the ring stretching modes. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

The Role of Hydrogen Bonds in Interactions between [PdCl4]2− Dianions in Crystal

[PdCl₄]²⁻ dianions are oriented within a crystal in such a way that a Cl of one unit approaches the Pd of another from directly above. Quantum calculations find this interaction to be highly repulsive with a large positive interaction energy. The placement of neutral ligands in their vicinity reduces the repulsion, but the interaction remains highly endothermic. When the ligands acquire a unit positive charge, the electrostatic component and the full interaction energy become quite negative, signalling an exothermic association. Raising the charge on these counterions to +2 has little further stabilizing effect, and in fact reduces the electrostatic attraction. The ability of the counterions to promote the interaction is attributed in part to the H-bonds which they form with both dianions, acting as a sort of glue.