MARTY - Modern ARtificial Theoretical phYsicist A C++ framework automating theoretical calculations Beyond the Standard Model

Abstract Studies Beyond the Standard Model (BSM) will become more and more important in the near future with the rapidly increasing amount of data from different experiments around the world. The full study of BSM models is in general an extremely time-consuming task involving long and difficult calculations. It is in practice not possible to do exhaustive predictions in these models by hand. Here we present MARTY, a new C++ framework that fully automates calculations from the Lagrangian to physical quantities such as amplitudes or cross-sections. It can fully simplify, automatically and symbolically, physical quantities in a very large variety of models and compute Wilson coefficients in effective theories. This will considerably facilitate BSM studies in flavour physics. Contrary to the existing public codes in this field MARTY aims at providing a unique, free, open-source, powerful and user-friendly tool for high-energy physicists studying predictive BSM models, in effective or full theories up to the one-loop level, which does not rely on any external package. With a few lines of code one can gather final expressions that may be evaluated numerically for statistical analysis. Program summary Program Title: MARTY CPC Library link to program files: https://doi.org/10.17632/y7n82jcj8d.1 Developer’s repository link: https://marty.in2p3.fr , https://fr.overleaf.com/project/5fb7c30c3175f910d2392059 Licensing provisions: GPLv3 Programming language: C++ Supplementary material: GNU-GSL [1], LoopTools [2] Nature of problem: From a given Lagrangian, Beyond the Standard Model phenomenology often requires one-loop theoretical calculations. Those calculations must be done analytically because a numerical evaluation of such quantities is not possible, the number of terms being too large for any computer to handle. Simplifications must be done analytically, and are very long and error prone justifying to automate them. This is why we need a computer algebra system dedicated to this task. Some packages written with Mathematica [3] - a private and closed software for symbolic manipulations - implement those calculations. There are multiple packages, the user needs to pay for Mathematica and the calculation of Wilson coefficients for BSM and at the one-loop level is actually hard to automate. Solution method: MARTY is a code written for BSM phenomenology that comes with its own computer algebra system, CSL, and automates the calculation of theoretical quantities at the one-loop level, in a very large variety of BSM model. Physicists now can have a unique C++ code, free and open-source, that can be used to do model building BSM and calculate any kind of amplitude, squared amplitude or Wilson coefficient BSM from the Lagrangian. From the developers point of view, having such a code, independent of any other framework using only the C++ standard library (C++17 standard), is a great opportunity to develop it even further, implementing new simplifications procedures, new calculation procedures and more. References: [1] GSL - GNU Scientific Library. https://www.gnu.org/software/gsl/ . [2] T. Hahn, M. Perez-Victoria, Automatized one loop calculations in four-dimensions and D-dimensions, Comput. Phys. Commun. 118 (1999) 153–165. arXiv:hep-ph/9807565 , doi: http://dx.doi.org/10.1016/S0010-4655(98)00173-8 . [3] Wolfram Research Inc., Mathematica. https://www.wolfram.com/mathematic .

[1]  Florian Staub,et al.  Exploring New Models in All Detail with SARAH , 2015, 1503.04200.

[2]  Matthew D. Schwartz,et al.  Quantum Field Theory and the Standard Model , 2013 .

[3]  A closed expression for the UV-divergent parts of one-loop tensor integrals in dimensional regularization , 2006, hep-ph/0609282.

[4]  G. Zanderighi,et al.  One-loop calculations in quantum field theory: From Feynman diagrams to unitarity cuts , 2011, 1105.4319.

[5]  Claude Duhr,et al.  UFO - The Universal FeynRules Output , 2011, Comput. Phys. Commun..

[6]  David Shih,et al.  FormFlavor Manual , 2016, 1606.00003.

[7]  A. Denner,et al.  Compact Feynman rules for Majorana fermions , 1992 .

[8]  A. Denner Techniques for the Calculation of Electroweak Radiative Corrections at the One‐Loop Level and Results for W‐physics at LEP 200 , 2007, 0709.1075.

[9]  Farvah Mahmoudi,et al.  SuperIso: A program for calculating the isospin asymmetry of B -> K*gamma in the MSSM , 2007, Comput. Phys. Commun..

[10]  A. Sherstnev,et al.  CompHEP 4.4 - Automatic Computations from Lagrangians to Events , 2004, hep-ph/0403113.

[11]  T. V. Ritbergen,et al.  GROUP THEORY FACTORS FOR FEYNMAN DIAGRAMS , 1998, hep-ph/9802376.

[12]  Benjamin C. Allanach,et al.  The inclusion of two-loop SUSYQCD corrections to gluino and squark pole masses in the minimal and next-to-minimal supersymmetric standard model: SOFTSUSY3.7 , 2016, Comput. Phys. Commun..

[13]  T. Hahn,et al.  Automatized One-Loop Calculations in 4 and D dimensions , 1998 .

[14]  Farvah Mahmoudi,et al.  SuperIso v3.0, flavor physics observables calculations: extension to NMSSM , 2009, Comput. Phys. Commun..

[15]  Claude Duhr,et al.  FeynRules 2.0 - A complete toolbox for tree-level phenomenology , 2013, Comput. Phys. Commun..

[16]  P. Cvitanović Group theory for Feynman diagrams in non-Abelian gauge theories , 1976 .

[17]  B. C. Allanach,et al.  SOFTSUSY: A program for calculating supersymmetric spectra☆ , 2001, hep-ph/0104145.

[18]  Farvah Mahmoudi,et al.  SuperIso Relic: A program for calculating relic density and flavor physics observables in Supersymmetry , 2009, Comput. Phys. Commun..

[19]  M. Ciuchini,et al.  Next-to-leading QCD corrections to B → Xsγ in supersymmetry , 1998, hep-ph/9806308.

[20]  Werner Porod,et al.  A flavor kit for BSM models , 2014, 1405.1434.

[21]  Mikhail Dubinin,et al.  CompHEP - a package for evaluation of Feynman diagrams and integration over multi-particle phase space. User's manual for version 33 , 1999, hep-ph/9908288.

[22]  Farvah Mahmoudi,et al.  SuperIso v2.3: A program for calculating flavor physics observables in supersymmetry , 2008, Comput. Phys. Commun..

[23]  G. C. Wick The Evaluation of the Collision Matrix , 1950 .

[24]  Alexandre Arbey,et al.  SuperIso Relic v3.0: A program for calculating relic density and flavour physics observables: Extension to NMSSM , 2011, Comput. Phys. Commun..

[25]  A. V. Semenov,et al.  LanHEP - a package for the automatic generation of Feynman rules in field theory. Version 3.0 , 2002, Comput. Phys. Commun..

[26]  Farvah Mahmoudi,et al.  SuperIso Relic v4: A program for calculating dark matter and flavour physics observables in supersymmetry , 2018, Comput. Phys. Commun..

[27]  Andrzej J. Buras Weak Hamiltonian, CP Violation and Rare Decays , 1998 .