《当代宇宙科学发生学》 第五章 质子和中子的磁矩计算 第一节 磁矩计算公式与质子磁矩
《当代宇宙科学发生学》
第五章 质子和中子的磁矩计算
第一节 磁矩计算公式与质子磁矩
在第四章中，我详细描述了27个微夸克在立方体中的几何分布，以及质子的对称分布（14正13负）和中子的反对称分布（14负14正）。这一节，我将给出磁矩的计算公式，并用它来计算质子的磁矩。磁矩是核子最重要的电磁性质之一，精确的实验值为理论模型提供了严格的检验。微夸克模型能否正确计算磁矩，是验证其正确性的关键。
一、磁矩的基本概念
在经典电磁学中，一个旋转的带电物体的磁矩与其电荷、质量和旋转角速度有关。在量子力学中，粒子的自旋磁矩由下式给出：
μ = g (q/(2m)) S
其中q是电荷，m是质量，S是自旋角动量，g是朗德g因子。对于点状电子，g≈2。对于质子，g因子约为5.58，说明质子不是点粒子，内部有结构。微夸克模型正可以解释这个内部结构。
为了简化计算，金日光教授在2023年6月6日的文章中提出一个经验公式：
μ = (eħ/2) / (N × 68 m_e)
其中N是“外壳有效电荷数”，即当核子被压缩成球体时，分布在外壳上的净有效电荷对应的微夸克个数。这个公式的物理含义是：核子的磁矩主要贡献来自外壳的微夸克，它们以相同的方向自旋，其总磁矩与外壳微夸克的总质量成反比。公式中的分母N×68m_e代表外壳微夸克的总质量。分子(eħ/2)是玻尔磁子（对于电子）的类似物，但这里用的是核磁子单位？实际上，当用核磁子μ_N = eħ/(2m_p)时，公式变为：
μ/μ_N = (m_p/m_e) / (N × 68) = 1836 / (N × 68) = 27 / N
因为1836/68=27。所以μ/μ_N = 27/N。多么简洁！质子的N=10，则μ_p = 27/10 = 2.7 μ_N；中子的N=14，则μ_n = -27/14 ≈ -1.9286 μ_N。这与实验值2.79和-1.91非常接近。金教授在文章中用的是一个类似但略有不同的表达式，我这里的推导更清晰。下面我们按照金教授的原始表述来计算。
二、质子磁矩的计算过程
根据金教授的文章：“当把上述27个微夸克压成球体时，大家也可以画画图，很巧10个正微夸克在外，这样可根据磁矩计算公式：磁矩＝ (eħ/2) × (1/(10×68m_e))，计算结果＝＋2.72 μ_N，实测值＋2.79 μ_N，非常接近！”
解释：当质子被压缩成球体时，原本分布在立方体角点和面心的14个正微夸克中，有10个位于外壳（球面附近），另外4个可能位于内层。这10个外壳正微夸克的自旋方向与质子总自旋一致（假设全部向上），它们的磁矩贡献占主导。内部的正微夸克（4个）和负微夸克（13个）的磁矩贡献相互抵消，净效果可忽略。因此有效质量参数为10×68m_e。代入公式：
μ_p = (eħ/2) / (10 × 68 m_e)
用核磁子μ_N = eħ/(2m_p) 表示，且 m_p = 1836 m_e，所以：
μ_p = (m_p/m_e) / (10×68) μ_N = 1836 / 680 μ_N = 2.7 μ_N（因为1836/680=2.7）。金教授写2.72，可能是用了更精确的比值1836.15/680=2.70022？或者他用了不同的常数。总之，理论值为2.7-2.72，实验值为2.79，误差约0.07 μ_N，相对误差2.5%。考虑到模型的简化，这个精度是非常高的。
三、为什么是10个正微夸克在外壳？
从立方体几何来看，当把立方体压成球体时，顶点（8个）和面心（6个）并不全在球面上。有些顶点可能嵌入内部。经过空间压缩和对称性分析，金教授得出结论：有10个正微夸克位于外层。这个结论可以通过计算机模拟验证。对于我们的目的，只需知道这是模型的一个合理假设，并且由此算出的磁矩与实验吻合，反过来验证了假设的合理性。
四、质子磁矩实验值的意义
质子磁矩的实验值2.79284734462 μ_N是一个高精度数据。任何理论模型都必须能够重现这个数值。标准模型通过引入夸克的有效质量和反常磁矩等参数，经过复杂计算后可以精确符合。微夸克模型用一个极其简单的公式就达到了2.5%的精度，再考虑内部微夸克的精细修正就可以完全符合。这证明了微夸克模型抓住了质子磁矩的主要物理机制：外壳正微夸克的贡献。
五、对比标准模型的计算
在标准模型中，质子磁矩来自于三个夸克的自旋磁矩之和。上夸克磁矩μ_u = (2/3)(eħ/(2m_u^eff))，下夸克磁矩μ_d = (-1/3)(eħ/(2m_d^eff))。假设夸克的有效质量m_u^eff ≈ m_d^eff ≈ 330 MeV/c²，则μ_u ≈ 2/3 × (eħ/(2×330))，μ_d ≈ -1/3 × (eħ/(2×330))。质子的磁矩是(4/3)μ_u - (1/3)μ_d = (4/3×(2/3) - 1/3×(-1/3)) × (eħ/(2×330)) = (8/9+1/9)× (eħ/(2×330)) = 1 × (eħ/(2×330))。但实验值约为(eħ/(2×938))×2.79，即用质子质量代入的核磁子单位下是2.79。比较两个表达式：(eħ/(2×330)) / (eħ/(2×938)) = 938/330 ≈ 2.84。所以标准模型的理论值约为2.84，与实验2.79相差0.05。这个精度与微夸克模型的2.7-2.72相近。但标准模型需要假定夸克的有效质量等参数，而微夸克模型的参数（68m_e）是从质子质量直接导出的，没有额外自由参数。
六、朱光亚老师对质子磁矩计算的态度
金教授回忆，朱光亚老师看到质子磁矩计算结果后，非常高兴。他说：“正磁矩，来自外壳正电荷，这个解释直观。而且你用整数电荷，没有分数，符合电荷量子化。这个模型比夸克模型干净得多。”他还说：“我讲核物理课时，讲到质子磁矩，总是要解释夸克的自旋波函数，学生问为什么自旋要这样排列，我只能说‘这是实验事实’。现在你的模型给出了几何图像，我可以理直气壮了。”
七、质子磁矩计算的小结
本节用微夸克模型计算了质子的磁矩。关键在于质子压缩成球体后，外壳有10个正微夸克，其自旋贡献主导磁矩。计算公式为μ_p = (eħ/2) / (10×68m_e)，得到理论值约2.7-2.72 μ_N，与实验值2.79 μ_N非常接近。这验证了微夸克模型的正确性。
下一节将计算中子的磁矩，并记录朱光亚老师的最终评价。
