查看: 400  |  回复: 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的回帖

chenkunkun

铜虫 (小有名气)

请问沥青分子如何建立模型

发自小木虫Android客户端
努力
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的回帖

jluchemlib

铁杆木虫 (正式写手)

【答案】应助回帖

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

luyaobao

金虫 (正式写手)

5楼2018-06-14 07:49:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

邹卓民

新虫 (小有名气)

引用回帖:
2楼: Originally posted by chenkunkun at 2018-06-12 23:39:08
请问沥青分子如何建立模型

MS 建模
6楼2018-06-14 20:58:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

邹卓民

新虫 (小有名气)

引用回帖:
3楼: Originally posted by 无远不界 at 2018-06-13 06:35:52
你的步长已经很小了,如果还有原子丢失:1.thermo modify lost ignore试试,如果愿意ignore的话。2.增大cut-off。3.我也做沥青,建议换换力场pcff试试,cvff之前试过,不怎么成功
...

好的,谢谢我试试看,您方便留个联系方式(QQ)吗,想请教一下您
7楼2018-06-14 21:02:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

邹卓民

新虫 (小有名气)

引用回帖:
3楼: Originally posted by 无远不界 at 2018-06-13 06:35:52
你的步长已经很小了,如果还有原子丢失:1.thermo modify lost ignore试试,如果愿意ignore的话。2.增大cut-off。3.我也做沥青,建议换换力场pcff试试,cvff之前试过,不怎么成功
...

我是用12-component的沥青分子模型,请问一下您的步长一般是设置多少,
8楼2018-06-14 22:10:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

无远不界

铁虫 (正式写手)

引用回帖:
8楼: Originally posted by 邹卓民 at 2018-06-14 22:10:58
我是用12-component的沥青分子模型,请问一下您的步长一般是设置多少,...

我们应该有加QQ,你加我微信吧18729909957

发自小木虫Android客户端
9楼2018-06-16 09:47:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 邹卓民 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 博士面试完了 +29 1204081141 2019-01-14 45/2250 2019-01-17 00:28 by 1204081141
[博后之家] 想做博后但是年纪超龄了 +14 rose_1844100 2019-01-14 19/950 2019-01-16 21:26 by 寒夜微曙
[硕博家园] 北大英语考试又没过,心塞了 +9 yizhidagou 2019-01-14 13/650 2019-01-16 20:49 by yizhidagou
[留学DIY] 德国HZG生物材料研究所Lendlein教授招收一名博士生 +3 kellysun 2019-01-15 8/400 2019-01-16 20:29 by 771731828
[考博] 境外硕士回来考博问题咨询 +4 funnn 2019-01-16 4/200 2019-01-16 18:48 by kingmax18996
[公派出国] 初来荷兰,少许感想。 +18 Chuck LIU 2019-01-13 19/950 2019-01-16 17:32 by lyf1850
[教师之家] 现在还有自己写论文的大学老师吗? +20 玉蝴蝶 2019-01-15 20/1000 2019-01-16 15:39 by nenugdi
[公派出国] 2019西班牙巴塞罗那自治大学招收CSC奖学金博士生最新招生列表 (金币+5) +3 UAB Shanghai 2019-01-15 5/250 2019-01-16 13:51 by UAB Shanghai
[论文投稿] 学术不端问题求助? +9 xuliao2379 2019-01-10 27/1350 2019-01-15 17:09 by xuliao2379
[公派出国] 留学基金委为什么还不出留学规定和简章 +4 dqm_502106 2019-01-14 7/350 2019-01-15 12:31 by 爱睡觉小青年
[硕博家园] 拥有少女心,岁月不饶人 +13 更好88 2019-01-12 27/1350 2019-01-15 09:08 by 永不止步小狼
[硕博家园] SPSS超详细教程求大牛分享 +3 Today03 2019-01-15 3/150 2019-01-15 06:44 by zhangwh2015
[考博] 怎么办?现在还没找到确定博导 +5 苦范儿 2019-01-14 12/600 2019-01-14 15:28 by 苦范儿
[硕博家园] 环境科学考博学校推荐 +3 calp 2019-01-13 5/250 2019-01-14 14:05 by hwysi
[论文投稿] 论文投稿 +7 hgswei 2019-01-12 11/550 2019-01-14 10:57 by 老者归来
[教师之家] 余额宝会亏吗? +13 asdzxc2004 2019-01-12 15/750 2019-01-13 09:57 by toochen
[硕博家园] 天天起床第一件事情就是加油 +8 嘻哈胡杨 2019-01-11 14/700 2019-01-13 01:14 by 嘻哈胡杨
[有机交流] 怎么下载专利 +6 微笑啊你 2019-01-12 6/300 2019-01-12 18:00 by 飘雪孤立
[论文投稿] 发几个2019年光学类期刊影响因子预测 +4 哒哒哒。 2019-01-11 9/450 2019-01-12 16:24 by 哒哒哒。
[有机交流] 柱层析后产物结晶出来去干燥,结果融化了怎么回事 +6 diestars 2019-01-10 10/500 2019-01-12 06:19 by trywpl
信息提示
请填处理意见