差分

この文書の現在のバージョンと選択したバージョンの差分を表示します。

seminar:al_eos 2014/02/01 10:40 — 現在
ライン 1: ライン 1:
-====== アルミニウム結晶の状態方程式(エネルギーの体積依存性) ====== 
-<code python> 
-def get_fcc_unitcell(p): 
-  c = [(  0.0, p/2.0, p/2.0), 
-      (p/2.0,  0.0, p/2.0), 
-      (p/2.0, p/2.0,  0.0)] 
-  return c 
- 
-b = 4.05 
- 
-from ase import Atom, Atoms 
-from ase.calculators.jacapo import Jacapo 
-bulk = Atoms([Atom('Al', (0, 0, 0))], pbc = True) 
- 
-para = {} 
-para.update(kpts = (4, 4, 4)) 
-para.update(nbands = 6) 
-para.update(ft = 0.01) 
-     
-print 'Volume(A^3) Energy(eV)' 
- 
-for fraction in [0.9, 0.95, 1.0, 1.05, 1.1]: 
-    ncfile = 'bulkAl_%1.2f.nc' % fraction 
-    new_cell = get_fcc_unitcell(b*(fraction**(1./3.))) 
-    bulk.set_cell(new_cell, scale_atoms = True) 
-    para.update(atoms = bulk) 
-    solver = Jacapo(ncfile, **para) 
-    vol = bulk.get_volume() 
-    eng = bulk.get_potential_energy() 
-    print vol, eng 
-</code> 
-<code sh> 
-$ python Al_equation_of_state.py  
-Volume(A^3) Energy(eV) 
-14.946778125 -56.4090488622 
-15.7771546875 -56.4473195776 
-16.60753125 -56.4583266527 
-17.4379078125 -56.4483594784 
-18.268284375 -56.4233507908 
-$  
-</code> 
-{{:seminar:al.png?512|アルミニウム結晶のエネルギーの体積依存性}} 
- 
-[[seminar:jacapo_sample|目次へもどる]] 
seminar/al_eos.1391218840.txt.gz · 最終更新: 2014/02/01 10:40 by kimi
www.chimeric.de Creative Commons License Valid CSS Driven by DokuWiki do yourself a favour and use a real browser - get firefox!! Recent changes RSS feed Valid XHTML 1.0