跳到主要内容
Sandbox Physics
分子云的引力与磁支撑

分子云 · L2

磁场支撑为什么要看方向

气体沿磁场方向聚集时,可以少挤压磁场;横过磁场聚集时,会把磁场线挤紧,磁力便参与抵抗。磁场有多强还不够,还要说明压缩沿哪个方向。把它统一加成一个“更大的声速”,会丢掉这个区别。

沿着场走,还是把场挤紧

把均匀磁场想成一束有方向的线。每个位置都有一个磁场矢量,线把这些方向连起来,方便想象。理想磁流体近似让磁通随气体运动,横向压缩会改变场强,弯曲也会带来恢复作用。

挤紧磁场后,向两侧推回的作用称为磁压。磁压差可以把磁场较强的区域推向较弱的区域;沿弯曲磁场的张力则试图把它拉直。两种作用合起来,才是磁场对气体的力。沿场压缩和横向压缩,触发它们的方式不同。

先选一段密度波,再旋转背景磁场,比较同一个波向上的响应。这里的角度始终指背景磁场与选定波向的夹角。整片云还可以允许其他方向的运动。

对应阅读:Chandrasekhar & Fermi 1953, Problems of Gravitational Stability in the Presence of a Magnetic Field · Molecular-cloud Jeans standalone scientific reproduction

把六个扰动量放在同一套坐标里

令波向量k=kx^\mathbf k=k\hat{\mathbf x},背景磁场位于x,yx,y平面,与xx轴夹角θ\theta。用国际单位制的磁场除以μ0ρ0\sqrt{\mu_0\rho_0},结果带速度单位;μ0\mu_0是真空磁导率,ρ0\rho_0是背景体密度。归一化场强为Alfvén速度vAv_A。

背景分量为Bx=vAcos⁡θB_x=v_A\cos\theta和By=vAsin⁡θB_y=v_A\sin\theta。这里的Bx,ByB_x,B_y已经用速度单位表示;扰动by,bzb_y,b_z也作同样归一化。密度差δ\delta与三向速度vx,vy,vzv_x,v_y,v_z合起来,共有六个随时间变化的量。

磁场没有源头,扰动满足∇⋅δB=0\nabla\cdot\delta\mathbf B=0。对沿xx变化的一个非零波数,便有bx=0b_x=0。这并没有把恒定的背景BxB_x设成零。

用eikxe^{ikx}记空间变化,最后取实部;i2=−1i^2=-1。乘一个ii相当于让正弦、余弦换一个相位。这种写法让六个变量都能存成复数幅度,同时保留同一幅空间图。

对应阅读:Chandrasekhar & Fermi 1953, Problems of Gravitational Stability in the Presence of a Magnetic Field · Molecular-cloud Jeans standalone scientific reproduction

磁力和感应一起决定运动

记ωJ2=4πGρ0\omega_J^2=4\pi G\rho_0。连续性仍把密度变化连到纵向流动;纵向动量里,压力和引力之外多了由横向磁场扰动带来的力。横向速度又会受背景纵向磁场牵动。

气体运动通过感应方程改变磁扰动,改变后的磁场再推气体。把两组方程一起积分,密度图、速度箭头和磁场方向便能对应同一个过程。

纵向磁项的负号也有直观来源。归一化后,磁压除以背景密度的一阶变化是BybyB_yb_y;沿xx取负梯度,就得到−ikByby-ikB_yb_y。磁场较挤的一侧推力较强,因而倾向把这段横向磁场推散。背景纵向场另外连接横向张力,这正是其余分量方程保留的作用。

先用一个很小的纵向位移ξx\xi_x准备初态。某段被压短时,连续性给δ=−∂xξx\delta=-\partial_x\xi_x;同一位移挤紧横向磁通,给by=−By∂xξxb_y=-B_y\partial_x\xi_x,因此初态满足by=Byδb_y=B_y\delta。斜向演化随后还会激起横向速度,要继续用完整感应方程,不能永远把这一初始关系当作代数替代。

δ˙=−ikvx,v˙x=−ik(cs2−ωJ2k2)δ−ikByby\dot\delta=-ikv_x,\qquad \dot v_x=-ik\left(c_s^2-\frac{\omega_J^2}{k^2}\right)\delta-ikB_yb_y
密度与纵向速度相连;引力项已经用Poisson方程消去势。
v˙y=ikBxby,v˙z=ikBxbz\dot v_y=ikB_xb_y,\qquad \dot v_z=ikB_xb_z
横向磁力依赖背景磁场沿波向的分量。
b˙y=ik(Bxvy−Byvx),b˙z=ikBxvz\dot b_y=ik(B_xv_y-B_yv_x),\qquad \dot b_z=ikB_xv_z
感应方程让磁扰动随速度变化,和上面的动量方程组成六变量系统。

对应阅读:Chandrasekhar & Fermi 1953, Problems of Gravitational Stability in the Presence of a Magnetic Field

压缩模和横向Alfvén模要分别看

设每个变量随时间按e−iωte^{-i\omega t}变化,消去各分量,含密度压缩的频率满足下面的四次式。它把两个ω2\omega^2分支放在一起;负的频率平方对应增长与衰减。

还有一对Alfvén模。它主要是横向速度和磁扰动来回交换,频率取决于磁场沿波向的分量。把这一对遗漏,六变量矩阵就会少两个本征值。

在θ=0\theta=0时,压缩分支恢复L1的ω2=cs2k2−ωJ2\omega^2=c_s^2k^2-\omega_J^2。磁场再强,沿场这段压缩仍保留原Jeans门槛。对精确横向的θ=90∘\theta=90^\circ,非零压缩频率则含cs2+vA2c_s^2+v_A^2,横向磁压参与恢复。

把场强逐渐降到零,也是一个可检查的极限。四次式分出ω2=0\omega^2=0与普通Jeans分支,Alfvén频率也降为零。横向速度可以保持常数,密度压缩则回到没有磁场的L1方程。这个极限能检查代码有没有把场强额外加错一次。

ω4−[(cs2+vA2)k2−ωJ2]ω2+vA2k2cos⁡2θ(cs2k2−ωJ2)=0\omega^4-\left[(c_s^2+v_A^2)k^2-\omega_J^2\right]\omega^2+v_A^2k^2\cos^2\theta(c_s^2k^2-\omega_J^2)=0
压缩模色散关系。角度明确出现,不能对所有方向统一换成一个标量声速。
ωA2=k2vA2cos⁡2θ\omega_A^2=k^2v_A^2\cos^2\theta
另行存在的Alfvén因子,与压缩四次式合起来给出全部六个时间本征值。

对应阅读:Chandrasekhar & Fermi 1953, Problems of Gravitational Stability in the Presence of a Magnetic Field · Molecular-cloud Jeans standalone scientific reproduction

算一组强磁场下的不同方向

沿用J=1.5J=1.5和整盒基波,取vA/cs=4v_A/c_s=4。把四次式当作关于ω2\omega^2的二次方程求根:沿场增长率为0.745356/t00.745356/t_0,夹角45∘45^\circ时降为0.537226/t00.537226/t_0,到89∘89^\circ时只有0.0135481/t00.0135481/t_0。

最后一个数很小,但仍为正。接近横向时,一段有限时长的记录可能看起来几乎没有慢增长。精确90∘90^\circ端点含中性分支和振荡分支;浮点计算留下的极小正数,需要结合数值尺度分类。

画面中的初态来自一个纵向位移:by=Byδb_y=B_y\delta,速度起初为零,密度幅度取0.020.02。它通常是多个模的组合,因此不能把整条密度曲线都当成上面某一个纯指数。保存的六分量状态决定了实际叠加。

三个场强的增长模速率随方向改变,沿场增长率保持相同。
按本章色散关系计算,Jeans比固定为1.5。强横向场能阻止这段压缩,沿场增长仍然存在。图中给出单个本征模的增长率;Lab从混合初态起步,密度不必沿一条纯指数曲线变化。

对应阅读:Molecular-cloud Jeans standalone scientific reproduction

保持线性,需要检查整个扰动

只看密度小不小还不够。密度可能暂时很小,速度或磁扰动却已积累起来。示例同时限制密度差、相对场扰动和归一化速度,达到指定小扰动上限前停止。

速度用cs2+vA2\sqrt{c_s^2+v_A^2}作归一化尺度,相对扰动上限取0.150.15。这个速度尺度用于判断何时停止;传播和增长仍由六个方程计算。它没有替代带方向的磁力。

非磁气体片段解的是有限幅度的质量与动量方程;这一章解的是均匀背景上的线性磁扰动。两组方程分别适合有限幅度气体和小幅磁扰动,选择时先看自己要解释的运动。

对应阅读:Molecular-cloud Jeans standalone scientific reproduction

方向支撑与整片云的磁临界性

要问整片云能否被磁场托住,还得指定云的形状、磁通穿过哪张面、外界怎样固定磁场。文献中的质量磁通临界值带着这些几何条件。单个周期密度波没有提供一个通用的整云临界质量。

在这里的均匀模型里,只要允许选择沿场的波向,就仍能找到原本Jeans不稳定的压缩方向。后面的论文导读再看三维湍流计算:它既有不同方向的流,也有另一套质量与分辨率问题。

对应阅读:Guszejnov et al. 2020, Can magnetized turbulence set the mass scale of stars? arXiv:2002.01421v1 · Federrath & Klessen 2012, arXiv:1209.2856v1

练习与解释性答案

  1. 把磁场增强到声速对应尺度的四倍,能保证所有方向的Jeans压缩都稳定吗?

    展开解释

    不能作这个保证。沿场分支仍满足ω2=cs2k2−ωJ2\omega^2=c_s^2k^2-\omega_J^2,因为这段压缩没有挤紧横向磁场。必须给定波向;强磁场能显著改变横向和斜向响应。

  2. 89∘89^\circ的有限记录几乎平坦,为什么还不能把它归入精确横向的中性端点?

    展开解释

    对应慢分支的增长率仍约0.0135481/t00.0135481/t_0。有限观测时间内增长可很慢,且初值含多个模。先求本征值并声明数值中性容差,再解释时间曲线。精确端点的背景纵向磁场要设为零,不能由浮点余弦残差决定物理分类。

  3. 对沿xx变化的扰动,bx=0b_x=0是否意味着背景磁场不能沿xx?

    展开解释

    没有这个含义。无散条件限制的是这个非零波数的磁场扰动。恒定的背景BxB_x可以非零,还会通过磁张力和感应方程连接横向速度与磁扰动。

动手检验预测

看压力与引力怎样改变同一团气体

  1. 进入线性磁分支,依次选沿场、斜向与横向响应,保持Jeans比和场强比相同;先预测哪种压缩受到的恢复最强。
  2. 拖场方向柄,比较密度、速度与磁场箭头。角度改变的是背景方向,新的初态随后按同一组六变量方程演化。
  3. 保存一条记录,比较实际密度时间史与模式增长率。说明为什么混合初值不必给出一条纯指数。

来源与计算范围

  1. Chandrasekhar & Fermi 1953, Problems of Gravitational Stability in the Presence of a Magnetic Field

    第V节,第138至141页,式158至173。原文采用绝热声速导数和高斯磁单位;章节另行说明等温导数与国际单位归一化。链接原文,不重发受限许可论文。

  2. Federrath & Klessen 2012, arXiv:1209.2856v1

    固定v1,第2.2/2.4、3.1至3.3与7节:三维等温受驱磁流体、吸积粒子判据,以及统计磁压近似的适用范围。气体片段计算没有加入这些驱动、粒子或数值发电机条件。

  3. Guszejnov et al. 2020, Can magnetized turbulence set the mass scale of stars? arXiv:2002.01421v1

    固定v1,第2至4节、图8至10、式8至11/17/22至23。分别阅读质量加权百分位与按颗数分布。图8模型名配对和近似量纲例子按上下文保留,没有自行替作者勘误。

  4. Molecular-cloud Jeans standalone scientific reproduction

    独立有限体积与过采样伪谱计算、Jeans解析参照和完整线性磁流体参照,含十个支持病例、固定长时病例、官方小输入及原始失败。安装普通Python依赖后,默认重放会离线实际重算科学步骤。

反馈这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策