24小时热门版块排行榜    

查看: 3450  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 中青基了要发朋友圈吗? +3 349506619 2026-08-28 3/150 2026-08-29 10:49 by 木风9012
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 09:48 by G6APbkg8SA6w
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 09:39 by G6APbkg8SA6w
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见