| 查看: 1845 | 回复: 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从结构优化到光学性质 |
» 猜你喜欢
车间生产，真空度很高，温度很高，但减压蒸馏速度很慢。
已经有32人回复
-
为什么我的博后基金还在流动站审核中？不会是学院给我卡了吧？
已经有22人回复
-
BO4的YQ答辩通知发布了吗？
已经有4人回复
-
Nature 11日发文，中国著名学者们称造假迫不得已
已经有3人回复
-
博后基金，以往的结果点不开，怎么回事呢？最后一次机会了，两次都没中前面。
已经有10人回复
-
请问事业编制和年薪制冲突吗？
已经有4人回复
-
累并快乐着
已经有8人回复
-
审稿问题：为什么荧光激发波长和紫外吸收波长差的大？
已经有7人回复
-
2024国社科通讯评审
已经有7人回复
-
34岁读博士晚吗
已经有22人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 结构优化参数设置
已经有6人回复
- 求回归方程参数估计的最优化算法matlab代码
已经有9人回复
- LC-MS 优化方法时用到的一些参数DP, EP, FP, CXP等分别起什么作用？
已经有7人回复
- 结构优化参数设置
已经有18人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
- 【求助】castep做几何优化时其参数如何设置？
已经有19人回复
- 【求助】vasp优化晶格参数a
已经有13人回复
- 【求助】关于复杂的结构参数的优化
已经有5人回复
- 【求助】Castep中结构优化后晶胞参数差别较大 vasp高手也帮我看一下
已经有22人回复
- 【求助】DMol3优化参数的选择
已经有17人回复
【答案】应助回帖
★ ★ ★
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 (小学生)
- 金币: 1431.8
- 散金: 1197
- 红花: 5
- 帖子: 636
- 在线: 186.9小时
- 虫号: 883307
- 注册: 2009-10-25
- 性别: GG
- 专业: 理论和计算化学
3楼2011-05-26 15:17:14
4楼2011-05-26 15:28:44
tider
木虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1950.4
- 散金: 40
- 红花: 5
- 帖子: 313
- 在线: 87.2小时
- 虫号: 1190614
- 注册: 2011-01-16
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2011-05-26 18:47:07
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值。 |
» 本帖已获得的红花（最新10朵）
洋洋19846楼2011-05-26 20:21:07
nicolemikeyi
禁虫 (著名写手)
navyzhang(金币+10): 谢谢 2011-06-02 15:58:59
|
本帖内容被屏蔽 |
7楼2011-05-26 20:28:24
navyzhang
木虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 1431.8
- 散金: 1197
- 红花: 5
- 帖子: 636
- 在线: 186.9小时
- 虫号: 883307
- 注册: 2009-10-25
- 性别: GG
- 专业: 理论和计算化学
8楼2011-06-02 15:58:41
liliangfang(金币+1): 多谢交流 2012-02-26 08:52:35
|
本帖内容被屏蔽 |
9楼2012-02-26 00:16:07
洋洋1984
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1791.1
- 散金: 143
- 红花: 2
- 帖子: 474
- 在线: 113.1小时
- 虫号: 1303656
- 注册: 2011-05-23
- 性别: GG
- 专业: 数值预报与数值模拟
10楼2012-05-09 11:24:59
5