24小时热门版块排行榜    

CyRhmU.jpeg
查看: 563  |  回复: 3
当前主题已经存档。

zxinxin

金虫 (小有名气)

[交流] 【求助】frophon 计算声子为啥总是在G点有一小点虚频呢

有用frophon 计算声子的么
为什么用这个作出的声子 常常发现在Gama点处有一小点虚频呢 很小很小的
是因为计算中得到的力不准确么
怎么能消除这种假的虚频呢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)

zxinxin(金币+2): 2010-03-30 10:50
zxinxin(金币+1):谢谢 大约有 0.01 cm-1 一定是假虚 可能是力算的不准 2010-03-30 10:57
有多小?
如果只有几个波数
可以忽略
有可能是因为数值计算误差导致的
2楼2010-03-30 08:16:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

zxinxin(金币+2):我加了负方向的 看来还是软件本身的问题 我可以试试phonopy 谢谢 2010-03-30 12:31
主要的原因其实还是软件本身的问题,但是也能解决。
解决的办法,主要是加负方向的计算。
you can specify displacements more than independent freedome to obtain precise eigenvalus. This treatment works especially near G-point. An effective example is + - displacements. if your DISP file is,
1 1.00 0.00 0.00
1 0.00 1.00 0.00
1 0.00 0.00 1.00
you can create displacements by hand (or use sed), e.g.,
1  1.00  0.00  0.00
1 -1.00  0.00  0.00
1  0.00  1.00  0.00
1  0.00 -1.00  0.00
1  0.00  0.00  1.00
1  0.00  0.00 -1.00
This setting may be called 'plus-minus'. The plus-minus setting works to increase stability of phonon band structure when the initial atomic geometry is not enough converged to the equilibrium position. Especially frequency to approximately zero is sensitive to this treatment. Sometimes fails when forecs are converged too much.
毕竟能量的收敛是二阶的话,那么力收敛就是一阶。
frophon好像现在已经不升级了。他们编的现在是phonopy,全是用python写的,包括核心计算部分。这个还在不断升级中,而且在处理plus-minus的时候,软件会自己判断,因为有不需要plus-minus就能算准G point,它就不加,减少了计算量。有需要加plus-minus才能算准的,软件会自动加。我觉得比较好。
真理是一点点接近的!
3楼2010-03-30 12:00:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)

★ ★
gavinliu7390(金币+2):indeed.大部分的时候都是这样! 2010-04-03 16:13
0.01 cm-1 就是0
4楼2010-03-30 13:02:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zxinxin 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见