由新Hession矩阵所构造的简谐势能函数及其应用

经过简单的矩阵变换，得到了新Ｈｅｓｓｉｏｎ矩阵。由它构造的多原子分子简谐振势函数可以准确地拟合实验频率，在此基础上，采用准经典轨迹法，模拟了ＳｉＨ４、Ｈ２Ｏ等体系分子内能量转移的动力学过程，得到的结果与实验及理论一致。     

生物体系中电子转移机理的理论研究Ⅰ: 一种计算电子转移 的方法及在螺环化 合物中的应用

应用"相应轨道变换"和"广义"密度矩阵的方法,向MOPAC程序包中加入了新的功能,使其能处理电子转移反应中的部分参数。然后用此程序包中AM1方法对具有螺环结构的分子进行处理,计算了化合物在不同外加电场强度下的势能面、反应热ΔE,重组能λ及电子转移矩阵元V~A~B,结果表明,λ,V~A~B受外加电场的影响较小,而ΔE则与之成正比。对标题化合物1的计算结果也同abinitio法的结果进行了比较,发现其变化趋势完全一致,这说明本方法在计算电子转移方面是可靠的。与abinitio方法相比,本程序不仅适用于计算较大体系(如生物大分子),而且还具有速度快,耗时少的优点。     

Formamide dimers: a computational and matrix isolation study

The dimerization of formamide (FMA) has been investigated by matrix isolation spectroscopy, static ab initio calculations, and ab initio molecular dynamics (AIMD) simulations. Comparison of the experimental matrix IR spectra with the ab initio calculations reveals that two types of dimers A and C are predominantly formed, with two and one strong NH...O hydrogen bonds, respectively. This is in accordance with previously published experiments. In addition, there is also experimental evidence for the formation of the thermally labile dimer B after deposition of high concentrations of FMA in solid xenon. The AIMD simulations of the aggregation process show that in all cases dimer C is initially formed, but rearrangement to the more stable doubly hydrogen-bonded structures A or B occurs for a fraction of collisions on the sub-picosecond time scale.     

Predicting shielding constants in solution using gauge invariant atomic orbital theory and the effective fragment potential method

A method to approximate ab initio shielding constants is presented, in which the ab initio density matrix is replaced in the gauge invariant atomic orbital formalism with the density matrix resulting from an effective fragment potential calculation. The resulting first-order density matrix is then iterated to self-consistency. The method is compared with fully ab initio gauge invariant atomic orbital restricted Hartree-Fock calculations on hydrogen chloride, water, and ammonia solutes with up to nine solvent water molecules using the 6-31G, 6-31G(d,p), and 6-31+G(d,p) basis sets. Using the 6-31G(d,p) basis sets, the average of the average absolute deviations for the three environments tested is 0.34 ppm. This is sufficiently accurate to allow for the identification of specific (1)H nuclei in a solvated molecule when the chemical shift between nuclei is not less than 1 ppm. The success of the method at this level of approximation is due to a cancellation of errors between the paramagnetic and diamagnetic terms of the shielding constant: the diamagnetic term is underestimated by roughly the same amount that the paramagnetic term is overestimated.     

Electronic coupling matrix elements from charge constrained density functional theory calculations using a plane wave basis set

We present a plane wave basis set implementation for the calculation of electronic coupling matrix elements of electron transfer reactions within the framework of constrained density functional theory (CDFT). Following the work of Wu and Van Voorhis [J. Chem. Phys. 125, 164105 (2006)], the diabatic wavefunctions are approximated by the Kohn-Sham determinants obtained from CDFT calculations, and the coupling matrix element calculated by an efficient integration scheme. Our results for intermolecular electron transfer in small systems agree very well with high-level ab initio calculations based on generalized Mulliken-Hush theory, and with previous local basis set CDFT calculations. The effect of thermal fluctuations on the coupling matrix element is demonstrated for intramolecular electron transfer in the tetrathiafulvalene-diquinone (Q-TTF-Q(-)) anion. Sampling the electronic coupling along density functional based molecular dynamics trajectories, we find that thermal fluctuations, in particular the slow bending motion of the molecule, can lead to changes in the instantaneous electron transfer rate by more than an order of magnitude. The thermal average, (<|H(ab)|(2)>)(1/2)=6.7 mH, is significantly higher than the value obtained for the minimum energy structure, |H(ab)|=3.8 mH. While CDFT in combination with generalized gradient approximation (GGA) functionals describes the intermolecular electron transfer in the studied systems well, exact exchange is required for Q-TTF-Q(-) in order to obtain coupling matrix elements in agreement with experiment (3.9 mH). The implementation presented opens up the possibility to compute electronic coupling matrix elements for extended systems where donor, acceptor, and the environment are treated at the quantum mechanical (QM) level.     

分子间相互作用的研究——-价轨道模型势从头算能量分解法

研究分子间相互作用是了解液体,固体性质和结构以及气体性质的关键,也是研究化学和生物化学催化机制及化学反应途径的重要方面.因此,近几年来这个领域的理论和实验研究引起了人们广泛的兴趣并取得了长足的进步.特别值得指出的是Morokuma等提出的基于单行列式从头算的能量分解法,较好地解决了很多体系中分子间相互作用的本质问题,受到了理论化学界的普遍重视.但是这种方法很费计算机时间,对于较大的体系,特别是含有重原子的体系,应用受到了限制.因此简化这种方法,使它能够比较容     

Limitations of the diatomics-in-molecules (DIM) method with neglect of overlap

The approximations made in the diatomics-in-molecules (DIM) method are examined by comparing the matrix elements of the DIM calculation with those of the ab initio calculation. It was found that if the overlap between atomic orbitals is small, the ab initio matrix elements can be expressed as the sums of DIM matrix elements and three-center terms ignored in the DIM method. The zero overlap of the atomic orbital approximation used in the DIM method makes the three-center terms simpler and apparently smaller which justifies the applicability of this approximation to the DIM method. Results of calculations show that deviations of the DIM values from accurate values are particularly large in the linear HXH configuration where X is a many-electron atom. It is suggested that a supplementary term be added to the DIM energy (DIM-SC) to take into account the neglected three-center terms.     

Free energy perturbation study of water dimer dissociation kinetics

An efficient approach is described for using accurate ab initio calculations to determine the rates of elementary condensation and evaporation processes that lead to nucleation of aqueous aerosols. The feasibility of the method is demonstrated in an application to evaporation rates of water dimer at 230 K. The method, known as ABC-FEP (ab initio/classical free energy perturbation), begins with a calculation of the potential of mean force for the dissociation (evaporation) of small water clusters using a molecular dynamics (MD) simulation with a model potential. The free energy perturbation is used to calculate how changing from the model potential to a potential calculated from ab initio methods would alter the potential of mean force. The difference in free energy is the Boltzmann-weighted average of the difference between the ab initio and classical potential energies, with the average taken over a sample of configurations from the MD simulation. In principle, the method does not require a highly accurate model potential, though more accurate potentials require fewer configurations to achieve a small sampling error in the free energy perturbation step. To test the feasibility of obtaining accurate potentials of mean force from ab initio calculations at a modest number of configurations, the free energy perturbation method has been used to correct the errors when some standard models for bulk water (SPC, TIP4P, and TIP4PFQ) are applied to water dimer. To allow a thorough exploration of sampling issues, a highly accurate fit to results of accurate ab initio calculations, known as SAPT-5s, as been used a proxy for the ab initio calculations. It is shown that accurate values for a point on the potential of mean force can be obtained from any of the water models using ab initio calculations at only 50 configurations. Thus, this method allows accurate simulations of small clusters without the need to develop water models specifically for clusters.     

烷基自由基β位裂解反应类反应势垒与速率常数的精确计算

提出反应类等键方法并用于高温燃烧机理中一类重要反应——烷基自由基β位裂解反应的反应势垒和速率常数的精确校正计算. 通过10种不同从头算水平对类反应中5个代表反应的反应势垒的计算发现, 用反应类等键反应方法和直接从头算方法获得的5 个代表反应的反应势垒最大绝对偏差的平均值分别为5.32 和16.16 kJ·mol^-1, 表明反应类等键反应方法计算的反应势垒对不同水平从头算方法的依赖性小, 可在较低从头算水平计算得到精确的反应势垒, 解决大分子体系反应势垒的精确计算问题. 此外应用反应类等键反应方法在BHandHLYP/cc-pVDZ 从头算水平计算了3 个代表反应的速率常数, 并与文献报道的实验值进行了比较, 其在500-2000 K温度区间内计算速率常数与实验速率常数中较大值与较小值的比值k_max/k_min的平均值为1.67, 最大值也仅有2.49. 表明应用反应类等键反应方法在较低从头算水平即可对同类反应的速率常数进行精确计算.最后在BHandHLYP/cc-pVDZ从头算水平用反应类等键反应方法计算了13个烷基自由基β位裂解反应的速率常数.     

Research in the method of large molecular calculations utilizing transferability of LMO

The localized molecular orbitals of some saturated hydrocarbons and their derivatives have been formed using ab initio method and M. P. [1–2] localization procedure. Two models, SLMO and ELMO , a set of parameters of LMO Fock matrix elements, and a technique called “Group Effect” are proposed. Based on these, we developed a procedure to simulate the ab initio calculations on large molecules. Some test calculations have been done. The results are compared with those of the ab initio method. In general, absolute errors of orbital energies are about 10^?3 a.u., and the relative errors of total energies are about 10^?4. For the original applications, we applied this procedure to some large systems of alkane and their derivatives as well as three Crown-ether compounds. Satisfactory results are obtained.     

Spin-orbit coupled potential energy surfaces and properties using effective relativistic coupling by asymptotic representation

A new method has been reported recently [H. Ndome, R. Welsch, and W. Eisfeld, J. Chem. Phys. 136, 034103 (2012)] that allows the efficient generation of fully coupled potential energy surfaces (PESs) including derivative and spin-orbit (SO) coupling. The method is based on the diabatic asymptotic representation of the molecular fine structure states and an effective relativistic coupling operator and therefore is called effective relativistic coupling by asymptotic representation (ERCAR). The resulting diabatic spin-orbit coupling matrix is constant and the geometry dependence of the coupling between the eigenstates is accounted for by the diabatization. This approach allows to generate an analytical model for the fully coupled PESs without performing any ab initio SO calculations (except perhaps for the atoms) and thus is very efficient. In the present work, we study the performance of this new method for the example of hydrogen iodide as a well-established test case. Details of the diabatization and the accuracy of the results are investigated in comparison to reference ab initio calculations. The energies of the adiabatic fine structure states are reproduced in excellent agreement with reference ab initio data. It is shown that the accuracy of the ERCAR approach mainly depends on the quality of the underlying ab initio data. This is also the case for dissociation and vibrational level energies, which are influenced by the SO coupling. A method is presented how one-electron operators and the corresponding properties can be evaluated in the framework of the ERCAR approach. This allows the computation of dipole and transition moments of the fine structure states in good agreement with ab initio data. The new method is shown to be very promising for the construction of fully coupled PESs for more complex polyatomic systems to be used in quantum dynamics studies.     

Nonempirical wave functions for very large molecules. I. The PRDDO/M method

We describe enhancements to the method of partial retention of diatomic differential overlap (PRDDO ). The new method, denoted PRDDO/M , employs a basis set of not quite orthogonal atomic orbitals (NQOAOS ) and utilizes sparse matrix techniques to greatly increase the computational efficiency for large molecules. Other modifications, including a complete reparametrization of the method against ab initio STO -3G calculations and implementation of integral screening/damping algorithms, are described. The method is an order of magnitude or more faster than are STO 3G single-point calculations using modern ab initio codes, with little loss in accuracy. © 1996 John Wiley & Sons, Inc.     

Validation of intermolecular transfer integral and bandwidth calculations for organic molecular materials

We present an interpretation of the intermolecular transfer integral that is independent from the origin of the energy scale allowing convergence studies of this important parameter of organic molecular materials. We present extensive numerical studies by using an ethylene pi dimer to investigate the dependence of transfer integrals on the level of theory and intermolecular packing. Transfer integrals obtained from semiempirical calculations differ substantially from one another and from ab initio results. The ab initio results are consistent across all the levels used including Hartree-Fock, outer valence Green's function, and various forms of density functional theory (DFT). Validation of transfer integrals and bandwidths is performed by comparing the calculated values with the experimental values of tetrathiafulvalene-tetracyanoquinodimethane (TTF-TCNQ), bis[1,2,5]thiadiazolo-p-quinobis(1,3-dithiole), (BTQBT) K-TCNQ, and hexagonal graphite. DFT in one of its presently popular forms, such as Perdew-Wang functionals (PW91), in combination with sufficient basis sets provides reliable transfer integrals, and therefore can serve as a basis for energy band calculations for soft organic materials with van der Waals gaps.     

A hierarchic sparse matrix data structure for large-scale Hartree-Fock/Kohn-Sham calculations

A hierarchic sparse matrix data structure for Hartree-Fock/Kohn-Sham calculations is presented. The data structure makes the implementation of matrix manipulations needed for large systems faster, easier, and more maintainable without loss of performance. Algorithms for symmetric matrix square and inverse Cholesky decomposition within the hierarchic framework are also described. The presented data structure is general; in addition to its use in Hartree-Fock/Kohn-Sham calculations, it may also be used in other research areas where matrices with similar properties are encountered. The applicability of the data structure to ab initio calculations is shown with help of benchmarks on water droplets and graphene nanoribbons.     

On the influence of semirigid environments on proton transfer along molecular chains

The dynamics of proton transfer along ammonia chains (chemical composition N(x)H(+)(3x+1), x=2, 4, and 6) in a constraining environment is investigated by ab initio molecular dynamics simulations. A carbon nanotube of defined length and diameter is used as an idealized constraining environment such that the ammonia chain is forced to maintain its quasilinear geometry. It is found that, although the energetics of proton transport shows considerable energetic barriers, proton translocation along the wire is possible at finite temperature for all chain lengths studied. The proton transport involves rotational reorientation of the proton-carrying ammonia molecule. High level ab initio calculations (MP2/aug-cc-pVTZ) yield barriers for internal rotation of 9.1 kcal/mol for NH(4) (+)-NH(3) and 11.7 kcal/mol for OH(3) (+)-OH(2), respectively. The infrared spectrum calculated from the dipole-dipole autocorrelation function shows distinct spectral features in the regions (2000-3000 cm(-1)) where the NHN proton transfer mode is expected to absorb. Assigning moderate opposite total charges between 0.002 and 0.2e to the carbon atoms at the end caps of the nanotube leads to a considerable speedup of the proton transfer.     

Ab initio ground and excited state potential energy surfaces for NO-Kr complex and dynamics of Kr solids with NO impurity

The intermolecular potentials for the NO(X 2Pi)-Kr and NO(A 2Sigma+)-Kr systems have been calculated using highly accurate ab initio calculations. The spin-restricted coupled cluster method for the ground 1 2A' state [NO(X 2Pi)-Kr] and the multireference singles and doubles configuration interaction method for the excited 2 2A' state [NO(A 2Sigma+)-Kr], respectively, were used. The potential energy surfaces (PESs) show two linear wells and one that is almost in the perpendicular position. An analytical representation of the PESs has been constructed for the triatomic systems and used to carry out molecular dynamics (MD) simulations of the NO-doped krypton matrix response after excitation of NO. MD results are shown comparatively for three sets of potentials: (1) anisotropic ab initio potentials [NO molecule direction fixed during the dynamics and considered as a point (its center of mass)], (2) isotropic ab initio potentials (isotropic part in a Legendre polynomial expansion of the PESs), and (3) fitted Kr-NO potentials to the spectroscopic data. An important finding of this work is that the anisotropic and isotropic ab initio potentials calculated for the Kr-NO triatomic system are not suitable for describing the dynamics of structural relaxation upon Rydberg excitation of a NO impurity in the crystal. However, the isotropic ab initio potential in the ground state almost overlaps the published experimental potential, being almost independent of the angle asymmetry. This fact is also manifested in the radial distribution function around NO. However, in the case of the excited state the isotropic ab initio potential differs from the fitted potentials, which indicates that the Kr-NO interaction in the matrix is quite different because of the presence of the surrounding Kr atoms acting on the NO molecule. MD simulations for isotropic potentials reasonably reproduce the experimental observables for the femtosecond response and the bubble size but do not match spectroscopic results. A general overall view of the results suggests that, when the Kr-NO interaction takes place inside the matrix, potentials are rather symmetric and less repulsive than those for the triatomic system.     

A distance-dependent parameterization of the extended Hubbard model for conjugated and aromatic hydrocarbons derived from stretched ethene

The Hubbard model, which is widely used in physics but is mostly unfamiliar to chemists, provides an attractive yet simple model for chemistry beyond the self consistent field molecular orbital approximation. The Hubbard model adds an effective electron-electron repulsion when two electrons occupy the same atomic orbital to the familiar Hückel Hamiltonian. Thus it breaks the degeneracy between excited singlet and triplet states and allows an explicit treatment of electron correlation. We show how to evaluate the parameters of the model from high-level ab initio calculations on two-atom fragments and then to transfer the parameters to large molecules and polymers where accurate ab initio calculations are difficult or impossible. The recently developed MS-RASPT2 method is used to generate accurate potential energy curves for ethene as a function of carbon-carbon bond length, which are used to parameterize the model for conjugated hydrocarbons. Test applications to several conjugated/aromatic molecules show that even though the model is very simple, it is capable of reasonably accurate predictions for bond lengths, and predicts molecular excitation energies in reasonable agreement with those from the MS-RASPT2 method.     

Inner Sphere Reorganization Energy for the Electron Transfer Reactions between M-C_6H_6 and M~ -C_6H_6 Complexes in Gaseous phase: An ab initio Computation

IntroductionTheelectrontransfer(ET)reactionisanfundamentalchemicalprocessandhasbecomethefocusoftheoreticalandexperimentalstudies1-3.AsuitableapproachtowardunderstandingtheETreactivitypatternsistoseparatethevariouscontributionsintoexponentialandpre-exponentialfactors.Theexponentialparts,i.e.,theFrank-Condonfactor,isacompositeofthestructuralreorganizationofthereactingsystemimmediate1yaftertheelectrontransfered.Basedonthecharacteristicsofthereactingsystem'sstructure,thereorganizationenergy(RE…     

The fluorescence of 5-cyano-2-(1-pyrrolyl)-pyridine (CPP) in different solvents and in solid argon: An experimental and theoretical study

The fluorescence spectrum of 5-cyano-2-(1-pyrrolyl)-pyridine (CPP) was measured in several solvents as well as in an argon matrix. Based on comparison with other compounds and on ab initio calculations it is proposed that the fluorescence in the argon matrix and in non-polar solvents is due to two electronic excited states: one is of locally excited nature, the other a charge transfer (CT) state. In polar solvents the spectrum is dominated by the CT emission. The photo-physical behavior of CPP is discussed by comparison with that of other molecules exhibiting dual fluorescence and in view of a recent model developed for the benzene analog.