欢迎光临散文网 会员登陆 & 注册

高维球与分子速率分布||Pi Day算个Pi

2021-03-14 00:27 作者:湮灭的末影狐  | 我要投稿

//今天是2021年3月14日,传说中的%5Cpi%0A日。

//同时也是物理学之神A.Einstein的142周年诞辰。

//以及是玩黑洞的S.J.Hawking逝世3周年。

//在这一天我想写一些不太一样的内容。我告诉你n维球的表面积不仅与%5Cpi有关,还和气体分子速率分布有关,你是否会觉得奇怪?那么往下看下去吧...

常数%5Cpi,被定义为圆的周长与直径之比,广泛出现在各种计算、公式、模型中,是最广为人知的无理数、超越数

而在热学中我们会了解到,单组分、平衡态的气体,分子在某个方向的速度分布满足以下麦克斯韦速度分布(这是一个概率分布函数):

f(v_i)%3D%5Csqrt%7B%5Cfrac%7Bm%7D%7B2%5Cpi%20kT%7D%7De%5E%7B-%5Cfrac%7Bmv_i%5E2%7D%7B2kT%7D%7D

其中,m是气体分子的质量,T是温度,而k是玻尔兹曼常数。i可以取x%2Cy%2Cz.

我们可以看到这是一个正态分布的形式。这很合理对吧。所以我们遇到了第一个关于%5Cpi%0A的问题:速率分布函数里面的%5Cpi是怎么出现的?

我们知道,作为一个概率分布函数,是必须满足归一化的:

%5Cint_%7B-%5Cinfty%7D%5E%7B%5Cinfty%7Df(v)%5Cmathrm%7Bd%7Dv%3D1%20

也就是说,一个分子出现这一方向所有可能速度分量的概率总和为1,这很显然。

在求积分%5Cint_%7B-%5Cinfty%7D%5E%7B%5Cinfty%7D%20e%5E%7B-%5Calpha%20x%5E2%7D%5Cmathrm%7Bd%7Dx 时,被积函数是没有初等形式的原函数的,但是可以用到一个巧妙的操作,将它化为对整个二维平面的积分,我在费曼讲义上第一次见到这种操作:

I%3D%5Cint_%7B-%5Cinfty%7D%5E%7B%5Cinfty%7D%20e%5E%7B-%5Calpha%20x%5E2%7D%5Cmathrm%7Bd%7Dx%3D%5Cint_%7B-%5Cinfty%7D%5E%7B%5Cinfty%7D%20e%5E%7B-%5Calpha%20y%5E2%7D%5Cmathrm%7Bd%7Dy%20

%20%5CRightarrow%20I%5E2%20%3D%20%5Ciint_%7B%5Crm%20All%7De%5E%7B-%5Calpha%20x%5E2%7De%5E%7B-%5Calpha%20y%5E2%7D%20%5Cmathrm%7Bd%7Dx%5Cmathrm%7Bd%7Dy%3D%5Cint_0%5E%5Cinfty%20e%5E%7B-%5Calpha%20r%5E2%7D%5Ccdot2%5Cpi%20r%20%5Cmathrm%7Bd%7Dr%3D%5Cfrac%7B%5Cpi%7D%7B%5Calpha%7D

首先将积分平方化为对整个二维平面的积分,再把积分改写为极坐标,从而严格求出了结果。%5Cpi的出现就是由于取极坐标这一操作。

有了这个积分的值,再代入%5Calpha%3D%5Cfrac%7Bm%7D%7B2kT%7D就不难验证前面的速度分量分布满足归一化了。这里再提一句,通过求得的这个积分式对%5Calpha%20偏导,我们可以求得高斯积分

%20%5Cmathscr%20G_n(%5Calpha)%20%3D%20%5Cint_0%5E%5Cinfty%20x%5E%7Bn-1%7De%5E%7B-%5Calpha%20x%5E2%7D%5Cmathrm%7Bd%7Dx%3D%5Cbeta_n%20%5Calpha%5E%7B-n%2F2%7D%20

其中,系数%5Cbeta需要分奇偶讨论:

%5Cbeta_n%20%3D%20%5Cleft%5C%7B%5Cbegin%7Barray%7D%7Bc%7D%0A%20%20%20%20%5Cfrac%7B%5Csqrt%7B%5Cpi%7D(2k-2)!%7D%7B2%5En%20(k-1)!%7D%2Cn%3D2k-1%5C%5C%0A%20%20%20%20%5Cfrac%7B(k-1)!%7D%7B2%7D%2Cn%3D2k%0A%5Cend%7Barray%7D%5Cright.

记住这个结果,后面要用到它。

我们简单由速度分布推到速率分布:速率为v的分子将分布在速度空间里以v为半径的球表面,所以速率分布函数的形式应该是速度分布函数乘球表面积的形式:(注意是不带方向的速率!)

f(v)%3D4%5Cpi%20v%5E2f(v_x)f(v_y)f(v_z)%3D4%5Cpi%20v%5E2%20(%5Cfrac%7Bm%7D%7B2%5Cpi%20kT%7D)%5E%7B3%2F2%7De%5E%7B-%5Cfrac%7Bmv%5E2%7D%7B2kT%7D%7D

接下来我将告诉你一个比较奇怪的事实:我可以由热力学中的能均分定理和麦克斯韦速率分布推出n维球的表面积公式!

能均分定理是热力学中的一条重要定理,说的是热平衡态的气体每个自由度的平均动能是相同的,每个自由度平均每个分子的动能均为%5Cfrac12kT自由度,大致可以理解为允许自由运动的方向,包括平动与转动:一个粒子,可以横着动,竖着动,前后动,这是3个自由度;如果这是一个分子,还可以在3个不同方向转动,也是3个自由度。

那么考虑n维空间的情况,如果只考虑平动,则在n维空间共有n个自由度,分子平均平动动能

%5Cfrac%7B1%7D%7B2%7Dm%5Cbar%7Bv%5E2%7D%3D%5Cfrac%7Bn%7D%7B2%7DkT

要计算速率分布函数,先考虑若n维球的表面积计算公式为

S_n%3D%20c_n%20r%5E%7Bn-1%7D

其中c_n是待定系数。显然由对称性任意维空间任意方向的速度分量分布都和三维x方向一样,所以n维的气体速率分布

f_n(v)%3Dc_n%20v%5E%7Bn-1%7D(%5Cfrac%7Bm%7D%7B2%5Cpi%20kT%7D)%5E%7Bn%2F2%7De%5E%7B-%5Cfrac%7Bmv%5E2%7D%7B2kT%7D%7D

接下来又根据前面讨论的能均分定理,

%5Cbar%7Bv_n%5E2%7D%3D%5Cfrac%7BnkT%7D%7Bm%7D%3D%5Cint_0%5E%5Cinfty%20v%5E2f_n(v)%5Cmathrm%7Bd%7Dv

再利用前面高斯积分的结果:

%5Cint_0%5E%5Cinfty%20v%5E2f_n(v)%5Cmathrm%7Bd%7Dv%3Dc_n(%5Cfrac%7Bm%7D%7B2%5Cpi%20kT%7D)%5E%7Bn%2F2%7D%5Cmathscr%20G_%7Bn%2B2%7D%20(%5Cfrac%7Bm%7D%7B2kT%7D)%3D%5Cfrac%7Bc_n%7D%7B%5Cpi%5E%7Bn%2F2%7D%7D%5Ccdot%5Cfrac%7B2kT%7D%7Bm%7D%5Ccdot%20%5Cbeta_%7Bn%2B2%7D

得到

c_n%20%3D%5Cfrac%7Bn%5Cpi%5E%7Bn%2F2%7D%7D%7B2%5Cbeta_%7Bn%2B2%7D%7D

所以n维球的表面积公式是

S_n(r)%3D%5Cfrac%7Bn%5Cpi%5E%7Bn%2F2%7D%7D%7B2%5Cbeta_%7Bn%2B2%7D%7Dr%5E%7Bn-1%7D

还记得当初见到这种神奇操作是在物竞书上(好像是郑永令的国培)。当时一直感到疑惑,高维球表面积公式这个理应是纯数学的东西却可以由两条热力学规律导出。现在回过头来思考,其实能均分定理终究是一个定理,它是可以被证明的。而且从另一个角度考虑,其实我们也可以选择不使用能均分定理。直接利用速率分布函数的归一化

%5Cint_0%5E%5Cinfty%20f_n(v)%5Cmathrm%7Bd%7Dv%3D%20c_n(%5Cfrac%7Bm%7D%7B2%5Cpi%20kT%7D)%5E%7Bn%2F2%7D%5Cmathscr%20G_%7Bn%7D%20(%5Cfrac%7Bm%7D%7B2kT%7D)%20%3D%201

从而直接得到

c_n%20%3D%20%5Cfrac%7B%5Cpi%5E%7Bn%2F2%7D%7D%7B%5Cbeta_n%7D

虽然这个结果和上面有一点点不一样,但是回头仔细看一下:

%5Cbeta_n%20%3D%20%5Cleft%5C%7B%5Cbegin%7Barray%7D%7Bc%7D%0A%20%20%20%20%5Cfrac%7B%5Csqrt%7B%5Cpi%7D(2k-2)!%7D%7B2%5En%20(k-1)!%7D%2Cn%3D2k-1%5C%5C%0A%20%20%20%20%5Cfrac%7B(k-1)!%7D%7B2%7D%2Cn%3D2k%0A%5Cend%7Barray%7D%5Cright.

你会发现恰好有%5Cbeta_%7Bn%2B2%7D%3D%5Cfrac%20n2%5Cbeta_n,所以前面两次的结果是完全相等的!球表面积公式应该是

S_n(r)%20%3D%20%5Cfrac%7B%5Cpi%5E%7Bn%2F2%7D%7D%7B%5Cbeta_n%7Dr%5E%7Bn-1%7D

其中,%5Cbeta_n%20%3D%20%5Cint_0%5E%5Cinfty%20t%5E%7Bn-1%7De%5E%7B-t%5E2%7D%5Cmathrm%20d%20t.

写到这里再看,其实能均分定理根本没有起到决定性作用。我有一种感觉,我们在默认各个方向的分子速度分量分布函数具有相同形式时,似乎也就默认以能均分定理作为大前提了。从而在后面步骤中根本没有使用能均分定理的必要。

另外,n维球的表面积公式一定可以通过纯几何导出,这一点是必然正确的,因此我们前面的操作,某种程度上可以说与能均分定理的证明的某些步骤在数学上等价

那么这就是在这个 %5Cpi day我想说的一些与%5Cpi、球和热力学有关的故事。

参考文献

[1] Richard Feynman. The Feynman Lectures on Physics (The New Millenium Edition, Volume I). 上海:上海科技出版社,2013.4. 415~418.

[2] 忘了是哪里看的。

[3] 我自己远古时期的笔记。

高维球与分子速率分布||Pi Day算个Pi的评论 (共 条)

分享到微博请遵守国家法律