查看: 6503  |  回复: 8

largehoo

木虫 (小有名气)

[交流] 【搬运】系综小结 已有7人参与

系综(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. 厦门大学固体表面物理化学国家重点实验室.网络下载

转自量子化学网http://www.quantumchemistry.net/node/123
原作者 Elizerbeth
回复此楼

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

第一性原理相关文档 分子动力学模拟 分子动力学 分子动力学模拟过程

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

追梦2016

铜虫 (小有名气)

有所收货,谢谢分享。
只要肯学,没有学不会的
2楼2016-05-26 14:53:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
谢谢分享!
3楼2016-05-26 15:02:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

一次性保温杯

金虫 (小有名气)

送红花一朵
感谢
4楼2016-05-26 17:26:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
很老的内容了,是@zhangsolomon 总结的。
5楼2016-05-27 08:40:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paul121

银虫 (初入文坛)

学习啦
6楼2016-06-04 17:54:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tina张

新虫 (小有名气)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
5楼: Originally posted by abinitio at 2016-05-27 08:40:09
很老的内容了,是zhangsolomon 总结的。

您好,我想问一下,我在模拟水的时候,用nve系综,用fix  9 thermup langevin 300.0 300.0 0.1 89994827  tally yes控制温度,但是温度还是往上升,已经升到350了
不忘初心
7楼2016-11-29 21:54:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
7楼: Originally posted by tina张 at 2016-11-29 21:54:31
您好,我想问一下,我在模拟水的时候,用nve系综,用fix  9 thermup langevin 300.0 300.0 0.1 89994827  tally yes控制温度,但是温度还是往上升,已经升到350了...

CODE:
2. 微正则系综(micro-canonical ensemble),简写为NVE,即表示具有确定的粒子数(N)、体积(V)、总能量(E)。微正则系综广泛被应用在分子动力学模拟中。假定N个粒子处在体积为V的盒子内,并固定总能量(E)。此时,系综的温度(T)和系统压强(P)可能在某一平均值附近起伏变化。

所以,如果你想控制温度,应该用nvt或npt系综
8楼2016-11-30 08:30:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jikaiheaven

新虫 (初入文坛)

楼主,感谢了
9楼2017-09-07 21:34:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 largehoo 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 申请科研项目纯靠实力能成功吗? 200+7 David1982 2026-10-13 11/550 2026-10-15 10:30 by zyqchem
[职场人生] 40+中年人失业后,还能找到好工作吗? +5 xqy2009 2026-10-13 5/250 2026-10-15 09:49 by 淘淘塑物性通
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见