| 查看: 939 | 回复: 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的胖友帮忙看下是哪里出了问题。 |
回复此楼» 猜你喜欢
2027广东省杰青
已经有5人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有72人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
get-it
新虫 (小有名气)
- 应助: 23 (小学生)
- 金币: 1825.8
- 散金: 26
- 红花: 4
- 帖子: 273
- 在线: 448.3小时
- 虫号: 3212202
- 注册: 2014-05-16
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2021-04-23 16:19:04
