Structure and stability of various states of the EuC and EuC2 molecules

B3LYP level density functional theory （DFT） and multiconfiguration self-consistent-field （MCSCF） level ab initio method calculations have been performed on the basis of relativistic effective core potentials to investigate the nature of EuC and EuC2 molecules. The computed results indicate that the ground states of EuC and EuC2 are ^12∑^＋ and SA2, respectively. Dissociation potential energy curves of the low-lying electronic states of EuC have been calculated using the MCSCF method, and the same level calculation on EuC2 indicates that the dissociation energy of EuC2 of ground state compares well with the available experimental data. The bond characteristic is also discussed using Mulliken populations.     

Pu_3体系的结构与势能函数

用相对论有效原子实势 (RECP)和密度泛函 (B3LYP)方法对Pun(n =2 ,3)体系的结构进行了优化 ,得到了Pu2 和Pu3分子的几何构型分别为D∞h,D3h,其基态分别为 13和 19重态 .在B3LYP RECP水平上得到Pu2 分子的光谱常数ωe=5 2 .3845cm- 1 ,ωe 　χe=0 .0 2 0 1cm- 1 和Pu3分子的谐振频率 (ν1 =5 6 .90 0 7cm- 1 ,ν2 =5 7.1816cm- 1 ,ν3=6 4 0 785cm- 1 )等性质 ,并通过正规方程组和多体展式理论 ,得到了Pu2 ,Pu3的分析势能函数 .     

SeF_n(n=0,+1,+2)基态及SeF低激发态的势能函数与微观性质

为了弄清氟与硒相互拮抗作用的动力学特征,了解SeFn(n=0, 1, 2)分子离子基态和激发态势能函数及微观结构、性质的基本信息是必要的.本文采用6-311 G**基组、B3LYP方法对SeFn(n=0, 1, 2)分子离子基态进行了结构优化和频率计算,用TDB3LYP/6-311 G**含时方法对激发态B12Π、B22Π进行了计算,得到SeFn(n=0, 1, 2)分子离子基态和SeF分子激发态B12Π、B22Π的平衡几何结构、电子状态、谐振频率、偶极矩、离解能De等性质,并在计算出来的一系列单点势能基础上,用正规方程组拟合Murrell-Sorbie(M-S)势能函数,得到相应态的解析势能函数,光谱参数Be、αe、ωe、和ωeχe,由此计算对应的光谱参数和力学性质.理论计算值与相关文献值吻合较好,说明用B3LYP(TDB3LYP)/6-311 G**方法计算SeFn(n=0, 1, 2)分子离子基态和激发态微观结构性质是可行的.     

Pu2分子的结构与势能函数

用相对论有效原子实势(RECP)和密度泛函(B3LYP)方法对Pu2分子的结 构进行了优化,对较高多重性优化得到两个平衡结构,并用Murrell-Sorbie函数导出了基 态两种结构的势能函数和光谱数据。     

AT1受体拮抗剂V8、V12和BZI8结构的量子化学研究

基于Mavromoustakos等人的分子动力学和蒙特卡罗方法结果,本文利用量子化学计算的方法研究了SARTANS(沙坦)类的3种AT1受体拮抗剂(V8、V12和BZI8).在B3LYP/6-31G(d)方法下计算得到了3种药物分子的平衡几何结构,并且采用B3LYP-GIAO/6-311G(d,p)方法计算了这些药物分子中的质子化学位移.通过研究结果可以看出,本文采用量子化学方法计算得到的3种分子的结构与Mavromoustakos等人通过分子动力学和蒙特卡罗方法得到的结果比较接近,计算结果与核磁共振实验的ROE数据和1H谱相比吻合得较好,说明所获得的结构较为合理.另外通过对所获得分子结构进行叠合分析发现,同属SARTANS类的3种分子在联苯咪唑环区域具有很大的相似性.     

Pu3体系的结构与势能函数

用相对论有效原子实势(RECP)和密度泛函(B3LYP)方法对Pu_n(n=2,3)体系的结构进行了优化,得到了Pu₂和Pu₃分子的几何构型分别为D_∞h,D_3h,其基态分别为13和19重态.在B3LYP/RECP水平上得到Pu₂分子的光谱常数ω_e=52.3845cm^-1,ω_e　χ_e=0.02 关键词： 2')" href="#">Pu₂ 3')" href="#">Pu₃ 分析势能函数     

B2H6分子结构和振动频率的DFT研究

用密度泛函理论(DFT)方法研究B2H6分子的几何结构、空间电子密度分布和轨道动量分布以及振动频率等基本性质.取DFT理论中的B3LYP方法和6-311++G* *基组对乙硼烷分子的几何构型进行全优化计算后分析其空间电子密度分布.利用不同方法和基组得到了轨道动量分布.在B3LYP/6-311++G* *基础上对优化后的结构进行正则振动频率分析.     

SiF2自由基与异硫氰酸反应机理的理论研究

采用密度泛函理论B3LYP 方法研究了CF2 自由基与HNCS 的反应机理, 并在B3LYP/6-311＋＋G**水平上对反应物、中间体、过渡态进行了全几何参数优化. 为了得到更精确的能量值,又用CCSD(T)/6-311＋＋G**方法计算了在B3LYP/6-311＋＋G**水平优化后的各个驻点的相对能量. 根据统计热力学及用Winger校正的Eyring过渡态理论,利用自编程序,计算不同温度下低势垒反应的平衡常数和速率常数.计算结果表明单重态的CF2 自由基与HNCS 的反应有6条反应通道,三重态的CF2 自由基与HNCS 的反应有1条反应通道.其中单重态反应通道CF2＋HNCS→IM1→TS1→IM2 HCF2NCS(P1)为主反应通道.     

LaO分子基态X2Σ+势能和平衡几何构型的比较研究

采用三种相对论有效核芯势(LanL2DZ,AREP,RCEP)和价电子基组在MP2、B3LYP、CCSD(T)和QCISD(T)水平上计算了LaO分子基态平衡几何构型和能量,并对各种计算方法在不同基组下的计算效率和精度等进行了综合比较.在此基础上,选用B3LYP方法进行能量扫描得到了LaO分子势能曲线.计算得到的LaO分子基态平衡几何,振动频率、解离能和力常数等均与实验结果吻合.最后对赝势方法计算的误差原因进行了分析.     

AT1受体拮抗剂V8、V12和BZI8结构的量子化学研究

基于Mavromoustakos等人的分子动力学和蒙特卡罗方法结果，本文利用量子化学计算的方法研究了SARTANS（沙坦）类的3种AT1受体拮抗剂（V8、V12和BZI8）. 在B3LYP/6-31G(d)方法下计算得到了3种药物分子的平衡几何结构，并且采用B3LYP-GIAO/6-311G(d，p)方法计算了这些药物分子中的质子化学位移. 通过研究结果可以看出，本文采用量子化学方法计算得到的3种分子的结构与Mavromoustakos等人通过分子动力学和蒙特卡罗方法得到的结果比较接近，计算结果与核磁共振实验的ROE数据和¹H谱相比吻合得较好，说明所获得的结构较为合理. 另外通过对所获得分子结构进行叠合分析发现，同属SARTANS类的3种分子在联苯咪唑环区域具有很大的相似性.     

PPP Hamiltonian for polar polycyclic aromatic hydrocarbons

Total energies of charged states and configurations of different spin multiplicity of two polar non-alternant polycyclic aromatic hdrocarbons (PAH), namely, pentaheptafulvalene and azulene, calculated by means of a Multi-Configurational (MCSCF) method which includes correlation only amongst π orbitals, have been fitted by exact solutions of the Pariser-Parr-Pople (PPP) and the Hubbard Hamiltonians for π electrons. As both molecules are planar, such an approach is in principle feasible. As found in our previous analysis of PAH, PPP fittings are significantly better than those attained with the Hubbard Hamiltonian. In addition, parameters for the Hubbard Hamiltonian are around twice those derived for the PPP model, indicating that parameters are not model independent. Fitted PPP parameters are close to those derived from a similar study of the PAH 2, 5, 8-trihydrogenated phenalene and those originally proposed by Pariser et al. providing further support to a wide applicability of the fitted parameters. Fittings obtained for a MCSCF method that MM298001also includes σσ and σπ correlations (MCSCF/MP2) are slightly less accurate giving an on-site repulsion 10–15% smaller. The accuracy of the fittings further diminishes when parameters are derived from energies obtained by means of a DFT method (B3LYP) with an additional decrease in U of 5–25%. In the latter two cases, parameters have to be considered as effective, accounting for effects of σ orbitals not explicitly included in the model Hamiltonians. Electron affinities, ionization energies and dipole moments, calculated by means of the model Hamiltonians, are compared to those derived from DFT and ab initio methods and, whenever available, to experimental data.     

The treatment of the spin coupling in the bonding of NO to the Ni-doped MgO (100) surface

The results are reported of density functional theory (DFT) and explicitly correlated wave-function (CASSCF and CASPT2) calculations on the bonding of NO with the Ni-doped MgO(100) surface. The surface is represented by means of a cluster of ions embedded in point charges. A comparison is made between unrestricted (spin polarized) and spin restricted approaches. While the geometry of the surface complex is described in quite an accurate way by a spin unrestricted DFT approach, e.g., using the B3LYP functional, the spin distribution does not correspond to that of the real physical situation. In fact, the spin polarized DFT treatment shows three unpaired electrons, two with spin up and one with spin down, while EPR experiments show clearly the existence of a single spin localized on an Ni 3d shell. A spin restricted B3LYP treatment, on the other hand, gives a correct spin distribution and geometry but fails in reproducing the adsorption energy. Other exchange-correlation functionals behave in a similar or even worse way. The CASPT2 results, by contrast, are in substantial agreement with the experiment, showing the importance of treating on the same footing the spin and electron correlation as well as the multi-configuration character of the wavefunction.     

密度泛函理论的不同方法研究小钨团簇的结构与稳定性

分别利用密度泛函理论的LSDA方法和B3LYP 方法在有效核势基组(LanL2DZ)水平上系统地研究了Wn(n=2-14)团簇的结构和稳定性.Wn(n=2-14)团簇全局能量最小的最稳定结构被确定.结果表明:在所讨论的尺寸范围内,LSDA方法比B3LYP 方法得到的平均结合能与一阶能量差分高;大部分Wn(n=2-14) 团簇在两种方法计算下具有相同的几何结构.两种方法计算结果均表明:除了W3自旋多重度为3重态外,其余尺寸团簇的自旋多重度均为单重态;n=7的最低能量结构为双戴帽五边环,而非正十面体,n=8-10为笼型结构,n=11-14为扁长型结构.综合平均结合能、一阶能量差分和二阶能量差分,两种计算方法均表明Wn(n=2-14)团簇具有相同的幻数序列(2,6,9 和 13).     

Multiplet splittings and other properties from density functional theory: an assessment in iron–porphyrin systems

In transition metal compounds with spin states close in energy, the magnitude and sign of the energy splitting calculated with density functional theory depends strongly on the functional used. Therefore we must turn to additional criteria to assess the level of accuracy and reliability of predictions based on this level of theory. We report optimized geometries, total energies, and Mössbauer quadrupole splitting values for low-spin and high-spin, ferric and ferrous model hemes using a variety of gradient-corrected and hybrid functionals. In one model, the iron–porphyrin is axially ligated by two strong-field imidazole ligands [FeP(Im)₂] and has a low-spin ground state. In the other model complex the axial ligands are two weak-field, water molecules [FeP(H₂O)₂], and have a high-spin ground state. Among all the functionals used (UHF, B3LYP, B3LYP*, BLYP, half-and-half, LSDA), the B3LYP hybrid functional most consistently reproduced the experimental geometry, Mössbauer, and spin state data for the two model hemes. Simply gradient-corrected functionals exhibit strong biases towards low spin states, while Hartree–Fock favours strongly high spin states. These findings suggest that for systems with similar characteristics of several accessible electronic spin configurations, it is imperative to include properties other than just the energy in the assessment of the DFT predictions.     

CH3S CH2SH异化反应的理论研究

利用密度泛函理论(DFT)和从头算(ab initio)研究了CH3S←→CH2SH互异化的反应机理．采用HF、B3LYP、MP2理论水平和中等基组6-31(d)，计算了CH3S、CH2SH及其过渡态的结构参数、谐振频率、零点能(ZPF)、总能量和相对能量，并利用B3LYP／6-31(d)的方法计算了反应的内禀反应坐标(IRC)，给出了分子构型和自旋污染沿反应坐标的变化曲线，以及最小能量曲线(MEP)、绝热能量曲线．此外，利用传统过渡态理论(CTST)研究了该互异化反应的速率常数和平衡常数在200～1000K的变化．     

密度泛函理论研究Nb_2Ge_n(n=1～4)团簇

运用密度泛函方法在(U)B3LYP/Lan L2DZ水平上对Nb_2Ge_n(n=1~4)团簇进行了系统的理论研究,得到Nb_2Ge_n(n=1~4)团簇的最低能结构的几何构型和电子性质.优化结果表明:Nb_2Ge_n(n=1~4)团簇最低能结构的自旋多重度均为单重态.团簇最低能结构的电子态与团簇的大小有关.当n为奇数时,团簇的电子态为~1A~1,n为偶数时电子态为1A.通过对计算平均束缚能和分裂能发现:Nb_2Ge_n(n=1~4)团簇中热力学稳定性最强的是Nb_2Ge_2团簇;最弱的是Nb_2Ge_4团簇.自然电荷分布的结果说明Nb_2Ge_n(n=1~4)团簇中当n=1-2时,电子转移正常,而当n=3-4时出现电荷反转现象.同时还研究了HOMOLUMO能隙、磁性和红外光谱.     

铝团簇的电子自旋极化与核自旋-自旋偶合的理论研究

利用传统的密度泛函理论在B3LYP/6-31 G(d)水平上优化了铝簇(Aln ,Aln与Aln-,n=2~9)的几何结构,并利用偶合的微扰的密度泛函理论在B3LYP/6-311 G(3df)水平上计算了核自旋-自旋偶合常数.优化结果表明Aln(n=2~9)中的电子是自旋极化的,与早期的质谱实验一致.核自旋-自旋偶合常数的计算结果表明电子的自旋极化与原子核的自旋取向有密切关系.     

Spin–spin coupling constants in linear substituted HCN clusters

The indirect nuclear spin–spin coupling constants of homogeneous hydrogen-bonded HCN clusters are compared with those of inhomogeneous HCN clusters where one of the terminal HCN molecules is substituted by its isomer HNC and by LiCN. Both the intra- and intermolecular (across the hydrogen bond) coupling constants are calculated for the linear form of the clusters containing up to three molecular monomers using different hybrid DFT functionals. The geometry of the monomers and clusters is optimised at the B3LYP/6-311++G(d,p) level. The effect of substitution by the ionic compound LiCN on the coupling constants of HCN is found to be more pronounced than that by HNC. The Ramsey parameters that form the total spin–spin coupling constants are also analysed individually. Among the four Ramsey parameters, the Fermi Contact term is found to be the dominant contributor to the total coupling constants in most cases. The presence of LiCN in the cluster tends to decrease the intramolecular Fermi Contact values, while HNC increases the same in all dimers and trimers. The contributions of localised molecular orbitals have been analysed for the HCN–HNC cluster to obtain some additional insight about the SSCC transmission mechanism along the coupling pathway.     

The calculation of the energies of phonon normal modes in the alpha-quartz crystal by analyzing the change of forces on atoms in the hydrogen-passivated small cluster of H16Si7O6

The lattice dynamics of quartz have been the subject of considerable experimental and theoretical study. The energies of phonon normal modes in the alpha-quartz crystal calculated by analyzing the change of forces on atoms with their displacements from their equilibrium positions in the hydrogen-passivated small cluster of H₁₆Si₇O₆ in other viewpoint of crystal symmetry. The calculation is done by using the standard density functional theory with the exchange–correlation potential in the B3LYP Hamiltonian form in 6-311+G (2d) basis set. In order to simulate the actual situation of the crystal, we propose an idea to make the atom displacements in accordance with the direction of each actual lattice vibration in the crystal. The normal mode frequencies and the corresponding normal modes assignments were theoretically studied using the Gaussian 03 package. Finally, we obtain that the calculated normal mode frequencies with the B3LYP density functionals are generally consistent with the observed spectra. Results obtained from calculations basing on B3LYP have been reported. The final agreement between experiment and theory appears convincing.