# Generic calculation of two-body partial decay widths at the full one-loop level

## TL;DR

This paper presents a comprehensive framework for calculating two-body partial decay widths at the full one-loop level across various models, incorporating virtual, real, and higher-order corrections with IR divergence cancellation.

## Contribution

It introduces a fully generic implementation in SARAH and SPheno for one-loop decay width calculations, supporting diverse decay channels and model-specific renormalization schemes.

## Key findings

- Validated against existing codes for MSSM decay processes.
- Reproduced results for R-parity violating neutralino decays.
- Compared and contrasted with other public decay codes.

## Abstract

We describe a fully generic implementation of two-body partial decay widths at the full one-loop level in the SARAH and SPheno framework compatible with most supported models. It incorporates fermionic decays to a fermion and a scalar or a gauge boson as well as scalar decays into two fermions, two gauge bosons, two scalars or a scalar and a gauge boson. We present the relevant generic expressions for virtual and real corrections. Whereas wavefunction corrections are determined from on-shell conditions, the parameters of the underlying model are by default renormalised in a DR (or MS) scheme. However, the user can also define model-specific counter-terms. As an example we discuss the renormalisation of the electric charge in the Thomson limit for top-quark decays in the standard model. One-loop induced decays are also supported. The framework additionally allows the addition of mass and mixing corrections induced at higher orders for the involved external states. We explain our procedure to cancel infra-red divergences for such cases, which is achieved through an infra-red counter-term taking into account corrected Goldstone boson vertices. We compare our results for sfermion, gluino and Higgs decays in the minimal supersymmetric standard model (MSSM) against the public codes SFOLD, FVSFOLD and HFOLD and explain observed differences. Radiative induced gluino and neutralino decays are compared against the original implementation in SPheno in the MSSM. We exactly reproduce the results of the code CNNDecays for decays of neutralinos and charginos in R-parity violating models. The new version SARAH 4.11.0 by default includes the calculation of two-body decay widths at the full one-loop level. Current limitations for certain model classes are described.