JPH10334076A - 分子動力学法における温度制御方法 - Google Patents

分子動力学法における温度制御方法

Info

Publication number
JPH10334076A
JPH10334076A JP14038897A JP14038897A JPH10334076A JP H10334076 A JPH10334076 A JP H10334076A JP 14038897 A JP14038897 A JP 14038897A JP 14038897 A JP14038897 A JP 14038897A JP H10334076 A JPH10334076 A JP H10334076A
Authority
JP
Japan
Prior art keywords
temperature
molecular dynamics
control method
state
temperature control
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.)
Pending
Application number
JP14038897A
Other languages
English (en)
Inventor
Haruo Kato
治男 加藤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NEC Corp
Original Assignee
NEC Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by NEC Corp filed Critical NEC Corp
Priority to JP14038897A priority Critical patent/JPH10334076A/ja
Publication of JPH10334076A publication Critical patent/JPH10334076A/ja
Pending legal-status Critical Current

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Testing Or Measuring Of Semiconductors Or The Like (AREA)
  • Complex Calculations (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

(57)【要約】 【課題】分子動力学計算において、温度設定を変更した
場合に、解析状態が迅速に安定化し、所望の設定温度T
oでの状態を正確に得ることができるようにする。 【解決手段】系の温度をTとし、温度Tに対応する各粒
子の速度をvとし、定数γを0<γ<1とし、 【数1】 にしたがって速度vを更新し、更新された速度vに応じ
て温度Tを更新し(ステップ102)、予め定められた
しきい値εに対して|To−T|<εを満たすまで、速
度vと温度Tの更新を繰り返し(ステップ103)、そ
の後、設定温度T oでの分子動力学計算を実行する(ス
テップ104)。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、各種のシミュレー
ションのための計算物理学的手法である分子動力学法に
関し、特に、分子動力学法によるシミュレーション計算
において計算モデルの温度を制御する方法に関する。
【0002】
【従来の技術】各種のシミュレーションを行うための手
法の一つとして、分子動力学(MD;Molecular Dynami
cs)法がある。分子動力学法は、液相や凝集系などの系
を構成する多数の分子、原子の運動をコンピュータによ
る数値計算でシミュレーションすることによって、その
系の統計力学的な量などを求める技術である。例えば、
半導体装置製造の分野において、金属配線におけるエレ
クトロマイグレーション(EM;Electromigration)現
象のシミュレーションや解析などに、この分子動力学法
は用いられている。
【0003】さて、分子動力学法によるシミュレーショ
ンを実行する場合、解析対象の系である計算モデルの設
定温度をある温度から別の温度に変化させたいことがあ
る。ところで、よく知られているように、粒子の運動エ
ネルギーと系の温度とは線形関係にある。つまり、粒子
の運動速度の2乗が温度と線形の関係を有する。分子動
力学法では、解析対象の温度を初期温度Tから設定温度
oに変更する場合、従来、例えば、河村雄行著「パソ
コン分子シミュレーション」第44頁(海文堂)にある
ように、粒子の初期速度vに(To/T)1/2を乗じて粒
子の速度を変更し、これによって運動エネルギーを変化
させることにより、温度を変更している。
【0004】図4は、この従来の温度制御方法の手順を
示すフローチャートである。系の初期温度T、粒子の初
期速度vとして分子動力学計算を行っていたとして(ス
テップ111)、
【0005】
【数2】 にしたがって速度を更新して設定速度voとし、温度を
更新して設定温度Toとする(ステップ112)。そし
て、このような設定速度vo及び設定温度Toによって、
分子動力学計算を再開する(ステップ113)。
【0006】
【発明が解決しようとする課題】上述した従来の温度制
御方法により各粒子の速度を変更し、分子動力学による
計算シミュレーションでの粒子間の運動を再開させた場
合、設定温度と初期温度との差が大きいときなどに、各
粒子間に働く粒子間ポテンシャルによる相互作用によっ
て受ける力のために解析状態が安定せず、所望の温度
(設定温度)での状態に迅速に収束することができない
という問題点がある。
【0007】本発明の目的は、温度設定を変更した場合
に、解析状態が迅速に安定化し、所望の設定温度での状
態を正確に得ることができる温度制御方法を提供するこ
とにある。
【0008】
【課題を解決するための手段】本発明の第1の温度制御
方法は、分子動力学法による計算を実行する際に、初期
状態から指定された設定温度Toに解析対象の系の状態
を変化させる温度制御方法において、初期状態の温度と
設定温度Toとの間に中間の温度段階を設定し、中間の
温度段階を経由して系の状態を変化させる。
【0009】本発明の第2の温度制御方法は、分子動力
学法による計算を実行する際に、初期状態から解析に必
要な指定された設定温度Toに系の状態を変化させる温
度制御方法において、系の温度をTとし、温度Tに対応
する各粒子の分子動力学法での速度をvとし、γを0<
γ<1を満たす予め定められた値とし、
【0010】
【数3】 にしたがって速度vを更新し、更新された速度vに応じ
て温度Tを更新し、予め定められたしきい値εに対して
|To−T|<εを満たすまで、速度vと温度Tの更新
を繰り返し、その後、設定温度Toでの分子動力学計算
を実行する。この温度制御方法において、初期状態の温
度Tが絶対零度か絶対零度に近い場合には、上述の式が
発散するので、予め定めた仮の温度T'(T'>0)を温
度Tに代入して速度vの更新を行うようにすることが好
ましい。
【0011】本発明では、解析対象の系の温度を、初期
状態の温度から一気に設定温度に変化させるのではな
く、これら温度の中間の温度段階を経て、これら中間の
温度段階ごとに分子動力学計算を行って解析対象の系を
安定化させるようにする。このように構成することによ
り、1回ごとの分子動力学計算での温度変化が小さいた
め、急激な運動エネルギーの変化が抑えられ、そのた
め、設定温度で系が安定化するまでのトータルの時間
は、一気に設定温度まで温度を変化させる場合に比べ、
大幅に短縮される。
【0012】
【発明の実施の形態】次に、本発明の実施の形態につい
て、図面を参照して説明する。図1は本発明の実施の一
形態の分子動力学法における温度制御方法の手順を示す
フローチャートである。ここでは、粒子の初期速度をv
とし、この初期速度に対応する初期温度をTとして、解
析対象の系の温度を設定温度Toに設定するものとす
る。設定温度Toに対応する速度を設定速度voとする。
【0013】初期温度T、初期速度vでの分子動力学計
算を行っていたとする(ステップ101)。ここで、0
<γ<1を満たす定数γを考える。例えば、γ=0.5
とする。そして、
【0014】
【数4】 にしたがって各粒子の速度vを更新し、更新された速度
vに対応する温度Tを求め、分子動力学計算を実行する
(ステップ102)。更新後の温度Tと設定温度To
の差を求め、この差の絶対値が予め定めたしきい値εよ
り小さいかどうかを判別する(ステップ103)。ここ
でしきい値εは、設定温度Toに応じて変化させるよう
にしてもよい。|To−T|≧εの場合には、温度Tが
設定温度Toにまだ十分に接近していないので、温度T
をさらに設定温度Toに近づけるために、ステップ10
2に戻る。|To−T|<εであれば、ステップ102
での更新後の温度Tが十分に設定温度Toに接近してい
る場合であるので、設定温度To、設定速度voでの分子
動力学計算を行い(ステップ104)、所望の結果を得
て、処理を終了する。
【0015】本実施の形態では、初期温度から設定温度
へ一気に系の温度を変化させるのではなく、式(2)にし
たがって徐々に温度を変化させ、温度変化の各段階でそ
れぞれ分子動力学計算を実行する。温度変化の各段階で
は、初期温度から設定温度へ一気に系の温度を変化させ
る場合に比べ、前後の段階との運動エネルギーの差が小
さいので、解析系の変動が小さく、分子動力学計算にお
いて解析対象の系の安定化に要する時間が短い。温度変
化の各段階ごとに分子動力学計算を行うとしても、それ
ぞれの段階での計算時間が短くてすむので、本実施の形
態によれば、初期温度から設定温度へ一気に系の温度を
変化させる場合に比べ、トータルとしての計算時間を短
くすることができ、設定温度での安定状態により早く到
達することができる。
【0016】
【実施例】以下、具体的な実施例によって本発明をさら
に詳しく説明する。
【0017】《実施例1》図2は、実施例1における解
析領域を説明する図である。ここでは、初期状態では解
析領域1内の各粒子2が静止状態、すなわち温度Tが0
K(絶対零度)であるとして、この状態からある設定温
度Toに系の温度を変化させるものとする。この場合、
図1に示すフローチャートにしたがい、速度vを少しず
つ変更することを繰り返して所望の設定温度Toに至ら
しめるのであるが、T=0を代入すると、式(2)の平方
根中の分母が0となり、この式が発散する。この困難を
避けるため、仮に温度Tを1Kぐらいの小さい正の値
T'として、式(2)に代入するようにすればよい。
【0018】《実施例2》図3は、実施例2における解
析領域を説明する図である。ここでは、初期状態とし
て、解析領域1内の各粒子2がそれぞれ速度3を持つ状
態、すなわち、その速度に対応する温度T(T>0)の
状態を考え、この状態から別の設定温度Toに系の温度
を設定するものとする。この場合も、図1に示すフロー
チャートにしたがい、速度vを少しずつ変更することを
繰り返して所望の設定温度Toに至らしめる。
【0019】
【発明の効果】以上説明したように本発明は、解析対象
の系の温度を、初期状態の温度から一気に設定温度に変
化させるのではなく、これら温度の中間の温度段階(ス
テップ)を経るとともにこれら中間の温度段階ごとに分
子動力学計算を行うことにより、解析系の安定化に要す
る時間が各ステップで短く、かつ一気に温度を上下させ
る場合に比べて各ステップでの運動エネルギーの変化量
が少なくて解析系の変動が小さいので、結局、より短時
間で系が安定化するようになるという効果がある。
【図面の簡単な説明】
【図1】本発明の実施の一形態の分子動力学法における
温度制御方法の手順を示すフローチャートである。
【図2】実施例1における解析領域を説明する図であ
る。
【図3】実施例2における解析領域を説明する図であ
る。
【図4】従来の分子動力学法における温度制御方法の手
順を示すフローチャートである。
【符号の説明】
1 解析領域 2 粒子 3 粒子の速度 101〜104,111〜113 ステップ

Claims (3)

    【特許請求の範囲】
  1. 【請求項1】 分子動力学法による計算を実行する際
    に、初期状態から指定された設定温度Toに解析対象の
    系の状態を変化させる温度制御方法において、 初期状態の温度と設定温度Toとの間に中間の温度段階
    を設定し、中間の温度段階を経由して系の状態を変化さ
    せることを特徴とする温度制御方法。
  2. 【請求項2】 分子動力学法による計算を実行する際
    に、初期状態から指定された設定温度Toに系の状態を
    変化させる温度制御方法において、 系の温度をTとし、温度Tに対応する各粒子の分子動力
    学法での速度をvとし、γを0<γ<1を満たす予め定
    められた値とし、 【数1】 にしたがって速度vを更新し、更新された速度vに応じ
    て温度Tを更新し、 予め定められたしきい値εに対して|To−T|<εを
    満たすまで、速度vと温度Tの更新を繰り返し、 その後、設定温度Toでの分子動力学計算を実行するこ
    とを特徴とする温度制御方法。
  3. 【請求項3】 温度Tが絶対零度か絶対零度に近い場合
    には、予め定めた仮の温度T'(T'>0)を温度Tに代
    入して速度vの更新を行う請求項2に記載の温度制御方
    法。
JP14038897A 1997-05-29 1997-05-29 分子動力学法における温度制御方法 Pending JPH10334076A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP14038897A JPH10334076A (ja) 1997-05-29 1997-05-29 分子動力学法における温度制御方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP14038897A JPH10334076A (ja) 1997-05-29 1997-05-29 分子動力学法における温度制御方法

Publications (1)

Publication Number Publication Date
JPH10334076A true JPH10334076A (ja) 1998-12-18

Family

ID=15267658

Family Applications (1)

Application Number Title Priority Date Filing Date
JP14038897A Pending JPH10334076A (ja) 1997-05-29 1997-05-29 分子動力学法における温度制御方法

Country Status (1)

Country Link
JP (1) JPH10334076A (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3467836A2 (en) 2017-10-06 2019-04-10 Sumitomo Heavy Industries, Ltd. Particle simulation method, simulation apparatus, and program
RU2727631C1 (ru) * 2018-08-23 2020-07-22 Тата Консалтенси Сервисес Лимитед Системы и способы прогнозирования структуры и свойств атомарных веществ и их сплавов

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3467836A2 (en) 2017-10-06 2019-04-10 Sumitomo Heavy Industries, Ltd. Particle simulation method, simulation apparatus, and program
EP3467836A3 (en) * 2017-10-06 2019-06-05 Sumitomo Heavy Industries, Ltd. Particle simulation method, simulation apparatus, and program
US10691147B2 (en) 2017-10-06 2020-06-23 Sumitomo Heavy Industries, Ltd. Simulation method, simulation apparatus, and program
RU2727631C1 (ru) * 2018-08-23 2020-07-22 Тата Консалтенси Сервисес Лимитед Системы и способы прогнозирования структуры и свойств атомарных веществ и их сплавов

Similar Documents

Publication Publication Date Title
Tan et al. Stable proportional-derivative controllers
JP4223894B2 (ja) Pidパラメータ調整装置
Nguyen-Tuong et al. Incremental online sparsification for model learning in real-time robot control
JP2019537132A (ja) アクション選択ニューラルネットワークをトレーニングすること
Han et al. Cellular automata modeling of static recrystallization based on the curvature driven subgrain growth mechanism
JP7215077B2 (ja) 予測プログラム、予測方法及び予測装置
Marshall et al. Autonomous Development of a Machine-Learning Model for the Plastic Response of Two-Phase Composites from Micromechanical Finite Element Models: Marshall and Kalidindi
US10691147B2 (en) Simulation method, simulation apparatus, and program
CN110554605A (zh) 一种基于约束跟踪的复杂机械系统自适应鲁棒控制方法
JPH0793290A (ja) カオスフィードバックシステム
KR20230157426A (ko) 노이즈 제거 목적을 사용하는 그래프 신경 네트워크 트레이닝
Ulz Coupling the finite element method and molecular dynamics in the framework of the heterogeneous multiscale method for quasi-static isothermal problems
Kunz et al. Atomic clusters and nanoscale particles: from coarse-grained dynamics to optimized annealing schedules
Denniston et al. LAMMPS lb/fluid fix version 2: Improved hydrodynamic forces implemented into LAMMPS through a lattice-Boltzmann fluid
Xu et al. Meta-learning via weighted gradient update
CN116595888A (zh) 基于傅里叶神经网络的生物动力学模型求解方法及系统
JPH10334076A (ja) 分子動力学法における温度制御方法
Nguyen–Tuong et al. Incremental sparsification for real-time online model learning
CN114839874B (zh) 一种面向系统模型部分未知的平行控制方法及系统
Chao et al. Dictionary-based fidelity measure for virtual traffic
Wolfrum et al. Complexity reduction of a thin film deposition model using a trajectory based nonlinear model reduction technique
Williamson et al. Granular synthesis for display of time-varying probability densities
Yoshioka et al. Fast Gaussian process regression using representative data
Mühlbacher et al. Electron transfer rates for asymmetric reactions
Van der Velden High-fidelity simulation surrogate models for systems engineering