首页 > 最新文献

Computer Physics Communications最新文献

英文 中文
Combustion Toolbox: An open-source thermochemical code for gas- and condensed-phase problems involving chemical equilibrium 燃烧工具箱:一个开放源代码的热化学代码,用于解决涉及化学平衡的气相和冷凝相问题
IF 3.4 2区 物理与天体物理 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2026-03-01 Epub Date: 2025-12-18 DOI: 10.1016/j.cpc.2025.110004
Alberto Cuadra, César Huete, Marcos Vera
The Combustion Toolbox (CT) is a newly developed open-source thermochemical code designed to solve problems involving chemical equilibrium for both gas- and condensed-phase species. The kernel of the code is based on the theoretical framework set forth by NASA’s computer program CEA (Chemical Equilibrium with Applications) while incorporating new algorithms that significantly improve both convergence rate and robustness. The thermochemical properties are computed under the ideal gas approximation using an up-to-date version of NASA’s 9-coefficient polynomial fits. These fits use the Third Millennium database, which includes the available values from Active Thermochemical Tables. Combustion Toolbox is programmed in MATLAB with an object-oriented architecture composed of three main modules: CT-EQUIL, CT-SD, and CT-ROCKET. The kernel module, CT-EQUIL, minimizes the Gibbs/Helmholtz free energy of the system using the technique of Lagrange multipliers combined with a multidimensional Newton-Raphson method, upon the condition that two state functions are used to define the mixture properties (e.g., enthalpy and pressure). CT-SD solves processes involving strong changes in dynamic pressure, such as steady shock and detonation waves under normal and oblique incidence angles. Finally, CT-ROCKET estimates rocket engine performance under highly idealized conditions. The new tool is equipped with a versatile Graphical User Interface and has been successfully used for teaching and research activities over the last six years. Results are in excellent agreement with CEA, Cantera within Caltech’s Shock and Detonation Toolbox (SD-Toolbox), and the Thermochemical Equilibrium Abundances (TEA) code. CT is available under an open-source GPLv3 license via GitHub https://github.com/CombustionToolbox/combustion_toolbox, and its documentation can be found in https://combustion-toolbox-website.readthedocs.io.
燃烧工具箱(CT)是一个新开发的开源热化学代码,旨在解决涉及气相和冷凝相物质的化学平衡问题。代码的核心是基于NASA计算机程序CEA(化学平衡与应用)提出的理论框架,同时结合了显着提高收敛速度和鲁棒性的新算法。热化学性质在理想气体近似下计算,使用最新版本的NASA 9系数多项式拟合。这些拟合使用第三个千年数据库,其中包括活性热化学表中的可用值。燃烧工具箱是在MATLAB编程与面向对象的体系结构组成的三个主要模块:CT-EQUIL, CT-SD,和CT-ROCKET。核心模块CT-EQUIL在使用两个状态函数定义混合性质(如焓和压力)的条件下,使用拉格朗日乘子技术结合多维牛顿-拉夫森方法最小化系统的吉布斯/亥姆霍兹自由能。CT-SD解决了涉及动压剧烈变化的过程,例如在正入射和斜入射下的稳定激波和爆震波。最后,CT-ROCKET在高度理想化的条件下估计火箭发动机的性能。这个新工具配备了一个多功能的图形用户界面,在过去的六年里已经成功地用于教学和研究活动。结果与CEA,加州理工学院的冲击和引爆工具箱(SD-Toolbox)中的Cantera以及热化学平衡丰度(TEA)代码非常一致。CT在开源GPLv3许可下可通过GitHub https://github.com/CombustionToolbox/combustion_toolbox获得,其文档可在https://combustion-toolbox-website.readthedocs.io找到。
{"title":"Combustion Toolbox: An open-source thermochemical code for gas- and condensed-phase problems involving chemical equilibrium","authors":"Alberto Cuadra,&nbsp;César Huete,&nbsp;Marcos Vera","doi":"10.1016/j.cpc.2025.110004","DOIUrl":"10.1016/j.cpc.2025.110004","url":null,"abstract":"<div><div>The Combustion Toolbox (CT) is a newly developed open-source thermochemical code designed to solve problems involving chemical equilibrium for both gas- and condensed-phase species. The kernel of the code is based on the theoretical framework set forth by NASA’s computer program CEA (Chemical Equilibrium with Applications) while incorporating new algorithms that significantly improve both convergence rate and robustness. The thermochemical properties are computed under the ideal gas approximation using an up-to-date version of NASA’s 9-coefficient polynomial fits. These fits use the Third Millennium database, which includes the available values from Active Thermochemical Tables. Combustion Toolbox is programmed in MATLAB with an object-oriented architecture composed of three main modules: CT-EQUIL, CT-SD, and CT-ROCKET. The kernel module, CT-EQUIL, minimizes the Gibbs/Helmholtz free energy of the system using the technique of Lagrange multipliers combined with a multidimensional Newton-Raphson method, upon the condition that two state functions are used to define the mixture properties (e.g., enthalpy and pressure). CT-SD solves processes involving strong changes in dynamic pressure, such as steady shock and detonation waves under normal and oblique incidence angles. Finally, CT-ROCKET estimates rocket engine performance under highly idealized conditions. The new tool is equipped with a versatile Graphical User Interface and has been successfully used for teaching and research activities over the last six years. Results are in excellent agreement with CEA, Cantera within Caltech’s Shock and Detonation Toolbox (SD-Toolbox), and the Thermochemical Equilibrium Abundances (TEA) code. CT is available under an open-source GPLv3 license via GitHub <span><span>https://github.com/CombustionToolbox/combustion_toolbox</span><svg><path></path></svg></span>, and its documentation can be found in <span><span>https://combustion-toolbox-website.readthedocs.io</span><svg><path></path></svg></span>.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"320 ","pages":"Article 110004"},"PeriodicalIF":3.4,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145880067","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Automation of a matching on-shell calculator 一个匹配的壳上计算器的自动化
IF 3.4 2区 物理与天体物理 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2026-03-01 Epub Date: 2025-11-19 DOI: 10.1016/j.cpc.2025.109935
Javier López Miras, Fuensanta Vilches
We introduce mosca, a Mathematica package designed to facilitate on-shell calculations in effective field theories (EFTs). This initial release focuses on the reduction of Green’s bases to physical bases, as well as transformations between arbitrary operator bases. The core of the package is based on a diagrammatic on-shell matching procedure, grounded in the equivalence of physical observables derived from both redundant and non-redundant Lagrangians. mosca offers a complete set of tools for performing basis transformations, diagram isomorphism detection, numerical substitution of kinematic configurations, and symbolic manipulation of algebraic expressions. Planned future developments include extension to one-loop computations, thus providing support for EFT renormalization directly in a physical basis and automated computation of one-loop finite matching, including contributions from evanescent operators.
PROGRAM SUMMARY Program Title: mosca CPC Library link to program files: (to be added by Technical Editor) Developer’s repository link: https://gitlab.com/matchingonshell/mosca Licensing provisions: GPLv3 Programming language: Mathematica Nature of problem: Matching calculations in effective field theories are traditionally performed off-shell, involving complicated basis reductions through non-trivial field redefinitions to eliminate redundant operators. This process is algebraically intensive and prone to errors. Although on-shell matching, which focuses directly on physical observables, could simplify these steps by avoiding field redefinitions, it has been considered impractical due to the presence of apparent non-localities that must cancel precisely. Automating on-shell matching has therefore been a long-standing challenge. Solution method: Our approach is based on a numerical solution of the on-shell matching equations, which naturally and effortlessly enforces the delicate cancellation of non-local terms between the full theory and the effective theory. By employing rational on-shell kinematics, the method achieves an exact analytic solution despite using numerical techniques. This allows the matching to be performed entirely within a physical operator basis. Additional comments including restrictions and unusual features: The workflow for handling Lagrangians and Feynman diagrams in mosca is based on the integration of FeynArts and FeynCalc. Consequently, users need to provide specific FeynArts model files patched for compatibility with FeynCalc. Additionally, a specialized input format is required to define Wilson coefficients along with their corresponding EFT order (EFTOrder). These requirements ensure the correct processing of models and coefficients.
我们介绍mosca,一个旨在促进有效场论(eft)的壳上计算的Mathematica软件包。这个初始版本的重点是将Green的基减少到物理基,以及任意算子基之间的转换。该包的核心是基于一个图解的壳匹配过程,基于从冗余和非冗余拉格朗日导出的物理可观测值的等效性。Mosca提供了一套完整的工具,用于执行基变换、图同构检测、运动学构型的数值替换和代数表达式的符号操作。计划的未来发展包括扩展到单回路计算,从而在物理基础上直接支持EFT重规范化和单回路有限匹配的自动计算,包括来自倏逝算子的贡献。程序摘要程序标题:mosca CPC库链接到程序文件:(由技术编辑添加)开发人员的存储库链接:https://gitlab.com/matchingonshell/mosca许可条款:GPLv3编程语言:Mathematica问题的性质:有效场理论中的匹配计算传统上是off-shell执行的,涉及通过非琐碎的字段重新定义来消除冗余运算符的复杂基约简。这个过程需要大量的代数运算,而且容易出错。尽管直接关注物理可观察对象的壳上匹配可以通过避免字段重新定义来简化这些步骤,但由于存在必须精确取消的明显非位置,它被认为是不切实际的。因此,在shell上的自动匹配一直是一个长期的挑战。求解方法:我们的方法是基于壳上匹配方程的数值解,它自然而毫不费力地强制在满理论和有效理论之间微妙地消去非局部项。通过采用有理壳上运动学,该方法实现了精确的解析解,尽管使用数值技术。这允许完全在物理操作符的基础上执行匹配。在mosca中处理拉格朗日图和费曼图的工作流程是基于feynnts和FeynCalc的集成。因此,用户需要提供特定的feynnts模型文件补丁,以便与FeynCalc兼容。此外,需要一种专门的输入格式来定义Wilson系数及其相应的EFT顺序(EFTOrder)。这些要求保证了模型和系数的正确处理。
{"title":"Automation of a matching on-shell calculator","authors":"Javier López Miras,&nbsp;Fuensanta Vilches","doi":"10.1016/j.cpc.2025.109935","DOIUrl":"10.1016/j.cpc.2025.109935","url":null,"abstract":"<div><div>We introduce <span>mosca</span>, a <span>Mathematica</span> package designed to facilitate on-shell calculations in effective field theories (EFTs). This initial release focuses on the reduction of Green’s bases to physical bases, as well as transformations between arbitrary operator bases. The core of the package is based on a diagrammatic on-shell matching procedure, grounded in the equivalence of physical observables derived from both redundant and non-redundant Lagrangians. <span>mosca</span> offers a complete set of tools for performing basis transformations, diagram isomorphism detection, numerical substitution of kinematic configurations, and symbolic manipulation of algebraic expressions. Planned future developments include extension to one-loop computations, thus providing support for EFT renormalization directly in a physical basis and automated computation of one-loop finite matching, including contributions from evanescent operators.</div><div><strong>PROGRAM SUMMARY</strong> <em>Program Title:</em> mosca <em>CPC Library link to program files:</em> (to be added by Technical Editor) <em>Developer’s repository link:</em> <span><span>https://gitlab.com/matchingonshell/mosca</span><svg><path></path></svg></span> <em>Licensing provisions:</em> GPLv3 <em>Programming language:</em> Mathematica <em>Nature of problem:</em> Matching calculations in effective field theories are traditionally performed off-shell, involving complicated basis reductions through non-trivial field redefinitions to eliminate redundant operators. This process is algebraically intensive and prone to errors. Although on-shell matching, which focuses directly on physical observables, could simplify these steps by avoiding field redefinitions, it has been considered impractical due to the presence of apparent non-localities that must cancel precisely. Automating on-shell matching has therefore been a long-standing challenge. <em>Solution method:</em> Our approach is based on a numerical solution of the on-shell matching equations, which naturally and effortlessly enforces the delicate cancellation of non-local terms between the full theory and the effective theory. By employing rational on-shell kinematics, the method achieves an exact analytic solution despite using numerical techniques. This allows the matching to be performed entirely within a physical operator basis. <em>Additional comments including restrictions and unusual features:</em> The workflow for handling Lagrangians and Feynman diagrams in mosca is based on the integration of FeynArts and FeynCalc. Consequently, users need to provide specific FeynArts model files patched for compatibility with FeynCalc. Additionally, a specialized input format is required to define Wilson coefficients along with their corresponding EFT order (EFTOrder). These requirements ensure the correct processing of models and coefficients.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"320 ","pages":"Article 109935"},"PeriodicalIF":3.4,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145681672","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
SpectraMatcher: A python program for interactive analysis and peak assignment of vibronic spectra SpectraMatcher:一个用于交互式分析和频谱峰分配的python程序
IF 3.4 2区 物理与天体物理 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2026-03-01 Epub Date: 2025-11-11 DOI: 10.1016/j.cpc.2025.109945
Johanna Langner, Isabelle Weber, Henryk A. Witek, Yuan-Pern Lee
SpectraMatcher is a cross-platform desktop application for interactive comparison of experimental and computed vibronic spectra, designed to assist in the recognition and assignment of spectral patterns. It provides an intuitive graphical interface — with no coding or scripting required — for importing experimental spectra, visualizing them alongside the corresponding theoretical spectra constructed from Gaussian frequency calculations, and adjusting key parameters such as peak width, intensity scaling factors, and vibration-type-specific anharmonic corrections. SpectraMatcher features an automated peak-matching algorithm that assigns experimental and computed peaks based on their intensity ratio and proximity. Assignments and spectra can be exported in multiple formats for publication or for further analysis. The software remains responsive even for large datasets, and supports efficient and reproducible interpretation of vibronic spectra.
SpectraMatcher是一个跨平台的桌面应用程序,用于实验和计算振动谱的交互式比较,旨在帮助识别和分配光谱模式。它提供了一个直观的图形界面-无需编码或脚本-用于导入实验光谱,将它们与从高斯频率计算构建的相应理论光谱一起可视化,并调整关键参数,如峰宽,强度缩放因子和振动类型特定的非谐波校正。SpectraMatcher具有自动峰值匹配算法,根据它们的强度比和接近度分配实验和计算的峰值。分配和光谱可以以多种格式导出,以供发布或进一步分析。该软件即使对大型数据集也保持响应,并支持有效和可重复的振动谱解释。
{"title":"SpectraMatcher: A python program for interactive analysis and peak assignment of vibronic spectra","authors":"Johanna Langner,&nbsp;Isabelle Weber,&nbsp;Henryk A. Witek,&nbsp;Yuan-Pern Lee","doi":"10.1016/j.cpc.2025.109945","DOIUrl":"10.1016/j.cpc.2025.109945","url":null,"abstract":"<div><div><em>SpectraMatcher</em> is a cross-platform desktop application for interactive comparison of experimental and computed vibronic spectra, designed to assist in the recognition and assignment of spectral patterns. It provides an intuitive graphical interface — with no coding or scripting required — for importing experimental spectra, visualizing them alongside the corresponding theoretical spectra constructed from Gaussian frequency calculations, and adjusting key parameters such as peak width, intensity scaling factors, and vibration-type-specific anharmonic corrections. <em>SpectraMatcher</em> features an automated peak-matching algorithm that assigns experimental and computed peaks based on their intensity ratio and proximity. Assignments and spectra can be exported in multiple formats for publication or for further analysis. The software remains responsive even for large datasets, and supports efficient and reproducible interpretation of vibronic spectra.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"320 ","pages":"Article 109945"},"PeriodicalIF":3.4,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145616217","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Fantômas unconfined: global QCD fits with Bézier parameterizations Fantômas unconstrained:全局QCD与bsamizier参数化拟合
IF 3.4 2区 物理与天体物理 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2026-03-01 Epub Date: 2025-12-14 DOI: 10.1016/j.cpc.2025.109969
Lucas Kotz , Aurore Courtoy , T.J. Hobbs , Pavel Nadolsky , Fredrick Olness , Maximiliano Ponce-Chavez , Varada Purohit
Fantômas is a C++ toolkit for exploring the parametrization dependence of parton distribution functions (PDFs) and other correlator functions in quantum chromodynamics (QCD). Fantômas facilitates the generation of adaptable polynomial parametrizations for PDFs, called metamorphs, to find best-fit PDF solutions and quantify the epistemic uncertainty associated with the parametrizations during their fitting. The method employs Bézier curves as universal approximators for a variety of PDF shapes. Integrated into the xFitter framework for the global QCD analysis, Fantômas provides a foundation for general models of PDFs, while reducing the computational time compared to the approaches utilizing traditional polynomial parametrizations as well as providing an interpretable alternative to neural-network-based models. This paper outlines the structure and practical usage of the Fantômas toolkit, including its inputs, outputs, and implementation within xFitter. It also provides a practical example of using Fantômas for uncertainty quantification as well as the combination of PDF fits into a single ensemble.
Fantômas是一个c++工具包,用于探索量子色动力学(QCD)中部子分布函数(pdf)和其他相关函数的参数化依赖。Fantômas有助于生成PDF的自适应多项式参数化,称为变形,以找到最适合的PDF解,并在拟合期间量化与参数化相关的认知不确定性。该方法采用bsamizier曲线作为各种PDF形状的通用逼近器。集成到xFitter框架中用于全局QCD分析,Fantômas为pdf的一般模型提供了基础,同时与使用传统多项式参数化的方法相比减少了计算时间,并提供了基于神经网络的模型的可解释替代方案。本文概述了Fantômas工具包的结构和实际用法,包括它的输入、输出和xFitter中的实现。它还提供了一个使用Fantômas进行不确定度量化的实际示例,以及PDF适合于单个集成的组合。
{"title":"Fantômas unconfined: global QCD fits with Bézier parameterizations","authors":"Lucas Kotz ,&nbsp;Aurore Courtoy ,&nbsp;T.J. Hobbs ,&nbsp;Pavel Nadolsky ,&nbsp;Fredrick Olness ,&nbsp;Maximiliano Ponce-Chavez ,&nbsp;Varada Purohit","doi":"10.1016/j.cpc.2025.109969","DOIUrl":"10.1016/j.cpc.2025.109969","url":null,"abstract":"<div><div><span>Fantômas</span> is a C++ toolkit for exploring the parametrization dependence of parton distribution functions (PDFs) and other correlator functions in quantum chromodynamics (QCD). <span>Fantômas</span> facilitates the generation of adaptable polynomial parametrizations for PDFs, called metamorphs, to find best-fit PDF solutions and quantify the epistemic uncertainty associated with the parametrizations during their fitting. The method employs Bézier curves as universal approximators for a variety of PDF shapes. Integrated into the <span>xFitter</span> framework for the global QCD analysis, <span>Fantômas</span> provides a foundation for general models of PDFs, while reducing the computational time compared to the approaches utilizing traditional polynomial parametrizations as well as providing an interpretable alternative to neural-network-based models. This paper outlines the structure and practical usage of the <span>Fantômas</span> toolkit, including its inputs, outputs, and implementation within <span>xFitter</span>. It also provides a practical example of using <span>Fantômas</span> for uncertainty quantification as well as the combination of PDF fits into a single ensemble.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"320 ","pages":"Article 109969"},"PeriodicalIF":3.4,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145880064","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
PWACG: Partial wave analysis code generator supporting Newton-conjugate gradient method 部分波分析代码生成器支持牛顿共轭梯度法
IF 3.4 2区 物理与天体物理 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2026-03-01 Epub Date: 2025-12-02 DOI: 10.1016/j.cpc.2025.109963
Xiang Dong , Yu-Chang Sun , Chu-Cheng Pan , Ao-Yan Cheng , Ao-Bo Wang , Hao Cai , Kai Zhu
This paper introduces a novel Partial Wave Analysis Code Generator (PWACG) that automatically generates high-performance partial wave analysis codes. This is achieved by leveraging the JAX automatic differentiation library and the jinja2 template engine. The resulting code is constructed using the high-performance API of JAX, and includes support for the Newton’s Conjugate Gradient optimization method, as well as the full utilization of parallel computing capabilities offered by GPUs. By harnessing these advanced computing techniques, PWACG demonstrates a significant advantage in efficiently identifying global optimal points compared to conventional partial wave analysis software packages.
PROGRAM SUMMARY
Program Title: PWACG: Partial Wave Analysis Code Generator
CPC Library link to program files: https://doi.org/10.17632/47ppcnx77x.1
Licensing provisions: This software is distributed under the MIT License.
Programming language: Python
External routines/libraries: jaxlib, jax, jinja2, matplotlib, numpy, scipy
Nature of problem: The program addresses the need for high-performance computational tools in high-energy physics partial wave analysis (PWA). It introduces the Newton-conjugate gradient method for optimization, enhancing the accuracy and stability of fits.
Solution method: PWACG employs code generation and automatic differentiation techniques to automate the creation of PWA code. It leverages the computational capabilities of JAX for efficient execution and supports advanced features such as multi-GPU computation.
GitHub repository: https://github.com/caihao/PWACG
介绍了一种能自动生成高性能部分波分析代码的部分波分析代码发生器(pacg)。这是通过利用JAX自动区分库和jinja2模板引擎实现的。生成的代码使用JAX的高性能API构建,并包括对牛顿共轭梯度优化方法的支持,以及对gpu提供的并行计算能力的充分利用。通过利用这些先进的计算技术,与传统的部分波分析软件包相比,pwag在有效识别全局最优点方面具有显着优势。•程序标题:PWACG:部分波分析代码生成器•CPC库链接到程序文件:https://doi.org/10.17632/47ppcnx77x.1•许可条款:本软件在麻省理工学院许可下分发。•外部例程/库:jaxlib, jax, jinja2, matplotlib, numpy, scipy•问题性质:该程序解决了高能物理分波分析(PWA)中高性能计算工具的需求。引入牛顿共轭梯度法进行优化,提高了拟合的精度和稳定性。•解决方法:pwag采用代码生成和自动区分技术,自动创建PWA代码。它利用JAX的计算能力进行高效执行,并支持多gpu计算等高级特性。•GitHub存储库:https://github.com/caihao/PWACG
{"title":"PWACG: Partial wave analysis code generator supporting Newton-conjugate gradient method","authors":"Xiang Dong ,&nbsp;Yu-Chang Sun ,&nbsp;Chu-Cheng Pan ,&nbsp;Ao-Yan Cheng ,&nbsp;Ao-Bo Wang ,&nbsp;Hao Cai ,&nbsp;Kai Zhu","doi":"10.1016/j.cpc.2025.109963","DOIUrl":"10.1016/j.cpc.2025.109963","url":null,"abstract":"<div><div>This paper introduces a novel Partial Wave Analysis Code Generator (PWACG) that automatically generates high-performance partial wave analysis codes. This is achieved by leveraging the JAX automatic differentiation library and the jinja2 template engine. The resulting code is constructed using the high-performance API of JAX, and includes support for the Newton’s Conjugate Gradient optimization method, as well as the full utilization of parallel computing capabilities offered by GPUs. By harnessing these advanced computing techniques, PWACG demonstrates a significant advantage in efficiently identifying global optimal points compared to conventional partial wave analysis software packages.</div><div><strong>PROGRAM SUMMARY</strong></div><div>• <strong>Program Title:</strong> PWACG: Partial Wave Analysis Code Generator</div><div>• <strong>CPC Library link to program files:</strong> <span><span>https://doi.org/10.17632/47ppcnx77x.1</span><svg><path></path></svg></span></div><div>• <strong>Licensing provisions:</strong> This software is distributed under the MIT License.</div><div>• <strong>Programming language:</strong> Python</div><div>• <strong>External routines/libraries:</strong> jaxlib, jax, jinja2, matplotlib, numpy, scipy</div><div>• <strong>Nature of problem:</strong> The program addresses the need for high-performance computational tools in high-energy physics partial wave analysis (PWA). It introduces the Newton-conjugate gradient method for optimization, enhancing the accuracy and stability of fits.</div><div>• <strong>Solution method:</strong> PWACG employs code generation and automatic differentiation techniques to automate the creation of PWA code. It leverages the computational capabilities of JAX for efficient execution and supports advanced features such as multi-GPU computation.</div><div>• <strong>GitHub repository:</strong> <span><span>https://github.com/caihao/PWACG</span><svg><path></path></svg></span></div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"320 ","pages":"Article 109963"},"PeriodicalIF":3.4,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145786464","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
JANC: A cost-effective, differentiable compressible reacting flow solver featured with JAX-based adaptive mesh refinement JANC:一个具有成本效益,可微的可压缩反应流求解器,具有基于jax的自适应网格细化功能
IF 3.4 2区 物理与天体物理 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2026-02-01 Epub Date: 2025-10-25 DOI: 10.1016/j.cpc.2025.109915
Haocheng Wen , Faxuan Luo , Sheng Xu, Bing Wang
<div><div>The compressible reacting flow numerical solver is an essential tool in the study of combustion, energy disciplines, as well as in the design of industrial power and propulsion devices. We have established the first JAX-based (a Python library developed by Google for accelerator-oriented array computation and high-performance numerical computing) block-structured adaptive mesh refinement (AMR) framework, called JAX-AMR, and then developed a fully-differentiable solver for compressible reacting flows, named JANC. JANC is implemented in Python and features automatic differentiation capabilities, enabling an efficient integration of the solver with machine learning. Furthermore, benefited by multiple acceleration features such as accelerated linear algebra (XLA)-powered Just-In-Time (JIT) compilation, GPU/TPU computing, parallel computing, and AMR, the computational efficiency of JANC has been significantly improved. In a comparative test of a two-dimensional detonation tube case, the computational cost of the JANC core solver, running on a single A100 GPU, was reduced to 1% of that of OpenFOAM, which was parallelized across 384 CPU cores. When the AMR method is enabled for both solvers, JANC’s computational cost can be reduced to 1-2% of that of OpenFOAM. The core solver of JANC has also been tested for parallel computation on a 4-card A100 setup, demonstrating its convenient and efficient parallel computing capability. JANC also shows strong compatibility with machine learning by combining adjoint optimization to make the whole dynamic trajectory efficiently differentiable. JANC provides a new generation of high-performance, cost-effective, and high-precision solver framework for large-scale numerical simulations of compressible reacting flows and related machine learning research.</div><div>Program summary</div><div><em>Program title</em>: JAX-AMR and JANC</div><div><em>CPC Library link to program files</em>: <span><span>https://doi.org/10.17632/pkbxp5tm8w.1</span><svg><path></path></svg></span></div><div><em>Developer’s repository link</em>: <span><span>https://github.com/JA4S/JAX-AMR</span><svg><path></path></svg></span>, <span><span>https://github.com/JA4S/JANC</span><svg><path></path></svg></span></div><div>Licensing provisions: MIT</div><div>Programming language: Python</div><div><em>Nature of problem</em>: The numerical solution of compressible reactive flows plays a crucial role in combustion, energy utilization, and the design and manufacturing of propulsion systems. However, the multi-species nature, highly transient behavior, and strong numerical stiffness of reactive flows lead to significantly higher computational costs compared to conventional flow problems. In addition, conventional reactive flow solvers are typically built on Fortran or C++ frameworks, making them difficult to integrate with data-driven methods based on existing Python ecosystems—particularly gradient-based optimization techniques such as machine learning
可压缩反应流数值求解器是研究燃烧、能源学科以及设计工业动力和推进装置的重要工具。我们建立了第一个基于jax(谷歌开发的用于面向加速器的阵列计算和高性能数值计算的Python库)的块结构自适应网格细化(AMR)框架,称为JAX-AMR,然后开发了一个可压缩反应流的完全可微求解器,称为JANC。JANC是在Python中实现的,并具有自动区分功能,从而实现求解器与机器学习的有效集成。此外,得益于加速线性代数(XLA)驱动的JIT (Just-In-Time)编译、GPU/TPU计算、并行计算和AMR等多种加速特性,JANC的计算效率得到了显著提高。在二维爆爆管案例的对比测试中,在单个A100 GPU上运行的JANC核心求解器的计算成本降低到OpenFOAM的1%,OpenFOAM在384个CPU核心上并行化。当两个求解器都启用AMR方法时,JANC的计算成本可以降低到OpenFOAM的1-2%。该核心求解器在4卡A100装置上进行了并行计算测试,证明了其方便高效的并行计算能力。通过结合伴随优化使整个动态轨迹高效可微,JANC与机器学习具有较强的兼容性。JANC为可压缩反应流的大规模数值模拟和相关的机器学习研究提供了新一代高性能、经济高效、高精度的求解器框架。项目摘要项目标题:JAX-AMR和JANCCPC库链接到程序文件:https://doi.org/10.17632/pkbxp5tm8w.1Developer的存储库链接:https://github.com/JA4S/JAX-AMR, https://github.com/JA4S/JANCLicensing条款:mit编程语言:python问题性质:可压缩反应流的数值解在燃烧、能量利用和推进系统的设计和制造中起着至关重要的作用。然而,与传统流动问题相比,反应流动的多物种特性、高瞬态特性和强数值刚度导致计算成本显着增加。此外,传统的响应式流求解器通常是在Fortran或c++框架上构建的,这使得它们很难与基于现有Python生态系统的数据驱动方法集成,特别是基于梯度的优化技术,如机器学习。尽管已经开发了几个基于JAX的微分求解器原型来解决这些问题,利用gpu的多核计算能力和JAX的自动微分框架来进行高效和可微分的流模拟,但还没有开发出基于JAX的可压缩反应流求解器。更关键的是,由于反应流的多尺度特性,目前还没有与JAX兼容的自适应网格细化(AMR)框架,这使得大规模、高分辨率的模拟非常昂贵且难以执行。解决方法:我们已经开发了JANC:一个具有成本效益,可微的可压缩反应流求解器,具有基于jax的自适应网格细化的特点。JANC是建立在JAX框架上的求解器,充分利用了JAX的XLA高性能计算编译和自动区分功能,与传统的响应式流求解器(如OpenFOAM)相比,具有显著的加速。通过结合基于伴随的优化算法,我们的求解器有效地兼容现有的可微优化方法,提供更低的内存消耗和更快的可微计算。在求解器的基础上,提出了一种基于固定位置和形状的动态更新多级块的自适应网格细化框架JAX-AMR。该框架与JAX的即时(JIT)编译和矢量化(vmap)操作完全兼容,能够将AMR功能有效地集成到流求解器中。此外,我们的反应流求解器支持多设备并行,以满足大规模问题的需求。JANC为定义边界条件、源项和其他特定于问题的功能提供了高度自定义和用户友好的界面。JANC可以通过pip install作为Python包轻松安装,并且可以在单个Jupyter Notebook中执行响应式流程模拟的整个工作流。附加注释,包括限制和不寻常的功能:JANC的一些功能依赖于开源的第三方Python库,这些库被设置为自动安装。
{"title":"JANC: A cost-effective, differentiable compressible reacting flow solver featured with JAX-based adaptive mesh refinement","authors":"Haocheng Wen ,&nbsp;Faxuan Luo ,&nbsp;Sheng Xu,&nbsp;Bing Wang","doi":"10.1016/j.cpc.2025.109915","DOIUrl":"10.1016/j.cpc.2025.109915","url":null,"abstract":"&lt;div&gt;&lt;div&gt;The compressible reacting flow numerical solver is an essential tool in the study of combustion, energy disciplines, as well as in the design of industrial power and propulsion devices. We have established the first JAX-based (a Python library developed by Google for accelerator-oriented array computation and high-performance numerical computing) block-structured adaptive mesh refinement (AMR) framework, called JAX-AMR, and then developed a fully-differentiable solver for compressible reacting flows, named JANC. JANC is implemented in Python and features automatic differentiation capabilities, enabling an efficient integration of the solver with machine learning. Furthermore, benefited by multiple acceleration features such as accelerated linear algebra (XLA)-powered Just-In-Time (JIT) compilation, GPU/TPU computing, parallel computing, and AMR, the computational efficiency of JANC has been significantly improved. In a comparative test of a two-dimensional detonation tube case, the computational cost of the JANC core solver, running on a single A100 GPU, was reduced to 1% of that of OpenFOAM, which was parallelized across 384 CPU cores. When the AMR method is enabled for both solvers, JANC’s computational cost can be reduced to 1-2% of that of OpenFOAM. The core solver of JANC has also been tested for parallel computation on a 4-card A100 setup, demonstrating its convenient and efficient parallel computing capability. JANC also shows strong compatibility with machine learning by combining adjoint optimization to make the whole dynamic trajectory efficiently differentiable. JANC provides a new generation of high-performance, cost-effective, and high-precision solver framework for large-scale numerical simulations of compressible reacting flows and related machine learning research.&lt;/div&gt;&lt;div&gt;Program summary&lt;/div&gt;&lt;div&gt;&lt;em&gt;Program title&lt;/em&gt;: JAX-AMR and JANC&lt;/div&gt;&lt;div&gt;&lt;em&gt;CPC Library link to program files&lt;/em&gt;: &lt;span&gt;&lt;span&gt;https://doi.org/10.17632/pkbxp5tm8w.1&lt;/span&gt;&lt;svg&gt;&lt;path&gt;&lt;/path&gt;&lt;/svg&gt;&lt;/span&gt;&lt;/div&gt;&lt;div&gt;&lt;em&gt;Developer’s repository link&lt;/em&gt;: &lt;span&gt;&lt;span&gt;https://github.com/JA4S/JAX-AMR&lt;/span&gt;&lt;svg&gt;&lt;path&gt;&lt;/path&gt;&lt;/svg&gt;&lt;/span&gt;, &lt;span&gt;&lt;span&gt;https://github.com/JA4S/JANC&lt;/span&gt;&lt;svg&gt;&lt;path&gt;&lt;/path&gt;&lt;/svg&gt;&lt;/span&gt;&lt;/div&gt;&lt;div&gt;Licensing provisions: MIT&lt;/div&gt;&lt;div&gt;Programming language: Python&lt;/div&gt;&lt;div&gt;&lt;em&gt;Nature of problem&lt;/em&gt;: The numerical solution of compressible reactive flows plays a crucial role in combustion, energy utilization, and the design and manufacturing of propulsion systems. However, the multi-species nature, highly transient behavior, and strong numerical stiffness of reactive flows lead to significantly higher computational costs compared to conventional flow problems. In addition, conventional reactive flow solvers are typically built on Fortran or C++ frameworks, making them difficult to integrate with data-driven methods based on existing Python ecosystems—particularly gradient-based optimization techniques such as machine learning","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"319 ","pages":"Article 109915"},"PeriodicalIF":3.4,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145517651","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Particle collisionality in scaled kinetic plasma simulations 尺度动力学等离子体模拟中的粒子碰撞
IF 3.4 2区 物理与天体物理 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2026-02-01 Epub Date: 2025-10-18 DOI: 10.1016/j.cpc.2025.109899
S.R. Totorica , K.V. Lezhnin , W. Fox
Kinetic plasma processes, such as magnetic reconnection, collisionless shocks, and turbulence, are fundamental to the dynamics of astrophysical and laboratory plasmas. Simulating these processes often requires particle-in-cell (PIC) methods, but the computational cost of fully kinetic simulations can necessitate the use of artificial parameters, such as a reduced speed of light and ion-to-electron mass ratio, to decrease expense. While these approximations can preserve overall dynamics under specific conditions, they introduce nontrivial impacts on particle collisionality that are not yet well understood. In this work, we develop a method to scale particle collisionality in simulations employing an artificial speed of light and/or an artificial ion-to-electron mass ratio. By introducing species-dependent scaling factors, we independently adjust inter- and intra-species collision rates to better replicate the collisional properties of the physical system. Our approach maintains the fidelity of electron and ion transport properties while preserving critical relaxation rates, such as energy exchange timescales, within the limits of weakly collisional plasma theory. We demonstrate the accuracy of this scaling method through benchmarking tests against theoretical relaxation rates and connecting to fluid theory, highlighting its ability to retain key transport properties. Existing collisional PIC implementations can be easily modified to include this scaling, which will enable deeper insights into the behavior of marginally collisional plasmas across various contexts.
动态等离子体过程,如磁重联、无碰撞冲击和湍流,是天体物理和实验室等离子体动力学的基础。模拟这些过程通常需要粒子池(PIC)方法,但是完全动力学模拟的计算成本可能需要使用人工参数,例如降低光速和离子与电子的质量比,以降低成本。虽然这些近似可以在特定条件下保持整体动力学,但它们对粒子碰撞产生的重要影响尚未得到很好的理解。在这项工作中,我们开发了一种在模拟中使用人工光速和/或人工离子与电子质量比来衡量粒子碰撞的方法。通过引入依赖物种的比例因子,我们独立调整物种间和物种内的碰撞率,以更好地复制物理系统的碰撞特性。我们的方法保持了电子和离子输运性质的保真度,同时在弱碰撞等离子体理论的限制内保持了临界弛豫率,如能量交换时间尺度。我们通过对理论松弛率和流体理论进行基准测试,证明了这种标度方法的准确性,强调了其保留关键输运性质的能力。现有的碰撞PIC实现可以很容易地修改,以包括这种缩放,这将使人们能够更深入地了解各种环境下边缘碰撞等离子体的行为。
{"title":"Particle collisionality in scaled kinetic plasma simulations","authors":"S.R. Totorica ,&nbsp;K.V. Lezhnin ,&nbsp;W. Fox","doi":"10.1016/j.cpc.2025.109899","DOIUrl":"10.1016/j.cpc.2025.109899","url":null,"abstract":"<div><div>Kinetic plasma processes, such as magnetic reconnection, collisionless shocks, and turbulence, are fundamental to the dynamics of astrophysical and laboratory plasmas. Simulating these processes often requires particle-in-cell (PIC) methods, but the computational cost of fully kinetic simulations can necessitate the use of artificial parameters, such as a reduced speed of light and ion-to-electron mass ratio, to decrease expense. While these approximations can preserve overall dynamics under specific conditions, they introduce nontrivial impacts on particle collisionality that are not yet well understood. In this work, we develop a method to scale particle collisionality in simulations employing an artificial speed of light and/or an artificial ion-to-electron mass ratio. By introducing species-dependent scaling factors, we independently adjust inter- and intra-species collision rates to better replicate the collisional properties of the physical system. Our approach maintains the fidelity of electron and ion transport properties while preserving critical relaxation rates, such as energy exchange timescales, within the limits of weakly collisional plasma theory. We demonstrate the accuracy of this scaling method through benchmarking tests against theoretical relaxation rates and connecting to fluid theory, highlighting its ability to retain key transport properties. Existing collisional PIC implementations can be easily modified to include this scaling, which will enable deeper insights into the behavior of marginally collisional plasmas across various contexts.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"319 ","pages":"Article 109899"},"PeriodicalIF":3.4,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145464292","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Acceleration of the CASINO quantum Monte Carlo software using graphics processing units and OpenACC 使用图形处理单元和OpenACC加速CASINO量子蒙特卡罗软件
IF 3.4 2区 物理与天体物理 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2026-02-01 Epub Date: 2025-10-26 DOI: 10.1016/j.cpc.2025.109900
B. Thorpe , M.J. Smith , P.J. Hasnip , N.D. Drummond
We describe how quantum Monte Carlo calculations using the CASINO software can be accelerated using graphics processing units (GPUs) and OpenACC. In particular we consider offloading Ewald summation, the evaluation of long-range two-body terms in the Jastrow correlation factor, and the evaluation of orbitals in a blip basis set. We present results for three- and two-dimensional homogeneous electron gases and ab initio simulations of bulk materials, showing that significant speedups of up to a factor of 2.5 can be achieved by the use of GPUs when several hundred particles are included in the simulations. The use of single-precision arithmetic can improve the speedup further without significant detriment to the accuracy of the calculations.
我们描述了使用CASINO软件的量子蒙特卡罗计算如何使用图形处理单元(gpu)和OpenACC加速。我们特别考虑了卸载Ewald和,Jastrow相关因子中远程二体项的评价,以及在一个小点基集中轨道的评价。我们给出了三维和二维均质电子气体和块状材料从头算模拟的结果,表明当模拟中包含数百个粒子时,使用gpu可以实现高达2.5倍的显著加速。使用单精度算法可以在不显著损害计算精度的情况下进一步提高加速。
{"title":"Acceleration of the CASINO quantum Monte Carlo software using graphics processing units and OpenACC","authors":"B. Thorpe ,&nbsp;M.J. Smith ,&nbsp;P.J. Hasnip ,&nbsp;N.D. Drummond","doi":"10.1016/j.cpc.2025.109900","DOIUrl":"10.1016/j.cpc.2025.109900","url":null,"abstract":"<div><div>We describe how quantum Monte Carlo calculations using the CASINO software can be accelerated using graphics processing units (GPUs) and OpenACC. In particular we consider offloading Ewald summation, the evaluation of long-range two-body terms in the Jastrow correlation factor, and the evaluation of orbitals in a blip basis set. We present results for three- and two-dimensional homogeneous electron gases and <em>ab initio</em> simulations of bulk materials, showing that significant speedups of up to a factor of 2.5 can be achieved by the use of GPUs when several hundred particles are included in the simulations. The use of single-precision arithmetic can improve the speedup further without significant detriment to the accuracy of the calculations.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"319 ","pages":"Article 109900"},"PeriodicalIF":3.4,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145464293","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Reconstructing Laurent expansion of rational functions using p-adic numbers 用p进数重建有理函数的劳伦展开式
IF 3.4 2区 物理与天体物理 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2026-02-01 Epub Date: 2025-10-25 DOI: 10.1016/j.cpc.2025.109908
Tianya Xia, Li Lin Yang
We propose a novel method for reconstructing Laurent expansion of rational functions using p-adic numbers. By evaluating the rational functions in p-adic fields rather than finite fields, it is possible to probe the expansion coefficients simultaneously, enabling their reconstruction from a single set of evaluations. Compared with the reconstruction of the full expression, constructing the Laurent expansion to the first few orders significantly reduces the required computational resources. Our method can handle expansions with respect to more than one variables simultaneously. Among possible applications, we anticipate that our method can be used to simplify the integration-by-parts reduction of Feynman integrals in cutting-edge calculations.
PROGRAM SUMMARY Manuscript Title: Reconstructing Laurent expansion of rational functions using p-adic numbers
Authors: Tianya Xia, Li Lin Yang
Program Title: LaurentExpPadicReconstruct
CPC Library link to program files: (to be added by Technical Editor)
Licensing provisions: GPLv3
Programming language: C++
External routines/libraries: FireFly, FLINT
Nature of problem: Reconstructing Laurent expansion of rational function arising in the IBP reuduction of Feynman Integrals.
Solution method: Uses p-adic numbers combined with rational function reconstruction over finite fields.
Running time: Typically ranges from several minutes to a few hours, depending on the size and algebraic complexity of the input.
提出了一种利用p进数重构有理函数劳伦展开式的新方法。通过计算p进域中的有理函数而不是有限域中的有理函数,可以同时探测展开系数,从而使它们能够从一组计算中重建。与完整表达式的重构相比,将Laurent展开构造到前几阶,大大减少了所需的计算资源。我们的方法可以同时处理多个变量的展开。在可能的应用中,我们期望我们的方法可以用于简化尖端计算中费曼积分的分部积分缩减。程序摘要手稿标题:用p进数重建有理函数的Laurent展开作者:夏天亚,杨李林程序标题:LaurentExpPadicReconstructCPC库链接到程序文件:(由技术编辑添加)许可条款:gplv3编程语言:c++外部例程/库:FireFly, flint问题性质:重建费曼积分IBP约简中有理函数的Laurent展开。求解方法:利用p进数结合有限域上的有理函数重构。运行时间:通常从几分钟到几个小时不等,取决于输入的大小和代数复杂度。
{"title":"Reconstructing Laurent expansion of rational functions using p-adic numbers","authors":"Tianya Xia,&nbsp;Li Lin Yang","doi":"10.1016/j.cpc.2025.109908","DOIUrl":"10.1016/j.cpc.2025.109908","url":null,"abstract":"<div><div>We propose a novel method for reconstructing Laurent expansion of rational functions using <span><math><mi>p</mi></math></span>-adic numbers. By evaluating the rational functions in <span><math><mi>p</mi></math></span>-adic fields rather than finite fields, it is possible to probe the expansion coefficients simultaneously, enabling their reconstruction from a single set of evaluations. Compared with the reconstruction of the full expression, constructing the Laurent expansion to the first few orders significantly reduces the required computational resources. Our method can handle expansions with respect to more than one variables simultaneously. Among possible applications, we anticipate that our method can be used to simplify the integration-by-parts reduction of Feynman integrals in cutting-edge calculations.</div><div><strong>PROGRAM SUMMARY</strong> <em>Manuscript Title:</em> Reconstructing Laurent expansion of rational functions using p-adic numbers</div><div><em>Authors:</em> Tianya Xia, Li Lin Yang</div><div><em>Program Title:</em> LaurentExpPadicReconstruct</div><div><em>CPC Library link to program files:</em> (to be added by Technical Editor)</div><div><em>Licensing provisions:</em> GPLv3</div><div><em>Programming language:</em> C++</div><div><em>External routines/libraries:</em> FireFly, FLINT</div><div><em>Nature of problem:</em> Reconstructing Laurent expansion of rational function arising in the IBP reuduction of Feynman Integrals.</div><div><em>Solution method:</em> Uses p-adic numbers combined with rational function reconstruction over finite fields.</div><div><em>Running time:</em> Typically ranges from several minutes to a few hours, depending on the size and algebraic complexity of the input.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"319 ","pages":"Article 109908"},"PeriodicalIF":3.4,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145464407","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Portal for high-precision atomic data and computation 用于高精度原子数据和计算的门户
IF 3.4 2区 物理与天体物理 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2026-02-01 Epub Date: 2025-11-14 DOI: 10.1016/j.cpc.2025.109951
Amani Kiruga , Charles Cheung , Dmytro Filin , Parinaz Barakhshan , Akshay Bhosale , Vipul Badhan , Bindiya Arora , Rudolf Eigenmann , Marianna S. Safronova
We’ve developed a scalable and sustainable online atomic data portal with an automated interface for easy update and addition of new data. The current portal provides energies, transition matrix elements, transition rates, radiative lifetimes, branching ratios, polarizabilities, hyperfine constants, and other data, for 28 atoms and ions. It also features an interactive polarizability plotting interface for neutral atoms and singly-charged ions. The data production is supported by recent developments of open-access atomic software based on our research codes, including new workflow algorithms, which allow large volumes of such data to be generated with automated accuracy assessments. This entails a new method of comparing our calculated values with data from the NIST Atomic Spectra Database. All calculated values include estimated uncertainties. Data for more systems will be added in the future. Experimental values are included with references, where high-precision data are available.
我们开发了一个可扩展的、可持续的在线原子数据门户,它带有一个自动界面,可以方便地更新和添加新数据。当前门户提供28个原子和离子的能量、跃迁矩阵元素、跃迁速率、辐射寿命、分支比、极化率、超精细常数和其他数据。它还具有中性原子和单电荷离子的交互极化率绘图界面。数据生产由基于我们研究代码的开放存取原子软件的最新发展提供支持,包括新的工作流程算法,它允许通过自动准确性评估生成大量此类数据。这需要一种新的方法来比较我们的计算值与NIST原子光谱数据库的数据。所有计算值都包括估计的不确定性。将来会添加更多系统的数据。参考文献中包含有高精度数据的实验值。
{"title":"Portal for high-precision atomic data and computation","authors":"Amani Kiruga ,&nbsp;Charles Cheung ,&nbsp;Dmytro Filin ,&nbsp;Parinaz Barakhshan ,&nbsp;Akshay Bhosale ,&nbsp;Vipul Badhan ,&nbsp;Bindiya Arora ,&nbsp;Rudolf Eigenmann ,&nbsp;Marianna S. Safronova","doi":"10.1016/j.cpc.2025.109951","DOIUrl":"10.1016/j.cpc.2025.109951","url":null,"abstract":"<div><div>We’ve developed a scalable and sustainable online atomic data portal with an automated interface for easy update and addition of new data. The current portal provides energies, transition matrix elements, transition rates, radiative lifetimes, branching ratios, polarizabilities, hyperfine constants, and other data, for 28 atoms and ions. It also features an interactive polarizability plotting interface for neutral atoms and singly-charged ions. The data production is supported by recent developments of open-access atomic software based on our research codes, including new workflow algorithms, which allow large volumes of such data to be generated with automated accuracy assessments. This entails a new method of comparing our calculated values with data from the NIST Atomic Spectra Database. All calculated values include estimated uncertainties. Data for more systems will be added in the future. Experimental values are included with references, where high-precision data are available.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"319 ","pages":"Article 109951"},"PeriodicalIF":3.4,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145576712","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
期刊
Computer Physics Communications
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1