op-closure-estimates

Operator closure estimates for proton decay (and neutrino masses)

https://github.com/johngarg/op-closure-estimates

Repository

Operator closure estimates for proton decay (and neutrino masses)

readme.org

#+title: Operator closure estimates for BNV nucleon decays

Our code contains the Mathematica code that implements our symbolic rewriting
system as well as the Python code in which we perform much of our analysis. If
you use any of the code contained in this repository, please cite our paper
"Model-independent estimates for loop-induced baryon-number-violating nucleon
decays" https://arxiv.org/abs/2401.04768.

* Rewriting system

Load up the definitions in the core files:
#+begin_src wolfram
Get["main.m"];
Get["operators.m"];
Get["rules.m"];
Get["decay.m"];
Get["pmatch.m"];
Get["symmetries.m"];
#+end_src

The main function of the package is =MatchingData=, which is a wrapper around a
number of auxiliary functions. =MatchingData= returns the matching estimates
associated with a specific field-string operator. Its inputs are the operator
label, provided as a string, and the flavour indices, a list of integers. For example
#+begin_src wolfram
MatchingData["11", {1, 1, 1, 1}]
#+end_src
The output is in the mass basis.

* Matching-estimate database

Our results can be accessed through the =Pandas= dataframes we provide in
=database.py=. We provide an example notebook showing how to access these data
in the =examples/= directory. We summarise the keys of the loop-level dataframe
below.

| Key                         | Description                                                                               |
|-----------------------------+-------------------------------------------------------------------------------------------|
| =smeft_op=                  | The `sympy` symbol representing the coefficient of the SMEFT operator being matched on to |
| =smeft_label=               | The label of the SMEFT operator                                                           |
| =smeft_flavour=             | The flavour indices of the SMEFT operator                                                 |
| =smeft_op_expr=             | The symbolic expression implementing the SMEFT-LEFT matching provided in the paper        |
| =process=                   | The decay process generated in LaTeX                                                      |
| =gamma=                     | The symbolic decay rate                                                                   |
| =left_dimension=            | The dimension of the LEFT operator generated                                              |
| =gamma_coeff_1=             | The decay rate with the dimensionless coefficient set to 1                                |
| =left_op=                   | The label of the LEFT operator driving the decay process                                  |
| =left_flavour=              | The flavour indices of the LEFT operator                                                  |
| =lambda_limit=              | The symbolic lower limit on the underlying new-physics scale                              |
| =lambda_limit_coeff_1=      | The numeric lower limit on the underlying new-physics scale with coefficient set to 1     |
| =lifetime_limit=            | The experimental lifetime limit in years                                                  |
| =lifetime_limit_ref=        | The reference from which the experimental bound is taken                                  |
| =matrix_elem=               | The nuclear matrix element relevant to the nucleon decay calculation                      |
| =fieldstring_label=         | The label of the field-string operator                                                    |
| =fieldstring_flavour=       | The flavour indices of the field string                                                   |
| =latex=                     | The LaTeX expression of the matching estimate                                             |
| =gamma_fieldstring_coeff=   | The decay rate with the matching estimate pasted in for the SMEFT coefficient             |
| =gamma_fieldstring_coeff_1= | Same as `gamma_fieldstring_coeff` but for field-string coefficient set to 1               |
| =is_future=                 | True if the row corresponds to a future sensitivity estimate, False otherwise             |

Owner

Citation (CITATION.cff)

cff-version: 1.2.0
message: "If you use this software, please cite arXiv:2401.XXXXX."
authors:
- family-names: "Gargalionis"
  given-names: "John"
  orcid: "https://orcid.org/0000-0002-0745-7121"
- family-names: "Herrero-Garcia"
  given-names: "Juan"
  orcid: "https://orcid.org/0000-0002-3300-0029"
- family-names: "Schmidt"
  given-names: "Michael"
  orcid: "https://orcid.org/0000-0002-8792-5537"
title: "Operator closure estimates for BNV nucleon decays"
version: 1.0.0
date-released: 2024-01-10
url: "https://github.com/johngarg/op-closure-estimates"

GitHub Events