应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】关于结构参数的优化
91/1返回列表
查看: 2002  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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的最优值.

如果这样做,那么和前面的那种优化方法有啥本质的区别吗?

搬个板凳,带着笔记,等大家拍砖.
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
Life, Love, Laugh.
1楼 2009-09-23 06:22:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxzj05

荣誉版主 (著名写手)
★ ★
qphll(金币+1,VIP+0):多谢参与. 9-23 09:38
wuchenwf(金币+1,VIP+0):谢谢 9-24 21:59
"固定b/a的值,看a变化的过程中"变化的是b和a，后一种方法需要更多次的循环才能达到第一种的效果。
赞一下(12人) 回复此楼
储氢家族欢迎储氢研究者！
2楼2009-09-23 09:32:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)
zxzj05的回复我有些不明白.

你是觉得这样做是合适的:
(1) 固定b/a的值,看a变化的过程中,体系总能量的变化,从而找到a的优化值;
(2) 固定a值为(1)中找到的最优值,然后调节b/a的值, (事实上这里就是等同于改变b的值), 然后找到b的最优值.
赞一下 回复此楼
Life, Love, Laugh.
3楼2009-09-23 09:37:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

y1ding

铁杆木虫 (著名写手)
★
freshgirl(金币+1,VIP+0):谢谢参与~ 9-23 14:10
pw的话,直接 vc-relax
赞一下(10人) 回复此楼
4楼2009-09-23 10:52:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)
★ ★ ★ ★
wuchenwf(金币+2,VIP+0):谢谢 9-24 22:00
qphll(金币+2,VIP+0): 9-25 11:42
引用回帖:
Originally posted by qphll at 2009-9-23 06:22:
一直有个问题不是很清楚,请有经验的大家解答一下.

在做结构参数优化的时候,如果需要优化的参数不止一个,那么应该怎样?

比如,如果需要对Orthorhombic(正交)盒子进行优化, 那么就会涉及到需要优化的参数有三个 ...

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的网址
赞一下(11人) 回复此楼
5楼2009-09-24 09:47:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)
多谢楼上三位的参与. 对于我自己的特定体系而已, 我发帖的里面提到的结构参数的优化问题已经想明白了.

对于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文件给我学习学习?
赞一下 回复此楼
Life, Love, Laugh.
6楼2009-09-25 12:02:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

apzhough

铁虫 (小有名气)
路过，踩踩
回复此楼
7楼2009-09-28 15:30:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxzj05

荣誉版主 (著名写手)
引用回帖:
Originally posted by qphll at 2009-9-23 09:37:
zxzj05的回复我有些不明白.

你是觉得这样做是合适的:
(1) 固定b/a的值,看a变化的过程中,体系总能量的变化,从而找到a的优化值;
(2) 固定a值为(1)中找到的最优值,然后调节b/a的值, (事实上这里就是等同于改变b ...

我是解释，问哪个合适，显然是后者
赞一下 回复此楼
储氢家族欢迎储氢研究者！
8楼2009-09-28 17:20:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)
引用回帖:
Originally posted by zxzj05 at 2009-09-28 17:20:17:
我是解释，问哪个合适，显然是后者

这是为什么，不都是改变b/a的比例吗？我现在也在尝试用这种方法优化…………
赞一下 回复此楼
9楼2011-03-31 12:35:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qphll 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭