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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】如何从算出的E-a数值拟和出晶格常数a0
191/2返回列表12下一页
查看: 856  |  回复: 18
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

yq2240711

金虫 (正式写手)

[交流] 【求助】如何从算出的E-a数值拟和出晶格常数a0

如何从算出的一系列E-a数值拟和出晶格常数a0
如何利用物态方程拟和晶格常数？

[ Last edited by yq2240711 on 2009-11-21 at 21:19 ]
回复此楼

» 猜你喜欢
1楼2009-11-11 19:19:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)
★ ★ ★ ★
yq2240711(金币+3,VIP+0): 11-11 19:35
yjmaxpayne(金币+1,VIP+0):谢谢解答 11-11 20:01
用pwscf可以直接拟合，参看：
http://blog.sina.com.cn/s/blog_5f15ead20100f5jx.html

可以直接用泰勒展开去拟合么？
请高人帮忙；）
赞一下(7人) 回复此楼
2楼2009-11-11 19:34:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdhlover

荣誉版主 (职业作家)

小木虫之莲花佛
★ ★ ★ ★ ★ ★
yq2240711(金币+1,VIP+0): 11-11 20:40
yjmaxpayne(金币+5,VIP+0):非常感谢提供信息 11-11 21:17
On the other hand, we may use a relatively simple approach to refine the determined equilibrium
value of the lattice constant by fitting our sampled data points to a second order
polynomial of the form

P(a)=A(a-a0)*(a-a0)+B(a-a0)+c
We may then choose the minimum of P(a)
amin=a0-B/(2A)
as the value of the lattice constant . We now set a0=2.999 (the central point of the interval) and let h= 0.01 Rydberg (the equidistant spacing between sample points). In addition, using a few algebraic manipulations leads to the following analytic expressions for the two constants A and B
:

A=(E(a0+h)+E(a0-h)-2E(a0))/(2h*2h)
and
B=(E(a0+h)-E(a0-h)))/2h

then you can get the lattice constant :amin=a0-B/(2A)
赞一下(6人) 回复此楼
3楼2009-11-11 20:06:57
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)
引用回帖:
Originally posted by zdhlover at 2009-11-11 20:06:
On the other hand, we may use a relatively simple approach to refine the determined equilibrium
value of the lattice constant by fitting our sampled data points to a second order
polynomial of th ...

您提供的信息太好了：）
请问，您从哪里获得的这些说明 ？
能给个链接么？
谢谢啦；）
赞一下 回复此楼
4楼2009-11-11 20:54:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdhlover

荣誉版主 (职业作家)

小木虫之莲花佛
★
yq2240711(金币+1,VIP+0): 12-17 18:35
引用回帖:
Originally posted by xirainbow at 2009-11-11 00:54:

您提供的信息太好了：）
请问，您从哪里获得的这些说明 ？
能给个链接么？
谢谢啦；）

这个就是二项式拟合，任何一个教材上都有的
赞一下 回复此楼
5楼2009-11-11 21:49:02
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdd2008

金虫 (职业作家)
★ ★
yjmaxpayne(金币+2,VIP+0):谢谢参与讨论 11-12 19:12
建议拟合时不要使用普通的多项式进行拟合！

当然如果你只为了得到平衡体积（平衡晶格常数）和对应的能量，你可以使用多项式拟合！

但使用多项式拟合没有物理意义,因此，建议使用物态方程拟合！这样拟合不仅可以得到平衡体积（平衡晶格常数）和对应的能量，以及体模量，及其一阶和二阶导数！呵呵！！
赞一下(6人) 回复此楼
6楼2009-11-11 22:15:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)
引用回帖:
Originally posted by zdd2008 at 2009-11-11 22:15:
建议拟合时不要使用普通的多项式进行拟合！

当然如果你只为了得到平衡体积（平衡晶格常数）和对应的能量，你可以使用多项式拟合！

但使用多项式拟合没有物理意义,因此，建议使用物态方程拟合！这样拟合不仅 ...

请问，如何使用物态方程拟合正交，单斜，三斜晶系的情况呢：）？

[ Last edited by xirainbow on 2009-11-11 at 22:38 ]
赞一下 回复此楼
7楼2009-11-11 22:35:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdd2008

金虫 (职业作家)
★ ★ ★ ★
yjmaxpayne(金币+4,VIP+0):谢谢分享 11-12 19:12
引用回帖:
Originally posted by xirainbow at 2009-11-11 22:35:


请问，如何使用物态方程拟合正交，单斜，三斜晶系的情况呢：）？

[ Last edited by xirainbow on 2009-11-11 at 22:38 ]

物态方程拟合的是体积与能量，体积与压力之间的关系！因此对于正交，单斜，三斜晶系也适用！！呵呵！！
我明白你的意思，比如拿正交来讲，拟合得到平衡体积后，根本不知道a，b，c分别是多少！此时可以采用全弛豫的a:b:c来确定！这是我的处理方法！呵呵！！不知道讲明白了没有！！
赞一下(6人) 回复此楼
8楼2009-11-11 23:49:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdd2008

金虫 (职业作家)
★ ★ ★ ★ ★
yjmaxpayne(金币+5,VIP+0):非常感谢细致的讲解 11-12 19:12
我做的一个图！！呵呵！
赞一下(6人) 回复此楼
9楼2009-11-12 13:10:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)
★ ★
fegg7502(金币+2,VIP+0):thank you very much! 11-15 17:23
引用回帖:
Originally posted by zdd2008 at 2009-11-11 23:49:



物态方程拟合的是体积与能量，体积与压力之间的关系！因此对于正交，单斜，三斜晶系也适用！！呵呵！！
我明白你的意思，比如拿正交来讲，拟合得到平衡体积后，根本不知道a，b，c分别是多少！此时可以采用 ...

非常感谢您的指点；）还给出了一张这么漂亮的图；）
您说的正是我想问的；）
对于三斜，不仅abc不定，连他们之间的夹角也不定；（

您所说的全弛豫是什么意思呢？
请您再说明一下；）
我用的软件无法优化晶胞，所以才面临这样的问题；）

[ Last edited by xirainbow on 2009-11-12 at 19:43 ]
赞一下(10人) 回复此楼
10楼2009-11-12 19:41:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yq2240711 的主题更新
191/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭