阿掖山
阿掖山

智力活动是一种生活态度 https://mountaye.github.io/blog/

非平衡态物理笔记01 | 朗之万方程和响应函数,热力学平衡和细致平衡条件

来来回回地协调了好几次时间之后,我们最终确定每周三下午进行一次讨论。刚刚过去的这一周讨论了第一讲,时长超过了一个小时。

在讨论会之前,我的主要精力放在了仅有的几张有公式的幻灯片上,结果这一部分完全没有出现在讨论当中。讨论的主要内容是组里本科生艾萨提出的问题:“细致平衡条件 (detailed balance) 和平衡 (equilibrium) 的关系是什么”。这个问题其实挺有质量的,而且我们确实也没有给出一个令人满意的回答。我觉得要想说清楚这个问题,其实需要的是逻辑学和科学哲学,我现在在这方面的水平还不够,这篇文章姑且尝试解释一下。不过不谈数学的原因恐怕不是在这方面没有疑问,毕竟一不用做作业,二没有考试,完全出于兴趣的情况下,直面数学是需要热情和勇气的。

朗之万方程和响应函数

朗之万方程

讲义第12页到15页,讲解了热力学平衡条件下的朗之万方程。所谓热力学平衡条件,真正用到的性质就只有能量均分定理。朗之万方程听起来高大上,实际上只不过是对牛顿第二定律进行了一下改写,明确指明了这个系统受到两种力的作用,下面的公式里,等号右边第一项是正比于速度的衰减项,防止方程的解在无限远处发散(”所以你承认你们科学家是在瞎编咯?”),另外一项是一个关于时间的随机函数,用来刻画不确定性的涨落。

朗之万是和爱因斯坦以及居里夫人(咳咳)同时代的人,依然可以用牛顿时代的东西青史标名,可见物理学果然还是内卷不足,广阔天地,大有可为。(逃)

粒子的位置是粒子之前受力作用的积累:响应函数是卷积的核

所谓牛顿第二定律,对于高中生来说是公式,对于大学生来说是方程,(“公式和方程有啥区别?”),对于刚上大学的高中毕业生来说是噩梦—— F=ma

是方程?这玩意需要解?这玩意也能解?这玩意怎么解?三个字母倒真不算方程,但是物理受到的力往往取决于其所在的位置 x(t) ,所以F其实是 F(x(t)) ;a 是加速度,位置的二阶导数,也就是 d^2 x(t)/dt^2 ,这里的未知量是 x(t)

,这是一个微分方程,常微分方程往往是本科阶段要上一个学期的课程。

响应函数用的是另一套思路,先不从微分方程的角度入手,直接来猜猜解应该长什么样子。此刻的位置与一段时间以前的受力以及当时的那个力对此刻残留的影响有关,这个“一段时间”究竟是多长时间呢?不知道,索性考虑所有时间长度,把这些影响全都加起来。学过微积分应该知道,对无限个无穷小量求和,再取个极限,这就是积分(“快把那个要先证明积分存在并收敛的数学家嘴堵上”),写作:

熟悉数学的朋友知道,等号右边的积分叫做卷积。一般卷积的定义里形式稍有不同,积分上下限分别是正负无穷,这对物理学家不是问题,只要令 τ>t 的时候 h(t−τ)=0 就行了。

一方面,卷积既是一种定义在两个函数(ξ 和 h)之间的运算,可以简写作 x=h∗ξ;另一方面,也可以看作是对 ξ 的变换 (transform), h(t−τ) 是这个变换的核 (kernel),也就是讲义中所讲的 response function。当一个具体问题确定的时候,方程的解就和响应函数一一对应。

在方程实际被求解之前,响应函数和待求解的函数一样,是一个未知量。因为牛顿第二定律是一个二阶微分方程,要想求解,所以很自然的思路就是将方程两边同时做傅里叶变换,将导数转化为含有 iw 的代数方程,然后再将解做傅里叶逆变换。因为两个函数的卷积的傅里叶变换,等于两个函数各自的傅里叶变换的乘积,所以在傅里叶倒空间里,方程的解、响应函数、随机涨落三项是很简单的乘积形式。

关于朗之万方程的解,讲义中不清楚的地方

第14页右侧蓝字似乎是 “spectral density”.

第15页 ~C(ω) 似乎是 ⟨~v(ω)2⟩

平衡 (equilibrium) 和细致平衡 (detailed balance)

先贴几个相关的链接:

目前直接在谷歌搜索 “equilibrium” 并进入维基百科的话会发现找不到和物理有关的词条,因为 “mechanical equilibrium”, “thermal equilibrium”, “thermaldynamic equilibrium” 是几个独立的词条,而且都没有链接到 “equilibrium” 的消歧义页面上。

所谓 “mechanical equilibrium”,就是受力平衡,分为三种:稳定平衡、不稳定平衡、随遇平衡。实际上就是牛顿第二定律作为一个动力学方程,其中的稳定、不稳定和中性三种驻点 (fixed point),和这一讲的关系不大。

热平衡简单。两个 系统处于热平衡,那么两个系统的温度相同。一个 系统处于热平衡,指的是这个系统温度时间上恒定,空间上均一。虽然这样说有点鸡贼了,盗用了读者在日常生活中对“温度”的理解。

热力学平衡的条件严格一些,不止热平衡,还要达到机械平衡,化学平衡,所有宏观上的热力学物理量(温度、质心坐标、粒子数……)都没有变化。

据此,大概可以归纳物理中一般性的“平衡”概念:某种物理性质在一个系统中的平衡,指的是这种物理量在系统内时间上的恒定,空间上的均匀;在多个系统之间的平衡,指的是这一物理量在这些系统之间没有流动。

其实上面这一段的草稿里,“指的是这个系统温度时间上恒定,空间上均一”后面想当然地来了一句“也就是任意两个子系统之间相互平衡”。当你听到了两个系统之间的平衡之后,再想到一个系统可以划分为多个子系统,思维的惯性就会用类比的方法来如此定义一个系统的平衡。

但是其实仔细一想,“时间上恒定,空间上均一”并不一定要求任意两个子系统之间相互平衡,但是任意两个子系统互相平衡是一个很有趣的状态,我们给他起名叫“细致平衡”。当任意两个子系统互相平衡的时候,可以确定地达到时间上恒定空间上均一的结果。

所以这是我对艾萨问题的回答:细致平衡条件是平衡的充分不必要条件。之所以常有这样的疑惑,是因为平衡的定义不能递归,一个系统的平衡和多个系统之间的平衡定义不同。

上面最后一个链接可以说是直接回答了艾萨的问题,但是这个页面的主标题是 “Equilibrium Means Detailed Balance”, 字面意思是“平衡意味着细致平衡”,也就是说平衡是细致平衡的充分条件,和我前面说的相反。但是文章的内容实际并没有讨论究竟是谁蕴含谁,所以姑且可以无视这个标题。

细致平衡条件在维基百科词条的正文中称作 “principle of detailed balance”,也就是说这是一个“原理”,和“定律”(law) 一样,这是科学 理论 这一逻辑体系的公理,是没办法在理论内部通过逻辑推导来证明或者证伪的。我们之所以采用这些原理(也就是认为他们是真命题),是因为现有的 实验 全都没有证明这些原理不正确。

多说两句,科学虽然存在理论,但归根结底是一种经验科学,必须落地于实验和观测。(这里的“经验”对应于“理性”,是有严格定义的哲学概念,不是我们日常用语中所说的那个经验。)而观测和实验总有认知的的边界,边界以外的事情科学家做不了保证,所以实验总没有办法一劳永逸地说某某是对的,只能说某某还没错。只要证而不伪,科学就说这是一个真命题。科学与非科学的界限,根据波普尔的科学哲学理论,就在于可证伪性

另外讨论课里完全没有涉及的一点是,这里的“子系统”必须依然是“宏观的”(macroscopic),各个子系统只能是小一号的总系统,(用软件工程的话说,定义的接口是一样的),不涉及组成这个系统的微观机制(不能是另一类对象)。比如说一盒理想气体,子系统只能是给盒子里加上隔板后的几小盒理想气体,不能精确到具体的气体原子分子,否则此时连机械平衡都满足不了。

这和经典力学里的“宏观无穷小”这一扭曲的概念有异曲同工之妙。归根结底,是因为宏观的科学理论建立得不够早,已经用上了微积分;但又太早,没考虑也不知道更低层次的物理规律。经典力学的更低层次已经形成了量子力学,而热力学的底层统计物理就只能和热力学挤在一起,和理论力学、电动力学、量子力学并列为四大力学,逃不出“N 大某某有 (N+1) 个”的魔咒。

CC BY-NC-ND 2.0 版权声明

喜欢我的文章吗?
别忘了给点支持与赞赏,让我知道创作的路上有你陪伴。

加载中…
加载中…

发布评论