24小时热门版块排行榜    

查看: 7693  |  回复: 27
本帖产生 1 个 模拟EPI ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lchunmei001

木虫 (正式写手)


[交流] 【求助】gulp中改变温度,晶格常数没变?

1.请教一下大家:我在gulp输入文件中改变温度,为什么相应的晶格常数没有变化?
输入文件如下:
conp opti  prop comp  
cell
4.000000   4.000000   4.000000  90.000000  90.000000  90.000000
fractional    4
Ni    core 0.0000000 0.0000000 0.0000000 0.00000000 1.0000 0.00000              
Ni    core 0.0000000 0.5000000 0.5000000 0.00000000 1.0000 0.00000              
Ni    core 0.5000000 0.0000000 0.5000000 0.00000000 1.0000 0.00000              
Ni    core 0.5000000 0.5000000 0.0000000 0.00000000 1.0000 0.00000              
species
Ni core  0.000
pressure  0.0
temperature  300.0
manybody
Ni core Ni core  0.0 12.0
eam_functional square_root
Ni core  1.000
eam_density power 6
Ni core   729.6928296
lennard 9 6
Ni core Ni core 1303.09841 0.0 0.0 12.0
dump example17.grs

发现temperature随便改为多少,得到的晶格常数都没变?望高手指点指点。我是新手,所以很多不明白!谢谢大家!

2.跑MD是出现以下错误,不知调哪些参数可以解决,试调整了MD相关参数,还是有此错误!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!! ERROR : temperature has exceeded maximum allowed set by mdmaxtemp
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
月只蓝: 金币+1, 鼓励交流! 2014-02-25 11:26:27
用准简谐近似(QHA)的方法考虑温度效应,注意温度不要达到1/2熔点,不然非谐效应主导,用准简谐近似无法算准:

conp free  prop comp
cell
4.000000   4.000000   4.000000  90.000000  90.000000  90.000000
fractional    4
Ni    core 0.0000000 0.0000000 0.0000000 0.00000000 1.0000 0.00000              
Ni    core 0.0000000 0.5000000 0.5000000 0.00000000 1.0000 0.00000              
Ni    core 0.5000000 0.0000000 0.5000000 0.00000000 1.0000 0.00000              
Ni    core 0.5000000 0.5000000 0.0000000 0.00000000 1.0000 0.00000              
species
Ni core  0.000
pressure  0.0
temperature  300.0
shrink 10 10 10
manybody
Ni core Ni core  0.0 12.0
eam_functional square_root
Ni core  1.000
eam_density power 6
Ni core   729.6928296
lennard 9 6
Ni core Ni core 1303.09841 0.0 0.0 12.0
dump example17.grs
21楼2014-02-25 06:02:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 28 个回答
lchunmei001(金币+5): 2011-01-07 17:24:47
lchunmei001(金币+5): 2011-01-07 19:55:37
你这样只是用势函数来优化,温度是不起作用的。
2楼2011-01-07 12:16:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lchunmei001

木虫 (正式写手)


不太明白,因为之前我试另外一个体系的时候就会变(输入控制参数同上,就体系和势不同),还请更详细地指点指点!或者说:针对上面这个输入,怎样修改才会变?谢谢!

[ Last edited by lchunmei001 on 2011-1-7 at 17:33 ]
3楼2011-01-07 17:27:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
zh1987hs(金币+2): 谢谢 2011-02-18 22:59:05
这样优化,晶格常数应该不会随温度而变化的。你可以用模拟的办法。
4楼2011-01-07 17:39:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:44 by XgCC8uTcwILl
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:42 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 GxESXJmptOtQ 2026-09-25 3/150 2026-09-27 06:32 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:08 by XgCC8uTcwILl
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:05 by XgCC8uTcwILl
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:56 by XgCC8uTcwILl
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:56 by XgCC8uTcwILl
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:56 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:53 by XgCC8uTcwILl
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:50 by XgCC8uTcwILl
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:44 by XgCC8uTcwILl
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:44 by XgCC8uTcwILl
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:44 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 6/300 2026-09-27 05:42 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 PGTSwSC3F6nQ 2026-09-24 5/250 2026-09-25 11:19 by yhq9807
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 PGTSwSC3F6nQ 2026-09-24 7/350 2026-09-25 11:18 by yhq9807
[硕博家园] 上海工程技术大学 激光智能制造课题组 2027级博士研究生招生 +4 oxidpl 2026-09-20 6/300 2026-09-24 20:40 by oxidpl
[教师之家] 课题组招2027级博士 上海工程技术大学 激光智能制造方向 +4 oxidpl 2026-09-20 6/300 2026-09-24 20:37 by oxidpl
[基金申请] 师弟论文见刊大半年才想起来申请专利,还能抢救一下吗? +4 13108017953 2026-09-21 4/200 2026-09-24 10:06 by kudofaye
[有机交流] 有机合成以后会不会被AI改变?做科研的虫友怎么看 +5 zapen 2026-09-20 8/400 2026-09-23 17:21 by tianxiaoxian
信息提示
请填处理意见