| 查看: 2019 | 回复: 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的最优值. 如果这样做,那么和前面的那种优化方法有啥本质的区别吗? 搬个板凳,带着笔记,等大家拍砖. |
回复此楼» 猜你喜欢
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有7人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有7人回复
-
磺酰氟产物，毕不了业了！
已经有5人回复
-
论文终于录用啦！满足毕业条件了
已经有25人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有3人回复
-
自荐读博
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有5人回复
-
不自信的我
已经有5人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
qphll
金虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
3楼2009-09-23 09:37:32
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