| 查看: 2483 | 回复: 8 | |||
qphll金虫 (正式写手)
|
[交流]
【求助】关于结构参数的优化
|
|
一直有个问题不是很清楚,请有经验的大家解答一下. 在做结构参数优化的时候,如果需要优化的参数不止一个,那么应该怎样? 比如,如果需要对Orthorhombic(正交)盒子进行优化, 那么就会涉及到需要优化的参数有三个. 如果是在pwscf中,那么就是a, b/a和c/a. 对于我而已,我是想研究小分子在表面的分解, 那么一般c取到一定大以后,就不去考虑c/a的优化了. 那么在实际的操作中怎么分别优化a和b/a? 我的理解和做法是这样的: 分两步来 (1) 固定b/a的值,看a变化的过程中,体系总能量的变化,从而找到a的优化值; (2) 固定a值为(1)中找到的最优值,然后调节b/a的值, (事实上这里就是等同于改变b的值), 然后找到b的最优值. 请大家评价一下,是不是就是这样做的? 还有,如果我是这样两步做的: (1) 固定b的值,看a变化的过程中,体系总能量的变化,从而找到a的优化值; (2) 固定a值为(1)中找到的最优值,然后改变b的值, 然后找到b的最优值. 如果这样做,那么和前面的那种优化方法有啥本质的区别吗? 搬个板凳,带着笔记,等大家拍砖. |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
fileCode有新解读？
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- MS 结构优化的结果不能理解 求高人指点
已经有10人回复
- 关于VASP一些琐事的求助（K点，吸附优化等）
已经有12人回复
- 结构优化参数设置
已经有6人回复
- 关于建模结构优化的一些问题
已经有11人回复
- 【活动or资源】vasp晶体结构优化
已经有142人回复
- 请大家帮帮忙，看看这CASTEP结构优化怎么失败了
已经有11人回复
- 【原创】关于结构优化的总结，体会
已经有33人回复
- 【求助】结构优化中，如何选择POTIM？
已经有7人回复
- 【求助】VASP怎么做晶格参数的结构优化
已经有17人回复
- 【求助】关于优化时对称性是否加入的问题 能量的选择问题 赝势的选择问题
已经有5人回复
- 【求助】关于优化结果
已经有17人回复
- 【求助】简单问题----关于vasp晶格优化
已经有7人回复
- 【求助】关于复杂的结构参数的优化
已经有5人回复
- 【求助】siesta中关于几何结构优化的问题【完结】
已经有17人回复
- 【求助】DMol3优化参数的选择
已经有17人回复
zxzj05
荣誉版主 (著名写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 贵宾: 1.242
- 金币: 5912.5
- 红花: 2
- 帖子: 1528
- 在线: 144.5小时
- 虫号: 406147
- 注册: 2007-06-18
- 性别: GG
- 专业: 金属功能材料
- 管辖: 第一性原理
2楼2009-09-23 09:32:08
qphll
金虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
3楼2009-09-23 09:37:32
y1ding
铁杆木虫 (著名写手)
- 1ST强帖: 1
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
4楼2009-09-23 10:52:22
xirainbow
木虫 (正式写手)
- 1ST强帖: 2
- 应助: 0 (幼儿园)
- 贵宾: 0.03
- 金币: 2132.6
- 红花: 2
- 帖子: 775
- 在线: 6.1小时
- 虫号: 431211
- 注册: 2007-08-11
- 专业: 光谱学
★ ★ ★ ★
wuchenwf(金币+2,VIP+0):谢谢 9-24 22:00
qphll(金币+2,VIP+0): 9-25 11:42
wuchenwf(金币+2,VIP+0):谢谢 9-24 22:00
qphll(金币+2,VIP+0): 9-25 11:42
|
PWscf的参数“cell_dofree”专门处理你提出的问题，只优化a,b而固定c 但是这个方法bfgs好像是暂时还没有实现，你需使用共轭梯度法 你可以看看这个：“http://blog.sina.com.cn/s/blog_5f15ead20100e1lu.html###”（Some fundamental concept and parameters of pwscf） 里面有关于此问题PWscf forum的网址 |
5楼2009-09-24 09:47:17
qphll
金虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
|
多谢楼上三位的参与. 对于我自己的特定体系而已, 我发帖的里面提到的结构参数的优化问题已经想明白了. 对于graphene而言,因为我自己编程生成graphene的原因,我的盒子是Orthorhombic的,a,b,c三个值各不相等,alpha, beta, gamma 都是90度. 如果想想结构优化的目的,我想对于我的体系而已,无非是因为在生成graphene的时候, C-C的键长使用了实验值1.42 那么这里结构优化的目的,就是通过盒子的'扩张'或者'收缩', 找到一个合适的C-C键长而已. 如果keep that in mind, 那么虽然我这里a和b的数值是不相等的,但是其实它们也不是互相独立的.如果我将a (X轴)拉伸,那么为了维持graphene中C-C键长的一致性, b(Y轴)必定得要按比例拉伸. 这里c (Z轴)在取到一定值以后,保证graphene面和它的Z轴镜像没有相互左右以后,就不需要优化啦. 所以我想,对于我的这个长方体的模拟盒子,放在pwscf里面考虑lattice parameter优化的话,是可以这么一步做的: (在in文件中,有celldm(1)=a, celldm(2)=b/a,celldm(3)=c/a 这三个参数) 固定b/a的值,固定c的值, 看a变化的过程中,体系总能量的变化,从而找到a的优化值 (事实上b/a固定的话,那么a的最优值找到了,那么b的值也相应找到了). 超级新手, pw里面的vc-relax没有尝试,看manual也似乎没说啥,例子也没有找到. y1ding可否提供一个in文件给我学习学习? |
6楼2009-09-25 12:02:36
7楼2009-09-28 15:30:47
zxzj05
荣誉版主 (著名写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 贵宾: 1.242
- 金币: 5912.5
- 红花: 2
- 帖子: 1528
- 在线: 144.5小时
- 虫号: 406147
- 注册: 2007-06-18
- 性别: GG
- 专业: 金属功能材料
- 管辖: 第一性原理
8楼2009-09-28 17:20:17
改王之王
金虫 (正式写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 金币: 1120.7
- 散金: 61
- 红花: 3
- 帖子: 817
- 在线: 113.5小时
- 虫号: 968169
- 注册: 2010-03-11
- 专业: 凝聚态物性 II ：电子结构
9楼2011-03-31 12:35:32
100