SCF perturbation theory for unimolecular reactions

A CNDO/2 SCF perturbation theory is presented for interpreting the form of CNDO/2 potential energy surfaces of unimolecular reactions. The analysis is performed by calculating the energy change E arising from a distortion of the molecular geometry along the reaction coordinate. E is decomposed into different perturbational contributions which are appropriate for an interpretation of the perturbation energy E. Moreover, E is resolved into energy parts arising from a single occupied orbital and contributions due to pairwise orbital interactions. In this way one evaluates numerically how the form of the occupied and unoccupied orbitals determines the magnitude of E. If the distortion occurs along a definite symmetry coordinate, group-theoretical arguments can be applied to discuss the magnitude of characteristic components of the perturbation energy. The SCF perturbation theory is used to analyze the isomerization of ethylene, cis-2-butene and cis-2-butenenitrile.This work was partially supported by Nato-Grant No. 1072     

Improving the efficiency and reliability of free energy perturbation calculations using overlap sampling methods

A challenge in free energy calculation for complex molecular systems by computer simulation is to obtain a reliable estimate within feasible computational time. In this study, we suggest an answer to this challenge by exploring a simple method, overlap sampling (OS), for producing reliable free-energy results in an efficient way. The formalism of the OS method is based on ensuring sampling of important overlapping phase space during perturbation calculations. This technique samples both forward and reverse free energy perturbation (FEP) to improve the free-energy calculation. It considers the asymmetry of the FEP calculation and features an ability to optimize both the precision and the accuracy of the measurement without affecting the simulation process itself. The OS method is tested at two optimization levels: no optimization (simple OS), and full optimization (equivalent to Bennett's method), and compared to conventional FEP techniques, including the widely used direct FEP averaging method, on three alchemical mutation systems: (a) an anion transformation in water solution, (b) mutation between methanol and ethane, and (c) alchemical change of an adenosine molecule. It is consistently shown that the reliability of free-energy estimates can be greatly improved using the OS techniques at both optimization levels, while the performance of Bennett's method is particularly striking. In addition, the efficiency of a calculation can be significantly improved because the method is able to (a) converge to the right answer quickly, and (b) work for large perturbations. The basic two-stage OS method can be extended to admit additional stages, if needed. We suggest that the OS method can be used as a general perturbation technique for computing free energy differences in molecular simulations.     

Direct calculation of ionization potentials of closed-shell atoms and molecules

Vertical ionization potentials, electron affinities and information about quasi-particles can be obtained by using the technique of the single-particle propagator. The expansion of the self-energy part up to third order perturbation theory can be evaluated numerically, but does not lead, in most cases, to satisfying results. A theoretical and numerical analysis of the diagrammatic expansion of the self-energy part requires the introduction of a renormalized interaction and renormalized hole and particle lines.     

Theoretical study of molecular conduction: I. Effective Green's function based on perturbation theory

There have been many experimental and theoretical studies on molecular conduction, as it is a fundamental parameter in the study of molecular‐scale electronics. We have investigated the features of molecular conduction using a Green's function method, which has often been used to solve problems in quantum transport and is also effective in elucidating electron transport in molecules. We have obtained the novel effective Green's functions, including the first‐order energy corrections, by accommodating the self‐energy of the electrodes as perturbation terms. Although these approximate Green's functions only provide information on the first‐order energy corrections, they can involve the elementary properties of molecular conduction. We propose a scheme for the analysis of the relations between molecular orbitals and their roles in molecular conduction and present analytical calculations for normal and cyclic polyenes. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

The electronic spectrum of VCr

The electronic spectrum of VCr has been studied using the complete-active-space self-consistent field complete-active-space second-order perturbation theory approach. Potential-energy curves for 12 electronic states have been computed. Transition energies, with respect to the X² ground state, for some of the calculated electronic states are (with possible experimental values within parentheses) 0.53 eV (0.56) for A²⁺, 1.03 eV (1.14) for A⁴, 1.20 eV (1.14) for B², 1.45 eV (1.51) for B⁴, 1.60 eV (1.51, 1.78) for C², and 1.61 eV (1.63) for A^4–. AcknowledgementsThe research reported in this communication was supported by a grant from the Swedish Natural Science Research Council (NFR). S. Alex and D.G. Leopold are acknowledged for providing results from their negative ion photoelectron spectroscopy on VCr.Contribution to Björn Roos Honorary Issue     

HCN(HNC)与NH3, H2O和HF分子间相互作用的理论研究

在MP2/aug-cc-pVTZ水平上, 对HCN(HNC)与NH3, H2O和HF分子间可能存在的氢键型复合物进行了全自由度能量梯度优化, 通过在相同水平上的频率验证分析发现了稳定的分子间相互作用形式是HCN(HNC)作为质子供体或作为质子受体形成的复合物. 基组重叠误差对总相互作用能的影响均小于3.34 kJ/mol. 通过自然键轨道(NBO)分析, 研究了单体和复合物中的原子电荷和电荷转移对分子间相互作用的影响. 对称性匹配微扰理论(SAPT, Symmetry Adapted Perturbation Theory)能量分解结果表明, 在分子间相互作用中, 静电作用与诱导作用占主导地位, 而诱导作用与复合物的电荷转移之间具有良好的正相关性.     

Theoretical study on the intramolecular proton transfer reactions of 3-methyl-5-hydroxyisoxazole and its water complexes

The optimized structures and proton transfer reactions of 3-methyl-5-hydroxyisoxazole and its water complexes (3-M-5-HIO · (H₂O)_n · (n = 0–3)) were computed at B3LYP and MP2 theoretical level. The results indicates that 3-M-5-HIO has four isomers (Ecis, Etrans, K1 and K2), and the keto tautomer, and K2 is the most stable isomer in the gas phase. Hydrogen bonding between 3-M-5-HIO and the water molecules can dramatically lower the barrier by the concerted transfer mechanism. Ecis · (H₂O)₃ → K1 · (H₂O)₃ and Ecis · (H₂O)₂ → K2 · (H₂O)₂ is found to be very efficient. Comparing with the proton transfer mechanism of 5-HIO shows that the methyl substitution prevents the intramolecular proton transfer.     

Estimates of the tail of the stationary density function of certain nonlinear autoregressive processes

We consider the Markov chainX _n+1=T(X _n )+ _n , where { _n ;n1} is a ^d -valued random sequence of independent identically distributed random variables, and the functionT: ^{d d} is measurable and satisfies a suitable growth condition. Under certain conditions involvingT and the probability distribution of _n , we show that this Markov chain is ergodic. Moreover, we obtain sharp upper bounds for the tail of the corresponding stationary probability density function. In our proofs, we make use of the Leray-Schauder fixed-point theorem.     

Novel perturbation expansion for the Langevin equation

We discuss the randomly driven systemdx/dt= -W(x) +f(t), wheref(t) is a Gaussian random function or stirring force withf(t)f(t)=2(t–t), andW(x) is of the formgx ¹⁺². The parameter is a measure of the nonlinearity of the equation. We show how to obtain the correlation functionsx(t)f(t)···x(t( n)) _f as a power series in. We obtain three terms in the expansion and show how to use Padé approximants to analytically continue the answer in the variable. By using scaling relations, we show how to get a uniform approximation to the equal-time correlation functions valid for allg and.