| 查看: 2116 | 回复: 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的最优值. 如果这样做,那么和前面的那种优化方法有啥本质的区别吗? 搬个板凳,带着笔记,等大家拍砖. |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
