Abstract
In recent works, the full set of three-loop diagrams (of order \(\alpha _{s}^{2}\) ) associated with the current-current operators O1 and O2 contributing to the decay amplitude for b → sγ were calculated. In this paper, one pair of these diagrams is calculated at various values of the charm-quark mass mc using different methods, and the numerical accuracy/ease of use of these methods are compared to the analytic results obtained by others. Using the programs AMFlow and DiffExp to solve the differential equations for the master integrals, precise numerical results were obtained across a range of charm-quark mass values. In addition, asymptotic and Taylor series expansions were performed around z = \({{m_{c}^{2}} \mathord{\left/ {\vphantom {{m_{c}^{2}} {m_{b}^{2}}}} \right. \kern-0em} {m_{b}^{2}}}\) = 0 and z = 1/10, respectively. A fully numerical evaluation of the master integrals using PySecDec was also carried out.
