| 查看: 4895 | 回复: 10 | ||||
fdd096030079新虫 (小有名气)
|
[交流]
vasp中，做截断能测试时，测试收敛的标准 已有5人参与
|
|||
|
在做截断能测试时，测试收敛的标准是什么呢？ 我看到候博士的手册上说的时总能量相差小于0.001eV时，就说是已经收敛了，但我最近做的NaAlH4的超胞，96个原子，很难收敛到0.001eV，因此，我想是不是候博士那个地方说的不够清楚，应该是平均每个原子的能量相差0.001eV，而不是总能量。后来，我问导师，导师说是按照总能量的，但我觉得这样不合理啊~如何一个大的体系，要求总能量的变化在0.001eV，那不是很困难吗？本来体系大的话就难以计算，还哪能要求这么严格？因此，带着这个疑问，我想在这里跟大家交流一下，到底是总能量还是平均每个原子的能量，希望有经验的虫友勇于发言啊~ 这是我计算结果曲线，虽然都计算到了800eV，但总能量还没有达到0.001eV。  [ Last edited by fdd096030079 on 2012-12-5 at 09:41 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 科研 |
» 猜你喜欢
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值
已经有7人回复
-
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有189人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- VASP计算离子步数设置太低，未能收敛，想继续算遇到的问题！
已经有8人回复
- DFPT超胞模型vasp算声子，超胞250个原子，在计算中心12个节点3天3夜还没算完
已经有8人回复
- 关于VASP中对ENCUT收敛性测试的疑问
已经有5人回复
- 测试截断能，k点，和结构优化，哪个在前？
已经有15人回复
- vasp阶段能测试求助，谢谢啊！急求。
已经有15人回复
- 对于收敛性检测
已经有24人回复
- 关于VASP计算收敛的判断
已经有25人回复
- 建模之后进行收敛性测试 主要是看哪些图行能看出是收敛的
已经有22人回复
- 关于截断能和K点测试的问题
已经有6人回复
- VASP静态自洽计算收敛特慢！
已经有31人回复
- VASP中截断能的选取对计算速度影响有多大？
已经有8人回复
- VASP计算表面吸附 但是能量总是不收敛/。。。
已经有9人回复
- k点和截断能的测试问题
已经有14人回复
- vasp计算：优化能收敛，但同样的设置静态计算不收敛
已经有3人回复
- 【求助成功】vasp测试K点收敛性的脚本问题
已经有20人回复
- 【原创】写了个集群下应用测试K点或截断能的脚本
已经有18人回复
- 【求助】vasp截断能不收敛
已经有7人回复
- 【求助】请问收敛测试的标准
已经有15人回复
- 【求助】收敛性测试的问题
已经有40人回复
- 【求助】思考：关于截断能
已经有11人回复
ifmc1234
木虫 (著名写手)
- 应助: 18 (小学生)
- 金币: 2797.9
- 散金: 6001
- 红花: 14
- 帖子: 1761
- 在线: 572.6小时
- 虫号: 825500
- 注册: 2009-08-10
- 专业: 理论和计算化学
2楼2012-12-05 15:36:19
10楼2015-09-02 16:33:14
fdd096030079
新虫 (小有名气)
- 应助: 8 (幼儿园)
- 金币: 481.4
- 散金: 20
- 红花: 5
- 帖子: 112
- 在线: 63小时
- 虫号: 1973034
- 注册: 2012-09-04
- 专业: 凝聚态物性 II ：电子结构
3楼2012-12-05 22:12:03
ifmc1234
木虫 (著名写手)
- 应助: 18 (小学生)
- 金币: 2797.9
- 散金: 6001
- 红花: 14
- 帖子: 1761
- 在线: 572.6小时
- 虫号: 825500
- 注册: 2009-08-10
- 专业: 理论和计算化学
4楼2012-12-06 21:16:09
fdd096030079
新虫 (小有名气)
- 应助: 8 (幼儿园)
- 金币: 481.4
- 散金: 20
- 红花: 5
- 帖子: 112
- 在线: 63小时
- 虫号: 1973034
- 注册: 2012-09-04
- 专业: 凝聚态物性 II ：电子结构
5楼2012-12-09 07:42:31
lrjhxh
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 3654.9
- 红花: 2
- 帖子: 267
- 在线: 371.8小时
- 虫号: 102835
- 注册: 2005-11-14
- 专业: 凝聚态物性 II ：电子结构
6楼2013-05-03 10:48:21
对马拉二胡
至尊木虫 (职业作家)
- 应助: 27 (小学生)
- 金币: 25915.6
- 散金: 244
- 红花: 13
- 帖子: 4065
- 在线: 938小时
- 虫号: 2471734
- 注册: 2013-05-19
- 性别: GG
- 专业: 无机非金属材料
7楼2014-03-15 09:04:51
fdd096030079
新虫 (小有名气)
- 应助: 8 (幼儿园)
- 金币: 481.4
- 散金: 20
- 红花: 5
- 帖子: 112
- 在线: 63小时
- 虫号: 1973034
- 注册: 2012-09-04
- 专业: 凝聚态物性 II ：电子结构
|
仅供参考哦： #PBS -N bulk-encut #PBS -j oe #PBS -l nodes=1:ppn=16 #PBS -l walltime=120:00:00 #PBS -q low cd $PBS_O_WORKDIR NP=`cat $PBS_NODEFILE|wc -l` for en in 300 350 400 450 500 550 600 ;do cat > INCAR <<! Startparameter for this run: # PREC = normal normal or accurate (medium, high low for compatibility) NPAR = 4 ISTART = 0 job : 0-new 1-cont 2-samecut ICHARG = 2 charge: 1-file 2-atom 10-const # ISPIN = 1 spin polarized calculation? # NGX = 26 # NGY = 26 # NGZ = 102 # NBANDS = 150 Electronic Relaxation 1 ENCUT = $en eV NELM = 60; NELMIN= 2; NELMDL=-12 # of ELM steps EDIFF = 0.1E-04 stopping-criterion for ELM # LREAL = F real-space projection Ionic relaxation EDIFFG = -.1E-01 stopping-criterion for IOM NSW = 0 number of steps for IOM IBRION = -1 ionic relax: 0-MD 1-quasi-New 2-CG # ISIF = 3 stress and relaxation # ISYM = 2 0-nonsym 1-usesym 2-fastsym # POTIM = 0.5000 time-step for ionic-motion DOS related values: # EMIN = 10.00; EMAX =-10.00 energy-range for DOS ISMEAR = 0; # SIGMA = 0.20 broadening in eV -4-tet -1-fermi 0-gaus Electronic relaxation 2 (details) IALGO = 48 algorithm Write flags LWAVE = F write WAVECAR LCHARG = F write CHGCAR ! ~/bin/mpirun -np $NP ~/bin/vasp > vasp.out cp CONTCAR POSCAR E=`tail -1 OSZICAR | awk '{printf "%12.6f \n", $5}'` echo $en $E >> en.dat done |
8楼2014-11-13 09:25:28
9楼2015-08-03 10:40:00
20