用离散-连续模型计算NH2-,NH3和NH4+的溶剂化自由能

通过理论计算推测NH2-,NH3和NH4+在水溶液第一溶剂化层中与之直接作用的水分子分别为2,4和4个,并采用离散-连续模型计算了NH2-,NH3,NH3和NH4+在水溶液中的溶剂化自由能.结果表明,由于离散-连续模型在从头算水平考虑了溶质分子与第一溶剂化层溶剂分子之间的作用,能更准确地描述溶剂化作用.此外,采用更加符合溶液中真实情况的溶剂化构型,能得到更准确的溶剂化性质.     

醌型D-π-A类化合物AMTOE在不同溶剂中吸收和发射光谱的蓝移现象

采用连续介质模型(PCM)以及明确/连续的混合溶剂模型,运用密度泛函理论(DFT)和含时密度泛函理论(TD-DFT)方法,研究了溶剂极性对D-π-A类两性离子化合物1-(4-aza-4-methylphenyl)-2-trans-(4-oxyphenyl)ethane(AMTOE)分子几何、电子结构以及光谱性质的影响.计算结果表明,随着溶剂极性和超分子簇大小的增加,AMTOE基态(S0)分子几何从醌式结构向芳式结构转化.从基态(S0)到激发态(S1),化合物AMTOE醌式结构增强.从弱极性的氯仿溶液到强极性的水溶液,AMTOE分子的吸收光谱和发射光谱均发生蓝移,并且吸收光谱的蓝移程度大于发射光谱的蓝移程度,与实验现象定性一致.吸收光谱和发射光谱发生蓝移的原因是随着溶剂极性增强,HOMO轨道能级与LUMO轨道能级之间的能隙增大.对于具有明显电荷转移的AMTOE分子的溶剂化显色效应,长程矫正的含时密度泛函TD-CAM-B3LYP方法比传统的含时密度泛函TD-B3LYP方法更为合理.另外,明确/连续的混合溶剂模型能更好的描述该类体系在强极性溶剂中的溶剂化显色效应.     

过氧化氢与苯乙烯环氧化的反应机理及溶剂效应

以过氧化氢与苯乙烯环氧化为模型反应, 采用Materials Studio软件中Dmol³模块, 模拟计算了过氧化氢与苯乙烯的环氧化反应机理. 并用连续介质-类导体屏蔽模型(COSMO)研究了反应体系分别在三种质子性溶剂(水, 乙醇, 叔丁醇)中的溶剂化效应. 为研究溶剂分子直接参与反应的微观过程, 用离散介质模型模拟了单个水分子、乙醇分子和叔丁醇分子分别对反应的影响. 两种溶剂模型所得的结果一致, 叔丁醇作溶剂时反应活性最好, 乙醇次之; 质子性溶剂能够促进过氧化氢分子的异裂, 形成活性氧物种, 从而使反应能垒降低.     

有机酸取代基和溶剂效应的静电模型

本文根据气相实验结果和以往的静电理论, 建立了一个描述有机酸电离取代基和溶剂效应的静电模型。该模型不仅满足气相边界条件, 而且其参数B具有较明确的物理意义, 即与酸官能团的溶剂化半径以及取代基与溶剂分子之间特殊相互作用有关, 因而是描述溶质-溶剂相互作用的重要参数。     

密度泛函理论研究Gd（H2O）n^3＋（n=8，9）化合物的结构及相对稳定性

用密度泛函理论方法研究了气相和水溶液中Gd（H2O）n^3＋（n=8，9）化合物的结构和相对稳定性，其中水溶剂效应利用极化连续介质方法结合多种溶质空腔模型进行模拟．气相计算得到的化合物结构与实验观察结果一致．计算结果表明，在气相中9配位Gd（H2O）9^3＋比8配位Gd（H2O）8^3＋稳定，而在水溶液中稳定顺序刚好相反，这一结果不依赖于计算中采用的空腔模型种类，而且也与实验结果吻合．最后，通过采用各种空腔模型计算Gd^3＋的水合自由能，并与实验值比较，发现当化合物只包含第一层配位水分子时，UA0、UAHF及UAKS空腔模型最适合研究Gd^3＋在水溶液中的性质．     

2-巯基苯并咪唑及其类似物互变异构的理论研究

采用B3LYP/6-311G**方法, 计算了2-巯基苯并咪唑及其类似物(2-巯基苯并噁唑、2-巯基苯并噻唑、2-羟基苯并咪唑、2-羟基苯并噁唑、2-羟基苯并噻唑以及2-巯基咪唑、2-巯基噁唑、2-巯基噻唑、2-羟基咪唑、2-羟基噁唑、2-羟基噻唑)的(硫)醇式与(硫)酮式结构进行质子迁移的3种可能途径: (a)分子内质子迁移; (b)水助质子迁移; (c)甲醇助质子迁移．结果表明, 途经b和c所需要的活化能较小, 氢键在降低反应活化能方面起重要作用．采用PCM方法研究了反应体系的溶剂化效应．结果表明孤立分子、一水合物和一甲醇合物的最稳定异构体相同, 都为(硫)酮式, 与气相结论一致．溶剂化效应对异构化能垒的影响较小.     

ZSM-5分子筛硅铝比与其催化性能关系的量子化学研究

利用Gaussian03计算软件包,在DFT-B3LYP/3-21G水平上对ZSM-5分子筛128T簇模型进行了研究.通过分析T—O—T基团键角、桥羟基键长、质子电荷和分子轨道能隙等数据,讨论了ZSM-5分子筛硅铝比与其Br nsted酸性的关系,预测了其催化性能.计算结果表明,ZSM-5分子筛硅铝比越低,酸性越强,相应的催化性能越好,与实验结果相符.     

PP体系双荧光及分子内电荷转移的理论研究

用含时密度泛函理论（TDDFF）对PP（5-phenyl-5H-phenanthridin-6-one）分子进行了量子化学计算．实验观察到该分子具有双荧光和分子内电荷转移（ICT）特征．理论计算采用B3LYP泛函和6-311＋G（2d,p）基组,并用可极化连续介质模型（PCM）处理溶剂效应．基态结构优化表明PP分子为非平面状,苯环和菲啶酮平面的二面角近似等于90．0°．气相和溶剂中低激发态的计算结果表明这些激发态均为局域激发（LE）．通过构建气相和溶液中LE—ICT反应和荧光发射的势能曲线,得到了ICT反应的焓变△H、能垒Ea以及由于结构弛豫所产生的能量差δEFC．势能曲线的分析和发射能的计算显示PP是一个双荧光体系,包含一个LE发射谱和一个红移的ICT谱．考虑了溶剂效应的计算结果表明,双荧光的产生是由于在ICT反应过程中分子构型的变化,即苯环和菲啶酮单元由相互垂直变成了共平面．     

质子水合结构的振动态密度分析

使用基于多态经验价键模型的分子动力学模拟, 对水溶液中质子的水合结构及其在质子传递过程中的动力学过程进行了研究. 在价键模型的方法下, 质子的水合结构主要以H₉O₄⁺(Eigen)以及过渡态的H₅O₂⁺(Zundel)结构形态存在, 且在这两种结构中以Eigen的形态表现明显. 通过对质子传递过程中不同水合结构的态密度频谱分析, 发现一个在2000~3000 cm^-1范围内的明显连续的宽吸收谱带, 主要归因于Eigen结构的贡献, 这些特征峰的出现与水合氢离子第一溶剂化层内的强氢键作用密切相关. 对于Zundel的结构, 在1760 cm^-1处出现一个较为明显的肩峰, 归属为质子传递模式的特征振动. 通过对质子水合结构态密度频谱的分析, 可望增强对于稀酸溶液红外光谱中的连续宽吸收带以及质子传递的微观动力学过程的理解.     

Al(H2O)63+水交换反应溶剂效应的量子化学团簇模型-密度泛函理论方法研究

在气相模型、极化连续模型、超分子模型和超分子-极化连续模型的基础上,采用量子化学团簇模型密度泛函理论方法,在B3LYP/6-311+G(d,p)基组水平下系统地开展了以下研究:优化得到Al(H2O)63+水交换反应的反应物、过渡态和产物构型,采用MP2方法在相同基组水平下计算得到相应的单点能,考虑零点振动能、热力学校正项和熵等参数的影响,计算得到Al(H2O)63+水交换反应的Gibbs自由能变和反应速率常数kex.计算结果表明:GP-SM//MP2-PCM和GP-SM-PCM//MP2-PCM模型得到的kex相近,并且与文献值相符,说明GP-SM//MP2-PCM模型可以充分考虑真实溶剂效应和主体溶剂效应,适用于Al(H3O)63+体系水交换反应的模拟.     

Clusters of hydrated methane sulfonic acid CH3SO3H.(H2O)n (n = 1-5): a theoretical study

Ab initio and density functional methods have been used to examine the structures and energetics of the hydrated clusters of methane sulfonic acid (MSA), CH3SO3H.(H2O)n (n = 1-5). For small clusters with one or two water molecules, the most stable clusters have strong cyclic hydrogen bonds between the proton of OH group in MSA and the water molecules. With three or more water molecules, the proton transfer from MSA to water becomes possible, forming ion-pair structures between CH3SO3- and H3O+ moieties. For MSA.(H2O)3, the energy difference between the most stable ion pair and neutral structures are less than 1 kJ/mol, thus coexistence of neutral and ion-pair isomers are expected. For larger clusters with four and five water molecules, the ion-pair isomers are more stable (>10 kJ/mol) than the neutral ones; thus, proton transfer takes place. The ion-pair clusters can have direct hydrogen bond between CH3SO3- and H3O+ or indirect one through water molecule. For MSA.(H2O)5, the energy difference between ion pairs with direct and indirect hydrogen bonds are less than 1 kJ/mol; namely, the charge separation and acid ionization is energetically possible. The calculated IR spectra of stable isomers of MSA.(H2O)n clusters clearly demonstrate the significant red shift of OH stretching of MSA and hydrogen-bonded OH stretching of water molecules as the size of cluster increases.     

The effect of hydrogen bonding on the excited-state proton transfer in 2-(2'-hydroxyphenyl)benzothiazole: a TDDFT molecular dynamics study

The dynamics of the excited-state proton transfer (ESPT) in a cluster of 2-(2'-hydroxyphenyl)benzothiazole (HBT) and hydrogen-bonded water molecules was investigated by means of quantum chemical simulations. Two different enol ground-state structures of HBT interacting with the water cluster were chosen as initial structures for the excited-state dynamics: (i) an intramolecular hydrogen-bonded structure of HBT and (ii) a cluster where the intramolecular hydrogen bond in HBT is broken by intermolecular interactions with water molecules. On-the-fly dynamics simulations using time-dependent density functional theory show that after photoexcitation to the S(1) state the ESPT pathway leading to the keto form strongly depends on the initial ground state structure of the HBT-water cluster. In the intramolecular hydrogen-bonded structures direct excited-state proton transfer is observed within 18 fs, which is a factor two faster than proton transfer in HBT computed for the gas phase. Intermolecular bonded HBT complexes show a complex pattern of excited-state proton transfer involving several distinct mechanisms. In the main process the tautomerization proceeds via a triple proton transfer through the water network with an average proton transfer time of approximately 120 fs. Due to the lack of the stabilizing hydrogen bond, intermolecular hydrogen-bonded structures have a significant degree of interring twisting already in the ground state. During the excited state dynamics, the twist tends to quickly increase indicating that internal conversion to the electronic ground state should take place at the sub-picosecond scale.     

Theoretical study of water-pyridine complexes using intermolecular potentials

Intermolecular potential functions have been used to determine the equilibrium structures of the water-pyridine complexes. The dimer and symmetrical 2:1 water pyridine systems have been studied. Three water models, ST2, TIPS2, and EMPWI have been combined with two different Lennard Jones nonbonded parameters and various charge distributions for the pyridine molecule to describe the systems. For the dimer, results show two distinguishable classes of preferential hydration sites, which are specific sites corresponding to hydrogen-bonded dimer and nonspecific sites located near the hydrophobic regions. Calculations performed on hydrogen-bonded symmetrical complexes show that the planar complex is generally less stable than the complex with water molecules perpendicular to the pyridine plane. For these complexes, the major factor that influences the hydrogen-bonded configurations is the choice of the water model. The importance of atomic charge distributions for the solute over the choice of potential parameters is pointed out. Finally, the effective lone pair representation on the aromatic nitrogen atom is shown to improve the hydrogen bond geometry and the stability of the complexes.     

Hydrogen-bonding and Proton Transfer Interactions between Harmane and Trifluoroethanol in the Ground and Excited Singlet States

A study of the hydrogen-bonding and proton transfer reactions of the ground and excited states of harmane (1-methyl-9H-pyrido/3,4-b/indole) and its N ₉-methyl derivative with 2,2,2-trifluoroethanol in cyclohexane is reported. Spectral measurements (UV–visible, Fourier trans-form IR, steady-state and time-resolved fluorescence) show the formation of fluorescent ground-state hydrogen-bonded complexes. The results have been interpreted assuming a tautomeric equilibrium between a 1:1 hydrogen-bonded complex and its 1:2 proton transfer tautomer (hydrogen-bonding ion pair). Upon excitation to its singlet excited state, the proton transfer tautomer of harmane reacts with an additional 2,2,2-trifluoroethanol molecule to give a zwitterionic exciplex, which fluoresces at longer wavelength.     

Excited-state proton transfer via hydrogen-bonded acetic acid (AcOH) wire for 6-hydroxyquinoline

Spectroscopic studies on excited-state proton transfer (ESPT) of hydroxyquinoline (6HQ) have been performed in a previous paper. And a hydrogen-bonded network formed between 6HQ and acetic acid (AcOH) in nonpolar solvents has been characterized. In this work, a time-dependent density functional theory (TDDFT) method at the def-TZVP/B3LYP level was employed to investigate the excited-state proton transfer via hydrogen-bonded AcOH wire for 6HQ. A hydrogen-bonded wire containing three AcOH molecules at least for connecting the phenolic and quinolinic -N- group in 6HQ has been confirmed. The excited-state proton transfer via a hydrogen-bonded wire could result in a keto tautomer of 6HQ and lead to a large Stokes shift in the emission spectra. According to the results of calculated potential energy (PE) curves along different coordinates, a stepwise excited-state proton transfer has been proposed with two steps: first, an anionic hydrogen-bonded wire is generated by the protonation of -N- group in 6HQ upon excitation to the S(1) state, which increases the proton-capture ability of the AcOH wire; then, the proton of the phenolic group transfers via the anionic hydrogen-bonded wire, by an overall "concerted" process. Additionally, the formation of the anionic hydrogen-bonded wire as a preliminary step has been confirmed by the hydrogen-bonded parameters analysis of the ESPT process of 6HQ in several protic solvents. Therefore, the formation of anionic hydrogen-bonded wire due to the protonation of the -N- group is essential to strengthen the hydrogen bonding acceptance ability and capture the phenolic proton in the 6HQ chromophore.     

Structure, energetics, and infrared spectra of weakly bound HC(2n+1)N···HCl complexes. A theoretical study

Three model systems, HCN···HCl, HC(3)N···HCl, and HC(5)N···HCl, have been investigated computationally with the use of the second-order M?ller-Plesset (MP2) and the coupled cluster (with single and double excitations and noniterative inclusion of triples) methods. The global minima are linear hydrogen-bonded structures with HCl as a proton donor. Bent structures are proton-side complexes with HCl as an electron donor, while the bifurcated hydrogen bonds are predicted for t-shape complexes. One of the most important findings in this paper is that, according to symmetry-adapted perturbation analysis, the induction-to-dispersion ratios are the biggest for linear complexes, and it is the most noticeable difference between linear, bent, and t-shape structures.     

Characterization of hydrated Na+(phenol) and K+(phenol) complexes using infrared spectroscopy

Hydrated alkali metal ion-phenol complexes were studied to model these species in aqueous solution for M=Na and K. IR predissociation spectroscopy in the O-H stretch region was used to analyze the structures of M+(Phenol)(H2O)n cluster ions, for n = 1-4. The onset of hydrogen bonding was observed to occur at n=4. Ab initio calculations were used to qualitatively explore the types of hydrogen-bonded structures of the M+(Phenol)(H2O)4 isomers. By combining the ab initio calculations and IR spectra, several different structures were identified for each metal ion. In contrast to benzene, detailed in a previous study of Na+(Benzene)n(H2O)m [J. Chem. Phys. 110, 8429 (1999)], phenol is able to bind directly to Na+ even in the presence of four waters. This is likely the result of the sigma-type interaction between the phenol oxygen and the ion. With K+, the dominant isomers are those in which the phenol O-H group is involved in a hydrogen bond with the water molecules, while with Na+, the dominant isomers are those in which the phenol O-H group is free and the water molecules are hydrogen-bonded to each other. Spectra and ab initio calculations for the M+(Phenol)Ar cluster ions for M=Na and K are reported to characterize the free phenol O-H stretch in the M+(Phenol) complex. While pi-type configurations were observed for binary M+(Phenol) complexes, sigma-type configurations appear to dominate the hydrated cluster ions.     

Proton transfer and the mobilities of the H+ and OH- ions from studies of a dissociating model for water

Hydrogen (H(+)) and hydroxide (OH(-)) ions in aqueous solution have anomalously large diffusion coefficients, and the mobility of the H(+) ion is nearly twice that of the OH(-) ion. We describe molecular dynamics simulations of a dissociating model for liquid water based on scaling the interatomic potential for water developed by Ojama?e-Shavitt-Singer from ab initio studies at the MP2 level. We use the scaled model to study proton transfer that occurs in the transport of hydrogen and hydroxide ions in acidic and basic solutions containing 215 water molecules. The model supports the Eigen-Zundel-Eigen mechanism of proton transfer in acidic solutions and the transient hyper-coordination of the hydroxide ion in weakly basic solutions at room temperature. The free energy barriers for proton transport are low indicating significant proton delocalization accompanying proton transfer in acidic and basic solutions. The reorientation dynamics of the hydroxide ion suggests changes in the proportions of hyper-coordinated species with temperature. The mobilities of the hydrogen and hydroxide ions and their temperature dependence between 0 and 50 °C are in excellent agreement with experiment and the reasons for the large difference in the mobilities of the two ions are discussed. The model and methods described provide a novel approach to studies of liquid water, proton transfer, and acid-base reactions in aqueous solutions, channels, and interfaces.     

Formation of hydrogen-bonded complexes of 3,3′,5,5′-tetrabromo-2,2′-biphenol with MTBD and triethylamine

The complexes of 3,3′,5,5′-tetrabromo-2,2′-biphenol (TBBPh) with 7-methyl-1,5,7-triazabicyclo[4.4.0]dec-5-ene (MTBD) and triethylamine (TEA) were studied by FTIR spectroscopy. In chloroform and in acetonitrile a proton transfer from TBBPh to N-bases (MTBD, TEA) occurs. In chloroform solution the protonated N-base molecules are hydrogen-bonded to the deprotonated TBBPh molecules, whereas in acetonitrile the complexes dissociate. The intra- as well as intermolecular hydrogen bonds within the chains show large proton polarizability.     

Cooperativity and proton transfer in hydrogen-bonded triads.

Ab initio MP2/6-311+G(3df,2pd) and MP2/aug-cc-pVTZ calculations have been carried out to investigate the structures and properties of AHXHYH(3) (A=F, Cl; X=F, Cl; Y=N, P) hydrogen-bonded complexes. Significant cooperative effects are observed in the XHYH3 dyads in the triads due to the presence of the polar near-neighbor AH. These effects are greater when the polar partner is HF, which is a better proton donor than HCl. Structural changes, red shifts of proton-donor stretching frequencies, nonadditive interaction energies, and electron density redistributions unambiguously demonstrate that the X--HY hydrogen bond (HB) is stronger in the triads than in the corresponding dyads, while the X--H bond of the proton donor becomes weaker. Even more pronounced cooperative effects are observed in the AHXH dyads due to the presence of the YH3 partner. These effects are weaker in complexes having PH3 rather than NH3 as the proton acceptor, since NH3 is a stronger base. Cooperativity also enhances the proton-donating ability of the YH3 moiety, with the result that all complexes except FHFHPH3 are cyclic. Cooperativity, together with the ease of breaking the Cl--H bond in ClHClHNH3 and FHClHNH3, leads to proton transfer (PT), so that these two complexes are better described as approaching hydrogen-bonded ClHCl- x +HNH3 and FHCl- x +HNH3 ion pairs.