[论文解读] soliton_solver: A GPU-based finite-difference PDE solver for topological solitons in two-dimensional non-linear field theories

Paul Leask|arXiv (Cornell University)|Mar 25, 2026

Scientific Research and Discoveries被引用 0

一句话总结

一个GPU加速、Python框架，具备面向理论无关的有限差分核心，用于跨多模型通过模块化理论组件进行二维拓扑孤子的仿真与可视化。

ABSTRACT

This paper introduces soliton_solver, an open-source GPU-accelerated software package for the simulation and real-time visualization of topological solitons in two-dimensional non-linear field theories. The software is structured around a theory-agnostic numerical core implemented using Numba CUDA kernels, while individual physical models are introduced through modular theory components. This separation enables a single computational framework to be applied across a broad class of systems, from nanoscale magnetic spin textures in condensed matter physics to cosmic strings spanning galaxies in high energy physics. The numerical backend provides finite-difference discretization, energy minimization, and GPU-resident evaluation of observables. A CUDA--PyOpenGL rendering pipeline allows direct visualization of evolving field configurations without staging full arrays through host memory. The package is distributed in Python via PyPI and supports both reproducible batch simulations and interactive exploration of metastable configurations, soliton interactions, and model-dependent initial states. We describe the software architecture, numerical workflow, and extensibility model, and we present representative example applications. We also outline how additional theories can be incorporated with minimal modification of the shared numerical infrastructure.

研究动机与目标

• 研究多样物理系统中的拓扑孤子的动机与需要一个可重用的计算工具。
• 提供一个理论无关的GPU后端，通过模块化理论组件驱动多种非线性场理论。
• 实现实时可视化和对中间稳态配置与孤子相互作用的交互式探索。
• 提供一个开源包（PyPI），支持批量仿真与新模型的快速原型设计。

提出的方法

• 用于二维场的GPU加速有限差分引擎，具有共享参数系统和halo增强的网格。
• 用于能量最小化以高效放松到静态孤子配置的 arrested Newton flow。
• 通过Numba CUDA实现的CUDA内核，配合理论注册表在运行时加载模型。
• 用于交互式、无主机内存的可视化的CUDA–PyOpenGL渲染管线，展示演化场。
• 模块化理论接口，每个模型提供场定义、能量泛函、参数、初始化以及可选的可视化工具。

Figure 1 : Typical usage of soliton_solver for a micromagnetic simulation of anti-skyrmions in Heusler compounds, modelled by a chiral ferromagnet with the Dzyaloshinskii–Moriya interaction (DMI), and including the effects of demagnetization.

实验结果

研究问题

• RQ1 一个单一计算框架如何支持具有拓扑孤子的多种二维非线性场理论？
• RQ2 与标准梯度下降相比，被抑制的牛顿流在孤子弛豫速度和鲁棒性方面的影响如何？
• RQ3 是否可以将实时GPU驻留可视化与求解器集成，以在不同理论之间实现交互式探索？
• RQ4 在对共享数值基础设施最少修改的情况下，新增理论有多容易被加入框架？

主要发现

• 该包提供一个可重用的计算框架，具备理论无关的核心和模块化理论接口。
• 内置理论覆盖Abelian Higgs宇宙字符串、Ginzburg–Landau超导体中的涡旋、anyon、在各种介质中的Skyrmion以及baby Skyrmion等。
• 该框架提供一个GPU驻留的CUDA–PyOpenGL可视化后端，适用于对演化配置的实时检查。
• arrested Newton flow通过利用虚构时间动力学并在能量增加时停止来加速放松，提高多孤子问题的效率。
• 通过实现紧凑的理论模块来添加新模型；求解器核心处理网格构建、微分、时间积分和渲染。
• 工作流支持交互式探索和批量仿真，安装和使用可通过PyPI或源码完成。

Figure 2 : A single anti-skyrmion in the bulk of a Heusler compound, with the associated Dzyaloshinskii–Moriya interaction and backreaction of the demagnetization field. All panels are in dimensionless units. The left panel shows the energy density $\mathcal{E}(\vec{x})$ , the middle shows the magne