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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】如何从算出的E-a数值拟和出晶格常数a0
191/2返回列表12下一页
查看: 979  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上和省基金B类撒花 +6 花田半亩～白 2026-04-21 6/300 2026-04-22 10:55 by 3126142009
[考博] 华师大读博 +3 xq83 2026-04-22 5/250 2026-04-22 10:42 by xq83
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +5 yexuqing 2026-04-19 5/250 2026-04-22 10:01 by easeheart
[考研] 0854求调剂 +24 门路摸摸 2026-04-15 28/1400 2026-04-22 09:54 by Sy199704！
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 8/400 2026-04-21 16:26 by 啃面包的小书虫
[考研] 295分求调剂 +6 ?要上岸? 2026-04-17 6/300 2026-04-21 08:18 by Equinoxhua
[考研] 一志愿中科大材料与化工，353分还有调剂学校吗 +11 否极泰来2026 2026-04-15 13/650 2026-04-20 22:31 by Equinoxhua
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[论文投稿] 期刊推荐 +3 材料研究生 2026-04-15 5/250 2026-04-20 16:02 by 豆豆7758
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
信息提示
关闭
请填处理意见
关闭