24小时热门版块排行榜    

查看: 4188  |  回复: 4

plafalines

捐助贵宾 (小有名气)

[求助] Gromacs是如何算溶液的介电常数的啊

求助~~

用的什么公式

计算步骤

多谢
回复此楼

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

分子动力学 模型 - 模拟 - 统计 - 热力学 TOMA'S INDEX gromacs

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)


【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
ghcacj(金币+8): 谢谢 2011-10-25 10:19:25
ghcacj(金币+7): 谢谢 2011-10-26 11:00:06
plafalines(金币+2): 很有用 2011-10-26 12:53:54
引用回帖:
1楼: Originally posted by plafalines at 2011-10-24 02:53:32:
求助~~

用的什么公式

计算步骤

多谢

Google "gromacs dielectric constant", I got:
http://www.icp.uni-stuttgart.de/ ... _II_T2_SS_10_11.pdf



M是体系的总的dipole moment.
是对trajectory中的每一个采样的dipole moment依次平方,加和后再平均;
^2是对所有采样的dipole moment平均后再平方;

- ^2求的是体系的dipole moment采样的均方偏差,即所有样本值在平均值的附近的分布的平均远近情况。如果均方偏差大,说明分的比较散,说明体系的dipole moment的波动很大,那么意味这施加一个外场的话,体系很容易被极化,也就意味着一个大的介电常数。
2楼2011-10-25 05:03:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

plafalines

捐助贵宾 (小有名气)

引用回帖:
2楼: Originally posted by ChemiAndy at 2011-10-25 05:03:26:
Google "gromacs dielectric constant", I got:
http://www.icp.uni-stuttgart.de/ ... _II_T2_SS_10_11.pdf

[eimg]90/24/836907_1319489726_226.png[ ...

那么对于电场是变化的

如何处理

貌似需要用到电动力学

数学真复杂啊
3楼2011-10-25 13:38:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)


★ ★ ★ ★ ★ ★ ★
ghcacj(金币+7): 谢谢 2011-10-26 11:00:12
体系对变化电场的响应就是介电弛豫谱了dielectric relaxation spectroscopy;是通过测定体系对交变电场(正弦变化)的响应,并从峰的情况,来探测体系中不同尺度的粒子的运动情况;电场变化的快慢,即频率不同,体系里面就有不同频率的运动与之相对应,响应,或共振。低频峰对应的是分子dipole moment, 或者cluster的dipole moment,比如一对缔合的离子对的转动驰豫时间;介电弛豫谱中的极高频对应的是体系的电子的瞬时响应,后者这个相对前两者的运动要快的多。

由上面分析,介电弛豫谱一般分为即时响应部分(高频)和弛豫响应部分,弛豫响应部分的变化要慢于电场的变化,落后半拍,所以弛豫谱是一个复数谱。

但我们平常我们谈论的介电常数却是一个常数,这是怎么回事呢?实际上对应的是体系在无限高频率下电场下的介电响应,即高频部分——电子的响应。为什么呢?因为测介电常数用的是不变电场,分子dipole moment受激发后不再变化。而极高频交变电场下,分子跟不上电场变化。都是与电场无关的自发运动造成的总dipole波动。(所谓与电场无关,即用不变电场测到极化响应也要外推到外场为0的情况下,去计算介电常数,否则不同外场的响应幅度不同,不为常数。)
4楼2011-10-26 03:27:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

郭倪倪倪倪

铁虫 (小有名气)

请问这个公式里面用到的数的单位是啥啊,得出来的介电常数有单位吗
5楼2019-10-08 14:58:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 plafalines 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +16 木水思豆 2026-08-25 19/950 2026-08-25 23:56 by dragonxp
[基金申请] 今天务委会开完了,明天出结果吗 +17 angus9576 2026-08-25 21/1050 2026-08-25 22:39 by 悠悠美女
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 2026国自然函评费到账 +20 羊腰板 2026-08-21 23/1150 2026-08-25 16:34 by zsna
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 今天基金会出结果吗?20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见