[论文解读] NuGrid: Toward High Precision Double-Degenerate Merger Simulations with SPH in 3D

Steven Diehl, Chris L. Fryer+12 more|ArXiv.org|Nov 28, 2008

Gamma-ray bursts and supernovae被引用 3

一句话总结

本文通过改进初始条件设置（采用加权Voronoi镶嵌法）并引入理想气体状态方程以包含激波，对双简并白矮星并合进行了高精度三维SPH模拟，强调了严格的验证与验证。关键结果表明，激波显著改变了并合动力学，形成一个高温、膨胀的包层并喷射携带角动量的质量，这可能影响此类并合作为Ia型超新星前身体的可行性。

ABSTRACT

We present preliminary results from recent high-resolution double-degenerate merger simulations with the Smooth Particle Hydrodynamics (SPH) technique. We put particular emphasis on verification and validation in our effort and show the importance of details in the initial condition setup for the final outcome of the simulation. We also stress the dynamical importance of including shocks in the simulations. These results represent a first step toward a suite of simulations that will shed light on the question whether double-degenerate mergers are a viable path toward type 1a supernovae. In future simulations, we will make use of the capabilities of the NuGrid collaboration in post-processing SPH particle trajectories with a complete nuclear network to follow the detailed nuclear reactions during the dynamic merger phase.

研究动机与目标

• 开展高精度、经验证和验证的三维SPH模拟，研究双简并白矮星并合作为Ia型超新星前身体的可行性。
• 研究初始条件设置（特别是粒子分布和共转）对模拟精度和动力学的影响。
• 通过采用不同状态方程比较含激波与不含激波的模拟，评估激波在并合动力学中的作用。
• 利用完整的核反应网络对SPH粒子轨迹进行后处理，模拟动态核反应，并与R Coronae Borealis型星的元素丰度进行比较。
• 在NuGrid合作框架内，为未来使用真实状态方程和恒星演化模型的模拟奠定基础。

提出的方法

• 采用修改后的SNSPH代码，利用高分辨率光滑粒子流体动力学（SPH）在三维空间中模拟双白矮星并合。
• 采用自洽场方法生成共转、质量不等的双白矮星系统的平衡初始条件，最大限度减少数值伪影。
• 提出一种基于加权Voronoi镶嵌（WVT）的粒子布置方法，在外层实现高分辨率，同时最小化粒子噪声和密度波动。
• 通过比较使用多相态状态方程（无激波）与理想气体状态方程（含激波）的模拟，隔离激波加热的动力学影响。
• 利用tppnp工具进行后处理，追踪SPH粒子通过完整核反应网络，实现核素合成，并与观测到的R Coronae Borealis型星元素丰度进行比较。
• 通过与结构网格模拟的代码对比和数值收敛性研究进行验证，并与R Coronae Borealis型星的观测约束进行对比。

实验结果

研究问题

• RQ1 不同的初始条件设置（特别是粒子分布和旋转速度分布）如何影响双简并白矮星并合模拟的结果？
• RQ2 在状态方程中引入激波对并合演化和质量喷射的动力学影响是什么？
• RQ3 并合残余物的结构和角动量损失情况，与R Coronae Borealis型星的观测约束相比如何？
• RQ4 并合过程中形成的激波加热、高熵包层在多大程度上改变了并合 timescale 和能量预算？
• RQ5 高精度SPH模拟结合真实状态方程和核反应网络后处理，能否证实双简并白矮星并合作为Ia型超新星前身体的可行性？

主要发现

• 通过理想气体状态方程引入激波后，形成了一个高温、膨胀的激波物质包层，包裹住两颗白矮星，与无激波的多相态模拟相比，根本性地改变了并合动力学。
• 激波气体从两颗恒星的后方被高效喷射，携带显著的角动量，可能加速轨道衰减。
• 并合残余物最终形成一个具有快速旋转核心的非均匀旋转、近似球对称的结构，热包层延伸至约100个太阳半径。
• 模拟显示，吸积流平均由超过1,000个SPH粒子解析，表明由于基于WVT的粒子布置方法，分辨率高且数值噪声低。
• 非共转系统因角动量向自转转移而表现出人为的轨道收缩，凸显了正确初始旋转设置的重要性。
• 结果表明，状态方程和初始条件质量是决定并合结果的关键因素，激波在动力学中起核心作用。