Two-loop tensor integral coefficients in OpenLoops

Stefano Pozzorini (Physik-Institut, Universität Zürich, Winterthurerstrasse 190, Zürich, CH-8057, Switzerland) ; Natalie Schär (Paul Scherrer Institut, Forschungsstrasse 111, Villigen PSI, CH-5232, Switzerland) ; Max Zoller (Paul Scherrer Institut, Forschungsstrasse 111, Villigen PSI, CH-5232, Switzerland)

We present a new and fully general algorithm for the automated construction of the integrands of two-loop scattering amplitudes. This is achieved through a generalisation of the open-loops method to two loops. The core of the algorithm consists of a numerical recursion, where the various building blocks of two-loop diagrams are connected to each other through process-independent operations that depend only on the Feynman rules of the model at hand. This recursion is implemented in terms of tensor coefficients that encode the polynomial dependence of loop numerators on the two independent loop momenta. The resulting coefficients are ready to be combined with corresponding tensor integrals to form scattering probability densities at two loops. To optimise CPU efficiency we have compared several algorithmic options identifying one that outperforms naive solutions by two orders of magnitude. This new algorithm is implemented in the OpenLoops framework in a fully automated way for two-loop QED and QCD corrections to any Standard Model process. The technical performance is discussed in detail for several 2 → 2 and 2 → 3 processes with up to order 105 two-loop diagrams. We find that the CPU cost scales linearly with the number of two-loop diagrams and is comparable to the cost of corresponding real-virtual ingredients in a NNLO calculation. This new algorithm constitutes a key building block for the construction of an automated generator of scattering amplitudes at two loops.

{
  "_oai": {
    "updated": "2022-10-08T00:32:17Z", 
    "id": "oai:repo.scoap3.org:70144", 
    "sets": [
      "JHEP"
    ]
  }, 
  "authors": [
    {
      "affiliations": [
        {
          "country": "Switzerland", 
          "value": "Physik-Institut, Universit\u00e4t Z\u00fcrich, Winterthurerstrasse 190, Z\u00fcrich, CH-8057, Switzerland", 
          "organization": "Universit\u00e4t Z\u00fcrich"
        }
      ], 
      "surname": "Pozzorini", 
      "email": "pozzorin@physik.uzh.ch", 
      "full_name": "Pozzorini, Stefano", 
      "given_names": "Stefano"
    }, 
    {
      "affiliations": [
        {
          "country": "Switzerland", 
          "value": "Paul Scherrer Institut, Forschungsstrasse 111, Villigen PSI, CH-5232, Switzerland", 
          "organization": "Paul Scherrer Institut"
        }
      ], 
      "surname": "Sch\u00e4r", 
      "email": "natalie.schaer@psi.ch", 
      "full_name": "Sch\u00e4r, Natalie", 
      "given_names": "Natalie"
    }, 
    {
      "affiliations": [
        {
          "country": "Switzerland", 
          "value": "Paul Scherrer Institut, Forschungsstrasse 111, Villigen PSI, CH-5232, Switzerland", 
          "organization": "Paul Scherrer Institut"
        }
      ], 
      "surname": "Zoller", 
      "email": "max.zoller@psi.ch", 
      "full_name": "Zoller, Max", 
      "given_names": "Max"
    }
  ], 
  "titles": [
    {
      "source": "Springer", 
      "title": "Two-loop tensor integral coefficients in OpenLoops"
    }
  ], 
  "dois": [
    {
      "value": "10.1007/JHEP05(2022)161"
    }
  ], 
  "publication_info": [
    {
      "page_end": "48", 
      "journal_title": "Journal of High Energy Physics", 
      "material": "article", 
      "journal_volume": "2022", 
      "artid": "JHEP05(2022)161", 
      "year": 2022, 
      "page_start": "1", 
      "journal_issue": "5"
    }
  ], 
  "$schema": "http://repo.scoap3.org/schemas/hep.json", 
  "acquisition_source": {
    "date": "2022-10-08T00:31:17.856573", 
    "source": "Springer", 
    "method": "Springer", 
    "submission_number": "57da281646a011ed91f9a2b26b995c7b"
  }, 
  "page_nr": [
    48
  ], 
  "license": [
    {
      "url": "https://creativecommons.org/licenses//by/4.0", 
      "license": "CC-BY-4.0"
    }
  ], 
  "copyright": [
    {
      "holder": "The Author(s)", 
      "year": "2022"
    }
  ], 
  "control_number": "70144", 
  "record_creation_date": "2022-05-25T06:30:18.169503", 
  "_files": [
    {
      "checksum": "md5:158037f59339ef9e95923308c1cbdd0f", 
      "filetype": "xml", 
      "bucket": "b859577c-b841-4a5c-bb3e-2b03bbde14a5", 
      "version_id": "48a1d382-ecb5-488a-a2c7-f7cad2248f47", 
      "key": "10.1007/JHEP05(2022)161.xml", 
      "size": 12131
    }, 
    {
      "checksum": "md5:b9e7e6cf664d564d3104c5d66f13500f", 
      "filetype": "pdf/a", 
      "bucket": "b859577c-b841-4a5c-bb3e-2b03bbde14a5", 
      "version_id": "70c82537-eac3-4ec0-81bd-a68d3dd9a42d", 
      "key": "10.1007/JHEP05(2022)161_a.pdf", 
      "size": 889592
    }
  ], 
  "collections": [
    {
      "primary": "Journal of High Energy Physics"
    }
  ], 
  "arxiv_eprints": [
    {
      "categories": [
        "hep-ph"
      ], 
      "value": "2201.11615"
    }
  ], 
  "abstracts": [
    {
      "source": "Springer", 
      "value": "We present a new and fully general algorithm for the automated construction of the integrands of two-loop scattering amplitudes. This is achieved through a generalisation of the open-loops method to two loops. The core of the algorithm consists of a numerical recursion, where the various building blocks of two-loop diagrams are connected to each other through process-independent operations that depend only on the Feynman rules of the model at hand. This recursion is implemented in terms of tensor coefficients that encode the polynomial dependence of loop numerators on the two independent loop momenta. The resulting coefficients are ready to be combined with corresponding tensor integrals to form scattering probability densities at two loops. To optimise CPU efficiency we have compared several algorithmic options identifying one that outperforms naive solutions by two orders of magnitude. This new algorithm is implemented in the OpenLoops framework in a fully automated way for two-loop QED and QCD corrections to any Standard Model process. The technical performance is discussed in detail for several 2 \u2192 2 and 2 \u2192 3 processes with up to order 105 two-loop diagrams. We find that the CPU cost scales linearly with the number of two-loop diagrams and is comparable to the cost of corresponding real-virtual ingredients in a NNLO calculation. This new algorithm constitutes a key building block for the construction of an automated generator of scattering amplitudes at two loops."
    }
  ], 
  "imprints": [
    {
      "date": "2022-05-24", 
      "publisher": "Springer"
    }
  ]
}
Published on:
24 May 2022
Publisher:
Springer
Published in:
Journal of High Energy Physics , Volume 2022 (2022)
Issue 5
Pages 1-48
DOI:
https://doi.org/10.1007/JHEP05(2022)161
arXiv:
2201.11615
Copyrights:
The Author(s)
Licence:
CC-BY-4.0

Fulltext files: