经典系综理论

article/2025/9/21 21:00:51

系综理论是统计力学的理论基础

本质上来说,统计热力学中只有一个问题,即给定能量$E$,如何分布在$N$个全同系统构成的系综上

                             --薛定谔


1.系综的定义

在下面的介绍中我们可以逐渐体会这句话的含义,此时我们暂且把这个定位我们解决的目标问题。之前的波尔兹曼分布/波色爱因斯坦分布/费米分布 已经取得了成功,但是它们针对的是近独立粒子组成的系统。如果粒子之间的相互作用不能忽略的时候,我们就需要系综理论。

在经典的波尔兹曼统计中曾引入单个粒子的相空间,我们称之为“$\mu$空间”。假设离子的自由度是r,则“$\mu$空间”是由r个广义坐标$q_i$和r个广义动量$p_i$组成,张成2r维的相空间。而对于N个(无相互作用的近独立)粒子,我们可以再定义一个“$\Gamma$空间”,那么“$\Gamma$空间”是由rN个广义坐标$q_i$和rN个广义动量$p_i$组成,张成2rN维的相空间。系统在某一时刻的状态,N个粒子的运动状态可以由$q_1,q_2,...,q_{rN};p_1,p_2,...,p_{rN}$,并可以用$\Gamma$空间上的一点来描述。此时可以用我们熟悉的哈密顿正则方程,其中哈密顿量可以表示为$H(q_1,q_2,...,q_{rN};p_1,p_2,...,p_{rN})$。对于保守系统,$H=E$。如果是孤立系统,那么总能量不变:

\[H(q_1,q_2,...,q_{rN};p_1,p_2,...,p_{rN})=E\]

其中独立的变量有$2rN-1$个,可以想象这是一个能量为E的等能面,而且这个面是非常对称的。可以想象一下$x^2+y^2+z^2=1$找找感觉。现在,我们已经拥有了一个2rN维的相空间,空间中的每一点代表了一个N粒子系统的状态。下一步我们还能做些什么?当我们调出一个宏观的参数,我们想知道的是,从微观角度,粒子的状态是怎么分布的。如果只有两个粒子,那只有一个或两个状态,如果考虑的是N个粒子,又如何刻画呢?对于N个粒子的系统,我们可以选择以下两种角度来描述:

a. N粒子系统经历了长时间的演化,每一时刻处于一种微观态,那么它演化的轨迹在$\Gamma$空间中就是一堆密集的点

b. 我们放置了大量的宏观性质完全相同的N粒子系统,它们在$\Gamma$空间中的分布,也还是一堆密集的点

这种思维类似于统计抛硬币得到正反面的概率,你可以用一枚硬币抛1000次,也可以同时抛1000枚硬币。如果我们不缺硬币,那么从效率上看,后者是更为明智的方法。这里我们按照思路b,先给出系综的定义:

大量的处于相同宏观条件、相同力学性质(初始条件可以不同),而各处于某一微观运动状态、并各自独立的系统(子体系)的集合称为系综。

简而言之:系综是系统的集合(系统宏观相同,微观不同)

 

上图是“非常严谨”的系综示意图,蓝色表示张成2rN维$\Gamma$空间,为了表现张成的含义,我画了伞状的高维空间。红色球体代表“系综”,由大量的红色点“系统”组成,系统的空间位置确定它自身的运动状态。

Gibbs的这些创意,都建立在他的一个认识上,力学量的系综平均=时间平均。波尔兹曼尝试用“各态经历假说”来解释。这个假说可以表述为:一个孤立系统从任一初态出发,只要我们给它足够时间,它能经历一切等能面上的微观状态。或者这样来解释假说:一个人在400m跑道上不停地奔跑,经过无穷时间,他可以跑过操场上的任何一点。但数学家证明:“人不可能踏入同一条跑道(河流)”。数学上虽然未被证明,不代表物理上就用错了。实际上,系统很难做到鼓励,外界的微扰会让系统具备跑满整条轨道的能力!

 

2. 微正则/正则/巨正则系统

三种系统的切换实际上是变分问题

 

 

 

转载于:https://www.cnblogs.com/xumh/p/6808889.html


http://chatgpt.dhexx.cn/article/0fLXma6e.shtml

相关文章

lammps教程:nve/nvt/npt系综设置方法

系综是指大量性质和结构完全相同的、处于各种运动状态的、各自独立的系统的集合。 在lammps模拟中,常用的系综有nve、nvt、npt三种系综,下面说说这三种系综的使用方法。 1. nve(微正则系综) nve系综保证系统中原子数量n、体系总…

数值分析(2)-多项式插值: 拉格朗日插值法

整理一下数值分析的笔记~ 目录: 1. 误差 2. 多项式插值与样条插值(THIS) 3. 函数逼近 4. 数值积分与数值微分 5. 线性方程组的直接解法 6. 线性方程组的迭代解法 7. 非线性方程求根 8. 特征值和特征向量的计算 9. 常微分方程初值问题的数值解 1. 拉格朗日基函数 定…

拉格朗日插值法求多项式系数 (附代码)

写在前面: 学了拉格朗日插值法之后发现大家都说可以在O(n^2)时间内得到多项式系数,但是没有找到代码,网上找了很多资料又因为我太弱了没能看懂,最后在emofunx学长的帮助下终于搞明白了。 由于太弱没能看懂的文章 引入 我们都知…

拉格朗日插值多项式的龙格现象

利用插值基函数可以很容易得到拉格朗日插值多项式,公式结构紧凑,在理论分析中甚为重要。但存在以下问题: 当插值节点增减时,计算要全部重新进行,甚为不便(当然现在我们有MATLAB,只需改变下参数…

拉格朗日多项式插值法

在数值分析中,拉格朗日常用于多项式插值。假定提供一组数据点[xi,yi],拉格朗日插值多项式就是由这些 数据的线性运算得到的。 其中基本的多项式有以下公式计算得到 注意 1.第一提供的xi应该是没有相同的,否则不能应用此算法 2.对于每一个xi&…

拉格朗日插值和求多项式系数

拉格朗日介绍 先说说拉格朗日是啥吧 首先 拉格朗日插值是给你 n1 个点 (x,y) 然后根据这n个点可以O(n^2)的求出多项式的系数。也就是解出这个多项式的答案。 假设给你一个多项式 ya0a1*xa2*x^2 然后给你3个解 (x1,y1)(x2,y2)(x3,y3)你第一个想法是怎么解?解方程啊…

数学建模准备 插值(拉格朗日多项式插值,牛顿多项式插值,分段线性插值,分段三次样条插值,分段三次Hermite插值)

文章目录 摘要(必看)0 基础概念什么是插值插值用途什么是拟合插值和拟合的相同点插值和拟合的不同点 1 常用的基本插值方法1.1 多项式插值法1.1.1 拉格朗日多项式插值法多项式插值并不是次数越大越好(龙格现象)分段低次线性插值以…

数模--拉格朗日多项式插值、matlab实现

拉格朗日多项式公式: matlab中插入一个值的代码 function yhlagrange(x,y,xh) %定义拉格朗日插入函数 nlength(x); %统计x和xh的长度 mlength(xh); yhzeros(1,m); %构建一行m列的zero矩阵 c1ones(n-1,1); …

计算方法实验(一):拉格朗日插值多项式

拉格朗日插值数学原理 给定平面上 n 1 n 1 n1个不同的数据点 ( x k , f ( x k ) ) (x_{k},f(x_{k})) (xk​,f(xk​)), k 0 , 1 , ⋯ , n k 0,1,\cdots,n k0,1,⋯,n, x i ≠ x j x_{i} \neq x_{j} xi​​xj​, i ≠ j i \neq j i​j&…

拉格朗日插值多项式在MATLAB中的实现

拉格朗日插值多项式在MATLAB中的实现 Hi! 这是我的一个CSDN博客 例 以下给出了针对题(2)的使用方法 输入: 1.节点值 2.需要插值的原函数 对于题(2)需要设置参数为: xx [0,0.25,0.5,1];% The nodes(n…

lagrange插值法:求拉格朗日插值多项式matlab实现(内附代码及例题)

lagrange插值法:求拉格朗日插值多项式matlab实现(内附代码及例题) 关于拉格朗日插值法相关理论知识,在这里小编不在赘述,请不明白的小伙伴自行百度。小编只负责给出matlab源码。 **例题:**看下面例题(如图): matlab代码: %%%% 求拉格朗日多项式及基函…

拉格朗日插值算法分析

这几天一直研究拉格朗日多项式,今天将自己对拉格朗日多项式的理解写在这里,方便大家交流。 在数值分析中,拉格朗日常用于多项式插值。假定提供一组数据点[xi,yi],拉格朗日插值多项式就是由这些 数据的线性运算得到的。 其中基本…

拉格朗日(Lagrange)插值多项式

题目: 基本原理: 拉格朗日(Lagrange)插值多项式python实现: # encoding: utf-8 from symtable import Symbol X[0.4,0.5,0.6,0.7,0.8] Y[-0.9163,-0.6931,-0.5108,-0.3567,-0.2231] print(X,Y) Lfloat(0.0) x0.54 len…

使用拉格朗日多项式(Lagrangian polynomials)的插值法(python,数值积分)

第三十五篇 拉格朗日多项式插值 插值多项式 首先考虑一个函数的推导,该函数精准地通过一系列np离散数据点。虽然有无限多的函数具备这个条件,但我们将专注于最简单的一个,一个n阶多项式,其中n np−1。我们称这个函数为“插值多…

三点估算法 PERT计划评审技术

三点估算也称PERT法,在计算每项活动的工期时都要考虑三种可能性,计算最悲观的工 期、最可能的工期、最乐观的工期,然后再计算出该活动的期望工期,PERT法计算的是 期望工期. 用PERT法计算工期,我们必须记住下面三个公…

计划评审技术PERT和关键路径法CP

PERT是利用网络分析制定计划以及对计划予以评价的技术。它能协调整个计划的各道工序,合理安排人力、物力、时间和资金,加速计划的完成。PERT网络是一种类似流程圈的箭线圈。它描绘出项目包含的各种活动的先后次序,标明每项活动的时间或相关的…

PERT(计划评审技术,Program Evaluation an Review Technique)

如果你对项目管理、系统架构有兴趣,请加微信订阅号“softjg”,加入这个PM、架构师的大家庭 PERT(计划评审技术,Program Evaluation an Review Technique) 的理论基础是假设项目持续时间以及整个项目完成时间是随机的,且服从某种概…

PERT网络分析法(计划评估和审查技术)

PERT网络分析法(计划评估和审查技术,Program Evaluation and Review Technique 什么是PERT网络分析? PERT(Program Evaluation and Review Technique)即计划评审技术,最早是由美国海军在计划和控制北极星导弹的研制时发展起来的。PERT技术使原先估计的、…

对项目工时的估算----( PERT “计划评审技术” ) 三点估算法

“三点估算法”也称“PERT”法,在计算每项活动的工期时都要考虑三种可能性:计算最悲观的工期、最可能的工期、最乐观的工期,然后再计算出该活动的期望工期,PERT法计算的是期望工期。 用PERT法计算工期,我们必须记住下面…

技术评审

在工作中,我们经常可以听到以下的声音: “我们不进行评审,是因为我们项目比较特殊,没有时间……”。 “我们的项目已经进行了测试,不需要再进行评审了”。 “评审都是在走过场,没有效果……”。 业界公认评…