应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 声子计算又负频的出现！！！
251/3返回列表123下一页
查看: 2129  |  回复: 24
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

hanyanli0475

铜虫 (小有名气)

[求助] 声子计算又负频的出现！！！

今天计算氮化铟的声子，输出后低频有-26cm-1的负频，而且其他低频的走势和文献也不太一样，我把输入文件粘出来，请大家帮我找找错误呗，看看有没有需要修改的地方，谢谢了！！！
自洽计算：
&CONTROL
calculation = 'scf' ,
restart_mode = 'from_scratch' ,
outdir = '/home/zhao/10.21/tmp/' ,
pseudo_dir = '/opt/espresso-4.2.1/pseudo/' ,
prefix = 'InN' ,
tprnfor= .true. ,
tstress= .true. ,
/
&SYSTEM
ibrav = 4,celldm(1) = 6.893,celldm(3) = 1.605,nat = 4,ntyp = 2,
ecutwfc = 45.0 ,ecutrho = 180 ,
/
&ELECTRONS
conv_thr = 1.0d-8 ,
mixing_beta = 0.7 ,
/
ATOMIC_SPECIES
In   114.820  In.pbe-d-rrkjus.UPF
N    14.007   N.pbe-rrkjus.UPF
ATOMIC_POSITIONS crystal
In       0.333333487   0.666666513   0.003058286    1   1   1
In       0.666666513   0.333333487   0.503058286    1   1   1
N        0.333333439   0.666666561   0.381941714    1   1   1
N        0.666666561   0.333333439   0.881941714    1   1   1
K_POINTS automatic
4 4 4      1 1 1
声子计算
phonon for InN
&inputph
outdir ='/home/zhao/10.21/tmp/',
tr2_ph =1.0d-12,
prefix ='InN',
fildvscf ='InNdv',
amass(1) =114.8200,
amass(2) =14.00700,
fildyn ='InN.dyn',
elph =.false.,
trans =.true.,
ldisp =.true.,
nq1=4,nq2=4,nq3=4
/
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
什么时候能算完！！！！
1楼2011-10-26 15:05:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★
贺仪(金币+2): 多谢指教！ 2011-10-26 21:59:43
4x4x4的k点也敢拿来做phonon计算啊？太大无畏了～
顺便，对于模守赝势，这点截断也忒低了～用超软做声子时这个截断还差不多。
另外，对于模守赝势，不要自己设置电荷密度截断能。
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
2楼2011-10-26 16:38:30
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by souledge at 2011-10-26 16:38:30:
4x4x4的k点也敢拿来做phonon计算啊？太大无畏了～
顺便，对于模守赝势，这点截断也忒低了～用超软做声子时这个截断还差不多。
另外，对于模守赝势，不要自己设置电荷密度截断能。

谢谢您的回答，我现在就去改改，重新跑一下。
今天下午把ph的源程序大概看看，真难啊，呵呵！
赞一下 回复此楼
什么时候能算完！！！！
3楼2011-10-26 17:24:28
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

【答案】应助回帖

★ ★
贺仪(金币+2): 多谢指教！ 2011-10-26 22:00:00
他用的貌似是超软的。
超软赝势ecutrho至少要是ecutwfc的8-10倍。对于rrkjus类型的，ecutwfc=45确实有些小。建议做一个能量收敛测试。
此外，4*4*4也有点小，也要做收敛测试。
赞一下(1人) 回复此楼
4楼2011-10-26 18:45:29
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)
引用回帖:
4楼: Originally posted by goldenfisher at 2011-10-26 18:45:29:
他用的貌似是超软的。
超软赝势ecutrho至少要是ecutwfc的8-10倍。对于rrkjus类型的，ecutwfc=45确实有些小。建议做一个能量收敛测试。
此外，4*4*4也有点小，也要做收敛测试。

请问如何做个测试呢，测试的标准怎么看呢？我以前没做过测试，请您把测试的过程说的仔细点呗，谢谢了。
赞一下 回复此楼
什么时候能算完！！！！
5楼2011-10-26 19:45:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

【答案】应助回帖

★ ★ ★
youzhizhe(金币+3): 谢谢交流。已计入专家考核存档。 2011-10-27 09:35:34
引用回帖:
5楼: Originally posted by hanyanli0475 at 2011-10-26 19:45:51:
请问如何做个测试呢，测试的标准怎么看呢？我以前没做过测试，请您把测试的过程说的仔细点呗，谢谢了。

测试很简单，就是把截断从小到达取不同数值进行自洽场迭代计算，看不同截断的时候的能量收敛精度。如果基本前后两步的数值基本一致，则可认为收敛。能量或者压强的收敛与截断、K点选择等有关系，你应该不同的参数分别测试。
PS：很多时候我们计算的时候都是不管收不收敛，直接给定数值就进行计算，这是完全不对的，不是科学的态度。我们应该一步一步的来，多看看别人给你的步骤和例子，里面都是有的。
自己想当然的算，然后靠别人帮你指出错误，尤其是很简单又讨论过无数次的错误，这对自己而言虽然轻松了，但少了一个自己探索的过程。而且简单重复的回答，也浪费了别人的时间。归根结底，是自身没有足够的学习和调研造成许多根本上的却很无谓的错误。
说多了。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

hanyanli0475 hanyanli0475
6楼2011-10-26 22:36:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
★ ★ ★
youzhizhe(金币+3): 谢谢交流。已计入专家考核存档。 2011-10-27 09:35:46
引用回帖:
6楼: Originally posted by goldenfisher at 2011-10-26 22:36:37:
测试很简单，就是把截断从小到达取不同数值进行自洽场迭代计算，看不同截断的时候的能量收敛精度。如果基本前后两步的数值基本一致，则可认为收敛。能量或者压强的收敛与截断、K点选择等有关系，你应该不同的参数 ...

你这么说也不管用啊，毕竟测试也是个比较耗时间精力的过程，并且还可能证明了之前看起来跟实验非常“符合”的结果有问题了~
事实上去年我在实验室里做组会报告时，做的就是如何测试收敛的问题，被组里几个也做计算的人说是很无聊和没有特殊的东西的报告~
结果，有人的计算就是低截断低K点得到了跟实验很相近的结果，而我在测试他的模型时，加大了K点和截断，于是就跟实验值差别非常大了~
不过，由于老板不懂计算，也没有这样的科学思想，所以，觉得我加大截断和K点的做法是没选对参数，而那个哥们做的结果符合实验，是好结果……
总的来说，工科的急功近利的思想啊~
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
7楼2011-10-26 23:26:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)
送鲜花一朵
引用回帖:
6楼: Originally posted by goldenfisher at 2011-10-26 22:36:37:
测试很简单，就是把截断从小到达取不同数值进行自洽场迭代计算，看不同截断的时候的能量收敛精度。如果基本前后两步的数值基本一致，则可认为收敛。能量或者压强的收敛与截断、K点选择等有关系，你应该不同的参数 ...

谢谢您的教导，学习了，您说的确实是我一直存在的问题，计算某个物理性质，把参数拿来就开始跑，总是抱着偶然结果很好这样的想法去做，现在知道错了，谢谢您，我会改正的。加油了！
赞一下 回复此楼
什么时候能算完！！！！
8楼2011-10-27 08:17:00
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)

【答案】应助回帖

★
youzhizhe(金币+1): 谢谢交流。 2011-10-27 23:30:30
引用回帖:
5楼: Originally posted by hanyanli0475 at 2011-10-26 19:45:51:
请问如何做个测试呢，测试的标准怎么看呢？我以前没做过测试，请您把测试的过程说的仔细点呗，谢谢了。

ecutwfc=45已经不小了，这里的单位是Ry, 45Ry就相当于600多eV 的截断，已经相当的高了
赞一下(1人) 回复此楼
9楼2011-10-27 10:52:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★
youzhizhe(金币+2): 谢谢交流。 2011-10-27 23:30:50
以他这么大的单胞，我认为8X8X8的K点就够了，振动出现虛频，有可能是结构没有优化好，也有可能是势能过软，换用较硬的赝势试试。另外，In原子较重，可以试试加U，再看看结果如何。
赞一下(1人) 回复此楼
10楼2011-10-27 10:59:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hanyanli0475 的主题更新
251/3返回列表123下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭