第一原理分子動力学プログラム STATE Senri Wiki

• 追加された行はこの色です。
• 削除された行はこの色です。
• STATE_tutorial/drho へ行く。
• STATE_tutorial/drho の差分を削除

#author("2024-02-14T10:45:51+09:00","default:StatE","StatE")
#author("2024-02-14T10:48:37+09:00","default:StatE","StatE")
* drho* [#lb12db6b]
差電荷密度とそのxy面平均を計算を実行するためのプログラム群で、以下にあります。
 state-5.6.14/util/chargeutil
ここでSTATEのバージョンは5.6.14を想定していますが、他のバージョンでも同様に動作します。~
主に使うのは以下のプログラムです。
- drho1~
$\Delta \rho = \rho_{\mathrm{A}}(\mathbf{r}) - \rho_{\mathrm{B}}(\mathbf{r})$
- drho2~
$\Delta \rho = \rho_{\mathrm{A}}(\mathbf{r}) - \rho_{\mathrm{B}}(\mathbf{r}) -  \rho_{\mathrm{C}}(\mathbf{r})$~
特にdrho2は吸着系における以下のような差電荷密度の計算をすることを目的に書かれています。~
$\Delta \rho = \rho_{\mathrm{mol/sub}}(\mathbf{r}) - \rho_{\mathrm{mol}}(\mathbf{r}) -  \rho_{\mathrm{sub}}(\mathbf{r})$~
ここで$\Delta \rho = \rho_{\mathrm{mol/sub}}(\mathbf{r})$、$\rho_{\mathrm{mol}}(\mathbf{r})$、$  \rho_{\mathrm{sub}}(\mathbf{r})$はそれぞれ吸着系の電荷密度、吸着分子の電荷密度、基板の電荷密度です。

**使い方 [#sdc04bfd]
ここではrho2を使用する場合を考えます。
先ず
 $ drho2
を実行すると