24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者Emanonvvv将赠送您 5 个金币

Emanonvvv

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 28.8
帖子: 8
在线: 3.8小时
虫号: 17213403
注册: 2019-09-04
性别: GG
专业: 物理电子学

[求助] 有用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的胖友帮忙看下是哪里出了问题。

回复此楼