| 查看: 851 | 回复: 1 | ||
| 【悬赏金币】回答本帖问题，作者Emanonvvv将赠送您 5 个金币 | ||
[求助]
有用GPAW的朋友吗，算声子一般怎么设置参数 已有1人参与
|
||
|
计算钙钛矿2X2X2超胞的声子谱，程序跑很久也不出结果，查看文件夹多了很多后缀为pckL的文件 phonon.1x+.pckl phonon.2z-.pckl phonon.4x+.pckl phonon.5z-.pckl phonon.7x+.pckl p4mm.cif phonon.1y-.pckl phonon.2z+.pckl phonon.4y-.pckl phonon.5z+.pckl phonon.7y-.pckl phonon.0x-.pckl phonon.1y+.pckl phonon.3x-.pckl phonon.4y+.pckl phonon.6x-.pckl phonon.7y+.pckl phonon.0x+.pckl phonon.1z-.pckl phonon.3x+.pckl phonon.4z-.pckl phonon.6x+.pckl phonon.7z-.pckl phonon.0y-.pckl phonon.1z+.pckl phonon.3y-.pckl phonon.4z+.pckl phonon.6y-.pckl phonon.7z+.pckl phonon.0y+.pckl phonon.2x-.pckl phonon.3y+.pckl phonon.5x-.pckl phonon.6y+.pckl phonon.eq.pckl phonon.0z-.pckl phonon.2x+.pckl phonon.3z-.pckl phonon.5x+.pckl phonon.6z-.pckl phonon.0z+.pckl phonon.2y-.pckl phonon.3z+.pckl phonon.5y-.pckl phonon.6z+.pckl phonon.1x-.pckl phonon.2y+.pckl phonon.4x-.pckl phonon.5y+.pckl phonon.7x-.pckl 这是我仿写ASE教程的计算程序： from ase import atoms from ase.io import read,write from gpaw import GPAW, FermiDirac import matplotlib.pyplot as plt from ase.phonons import Phonons at = read('mm.cif') calc = GPAW(kpts=(4,4,4), h=0.2, occupations=FermiDirac(0.)) ph = Phonons(at, calc, supercell = (1,1,1), delta = 0.05) ph.run() ph.read(method='frederiksen', acoustic=True) ph.clean() path = ph.cell.bandpath('GXMGRX,MR', npoints=100) bs = ph.get_band_structure(path) dos = ph.get_dos(spin=0, npts=1000) fig = plt.figure(1, figsize=(7, 4)) ax = fig.add_axes([.12, .07, .67, .85]) emax = 0.035 bs.plot(ax=ax, emin=0.0, emax=emax) dosax = fig.add_axes([.8, .07, .17, .85]) dosax.fill_between(dos.weights[0], dos.energy, y2=0, color='grey', edgecolor='k', lw=1) dosax.set_ylim(0, emax) dosax.set_yticks([]) dosax.set_xticks([]) dosax.set_xlabel("DOS", fontsize=18) fig.show() fig.savefig('phonon.png') 有没有用GPAW的胖友帮忙看下是哪里出了问题。 |
回复此楼» 猜你喜欢
281求调剂（0805）
已经有8人回复
-
环境领域全国重点实验室招收博士1-2名
已经有3人回复
-
材料专硕306英一数二
已经有10人回复
-
301求调剂
已经有6人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有7人回复
-
302求调剂
已经有6人回复
-
26博士申请
已经有3人回复
-
268求调剂
已经有3人回复
-
311求调剂
已经有10人回复
-
被我言中：新模板不强调格式了，假专家开始管格式了
已经有4人回复
get-it
新虫 (小有名气)
- 应助: 23 (小学生)
- 金币: 1825.8
- 散金: 26
- 红花: 4
- 帖子: 273
- 在线: 448.3小时
- 虫号: 3212202
- 注册: 2014-05-16
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2021-04-23 16:19:04
