6. 原子構造モデリング

サイドメニューの「Modeling」をクリックすると、モデリングのタイプを選択するメニューが現れます。 現在

  • Basic(セル変形、スーパーセル、原子削除、置換、移動)

  • Surface &Slabs(表面スラブモデル)

  • Interface(界面モデル)

  • Add Molecule(表面分子吸着)

  • Add Cell(孤立分子からスーパーセルモデルを作る)

というタイプのモデリングが可能です。以下で各モデリング機能の説明を行います。

6.1. 基本モデリング

../_images/%E3%82%B9%E3%82%AF%E3%83%AA%E3%83%BC%E3%83%B3%E3%82%B7%E3%83%A7%E3%83%83%E3%83%88%202025-03-25%20082300.png

モデリングの基本画面を示します。こちらはBasicモデリングに限らず、 他のモデリングでも基本的な操作は同様の画面で行います。表示されている原子モデルは、 標準的なマウス操作により、回転、拡大・縮小、平行移動が可能です。右側にあるツール群を使って、 原子の追加、削除、置換、スーパーセルの作成、表面スラブモデルの作成等を行い、 最後に左上にある「SAVE」ボタンを押すと、作成した原子モデルが一つのジョブとして登録されます。 以下で、モデリング画面にあるツール群やボタンの説明を行います。

原子をクリック(左クリック)すると、 クリックした原子の色が変わりマーキングされます。2 個以上の原子をマーキングすると、 マーキングされた原子間の距離が表示されます。原子上で右クリックをした場合、選択した原子を削除するか、 または別の元素に置換するかの選択メニューが現れます。 「削除 (Remove)」を選択した場合はその原子がモデルから消え、 「置換 (Substitute)」を選択した場合は周期律表が現れ、その中から選んだ元素に置換されます。 さらに、初期スピン差、構造最適化時の原子移動の拘束条件 (動く方向を制限したり、完全に動かさないようにしたり)といった設定も可能です。 また表面スラブモデル作成(後述)で付加できる終端用の擬水素原子については、上記に加えて、 結合長と電荷の設定も可能になります。「移動 (Move)」を選択すると、選択した原子の上に矢印が表示され、 マウスドラッグで、その矢印の方向に原子を動かすことができます。

画面右側にあるモデリングのツール群について説明します。ツールのメニューには3つのタブがあり、 それぞれ

  • BASIC: 基本的な操作

  • CELL: Primitive cell ⇔ Conventional Cellの変換を行う

  • VIEW: 原子球のボンドの太さなど、表示に関する設定

といった機能がまとめられています。以下は全てのタブについての機能の説明となっています。

UNIT CELL

単位胞の長さを変更できます。変更したい数値を入力した後「UPDATE」 ボタンを押すと新しい格子定数のモデルになります。モデリング開始時点の単位胞の長さの変更であり、 スーパーセルモデルの単位胞の長さではないことに注意してください。

SUPERCELL

モデリング開始時点の単位胞を基準にして、a1,a2,a3 方向の単位胞の数がそれぞれ n1,n2,n3 倍となるようなスーパーセルモデルを作成します。n1,n2,n3 を指定した後「UPDATE」 ボタンを押すとスーパーセルが生成されます。

ATOMS

  • sphere:原子の球のサイズ変更

  • label:元素名+通し番号の表示/非表示切り替え

  • distance:マーキングされた原子間の距離の表示/非表示切り替え

BONDS

  • thickness:ボンドの太さの変更

  • AXIS

  • show axis:座標軸の表示/非表示切り替え

RESET

  • VIEW:回転、拡大・縮小、平行移動されたモデルの見た目を初期値に戻します

  • PLANE:原子の移動および追加のためのガイド平面を初期位置に戻します

  • MARKS:原子のマーキングを全て解除します

6.1.1. ボタン

画面左上に「SAVE」「Undo」「Redo」の4つのボタンがあります。 SAVE は作成したモデルを最後に保存するボタンです。Undo および Redo はモデリング途中の操作を戻したり、 戻すのをやめたりするボタンです。例えば、原子を 1 個追加したとき、Undo すると追加前の状態に戻ります。 追加前の状態に戻った後 Redo を押すと、再び原子が1個追加された状態になります。Undo されていない状態で Redo を押した場合は何も起こりません。

6.2. 表面スラブモデル

このモデリングモードではツール群に「SURFACE」のタブが加わります。 これに切り替えると、表面スラブモデル作成のツールが現れます。制御可能なパラメータは以下の通りです。

  • ミラー指数

  • Slab size(元になるバルク構造のユニットセルを単位とする)

  • Vacuum size(単位はスラブサイズと同じ)

  • Orthogonal_c_slab(表面に垂直な軸を c 軸とするセルを取る)

  • Cener_slab(スラブ領域をセルの中央に配置)

  • Pseudo H(表面の上部または下部を水素終端する。Z=1 以外の電荷を設定したり、表面原子との結合距離を設定することが可能

  • Flip up-side down(スラブの上下を反転させる)

その他にもいくつかのパラメータが設定可能ですが、それらはすべて pymatgen(https://pymatgen.org)の SlabGenerator の機能に準じていますので詳しくはpymatgenのドキュメントをご参照ください。 各種設定を行なったのち「GENERATE SLAB」ボタンを押すと、 ミラー指数だけで一意に決まらない面の選択ができる他、指定したミラー指数に対応する、 スラブモデルの構成単位となるバルク構造のいずれか一つが選択できるようになっています。 選択後、「select」ボタンを押すと、通常のモデリングに戻ります。 終端用の擬水素原子の結合長や電荷の設定は、この通常モデリングのモードで行うことができます。

6.3. 界面モデリング

../_images/%E3%82%B9%E3%82%AF%E3%83%AA%E3%83%BC%E3%83%B3%E3%82%B7%E3%83%A7%E3%83%83%E3%83%88%202025-03-25%20084025.png

モデリングタイプ「Interface」を選択すると上図のような設定メニューが現れます。 ジョブ詳細画面で見ている物質を基板として、その上でフィルムとなる物質の結晶構造を、 外部ファイルよりアップロードするという作成法になります。 格子定数なるべくマッチさせるように界面構造の候補を探すため、 探索範囲を面積で指定する必要もあります。設定メニューでは

  • フィルムを構成する結晶構造ファイルのアップロードまたは既存のJOBあるいはMATERIALの選択

  • フィルムおよび基板のミラー指数、および膜厚の指定

  • 界面間の初期ギャップ

  • 真空領域のサイズ

  • 界面構造の候補を探索する面積の範囲

といった設定が可能です。設定後「GENERATE」ボタンを押すと界面構造の候補が一覧表示されます (サイズの小さいモデルから順に表示されます)。 この中からターゲットとする構造を選択し(複数選択可)「REGISTER」を押すと、 選択した界面構造が一つのジョブとして登録されます。

6.4. 分子挿入

モデリングタイプ「Add Molecule」 を選択するとバルクやスラブといった周期構造に孤立分子を挿入したモデルを作成することが可能です。 「Add Molecule」を選択すると図 48のようなメニューが現れ、 ここでは孤立分子の構造ファイルをアップロードする、あるいはQuloudに登録されている孤立分子の構造から選ぶ、 という方法で作成が可能です。分子構造を選択し「ADD」を押すとモデリング画面に遷移します。 図 49ではSiスラブモデルに窒素分子を追加するモデリングの例を示します。 追加した分子に矢印が表示されており、この矢印に沿って分子をモデル内で平行移動させることができます。 また画面左上の「Rotate Molecule」のボタンを押すと、今度は分子を回転させるモードになります。 「Translate Molecule」を押すと再び平行移動モードに戻ります。

6.5. 孤立分子のスーパーセル

モデリングタイプ「Add Cell」を選択するとモデリング画面が開きます。 ここでは、例えばxyz形式のような、ユニットセルの情報を持たない孤立分子の構造に、 ユニットセルの情報を追加してスーパーセルモデルを作成することができます。スーパーセルモデルにすると、 周期系専用の第一原理計算ソフト等を用いて、分子に対する計算が可能になります。