24小时热门版块排行榜

返回列表

本帖产生 1 个模拟EPI ，点击这里进行查看

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

lchunmei001

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 1738.5
帖子: 797
在线: 151.3小时
虫号: 863507

[交流] 【求助】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
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

回复此楼

» 收录本帖的淘帖专辑推荐

天地灵气，日月精华

gulp

» 猜你喜欢

职称评审没过，求安慰已经有41人回复
回收溶剂求助已经有7人回复
硝基苯如何除去已经有3人回复
A期刊撤稿已经有4人回复
垃圾破二本职称评审标准已经有17人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
求助文献已经有3人回复
三无产品还有机会吗已经有6人回复

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

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

查看全部散金贴

1楼 2011-01-05 22:25:01

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

老虎大王

木虫 (著名写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+10, 模拟EPI+1): 谢谢 2011-02-24 13:33:37

你好。昨晚我用你的输入文件作了测试，好像确实有问题。我另换了Johnson势函数，也改过几次控温参数，均发现有问题，有时候是温度，有时候是压力。但这样的问题对于晶体的模拟是不该出现的，我一时也看不出是程序有Bug还是输入有问题。

你用的Sutton-Chen势太长了，运行很慢，我没有时间继续测试下去，我建议在找到问题症结之前，你先把Gulp程序中的mxmdtemp参数增大，在initial.f90中找到
rmdmaxtemp = 100.0_dp，
把100改为1000，这样容错性就好一些。注意改后要重新编译生成新的可执行文件。
但如果这样还是有问题，那就真的是不好办了。另外，请你打开mdprop.f90这个文件，找到大约127行处，有如下代码：
  if (tempratio.gt.rmdmaxtemp) then
call outerror('temperature has exceeded maximum allowed set by mdmaxtemp',0_i4)
call gulpfinish
  endif

你还可以在call gulpfinish 这一句之前（即call outerror这一句之后），添加如下代码：
   write(ioout,'(''  ** Time : '',f18.5,'' ps :'')') nsteps*tstep(ncf)
   write(ioout,*) 'temperature, targettemperature, tempratio:'
   write(ioout,*) temperature, targettemperature, tempratio

改后要重新编译。这个是在出错的时候输出错在什么地方。这一句对改错是没什么用的，只是弄明白错的程度有多大，如果输出的tempratio不是特别大（比如你前面设的rmdmaxtemp是1000，这里输出的tempratio是1020）你就可以再加大rmdmaxtemp，比如设成1200之类，如果tempratio大得很，比如2000，3000之类，你就没有必要再搞了，肯定问题很大。

再一个你的输入文件太麻烦了。对FCC金属，可以这样写：
nosy cart md conp
cell
4.09  4.09  4.09  90.000000  90.000000  90.000000
frac
Ag 0.0  0.0  0.0
Ag 0.5  0.0  0.5
Ag 0.5  0.5  0.0
Ag 0.0  0.5  0.5
supercell  7 7 7
libr  suttonchen
下面就是pressure=0这些，是你的MD控制参数。
这样你换其它FCC体系也容易，只要改晶格参数、元素名就可以了。suttonchen势是在Gulp的势函数库当中的。用的时候把Libraries中的suttonchen.lib拷到当前运行目录下。用7*7*7的晶格也太大了，大概256或500个，对于晶格优化也就差不多了，864也可以。时间步数大概0.001ps就够了。你对金属Ag，应该在熔点以下运行，温度不要太高。

当然输入文件的写入是小事，关键还是温度控不住的情况比较烦人。而且一般是突然发生的错误。实在不行你可以给程序作者发Email，说说这个问题，这个作者Julian Gale是个很好、很有耐心的人。要不就换DLPOLY做MD吧，DLPOLY计算MD专业一点。呵呵呵。