SST Lab Dokuwiki Header header picture

@surface

トレース:

seminar:graphene

差分

このページの2つのバージョン間の差分を表示します。

この比較画面へのリンク

--- seminar:graphene [2019/12/03 15:33] – [原子の情報の設定] kimi
+++ seminar:graphene [2022/08/23 13:34] (現在) – 外部編集 127.0.0.1
@@ 行 1: / 行 1: @@
 ====== 総合演習 ======
-グラフェンを例に、より詳しい計算の演習を行おう。
+グラフェンを例に、より詳しい計算の演習を行う。
-===== 単位格子の設定 =====
+===== 計算の前に =====
+{{ :seminar:graphene.png.001.png?400 |}}
+===== プログラムの概要 =====
 <code python>
 a0 = 2.461 # Lattice constant of Graphene
@@ 行 11: / 行 15: @@
 c3 = (0.0, 0.0, a0*8.0)              # Height of 3D unit cell
-unit_cell = [c1, c2, c3] # Unit cell
-</code>
 from ase import Atom
 atom1 = Atom('C', (1.0/3.0, 2.0/3.0, 0.5))
 atom2 = Atom('C', (2.0/3.0, 1.0/3.0, 0.5))
-list_of_atoms = [atom1, atom2]
+from ase import Atoms
+box1 = Atoms([atom1, atom2], pbc = True)
+box1.set_cell([c1, c2, c3], scale_atoms = True)
+from ase.calculators.jacapo import Jacapo
+solver1 = Jacapo(nbands = 8, kpts = (2, 2, 1), pw = 150, dw = 150)
+box1.set_calculator(solver1)
+solver1.calculate()
+</code>
+===== 単位格子の設定 =====
+<code python>
+b0 = 1.42 # C-C bond length
+from math import sqrt
+a0 = b0*sqrt(3.0)
+ax = a0*sqrt(3.0)/2.0    # ax = b0*3.0/2.0
+ay = a0/2.0
+az = a0*8                # Height of 3D unit cell
+c1 = ( ax, ay, 0.0)
+c2 = (-ax, ay, 0.0)
+c3 = (0.0, 0.0, az)
+unit_cell = [c1, c2, c3] # Unit cell
+</code>
+===== 原子の位置と種類 =====
+<code python>
+v1 = (1.0/3.0, 2.0/3.0, 0.5)
+v2 = (2.0/3.0, 1.0/3.0, 0.5)
+positions = [v1, v2]
+from ase import Atom
+list_of_atoms = []
+for pos in positions:
+  atm = Atom('C', pos)
+  list_of_atoms.append(atm)
+# # list_of_atoms is a list of [atom1, atom2]
 </code>

seminar/graphene.1575354838.txt.gz · 最終更新: 2022/08/23 13:34 (外部編集)