| 查看: 3028 | 回复: 17 | ||||
qyshang1989新虫 (小有名气)
|
[交流]
请教DPD计算界面张力问题
|
|||
| 请大家指教:在MS中用DPD模拟出结果后,得到的结果都是简化后的数据,例如结果中的温度,压力,表面张力这些数据该如何转换成现实中的物理量呢? |
» 收录本帖的淘帖专辑推荐
xuexijisuan |
» 猜你喜欢
A期刊撤稿
已经有5人回复
临港实验室与上科大联培博士招生1名
已经有8人回复
26申博自荐
已经有7人回复
想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗?
已经有4人回复
带资进组求博导收留
已经有9人回复
求助大佬们,伤口沾上了乙腈
已经有6人回复
最近几年招的学生写论文不引自己组发的文章
已经有9人回复
» 本主题相关价值贴推荐,对您同样有帮助:
请教扭矩计算问题~~
已经有3人回复
请教计算隶属函数的问题
已经有11人回复
请教去除率计算问题
已经有14人回复
请教关于轴向刚度计算公式的问题
已经有8人回复
请教各位虫子,NH3-TPD中NH3的吸附量的计算问题
已经有10人回复
求助materials studio计算dpd 需要的摩尔体积和溶解度参数
已经有12人回复
请教comsol中网格设置的一个问题
已经有3人回复
Gaussian 热化学计算 问题请教
已经有8人回复
请教有关载药量的计算问题
已经有3人回复
表面张力与固液界面张力是不是相关的?
已经有9人回复
请教LSDA+U计算的问题
已经有24人回复
求助固液界面张力相关知识
已经有4人回复
界面张力的计算
已经有5人回复
在DPD中如何计算获得均方末端距?
已经有11人回复
【求助】求教DPD高手怎样用amorphous cell计算Flory huggins参数
已经有8人回复
【求助】请问金属熔体的界面张力与什么相关?谢谢
已经有16人回复
【求助】模拟如何计算 界面张力?MS
已经有21人回复
【求助】请教有效质量和载流子浓度计算方法问题
已经有7人回复
【交流】表面活性剂驱油剂测表面张力和界面张力有何区别
已经有9人回复
lsloneil
专家顾问 (正式写手)
-

专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
qyshang1989: 金币+2 2015-07-10 09:31:51
小木虫: 金币+0.5, 给个红包,谢谢回帖
qyshang1989: 金币+2 2015-07-10 09:31:51
|
你要搞清楚是用哪个密度泛函理论(DFT)算表面张力。人们通常所说的DFT是指量子化学里面的DFT,也就是把能量表示为电子密度的泛函,据我了解量子化学里的DFT很难用来计算表面张力。MS和高斯所算的DFT都是指这种DFT。 另外一种DFT是指液体理论里的DFT,基于统计力学原理,是把自由能表示为粒子密度函数的泛函,可以计算体系很多热力学性质,计算表面张力也是可行的。据我所了解很少有计算化学软件包含了这种DFT的功能,很多课题组都有自己的代码来进行这方面的计算,你可以尝试联系相关课题组看别人是否愿意给代码。 |
6楼2015-07-10 00:35:00
lsloneil
专家顾问 (正式写手)
-

专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
qyshang1989: 金币+6 2015-07-11 14:20:01
小木虫: 金币+0.5, 给个红包,谢谢回帖
qyshang1989: 金币+6 2015-07-11 14:20:01
|
你这人怎么老是做伸手党啊?自己搜一下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
chuntao118
新虫 (著名写手)
- 应助: 93 (初中生)
- 金币: 2590.9
- 散金: 1402
- 红花: 25
- 帖子: 1611
- 在线: 241.2小时
- 虫号: 1791567
- 注册: 2012-05-02
- 专业: 计算数学与科学工程计算
2楼2015-07-09 10:57:56
lsloneil
专家顾问 (正式写手)
-

专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
qyshang1989: 金币+2 2015-07-10 12:55:55
小木虫: 金币+0.5, 给个红包,谢谢回帖
qyshang1989: 金币+2 2015-07-10 12:55:55
|
"自由能密度泛函理论应该就是您说的液体理论里的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
qyshang1989
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 5259.8
- 散金: 29
- 帖子: 157
- 在线: 96.4小时
- 虫号: 1723573
- 注册: 2012-03-28
- 性别: MM
- 专业: 化工热力学和基础数据
3楼2015-07-09 11:54:36
chuntao118
新虫 (著名写手)
- 应助: 93 (初中生)
- 金币: 2590.9
- 散金: 1402
- 红花: 25
- 帖子: 1611
- 在线: 241.2小时
- 虫号: 1791567
- 注册: 2012-05-02
- 专业: 计算数学与科学工程计算
4楼2015-07-09 19:12:45
qyshang1989
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 5259.8
- 散金: 29
- 帖子: 157
- 在线: 96.4小时
- 虫号: 1723573
- 注册: 2012-03-28
- 性别: MM
- 专业: 化工热力学和基础数据
5楼2015-07-09 20:47:59
阿Q~~
至尊木虫 (文坛精英)
- 应助: 77 (初中生)
- 金币: 33343.6
- 红花: 107
- 帖子: 36186
- 在线: 1028.1小时
- 虫号: 3237212
- 注册: 2014-05-27
- 性别: MM
- 专业: 宇宙学

7楼2015-07-10 06:45:45
chuntao118
新虫 (著名写手)
- 应助: 93 (初中生)
- 金币: 2590.9
- 散金: 1402
- 红花: 25
- 帖子: 1611
- 在线: 241.2小时
- 虫号: 1791567
- 注册: 2012-05-02
- 专业: 计算数学与科学工程计算
8楼2015-07-10 08:23:52
qyshang1989
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 5259.8
- 散金: 29
- 帖子: 157
- 在线: 96.4小时
- 虫号: 1723573
- 注册: 2012-03-28
- 性别: MM
- 专业: 化工热力学和基础数据
9楼2015-07-10 09:31:29













回复此楼