BF分子X1∑+,A1∏和B1∑+电子态的势能函数

使用SAC/SAC-CI和D95++,6-311++g,6-311++g**及D95(d)基组,分别对BF分子的基态X1∑+、第一简并激发态A1∏和第二激发态B1∑+的平衡结构和谐振频率进行优化计算.对所有计算结果进行比较,得出6-311++g**基组为最优基组.运用6-311++g**基组和SAC方法对基态X1∑+,SAC-CI方法对激发态A1∏和B1∑+进行单点能扫描计算,并用正规方程组拟合Murrell-Sorbie函数,得到相应电子态的势能函数解析式,由得到的势能函数计算了与X1∑+,A1∏和B1∑+态相对应的光谱常数,结果与实验数据较为一致.     

LiBe+分子基态及其低电子激发态

采用Davidson校正的多参考组态相互作用方法(MRCI+Q)和多参考平均二次耦合团簇方法(MRAQCC)结合基组ROOS-ANO-TZP得到了LiBe+分子基态(x1∑+)以及六个低电子激发态(a3Σ+,b3∏,A1∑+,B1∏,c3∑+,C1∑+)的势能曲线(PECs).计算结果表明:X1∑+,a3∑+,b3∏和...     

LaH分子的X1∑+态的势能函数

应用原子分子反应静力学原理导出LaH分子的电子状态和可能的离解极限,考虑相对论紧致有效势RCEP(RelativisticCompactEffectivePotential)近似下,用QCISD方法计算了LaH分子基态X1∑+的平衡几何Re和离解能De为2.125A和2.623eV,并在计算出来的一系列单点势能基础上,用正规方程组拟合Murrell-Sorbie(M-S)势能函数,得到相应态的解析势能函数,由此计算对应的光谱参数,其Be、ae、ωe和ωexe的理论值,分别为:3.7333、0.0723、1461.73和21.383cm-1.     

基态UC2分子的结构和势能函数

采用密度泛函理论(DFT)的B3LYP方法和相对论有效原子实势理论模型(RECP),对UC2分子可能的结构进行优化计算,得到UC2分子稳定构型为角形C-U-C(C2v);由微观可逆性原理,判断了UC2分子的离解极限;并且导出了基态UC2分子(X 5B1)的多体项展式势能函数,其势能面等值图展现了C-U-C(C2v)稳定结构;根据势能面等值图,讨论了C+UC(X 3П)反应和U+C2(X 1∑+g)反应的势能面静态特征.     

KH(KD)分子基态(X~1Σ~+)的结构与势能函数研究

利用群论及原子分子反应静力学的有关原理,推导了KH(KD)分子基态X1Σ+的电子态和合理的离解极限.采用Gaussian03程序包中的多种方法和基组,对KH(KD)分子基态X1Σ+的平衡结构和谐振频率进行了优化计算.通过比较计算结果,发现B3LYP为最优方法,6-311g(3df,3pd)为最佳基组.运用优选出的方法和基组对KH(KD)分子基态进行了单点势能扫描,然后分别采用Murrell-Sorbie函数及修正的Murrell-Sorbie+c6函数进行了非线性最小二乘拟合,得到了KH(KD)分子基态的势能函数和相应的光谱常数.计算结果表明,利用修正的Murrell-Sorbie+c6函数计算所得的光谱常数与实验数据吻合得更好.     

PuN基态分子势能函数与热力学函数的理论计算

在Pu的相对论有效原子实势近似和N原子6-311G*全电子基函数下,用密度泛函B3LYP方法计算得到PuN分子基态X6∑+的结构与势能函数、力常数与光谱数据.同时计算得到PuN(g)分子在298 K时的标准生成热力学函数△fH0、△S0和△fG0,分别为-487.239 kJ/mol、95.345 J/mol K和-515.6661 kJ/mol.     

PuN和PuN_2基态分子的结构与势能函数

用相对论有效原子实势和密度泛函理论方法对PuN和PuN_2分子的结构进行优化 ，得到了其平衡几何构型和谐振频率。采用最小二乘法拟合出PuN基态分子的 Murrell-Sorbie解析势能函数，在此基础上推导出光谱数据和力常数，并用多体展 式理论导出PuN_2基态分子的解析势能函数，正确地反映了其平衡构型特性。     

BF分子X~1Σ~+,A~1Π和B~1Σ~+电子态的势能函数

使用SAC/SAC-CI和D95++,6-311++g,6-311++g**及D95(d)基组,分别对BF分子的基态X~1Σ~+、第一简并激发态A~1Π和第二激发态B~1Σ~+的平衡结构和谐振频率进行优化计算.对所有计算结果进行比较,得出6-311++g**基组为最优基组.运用6-311++g**基组和SAC方法对基态X~1Σ~+,SAC-CI方法对激发态A~1Π和B~1Σ~+进行单点能扫描计算,并用正规方程组拟合Murrell-Sorbie函数,得到相应电子态的势能函数解析式,由得到的势能函数计算了与X~1Σ+,A~1Π和B~1Σ~+态相对应的光谱常数,结果与实验数据较为一致.     

MRCI及其近似方法下C2分子X1∑+g, B1△g和B'1∑+g态的势能曲线

运用MRCI及其近似方法计算得到C2分子X1∑+g,B1△+g,和B'1∑+g 3个电子态的势能曲线,并以FCI方法为基准,考察比较各种方法的优缺点.MRCI计算与FCI计算几乎不存在差别,在FCI计算不现实的情况下,MRCI计算足以作为基准考验其他计算方法.多种近似方法均能达到化学精度,并提出一个有效的选择参考组态的方法,提高MRCI计算效率的同时,保证计算的精确性.     

基态Hse分子及Hsex(x=+1,-1) 离子的量子力学计算

用原子-分子反应静力学原理推导出了HSe分子、HSe^x（x=＋1，-1）离子的基态电子状态及其离解极限．在cc—pVDZ水平基础上，用B3LYP方法计算了它们的平衡核间距Re和离解能De．并在计算出来的一系列单点势能基础上，用正规方程组拟合Murrell—Sorbie（M—S）势能函数，得到相应态的解析势能函数，由此计算对应的光谱参数（Be，αe，ωe，ωe，χe）和力学性质，计算结果表明HSe分子、HSe^x（x=＋1，-1）离子可稳定存在．     

Smoothing techniques of global optimization: Distance scaling method in searches for most stable Lennard–Jones atomic clusters

Spatial averaging of the potential energy function facilitates the search for the most stable configuration of a molecular system. Recently some global optimization methods of this kind have been designed in the literature that rely on physical phenomena such as diffusion, wave function evolution in quantum mechanics, Smoluchowski dynamics, evolution in temperature of canonical ensembles, etc. In the present article we highlight the fact that all these methods, when applied to the Gaussian distributions of an ensemble, represent special cases of a set of differential equations involving the spatially averaged potential energy. Their structure suggests that the nature's strategy to cope with the global optimization is robust and differs only in the details in particular applications. The strategy consists of going downhill of the averaged potential energy, removing the barriers, and hunting for low energy regions by a selective increasing of the spatial averaging. In this study we explore the deformation of the potential rather than its averaging. The deformation comes from scaling of atomic distances and reduces the barriers even more effectively than the Gaussian averaging. The position and widths of the Gaussian distribution evolve similarly to the Gaussian density annealing (GDA), but we allow elliptical instead of spherical Gaussians as well as branching of the single trajectory of the system into multiple ones. When the temperature reaches 0 K, one has a number of independent Gaussian distributions, each corresponding to a structure and (usually low) energy of the system. The multiple elliptic-Gaussian distance scaling method has been applied to clusters of argon atoms (N=5,…,31), a system serving usually as a benchmark domain. The method found the global minima for all but three clusters (of very low energy). The procedure is 20 or more times less expensive than the GDA one. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 2040–2049, 1997     

Interaction energy contours of molecules as probed by a test atom

Potential energy contour curves for some simple molecules interacting with a test atom (a Li atom in the ground state) are calculated by theab initio method. Relative hardness of the repulsive molecular surface is found to be anisotropic. This indicates that the outer boundary surface of a molecule changes its form depending on the extent of applied external forces. The attractive nature of the molecular surface is demonstrated to be local. Surroundings of electronegative atoms having nonbonded electron pairs are found to be attractieve to the test atom having both electron-donating and electron-accepting abilities.     

An analytic potential energy function for the amide–amide and amide–water intermolecular hydrogen bonds in peptides

An analytic potential energy function is proposed and applied to evaluate the amide–amide and amide–water hydrogen‐bonding interaction energies in peptides. The parameters in the analytic function are derived from fitting to the potential energy curves of 10 hydrogen‐bonded training dimers. The analytic potential energy function is then employed to calculate the N? H…O?C, C? H…O?C, N? H…OH₂, and C?O…HOH hydrogen‐bonding interaction energies in amide–amide and amide–water dimers containing N‐methylacetamide, acetamide, glycine dipeptide, alanine dipeptide, N‐methylformamide, N‐methylpropanamide, N‐ethylacetamide and/or water molecules. The potential energy curves of these systems are therefore obtained, including the equilibrium hydrogen bond distances R(O…H) and the hydrogen‐bonding energies. The function is also applied to calculate the binding energies in models of β‐sheets. The calculation results show that the potential energy curves obtained from the analytic function are in good agreement with those obtained from MP2/6‐31+G** calculations by including the BSSE correction, which demonstrate that the analytic function proposed in this work can be used to predict the hydrogen‐bonding interaction energies in peptides quickly and accurately. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Measurement of real free energies of transfer of individual ions from water to other solvents with the jet cell

Resolution of the activities of solutions of electrolytes into the individual ionic contributions cannot be carried out rigorously and requires the introduction of extrathermodynamic assumptions which have inherent uncertainties. The most commonly used approaches are basically similar in that they are based on the assumed solvent independence of the difference in the enthalpy or Gibbs energy of transfer of pairs of model solutes, e.g., tetraphenylarsonium and tetraphenylborate ions, or ferricinium ion and ferrocene. In this work we follow an alternative approach pioneered by Parsons involving measurement in the jet (Kenrick) cell of outer-potential differences between solutions of the same electrolyte in two solvents. These potential differences provide the real free energies of transfer of individual ions which, in turn, differ from the usual Gibbs energies of transfer by the work required to transfer the ion through the dipolar layers at the two solvent-gas interfaces. One objective of this work was to improve the reliability of real free energy of transfer measurements, which are experimentally demanding, to within ca. ±0.5 kJ-mol^–1 in order to match typical uncertainties in Gibbs transfer energies of electrolytes. This goal was met, in most instances, by careful evaluation of experimental parameters (particularly jet pressure). A major improvement over previous measurements was made by adding a supporting electrolyte which allowed stable potentials to be obtained at test electrolyte concentrations as low as 10^–4M. Real free energy changes are reported for the transfer of silver ion from water to methanol, ethanol, acetonitrile, propylene carbonate and dimethyl sulfoxide, as well as for the transfer of chloride ion from water to methanol and ethanol. Reliable data of this kind may lead to improved understanding of either the properties of the surfaces of solvents or the interactions of model solutes with solvents, depending on which of the two fields develops most.     

Quantum‐based models of charge‐dependent potential energy surfaces: Three‐state models

Quantum‐based models of how potential energies depend on charge are developed from a three‐state model, at the level of neglecting state‐to‐state overlap. The energy as a function of charge is defined as proposed previously (Valone and Atlas, J Chem Phys 2004, 120, 7262). With this definition, addition of a third state smooths the derivatives of the energy model with respect to charge at integer values of charge that are in the interior of the allowed charge range. These derivatives are related to the chemical potential. At the dissociation limit, this model converges to established limits. Another dependence is proposed that uses two different charges simultaneously. The concepts are illustrated, with calculations on an OH molecule. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Exploring the potential energy surface of retinal, a comparison of the performance of different methods

The ground state structure of retinal has been investigated. We found that DFT and CASSCF produce different results for the bond length alternation in a model system of retinal. Quantum mechanics/molecular mechanics calculations including the closest surrounding amino acids have been performed, using DFT and CASSCF to calculate the structure of retinal in the protein cavity. The planarity of the retinal molecule is affected by the surrounding protein. DFT and CASSCF produce different twist angles. The difference between CASSCF and DFT appears to be related to the positively charged nitrogen of the Schiff base, which leads to different pi-bond orders produced by the two methods.     

Simple molecular wavefunctions with correlation corrections II

The potential curves for the ground state of Li₂ (¹ _g ⁺ ) and FH (¹ _g ⁺ ) are computed. The correlation energy is calculated using a functional of the one- and two-electron density matrices derived from an MC SCF reference wavefunction and is added to the reference energy to obtain a correlated potential curve.     

Theoretical Studies on the Potential Energy Surface and Vibrational Energy Levels of HXeBr

The potential energy surface for the electronic ground state of the HXeBr molecule is constructed from more than 4200 ab initio points calculated using the internally contracted multi-reference configuration interaction method with the Davidson correction （icMRCI ＋ Q）. The stabilities and dissociation barriers are identified from the potential energy surface. The three-body dissociation channel is found to be the dominant dissociation channel for HXeBr. Low-lying vibrational energy levels of HXeBr calculated using the Lanczos algorithm are found to be in good agreement with the available experimental band origins.