Hydrogen bond vs proton transfer in HZSM5 zeolite. A theoretical study

The interaction of a large set of bases covering a wide range of the basicity scale with HZSM5 medium-size zeolites has been investigated through the use of two model clusters, namely 5T and 7T:63T. The 5T cluster has been treated fully ab initio at the B3LYP level, whereas the 63T cluster has been treated with the ONIOM2 scheme using the B3LYP:MNDO combination for geometry optimizations and B3LYP:HF/3-21G for adsorption energies. The optimized geometries of the different hydrogen bond (HB) and ion pair (IP) complexes obtained with both models are rather similar. However, there are significant dissimilarities as far as the adsorption energies are concerned, in particular when dealing with IP clusters whose intrinsic stability is largely underestimated when the simpler 5T model is used. 5T clusters could be used to obtain reasonable estimates of adsorption energies provided these are scaled by a factor of 1.1 for HB complexes and 1.4 for IP complexes. The zeolite cavity favors the proton transfer process, similarly to that found by third polar partners in gas-phase HB trimers. The intrinsic basicity of the base and its adsorption energy within the zeolite are correlated. From this correlation, is possible to conclude that, in general, bases with proton affinities (PA) larger than 200 kcal mol(-1) should lead to the formation of IPs, whereas bases with PA smaller than this value should form HB complexes.     

The nicotinic pharmacophore: thermodynamics of the hydrogen-bonding complexation of nicotine, nornicotine, and models

The thermodynamics of the hydrogen-bonding complexation of the acetylcholine agonists nicotine and nornicotine and of model pyridines, pyrrolidines, and N-methylpyrrolidines has been measured in CCl(4) by FTIR spectrometry toward a reference hydrogen-bond donor, 4-fluorophenol. Various methods are devised for measuring separately the hydrogen-bond acceptor strength of each nitrogen of nicotine and nornicotine: variation of the stoichiometry of complexation; correlations with electrostatic potentials on nitrogens and with substituent constants in the series of 3-substituted pyridines, 2-substituted pyrrolidines, and 2-substituted N-methylpyrrolidines; and linear free energy relationships between 4-fluorophenol and hydrogen fluoride hydrogen-bonded complexes. It is consistently found that nicotine and nornicotine have two active hydrogen-bond acceptor sites, the pyridine and pyrrolidine nitrogens, and that ca. 90% (for nicotine) and 80% (for nornicotine) of the 1:1 hydrogen-bonded complexes are formed to the pyridine nitrogen, although the pyrrolidine nitrogen is the first protonation site of nicotine and nornicotine in water. The low hydrogen-bond basicity of the pyrrolidine nitrogen in nicotine is mainly explained by the inductive electron-withdrawing and steric effects of the 2-(3-pyridyl) substituent. The partition of the Gibbs energy of the isomerism of complexation (AH...Nsp(2) <==> AH...Nsp(3)) into enthalpic and entropic contributions shows that the selectivity in favor of the pyridine nitrogen is driven by entropy. It is important to recognize the bifunctionality of nicotine in hydrogen bonding for understanding its lipophilicity and molecular recognition in non protonic media. When monoprotonated on their sp(3) nitrogen, nicotine and nornicotine keep, through their sp(2) nitrogen, a significant hydrogen-bond basicity which is greater than that of the ester group of acetylcholine.     

Competition between nonclassical hydrogen-bonded acceptor sites in complexes of neutral AH2 Radicals (A = B, Al, and Ga): A theoretical investigation

An ab initio computational study of the properties of the neutral AH2 radicals (A = B, Al, Ga) as hydrogen-bond (HB) acceptors, with H-X (X = F, Cl, Br, CN, and CCH) as HB donors, is carried out at the UMP2/6-311++G(2d,2p) level. Two different minima have been found for each of the 15 possible dimers. One structure corresponds to a single-electron hydrogen-bonded complex (SEHB), with the A atom acting as an HB acceptor. The second corresponds to a dihydrogen bond complex between one of the hydrogen atoms of AH2 and the H-X molecule. Thus, all the atoms of the neutral AH2 molecule can act as HB acceptors and none as donors. The stability of the SEHB complexes decreases as BH2 > AlH2 > GaH2, while for the dihydrogen-bonded complexes the order is AlH2 > GaH2 > BH2. For the BH2 radical the SEHB complexes are stronger than the dihydrogen bonded ones, while the opposite is found for the AlH2 and GaH2 systems. Regarding the HB donors, the order found for the binding energy in the two types of complexes is H2A...HF > H2A...HCl > H2A...HBr > H2A...HCN > H2A...HCCH.     

Acidity constants from vertical energy gaps: density functional theory based molecular dynamics implementation

The question of how to compute acidity constants (pK(a)) treating solvent and solute at the same level of theory remains of some interest, for example in the case of high or low pH conditions. We have developed a density functional theory based molecular dynamics implementation of such a method. The method is based on a half reaction scheme computing free energies of dissociation from the vertical energy gaps for insertion or removal of protons. Finite system size effects are important, but largely cancel when half reactions are combined to full reactions. We verified the method by investigating a series of organic and inorganic acids and bases spanning a wide range of pK(a) values (20 units). We find that the response of the aqueous solvent to vertical protonation/deprotonation is almost always asymmetric and correlated with the strength of the hydrogen bonding of the deprotonated base. We interpret these observations in analogy with the picture of solvent response to electronic ionization.     

Energetic,Topological and Electric Field Analyses of Cation-Cation Nucleic Acid Interactions in Watson-Crick Disposition

A theoretical study of the effect of the diprotonation on the nucleic acid bases (A : U, A : T and G : C) in Watson-Crick conformation has been carried out by means of DFT computational methods in vacuum. In addition, the corresponding neutral and monoprotonated binary complexes have been considered. Most of the diprotonated species studied are stable, even though the binding energy is positive due to the overall repulsive electrostatic term. Local electrostatic attractive forces in the regions of hydrogen bonds (HBs) are responsible for equilibrium geometries, as shown by the electric field lines connecting the electrophilic and nucleophilic sites involved in the HB interactions. Secondary electrostatic effects also affect the assembling of the nucleic acid complexes in either neutral or cationic form. In particular, the electric field lines flowing from electrophilic sites in one base to nucleophilic sites in the other reinforce the linking between them. Hence, when the nucleophilic site concerns the free lone pair of the heteroatom involved in the HB interaction as acceptor, the HB distance shortens. However, if the free lone pair of the HB acceptor interacts with an electrophilic site in the same molecule, the HB distance elongates, weakening the HB interaction. The topological analysis of the electron density distribution in HB regions indicates that neutral, monoprotonated and diprotonated complexes show no differences in the nature of their HB's.     

PFOS在锐钛型TiO2表面吸附行为的理论研究

采用密度泛函理论(DFT)B3LYP方法对全氟辛烷磺酸(PFOS)在锐钛型TiO2表面的化学吸附和物理吸附行为进行了研究,其中化学吸附包含双齿双核(BB)和单齿单核(MM)在内的4种可能的吸附构型.吸附能(Eads)及反应吉布斯自由能(ΔGads)的计算结果表明,PFOS分子易于与TiO2表面发生氢键作用吸附;化学吸附表现为PFOS分子与TiO2表面的水分子(H2O)和羟基(—OH)反应,且与取代—OH相比,H2O取代相对更容易发生,其中,MM1构型(取代一个表面水分子)为化学吸附中的优势构型.PFOS在锐钛矿表面吸附的热力学稳定性和反应自发性顺序如下:H-Bonded(氢键吸附)>MM1(取代一个表面水分子)>BB1(取代两个表面水分子)>MM2(取代一个表面羟基)>BB2(取代一个表面水分子和一个表面羟基).成键结构分析表明,TiO2表面H2O/—OH官能团与PFOS上的磺酸基之间形成了中等强度的氢键;在化学吸附过程中,电荷从PFOS分子向TiO2表面发生转移,生成Ti—O—S化学键,电荷转移主要来自PFOS分子的O和F原子.     

Coupled cluster and density functional investigation of the hydrogen bond between halides,paraffines, olefins,and alkynes

By means of coupled cluster and density functional M06-2X calculations we have studied the directionality and energetics of the hydrogen bond (HB) between halides and the following H donors: methane, ethylene, and acetylene. The largest interaction energy was obtained for acetylene, the molecule in which the hydrogen atoms present the largest positive charge. For the complex [X···HCCH]^?, fluoride presented a different behavior since acetylene donates a hydrogen to fluoride and the [HF···CCH]^? complex was formed. We found that for ethylene there is a strong competition between the monodentate- and bidentate-like structures. Because of its small size, fluoride leads to a monodentate complex with C₂H₄. However, as we move down in the periodic table, the bidentate-like structure becomes more stable. In effect, for chloride and bromide they are nearly isoenergetic, while for iodide the bidentate-like complex is formed. These present results question previous investigations about the directionality of the HB in halide–ethylene complexes.     

Substituent effects on cobalt diglyoxime catalysts for hydrogen evolution

The design of efficient, robust, and inexpensive hydrogen evolution catalysts is important for the development of renewable energy sources such as solar cells. Cobalt diglyoxime complexes, Co(dRgBF(2))(2) with substituents R, are promising candidates for such electrocatalysts. The mechanism for hydrogen production requires a series of reduction and protonation steps for various monometallic and bimetallic pathways. In this work, the reduction potentials and pK(a) values associated with the individual steps were calculated for a series of substituents. The calculations revealed a linear relation between the reduction potentials and pK(a) values with respect to the Hammett constants, which quantify the electron donating or withdrawing character of the substituents. Additionally, the reduction potentials and pK(a) values are linearly correlated with each other. These linear correlations enable the prediction of reduction potentials and pK(a) values, and thus the free energy changes along the reaction pathways, to assist in the design of more effective cobaloxime catalysts.     

Gas-Phase and Solution-Phase Homolytic Bond Dissociation Energies of H-N(+) Bonds in the Conjugate Acids of Nitrogen Bases

The oxidation potentials of 19 nitrogen bases (abbreviated as B: six primary amines, five secondary amines, two tertiary amines, three anilines, pyridine, quinuclidine, and 1,4-diazabicyclo[2,2,2]octane), i.e., E(ox)(B) values in dimethyl sulfoxide (DMSO) and/or acetonitrile (AN), have been measured. Combination of these E(ox)(B) values with the acidity values of the corresponding acids (pK(HB)(+)) in DMSO and/or AN using the equation: BDE(HB)(+) = 1.37pK(HB)(+) + 23.1 E(ox)(B) + C (C equals 59.5 kcal/mol in AN and 73.3 kcal/mol in DMSO) gave estimates of solution phase homolytic bond dissociation energies of H-B(+) bonds. Gas-phase BDE values of H-B(+) bonds were estimated from updated proton affinities (PA) and adiabatic ionization potentials (aIP) using the equation, BDE(HB(+))(g) = PA + aIP - 314 kcal/mol. The BDE(HB)(+) values estimated in AN were found to be 5-11 kcal/mol higher than the corresponding gas phase BDE(HB(+))(g) values. These bond-strengthening effects in solution are interpreted as being due to the greater solvation energy of the HB(+) cation than that of the B(+*) radical cation.     

NMR parameters and geometries of OHN and ODN hydrogen bonds of pyridine-acid complexes

In this paper, equations are proposed which relate various NMR parameters of OHN hydrogen-bonded pyridine-acid complexes to their bond valences which are in turn correlated with their hydrogen-bond geometries. As the valence bond model is strictly valid only for weak hydrogen bonds appropriate empirical correction factors are proposed which take into account anharmonic zero-point energy vibrations. The correction factors are different for OHN and ODN hydrogen bonds and depend on whether a double or a single well potential is realized in the strong hydrogen-bond regime. One correction factor was determined from the known experimental structure of a very strong OHN hydrogen bond between pentachlorophenol and 4-methylpyridine, determined by the neutron diffraction method. The remaining correction factors which allow one also to describe H/D isotope effects on the NMR parameters and geometries of OHN hydrogen bond were determined by analysing the NMR parameters of the series of protonated and deuterated pyridine- and collidine-acid complexes. The method may be used in the future to establish hydrogen-bond geometries in biologically relevant functional OHN hydrogen bonds.     

Induction-driven stabilization of the anion-π interaction in electron-rich aromatics as the key to fluoride inclusion in imidazolium-cage receptors

Intermolecular interactions that involve aromatic rings are key processes in both chemical and biological recognition. It is common knowledge that the existence of anion-π interactions between anions and electron-deficient (π-acidic) aromatics indicates that electron-rich (π-basic) aromatics are expected to be repulsive to anions due to their electron-donating character. Here we report the first concrete theoretical and experimental evidence of the anion-π interaction between electron-rich alkylbenzene rings and a fluoride ion in CH(3)CN. The cyclophane cavity bridged with three naphthoimidazolium groups selectively complexes a fluoride ion by means of a combination of anion-π interactions and (C-H)(+)···F(-)-type ionic hydrogen bonds. (1)H NMR, (19)F NMR, and fluorescence spectra of 1 and 2 with fluoride ions are examined to show that only 2 can host a fluoride ion in the cavity between two alkylbenzene rings to form a sandwich complex. In addition, the cage compounds can serve as highly selective and ratiometric fluorescent sensors for a fluoride ion. With the addition of 1 equiv of F(-), a strongly increased fluorescence emission centered at 385 nm appears at the expense of the fluorescence emission of 2 centered at 474 nm. Finally, isothermal titration calorimetry (ITC) experiments were performed to obtain the binding constants of the compounds 1 and 2 with F(-) as well as Gibbs free energy. The 2-F(-) complex is more stable than the 1-F(-) complex by 1.87 kcal mol(-1), which is attributable to the stronger anion-π interaction between F(-) and triethylbenzene.     

Catalytic O[bond]O activation chemistry mediated by iron hangman porphyrins with a wide range of proton-donating abilities

[reaction: see text] A library of hanging porphyrin xanthene (HPX) compounds containing pendant groups with various proton-donating abilities (pK(a) ranging from approximately 2 to 25) has been synthesized. Their corresponding chloroiron(III) complexes promote the catalase-like disproportionation of hydrogen peroxide. The overall activity and turnover numbers (TONs) are maximal for iron HPX complexes bearing acidic hydrogen-bond pendants. These results establish that careful control of intramolecular proton inventory can dramatically influence the catalytic activation of O-O bonds.     

Double-cavity calix[4]pyrrole derivative with enhanced capacity for the fluoride anion

A double-cavity calix[4]pyrrole derivative, meso-tetramethyl-tetra[N-(2-phenoxyethyl)-N'-phenylurea]calix[4]pyrrole, 1, with enhanced hosting ability for the fluoride anion has been designed and characterized. Its interaction with anions (fluoride, chloride, bromide, iodide, dihydrogen phosphate, hydrogen sulfate, perchlorate, nitrate, and trifluoromethane sulfonate) was qualitatively and quantitatively assessed through 1H NMR, conductance, and calorimetric studies. The outcome of these investigations demonstrates that 1 interacts only with fluoride and dihydrogen phosphate anions in dipolar aprotic media. However, the composition of these complexes differs in that two units of fluoride are taken per unit of 1, while a 1:1 anion/ligand complex is formed with the dihydrogen phosphate anion. Results from the 1H NMR studies are striking in that these not only provide information about the active sites of the ligand-anion interaction but also allow the establishment of the sequence of events taking place during fluoride complexation. Thus, hydrogen-bond formation between the pyrrolic hydrogen and the fluoride anion is followed by the uptake of a second anion through the same type of interaction, but with the phenyl urea. It is also the latter group that is responsible for the interaction of 1 with the dihydrogen phosphate anion. Finally, this paper illustrates the importance of structural information for the interpretation of the thermodynamics associated with these systems.     

Assessing group-based cutoffs and the Ewald method for electrostatic interactions in clusters and in saturated,superheated, and supersaturated vapor phases of dipolar molecules

Monte Carlo simulations in the canonical, isobaric-isothermal, grand canonical, and Gibbs ensembles were used to assess whether the computationally expensive Ewald summation method for the computation of the first-order electrostatic energy can be replaced with a simpler truncation approach for accurate simulations of the saturated, superheated, and supersaturated vapor phases of dipolar and hydrogen-bonding molecules. Rotationally averaged hydrogen fluoride dimer and trimer energies, thermophysical properties and aggregation in the superheated vapor phase of hydrogen fluoride, nucleation free energy barriers for water, and the vapor–liquid coexistence properties of hydrogen fluoride and water were investigated over a wide range of state points. We find that for densities not too close to the critical density, results obtained from simulations using a spherical potential truncation based on neutral groups (molecules or fragments) for the Coulomb interactions are statistically identical to those obtained using the Ewald summation method. Use of the simpler spherical truncation results in a significant reduction of the computational effort for simulations employing molecular mechanics force fields and also allows for straightforward implementation of many-body expansion methods to compute the potential energy from electronic structure calculations of subsystems of the entire vapor-phase system.     

Titrametric characterization of pH-induced phase transitions in functionalized microgels

The chain and radial functional group distributions in carboxylic acid-functionalized poly(N-isopropylacrylamide)-based microgels have significant impacts on the types of swelling responses exhibited by the microgels upon the application of a temperature and/or pH stimulus. Potentiometric, conductometric, and calorimetric titration approaches are used in this work to characterize the chain distributions of -COOH groups in five microgels prepared using different -COOH-functionalized monomers. A direct correlation is observed between the kinetically predicted formation of functional monomer blocks within the microgel, the excess Gibbs free energy of ionization, and the apparent pK(a) versus degree of ionization profiles generated from potentiometric titration. Isothermal titration calorimetry (ITC) can be used to quantify the relative number of functional groups present in microgels prepared with the same functional monomer and/or identify differences between microgels with the same bulk -COOH content but different chain distributions. In particular, microgels prepared with diacid-functionalized monomers exhibit a characteristic two-step ITC profile. For microgels with the same bulk -COOH content, the heat of ionization measured via ITC increases systematically with the overall change in both the pK(a) and the excess Gibbs free energy for microgels prepared with monoacid-functionalized monomers. Diacid monomer-functionalized microgels have lower ionization enthalpies attributable to the break-up of hydrogen-bonded intramolecular ring complexes upon carboxylic acid ionization. The inferred chain functional group distributions can be used to understand differences in microgel swelling across the pH-induced phase transition.     

Reactions of cisplatin and glycine in solution with constant pH: a computational study

In this study interaction between glycine and cisplatin was explored. The structures were optimized in both gas phase and implicit water solution at the B3LYP/6-31G(p)/IEFPCM/UAKS computational level. Vibrational analysis confirms the corresponding character of all the complexes from examined reaction coordinates. The enthalpy and entropy contributions to Gibbs free energies were estimated from the vibrational frequencies and statistical canonical ensemble model. In the solution, several deprotonated forms were considered, and pK(a) values of the reactants and products were determined. The reaction Gibbs free energies and pK(a) values were used for determination of the transformed thermodynamic potential ΔG' with pH as a natural variable.     

Relaxation dynamics in (HF)x(H2O)1-x solutions

The high-frequency dynamics of (HF)(x)(H(2)O)(1-x) solutions has been investigated by inelastic x-ray scattering. The measurements have been performed as a function of the concentration in the range x = 0.20-0.73 at fixed temperature T = 283 K. The results have been compared with similar data in pure water (x = 0) and pure hydrogen fluoride (x = 1). A viscoelastic analysis of the data highlights the presence of a relaxation process characterized by a relaxation time and a strength directly related to the presence of a hydrogen-bond network in the system. The comparison with the data on water and hydrogen fluoride shows that the structural relaxation time continuously decreases at increasing concentration of hydrogen fluoride passing from the value for water to the one for hydrogen fluoride tau(alphaHF), which is three times smaller. This is the consequence of a gradual decreasing number of constraints of the hydrogen-bond networks in passing from one liquid to the other.     

Hydrogen-bond configuration in hydrogen-bond solutions. 1. Alcohol/alcohol solutions

A theory based on the hydrogen-bond configuration is proposed and applied to alcohol/alcohol binary solutions. The theory leads explicit expressions for the mixing Gibbs energy and reproduces the experiments on the crystal–liquid phase-diagrams of pure crystals and co-crystals and mixing heats with the parameters common to these experiments. The mixing entropy arises from the increase in the hydrogen-bonding availability of proton donors to approach hydrogen-bond-free proton acceptors. The mixing heat arises from a balance between the contribution from maintaining the original associations in pure liquids and the contribution from a construction of hydrogen bonds freely to hydrogen-bond-free acceptors. When hydrogen-bond associations between component-1 and component-2 are distinguished statistically from associations in each pure component, we call the solution as a cooperative solution that has at least one stoichiometric cooperative concentration point. Some shorter alcohol/alcohol solutions and some aromatic alcohol/aromatic alcohol solutions, however, have no cooperative point and we call those solutions as the ideal hydrogen-bond solutions of which properties are mainly governed by the ideal-gas-like mixing hydrogen-bond entropy. The hydrogen-bond energies of various combinations of the proton acceptor and the proton donor have been estimated consistently from the fittings of the theory, the shifts by hydrogen bonding of the OH stretching in the Raman or IR spectroscopy, and the sublimation energy of crystals. The present theory reveals the characteristics of hydrogen-bond solutions and gives some predictions.