电子转移溶剂重组能计算的自洽反应场新方法

基于非平衡溶剂化理论, 推导了用于非平衡溶剂化能数值计算的类导体屏蔽模型(COSMO)的相关公式. 在此基础上, 修改了HONDO99中COSMO模块, 并用以估算了[(CH₂)₂C]+—(CH₂)_n—C(CH₂)₂(n＝1～13)体系中的电子转移溶剂重组能. 结果表明, 溶剂重组能值与电子转移距离的倒数有很好的线性关系. 根据溶剂重组能数值解结果, 用新的双球模型给出了合理的给受体球半径.     

电子转移反应O2+O2·-→O2·-+O2的从头算研究

利用abinitio方法,在UHF,UMP2及不同基组3-21G,6-31G^*,6-311+G^*和UMP2(full)/6-311+G^*水平上,研究了O~2/O~2^.^-自交换电子转移反应。优化了电子转移前后反应物和产物的结构,研究了体系能量的变化,计算了自交换电子转移反应的内重组能。对UHF方法和UMP2方法的计算结果进行了比较,并与实验结果进行了对照。结果表明UHF方法由于没有考虑组态相互作用,计算结果存在较大偏差,UMP2(full)/6-311+G^*水平上的计算结果与实验值吻合较好。在UMP2(full)/6-311+G^*水平上计算了气相自交换电子转移反应速率常数。在优化了电子转移复合物结构的基础上考虑了溶剂效应的影响,计算了水溶液中的溶剂重组能。研究结果表明O~2/O~2^.^-体系电子转移反应的活化能主要来源于溶剂重组能的贡献。最后计算了该反应在水溶液中的反应速率常数。理论计算结果与实验值吻合得很好。     

电子转移反应O~2+O~2^.^-→O~2^.^-+O~2的从头算研究

利用abinitio方法,在UHF,UMP2及不同基组3-21G,6-31G^*,6-311+G^*和UMP2(full)/6-311+G^*水平上,研究了O~2/O~2^.^-自交换电子转移反应。优化了电子转移前后反应物和产物的结构,研究了体系能量的变化,计算了自交换电子转移反应的内重组能。对UHF方法和UMP2方法的计算结果进行了比较,并与实验结果进行了对照。结果表明UHF方法由于没有考虑组态相互作用,计算结果存在较大偏差,UMP2(full)/6-311+G^*水平上的计算结果与实验值吻合较好。在UMP2(full)/6-311+G^*水平上计算了气相自交换电子转移反应速率常数。在优化了电子转移复合物结构的基础上考虑了溶剂效应的影响,计算了水溶液中的溶剂重组能。研究结果表明O~2/O~2^.^-体系电子转移反应的活化能主要来源于溶剂重组能的贡献。最后计算了该反应在水溶液中的反应速率常数。理论计算结果与实验值吻合得很好。     

非平衡溶剂效应的约束平衡理论及在电子转移中的应用

基于经典热力学约束平衡方法,采用新非平衡溶剂化理论研究了Np O+2-Np O2+2体系电子转移反应的溶剂重组能,采用限制密度泛函理论实现电荷定域,用积分方程可极化连续介质模型获得水溶液中极化电荷.计算结果表明,新双球模型和数值解都给出了一致的溶剂重组能理论计算值.     

有机共轭体系电子转移反应的溶剂重组能

应用Marcus双球模型计算溶剂重组能λs时,在AM1法优化给受体几何构型基础上,提出了共轭体系电子云分布的扁球模型,并用统计的方法求出了rD/A.同时依照Miller等的处理办法,结合其他理论及实验证据将电子转移交叉反应中联苯分子的扭转能计入溶剂重组能λs中,从而用实验速率常数拟合出含扭转能的λs值.此实验拟合值与扁球法得到的λs计算值吻合得很好.通过比较理论值与实验值,发现了给受体间距的大小、受体分子的变化、溶剂的不同对λs计算值相对λs实验值的偏差的影响,直接证实了电子给受体的耦合作用,溶剂分子参与的超交换电子转移及溶质溶剂分子表面相互作用等量子因素造成的实际反应体系对溶剂经典连续介质模型的偏离.     

非平衡溶剂化的连续介质模型和超快过程溶剂效应

在连续介质理论基础上, 根据Jackson的能量积分公式导出非平衡态静电自由能和溶剂化能的正确表达式. 引入“弹簧能”概念, 对平衡态和非平衡态的静电能构成给出了合理解释, 即此能量由溶质自由电荷和溶剂极化电荷的自能、 两者之间的相互作用能和极化电荷的“弹簧能”构成. 对目前几种代表性的非平衡溶剂化理论进行了论证和比较, 指出其中存在的基本理论问题. 根据新的非平衡溶剂化能建立了电子转移反应溶剂重组能的双球模型、 光谱移动的单球孔穴点偶极模型, 多级展开方法和非平衡溶剂效应的数值解方法.在Poisson方程求解中引入类导体屏蔽模型, 建立了任意孔穴极化电荷数值解方法并应用到Closs-Miller电子转移体系, 得到与实验值吻合的溶剂重组能, 解决了传统非平衡溶剂化理论高估溶剂重组能的问题.     

含酪氨酸和色氨酸的肽内电子转移的溶剂重组能计算

基于连续介质模型，本文考察了多肽体系Trp-(Pro)n-Tyr (n=1,2) 从酪氨酸到色氨酸的分子内电子转移，并根据电荷定域的反应物和产物构型和线性反应坐标近似构造了电子转移的双势阱，通过势能曲线的交叉点确定了电子转移过渡态。本文重点讨论了电子转移溶剂重组能。根据作者的非平衡溶剂化理论和可极化连续介质模型编写了溶剂重组能计算程序并用于本文体系的计算。计算得到Trp-Pro-Tyr 和Trp-(Pro)2-Tyr.体系的溶剂重组能分别为20.89 kcal/mol和25.30 kcal/mol.     

联苯负离子自由基与萘分子间长程电子转移的从头算研究

分离处理π σ π体系中的给体 ,受体和σ桥体 ,在HF/4 31G和HF/DZP水平上优化了联苯 ,联苯负离子自由基 ,萘和萘负离子自由基的几何构型 ,计算了分子间电子转移的内重组能 .取线性反应坐标R =0 .5 ,在STO 3G水平上用变分原理和分子轨道跃迁能方法 ,计算了π σ π体系自交换反应的电子转移矩阵元 .对交叉反应体系 ,沿线性反应坐标搜寻最小轨道能级分裂Δmin,确定了电子转移矩阵元和过渡态构型 .用Marcus双球模型计算液相电子转移的溶剂重组能 ,结合半经典模型计算了几种以联苯负离子自由基为给体 ,联苯和萘为受体的π σ π体系分子内电子转移速率常数 .     

溶剂-导体界面电子转移溶剂重组能的球-界面模型

在连续介质理论基础上, 根据热力学基本原理, 用一个外加电场Eex将非平衡态2[Enon2, Dnon2]变成约束平衡态[E*2, D*2], 推导出了正确普适的溶剂重组能公式. 基于球-界面近似, 推导出了正确的溶剂-导体界面电子转移溶剂重组能公式. 和Marcus的公式相比, 本文的结果多了(εs-εop)/(εop(εs-1))因子. 对极性溶剂, 预测的溶剂重组能约为Marcus模型所得结果的一半. 以C343(Coumarin 343)-TiO2体系为算例, 计算了溶剂重组能并与实验值进行了比较.     

苯硝化反应中电子转移反应重组能的计算

提出了重组能的量子化学算法,在用CISD/6-31G基组水平上,得到苯硝化反应中反应物及过渡态的结构.并计算了各自交换电子转移反应以及交叉电子转移反应的重组能,同实验重组能进行了比较.计算用了Gaussian 94程序.从重组能的角度分析了苯硝化反应.结果表明,对于NO2++NO2→NO2+NO2+的自交换电子转移反应,重组能较大,结论为: 在芳烃硝化反应中,存在以NO2+为氧化剂的电子转移步骤的可能性很小,而从动力学的角度上,用NO+作反应的氧化剂更有可能.     

Pendant unit effect on electron tunneling in U-shaped molecules

The electron transfer reactions of three U-shaped donor-bridge-acceptor molecules with different pendant groups have been studied in different solvents as a function of temperature. Analysis of the electron transfer kinetics in nonpolar and weakly polar solvents provides experimental reaction Gibbs energies that are used to parameterize a molecular solvation model. This solvation model is then used to predict energetic parameters in the electron transfer rate constant expression and allow the electronic coupling between the electron donor and electron acceptor groups to be determined from the rate data. The U-shaped molecules differ by alkylation of the aromatic pendant group, which lies in the ‘line-of-sight’ between the donor and acceptor groups. The findings show that the electronic coupling through the pendant group is similar for these molecules.     

High pressure ESR studies of electron self-exchange reactions of organic radicals in solution

Simple electron self-exchange reactions are often used to study the role of the reaction medium on a chemical process, commonly implying the use of various solvents with different physical properties. In principle, similar studies may be conducted using a single solvent, changing its physical properties by application of elevated pressures, but so far only little information is available on pressure dependent exchange reactions. In this work, we have used a recently constructed high pressure apparatus for use with electron spin resonance (ESR) spectroscopy to investigate simple electron self-exchange reactions involving 2,3-dichloro-5,6-dicyanobenzoquinone (DDQ) and tetracyanoethylene (TCNE) and their respective radical anions as well as TMPPD and its radical cation in three different solvents. The self-exchange was observed by ESR line broadening experiments, yielding rate constants and volumes of activation. The experimental results were compared to theoretical calculations based on Marcus theory and taking into account solvent dynamic effects. The use of elevated pressures has enabled the study of solvent effects without commonly encountered problems like solubility issues or chemical reactions between solvent and solute which sometimes limit the range of useable solvents.     

A phenomenological model of dynamical arrest of electron transfer in solvents in the glass-transition region

A phenomenological model of electron transfer reactions in solvents undergoing glass transition is discussed. The reaction constant cuts off slow polarization modes from the spectrum of nuclear thermal motions active on the observation time scale. The arrest of nuclear solvation in turn affects the reaction activation barrier making it dependent on the rate. The resultant rate constant is sought from a self-consistent equation. The model describes well the sharp change in the solvent Stokes shift of optical lines in the glass-transition region. It is also applied to describe the temperature dependence of primary charge separation and reduction of primary pair in photosynthetic reaction centers. The model shows that a weak dependence of the primary charge separation rate on temperature can be explained by dynamical arrest of nuclear solvation on the picosecond time scale of electron transfer. For reduction of primary pair by cytochrome, the model yields a sharp turnover of the reaction kinetics at the transition temperature when nuclear solvation freezes in.     

Dynamical arrest of electron transfer in liquid crystalline solvents

We argue that electron transfer reactions in slowly relaxing solvents proceed in the nonergodic regime, making the reaction activation barrier strongly dependent on the solvent dynamics. For typical dielectric relaxation times of polar nematics, electron transfer reactions in the subnanosecond time scale fall into nonergodic regime in which nuclear solvation energies entering the activation barrier are significantly lower than their thermodynamic values. The transition from isotropic to nematic phase results in weak discontinuities of the solvation energies at the transition point and the appearance of solvation anisotropy weakening with increasing solute size. The theory is applied to analyze experimental kinetic data for the electron transfer kinetics in the isotropic phase of 5CB liquid crystalline solvent. We predict that the energy gap law of electron transfer reactions in slowly relaxing solvents is characterized by regions of fast change of the rate at points where the reaction switches between the ergodic and nonergodic regimes. The dependence of the rate on the donor-acceptor separation may also be affected in a way of producing low values for the exponential falloff parameter.     

Bidirectional and unidirectional PCET in a molecular model of a cobalt-based oxygen-evolving catalyst

The oxidation of water to molecular oxygen is a kinetically demanding reaction that requires efficient coupling of proton and electron transfer. The key proton-coupled electron transfer (PCET) event in water oxidation mediated by a cobalt-phosphate-based heterogeneous catalyst is the one-electron, one-proton conversion of Co(III)-OH to Co(IV)-O. We now isolate the kinetics of this PCET step in a molecular Co(4)O(4) cubane model compound. Detailed electrochemical, stopped-flow, and NMR studies of the Co(III)-OH to Co(IV)-O reaction reveal distinct mechanisms for the unidirectional PCET self-exchange reaction and the corresponding bidirectional PCET. A stepwise mechanism, with rate-limiting electron transfer is observed for the bidirectional PCET at an electrode surface and in solution, whereas a concerted proton-electron transfer displaying a moderate KIE (4.3 ± 0.2), is observed for the unidirectional self-exchange reaction.     

吩噻嗪与四甲基哌啶氧铵正离子的电子转移反应动力学研究

以Marcus半经典电子转移理论为基本框架, 改进了重组能的计算方法, 建立了一套研究自交换和交叉电子转移反应的理论方案。用密度泛函理论和半经验分子轨道理论具体研究了四甲基哌啶氧铵正离子与吩噻嗪在乙腈溶液中的交叉电子转移反应以及相应的2个自交换反应的动力学性质, 计算了反应的活化能、重组能、耦合矩阵元等有关参数,获得了和实验结果相一致的电子转移速率常数。     

构型对芳基桥体分子内电子转移矩阵元影响的理论研究

在HF/6-31G水平上,研究了有机化合物二甲氧基-4-甲苯-四甲苯基桥体-二甲氧基-4-甲苯正离子间的电子转移.用线性反应坐标确定电子转移的过渡态,用两态变分法计算了电子转移矩阵元V_AB,在考虑非平衡态溶剂化效应下,计算了电子转移速率常数.改变桥体与氧化还原中心的二面角,计算了相应的电子转移矩阵元.通过电子结构分析,将电子转移矩阵元分为通过空间的直接耦合和通过键的耦合,提出了将后者进一步分为通过σ键和π键的耦合.得出通过π键的耦合与二面角余弦的平方成正比的结论.     

Computational Analysis of Solute–Solvent Coupling Magnitude in the Z/E Isomerization Reaction of Nitroazobenzene and Benzylideneanilines

The dynamic solvent effect often arises in solution reactions, where coupling between chemical reaction and solvent fluctuation plays a decisive role in the reaction kinetics. In this study, the Z/E isomerization reaction of nitoroazobenzene and benzylideneanilines in the ground state was computationally studied by molecular dynamics simulations. The non-equilibrium solvation effect was analyzed using two approaches: (1) metadynamics Gibbs energy surface exploration and (2) solvation Gibbs energy evaluation using a frozen solvation droplet model. The solute–solvent coupling parameter (C_coupled) was estimated by the ratio of the solvent fluctuation Gibbs energy over the corresponding isomerization activation Gibbs energy. The results were discussed in comparison with the ones estimated by means of the analytical models based on a reaction–diffusion equation with a sink term. The second approach using a frozen solvation droplet reached qualitative agreement with the analytical models, while the first metadynamics approach failed. This is because the second approach explicitly considers the non-equilibrium solvation in the droplet, which consists of a solute at the reactant geometry immersed in the pre-organized solvents fitted with the solute at the transition state geometry.     

Anisotropic solvent model of the lipid bilayer. 1. Parameterization of long-range electrostatics and first solvation shell effects

A new implicit solvation model was developed for calculating free energies of transfer of molecules from water to any solvent with defined bulk properties. The transfer energy was calculated as a sum of the first solvation shell energy and the long-range electrostatic contribution. The first term was proportional to solvent accessible surface area and solvation parameters (σ(i)) for different atom types. The electrostatic term was computed as a product of group dipole moments and dipolar solvation parameter (η) for neutral molecules or using a modified Born equation for ions. The regression coefficients in linear dependencies of solvation parameters σ(i) and η on dielectric constant, solvatochromic polarizability parameter π*, and hydrogen-bonding donor and acceptor capacities of solvents were optimized using 1269 experimental transfer energies from 19 organic solvents to water. The root-mean-square errors for neutral compounds and ions were 0.82 and 1.61 kcal/mol, respectively. Quantification of energy components demonstrates the dominant roles of hydrophobic effect for nonpolar atoms and of hydrogen-bonding for polar atoms. The estimated first solvation shell energy outweighs the long-range electrostatics for most compounds including ions. The simplicity and computational efficiency of the model allows its application for modeling of macromolecules in anisotropic environments, such as biological membranes.     

卟啉-酞菁分子内能量传递和电子转移的溶剂效应

用稳态荧光光谱研究了以氧原子和哌嗪作为连接基的卟啉酞菁二元分子在不同溶剂中的分子内能量传递和电子转移过程结果表明;分子内的能量传递和电子转移是两个相互竞争的过程，在非极性溶剂中，激发单重态的能量传递是主要过程，而在极性溶剂中则以电子转移为主运用Rehm-Weller公式计算了两种二元化合物在不同溶剂中的电子转移反应的自由能变化△G0ET，表明溶剂的极性对电子转移反应的自由能变化△G0ET影响很大极性越大;体系中的电子转移反应的△G0ET、越负，电子转移反应越易进行由于电子转移过程较能量传递过程进行得快，所以表现为体系中能量传递效率降低而电子转移效率增大。两种二元化合物的能量传递效率（φEnT）利和电子转移效率（φET）随溶剂的极性的变化具有相同的变化趋势