SST Lab Dokuwiki Header
内容へ移動
@surface
ユーザ用ツール
ログイン
サイト用ツール
検索
ツール
文書の表示
以前のリビジョン
バックリンク
最近の変更
メディアマネージャー
サイトマップ
ログイン
>
最近の変更
メディアマネージャー
サイトマップ
トレース:
seminar:h_on_co
この文書は読取専用です。文書のソースを閲覧することは可能ですが、変更はできません。もし変更したい場合は管理者に連絡してください。
====== コバルト(0001)表面の水素原子吸着 ====== <code python> from math import sqrt a = 2.5 c = 1.622*a a1 = (2*a, 0, 0 ) a2 = (a, a*sqrt(3), 0 ) a3 = (0, 0, 3*c) from ase import Atom, Atoms from ase.calculators.jacapo import Jacapo from ase.dft.kpoints import cc18_1x1 slab = Atoms(pbc = True) slab.append(Atom('Co', (0, 0, 0), magmom=1.6)) slab.append(Atom('Co', (1/2., 0, 0), magmom=1.6)) slab.append(Atom('Co', (0, 1/2., 0), magmom=1.6)) slab.append(Atom('Co', (1/2., 1/2., 0), magmom=1.6)) slab.append(Atom('Co', (1/6., 1/6., -1/6.), magmom=1.6)) slab.append(Atom('Co', (2/3., 1/6., -1/6.), magmom=1.6)) slab.append(Atom('Co', (1/6., 2/3., -1/6.), magmom=1.6)) slab.append(Atom('Co', (2/3., 2/3., -1/6.), magmom=1.6)) slab.set_cell([a1, a2, a3], scale_atoms = True) slab.append(Atom('H', (0, 0, 1.4598))) para = {} para.update(kpts = cc18_1x1) # set the k-points (Chadi-Cohen) para.update(pw = 340) # planewavecutoff in eV para.update(dw = 340) para.update(nbands = 10 + 8*6 + 1*1) # set the number of electronic bands para.update(spinpol = True) # this calculation should be spinpolarized para.update(symmetry = True) # use symmetry to reduce the k-point set para.update(atoms = slab) calc = Jacapo('H_Co_ontop.nc', **para) energy = calc.get_potential_energy() print 'energy = ', energy </code> <code sh> energy = -8741.97691323 </code> <code sh> $ nc2cube H_Co_ontop.nc H_Co_ontop.cube </code> {{:seminar:h_co_ontop.png|H/Co(0001)}} [[seminar:jacapo_sample|目次へもどる]]
seminar/h_on_co.txt
· 最終更新: 2022/08/23 13:34 by
127.0.0.1
ページ用ツール
文書の表示
以前のリビジョン
バックリンク
文書の先頭へ