SLIDE 24 Grand canonical ensemble in electronic system
ｙ ｀ ６
４ ０ ０ ３ ５ ０ ３ ０ ０ ０ ２ ４ ６ ８ １ ０ １ ／ ｐ ｓ 図 ７ ． ３ 水 の 温 度 制 御 １＝５ｐｓにおいて，温度７；，を３００Ｋから３５０Ｋにあげた． Ｅ Ｓ に ｀ ６ ０ ０ ０ ４ ０ ０ ０ ２ ０ ０ ０ ０
２ ０ ０ ０
７ 。 ５ バ リ ネ ロ ・ ラ ー マ ン の 方 法 月 ア ０ ２ ４ ６ ８ １ ０ １ ／ ｐ ｓ 図７．４水のＪ１モカ制御 ＝５ｐｓにおいて，圧力／４を１０００ａｔｍから３０００ａｕ１１にあげた．
こ こ で は ， 立 方 休 セ ル に 対 し 休 秘 レ を 一 ・ 般 化 座 標 と し て 運 助 方 程 式 を 導 い た ． も ち ろ ん ， 体 秘 ド の 代 わ り に セ ル の 一 辺 の 長 さ 乙 ＝ 副 （ ７ ． ５ ３ ） を 変 数 と し て も か ま わ な い ． さ ら に は ， 直 方 体 セ ル に 対 し ， 各 辺 の 長 さ
回 口 し
を 一 般 化 座 標 と し て ， 三 つ の 方 向 の 圧 力 が そ れ ぞ れ
ｊ に 以
（ ７ ． ５ ４ ） （ ７ ． ５ ５ ）
となるように制御することもできる．たとえば，鳥・＝几＝ｈ１１１１１として，八＝ ０．５ａｔｍのようにとれば，表而張力-定のＭＤ計拿１：力４１呼能となる．次節ではこ れ を さ ら に 拡 張 し て ， セ ル の 傾 き ま で 自 山 度 を 与 え る こ と に よ り 圧 カ テ ン ソ ル Ｐ を 制 御 す る 方 法 を 示 す． ７ ． ５ パ リ ネ ロ ・ ラ ー マ ン の 方 法 ここでは，基本セルの各辺の長さが自山に変わると｜・り時に，セルの傾きも変 化 し ， 圧 力 の 平 均 値 が ハ と 一 定 と な る よ う な Ｍ Ｄ 計 算 に つ い て 説 明 す る ． こ れ に よ り ， 紘 品 系 が 界 な る ｊ 川 休 間 で の 相 転 移 や， 不 均 ・ 系 で の ひ ず み の な い ＭＤ計算が可能となる． 基 本 セ ル を 記 述 す る の に ， 図 ７ ． ５ に 示 すよ う に 平 行 六 而 体 セ ル の 各 辺 を 表 す ペ ク トル ａ ， ６ ， ｃ を 川 いて 行 列 １ を 定 義 す る ． £＝（ａろｃ）
Conventional NPT MD simulation If we can introduce a fictitious motion for amount of charge , we can realize NVTμe MD simulation μ (eV) t (ps)
ne
vcell ˙ vcell = Pvcell Mvcell ˙ Pvcell = P − Pext
: Cell volume
Mvcell : Fictitious mass for variable cell ˙ ncell = Pncell Mncell ˙ Pncell = µ − µext
from Virial theorem
23
