24小时热门版块排行榜    

CyRhmU.jpeg
查看: 4504  |  回复: 23

gemucai

木虫 (正式写手)

刚用occupation=fixed算了一遍,跟取比较大的smearing效果是一样一样滴,这是怎么回事呢?
11楼2012-05-26 11:17:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

姚海峰456

银虫 (初入文坛)

已经用第一性原理优化好的单胞,用gulp计算声子谱会出现虚频,怎么解决?——求大神帮忙啊!!!
12楼2012-08-01 21:42:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

引用回帖:
12楼: Originally posted by 姚海峰456 at 2012-08-01 21:42:17
已经用第一性原理优化好的单胞,用gulp计算声子谱会出现虚频,怎么解决?——求大神帮忙啊!!!

我记得gulp是经典的紧束缚方法吧?朝出现虚频的方向改变原子位置试试。
13楼2012-08-02 18:56:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

引用回帖:
13楼: Originally posted by gemucai at 2012-08-02 18:56:11
我记得gulp是经典的紧束缚方法吧?朝出现虚频的方向改变原子位置试试。...

GULP是紧束缚么……我怎么记得GULP是类似于经典的分子动力学的东西呢……记得实验室的人用这货跑上千个原子的分子动力学,也就一天的事~还是在一台三年前的四核电脑上……不过,至于他参数是不是足够收敛,我就不清楚了……当初那哥们就是为了凑毕业论文的数据的,用GULP出图快~
转个官方的说明:
CODE:
GULP is a program for performing a variety of types of simulation on materials using boundary conditions of 0-D (molecules and clusters), 1-D (polymers), 2-D (surfaces, slabs and grain boundaries), or 3-D (periodic solids). The focus of the code is on analytical solutions, through the use of lattice dynamics, where possible, rather than on molecular dynamics. A variety of force fields can be used within GULP spanning the shell model for ionic materials, molecular mechanics for organic systems, the embedded atom model for metals and the reactive REBO potential for hydrocarbons. Analytic derivatives are included up to at least second order for most force fields, and to third order for many.

思想重于技巧,内涵重于表象
14楼2012-08-02 20:39:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yinwei_li

木虫 (正式写手)

引用回帖:
3楼: Originally posted by gemucai at 2012-05-25 18:59:36
问题解决了。

跟大家分享一下我的解决方法,以后碰到有虚频的情况不妨增加smearing的宽度,即degauss的大小,效果很灵的。其实有时候根本不需要过度地优化,我优化到0.0000001ry/bohr也有很大虚频,但是degauss增 ...

你如何确保你增加smearing后计算的结果是精确的呢?
与其临渊羡鱼,不如退而结网!
15楼2012-08-02 20:49:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

内容已删除
16楼2012-08-03 09:12:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

引用回帖:
14楼: Originally posted by souledge at 2012-08-02 20:39:39
GULP是紧束缚么……我怎么记得GULP是类似于经典的分子动力学的东西呢……记得实验室的人用这货跑上千个原子的分子动力学,也就一天的事~还是在一台三年前的四核电脑上……不过,至于他参数是不是足够收敛,我就不 ...

其实这不是真正的紧束缚,只是一种经典的紧束缚,类似于phonon里面的小位移方法。它通过算势能对位移的二阶偏导得出近邻次近邻相互作用的动力学矩阵,解久期方程得到振动频率。貌似这个势能确定了之后算起来确实相当快,我去年手算一简单体系也只花了一个下午。
17楼2012-08-03 09:16:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yinwei_li

木虫 (正式写手)

引用回帖:
16楼: Originally posted by gemucai at 2012-08-03 09:12:12
问得非常好,这是一个我没搞清楚的问题。我曾经保持k点不变,改变展宽,结果是不收敛的。有大虾说一般是先定展宽,再增加k点达到收敛。当然对于我这个绝缘体系,fixed就足够了,但我也很想知道对于导体而言,展宽为 ...

其实,展宽的参数也是需要测试的,在其他参数固定的情况下,改变smearing数值,当声子频率收敛时的smearing就是精确地,
与其临渊羡鱼,不如退而结网!
18楼2012-08-03 14:11:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

内容已删除
19楼2012-08-03 16:18:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘小伟110

木虫 (小有名气)

楼主好人啊,我有个问题就是我要做零点能校正的话,按E=1/2hv公式,把所有频率的能量相加,大约就是0.7eV,按理说,零点能校正的话,不会是这么大的值,也就是0.02eV,是不是这中计算方法有误呢?(我计算的是单个分子的频率)
最近甚是苦恼啊,还望您能指点迷津,期待您的回复,在此,万分感谢。。。
慢慢的思考。。。
20楼2014-01-02 16:24:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gemucai 的主题更新
信息提示
请填处理意见