| 查看: 2869 | 回复: 13 | ||||
| 当前主题已经存档。 | ||||
pingguo802929木虫 (初入文坛)
|
[交流]
【讨论】讨论VASP里边ISYM参数的重要性！
|
|||
|
各位虫友： 据我的理解，对于VASP里边的ISYM参数，如果不对它进行设置， 让它保持默认值1或2，那么优化计算时，自始至终体系的初始 对称性将保持不变。 如果把它设置为0，则计算时将不考虑对称性， 优化出来的结构和初始结构的对称性有可能不一样。 我是用PSTRESS参数给体系做加压的优化计算，我做了一些测试，感觉 ISYM对优化的结果（如能量，构型）影响并不是很大，大家对ISYM参数 有没有什么其他看法？加压优化计算时，需要把ISYM设置为0吗？ [ Last edited by aylayl08 on 2010-4-21 at 08:59 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP |
» 猜你喜欢
085600，321分求调剂
已经有11人回复
-
复试调剂
已经有8人回复
-
材料调剂
已经有12人回复
-
301求调剂
已经有11人回复
-
283求调剂
已经有10人回复
-
0703化学
已经有11人回复
-
322求调剂：一志愿湖南大学 材料与化工（085600），已过六级。
已经有18人回复
-
070300化学学硕311分求调剂
已经有11人回复
-
材料化工306分找合适调剂
已经有12人回复
-
298求调剂
已经有3人回复
y1ding
铁杆木虫 (著名写手)
- 1ST强帖: 1
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
2楼2009-05-07 15:46:21
y1ding
铁杆木虫 (著名写手)
- 1ST强帖: 1
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
3楼2009-05-07 15:49:53
y1ding
铁杆木虫 (著名写手)
- 1ST强帖: 1
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
★ ★
spur(金币+2,VIP+0):辛苦！^_^ 5-8 13:38
spur(金币+2,VIP+0):辛苦！^_^ 5-8 13:38
|
9.6 Volume vs. energy, volume relaxations, Pulay Stress If you are doing energy-volume calculations or cell shape and volume relaxations you must understand the Pulay stress, and related problems. The Pulay stress arises from the fact that the plane wave basis set is not complete with respect to changes of the volume. Thus, unless absolute convergence with respect to the basis set has been achieved - the diagonal components of the stress tensor are incorrect. This error is often called ``Pulay stress''. The error is almost isotropic (i.e. the same for each diagonal component), and for a finite basis set it tends to decrease volume compared to fully converged calculations (or calculations with a constant energy cutoff). The Pulay stress and related problems affect the behavior of VASP and any plane wave code in several ways: First it evidently affects the stress tensor calculated by VASP, i.e. the diagonal components of the stress tensor are incorrect, unless the energy cutoff is very large (ENMAX=1.3 *default is usually a save setting to obtain a reliable stress tensor). In addition it should be noted that all volume/cell shape relaxation algorithms implemented in VASP work with a constant basis set. In that way all energy changes are strictly consistent with the calculated stress tensor, and this in turn results in an underestimation of the equilibrium volume unless a large plane wave cutoff is used. Keeping the basis set constant during relaxations has also some strange effect on the basis set. Initially all G-vectors within a sphere are included in the basis. If the cell shape relaxation starts the direct and reciprocal lattice vectors change. This means that although the number of reciprocal G-vectors in the basis is kept fixed, the length of the G-vectors changes, changing indirectly the energy cutoff. Or to be more precise the shape of cutoff region becomes an elipsoide. Restarting VASP after a volume relaxation causes VASP to adopt a new ``spherical'' cutoff sphere and thus the energy changes discontinuously (see section 7.10). One thing which is important to understand, is that problems due to the Pulay stress can often be neglected if only volume conserving relaxations are performed. This is because the Pulay stress is usually almost uniform and it therefore changes the diagonal elements of the stress tensor only by a certain constant amount (see below). In addition many calculations have shown that Pulay stress related problems can also be reduced by performing calculations at different volumes using the same energy cutoff for each calculation (this is what VASP does per default, see section 7.10), and fitting the final energies to an equation of state. This of course implies that the number of basis vectors is different at each volume. But calculations with many plane wave codes have shown that such calculations give very reliable results for the lattice constant and the bulk modulus and other elastic properties even at relatively small energy cutoffs. In a certain way constant energy cut-off are thus less prone to errors cause by the basis set incompleteness than constant basis set calculations. But it should be kept in mind that volume changes and cell shape changes must be rather large in order to obtain reliable results from this method, because in the limit of very small distortions the energy changes obtained with this method are equivalent with that obtained from the stress tensor and are therefore affected by the Pulay stress. Only large volume changes guarantee that the errors introduced by the basis set incompleteness are averaged out. |
4楼2009-05-07 15:50:09
pingguo802929
木虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 2185.3
- 散金: 4
- 帖子: 34
- 在线: 9.3小时
- 虫号: 708054
- 注册: 2009-02-25
- 性别: MM
- 专业: 会计与审计
5楼2009-05-07 16:35:24
pingguo802929
木虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 2185.3
- 散金: 4
- 帖子: 34
- 在线: 9.3小时
- 虫号: 708054
- 注册: 2009-02-25
- 性别: MM
- 专业: 会计与审计
6楼2009-05-07 16:36:48
y1ding
铁杆木虫 (著名写手)
- 1ST强帖: 1
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
7楼2009-05-08 09:13:37
hchfox
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3216.5
- 红花: 1
- 帖子: 609
- 在线: 101小时
- 虫号: 280207
- 注册: 2006-09-17
- 专业: 其他无机非金属材料
8楼2009-05-08 09:45:52
spur
金虫 (正式写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 贵宾: 1.21
- 金币: 4994.7
- 散金: 559
- 红花: 5
- 帖子: 866
- 在线: 374.9小时
- 虫号: 215683
- 注册: 2006-03-11
- 性别: GG
- 专业: 第一性原理计算
9楼2009-05-08 13:36:46
hchfox
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3216.5
- 红花: 1
- 帖子: 609
- 在线: 101小时
- 虫号: 280207
- 注册: 2006-09-17
- 专业: 其他无机非金属材料
10楼2009-05-08 15:36:27
