24小时热门版块排行榜     意得辑英文润色,师兄用过都说好

查看: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 核心为什么难发? +13 阿童木? 2018-10-17 14/700 2018-10-19 16:51 by 可爱的小帅锅
[找工作] 拜耳,默沙东,诺和诺德,罗氏 +4 牛奶维他命 2018-10-17 9/450 2018-10-19 16:36 by changray
[外语学习] 想提高口语,也许你可以进来看看? +3 不吃鱼丸 2018-10-18 7/350 2018-10-19 15:23 by 不吃鱼丸
[教师之家] 找工作 +5 马和平 2018-10-19 5/250 2018-10-19 13:41 by 放牛的小孩
[论文投稿] ASCE旗下期刊,投稿不到一个月,Submissions with a Decision,是不是不妙? +3 candypuck 2018-10-19 6/300 2018-10-19 13:14 by 国际科学编辑
[留学DIY] 各位道友,我想问下有关伯克利的申请 +5 皇甫大人 2018-10-19 7/350 2018-10-19 12:16 by 419765734
[有机交流] 核磁碳谱溶剂 10+3 yugw0214 2018-10-18 8/400 2018-10-19 10:21 by celan
[硕博家园] 垃圾导师 +41 念北LL 2018-10-18 49/2450 2018-10-19 10:01 by 1464736696
[论文投稿] elsevier投稿的一个问题 +4 qgking 2018-10-18 7/350 2018-10-19 07:42 by rohdy
[硕博家园] 工程博士 +5 我也搞学术 2018-10-16 5/250 2018-10-18 16:47 by skyking20
[硕博家园] 路边的女孩 +8 cxh90 2018-10-16 9/450 2018-10-18 14:31 by manbu2017
[考博] 老师什么意思,希望大吗 +7 爱7号的me 2018-10-17 12/600 2018-10-18 11:40 by 爱7号的me
[硕博家园] 女孩子,物理系硕士毕业可以做什么呢 +16 浅时言 2018-10-18 16/800 2018-10-18 11:13 by 博出位来
[博后之家] 湖南大学生物学院谭志凯课题组诚聘(生物材料/组织工程方向)博士后 +3 _耳东陈_ 2018-10-15 4/200 2018-10-18 10:17 by --LUCKY--
[硕博家园] 人应不应该一条道走到黑? +6 海上的小船 2018-10-14 9/450 2018-10-17 17:20 by Sponge丶
[硕博家园] 有没有想读博又怕毕业不了的 +19 18292454471 2018-10-15 23/1150 2018-10-16 20:57 by 18292454471
[分子生物] RNA提取 +4 差点是个镁铝 2018-10-14 4/200 2018-10-15 13:26 by 墨秦墨秦
[硕博家园] 急!有偿 寻求熟悉Material Studio的同学交流指导吸附扩散模拟 +3 Skyezj 2018-10-12 6/300 2018-10-15 10:22 by 风在飘丿
[版块工场] 【早起签到贴】2018年10月14日(q q 48626260 进群申请:小木虫——您个人ID) +56 8475 2018-10-14 58/2900 2018-10-15 07:02 by hnhbtj
[版块工场] 【早起签到贴】2018年10月13日(q q 48626260 进群申请:小木虫——您个人ID) +35 8475 2018-10-13 36/1800 2018-10-14 06:08 by hnhbtj
信息提示
请填处理意见