[论文解读] New Developments in FormCalc 4.1

Thomas Hahn|ArXiv.org|Jun 21, 2005

Parallel Computing and Optimization Techniques被引用 8

一句话总结

FormCalc 4.1 引入了高能物理中自动一阶矩阵元计算的重大改进，包括用于高效费米子振幅计算的外尔–范德瓦耳登旋量形式体系、与Cuba库集成以实现多维相空间积分、优化的Fortran代码生成，以及与FeynHiggs和MSSM中非最小味违反的新型接口。这些改进显著提升了性能、数值稳定性与可用性，适用于Higgs粒子产生和顶夸克衰变等过程的高精度计算。

ABSTRACT

FormCalc is a matrix-element generator that turns FeynArts amplitudes up to one loop into a Fortran code for computing the squared matrix element. The generated code can be run with FormCalc's own driver programs or used with other `frontends', e.g. Monte Carlos. Major new or enhanced features in Version 4.1 are: treatment of external fermions, phase-space integration, code-generation functions, extensions for the MSSM, the HadCalc frontend.

研究动机与目标

• 提高涉及外部费米子的一阶矩阵元计算的效率与数值稳定性。
• 通过与FeynHiggs 2.2集成，实现在MSSM中对Higgs粒子质量的高精度计算。
• 通过支持上型与下型标夸克的完整6×6混合矩阵，实现MSSM中非最小味违反（NMFV）的系统化处理。
• 通过公开的Fortran代码生成函数与新型前端界面，简化代码生成与集成工作流。
• 通过使用Cuba 1.2库，结合先进的蒙特卡洛与确定性算法，增强相空间积分能力。

提出的方法

• 在FormCalc中实现外尔–范德瓦耳登旋量形式体系，直接计算振幅而不使用迹技术，将n个费米子链的计算复杂度从O(n²)降低至O(n)。
• 集成Cuba 1.2库用于多维相空间积分，支持Vegas、Suave、Divonne与Cuhre算法，具备梅森尼采样与网格记忆化等特性。
• 引入公开的Mathematica函数，用于生成经过优化的模块化Fortran代码，包括公共子表达式消除与大表达式分拆至多个文件的支持。
• 通过预处理器标志启用FeynHiggs 2.2集成，以实现具有完整一阶精度的辐射修正Higgs质量计算。
• 通过扩展模型文件与数值计算，支持MSSM中NMFV，实现对上型与下型标夸克的完整6×6混合矩阵。
• 提供新型HadCalc前端，用于自动计算强子截面，支持PDF卷积、运动学截断及不变质量、快度与横动量的微分分布。

实验结果

研究问题

• RQ1 如何在不牺牲精度的前提下，降低包含外部费米子的平方矩阵元的计算成本？
• RQ2 使用先进的蒙特卡洛与确定性算法，相空间积分的性能与稳定性可实现哪些改进？
• RQ3 如何对大型表达式与模块化编译进行高能物理计算中的自动Fortran代码生成优化？
• RQ4 通过FeynHiggs实现完整一阶Higgs质量修正，对MSSM中高精度预测有何影响？
• RQ5 如何系统性地将MSSM中的非最小味违反（NMFV）整合到自动矩阵元生成器中？

主要发现

• 外尔–范德瓦耳登形式体系将所需振幅收缩次数从O(n²)减少至O(n)，实现了对外部费米子数量较多的振幅的高效计算。
• Cuba 1.2集成结合梅森尼采样与网格记忆化，显著提升了收敛性与鲁棒性，尤其在高维相空间中表现更优。
• 新型Fortran代码生成函数可生成经过优化的模块化代码，支持公共子表达式消除，并可将大型表达式分拆至多个文件。
• 与FeynHiggs 2.2的集成支持精确的一阶Higgs质量计算，用户可选择使用最精确的FeynHiggsFast近似或完整一阶结果。
• 通过预处理器标志与扩展模型文件，MSSM中的NMFV现已支持，可访问完整的6×6标夸克混合矩阵。
• HadCalc前端可实现自动、高精度的强子截面计算，支持PDF卷积、运动学截断及微分分布，相关功能可应要求提供。