24小时热门版块排行榜    

查看: 5353  |  回复: 2

h-bond

金虫 (正式写手)

[求助] 系综如何选择 已有2人参与

如何确定是选NVT还是NVE或者MS-MD-NVT,有什么区别?
如果模拟过程中有化学反应发生,那么选NVE是否合理?
回复此楼

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

晶体模型计算分析 学学学

» 猜你喜欢

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

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

abinitio

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
h-bond: 金币+5, ★★★很有帮助 2015-02-06 23:52:40
作者:solomon    文章来源:http://solomon.blog.edu.cn
系综(ensemble) 是指在一定的宏观条件下(约束条件),大量性质和结构完全相同的、处于各种运动状态的、各自独立的系统的集合。全称为统计系综。系综是用统计方法描述热力学系统的统计规律性时引入的一个基本概念;系综是统计理论的一种表述方式,系综理论使统计物理成为普遍的微观统计理论 ;系综并不是实际的物体,构成系综的系统才是实际物体。
约束条件是由一组外加宏观参量来表示。在平衡统计力学范畴下,可以用来处理稳定系综。

一、常用系综分类
根据宏观约束条件,系综被分为以下几种:
1. 正则系综(canonical ensemble),全称应为“宏观正则系综”,简写为NVT,即表示具有确定的粒子数(N)、体积(V)、温度(T)。正则系综是蒙特卡罗方法模拟处理的典型代表。假定N个粒子处在体积为V的盒子内,将其埋入温度恒为T的热浴中。此时,总能量(E)和系统压强(P)可能在某一平均值附近起伏变化。平衡体系代表封闭系统,与大热源热接触平衡的恒温系统。正则系综的特征函数是亥姆霍兹自由能F(N,V,T)。
2. 微正则系综(micro-canonical ensemble),简写为NVE,即表示具有确定的粒子数(N)、体积(V)、总能量(E)。微正则系综广泛被应用在分子动力学模拟中。假定N个粒子处在体积为V的盒子内,并固定总能量(E)。此时,系综的温度(T)和系统压强(P)可能在某一平均值附近起伏变化。平衡体系为孤立系统,与外界即无能量交换,也无粒子交换。微正则系综的特征函数是熵S(N,V,E)。
3. 等温等压(constant-pressure,constant-temperature),简写为NPT,即表示具有确定的粒子数(N)、压强(P)、温度(T)。一般是在蒙特卡罗模拟中实现。其总能量(E)和系统体积(V)可能存在起伏。体系是可移动系统壁情况下的恒温热浴。特征函数是吉布斯自由能G(N,P,T)。
4. 等压等焓(contant-pressure,constant- enthalpy),简写为NPH,即表示具有确定的粒子数(N)、压强(P)、焓(H)。由于由于H =E+PV,故在该系综下进行模拟时要保持压力与焓值为固定,其调节技术的实现也有一定的难度,这种系综在实际的分子动力学模拟中已经很少遇到了。
5. 巨正则系综(grand canonical ensemble),简写为VTμ,即表示具有确定的粒体积(V)、温度(T)和化学势(μ)。巨正则系综通常是蒙特卡罗模拟的对象和手段。此时、系统能量(E)、压强(P)和粒子数(N)会在某一平均值附近有一个起伏。体系是一个开放系统,与大热源大粒子源热接触平衡而具有恒定的T,。特征函数是马休(Massieu)函数J(μ,V,T)。

二、系综调节
系综调节主要是指在进行分子动力学计算过程中,对温度和压力参数的调节,分为调温技术和调压技术。
1.      调温技术
    在 NVT 系综或 NPT系综中,即使在 NVE系综模拟的平衡态中,也经常调整温度到期望值。如果希望知道系统的平衡态性质怎样依赖于温度,那么就必须在不同的温度下进行模拟。
目前实现对温度的调节有 4种方式:速度标度、Berendsen热浴、Gaussian热浴、Nose—Hoover热浴。

2.      调压技术
在等压模拟中,可以通过改变模拟原胞的三个方向或一个方向的尺寸来实现体积的变化.类似于温度控制的方法,也有许多方法用于压力控制,总的来说有以下 3种技术:Berendsen方法、Anderson方法、Parrinello-Rahman方法。

三、系综选择
原则上巨正则系综应用最广,但却不一定是最方便的。因为可以看到三种系综的演化过程既是约束解除的过程,却是以增加变量为代价的,这也就增加了数学上的复杂性。因此一般情况下如果不需求解。,则不必使用巨正则系综。
系综选择的基本原则为:
1.    微正则系综能够简单的求得近独立,全同,定域粒子系统,并且每个粒子只能有两个不同的可能状态,例如简单的铁磁,顺磁模型
2.    微正则系综难求的系统,可用正则系综求解
3.    当微正则和正则系综均难求时,可用巨正则系综求解

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

参考资料:
1. D.罗伯,计算材料学。化学工业出版社
2. 文玉华等。分子动力学模拟的主要技术。力学进展,2003,33:65
3. 潘娜娜。统计物理期中小结。网络下载
4. 维基百科。http://www.wiki.cn/wiki/%E7%B3%BB%E7%BB%BC
5. Statistical Thermodynamics and Chemical Kinetics. 厦门大学固体表面物理化学国家重点实验室.网络下载
2楼2015-02-06 09:04:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
有化学反应的话建议NVT或者NPT。否则用NVE的话能量不守衡,采样得到的结果没有任何意义。
3楼2015-02-07 15:02:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 h-bond 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 什么时候开奖? +4 CrisMessi 2026-08-18 4/200 2026-08-19 00:53 by gatelove
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 27/1350 2026-08-18 21:44 by Noways
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +3 瞬息宇宙 2026-08-15 3/150 2026-08-18 21:40 by xiaohuixin
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 9/450 2026-08-18 20:09 by 奶牛小黑
[基金申请] 朋友圈看到的 +4 wangzilk 2026-08-18 5/250 2026-08-18 19:25 by wangzilk
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +51 医学老男孩 2026-08-13 110/5500 2026-08-18 19:07 by icerianxia
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
请填处理意见