数值参数对潮汐触发条形形成的影响

我们研究了诸如粒子数 N、引力软化长度 $\epsilon$ 和时间步长 $\Delta t$ 等数值参数的联合影响，探讨其在星系模拟中的作用。在孤立星系模型中，我们进行了收敛性研究，估算了哪些数值参数范围内的模型不会显著偏离初始配置。固定 N 后，我们计算了沿星系盘半径的平均粒子间距 $\lambda(r)$ 的范围。我们发现，在 N=1310720 粒子的模拟中，$\lambda$ 变化约为 6 倍，而相应的最终托姆热参数 Q 仅改变约 5%。降低 N 后，$\lambda$ 和 Q 的范围会扩大。较大的 $\epsilon$ 和较小的 N 导致更早的条形形成。在给定参数集合下，采用普勒软化方法时的盘加热小于采用样条软化方法。在星系碰撞模型中，数值模拟表明形成条形的性质强烈取决于 N 和 $\epsilon$ 的选择。较大的引力软化参数和较少的粒子数会导致快速形成一个清晰且旋转速度较慢的条形。相反，较小的 $\epsilon$ 会产生小且快速旋转的盘，呈紧密卷曲的螺旋臂，随后形成较弱的条形。我们发现，随着 N 的增加，条形性质趋于收敛，软化参数的影响逐渐减弱。最后，在某些情况下，条形末端可观察到短暂的螺旋臂，这些螺旋臂定期从倒挂变为前绒，反之亦然——即“波动”。