| 查看: 2234 | 回复: 8 | ||||
剑魄诗魂木虫 (正式写手)
|
[求助]
求教表面吸附的频率计算问题
|
|
最近在尝试计算表面吸附结构的自由能变化，所以需要计算频率以得到ZPE。但是在实际的计算中发现，计算时间非常久（超过4天），而且CHG, CHGCAR, WAVECAR变得非常大（合计超过50G），所以不得不停止了计算。 在此想跟同行求教一下是否我的具体设置上存在什么问题所以导致不容易计算？我的体系是在金属团簇表面的吸附，具体策略是固定所有金属原子（坐标文件如下）而不固定吸附物。 POSCAR： 1.00000000000000 40.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 40.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 40.0000000000000000 M C O 79 1 2 Direct 0.4993046399203109 0.4992548039135158 0.4992924124318137 F F F 0.4990485543665031 0.5438863108801399 0.5435908768283824 F F F 0.4995033566807336 0.4546647224233595 0.4550508152545660 F F F .. (为节约空间省略若干行） 0.5877448409594663 0.4996622676310523 0.4992604420477491 F F F 0.4996940843718747 0.4109563411793573 0.4996687239395428 F F F 0.4989078996841698 0.5876795921551263 0.4988993667433306 F F F 0.5944357731254115 0.6390016216759307 0.5953946665121848 T T T 0.5823633152745984 0.6624175346966013 0.5822536135436280 T T T 0.6066869688169945 0.6157805664098532 0.6086636088716872 T T T INCAR如下： ISTART = 1 ICHARG = 1 PREC = Normal ALGO = Normal ISYM = 0 !0-no symmetry, 2-use symmetry NSW = 1 ! 0 ionic steps ENCUT = 400 ! energy cutoff for planewaves, ENMAX*1.3 when high ISPIN = 2 !1-nonspin, 2-spin GGA = PE ! PAW-PW91 LVDW = True ! DFT-D IVDW = 11 ! zero damping DFT-D3 method of Grimme EDIFF = 1E-03 ! Energy convergency EDIFFG = -0.03 LREAL = .True. ISMEAR = 1 ! Metallic SIGMA = 0.2 ! frequency POTIM = 0.02 IBRION = 5 ! frequency calculation NFREE = 2 还望有经验的同行不吝赐教，我这个体系如何设置能够顺利计算表面吸附分子的振动频率。 |
回复此楼» 猜你喜欢
【考研调剂】化学专业 281分，一志愿四川大学，诚心求调剂
已经有11人回复
-
一志愿深大，0703化学，总分302，求调剂
已经有4人回复
-
材料学硕333求调剂
已经有3人回复
-
材料求调剂
已经有3人回复
-
333求调剂
已经有5人回复
-
求助
已经有6人回复
-
278求调剂
已经有9人回复
-
297求调剂
已经有12人回复
-
0703化学297求调剂
已经有3人回复
-
材料学硕301分求调剂
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
漫天飘雪
版主 (知名作家)
认真做事，踏实做人
- 1ST强帖: 1
- 应助: 362 (硕士)
- 贵宾: 0.108
- 金币: 16867.2
- 散金: 13424
- 红花: 392
- 沙发: 2
- 帖子: 5766
- 在线: 2306.6小时
- 虫号: 914852
- 注册: 2009-11-27
- 性别: GG
- 专业: 材料物理化学
- 管辖: 第一性原理
2楼2017-08-04 12:49:14
剑魄诗魂
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 6 (幼儿园)
- 贵宾: 0.06
- 金币: 2685.1
- 散金: 590
- 红花: 7
- 帖子: 405
- 在线: 109.9小时
- 虫号: 1061894
- 注册: 2010-07-21
- 性别: GG
- 专业: 理论和计算化学
3楼2017-08-05 01:20:42
漫天飘雪
版主 (知名作家)
认真做事，踏实做人
- 1ST强帖: 1
- 应助: 362 (硕士)
- 贵宾: 0.108
- 金币: 16867.2
- 散金: 13424
- 红花: 392
- 沙发: 2
- 帖子: 5766
- 在线: 2306.6小时
- 虫号: 914852
- 注册: 2009-11-27
- 性别: GG
- 专业: 材料物理化学
- 管辖: 第一性原理
4楼2017-08-05 09:01:23
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3511.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
漫天飘雪: 金币+2, 谢谢交流 2017-08-06 09:46:55
剑魄诗魂(franch代发): 金币+10, 谢谢回帖交流 2017-08-27 13:57:10
漫天飘雪: 金币+2, 谢谢交流 2017-08-06 09:46:55
剑魄诗魂(franch代发): 金币+10, 谢谢回帖交流 2017-08-27 13:57:10
|
我不知道为啥你的格子会是40x40x40这么大，吸附的话格子很少是这样的整数，莫非你做的是团簇的吸附？如果是团簇的吸附你可以考虑放弃平面波的方案，考虑使用原子轨道，40x40x40的格子效率很低。 如果是slab体系的话，我也不建议轻易去降低精度，可能会导致引起虚频。如果体系过大，硬盘难以容纳的话我建议用phonopy引导vasp用有限差分算，这样每个有限位移（ 不考虑对称性的话是6N个结构，因为每个原子要移动三个方向的正负各一次）产生的结构可以分开到不同的服务器算，最后汇总记录力的文件vasp.xml，然后用phonopy读取做差分即可，缺点就是如果要对原子做固定需要一些技巧，而且需要写脚本。 |
5楼2017-08-05 09:41:03
剑魄诗魂
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 6 (幼儿园)
- 贵宾: 0.06
- 金币: 2685.1
- 散金: 590
- 红花: 7
- 帖子: 405
- 在线: 109.9小时
- 虫号: 1061894
- 注册: 2010-07-21
- 性别: GG
- 专业: 理论和计算化学
6楼2017-08-06 01:19:54
剑魄诗魂
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 6 (幼儿园)
- 贵宾: 0.06
- 金币: 2685.1
- 散金: 590
- 红花: 7
- 帖子: 405
- 在线: 109.9小时
- 虫号: 1061894
- 注册: 2010-07-21
- 性别: GG
- 专业: 理论和计算化学
。非常感谢帮助！ 7楼2017-08-06 01:21:31
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3511.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
8楼2017-08-06 10:36:15
9楼2017-08-25 22:09:27
40