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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » MS » 关于一篇粘性模拟文章的疑惑
451/1返回列表
查看: 2048  |  回复: 44
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

yeahhanpei

金虫 (正式写手)
已领完
关于一篇粘性模拟文章的疑惑
领取红包 (小木虫手机app专属红包)

扫一扫，下载小木虫客户端

有看过[compressive shear reactive molecular dynamics studies indicating that cocrystals of TNT/CL-20 decreasing sensitivity]这篇文章的么，文章里说模拟TNT的粘性，办法是优化后依次做npt,nvt,nve三个动力学，其中nve是在shear中做的，前两个都是在动力学中设置的。可是是在常温下做的，常温是固体啊，怎么会有粘性这个概念呢？

发自小木虫Android客户端
回复此楼

» 猜你喜欢
1楼 2016-10-18 10:33:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

abinitio

版主 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
39楼: Originally posted by yeahhanpei at 2016-10-18 14:55:28
刚才看了这个证明，但shear的help文件中，没有提到那两个公式，搞不懂了...

我不知道你所说的shear是什么软件，反正我们做动力学的基本上都用这两种方法来算黏度
赞一下(1人) 回复此楼
40楼2016-10-18 15:03:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yeahhanpei

金虫 (正式写手)
有懂得交流下嘛，只抢币不说话

发自小木虫Android客户端
赞一下 回复此楼
33楼2016-10-18 10:36:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
33楼: Originally posted by yeahhanpei at 2016-10-18 10:36:43
有懂得交流下嘛，只抢币不说话

下次要设置成手动发放金币。

------------------------------------

关于问题：在md模拟中，只要不是0K，就会有振动，就会有扩散，所以就会有黏度。固体的黏度要比液态黏度高2-3个数量级。
赞一下 回复此楼
35楼2016-10-18 14:11:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
35楼: Originally posted by abinitio at 2016-10-18 14:11:31
下次要设置成手动发放金币。

------------------------------------

关于问题：在md模拟中，只要不是0K，就会有振动，就会有扩散，所以就会有黏度。固体的黏度要比液态黏度高2-3个数量级。...

可以参考看看本版专家@brucefan的帖子

计算扩散系数的爱因斯坦公式与格林久保公式的等价性之数学证明 - 分子模拟 - 小木虫论坛 http://muchong.com/bbs/viewthread.php?tid=10716181&fpage=1
赞一下 回复此楼
36楼2016-10-18 14:16:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yeahhanpei

金虫 (正式写手)
引用回帖:
35楼: Originally posted by abinitio at 2016-10-18 14:11:31
下次要设置成手动发放金币。

------------------------------------

关于问题：在md模拟中，只要不是0K，就会有振动，就会有扩散，所以就会有黏度。固体的黏度要比液态黏度高2-3个数量级。...

按照标题里的那个文献的方法做了下，文章用的是CL-20晶胞然后建超晶胞做的。我换成了另两种物质建层，然后跑动力学，最后出来的结果不是很好，粘性只有0.7-0.8cP左右，而文献里提到的是约1.8cP.而且这做的是常温的。如果我把温度升到熔点以上，按照通常的理解，温度高粘性应该是会小的。可本来都零点几了，这样再做可能就出不了趋势了（不同质量分数的）。还是说，这种做法本来就有问题
赞一下 回复此楼
37楼2016-10-18 14:38:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yeahhanpei

金虫 (正式写手)
引用回帖:
35楼: Originally posted by abinitio at 2016-10-18 14:11:31
下次要设置成手动发放金币。

------------------------------------

关于问题：在md模拟中，只要不是0K，就会有振动，就会有扩散，所以就会有黏度。固体的黏度要比液态黏度高2-3个数量级。...

那如果我把体系升高温度设置成熔融态的话，粘度就会更小了吧
赞一下 回复此楼
38楼2016-10-18 14:53:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yeahhanpei

金虫 (正式写手)
引用回帖:
36楼: Originally posted by abinitio at 2016-10-18 14:16:15
可以参考看看本版专家@brucefan的帖子

计算扩散系数的爱因斯坦公式与格林久保公式的等价性之数学证明 - 分子模拟 - 小木虫论坛 http://muchong.com/bbs/viewthread.php?tid=10716181&fpage=1...

刚才看了这个证明，但shear的help文件中，没有提到那两个公式，搞不懂了
赞一下 回复此楼
39楼2016-10-18 14:55:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)
引用回帖:
38楼: Originally posted by yeahhanpei at 2016-10-18 14:53:03
那如果我把体系升高温度设置成熔融态的话，粘度就会更小了吧...

应该是。md方法肯定是。
赞一下 回复此楼
41楼2016-10-18 15:03:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yeahhanpei

金虫 (正式写手)
引用回帖:
40楼: Originally posted by abinitio at 2016-10-18 15:03:31
我不知道你所说的shear是什么软件，反正我们做动力学的基本上都用这两种方法来算黏度...

是我没说清楚，不好意思。shear是forcite模块的其中一个任务选项

发自小木虫Android客户端
赞一下 回复此楼
42楼2016-10-18 16:39:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yeahhanpei

金虫 (正式写手)
引用回帖:
40楼: Originally posted by abinitio at 2016-10-18 15:03:31
我不知道你所说的shear是什么软件，反正我们做动力学的基本上都用这两种方法来算黏度...

能不能推荐几篇这种方法的文献，谢谢哦
赞一下 回复此楼
43楼2016-10-18 17:26:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
42楼: Originally posted by yeahhanpei at 2016-10-18 16:39:19
是我没说清楚，不好意思。shear是forcite模块的其中一个任务选项
...

green-kubo 公式就是计算剪切黏度的方法，不知下面的对你是否有用。

The shear viscosity http://www.nyu.edu/classes/tucke ... cture_21/node6.html
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

yeahhanpei
44楼2016-10-18 17:30:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yeahhanpei

金虫 (正式写手)
送红花一朵
引用回帖:
44楼: Originally posted by abinitio at 2016-10-18 17:30:26
green-kubo 公式就是计算剪切黏度的方法，不知下面的对你是否有用。

The shear viscosity http://www.nyu.edu/classes/tuckerman/stat.mech/lectures/lecture_21/node6.html...

嗯，谢谢你，我看下

发自小木虫Android客户端
回复此楼
45楼2016-10-18 19:49:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
10492950392楼
2016-10-18 10:33   回复  
yeahhanpei(金币+1): 谢谢参与
已获得1个金币 发自小木虫IOS客户端
yl火焰3楼
2016-10-18 10:33   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
前方向右向左4楼
2016-10-18 10:34   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
齿轮初学者5楼
2016-10-18 10:34   回复  
yeahhanpei(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
徐二姐6楼
2016-10-18 10:34   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫IOS客户端
9723928787楼
2016-10-18 10:34   回复  
yeahhanpei(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
ammonup8楼
2016-10-18 10:34   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
haodong289楼
2016-10-18 10:34   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
ZHOUSHUNMIN10楼
2016-10-18 10:34   回复  
yeahhanpei(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
heqiang12311楼
2016-10-18 10:34   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
记忆中的风琴12楼
2016-10-18 10:34   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
hanhan20913楼
2016-10-18 10:34   回复  
yeahhanpei(金币+1): 谢谢参与
已获得1个金币 发自小木虫IOS客户端
wyklwyklxin14楼
2016-10-18 10:34   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
youyuxiu15楼
2016-10-18 10:34   回复  
yeahhanpei(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
隐蔽的树16楼
2016-10-18 10:34   回复  
yeahhanpei(金币+1): 谢谢参与
已获得1个金币 发自小木虫IOS客户端
shuaka17楼
2016-10-18 10:35   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫IOS客户端
mzdhgj18楼
2016-10-18 10:35   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
明年我在哪里19楼
2016-10-18 10:35   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
lichong123720楼
2016-10-18 10:35   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫IOS客户端
2010lfb21楼
2016-10-18 10:35   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
bkf9222楼
2016-10-18 10:35   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
kathrinyang23楼
2016-10-18 10:35   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
xiaomagescu24楼
2016-10-18 10:35   回复  
yeahhanpei(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
25499566325楼
2016-10-18 10:35   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
屋后梧桐26楼
2016-10-18 10:35   回复  
yeahhanpei(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
5612228427楼
2016-10-18 10:35   回复  
yeahhanpei(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
叶落了28楼
2016-10-18 10:35   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
cgljuny29楼
2016-10-18 10:35   回复  
yeahhanpei(金币+1): 谢谢参与
已获得1个金币 发自小木虫IOS客户端
ai小爱心30楼
2016-10-18 10:36   回复  
yeahhanpei(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
010308Jing31楼
2016-10-18 10:36   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
dd_123jnby32楼
2016-10-18 10:36   回复  
yeahhanpei(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
阿阿林仔34楼
2016-10-18 10:39   回复  
相关版块跳转 我要订阅楼主 yeahhanpei 的主题更新
451/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 23jxep3nCNZb 2026-09-14 6/300 2026-09-18 11:30 by HjTo14L7GlnE
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-18 10:18 by yS4CkeWDhD6V
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +8 s3fFTmArrBt6 2026-09-13 9/450 2026-09-18 10:04 by yS4CkeWDhD6V
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-18 09:18 by yS4CkeWDhD6V
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 C79jjtjAKjEn 2026-09-15 3/150 2026-09-18 03:40 by cITaGg3p5edr
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-18 00:28 by cITaGg3p5edr
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 7/350 2026-09-17 23:28 by cuPZTDXS3VOd
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 5BDX0d0WFp7t 2026-09-15 4/200 2026-09-17 16:29 by 0ngsZVSCTS7z
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 40ms4Wlo4umh 2026-09-14 5/250 2026-09-17 14:05 by ObVzqQQh1rrE
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 vZfe6xYu34yj 2026-09-14 4/200 2026-09-17 05:50 by k06BKNblrGRH
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-17 03:16 by Tql5LQhh5rLK
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-17 03:14 by Tql5LQhh5rLK
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-16 22:26 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-16 17:27 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
[教师之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
信息提示
关闭
请填处理意见
关闭