24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

qx08131989

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2374
帖子: 196
在线: 155.1小时
虫号: 659563

[交流] 【求助】E-V曲线晶体状态方程拟合问题

各位虫友，我对金属单质的hcp和bcc结构进行了结构优化，获得了最稳定的晶格常数a，b，c，然后逐渐增大和减少晶格常数a，b，c，变化幅度为0.02,进而获得一系列的E-V数据点，然后用这些数据点绘制成曲线，再用三阶BM-EOS的能量-体积关系进行拟合。我的问题如下：
1.这样做获得的E-V曲线是否合理。若不合理，应该如何获得E-V曲线？
2.晶体状态方程的三阶BM-EOS的能量-体积关系是否适合于任何一种结构。
3.状态方程参量的单位，例如，体积模量。

回复此楼

» 猜你喜欢

青椒八年已不青，大家都被折磨成啥样了？已经有11人回复
限项规定已经有9人回复
免疫学博士有名额，速联系已经有4人回复
交叉科学部支持青年基金，对三无青椒是个机会吗？已经有5人回复
国家基金申请书模板内插入图片不可调整大小? 已经有6人回复
国家级人才课题组招收2026年入学博士已经有5人回复
Fe3O4@SiO2合成已经有6人回复
青年基金C终止已经有4人回复
26申博求博导推荐-遥感图像处理方向已经有4人回复
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

四条曲线拟合成一个设定的方程已经有7人回复
origin中如何自定义Birch–Murnaghan三阶拟合已经有40人回复
急切求助EOS拟合已经有21人回复
关于晶体优化问题已经有6人回复
请教如何用VASP计算氧化物的晶格常数？以alpha-Al2O3为例已经有7人回复
求助Freundlich拟合方程那个是Kf值，哪个是n值已经有5人回复
VASP 计算E-V曲线求助。。已经有10人回复
求助用Oringin拟合转状态方程已经有6人回复
大家看看我的excell拟合方程，怎么能合理的解释相关性呢？已经有8人回复
求绳子的曲线方程已经有6人回复
关于药时曲线的合理拟合问题已经有28人回复
如何用VTF方程拟合电导率与温度曲线已经有8人回复
关于方程式创建曲线的问题已经有6人回复
请教下已知曲线以及曲线上的各个散点如何求出这条曲线的方程表达式已经有10人回复
实验数据的数学模型构建、曲线拟合方面问题求助已经有6人回复
【讨论】我收集到的状态方程编程已经有7人回复
【求助】状态方程已经有6人回复
【求助】拟合晶格常数已经有14人回复
【原创】elk中的状态方程拟合工具eos 已经有128人回复
【讨论】vasp结构优化之前的测试已经有11人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-12-31 11:27:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

babaleo

木虫 (正式写手)

应助: 18 (小学生)
金币: 2391.1
帖子: 375
在线: 477.1小时
虫号: 981394

引用回帖:

Originally posted by zhang668 at 2010-12-31 12:54:23:
cname             : name of crystal up to 256 characters
natoms             : number of atoms in unit cell
etype             : equation of state type (see below)
vplt1, vplt2, nvplt : volum ...

" ....hcp有两个自由度a，c。最好保持c/a比例不变来增加体积...". 然后是不是还要在每个固定体积下变c/a在找能量最低点以后才算找到用来拟合的E-V值?

赞一下

回复此楼

4楼2011-01-17 09:41:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 9 个回答

zhang668

木虫 (正式写手)

★ ★ ★
ellsaking(金币+3):谢谢详细说明！ 2010-12-31 13:58:49
qx08131989(金币+3): 2010-12-31 15:16:02
qx08131989(金币+3):再次感谢版主！ 2010-12-31 15:53:30

cname             : name of crystal up to 256 characters
natoms             : number of atoms in unit cell
etype             : equation of state type (see below)
vplt1, vplt2, nvplt : volume interval over which to plot energy, pressure etc.
                     as well as the number of points in the plot
nevpt             : number of energy-volume points to be inputted
vpt(i) ept(i)    : energy-volume points (atomic units)

Note that the input units are atomic - Bohr and Hartree (NOT Rydbergs).

The equations of state currently implemented are:
1. Universal EOS (Vinet P et al., J. Phys.: Condens. Matter 1, p1941 (1989))
2. Murnaghan EOS (Murnaghan F D, Am. J. Math. 49, p235 (1937))
3. Birch-Murnaghan 3rd-order EOS (Birch F, Phys. Rev. 71, p809 (1947))
4. Birch-Murnaghan 4th-order EOS
5. Natural strain 3rd-order EOS (Poirier J-P and Tarantola A, Phys. Earth
Planet Int. 109, p1 (1998))
6. Natural strain 4th-order EOS
7. Cubic polynomial in (V-V0)

以上是eos的说明，对于bcc，只有一个晶格常数自由度，因此没有问题。hcp有两个自由度a，c。最好保持c/a比例不变来增加体积。关于选用的物态方程是否合理最好看以上每个物态方程的原始文献。eos用的是原子单位，因此体弹模量单位是Ha/au3

赞一下(2人)

回复此楼

2楼2010-12-31 12:54:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖