| 查看: 2351 | 回复: 24 | ||||
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 / |
回复此楼» 猜你喜欢
复试调剂
已经有14人回复
-
化学070300 求调剂
已经有5人回复
-
一志愿085502，267分求调剂
已经有13人回复
-
277求调剂
已经有22人回复
-
268分085602化学工程调剂
已经有29人回复
-
生物学308求调剂
已经有5人回复
-
一志愿哈工大 085600 277 12材科基求调剂
已经有29人回复
-
2本，初试303，0860求调剂
已经有8人回复
-
279学硕食品专业求调剂院校
已经有16人回复
-
材料工程085601，270求调剂
已经有34人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 为何用晶格常数的实验值计算声子？
已经有18人回复
- 关于用超胞方法计算声子的一些经验
已经有21人回复
- Abinit计算声子问题
已经有10人回复
- ZnO的声子计算
已经有18人回复
- MS计算晶体声子谱的时候，横坐标有负能量值，为什么？
已经有11人回复
- 自旋体系的声子计算应该用什么软件呢？
已经有5人回复
- 【求助】castep计算声子谱，为什么不能进行振动分析
已经有9人回复
- 【求助】使用vasp进行带心冻结声子计算的几个问题~~~
已经有9人回复
- 【求助】请教CASTEP计算声子谱后出现虚频
已经有6人回复
- 【求助】有人用有限差分方法计算过G点声子谱吗？
已经有10人回复
- 【求助】关于PWSCF介电常数和声子谱的计算
已经有7人回复
- 【求助】请问大家在做声子计算时，这个 qpt是按照什么设置的？？？？
已经有9人回复
- 【求助】用MS中的castep计算声子振动模频率，具体怎么算
已经有10人回复
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
2楼2011-10-26 16:38:30
hanyanli0475
铜虫 (小有名气)
- 应助: 11 (小学生)
- 金币: 173.1
- 散金: 32
- 红花: 3
- 帖子: 296
- 在线: 130.5小时
- 虫号: 1160645
- 注册: 2010-12-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2011-10-26 17:24:28
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
4楼2011-10-26 18:45:29
hanyanli0475
铜虫 (小有名气)
- 应助: 11 (小学生)
- 金币: 173.1
- 散金: 32
- 红花: 3
- 帖子: 296
- 在线: 130.5小时
- 虫号: 1160645
- 注册: 2010-12-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2011-10-26 19:45:51
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
【答案】应助回帖
★ ★ ★
youzhizhe(金币+3): 谢谢交流。已计入专家考核存档。 2011-10-27 09:35:34
youzhizhe(金币+3): 谢谢交流。已计入专家考核存档。 2011-10-27 09:35:34
|
测试很简单，就是把截断从小到达取不同数值进行自洽场迭代计算，看不同截断的时候的能量收敛精度。如果基本前后两步的数值基本一致，则可认为收敛。能量或者压强的收敛与截断、K点选择等有关系，你应该不同的参数分别测试。 PS：很多时候我们计算的时候都是不管收不收敛，直接给定数值就进行计算，这是完全不对的，不是科学的态度。我们应该一步一步的来，多看看别人给你的步骤和例子，里面都是有的。 自己想当然的算，然后靠别人帮你指出错误，尤其是很简单又讨论过无数次的错误，这对自己而言虽然轻松了，但少了一个自己探索的过程。而且简单重复的回答，也浪费了别人的时间。归根结底，是自身没有足够的学习和调研造成许多根本上的却很无谓的错误。 说多了。 |
» 本帖已获得的红花（最新10朵）
hanyanli04756楼2011-10-26 22:36:37
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
★ ★ ★
youzhizhe(金币+3): 谢谢交流。已计入专家考核存档。 2011-10-27 09:35:46
youzhizhe(金币+3): 谢谢交流。已计入专家考核存档。 2011-10-27 09:35:46
|
你这么说也不管用啊，毕竟测试也是个比较耗时间精力的过程，并且还可能证明了之前看起来跟实验非常“符合”的结果有问题了~ 事实上去年我在实验室里做组会报告时，做的就是如何测试收敛的问题，被组里几个也做计算的人说是很无聊和没有特殊的东西的报告~ 结果，有人的计算就是低截断低K点得到了跟实验很相近的结果，而我在测试他的模型时，加大了K点和截断，于是就跟实验值差别非常大了~ 不过，由于老板不懂计算，也没有这样的科学思想，所以，觉得我加大截断和K点的做法是没选对参数，而那个哥们做的结果符合实验，是好结果…… 总的来说，工科的急功近利的思想啊~ |
7楼2011-10-26 23:26:21
hanyanli0475
铜虫 (小有名气)
- 应助: 11 (小学生)
- 金币: 173.1
- 散金: 32
- 红花: 3
- 帖子: 296
- 在线: 130.5小时
- 虫号: 1160645
- 注册: 2010-12-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
8楼2011-10-27 08:17:00
huazhorg
铁杆木虫 (正式写手)
- 应助: 21 (小学生)
- 金币: 4953.8
- 散金: 5
- 红花: 2
- 帖子: 381
- 在线: 363.2小时
- 虫号: 629899
- 注册: 2008-10-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
9楼2011-10-27 10:52:50
huazhorg
铁杆木虫 (正式写手)
- 应助: 21 (小学生)
- 金币: 4953.8
- 散金: 5
- 红花: 2
- 帖子: 381
- 在线: 363.2小时
- 虫号: 629899
- 注册: 2008-10-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
10楼2011-10-27 10:59:09
