随着人类对宇宙认知的不断深入，系外行星研究已成为天文学领域最引人入胜的前沿之一。自1992年首次发现围绕脉冲星PSR1257+12运行的行星以来，已知系外行星数量已增长至6000多颗，围绕4500多颗宿主恒星运行。这些行星在质量、大小和轨道方面展现出极端的多样性，从与地球一样小而岩石的天体，到比木星更庞大的气态巨行星，甚至接近褐矮星的质量极限。这种物理多样性不仅体现在行星本身，也反映在它们的宿主恒星上——恒星的质量、年龄、化学成分和演化阶段各不相同。
面对这一浩瀚的行星种群，欧洲空间局（ESA）的Ariel空间任务肩负着重要使命：通过观测数百颗系外行星的大气，特别是热和温气体巨行星，来表征其化学成分，从而建立首个连贯的行星分类系统并理解行星多样性的根源。然而，当前试图完全解码詹姆斯·韦伯空间望远镜（JWST）提供的复杂光谱信息的努力正受到我们模型和理解局限性的阻碍。
正是在这一背景下，OPAL（Origins of Planets for Ariel）关键科学项目应运而生。这项由意大利超级计算中心（ICSC）支持的研究，利用CINECA的LEONARDO预百亿亿级EuroHPC超级计算机，旨在创建具有JWST级别复杂度的详细合成系外行星光谱和大气模型库。这些模型将在受控环境下测试社区的工具和方法，从而最大化Ariel任务的科学产出。该研究发表于《Astronomy and Computing》期刊，展示了如何通过高性能计算方法解决行星形成这一高维度问题。
为开展这项研究，作者团队开发了一套完整的技术流程。他们利用GGChem基于宿主恒星成分量化原行星盘尘埃和岩石星子的平衡化学组成。通过Jade（Joint Astrochemistry and Disc Evolution）代码模拟原行星盘随时间的演化，追踪物理和化学过程。采用GroMiT（行星生长和迁移轨迹）代码通过蒙特卡洛方法约束可能的行星形成历史。使用Mercury-Arχes这一并行N体代码模拟生长和迁移行星与嵌入原行星盘中的星子盘之间的相互作用。借助Hephaestus组合不同盘成分模型与行星形成和迁移模拟，构建行星整体成分。利用FastChem计算化学平衡条件下的行星大气成分。通过Vulcan研究化学非平衡条件下的行星大气过程。所有模拟均在LEONARDO预百亿亿级基础设施上运行，充分利用其数据中心通用计算模块的高性能计算能力。
2. The OPAL project
OPAL项目作为Ariel联盟的关键科学项目，其主要任务是创建真实的合成大气光谱库，以测试和验证Ariel任务的数据处理工具。项目利用Arχes行星形成代码套件进行端到端的模拟活动，追踪形成中的巨行星从星周盘成分到其大气成分的天体化学路径。通过与Exoclimes套件中的FastChem和Vulcan等大气建模工具结合，OPAL能够产生平衡和非平衡条件下的现实大气成分。
2.1. The OPAL sample
OPAL模拟活动目前聚焦于三个行星系统：WASP-69、HD 209458和HIP 67522。这些系统被选入研究是因为它们被列入Ariel任务的目标列表，且Ariel恒星表征工作组已为其提供了均匀的物理和成分参数。这三个系统的行星平衡温度和辐照条件代表了Ariel将观测到的范围，使其成为理想的研究案例。
2.2. Computational aspects of the simulation campaign
OPAL模拟活动在LEONARDO预百亿亿级基础设施上运行，利用其数据中心通用计算模块。每个节点配备两个56核Intel Xeon Platinum 8480+ CPU和512 GB DDR5-4800 RAM。项目中的两个核心代码——Jade和Mercury-Arχes——具有不同的并行特性，需要针对性地优化资源分配策略，以最大化模拟活动的整体效率。
3. End-to-end modelling of planet formation: the OPAL pipeline
OPAL管道将Arχes套件中的代码与GGChem和Exoclimes套件结合，模拟宿主恒星、可能的原行星盘及其行星之间遗传联系的每个环节。该管道包括多个专门代码：Jade追踪原行星盘的物理和天体化学演化；GroMiT基于卵石吸积情景追踪行星形成如何受原生原行星盘演化特征影响；Mercury-Arχes研究嵌入原生盘中的行星系统的形成和演化；Hephaestus允许将来自Jade的盘成分天体化学信息与来自Mercury-Arχes的吸积历史结合，量化最终的行星成分。
3.1. GGChem: initial composition of the planetary building blocks
研究使用GGChem基于中平面的压力和温度剖面以及宿主恒星的成分来量化原生原行星盘中尘埃和岩石星子的平衡化学组成。通过结合来自分子云的挥发性分子丰度数据，该方法能够真实近似每颗恒星原行星盘储层的初始化学性质。
3.2. Jade: protoplanetary disc composition over space and time
Jade模拟由于动力学化学和物质传输相互作用导致的原行星盘随时间演化。该代码采用算子分裂方案来演化盘，在每个时间步顺序解决盘的演化物理和化学过程。模型聚焦于包含气体、尘埃和星子的II类原行星盘的中平面，探索不同的化学场景、尘埃颗粒种群和盘参数。
3.3. GRoMiT: constraining the possible planet formation histories
GroMiT代码通过卵石吸积模拟固体行星的生长和迁移，从月球质量的初始行星胚胎开始。一旦生长中的行星达到卵石隔离质量，它们开始扰动轨道外部的气体流动，产生局部压力最大值，有效地隔离行星与固体物质流。种群合成版本的GroMiT允许蒙特卡洛提取盘的多个初始特性和行星种子的初始条件，以识别研究行星的潜在形成路径。
3.4. Mercury-Arχes: modelling gas and planetesimal accretion by forming planets
Mercury-Arχes这一并行N体代码用于模拟生长和迁移行星与为三个研究恒星系统重建的原行星盘中嵌入的星子盘之间的相互作用。该代码允许模拟：形成巨行星的质量增长和行星半径演化；行星形成不同阶段的轨道迁移（I型和II型迁移）；盘重力对行星体运动的激发效应；盘气体对星子运动的空气动力学阻尼效应。
3.5. Hephaestus: building the bulk planetary compositions
Hephaestus组成后处理代码用于将不同的盘成分模型与不同的行星形成和迁移模拟结合。该代码存储来自Jade的原行星盘空间化学结构及其随时间演化的信息，以及来自Mercury-Arχes的生长行星对气体和星子的吸积信息，并利用它们确定行星随时间从不同源吸积的物质的元素组成。
3.6. FastChem: planetary atmospheres in chemical equilibrium
为推导选定巨行星在化学平衡条件下可能的大气，研究使用FastChem。该代码实现了一种半分析方法，用于最小化约396种中性物种和114种带电物种的吉布斯能量。通过将Hephaestus计算的元素丰度与行星大气的垂直温度-压力剖面结合，FastChem返回给定系外行星大气的垂直成分结构。
3.7. Vulcan: the role of chemical disequilibrium in planetary atmospheres
为研究大气动力学对最终垂直成分剖面的影响，研究使用Vulcan代码处理Hephaestus丰度的子集。该化学动力学模型直接积分一组质量连续性方程，直到通过Rosenbrock方法达到稳态。与FastChem等化学平衡模型相比，Vulcan允许包含各种非平衡诱导现象，如光化学和离子化学、通过平流、涡流扩散和分子扩散的垂直传输、过饱和物种的凝结和重力沉降以及内部脱气或大气逃逸。
4. Emerging complexity: early results of the OPAL campaign
OPAL模拟活动的早期结果已经展示了围绕单个恒星运行的个体行星可能具有的巨大多样性的整体和大气成分。这种多样性直接与行星形成过程初始条件的不确定性相关，并清楚地强调了一个配方并不适用于所有情况。以WASP-69b为例，研究显示广泛使用的C/O比在大多数情况下不能单独明确追踪特定的形成历史。C/N比很好地追踪了吸积到生长行星上的物质性质，但本身不允许轻易约束迁移。只有联合使用涉及不同挥发性元素的多个归一化元素比，才能打破简并性。
4.1. Computational lessons learned
在OPAL项目运行期间，研究团队确定了几个关键领域，优化管道及其代码将提高未来模拟活动的效率并增强其采样完整性。这些关键领域包括：管道自动化增强、代码设计优化和模拟活动的维度管理。当前管道版本中，每个单独代码的输出已经格式化为与后续管道步骤的互操作性，但模拟活动的管理仍由手工执行。此外，两个核心代码Mercury-Arχes和Jade的不同优化设计也限制了性能的可扩展性。高维端到端建模问题导致的简并形成路径是另一个瓶颈，需要识别决定最终结果的问题维度。
OPAL项目通过将多个专门模型结合到由LEONARDO超级计算集群高性能能力支持的集成大规模模拟活动中，应对系外行星研究中的高维挑战。这一开创性工作提供了前所未有的真实且高度详细的大气模型库，目前正被Ariel联盟用于构建具有JWST级别复杂度的合成光谱。这些合成光谱将在受控环境下测试和验证Ariel任务的代码、方法和管道，确保我们在面对Ariel观测将呈现的挑战时做好充分准备。该研究不仅为系外行星大气表征建立了新标准，也展示了高性能计算方法在解决天体物理学复杂问题中的变革性潜力，为未来系外行星研究奠定了坚实基础。
打赏
