ab initio法と融合したポテンシャル面自動生成に関する研究

结合从头算法自动势面生成研究

• 批准号: 12042232
• 负责人: 石田 俊正
• 金额: $ 1.09万
• 依托单位: Shizuoka University
• 依托单位国家: 日本
• 项目类别: Grant-in-Aid for Scientific Research on Priority Areas (A)
• 财政年份: 2000
• 资助国家: 日本
• 起止时间: 2000 至 2001
• 项目状态: 已结题
• 来源: https://kaken.nii.ac.jp/grant/KAKENHI-PROJECT-12042232/
• 关键词: ポテンシャル画; ab initio法; 反応動力学; 局所内挿法; 最小二乗法; H_4; ポテンシャル面; 分子動力学; H_2O; H_3O

ポテンシャル面は、電子の運動を論じる量子化学と原子核の運動を論じる反応動力学の橋渡しをする重要な役割を果たす。われわれは、エネルギーに関する座標微分を全く使わない方法として、内挿移動最小二乗法(interpolant moving least squares, IMLS)法とShepard内挿法を併用した方法(IMLS/Shepard法)を開発し、いくつかの系に応用してきた。本年度は、H_4系に対してcc-pVTZ基底による4電子15軌道に対するCAS SCF計算でポテンシャルの内挿点のエネルギーを得た。はじめに、28点をとり、100点を加えた。原子核の置換によって24の等価な点があることを考慮すると、点の総数は3,072となる。100点を加えるにあたっては、すべての原子間距離が0.8-7.0auになる配置をランダムに発生し、内挿ポテンシャル面上で、H_2同士を無限に離したときよりも0.3eVだけ高いエネエギーよりも低いエネルギーとなる点をサンプルした。得られたポテンシャル面は、低エネルギーでBoothroydらがH_4に対する、6101点にのぼる多参照CI計算を行い、そのデータに対して、Aguadoらが大局的ポテンシャルフィットを行った結果をよく再現している。また、H_2同士を無限に離したときよりも0.5eVだけ高いエネエギーをサンプルの閾値に使ったにもかかわらず、大局的に6eV程度のH_2-H_2の障壁を定性的に再現している。Aguadoらのフィットにおいて誤差が1kcal/mol程度と報告されているが、われわれの方法は局所内挿法に基づいているため、同じデータを用いた場合、データに用いた計算点における誤差は全くなくなるので、高精度のポテンシャルフィットにも適用性は高いと考えられる。以上のほかにBayesianの定理に基づいた重み関数を用いた場合について計算を行った。

ポ テ ン シ ャ ル は, electronic の movement theory of を じ る と nuclei を の movement theory of quantum chemistry じ る anti 応 dynamics の crossing し を す る important な "を cut fruit た す. わ れ わ れ は, エ ネ ル ギ ー に masato す る coordinate differential を the く has わ な い method と し て, scions in moving least squares method (interpolant moving further squares, IMLS) method と Shepard を scions method and used in し た method (IMLS/Shepard method) を open 発 し, い く つ か の is に 応 with し て き た. This year's は, H_4 に し seaborne て cc - pVTZ basal に よ る 4 electronic 15 tracks に す seaborne る CAS SCF calculation で ポ テ ン シ ャ ル の の scions point in エ ネ ル ギ ー を た. Youdaoplaceholder0, 28 points をと をと, 100 points を plus えた. Nuclear <s:1> permutation によって24 価な and other 価な points がある がある とを とを consider すると, point <s:1> 総 number <e:1> 3,072となる. 100 points of を え る に あ た っ て は, す べ て の atomic distance between 0.8 7.0 au が に な る configuration を ラ ン ダ ム に 発 し, within the scions ポ テ ン シ ャ で, H_2 ル surface with the infinite "を に from し た と き よ り も 0.3 eV だ け high い エ ネ エ ギ ー よ り も low い エ ネ ル ギ ー と な る point を サ ン プ ル し た. Have ら れ た ポ テ ン シ ャ ル は, low エ ネ ル ギ ー で Boothroyd ら が H_4 に す seaborne る, 6101 に の ぼ る computed を い, with reference to CI more そ の デ ー タ に し seaborne て, Aguado ら が the general situation of ポ テ ン シ ャ ル フ ィ ッ ト を line っ た results を よ く reappearance し て い る. ま た, H_2 with infinite "を に し た と き よ り も 0.5 eV だ け high い エ ネ エ ギ ー を サ ン プ ル の threshold numerical に make っ た に も か か わ ら ず 6 eV degree, the general situation of に の H_2 - H_2 の barrier を qualitative に reappearance し て い る. Aguado ら の フ ィ ッ ト に お い て error が と degree 1 kcal/mol さ れ て い る が, わ れ わ れ の way は bureau within the scions subgrade に づ い て い る た め, with じ デ ー タ を with い た occasions, デ ー タ に with い た calculation point に お け る error は full く な く な る の で, high precision の ポ テ ン シ ャ ル フ ィ ッ ト に も applicability Youdaoplaceholder7 gao と exam えられる. Above の ほ か に Bayesian の theorem に base づ い た heavy み masato number を with い た occasions に つ い て count を っ た.

项目成果

J.Aihara and T.Ishida: "Kinetic Stability of Azaborafullerene C_<58>BN as Predicted by the Bond Resonance Energy Model"Bull.Chem.Soc.Jpn.. 73(8). 1791-1795 (2000)

J.Aihara 和 T.Ishida：“键合共振能量模型预测的 Azaborafullerene C_ 58 > BN 的动力学稳定性”Bull.Chem.Soc.Jpn.. 73(8)。

T.Ishida, K.Ohno: "Valence Electronic Force and Core Electronic Force. Significance of Core Electrons on Chemical Binding"J.Molecular Structure(THEOCHEM). 574. 145-152 (2001)

T.Ishida、K.Ohno：“价电子力和核心电子力。核心电子对化学结合的意义”J.分子结构（THEOCHEM）。

Y.Tsuchiya,T.Ishida,H.Katagiri: "A full configuration interaction calculation based on Slater determinants. Application to AIH spectroscopic constants"JCPE Journal. 12(4). 301-308 (2000)

Y.Tsuchiya，T.Ishida，H.Katagiri：“基于 Slater 行列式的完整构型相互作用计算。在 AIH 光谱常数中的应用”JCPE 期刊。

T.Ishida, H.Katagiri: "An ab initio molecular orbital and dynamics study on Penning ionization of Ar with He metastables(He^*(2^1S,2^3S))into the spin-orbit Ar^+(^2P_<3/2>, ^2P_<1/2>) states"J.Phys.Chem. A. 105. 9379-9387 (2001)

T.Ishida、H.Katagiri：“从头开始分子轨道和动力学研究 Ar 与 He 亚稳态 (He^*(2^1S,2^3S)) 进入自旋轨道 Ar^ (^2P_<

K.Fujine, T.Ishida, J.Aihara: "Localization Energies for Graphite and Fullerenes"Phys.Chem.Chem.Phys. 3. 3917-3919 (2001)

K.Fujine、T.Ishida、J.Aihara：“石墨和富勒烯的局域化能量”Phys.Chem.Chem.Phys。