24小时热门版块排行榜    

查看: 3395  |  回复: 17

qyshang1989

新虫 (小有名气)

[交流] 请教DPD计算界面张力问题

请大家指教:在MS中用DPD模拟出结果后,得到的结果都是简化后的数据,例如结果中的温度,压力,表面张力这些数据该如何转换成现实中的物理量呢?
回复此楼

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

xuexijisuan

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

lsloneil

专家顾问 (正式写手)

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
qyshang1989: 金币+2 2015-07-10 09:31:51
引用回帖:
5楼: Originally posted by qyshang1989 at 2015-07-09 00:47:59
这本书有没有中文版本?还有个问题,用密度泛函计算界面张力,MS和高斯这两个软件哪个更容易学习啊?多谢!...

你要搞清楚是用哪个密度泛函理论(DFT)算表面张力。人们通常所说的DFT是指量子化学里面的DFT,也就是把能量表示为电子密度的泛函,据我了解量子化学里的DFT很难用来计算表面张力。MS和高斯所算的DFT都是指这种DFT。

另外一种DFT是指液体理论里的DFT,基于统计力学原理,是把自由能表示为粒子密度函数的泛函,可以计算体系很多热力学性质,计算表面张力也是可行的。据我所了解很少有计算化学软件包含了这种DFT的功能,很多课题组都有自己的代码来进行这方面的计算,你可以尝试联系相关课题组看别人是否愿意给代码。
6楼2015-07-10 00:35:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
qyshang1989: 金币+6 2015-07-11 14:20:01
引用回帖:
13楼: Originally posted by qyshang1989 at 2015-07-10 13:18:25
谢谢,我要算的体系文献上的方法就是DFT和状态方程结合,或着Monte Carlo模拟,或者简单的数据回归(这个老师说理论深度不够),我的实验数据已经采集完了,现在就是要计算了。
您讲的这个“压力张量在xx,yy,zz上 ...

你这人怎么老是做伸手党啊?自己搜一下material studio surface tension就会有东西出来。比如这个

http://accelrys.com/events/webin ... 20(NXPowerLite).pdf

参见第29张幻灯片。

或者参考这篇文献
http://pubs.acs.org/doi/abs/10.1021/la702286h

Method部分中的公式。

简而言之如果垂直于界面的方向定为z轴,需要计算的量就是压力张量在xx, yy, zz 方向上的分量,然后表面张力
\gamma = (P_zz - 0.5*(P_xx + P_yy ) )*L_z

L_z is the box length in z direction.

至于你怎么从MS的输出中提取压力张量在xx, yy, zz 方向上的分量请你自己研究吧,我没用过MS,只能帮你到这里。
17楼2015-07-11 11:33:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuntao118

新虫 (著名写手)

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
qyshang1989: 金币+2 2015-07-09 11:54:49
通过无量纲逆变换即可。关键是你要搞明白模拟的初始的无量纲值是怎么得到的,然后你就会换回来了。
2楼2015-07-09 10:57:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
qyshang1989: 金币+2 2015-07-10 12:55:55
引用回帖:
9楼: Originally posted by qyshang1989 at 2015-07-09 13:31:29
自由能密度泛函理论应该就是您说的液体理论里的DFT吧?
我是个初学者,看文献使用密度泛函理论计算的界面张力,以前没做过相关工作,还请多多指教,非常感谢!
我看小木虫上有人说Gromacs可计算界面张力,这个软件 ...

"自由能密度泛函理论应该就是您说的液体理论里的DFT吧?"
是。

我不清楚你到底想算什么,一会儿在MS里跑DPD,一会儿又想着DFT,一会儿又想用gromacs。

流体DFT不属于分子模拟范畴,而是计算液体结构的一种数值方法,所以gromacs这种分子模拟软件是不包含流体DFT的相关计算的。网上有时看到的gromacs和DFT结合是指gromacs和量化计算中的DFT结合进行QM/MM计算。

看你之前的帖子,你似乎是用MS做了个DPD模拟,想求表面张力,那你就应该尽量想办法从你现有的MS的输出文件中求出表面张力,如果像楼上所说只是单位转换的话,那难度很小。

另外使用别人的DFT代码或者使用gromacs重新模拟计算表面张力都需要额外的时间去学习,你要考虑你是否有这么多时间。如果你花了很长时间去学习另外一个软件最后得到的效果和你处理MS数据得到的效果是一样的,那只能说是得不偿失。
10楼2015-07-10 10:57:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

qyshang1989

新虫 (小有名气)

引用回帖:
2楼: Originally posted by chuntao118 at 2015-07-09 10:57:56
通过无量纲逆变换即可。关键是你要搞明白模拟的初始的无量纲值是怎么得到的,然后你就会换回来了。

有相关文献或者详细的教程资料吗?非常感谢!!
3楼2015-07-09 11:54:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuntao118

新虫 (著名写手)

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
qyshang1989: 金币+2 2015-07-09 20:48:27
引用回帖:
3楼: Originally posted by qyshang1989 at 2015-07-09 11:54:36
有相关文献或者详细的教程资料吗?非常感谢!!...

当然有,应该做DPD模拟的比较好的文献都会介绍这些的。如《Self-consistent dissipative particle dynamics algorithm》或者《Introduction to Practice of Molecular Simulation》
4楼2015-07-09 19:12:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qyshang1989

新虫 (小有名气)

引用回帖:
4楼: Originally posted by chuntao118 at 2015-07-09 19:12:45
当然有,应该做DPD模拟的比较好的文献都会介绍这些的。如《Self-consistent dissipative particle dynamics algorithm》或者《Introduction to Practice of Molecular Simulation》...

这本书有没有中文版本?还有个问题,用密度泛函计算界面张力,MS和高斯这两个软件哪个更容易学习啊?多谢!
5楼2015-07-09 20:47:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿Q~~

至尊木虫 (文坛精英)


小木虫: 金币+0.5, 给个红包,谢谢回帖
路过的看了一下,帮顶,祝福你好运!~~~~~~~~~~~~~~~~~~
自强不息,厚德载物;独立精神,自由思想。
7楼2015-07-10 06:45:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuntao118

新虫 (著名写手)

引用回帖:
5楼: Originally posted by qyshang1989 at 2015-07-09 20:47:59
这本书有没有中文版本?还有个问题,用密度泛函计算界面张力,MS和高斯这两个软件哪个更容易学习啊?多谢!...

没用过软件。
8楼2015-07-10 08:23:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qyshang1989

新虫 (小有名气)

引用回帖:
6楼: Originally posted by lsloneil at 2015-07-10 00:35:00
你要搞清楚是用哪个密度泛函理论(DFT)算表面张力。人们通常所说的DFT是指量子化学里面的DFT,也就是把能量表示为电子密度的泛函,据我了解量子化学里的DFT很难用来计算表面张力。MS和高斯所算的DFT都是指这种DFT ...

自由能密度泛函理论应该就是您说的液体理论里的DFT吧?
我是个初学者,看文献使用密度泛函理论计算的界面张力,以前没做过相关工作,还请多多指教,非常感谢!
我看小木虫上有人说Gromacs可计算界面张力,这个软件里包含的DFT理论是液体理论中的DFT吗?
9楼2015-07-10 09:31:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qyshang1989 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 9/450 2026-08-04 03:56 by kf2781974
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 03:35 by k7OM8YghWbkC
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 03:34 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 03:30 by k7OM8YghWbkC
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 03:27 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 4/200 2026-08-04 03:25 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 12/600 2026-08-04 03:17 by k7OM8YghWbkC
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 02:10 by Ocx3OKOzGU4q
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 16/800 2026-08-03 12:23 by fuweiguochen
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 12/600 2026-08-03 08:22 by syl200707
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见