この記事には参考文献や外部リンクの一覧が含まれていますが、脚注によって参照されておらず、情報源が不明瞭です。脚注を導入して、記事の信頼性向上にご協力ください。(2015年11月)
能勢=フーバー・サーモスタット(のせ=フーバー・サーモスタット、英: Nose?Hoover thermostat)は、等温分子動力学シミュレーションのための決定的アルゴリズムである。初めに能勢修一によって開発され[1]、フーバー
(英語版)によってさらに改良された[2]。能勢=フーバー・サーモスタットの熱浴はただ一つの想像上の粒子から成るが、シミュレーション系は現実的な定温条件(正準集団)を満たす。そのため、能勢=フーバー・サーモスタットは定温分子動力学シミュレーションのための最も正確かつ効率的な手法の一つとして一般的に使用されている。古典的分子動力学において、シミュレーションはもっとも単純には小正準集団(粒子の数、容積、エネルギー (N, V, E) が一定値を有する)中で行われる。しかしながら、一般的に実際の実験条件では、エネルギーの代わりに温度が制御される(N, V, T 一定)。この実験条件を記述するアンサンブル(統計集団)は正準集団と呼ばれる。重要なことに、統計力学の観点からすれば、正準集団は小正準集団とは全く異なる。そこで、小正準集団を用いつつも温度を一定に保つための複数の手法が発表されている。温度を制御するための人気のある手法としては、速度リスケーリング(V-rescale)、アンダーセン・サーモスタット、能勢=フーバー・サーモスタット、能勢=フーバー・チェイン法、ベレンゼン・サーモスタット、ランジュバン動力学がある。
中心となる考えは正準分布を得るようなやり方でシミュレーションを行うことである。これはシミュレーション下の系の平均温度を固定することを意味するが、同時に正準分布に典型的な分布に従った温度の変動を許す。 能勢のアプローチでは、追加の自由度 s を持つ熱浴
能勢=フーバー・サーモスタット
H ( P , R , p s , s ) = ∑ i p i 2 2 m i s 2 + U ( R ) + p s 2 2 Q + g k T ln ( s ) {\displaystyle {\mathcal {H}}(P,R,p_{s},s)=\sum _{i}{\frac {{\boldsymbol {p}}_{i}^{2}}{2m_{i}s^{2}}}+U(R)+{\frac {p_{s}^{2}}{2Q}}+gkT\ln \left(s\right)}
上式において、g は系の独立した運動量自由度の数であり、R および P は全ての座標 {ri} および {pi} を表わし、Q は系にあわせて慎重に選ぶべき仮想質量である。このハミルトニアンにあらわれる座標 R, P, t は仮想のものであり、実座標とは以下のような関係式で結び付いている。
R ′ = R , P ′ = P s and t ′ = ∫ t d τ s {\displaystyle R'=R,~P'={\frac {P}{s}}~{\text{and}}~t'=\int ^{t}{\frac {\mathrm {d} \tau }{s}}}
上式において、プライムが付いた座標が実座標である。注目すべきは、g = 3N における上記のハミルトニアンの集団平均が小正準集団平均と等しいことである。
フーバーは能勢の方法に位相空間における連続条件、一般化リウヴィル方程式を導入して改良し、現在能勢=フーバー・サーモスタットと呼ばれている手法を確立した[2]。この手法は s による時間スケーリングを必要としない。
脚注^ Nose 1984.
^ a b Hoover 1985.
^ Nose 1984, 式2.5.
参考文献
Nose, Shuichi (1984). “A unified formulation of the constant temperature molecular-dynamics methods”. Journal of Chemical Physics 81 (1): 511?519. Bibcode: 1984JChPh..81..511N. doi:10.1063/1.447334.
Hoover, William G. (Mar 1985). “Canonical dynamics: Equilibrium phase-space distributions”. Phys. Rev. A (American Physical Society) 31 (3): 1695?1697. Bibcode: 1985PhRvA..31.1695H. doi:10.1103/PhysRevA.31.1695.
Thijssen, J. M. (2007). Computational Physics (2nd ed.). Cambridge University Press. pp. 226?231. .mw-parser-output cite.citation{font-style:inherit;word-wrap:break-word}.mw-parser-output .citation q{quotes:"\"""\"""'""'"}.mw-parser-output .citation.cs-ja1 q,.mw-parser-output .citation.cs-ja2 q{quotes:"「""」""『""』"}.mw-parser-output .citation:target{background-color:rgba(0,127,255,0.133)}.mw-parser-output .id-lock-free a,.mw-parser-output .citation .cs1-lock-free a{background:url("//upload.wikimedia.org/wikipedia/commons/6/65/Lock-green.svg")right 0.1em center/9px no-repeat}.mw-parser-output .id-lock-limited a,.mw-parser-output .id-lock-registration a,.mw-parser-output .citation .cs1-lock-limited a,.mw-parser-output .citation .cs1-lock-registration a{background:url("//upload.wikimedia.org/wikipedia/commons/d/d6/Lock-gray-alt-2.svg")right 0.1em center/9px no-repeat}.mw-parser-output .id-lock-subscription a,.mw-parser-output .citation .cs1-lock-subscription a{background:url("//upload.wikimedia.org/wikipedia/commons/a/aa/Lock-red-alt-2.svg")right 0.1em center/9px no-repeat}.mw-parser-output .cs1-ws-icon a{background:url("//upload.wikimedia.org/wikipedia/commons/4/4c/Wikisource-logo.svg")right 0.1em center/12px no-repeat}.mw-parser-output .cs1-code{color:inherit;background:inherit;border:none;padding:inherit}.mw-parser-output .cs1-hidden-error{display:none;color:#d33}.mw-parser-output .cs1-visible-error{color:#d33}.mw-parser-output .cs1-maint{display:none;color:#3a3;margin-left:0.3em}.mw-parser-output .cs1-format{font-size:95%}.mw-parser-output .cs1-kern-left{padding-left:0.2em}.mw-parser-output .cs1-kern-right{padding-right:0.2em}.mw-parser-output .citation .mw-selflink{font-weight:inherit}ISBN 978-0-521-83346-2
外部リンク
⇒Berendsen and Nose-Hoover thermostats
⇒A simple (c++) implementation of the Nose-Hoover chains thermostat