24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 2995  |  回复: 12

jmforever123

新虫 (初入文坛)

[求助] 石墨烯片上的水滴的模拟 已有1人参与

各位好,我做的是石墨烯衬底上水滴的模拟,现在做的是平衡态的,在NVT系综下进行,当T到400K的时候,理论上应该有水分子蒸发,但是,我的结果水滴在石墨烯上动来动去,却没有蒸发,水滴为什么会动来动去,会跑到石墨烯片的边上。请各位大神、老师帮忙解答,谢谢。
units          metal
variable       T equal 500
variable       V equal vol
variable       dt equal 0.001
#**************************************************************
dimension      3
boundary       p p p
atom_style     full
read_data      h2o.data
group          water type 1 2
group          C    type 3

pair_style     hybrid lj/cut/tip4p/long 2 1 1 1 0.1577 10.0 8.5 lj/cut 10.0  rebo
pair_coeff     1 1 lj/cut/tip4p/long 0.0000 0.0000
pair_coeff     1 2 lj/cut/tip4p/long 0.0000 0.0000
pair_coeff     2 2 lj/cut/tip4p/long 0.009143 3.1668
pair_coeff     1 3 lj/cut 0.0000 0.0000
pair_coeff     2 3 lj/cut 0.004 3.262
pair_coeff     * * rebo CH.airebo NULL NULL C

kspace_style   pppm/tip4p 1.0e-4

bond_style     harmonic
bond_coeff     1 19.514 0.9572

angle_style    harmonic
angle_coeff    1 2.385 104.52

timestep       0.001
thermo         1000

neighbor       2.0 bin
neigh_modify   delay 0 every 1 check yes
#*********************************************************

min_style      cg
minimize       1.0e-12 1.0e-12 1000 1000

#****************************************************
velocity      all create $T 4928459 mom yes rot yes dist gaussian
thermo_style   custom step temp vol etotal press
#**************************************************************

fix            fixshake water shake 0.0001 20 0 b 1 a 1
fix            NVT all nvt temp $T $T 0.1
dump 1 all custom 1000 dump.myforce.lammpstrj id type x y z
run           500000@oxox6085
回复此楼

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

仿真建模与计算

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

alundilong

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
对真个droplet受力分析吧, F=ma, s = v0+at。
水珠受力不平衡。由于water-C的作用,在给定的初始构型就不能保证平衡的(楼主可以check,吧所有的水分子上的力求和)。
再加上初始速度也不为零吧,所以动来动去的(楼主可以check, 质心的速度到底多少)。

p.s.-:
话说回来,动来动去很正常的,如果不想动,用fix recenter就好了(把整个体系的中心固定起来)。
3楼2016-07-28 02:39:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alundilong

木虫 (正式写手)

【答案】应助回帖

引用回帖:
6楼: Originally posted by jmforever123 at 2016-07-28 01:54:20
还有一个问题就是水分子在400K的时候基本都不蒸发,蒸发是不是需要某个语句。或者是我模拟的步数太短了还是怎么回事。...

不需要其他语句的。
400K蒸不蒸发取决于力场的准确程度(力场是个大问题)
步数少也可能是个问题,但是力场dominate
7楼2016-07-28 23:33:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jmforever123

新虫 (初入文坛)

引用回帖:
9楼: Originally posted by fengguanfeng at 2016-08-12 08:43:40
想请教一下您的石墨烯衬底上的水滴模型是用的什么软件,如何得到的data文件?

你好,我是用的Material Studio,将坐标文件导出,自己写的data文件

发自小木虫Android客户端
10楼2016-08-28 21:24:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

李加伟

荣誉版主 (文坛精英)

不过凡人

优秀版主优秀版主

不懂,觉得很有意思啊
人生海海,在于内心
2楼2016-07-27 17:42:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jmforever123

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by alundilong at 2016-07-28 02:39:06
对真个droplet受力分析吧, F=ma, s = v0+at。
水珠受力不平衡。由于water-C的作用,在给定的初始构型就不能保证平衡的(楼主可以check,吧所有的水分子上的力求和)。
再加上初始速度也不为零吧,所以动来动去的( ...

谢谢你的回复,这个怎么check在lammps说明书有吗,那我好好看看,还有就是,如果check完,受力不平衡怎么处理,谢谢老师

发自小木虫Android客户端
4楼2016-07-28 08:17:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jmforever123

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 李加伟 at 2016-07-27 17:42:31
不懂,觉得很有意思啊

5楼2016-07-28 09:30:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jmforever123

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by alundilong at 2016-07-28 02:39:06
对真个droplet受力分析吧, F=ma, s = v0+at。
水珠受力不平衡。由于water-C的作用,在给定的初始构型就不能保证平衡的(楼主可以check,吧所有的水分子上的力求和)。
再加上初始速度也不为零吧,所以动来动去的( ...

还有一个问题就是水分子在400K的时候基本都不蒸发,蒸发是不是需要某个语句。或者是我模拟的步数太短了还是怎么回事。
6楼2016-07-28 15:54:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jmforever123

新虫 (初入文坛)

引用回帖:
7楼: Originally posted by alundilong at 2016-07-28 23:33:27
不需要其他语句的。
400K蒸不蒸发取决于力场的准确程度(力场是个大问题)
步数少也可能是个问题,但是力场dominate...

哦哦,那我再看看,我觉得可能是受到碳的影响,因为我把水滴单独拿出来模拟,是没有问题的,谢谢你

发自小木虫Android客户端
8楼2016-07-28 23:35:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengguanfeng

木虫 (初入文坛)

想请教一下您的石墨烯衬底上的水滴模型是用的什么软件,如何得到的data文件?
never give up
9楼2016-08-12 08:43:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jmforever123 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 药学求调剂 +11 喽哈加油 2026-04-14 13/650 2026-04-14 21:14 by qingfeng258
[考研] 一志愿沪9,326求生物学调剂 +10 刘墨墨 2026-04-13 10/500 2026-04-14 15:16 by zs92450
[考研] 284求调剂 +17 让我上岸吧阿西 2026-04-09 17/850 2026-04-14 14:44 by 不我拉绿卡
[考研] 生物学调剂 +7 纸扇zhishan 2026-04-13 7/350 2026-04-14 14:21 by jyl0317
[考研] 一志愿哈工大 085600 277 12材科基求调剂 5+5 chenny174 2026-04-10 37/1850 2026-04-14 07:39 by Abskk
[考研] 求调剂288 +7 ioodiiij 2026-04-10 9/450 2026-04-13 08:33 by Hayaay
[考研] 化工调剂求导师收留!一志愿失利,踏实肯干,有植物提取科研经历 +20 yzyzx 2026-04-09 21/1050 2026-04-12 00:12 by 小小小小啦啦啦
[考研] 086003调剂求助 +21 苏弋万 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[考研] 0854调剂 +5 音像店听花鼓戏 2026-04-10 5/250 2026-04-11 10:49 by qingpingzhu
[考研] 282,求调剂 +12 jggshjkkm 2026-04-09 14/700 2026-04-11 09:39 by 猪会飞
[考研] 342电子信息专硕求调剂 +9 你让我怎么荔枝 2026-04-10 10/500 2026-04-11 08:33 by zhq0425
[考研] 083200 305分 求二轮调剂 不接受跨专业 +9 Claireyyyy 2026-04-09 10/500 2026-04-10 21:21 by Claireyyyy
[考研] 一志愿矿大,材料工程专硕314分,0856可调都可以 +15 无懈可击的巨人 2026-04-09 15/750 2026-04-10 18:10 by hmn_wj
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[考研] 332,085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
[考研] 349学科化学045106求调剂,化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 086004 求调剂 309 +7 Yin DY 2026-04-08 7/350 2026-04-09 13:59 by Delta2012
[考研] 085404,334分,求调剂 +5 sunjie8888 2026-04-08 8/400 2026-04-09 07:26 by sunjie8888
[考研] 软件工程求调剂22软工296分求调剂,接受跨调 +4 yangchen2017 2026-04-08 5/250 2026-04-08 21:56 by 土木硕士招生
信息提示
请填处理意见