24小时热门版块排行榜    

查看: 3147  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

peter10373

新虫 (初入文坛)

[求助] 新手计算晶格常数上的问题 已有1人参与

С???????????VASP???ú???????????????QE?????????????Ц?????λ???????????
С????VASP????????????????????:
???????????趨??????????????

??????????????????????vc-relax????????????????????????????????????????????????input?????鷳??λ????械?????
&control                                                            
    calculation = 'vc-relax'                                       
    restart_mode = 'from_scratch'                                   
    prefix = 'MoO3'                                                
    outdir = '\Documents\peter\QE\test\MoO3'                        
    pseudo_dir = '/aracbox/app/espresso-5.1/pseudo/'               
/                                                                  
&system                                                            
    ibrav= 8                                                        
    a = 7.5371999741                                                
    b = 7.7888998985                                                
    c = 18.2537994385                                               
    nat  = 32                                                      
    ntyp = 2                                                        
    ecutwfc = 30                                                   
/                                                                  
&electrons                                                         
    mixing_mode = "plain",                                          
    electron_maxstep = 100,                                         
    conv_thr = 1D-6,                                                
/                                                                  
&IONS                                                               
    ion_dynamics="bfgs",                                            
/                                                                  
&CELL                                                               
    cell_factor = 4.0D0                                             
    cell_dofree = xy                                                
/                                                                  
ATOMIC_SPECIES                                                      
Mo   95.94   Mo.pbe-sp-hgh.UPF                                    
O    16      O.pbe-hgh.UPF                                         
ATOMIC_POSITIONS {angstrom}                                         
Mo                 2.82645000    2.18216938    1.63497461         
Mo                 0.94215000    1.71265449    4.61822965         
Mo                 6.59505000    2.18216938    1.63497461         
Mo                 4.71075000    1.71265449    4.61822965         
Mo                 2.82645000    6.07661938    1.63497461         
Mo                 0.94215000    5.60711228    4.61822965         
Mo                 6.59505000    6.07661159    1.63497461         
Mo                 4.71075000    5.60711228    4.61822965         
O                  0.94215000    1.95893951    2.17501329         
O                  2.82645000    1.93554165    4.07939573         
O                  2.82645000    0.07471892    1.85987968         
O                  0.94215000    3.81871074    4.39540552         
O                  2.82645000    2.06826451   18.20071795         
O                  0.94215000    1.82679303    6.30535537         
O                  4.71075000    1.95893951    2.17501329         
O                  6.59505000    1.93554165    4.07939573         
O                  6.59505000    0.07471892    1.85987968         
O                  4.71075000    3.81871074    4.39540552         
O                  6.59505000    2.06826451   18.20071795         
O                  4.71075000    1.82679303    6.30535537         
O                  0.94215000    5.85338951    2.17501329         
O                  2.82645000    5.82999165    4.07941399         
O                  2.82645000    3.96916892    1.85987968         
O                  0.94215000    7.71316853    4.39540552         
O                  2.82645000    5.96271451   18.20071795         
O                  0.94215000    5.72124303    6.30535537         
O                  4.71075000    5.85338951    2.17501329         
O                  6.59505000    5.82999165    4.07941399         
O                  6.59505000    3.96916113    1.85987968         
O                  4.71075000    7.71316853    4.39540552         
O                  6.59505000    5.96271451   18.20071795         
O                  4.71075000    5.72124303    6.30537363         
K_POINTS {automatic}                                                
3 1 1 0 0 0                                                         



output?????????CELL_PARAMETER???????????????
CELL_PARAMETERS (alat= 14.24324376)         
   0.693664469   0.000000000   0.000000000   
   0.000000000   0.768702659   0.000000000   
   0.000000000   0.000000000   2.421827668   
???????????????д?
a = 0.693664469*14.24324376*0.529177 = 5.22828 angstrom  (??????7.5371999741 ??Щ???)
???????b,c


????韩??λ????????!!С???????????????]???????????????λ??
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

KalaShayminS

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
peter10373: 金币+5, ★★★很有帮助 2014-08-25 22:35:47
引用回帖:
3楼: Originally posted by peter10373 at 2014-08-07 10:29:26
你好,谢谢回应,小弟用这套算完之后发现结构x和y方向都被压缩很多,使得结构有点变形,原子排列也与先前的差异有点多,这边不晓得是不是dofree这个参数的设置问题。
请问要去哪搜寻hgh的赝势大概是用多少Ry呢?
我 ...

你这个模型是slab模型吧,我的问题是你的体相是否优化过?
需要测试ecutwfn和k点在体相优化时收敛,做slab计算才有意义。
5楼2014-08-08 01:54:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

KalaShayminS

铁杆木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
这套参数用来优化体相结构的效果怎么样?
hgh的赝势我记得ecut=30是不够的
另外这个3 1 1 的k点也很奇怪,xy方向的k点应该一样多才对。
2楼2014-08-05 21:21:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peter10373

新虫 (初入文坛)

你好,谢谢回应,小弟用这套算完之后发现结构x和y方向都被压缩很多,使得结构有点变形,原子排列也与先前的差异有点多,这边不晓得是不是dofree这个参数的设置问题。
请问要去哪搜寻hgh的赝势大概是用多少Ry呢?
我在算VASP时是用PAW当时是用400ev大概是30Ry

K点问题感谢提醒,目前已修正但算出来还是有些怪异

谢谢!
3楼2014-08-07 10:29:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pangrui1985

铜虫 (小有名气)

引用回帖:
3楼: Originally posted by peter10373 at 2014-08-07 10:29:26
你好,谢谢回应,小弟用这套算完之后发现结构x和y方向都被压缩很多,使得结构有点变形,原子排列也与先前的差异有点多,这边不晓得是不是dofree这个参数的设置问题。
请问要去哪搜寻hgh的赝势大概是用多少Ry呢?
我 ...

一般YS文件里面有,有PAW的还是选PAW比较好。hgh是模守恒的,截断200多差不多。
4楼2014-08-07 23:01:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +3 wangyl_123 2026-10-10 6/300 2026-10-11 15:04 by Abin0406
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 5/250 2026-10-11 09:00 by 好运来~9527
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +5 13108017953 2026-10-08 5/250 2026-10-10 11:15 by wdwd123
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见