| 查看: 2509 | 回复: 11 | ||||||
| 本帖产生 2 个 1ST强帖 ，点击这里进行查看 | ||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||||
navyzhang木虫 (正式写手)
|
[求助]
最复杂情况下的参数优化
|
|||||
|
看了侯老师的VASP教程，关于晶格参数优化存在一些疑问！ 对于a,b,c都不相等的体系，是不是可以这样 1、先优化形状(即c/a,b/a)和位置，(将该计算得到的CONTCAR作为下一步的POSCAR) 2、仅优化位置，得到该情况下的能量E1 3、改变a参数，按上面两步，先优化形状和位置，再单优化位置，得到该情况下的能量E2 4、如此往复共改变十次a参数，得到是个能量E1~E10，拟合E-V曲线，得到能量最低时的V0 那么下面参数a到底该怎么求呢？将V0/(b*c)？如果是这样，b和c从哪里来？是能量最低情况下第一步对应的结果？ |
回复此楼» 收录本帖的淘帖专辑推荐
 第一原理资料汇编 | 计算-vasp | 分子模拟专题 | vasp从结构优化到光学性质 |
» 猜你喜欢
0856材料专硕353求调剂
已经有8人回复
-
材料与化工085600，总分304，本科有两篇sci参与，求调剂
已经有6人回复
-
南昌大学材料专硕311分求调剂
已经有9人回复
-
材料与化工考研调剂
已经有12人回复
-
340求调剂
已经有3人回复
-
安徽大学专硕生物与医药专业(086000)324分，英语已过四六级，六级521，求调剂
已经有3人回复
-
292求调剂
已经有12人回复
-
085600材料与化工306
已经有10人回复
-
359求调剂
已经有4人回复
-
求调剂
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 结构优化参数设置
已经有6人回复
- 求回归方程参数估计的最优化算法matlab代码
已经有9人回复
- LC-MS 优化方法时用到的一些参数DP, EP, FP, CXP等分别起什么作用？
已经有7人回复
- 结构优化参数设置
已经有18人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
- 【求助】castep做几何优化时其参数如何设置？
已经有19人回复
- 【求助】vasp优化晶格参数a
已经有13人回复
- 【求助】关于复杂的结构参数的优化
已经有5人回复
- 【求助】Castep中结构优化后晶胞参数差别较大 vasp高手也帮我看一下
已经有22人回复
- 【求助】DMol3优化参数的选择
已经有17人回复
nicolemikeyi
禁虫 (著名写手)
- 1ST强帖: 1
- 应助: 6 (幼儿园)
- 金币: 5725.1
- 散金: 94
- 红花: 15
- 帖子: 2539
- 在线: 616.8小时
- 虫号: 1241122
- 注册: 2011-03-22
- 性别: GG
- 专业: 固体力学
【答案】应助回帖
★ ★ ★
gzqdyouxia(金币+3, 1ST强帖+1): 谢谢详细指教 2011-05-27 10:43:47
gzqdyouxia(金币+3, 1ST强帖+1): 谢谢详细指教 2011-05-27 10:43:47
|
那里头的脚本分为两部分，一部分是固定晶胞体积的弛豫，得到CONTCAR后复制到POSCAR中进行另一部分静态计算，得出能量。所以侯的优化思路是这样的： 1.取不同的a值，首先进行ISIF=5的固定体积但不固定形状的弛豫，然后静态计算出对应的能量。 2.每一个a值对应一个体积和能量，多个a值可以得到一条E~V曲线，通过拟合可以得到最小的那个E对应的体积和a值。 3.注意，此时如果是立方晶格，a值即为优化结果。但例如是四方或其他，此时你只是得到了最佳的形状（a/c值），你再取不同的a进行下一次循环，以得到最佳的a值与c值。 |
洋洋19846楼2011-05-26 20:21:07
【答案】应助回帖
★ ★ ★
youzhizhe(金币+3): 谢谢交流。 2011-05-26 12:27:49
youzhizhe(1ST强帖+1): 谢谢交流。 2011-05-26 16:41:46
youzhizhe(金币+3): 谢谢交流。 2011-05-26 12:27:49
youzhizhe(1ST强帖+1): 谢谢交流。 2011-05-26 16:41:46
|
首先，结构优化是为了寻找晶体的基态。 在侯博士的那片文章中提到的对复杂晶体的优化，再脚本中的a首先是在实验值的左右各取若干个，然后对晶格采用固定体积的方法进行优化，这样就会得到一系列的a-E图像，做出曲线，这样的话，你就能看到这个曲线呈现U型。我们知道基态的能量应该是最低的，所以a的取值应该是对应于E最小的a，但是一次并不能精确的确定a的值，只能错略的估计a的范围，但是通过多次的缩小这个区间，就能得到较为准确的a值，而由于是允许晶格形可变，且原子位置驰豫的，所以就可以得到基态结构。 |
2楼2011-05-26 11:37:46
navyzhang
木虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 1411.8
- 散金: 1197
- 红花: 5
- 帖子: 640
- 在线: 187.1小时
- 虫号: 883307
- 注册: 2009-10-25
- 性别: GG
- 专业: 理论和计算化学
3楼2011-05-26 15:17:14
4楼2011-05-26 15:28:44
