作者：李文龙

摘要： 星系间的相互作用是驱动宇宙大尺度结构演化的常见现象。本文以NASA在2026年4月10日APOD上公布的“触角星系”（NGC 4038/4039）图像为切入点，梳理了星系交汇过程中的物理机制。动力学分析表明，在星系碰撞期间，恒星系统表现为典型的“无碰撞流体”；而大面积分布的分子气体云则因物理挤压产生强激波，进而触发局部的星暴活动。此外，非对称的引力潮汐力剥离了星系外围物质，形成大尺度潮汐尾。为了直观展示这一天体力学过程，本文利用前端Web技术搭建了一个受限N体模拟沙盘，加入引力平滑算法，避免星体近距离引力计算失常，平稳还原星系引力运动过程，并通过渲染帧叠加机制复现了潮汐尾的动态演化，为天体物理教学提供了一种低门槛的可视化方案。

关键词： 星系碰撞；引力潮汐；星暴；受限N体模拟；交互式可视化；NGC 4038

一、引言

在平时看科普新闻的时候，大众往往认为星系是静止且孤立的，但实际在天体观测中，星系之间的引力拖拽和碰撞非常普遍。2026年4月10日，NASA的每日天文一图（APOD）发布了一张“触角星系”（NGC 4038和NGC 4039）的高清图像[1]。照片里这两个大型星系正处于剧烈的交汇阶段。

很多人一听到“碰撞”，第一反应就是两块固体直接撞碎。其实不然，星系级别的碰撞涉及复杂的流体力学与轨道动力学，实际过程要反直觉得多。为了把照片背后的物理机制讲清楚，我在这篇报告里梳理了相关的动力学特征，顺便分享一个我用JavaScript写的网页沙盘，带大家直观地推演一下星系并合时的演化过程。

二、星系碰撞的物理机制

要弄懂这张照片，我们不能只看表面，得把星系拆解成恒星、星际气体和宏观引力场三个层面来分析。

2.1 形貌特征与“触角”的由来

先看照片中间那团高亮且扭曲的区域，这是两个星系NGC 4038和NGC 4039的核心，它们目前距离地球大约6000万光年。顺着交汇区往外围看，视野里有两条往两边延伸得很长的微弱光带。这种结构在天文上叫“潮汐尾”。因为它的整体轮廓长得特别像昆虫头顶上的两根长触须，所以这对星系（编号Arp 244）在观测圈里一般被俗称为“触角星系”^[1]。

   

图 1：触角星系（NGC 4038/4039）的光学形貌。图中清晰展示了中心交汇区与向外延伸的微弱引力潮汐尾（图片来源：NASA/APOD [1]）。

2.2 恒星的“无碰撞”穿透特性

刚才提到星系碰撞不是两块石头硬碰硬，这主要是因为太空实在是太宽广了。即使两个星系相互穿透，里面的恒星（比如像太阳这种主序星）也极少发生真实的物理相撞。

我们可以简单从尺度方面来看：恒星的自身物理半径比起它们之间的平均距离（通常在几光年左右）实在不值一提。这就好比把两把沙子分别撒在太平洋两岸，指望它们在海中心刚好撞上，概率基本为零。所以，在长达数亿年的交汇期里，这两个星系里的恒星群主要受整个系统的宏观引力势阱支配，互相穿插而过。在轨道动力学模型中，我们通常把这种现象处理为“无碰撞流体”（Collisionless fluid）^[2]。

2.3 分子气体相撞与星暴现象

既然恒星碰不上，照片中间为什么会爆发出那么亮的光芒？这其实是星系里大面积分布的冷分子气体和尘埃云造成的。

和散落的恒星不同，气体云体积庞大，属于连续分布的星际介质。当星系交汇时，恒星能躲开，但这些气体云是实打实地会撞在一起的。气流的对撞在局部产生了极高的压力和强烈的激波（Shockwaves）^[3]。这种剧烈的流体挤压瞬间打破了气体云原本的流体静力学平衡，导致它们迅速发生引力坍缩。根据天体物理中的肯尼卡特-施密特定律（Kennicutt-Schmidt Law），气体面密度的急剧上升会直接触发极高效率的恒星形成活动^[4]。这就是大家常听到的“星暴”（Starburst）。照片中心那些蓝粉交织的明亮斑块，正是这么批量出生的年轻超级星团^[5]。

                           

图 2：冷气体面密度与恒星形成率的经验关系图（Kennicutt-Schmidt Law）。星际气体的物理挤压是星系交汇中触发极端星暴活动的核心机制（图片来源：参考文献 [4]）。

2.4 引力潮汐与物质剥离

关于那两条超过5万光年长的“触角”是怎么被拉扯出来的。这完全是空间中的引力梯度导致的。

当两个大质量星系互相靠近并绕转时，引力场在空间上的分布是不均匀的。星系靠近伴星系的那一侧受到的引力，远大于背对侧受到的引力。这种空间上的引力差，也就是“引力潮汐力”。星系外围的恒星和气体原本受到的自引力束缚就比较弱，在潮汐力拉扯和轨道离心力的共同作用下，它们直接被剥离出了星系主盘^[6]。随着核心的持续绕转，这些被剥离的物质被抛射到太空深处，慢慢拉伸成了我们看到的巨大弧形潮汐尾^[1]。

三、可视化模拟器的参数设计

前半部分主要从物理的角度推导潮汐尾的成因，接下来我将会从计算机的角度带大家直观感受两星系相撞的过程。为了验证潮汐尾的成因，我在本地用HTML5 Canvas绘图技术，以JavaScript为程序主体绘制了一个简单的二维引力沙盘。界面去掉了多余的UI，采取轻量化的界面风格，只留了两个控制轨道初始条件的滑块：

1.接近速度（Velocity）： 决定了系统的初始平动动能。动能太大，星系交汇一次就飞走了；动能降下来，它们才会被彼此的引力捕获进入轨道锁定。

2.碰撞角度（Offset）： 天体力学里叫碰撞参数（Impact parameter）。如果把角度调成完全的正面相撞，系统很快就会混合成一团无规则的椭球体；但如果稍微调大点，模拟“擦肩而过”的掠射状态，就能非常清楚地看到非对称的引力梯度是怎么把潮汐尾甩出来的。

                           

图 3：本文自行开发的 Web 交互式星系引力沙盘界面截图。图中展示了在“掠射”初始参数下，模拟器成功重现的引力潮汐尾演化现象（来源：作者自制）。

四、底层物理算法与渲染逻辑

由于是在普通浏览器里跑代码，我没有调用庞大的第三方物理引擎，底层全是基于牛顿力学公式。为了保证程序既不卡顿又不违背物理常理，代码里做了几个关键的数值处理。

4.1 受限N体问题近似

真实的星系有上千亿颗恒星，要是算完整的N体引力交互，时间复杂度是 O(N^2)，前端浏览器的算力根本撑不住。既然前面在2.2节论证了恒星可以看作无碰撞流体，我在代码里就借用了经典的“受限N体模型”（Restricted N-body model）^[7]。

程序区分了大质量的星系核心（Core）和无质量的测试粒子（Particle，代表恒星）。粒子只受两个核心的引力吸引，粒子之间互不干涉。这样一来，计算复杂度直接降到了O(N)，不仅保证了60帧的流畅运行，在物理逻辑上也完全站得住脚。

4.2 引力软化处理

写过引力数值积分代码的同学都知道一个通病：万有引力公式里，距离r在分母上。如果积分步长内有粒子太靠近核心，r趋近于0，算出来的引力会无穷大，直接导致数值溢出，粒子瞬间飞出屏幕。

为了解决这个“奇点”Bug，我在更新粒子位置的代码里加了个限制条件：let distanceSq = Math.max(force.magSq(), 150);。通过强制给距离的平方卡一个150的下限，人为避免了分母趋于0的情况。在天体物理计算中，这叫作“引力软化”（Gravitational Softening）^[8]，属于保障程序数值稳定性的常规操作。

4.3 视觉暂留与拖影渲染

最后是视觉层面的处理。为了让几千个离散的几何粒子跑出星际气体那种连绵不断的流体质感，我利用了Canvas渲染的一个小技巧。

在动画帧循环里，我没有用 clearRect 把上一帧全擦掉，而是每次给整块画布蒙一层透明度15%的暗色背景（rgba(3, 3, 8, 0.15)）。靠着简单的视觉暂留效应，粒子运动时自然就拖拽出了渐暗的尾迹。没用复杂的流体力学方程，仅靠前端渲染技巧就在屏幕上还原出了潮汐尾的动态张力。

五、结论

结合APOD的光学观测和程序跑出来的数值结果，我们可以很清楚地看到：NGC 4038/4039的碰撞是一场复杂的宏观动力学演化。恒星群的无碰撞穿透、分子气体的受压星暴，以及引力潮汐对星系盘的外围剥离，共同塑造了其现有的宏观形态。通过受限N体沙盘，我们不仅在数值层面验证了潮汐尾的力学成因，也为天体物理理论的科普展示提供了一个低门槛、互动性强的实用工具。

参考文献

[1]NASA. Astronomy Picture of the Day: The Antennae Galaxies[EB/OL]. (2026-04-10)[2026-05-19]. https://apod.nasa.gov/apod/ap260410.html.

[2] STRUCK C. Galaxy collisions[J]. Physics Reports, 1999, 321(1-3): 1-137.

[3] UEDA J, WATANABE Y, IONO D, et al. ALMA observations of the dense and shocked gas in the nuclear region of NGC 4038 (Antennae galaxies) [J]. Publications of the Astronomical Society of Japan, 2017, 69(1): 6.

[4] KENNICUTT R C. Star formation in galaxies along the Hubble sequence[J]. Annual Review of Astronomy and Astrophysics, 1998, 36(1): 189-231.

[5] WHITMORE B C, ZHANG Q, LEITHERER C, et al. The star clusters in the Antennae galaxies: observations with the Hubble Space Telescope[J]. The Astronomical Journal, 1999, 118(4): 1551-1576.

[6] BARNES J E, HERNQUIST L. Dynamics of interacting galaxies[J]. Annual Review of Astronomy and Astrophysics, 1992, 30(1): 705-742.

[7] TOOMRE A, TOOMRE J. Galactic bridges and tails[J]. The Astrophysical Journal, 1972, 178: 623-666.

[8] AARSETH S J. Gravitational N-body simulations: tools and algorithms[M]. Cambridge: Cambridge University Press, 2003.

作者：李文龙，华中科技大学基础医学院本科生，儿科学2501班

编校：雷卫华

文章编号：华中大天文260526A