界面の解析 -NPzT計算-

LAMMPSにはfix npt コマンドなど、もとからNPT計算用のコマンドは用意されています。しかし、多くのコマンドではxyz方向全てで周期境界条件であることが前提になります。
例えば、気液界面の計算をしたい時などはZ方向に真空領域を用意しますが、このとき周期境界条件で計算された解析結果を使うのは難しいです。

周期境界条件を使用しないための計算の手順を書いていきます
①液体を計算空間に配置する
②原子IDを一つ増やし、相互作用の値を設定する
③増やしたIDの原子を液体の上下にxyの2次元平面となるよう配置

									def hello_world():
										print("Hello, World!")
											
										hello_world()
									

Interdum et malesuada fames ac ante ipsum primis in faucibus. Pellentesque venenatis dolor imperdiet dolor mattis sagittis. Praesent rutrum sem diam, vitae egestas enim auctor sit amet. Pellentesque leo mauris, consectetur id ipsum sit amet, fersapien risus, commodo eget turpis at, elementum convallis elit. Pellentesque enim turpis, hendrerit tristique lorem ipsum dolor.

Phone

(000) 000-0000 x12387

Address

1234 Somewhere Road #5432
Nashville, TN 00000
United States of America