版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (3498)
>虫友互识 (199)
>学术会议 (187)
>基金申请 (178)
>文献求助 (133)
>导师招生 (86)
>论文投稿 (66)
>休闲灌水 (63)
>论文道贺祈福 (33)
>硕博家园 (31)
>招聘信息布告栏 (30)
>考博 (28)
>博后之家 (27)
>教师之家 (21)
>找工作 (21)
>考研 (21)

返回列表

leila3409

捐助贵宾 (正式写手)

应助: 0 (幼儿园)
金币: 616.2
帖子: 310
在线: 71.5小时
虫号: 985885

[交流] 【求助】能量与region的设置有如此大的联系吗？

在设置两个任务时，只将region在z方向的长度差了62.816864-62.81682=0.000044，计算的原子数也相同，但能量的结果出现了巨大的差异(738.8918136-501.4444986=237.447315eV)。有可能是什么原因造成的呢？灰常谢谢啦~~~

# LAMMPS Input File for Grain Boundaries
# Mark Tschopp, Dec2009
# This file will generate numerous input files for LAMMPS
# using a large number of grain boundaries

# ---------- Setup Variables ---------------------
variable etol equal 0
variable ftol equal 1.0e-25
variable maxiter equal 10000
variable maxeval equal 100000
variable latparam equal 2.855312
variable minimumenergy equal -4.122435
variable overlapdist equal 0.275
# ---------- Initialize Simulation ---------------------
clear
units metal
dimension 3
boundary p p p
atom_style atomic

# ---------- Create Atomistic Structure ---------------------
lattice bcc ${latparam}
region whole block 0 63.20502545 -62.302096515  63.20502545 0.000000 62.816864  units box    ####另一个是region whole block 0 63.20502545 -62.302096515  63.20502545 0.000000 62.81682  units box
create_box 2 whole
region upper block INF INF 0.000000 63.20502545 INF INF units box
lattice bcc ${latparam} orient x  3 1 0 orient y  1 -3 0 orient z  0 0 -1
create_atoms 1 region upper
region lower block INF INF -62.302096515  0.000000 INF INF units box
lattice bcc ${latparam} orient x 3 1 0 orient y  -1 3 0  orient z  0 0 1
create_atoms 1 region lower
#region upper2 block  0 9.0  0.000000 9.0 0 2.8  units box
#region lower2 block  0 9.0  -8.9 0.0  0 2.8 units box
group upper region upper
group lower region lower
# ---------- Define Interatomic Potential ---------------------
pair_style  hybrid/overlay eam/alloy  eam/fs table linear 12001
pair_coeff    * * eam/alloy  pot.fe.eam.fs Fe NULL
pair_coeff    * * eam/fs pot.fehe.eam.fs Fe He
pair_coeff    2 2 table pot.he He_He

neighbor 1.0 bin
neigh_modify every 1 delay 0 check yes
# ---------- Displace atoms and delete overlapping atoms ---------------------
delete_atoms overlap ${overlapdist} lower upper
dump       id all atom  100 lammps.lammpstrj
dump_modify id scale no

# ---------- Run Minimization ---------------------
velocity all create 0.0  87287  rot yes mom yes dist gaussian
fix time all dt/reset 1  1.0e-5 1.0e-3 0.2 units box
thermo 10
fix          1 all nve
min_style cg
minimize ${etol} ${ftol} ${maxiter} ${maxeval}

#----------LXC------------------------------------
mass                1    55.845
mass                2    4.00260
group featoms type 1
group heatoms type 2
variable feec equal 4.012986464
variable heec equal 0.007072256332
variable pe equal pe
variable ec equal pe/atoms
variable nfe equal count(featoms)
variable nhe equal count(heatoms)
variable ef equal ${nfe}*${feec}+${nhe}*${heec}+pe
print "This system contains ${nfe} Fe and ${nhe} He atoms"
print "Pe = ${pe}"
print "Ec = ${ec}"
print "Ef = ${ef}"
############################################################

回复此楼

» 猜你喜欢

关于代码变化问题，想知道的进来已经有14人回复
这样的filecode谁见过已经有5人回复
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率已经有4人回复
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有4人回复
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有6人回复
fileCode有新解读？已经有16人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有5人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有4人回复
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有4人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有3人回复

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

查看全部散金贴

1楼 2011-03-31 13:47:51

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

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

H原子在Al2O3 晶体中掺杂， Kpoints 设置已经有6人回复
研究所里的小硕共派出国机会大吗？已经有10人回复
申论答题写错了用黑笔划掉，会被当作弊吗？已经有12人回复
油水界面，加入表面活性剂后，界面会减弱或消失吗已经有10人回复
求助关于荧光共振能量转移的问题？已经有8人回复
请教NEB收敛判据中能量的设置已经有6人回复
【求助】可以通过能量差排除某种物质对反应的影响吗已经有5人回复
【求助】同一个结构，G09和G03 opt计算的能量差别大吗？已经有7人回复
【求助】合成出了380度下高离子电导率的材料, 怎么联系投资商啊已经有13人回复
【求助】镁在常温下能与甲醇或乙醇反应吗？已经有12人回复

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 leila3409 的主题更新