稠环共轭分子的能级

我们对分子轨道的差分方程法作了改进,提出了计算稠环共轭分子能级的简化的分子轨道差分方程法。该法主要包括以下几部分:1.不同子集合碳原子的划分:划分的原则是先找出整个分子的对称中心,然后以它为圆心,以中心到不同碳原子的距离为半径画圆,凡在同一圆上的碳原子,属于同一子集合。     

Hamiltonian for rovibrational spectra of “hot” molecules

A variational method based on the determination of rotation as a state with a constant angular momentum has been proposed for calculating rovibrational energy levels of a polyatomic molecule. By using this method, energy level values can be determined for any vibrational state with arbitrary values of vibrational quantum numbers. This makes it possible to calculate the rovibrational energy levels of the “hot” molecules.     

Computation of thermodynamic plasma properties: A simplified approach

An approximate method for calculating the electronic partition functions of atomic systems is reported. The method is based on the idea of combining a multitude of atomic energy levels into two or three grouped levels. Dimensionless formulation suitable for practical calculations is presented. Application to real systems shows that two grouped levels are enough to model hydrogen-like and noble gas atoms, whereas three grouped levels are required to describe atoms with low lying excited states. Two methods for the calculation of degeneracy and energy values of the grouped levels are investigated. Representative mono-atomic LTE plasma properties calculations are reported. The results agree with accurate computations using partition functions that include several thousands energy levels.     

A simple scheme for calculating the energy levels and the spectra of molecular crystals

A simple matrix method for calculating the vibrational and electron energy levels and spectra in molecular crystals is proposed.     

Including anharmonicity in the calculation of rate constants. 1. The HCN/HNC isomerization reaction

A method for calculating anharmonic vibrational energy levels in asymmetric top and linear systems that is based on second-order perturbation theory in curvilinear coordinates is extended to the bound generalized normal modes at nonstationary points along a reaction path. Explicit formulas for the anharmonicity coefficients, x(ij), and the constant term, E0, are presented, and the necessary modifications for resonance cases are considered. The method is combined with variational transition state theory with semiclassical multidimensional tunneling approximations to calculate thermal rate constants for the HCN/HNC isomerization reaction. Although the results for this system are not very sensitive to the choice of coordinates, we find that the inclusion of anharmonicity leads to a substantial improvement in the vibrational energy levels. We also present detailed comparisons of rate constants computed with and without anharmonicity, with various approximations for incorporating tunneling along the reaction path, and with a more practical approach to calculating the vibrational partition functions needed for larger systems.     

势阱中粒子能级与波函数微扰计算的代数递推公式

利用超位力定理(HVT)和Hellmann-Feynman 定理(HFT),导出了由有精确解的势阱的能级值用微扰法直接计算一维势阱的各级近似能级的普遍代数公式,并导出了由能级近似值计算定态波函数近似表达式的代数公式.给出了代数公式具体应用的几个典型一维势阱实例.此法可推广到二维势阱与三维势阱的情形.     

镧系元素f-f跃迁光谱计算

运用不可约张量算符方法,建立了f-f超灵敏跃迁光谱的计算模型,编制了计算f离子^2s+1L~J能级、约化矩阵元、晶体场Stark能级、晶体场态-态跃迁振子强度及模拟吸收光谱的计算机程序.其中,自由f离子能级计算采用包括双电子与叁电子组态相互作用的二三参数模型,晶体场Stark分裂计算采用单电子轨道近似,用Newman角叠加模型计算晶体场参数.f-f跃迁振子强度计算包括静电诱导偶极跃、配体极化偶极跃、振动诱导电偶极跃迁及磁偶极跃迁的贡献.     

双原子分子振转激发态

发展了一种较简便的方法来计算双原子分子的振转激发态。原则上,该方法可以给出同数值积分方法一样精确的结果。作为例子, 给出了H_2的162个振转激发态的准确能级和波函数。     

计算高聚物增塑体系混合自由能新方法

高聚物增塑体系混合自由能与浓度关系 ΔGm=f(ω1) ,不仅反映了增塑体系热力学稳定性好坏 ,而且可以预测溶液各组分间的作用能力 [1,2 ] .通过实验可直接得到增塑体系混合自由能的相关数据 ,其中包括偏蒸气压值测量法、渗透压法、反相气相色谱法及小角 X射线散射法等[3～ 8] .这些方法可以测量增塑体系挥发性组分的偏摩尔化学位值 ,利用下面公式计算Δμ1/RT =ln(p1/p1,S)式中 p1,p1,S是某温度下增塑体系增塑剂的偏蒸气压值和纯增塑剂的饱和蒸气压 .1　现有的高聚物混合体系自由能 ηGm计算方法对于低分子溶液体系各组分的偏摩尔化学位Δ…     

统一统计理论多分割面热反应速率的计算程序

应用GAUSS积分方法，设计了建立在统一统计理论基础上的多分割面量子方法计算热反应速率的计算程序．通量积分N（E，S）基于内禀反应坐标（IRC）的信息采用阶梯函数量子数态法计算，由其微正则变分确定分割面．IRC信息由量子化学从头算方法得到．程序用FORTRAN语言编写．详细阐述了基本理论和程序算法及对H2＋H反应体系的成功应用．该程序可进一步适用于多原子分子的较大体系．     

Intermolecular electrostatic energies using density fitting

A method is presented to calculate the electron-electron and nuclear-electron intermolecular Coulomb interaction energy between two molecules by separately fitting the unperturbed molecular electron density of each monomer. This method is based on the variational Coulomb fitting method which relies on the expansion of the ab initio molecular electron density in site-centered auxiliary basis sets. By expanding the electron density of each monomer in this way the integral expressions for the intermolecular electrostatic calculations are simplified, lowering the operation count as well as the memory usage. Furthermore, this method allows the calculation of intermolecular Coulomb interactions with any level of theory from which a one-electron density matrix can be obtained. Our implementation is initially tested by calculating molecular properties with the density fitting method using three different auxiliary basis sets and comparing them to results obtained from ab initio calculations. These properties include dipoles for a series of molecules, as well as the molecular electrostatic potential and electric field for water. Subsequently, the intermolecular electrostatic energy is tested by calculating ten stationary points on the water dimer potential-energy surface. Results are presented for electron densities obtained at four different levels of theory using two different basis sets, fitted with three auxiliary basis sets. Additionally, a one-dimensional electrostatic energy surface scan is performed for four different systems (H2O dimer, Mg2+-H2O, Cu+-H2O, and n-methyl-formamide dimer). Our results show a very good agreement with ab initio calculations for all properties as well as interaction energies.     

Term crossing in the auxiliary Hamiltonian in calculations of energy-level lower bounds in quantum mechanical systems

We study a method for calculating lower bounds on energy levels that uses an auxiliary Hamiltonian with projection operators. We show that a change in the physical parameters of the problem can lead to the crossing of terms of the auxiliary Hamiltonian with the same symmetry. Such crossings, if not accounted for, will lead to incorrect estimates of the energy levels of quantum systems. As an illustration we consider the lower bound for the ground-state energy of the ion H^–.Translated from Teoreticheskaya i Éksperimental'naya Khimiya, Vol. 24, No. 5, pp. 513–519, September–October, 1988.     

Mixed quantum-classical reaction path dynamics of C2H5F --> C2H4 + HF

A mixed quantum-classical method for calculating product energy partitioning based on a reaction path Hamiltonian is presented and applied to HF elimination from fluoroethane. The goal is to describe the effect of the potential energy release on the product energies using a simple model of quantized transverse vibrational modes coupled to a classical reaction path via the path curvature. Calculations of the minimum energy path were done at the B3LYP/6-311++G(2d,2p) and MP2/6-311++G** levels of theory, followed by energy-partitioning dynamics calculations. The results for the final HF vibrational state distribution were found to be in good qualitative agreement with both experimental studies and quasiclassical trajectory simulations.     

Molecular-orbital coefficients for dinuclear polymethyne dyes in the effective additive parameter method

A method previously proposed for determining the energy structure of a polymethyne dye with any terminal groups is used in considering formulas for the molecularorbital coefficients and the differences in the distribution on the atoms in the polymethyne chain for localized and delocalized energy levels, as well as the accuracy in calculating the molecular-orbital coefficients by means of a finite number of effective additive parameters. It is found that the localized states are important to the electron-density distribution on the chain atoms characteristic of the polymethyne dyes.Translated from Teoreticheskaya i Éksperimental'naya Khimiya, Vol. 22, No. 2, pp. 137–145, March–April, 1986.     

Linear scaling calculation of band edge states and doped semiconductors

Linear scaling methods provide total energy, but no energy levels and canonical wave functions. From the density matrix computed through the density matrix purification methods, we propose an order-N [O(N)] method for calculating both the energies and wave functions of band edge states, which are important for optical properties and chemical reactions. In addition, we also develop an O(N) algorithm to deal with doped semiconductors based on the O(N) method for band edge states calculation. We illustrate the O(N) behavior of the new method by applying it to boron nitride (BN) nanotubes and BN nanotubes with an adsorbed hydrogen atom. The band gap of various BN nanotubes are investigated systematically and the acceptor levels of BN nanotubes with an isolated adsorbed H atom are computed. Our methods are simple, robust, and especially suited for the application in self-consistent field electronic structure theory.     

A new method for estimating the relative binding free energy,derived from a free energy variational principle for the Pim-1-kinase–ligand and FKBP–ligand systems

Journal of Computer-Aided Molecular Design - In this study, a new method is proposed for calculating the relative binding free energy between a ligand and a protein, derived from a free energy...     

过渡金属化合物能谱的理论计算

本文提出了在限制Hartree-Fock方法计算的单电子分子轨道基础上计算分子价电子体系能谱和多电子状态的不可约张量方法, 并利用这种方法, 计算了过渡金属化合物CrCl_6~(3-), CrF_6~(3-), CrO_6~(9-), Cr(H_2O)_6~(3+), Cr(NH_3)_6~(3+), Cr(CN)_6~(3-)的d-d跃迁能谱, 计算结果与实验基本符合, 误差一般为1～2 kcm~(-1), 只有个别情况为3.9 kcm~(-1)。     

Time-dependent treatment of intramolecular energy transfer for HeICl complex

The time-dependent golden wave packet method has been used for calculating the decay widths of vibrational predissociation for HeICl complex in the B state with total angular momentum J=0. This is a good example of intramolecular energy transfer. We examine the dependence of the final rotational distribution (partial decay width) of ICl fragment on the stretching excitation. It is found that computed final rotational distributions are weakly dependent on the vibrational level being excited. Unlike the smoothly varying rotational distribution for lower initial vibrational levels, for higher initial vibrational levels the rotational distribution indicates the very pronounced oscillatory structure. The analysis of the rotational distribution as a function of propagation time reveals the predominant role of the final states interaction in determining the final rotational distribution.     

正八面体过渡金属化合物能谱的理论计算Ⅱ.

本文提出了在非限制Hartree-Fock方法计算的分子轨道基础上,计算分子价电子体系能谱和多电子状态的不可约张量方法。使用这种方法,计算了具有近似正八面体对称性的过渡金属化合物MnBr_6~(4-)、MnCl_6~(4-)、FeCl_6~(3-)、MnF_6~(4-)、Mn(H_2O)_6~(2+)、Co(NH_3)_6~(3+)、NiF_6~(4-)、Ni(H_2O)_6~(2+)、Ni(NH_3)_6~(2+)的d-d跃迁能谱。理论计算结果与实验符合较好,误差为若干k cm~(-1)     

Including anharmonicity in the calculation of rate constants. II. The OH+H2-->H2O+H reaction

A recently developed method for calculating anharmonic vibrational energy levels at nonstationary points along a reaction path that is based on second-order perturbation theory in curvilinear coordinates is combined with variational transition state theory with semiclassical multidimensional tunneling approximations to calculate thermal rate constants for the title reaction. Two different potential energy surfaces were employed for these calculations, an improved version of the author's surface 5 and the WSLFH surface of Wu et al. [J. Chem. Phys. 113, 3150 (2000)]. We present detailed comparisons of rate constants computed for the two surfaces with and without anharmonicity and with various approximations for incorporating tunneling along the reaction path. The results for this system are quite sensitive to the surface employed, the choice of coordinates (curvilinear versus rectilinear), and the inclusion of anharmonicity. A comparison with experiment provides information on the accuracy of these surfaces.