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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 石墨烯片上的水滴的模拟
汕头大学海洋科学接受调剂
131/2返回列表12下一页
查看: 2989  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

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

仿真建模与计算

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2016-07-27 15:41:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

alundilong

木虫 (正式写手)

【答案】应助回帖

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

p.s.-:
话说回来，动来动去很正常的，如果不想动，用fix recenter就好了（把整个体系的中心固定起来）。
赞一下(1人) 回复此楼
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
赞一下(1人) 回复此楼
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客户端
赞一下(1人) 回复此楼
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 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 2026硕士调剂_能动_河南农业大学 +4 河南农业大学-能 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[考研] 293求调剂 +16 我爱高数高数爱� 2026-04-12 18/900 2026-04-13 21:47 by 学员JpLReM
[考研] 材料复试求调剂 +24 xhhdjdjsjks 2026-04-09 24/1200 2026-04-13 15:49 by 幸免 ..
[考研] 生物学调剂，一志愿西南大学348，Top期刊一区二作、二区三作，三等奖学金三次 +5 candyyyi 2026-04-09 5/250 2026-04-13 09:02 by 可淡不可忘
[考研] 。 +10 李多米lee. 2026-04-12 11/550 2026-04-12 22:58 by yuyin1233
[考研] 电气工程专硕320求调剂 +5 小麻子111 2026-04-10 5/250 2026-04-12 10:47 by zhouyuwinner
[考研] 求调剂 +18 MAX怅惘 2026-04-09 20/1000 2026-04-11 23:31 by zhen～
[考研] 352 求调剂 +6 yzion 2026-04-11 8/400 2026-04-11 16:24 by 明月此时有
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 283求调剂，工科！ +12 苏打水7777 2026-04-08 12/600 2026-04-11 10:28 by 逆水乘风
[考研] 297求调剂 +9 Kwgyz 2026-04-09 9/450 2026-04-11 10:09 by zhq0425
[考研] 一志愿沪9，326生物学求相关专业调剂 +4 刘墨墨 2026-04-09 4/200 2026-04-10 12:07 by pengliang8036
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +20 慕绝cc 2026-04-09 24/1200 2026-04-10 10:22 by xujun0624
[考研] 已调剂 +18 柴郡猫_ 2026-04-09 19/950 2026-04-09 22:10 by 柴郡猫_
[考研] 085400电子信息类（川大控制工程）求调剂可跨专业 求老师联系 +3 626776879 2026-04-08 3/150 2026-04-09 16:05 by 猪会飞
[考研] 265求调剂 +4 风说她早忘了 2026-04-07 4/200 2026-04-09 13:59 by only周
[考研] 0860004 求调剂 309分 +6 Yin DY 2026-04-09 6/300 2026-04-09 10:19 by 啊李999
[考研] 生物学学硕，初试351分，求调剂 +4 …～、王…～ 2026-04-08 5/250 2026-04-08 21:49 by limeifeng
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 22:00 by hemengdong
[考研] 318求调剂 +5 李青山山山 2026-04-07 5/250 2026-04-07 18:24 by 蓝云思雨
信息提示
关闭
请填处理意见
关闭