24小时热门版块排行榜

>论坛更新日志 (3591)
>导师招生 (436)
>文献求助 (374)
>休闲灌水 (294)
>虫友互识 (136)
>考博 (123)
>招聘信息布告栏 (92)
>论文道贺祈福 (88)
>硕博家园 (84)
>公派出国 (77)
>博后之家 (75)
>论文投稿 (75)
>基金申请 (63)
>教师之家 (59)
>考研 (41)
>找工作 (35)

返回列表

mansamsung8113

木虫 (正式写手)

应助: 2 (幼儿园)
金币: 3082.1
红花: 1
帖子: 562
在线: 208.6小时
虫号: 712068
注册: 2009-03-01
性别: GG
专业: 金属材料的力学行为

[求助] 关于结构参数的精修

看了侯老师的教程，有一个是关于Mg的参数精修的例子，最后一步教程上说“在该体积下，重复上面第2和3步，得到平衡状态下的其他晶胞参数。这一步就是：在得到了E(V)曲线后，通过拟合状态方程得到平衡状态下的体积，计算出上面脚本中的变量$i的值，并改变$i的循环值，再运行run_cell计算一次，得到其他的结构参数c和位置u.。
我的理解是在第一次运行完run_cell，用murn.x拟合E-V曲线，得到平衡状态的体积V0和晶格常数a0后，然后在把run_cell的 $i 值就取为a0后，再运行run_cell 吧?
第二次运行的run_cell时，对应的脚本开头是不是应该这样:
#!/bin/sh
rm WAVECAR
for i in a0
do
cat > POSCAR < Mg-hex
$i
0.0 -1.0 0.0
0.8660254037844 0.5 0.0
0.0 0.0 1.6230529595
.....................
还请各位大侠赐教

回复此楼