【答案】应助回帖

» 猜你喜欢

售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急已经有13人回复
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急已经有12人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有22人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有19人回复
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急已经有7人回复
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急已经有4人回复
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急已经有7人回复
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有9人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有10人回复
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急已经有7人回复

1楼 2016-08-02 15:28:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

回帖置顶 ( 共有3个 )

lfhuang

木虫 (著名写手)

应助: 39 (小学生)
金币: 3660.6
散金: 2321
红花: 37
帖子: 1432
在线: 681.5小时
虫号: 505824
注册: 2008-02-17
性别: GG
专业: 金属材料表面科学与工程

【答案】应助回帖

★
感谢参与，应助指数 +1
蛋蛋小童鞋: 金币+1, 谢谢参与 2016-08-11 08:37:38
漫天飘雪: 回帖置顶 2018-06-24 13:12:01

按照如下文章里的参数，肯定能得到好的声子谱：
http://journals.aps.org/prb/abstract/10.1103/PhysRevB.91.205433

尤其你的conv_thr还太大，建议调到d-12。nq应该还太小，一般在8以上。

先进使役材料计算（https://aidme.nimte.ac.cn）

14楼2016-08-04 00:21:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lfhuang

木虫 (著名写手)

应助: 39 (小学生)
金币: 3660.6
散金: 2321
红花: 37
帖子: 1432
在线: 681.5小时
虫号: 505824
注册: 2008-02-17
性别: GG
专业: 金属材料表面科学与工程

【答案】应助回帖

漫天飘雪: 回帖置顶 2018-06-24 13:12:03

另外，声子计算用什么参数，结构优化也最好用一样严格的参数，另外，力的收敛最好达到d-4。

先进使役材料计算（https://aidme.nimte.ac.cn）

15楼2016-08-04 00:24:05

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lfhuang

木虫 (著名写手)

应助: 39 (小学生)
金币: 3660.6
散金: 2321
红花: 37
帖子: 1432
在线: 681.5小时
虫号: 505824
注册: 2008-02-17
性别: GG
专业: 金属材料表面科学与工程

【答案】应助回帖

漫天飘雪: 回帖置顶 2018-06-24 13:12:05

当你nq是3的倍数或足够密的的时候，你自然会得到简并了。

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

坏逊平儿

先进使役材料计算（https://aidme.nimte.ac.cn）

16楼2016-08-04 00:26:00

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

回帖支持 ( 显示支持度最高的前 50 名 )

liqizuiyang

木虫 (著名写手)

应助: 59 (初中生)
金币: 6859.5
散金: 261
红花: 94
帖子: 2093
在线: 708.9小时
虫号: 2366686
注册: 2013-03-21
性别: GG
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★
蛋蛋小童鞋: 金币+2, 应助指数+1, 鼓励回答问题详细者 2016-08-03 14:34:12

重新优化结构，力收敛判据调至0.0001eV/A。

改用MP型展宽，做严格收敛测试。

K点改成18*18*1.

diagonalization改成david。cg只在空带数很大的时候才有优势，速度上并不占优。

6楼2016-08-03 09:37:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

c19940408

银虫 (初入文坛)

应助: 7 (幼儿园)
金币: 257.4
帖子: 42
在线: 29.9小时
虫号: 2486278
注册: 2013-05-28
专业: 理论和计算化学

【答案】应助回帖

★
蛋蛋小童鞋: 金币+1, 谢谢回答 2016-08-11 08:37:15

建议使用和优化时相近的参数来计算。。比如相同的能量截断，k点等等

20楼2016-08-08 02:14:24

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

普通回帖

坏逊平儿

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 282
帖子: 20
在线: 32.7小时
虫号: 3871845
注册: 2015-05-15
性别: GG
专业: 电磁场与电路

scf.in参数设置

&control
calculation='scf',
restart_mode='from_scratch',
prefix='si'
pseudo_dir = '/home/bxu/zjp/2sixi/4/si/',
outdir='/home/bxu/zjp/2sixi/4/si/tmp/'
/
&system
ibrav = 4, celldm(1) =7.3052690583, celldm(3) =3.0593725672, nat= 2, ntyp= 1,
ecutwfc = 50,ecutrho=450
occupations='smearing',degauss=0.02,smearing='gaussian'
/
&electrons
electron_maxstep=200
mixing_beta = 0.5
conv_thr = 1.0d-10
diagonalization='cg'
/
ATOMIC_SPECIES
Si 28.086 Si.pbe-n-rrkjus_psl.1.0.0.UPF
ATOMIC_POSITIONS {crystal}
Si 0.333330005 0.666660011 0.480724089
Si 0.666660011 0.333330005 0.519265927
K_POINTS {automatic}
20 20 1 0 0 0

2楼2016-08-02 15:31:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

坏逊平儿

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 282
帖子: 20
在线: 32.7小时
虫号: 3871845
注册: 2015-05-15
性别: GG
专业: 电磁场与电路

ph.in参数设置

&inputph
tr2_ph=1.0d-16,
prefix='si',
amass(1)=28.086,
outdir='/home/bxu/zjp/2sixi/4/si/tmp',
fildyn='si.dyn',
elph=.false.,
trans=.true.,
ldisp=.true.
nq1=4, nq2=4,nq3=1
/

3楼2016-08-02 15:32:17

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

坏逊平儿

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 282
帖子: 20
在线: 32.7小时
虫号: 3871845
注册: 2015-05-15
性别: GG
专业: 电磁场与电路

q2r.in 参数设置

&input
zasr='simple',fildyn='si.dyn',flfrc='si.fc',la2F=.false.

4楼2016-08-02 15:32:27

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

golddoushi

木虫 (正式写手)

应助: 36 (小学生)
金币: 4007.2
散金: 25
红花: 15
帖子: 337
在线: 447.7小时
虫号: 2546370
注册: 2013-07-15
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

蛋蛋小童鞋: 应助指数+1, 谢谢参与 2016-08-03 08:42:11

gamma 点本来就有虚频，对应长波极限。一般直接忽略不计

5楼2016-08-03 07:04:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

坏逊平儿

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 282
帖子: 20
在线: 32.7小时
虫号: 3871845
注册: 2015-05-15
性别: GG
专业: 电磁场与电路

引用回帖:

5楼: Originally posted by golddoushi at 2016-08-03 07:04:41
gamma 点本来就有虚频，对应长波极限。一般直接忽略不计

大神，你应该搞错了。。。问题不是在gamma点，是在gamma点附近的出现较大虚频，这跟文献上明显不对应啊。其次，在另一个好对称点处两光学支没有简并。是什么原因呢

7楼2016-08-03 11:01:25

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

坏逊平儿

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 282
帖子: 20
在线: 32.7小时
虫号: 3871845
注册: 2015-05-15
性别: GG
专业: 电磁场与电路

引用回帖:

6楼: Originally posted by liqizuiyang at 2016-08-03 09:37:36
重新优化结构，力收敛判据调至0.0001eV/A。

改用MP型展宽，做严格收敛测试。

K点改成18*18*1.

diagonalization改成david。cg只在空带数很大的时候才有优势，速度上并不占优。

请问改用mp展宽，做测试，能具体讲下细节吗，刚入门，诸多不太懂。。

发自小木虫Android客户端

8楼2016-08-03 11:27:34

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)

应助: 59 (初中生)
金币: 6859.5
散金: 261
红花: 94
帖子: 2093
在线: 708.9小时
虫号: 2366686
注册: 2013-03-21
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

8楼: Originally posted by 坏逊平儿 at 2016-08-03 11:27:34
请问改用mp展宽，做测试，能具体讲下细节吗，刚入门，诸多不太懂。。
...

可以参考vasp手册。测试流程就是不断增大展开，同时平均到每个原子上的TS项低于某个数值。

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

坏逊平儿

9楼2016-08-03 14:21:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

坏逊平儿

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 282
帖子: 20
在线: 32.7小时
虫号: 3871845
注册: 2015-05-15
性别: GG
专业: 电磁场与电路

送红花一朵

引用回帖:

9楼: Originally posted by liqizuiyang at 2016-08-03 14:21:59
可以参考vasp手册。测试流程就是不断增大展开，同时平均到每个原子上的TS项低于某个数值。...

谢谢耐心指导