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.
Load up the definitions in the core files:
Get["main.m"];
Get["operators.m"];
Get["rules.m"];
Get["decay.m"];
Get["pmatch.m"];
Get["symmetries.m"];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
MatchingData["11", {1, 1, 1, 1}]The output is in the mass basis.
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 |