24小时热门版块排行榜    

【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 3796  |  回复: 8
【悬赏金币】回答本帖问题,作者邹卓民将赠送您 5 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

邹卓民

新虫 (小有名气)

[求助] lammps采用Green-Kubo计算沥青分子的粘度已有1人参与

lammps采用Green-Kubo计算沥青分子的粘度
出现以下问题:ERROR on proc 0: Bond atoms 3078 3115 missing on proc 0 at step 4275849 (../ntopo_bond_all.cpp:63)

#采用Green-Kubo方法计算沥青分子粘度

#采用real单位系统(记住:长度:A;时间:fs)

#定义几个参数
variable T equal 533.15       #温度为533.15K
variable DT equal 10          #积分步长为10fs
variable V equal vol
#convert from LAMMPS real units to SI

variable kB equal 1.3806504e-23     #J/K Boltzman
variable atm2Pa equal 101325.0
variable A2m equal 1.0e-10
variable fs2s equal 1.0e-15
variable convert equal ${atm2Pa}*${atm2Pa}*${fs2s}*${A2m}*${A2m}*${A2m}


#初始化
dimension 3
units real
boundary p p p
atom_style full
bond_style harmonic
angle_style harmonic
dihedral_style harmonic
improper_style cvff

pair_style lj/cut/coul/long 10.0
kspace_style ewald 1.0e-4
neighbor 2.0 bin
read_data data.vigrin
region box block -50 100 -50 100 -50 100
velocity all create ${T} 12345 mom yes rot yes dist gaussian  

#在NVT系综下平衡体系
fix NVT all nvt temp ${T} ${T} 10  #Nose-Hoover热浴
timestep 0.001
thermo 1000
run 100000

#在产出阶段将用NVE系综
unfix NVT
fix NVE all nve

#在产出阶段计算自关联函数
variable Ns equal 100              #取样间隔
variable Nc equal 100000             #关联数据量
variable Np equal ${Ns}*${Nc}  #产出步数
variable pxy equal pxy
variable pxz equal pxz
variable pyz equal pyz
fix SS all ave/correlate ${Ns} ${Nc} ${Np} &
    v_pxy v_pxz v_pyz type auto file SS.txt ave running

variable M equal 1     #粘度自关联函数的个数
variable N equal ${M}*${Np}

#跑产出阶段
run ${N}
variable scale equal ${convert}/(${kB}*$T)*$V*${Ns}*${DT}
variable v11 equal trap(f_SS[3])*${scale}
variable v22 equal trap(f_SS[4])*${scale}
variable v33 equal trap(f_SS[5])*${scale}

thermo_style   custom step temp press v_pxy v_pxz v_pyz v_v11 v_v22 v_v33
variable v equal (v_v11+v_v22+v_v33)/3.0
variable ndens equal count(all)/vol
print    "average viscosity: $v [Pa.s] @ $T K, ${ndens} /A^3"
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

jluchemlib

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
楼主你试一下用单核运行的时候是否出错
4楼2018-06-13 09:35:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

chenkunkun

铜虫 (小有名气)

努力
2楼2018-06-12 23:39:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

无远不界

铁虫 (正式写手)

你的步长已经很小了,如果还有原子丢失:1.thermo modify lost ignore试试,如果愿意ignore的话。2.增大cut-off。3.我也做沥青,建议换换力场pcff试试,cvff之前试过,不怎么成功

发自小木虫Android客户端
3楼2018-06-13 06:35:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luyaobao

木虫 (著名写手)

5楼2018-06-14 07:49:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 299求调剂 +8 LChogD 2024-03-27 12/600 2024-03-29 14:58 by fyyybin
[考研] 304材料专硕求调剂 +3 烟云墨雨飞啊 2024-03-28 3/150 2024-03-29 10:45 by feizhang1122
[考研] 312求调剂 +3 hffg. 2024-03-27 6/300 2024-03-29 09:37 by fenggou82
[考研] 287求调剂 +6 重获新的我 2024-03-27 6/300 2024-03-29 09:23 by fenggou82
[论文投稿] 一个月了,还没proof,需要催一下吗? +5 214356 2024-03-26 6/300 2024-03-29 08:42 by 莱茵润色
[论文投稿] 请教如何恢复编辑部的意见 30+3 chungea 2024-03-26 5/250 2024-03-28 23:06 by chungea
[考研] 286求调剂 +7 白亦安安 2024-03-24 9/450 2024-03-28 22:36 by 星空星月
[考研] 299求调剂 +3 LChogD 2024-03-28 5/250 2024-03-28 21:43 by LChogD
[考研] 哈工大初试315,有大创经历,科研动手能力强!求调剂!!! +3 隔壁的狗 2024-03-27 3/150 2024-03-28 14:51 by renwenfeng
[考研] 311求调剂 +6 林幸运 2024-03-28 6/300 2024-03-28 14:16 by andm
[考研] 环境科学与工程初试303分求调剂 +4 来瓣儿柚子 2024-03-26 4/200 2024-03-28 11:04 by 小雨哥
[考研] 环境学硕334求调剂 +6 太阳553 2024-03-23 8/400 2024-03-27 20:50 by A我要上学
[考研] 334 化工学硕 求调剂 +6 奋23 2024-03-26 9/450 2024-03-27 16:27 by siwuhen
[考研] 309求调剂 +8 小笨熊12138 2024-03-23 8/400 2024-03-27 16:26 by tkdmxl
[考研] 279求调剂一志愿华东理工材料与化工085602专硕 +4 yt0312 2024-03-24 7/350 2024-03-27 16:21 by barlinike
[考研] 335求调剂 +6 云怎么挽留 2024-03-26 7/350 2024-03-26 18:29 by Whisper890
[考研] 化学工程与技术320求调剂 +12 Wjrjxmt 2024-03-23 14/700 2024-03-26 13:01 by hello呀6
[考研] 294求调剂 +7 迎秋87 2024-03-23 8/400 2024-03-26 11:29 by 201531180024
[考研] 一志愿中国农业大学工学院农业工程350 求调剂 +3 24考研仔 2024-03-24 5/250 2024-03-26 00:28 by 24考研仔
[考研] 329求调剂 +3 YYRr21 2024-03-22 3/150 2024-03-24 13:44 by 旋风马
信息提示
请填处理意见