1700924042
1700924043
在继续之前,请注意,量子抖动产生的压强,并不依赖于恒星的温度。它只与恒星受挤压的程度有关。如果考虑到温度,电子会因此而“正常”地颤动,并贡献额外的压强;而恒星愈热,它就颤动得越厉害。我们没有讨论来自这一部分的压强,因为时间不够;并且如果真的去计算,会发现它在大得多的量子压强面前相形见绌。
1700924044
1700924045
最后,我们已经准备好将量子压强的表达式代回关键的(1)式,这里值得再写一遍:
1700924046
1700924047
1700924048
1700924049
1700924050
(1)
1700924051
1700924052
但这并不像听起来那么容易,因为还需要知道小立方体上下两面的压强差。可以完全用恒星内部的密度来重写(1)式,而密度本来就是随着恒星内的位置不同而变化的(一定得是这样,否则就不会有压强差了),然后就可以尝试通过解方程来确定密度是如何随着到恒星中心的距离而变化的。这么做就要解一个微分方程(differential equation),而笔者希望避免这种程度的数学。相反,我们要更讲策略,想得更努力(并且计算得更少),以便利用(1)式来推导出白矮星的质量和半径之间的关系。
1700924053
1700924054
显然,这个小立方体的尺寸和它在恒星中的位置都是完全任意的;我们要得出的关于整体恒星的结论不能依赖于小立方体的细节。我们先来做一些可能看似毫无意义的事情。完全可以用恒星的尺寸来表示立方体的位置和大小。如果R是恒星的半径,则可以将立方体到恒星中心的距离写作r=aR[281],其中a是一个介于0和1之间的无量纲(量纲是dimension)数。所谓无量纲,就是说它是纯粹的数,不含有单位。如果a=1,立方体就在恒星的表面;而如果a=1/2,则它在中心到表面一半的位置。类似地,可以用恒星的半径来表示小立方体的尺寸。如果L是立方体的边长,则可以写成L=bR,其中b又是一个纯数;如果希望立方体相对恒星很小,则b也会很小。这里面绝对没有什么深奥的东西;在这个阶段,一切都很明显,乃至看似没有意义。唯一值得注意的是,R是一个十分自然的长度单位,因为没有其他任何跟白矮星有关的长度可以合理地替代它。
1700924055
1700924056
我们可以继续用类似的方法困扰自己,用恒星的平均密度来表示立方体所处位置的恒星密度,即写成ρ=fρ,其中f还是一个纯数,而是恒星的平均密度。如前所述,立方体的密度取决于它在恒星内部的位置:如果更接近中心,它的密度就会更大。既然ρ平均密度与小立方体的位置无关,那么f就得与之相关,即f得由距离r所决定,这显然意味着它取决于乘积aR。现在这里是决定剩余计算的关键信息:f是一个纯数,但R不是(因为它度量距离)。这意味着f只能取决于a,而和R毫无关系。这是一个非常重要的结论,因为它告诉我们,白矮星的密度分布曲线是“尺度不变”(scale invariant)的。这就是说,无论恒星的半径是多少,它的密度都会以相同的方式随半径而变化。例如,在距离中心3/4的地方,无论恒星尺寸如何,那里的密度与白矮星平均密度的比值对于所有白矮星都一样。有两种方法可以看出这个结果的关键性,笔者认为应该把它们都展示出来。我们中的一人是这样说的:“这是因为,任何函数,如果是无量纲的,它的变量就必须也是无量纲的;而对于函数f,它依赖于有量纲的r,则能作为其变量的、唯一的无量纲组合是r/R=a,因为R是我们所知唯一具有长度量纲的量。”
1700924057
1700924058
另一位笔者觉得下面的话更清楚:“一般而言,f可以由r以复杂的方式决定,而后者即小立方体与恒星中心距离。但为了本段的说明,我们假设它们是简单的成正比,即f∝r。换言之,f=Br,其中B是一个常数。这里的关键是,我们希望f是一个纯数,而r是以(比如说)米为单位。这意味着B必须以1/米为单位,这样长度单位才能互相抵消。那么从结果来看,B会是什么呢?我们不能随意选一些东西,比如‘1米的倒数’,因为这没有意义,与恒星毫无关系。比如说,为什么不选1光年的倒数,并且得到截然不同的正比关系呢?我们手头唯一的长度是R,即恒星的物理半径,所以我们被迫用它来确保f永远是一个纯数。这意味着f只取决于r/R。应该可以看到,如果我们开始时的假设是,比如说f∝r2,也会得出相同的结论。”这些其实和前一位笔者说的一样,只是更长。
1700924059
1700924060
这个结论意味着,可以将尺寸为L、体积为L3的小立方体的质量,表示为Mcube=f(a)L3ρ。我们把f写成f(a),是为了提醒你,f实际上只取决于我们对a=r/R的选择,而与恒星大规模的性质无关。用相同的论证还可以写下Min=g(a)M这个式子,其中g(a)也只是a的函数。例如,函数g(a)在a=1/2时的值就会告诉我们,具有一半恒星半径的同心球体,其质量占恒星质量的比例;并且它对于所有白矮星都是一样的,与具体的白矮星半径无关,这在前一段已经论证过了[282]。你可能已经注意到,我们在稳步处理(1)式中出现的各种符号,用无量纲的量(a、b、f和g)乘以只和恒星质量及半径有关的量(恒星的平均密度也由M和R决定,因为ρ=M/V,而根据球体积公式,V=4πR3/3)来代替它们。要完成任务,只需对压强差进行相同的处理;可以写作Pbutton-Ptop=h(a,b)κρ5/3,其中h(a,b)是一个无量纲的量。h(a,b)由a和b共同决定,是因为压强差不仅取决于小立方体的位置(由a表示),还取决于立方体的大小(由b表示):较大的立方体,对应的压强差也较大。关键是,和f(a)与g(a)类似,h(a,b)也不能单单与恒星的半径有关。
1700924061
1700924062
可以用刚刚推导出的表达式,重新写出(1)式:
1700924063
1700924064
1700924065
(hκρ5/3)×(b2R2)=
1700924066
1700924067
这看起来有点乱,不太像是在一页内就能中大奖。关键之处在于,这是在表达恒星的质量与半径的关系;两者间的具体关系已经近在眼前(或者远在天边,这取决于你的数学能力)。在代入恒星的平均密度[即ρ=M/(4πR3/3)]之后,这个乱七八糟的式子可以整理成
1700924068
1700924069
RM1/3=κ/(λG)
1700924070
1700924071
(5)
1700924072
1700924073
其中,
1700924074
1700924075
1700924076
1700924077
1700924078
现在的λ只依赖于无量纲的量a、b、f、g和h,这意味着它不由描述恒星整体的量M和R所决定,因此这个式子对于所有的白矮星都是一样的。
1700924079
1700924080
如果你担心,当改变a和/或b(这意味着改变小立方体的位置和/或尺寸)时会怎么样,那么你肯定没有领会到这段论证的威力。从表面上看,改变a和b当然会改变λ,这样我们就会得到RM1/3的不同结果。但这是行不通的,因为我们知道,RM1/3是取决于恒星整体的,而不是其内部一个、我们可能关心也可能不关心的小立方体的具体性质。这意味着a和b的任何变化,都必须由f、g和h的相应变化所完全补偿,使得λ保持不变。
1700924081
1700924082
(5)式相当明确地指出,白矮星可以存在;这是因为我们已经成功地解出了引力-压强的平衡方程[(1)式]。这不是一件平凡的事——因为有可能会发现,平衡方程(1)对于M和R的任何组合都无法被满足。(5)式还作出预言,RM1/3一定是常数。换句话说,如果我们仰望星空,测量白矮星的半径和质量,就应该会发现,每颗白矮星的半径乘以其质量的立方根,都会得出相同的数。这实在是一个大胆的预言。
1700924083
1700924084
笔者刚才展示的论证还可以改进,因为可以准确地计算出λ的数值,但这需要解出一个关于密度的二阶微分方程,而所需的数学技巧对于本书太遥远了。请记住,λ是一个纯数;“它就是它”,可以用一点更高级的数学去算出来。我们没有在这里把它算出来,但这对我们的成就没有丝毫影响:我们已经证明,白矮星可以存在,并且对于其质量与半径的关系做出了预言。在计算出λ之后(可以在家用电脑上完成),并代入κ和G的值,预言的结果是:
1700924085
1700924086
RM1/3=(3.5×1017kg1/3m)×(Z/A)5/3
1700924087
1700924088
对于纯氦、碳或氧(Z/A=1/2)组成的核心,等于1.1×1017kg1/3m。对于铁核心,Z/A=26/56,式子中的1.1略微变小到1.0。笔者翻阅了学术文献,收集了散布在我们的银河系中的16颗白矮星的质量和半径。对于每一颗,我们都计算出了RM1/3的值,结果是天文观测给出RM1/3≈0.9×1017kg1/3m。令人激动的是观测和理论比较一致;我们成功地利用泡利不相容原理、海森伯不确定性原理和牛顿引力定律,预言了白矮星的质量-半径关系。
1700924089
1700924090
当然,这些数字有一些不确定性(理论值为1.0或1.1,观测值等于0.9)。如果是正确的科学分析方法,现在就会开始讨论理论和实验一致的可能性有多大,但对我们的目的而言,这种层次的分析是不必要的,因为一致得惊人了。我们得以算出这一切,并且误差率仅有约10%,是十分奇妙的。这也是有力的证据,证明我们对恒星和量子力学有了相当的了解。
1700924091
[
上一页 ]
[ :1.700924042e+09 ]
[
下一页 ]