JP2008210306A - 化合物のシミュレーション方法 - Google Patents
化合物のシミュレーション方法 Download PDFInfo
- Publication number
- JP2008210306A JP2008210306A JP2007048381A JP2007048381A JP2008210306A JP 2008210306 A JP2008210306 A JP 2008210306A JP 2007048381 A JP2007048381 A JP 2007048381A JP 2007048381 A JP2007048381 A JP 2007048381A JP 2008210306 A JP2008210306 A JP 2008210306A
- Authority
- JP
- Japan
- Prior art keywords
- calculation
- molecular orbital
- calculating
- calculated
- molecular
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 150000001875 compounds Chemical class 0.000 title claims description 21
- 238000004088 simulation Methods 0.000 title description 7
- 238000004364 calculation method Methods 0.000 claims abstract description 90
- 238000004776 molecular orbital Methods 0.000 claims abstract description 66
- 238000004965 Hartree-Fock calculation Methods 0.000 claims abstract description 18
- 230000008569 process Effects 0.000 claims abstract description 13
- 238000006073 displacement reaction Methods 0.000 claims description 17
- 238000007796 conventional method Methods 0.000 claims description 10
- 238000004836 empirical method Methods 0.000 claims 3
- 238000005182 potential energy surface Methods 0.000 abstract 1
- 239000011159 matrix material Substances 0.000 description 12
- 238000003078 Hartree-Fock method Methods 0.000 description 4
- 230000004913 activation Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 125000001153 fluoro group Chemical group F* 0.000 description 2
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000004057 DFT-B3LYP calculation Methods 0.000 description 1
- 238000003775 Density Functional Theory Methods 0.000 description 1
- 238000000862 absorption spectrum Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010494 dissociation reaction Methods 0.000 description 1
- 230000005593 dissociations Effects 0.000 description 1
- 230000005283 ground state Effects 0.000 description 1
- 150000002500 ions Chemical class 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
単一参照関数理論による非平衡構造の計算を可能とすると共に、ポテンシャル曲面計算の高速化も可能とする。
【解決手段】
対象となる系の分子軌道の初期値として、SCF計算過程が収束する構造のうち、対象となる系に比較的近い構造を持つものの分子軌道を計算し、計算された分子軌道を初期値として、対象系に構造を近づけた系について再度分子軌道を計算することを続けていき、対象系の分子軌道計算について対象系に比較的構造が近いものから計算された分子軌道を初期値として計算する。
【選択図】図1
Description
Consistent field)過程における初期分子軌道作成に関するものである。
Configuration SCF(MCSCF)法を利用し、さらに電子相関を考慮する補正を加えていくが、計算コストは密度汎関数法よりもはるかに高い。
B3LYPを利用し、基底関数には、内殻軌道を1つの関数で、最外殻軌道を2つの関数で記述するVDZ(Valence Triple Zeta)を利用した。
20%程度短くなった。
Claims (9)
- 計算機で目的とする系の状態を求める第一原理計算において、
SCF計算の収束が可能な構造で、かつ前記目的とする系に近い構造を持つ第一の系を特定し、
前記第一の系の第一原理計算により第一の系の分子軌道を算出し、
前記目的とする系の分子軌道の初期値として前記第一の系の分子軌道を用いて、前記目的とする系の分子軌道を計算することを特徴とする第一原理計算法。 - 請求項1に記載された第一原理計算法であって、
前記第一の系は、前記目的とする系の構造との変位差の二乗和平方根が1.0Å 以内であることを特徴とする第一原理計算法。 - 計算機で目的とする系の状態を求める第一原理計算において、
SCF計算の収束が可能な構造で、かつ前記目的とする系に近い構造を持つ第一の系を特定し、
前記第一の系の第一原理計算により第一の系の分子軌道を算出し、
前記第一の系よりも前記目的とする系の構造との変位差が小さく、かつ第一の系の構造との変位差の二乗和平方根が1.0Å 以内の構造を持つ第二の系を特定し、
前記第二の系の分子軌道の初期値として前記第一の系の分子軌道を用いて、前記第二の系の分子軌道を算出し、
前記第二の系の特定と同様に、第N−1の系より第Nの系を特定し、第N−1の系の分子軌道を初期値として用いて第Nの系の分子軌道を第一原理計算により計算し、
前記第Nの系の計算を複数回繰り返すことにより目的とする系の分子軌道を計算することを特徴とする第一原理計算法。 - 第一原理計算を用いた計算において、平衡構造を持つ系E、もしくは平衡構造との変位差の二乗和平方根が1.0Å 以内となる構造を持つ系E′について初期分子軌道生成過程に半経験的な手法や経験的手法を利用した従来の第一原理計算法により計算し、系Eの構造か系E′の構造との変位差の二乗和平方根が1.0Å 以内で、かつ系Eの構造との変位差が系E′の構造より大きい構造を持つ系Aの第一原理計算について、分子軌道の初期値として系Eか系E′の計算で得られた分子軌道を初期値として系Aについて計算し、系Aの計算で得た分子軌道を初期値として、系Eの構造か系E′の構造との変位差が系Aの構造より大きく、かつ系Aの構造との変位差の二乗和平方根が1.0Å 以内である構造を持つ系Bを計算し、この手順を繰り返すことで、求めるポテンシャル曲面を計算することを特徴とする第一原理計算法。
- 計算機で第一原理計算法を用いてポテンシャル曲面を算出する方法であって、
前記ポテンシャル曲面を構成する複数の構造を特定し、
前記複数の構造のそれぞれについて、第一原理計算法により所定時間内のSCF計算を実施して分子軌道を計算し、
前記複数の構造のうち、SCF計算が収束した一または複数の構造A、及びSCF計算が収束しない一または複数の構造Zを特定し、
前記構造Zのうちの各構造zについて、前記構造zとの最も変異差の小さい構造Aのうちの一の構造aを特定し、前記構造aとの変位差の二乗和平方根が1.0Å 以下である構造bを特定し、前記構造aの分子軌道を初期値として用いて前記構造bの分子軌道を計算し、前記構造bの分子軌道に基づく分子軌道計算を繰り返して構造Zの分子軌道を計算することを特徴とするポテンシャル曲面の計算方法。 - 第一原理計算を用いたポテンシャル曲面の計算において、まず初期分子軌道生成過程に半経験的な手法や経験的手法を利用した従来の第一原理計算法を用いて全ての構造について計算を実施し、SCF計算が収束しなかった系Zについては、SCF計算が収束していて、かつ系Zの構造との変位差が最も小さい構造を持つ系Aの分子軌道を初期値として、系Zの構造との変位差が系Aの構造より小さくなり、かつ系Aの構造系Bの構造を計算し、系Bの構造にて計算された分子軌道を初期値として、系Bの構造より系Zの構造との変位差が小さく、かつ系Bの構造との変位差の二乗和平方根が1.0Å 以下である構造を持つ系Cを計算し、この手順を繰り返すことで、従来法では収束しなかった系Zにおける計算を収束させ、求めるポテンシャル曲面を計算することを特徴とする第一原理計算法。
- 目的化合物のポテンシャル曲面を第一原理計算により算出する計算方法であって、
前記目的化合物と異なる第一の化合物と、前記第一の化合物よりも非平衡構造で、かつ前記第一の化合物よりも目的化合物に類似する第二の化合物を選択する第一のステップと、前記第一の化合物の分子軌道を計算する第二のステップと、前記計算された分子軌道を初期値に用いて前記第二の化合物の分子軌道を計算する第三のステップを有し、
前記第一ないし第三のステップを複数回繰り返して、前記目的化合物の分子軌道を計算することを特徴とするポテンシャル曲面の計算方法。 - 目的化合物のポテンシャル曲面を第一原理計算により算出する計算方法であって、
前記目的化合物よりも平衡状態に近い構造の基準化合物の分子軌道を計算し、
前記計算された分子軌道を初期値として前記第一原理計算に使用し、
前記目的化合物の分子軌道を計算することを特徴とするポテンシャル曲面の計算方法。 - 請求項8に記載されたポテンシャル曲面の計算方法において、
前記基準化合物は、前記目的化合物との変位差の二乗和平方根が1.0Å 以内となる構造であることを特徴とするポテンシャル曲面の計算方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007048381A JP5145731B2 (ja) | 2007-02-28 | 2007-02-28 | 化合物のシミュレーション方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007048381A JP5145731B2 (ja) | 2007-02-28 | 2007-02-28 | 化合物のシミュレーション方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2008210306A true JP2008210306A (ja) | 2008-09-11 |
JP5145731B2 JP5145731B2 (ja) | 2013-02-20 |
Family
ID=39786523
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2007048381A Expired - Fee Related JP5145731B2 (ja) | 2007-02-28 | 2007-02-28 | 化合物のシミュレーション方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5145731B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015132865A1 (ja) * | 2014-03-03 | 2015-09-11 | 富士通株式会社 | シミュレーション方法、記憶媒体、及び装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04332864A (ja) * | 1991-05-09 | 1992-11-19 | Nec Software Ltd | 分子軌道初期値生成装置 |
JP2000011025A (ja) * | 1998-06-23 | 2000-01-14 | Toshiba Corp | 材料設計支援方法 |
JP2002109446A (ja) * | 2000-09-29 | 2002-04-12 | Toshiba Corp | シミュレーション方法、シミュレーション・システムおよびシミュレーション・プログラムを記録した記録媒体 |
JP2005050321A (ja) * | 2003-07-17 | 2005-02-24 | Toray Ind Inc | プログラム及びそれを用いた計算方法 |
JP2006215804A (ja) * | 2005-02-03 | 2006-08-17 | National Institute Of Information & Communication Technology | 最適値を求めるためのコンピュータシステム |
JP2006277582A (ja) * | 2005-03-30 | 2006-10-12 | Nec Corp | 分子軌道法計算方法及び計算装置 |
-
2007
- 2007-02-28 JP JP2007048381A patent/JP5145731B2/ja not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04332864A (ja) * | 1991-05-09 | 1992-11-19 | Nec Software Ltd | 分子軌道初期値生成装置 |
JP2000011025A (ja) * | 1998-06-23 | 2000-01-14 | Toshiba Corp | 材料設計支援方法 |
JP2002109446A (ja) * | 2000-09-29 | 2002-04-12 | Toshiba Corp | シミュレーション方法、シミュレーション・システムおよびシミュレーション・プログラムを記録した記録媒体 |
JP2005050321A (ja) * | 2003-07-17 | 2005-02-24 | Toray Ind Inc | プログラム及びそれを用いた計算方法 |
JP2006215804A (ja) * | 2005-02-03 | 2006-08-17 | National Institute Of Information & Communication Technology | 最適値を求めるためのコンピュータシステム |
JP2006277582A (ja) * | 2005-03-30 | 2006-10-12 | Nec Corp | 分子軌道法計算方法及び計算装置 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015132865A1 (ja) * | 2014-03-03 | 2015-09-11 | 富士通株式会社 | シミュレーション方法、記憶媒体、及び装置 |
JPWO2015132865A1 (ja) * | 2014-03-03 | 2017-03-30 | 富士通株式会社 | シミュレーション方法、記憶媒体、及び装置 |
Also Published As
Publication number | Publication date |
---|---|
JP5145731B2 (ja) | 2013-02-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Quirynen et al. | Autogenerating microsecond solvers for nonlinear MPC: a tutorial using ACADO integrators | |
Minion | Semi-implicit projection methods for incompressible flow based on spectral deferred corrections | |
Förster et al. | Low-order scaling quasiparticle self-consistent GW for molecules | |
KR20130053374A (ko) | 자동화된 멀티레벨 부분구조화를 이용한 부분구조 생성 | |
Ohtsuka et al. | Selected configuration interaction method using sampled first-order corrections to wave functions | |
EP2680157B1 (en) | Co-simulation procedures using full derivatives of output variables | |
Akinaga et al. | Two-component relativistic equation-of-motion coupled-cluster methods for excitation energies and ionization potentials of atoms and molecules | |
Schmitz et al. | Comparison of time splitting and backsubstitution methods for integrating Vlasov's equation with magnetic fields | |
Park et al. | A mode selection criterion based on flexibility approach in component mode synthesis | |
JP5145731B2 (ja) | 化合物のシミュレーション方法 | |
Arjmandi et al. | Computing mode shapes of fluid-structure systems using subspace iteration methods | |
Iserles et al. | Applications of Magnus expansions and pseudospectra to Markov processes | |
US8510084B2 (en) | Method, apparatus and computer program for multiple time stepping simulation of a thermodynamic system using shadow hamiltonians | |
Zheng et al. | Modulus-based successive overrelaxation method for pricing American options | |
Cross et al. | Level set methods for computing reachable sets of systems with differential algebraic equation dynamics | |
Kim et al. | A mode selection algorithm for the flexibility based component mode synthesis | |
Chiew et al. | Stability analysis of dual-time stepping | |
Irisa et al. | Solvation free energies of non-polar and polar solutes reproduced by a combination of extended scaled particle theory and the Poisson-Boltzmann equation | |
Wallraff et al. | Higher order multigrid algorithms for a discontinuous Galerkin RANS solver | |
Kanal et al. | Theory and algorithm of the inversion method for pentadiagonal matrices | |
Ma et al. | MATLAB Simulink modeling and simulation of Zhang neural networks for online time-varying Sylvester equation solving | |
Tseng et al. | A novel design space reduction method for efficient simulation-based optimization | |
Anto-Sztrikacs et al. | Effective Hamiltonian theory: An approximation to the equilibrium state of open quantum systems | |
Huang et al. | Novel approach to multibody system modeling: Cascading and clustering | |
Sabet et al. | Forward kinematic analysis of non-deterministic articulated multibody systems with kinematically closed-loops in polynomial chaos expansion scheme |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20090218 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20111101 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20111228 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20120619 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20120809 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20121030 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20121112 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20151207 Year of fee payment: 3 |
|
LAPS | Cancellation because of no payment of annual fees |