|
|
---
|
|
|
title: RGE Running with Mathematica
|
|
|
permalink: /RGE_Running_with_Mathematica/
|
|
|
---
|
|
|
|
|
|
[Category:Calculations](/Category:Calculations "wikilink") [Category:Output](/Category:Output "wikilink") SARAH writes the RGEs into a file which can be used for a numerical evaluation of the RGEs in <span>Mathematica</span>. This file is stored in the same directory and called <span>RunRGEs.m</span> and provides a function called <span>RunRGEs</span>. The syntax of this function is
|
|
|
|
|
|
RunRGEs[values, start, finish, Options]
|
|
|
|
|
|
The input is:
|
|
|
|
|
|
- `values`: all non-zero values for parameters at the scale where the running starts
|
|
|
- `start`: logarithm of the scale where the running starts
|
|
|
- `finish`: logarithm of the scale where the running should end
|
|
|
- `Options`: optionally two-loop contributions can be turned off (`TwoLoop->False` )
|
|
|
|
|
|
### Example
|
|
|
|
|
|
One can load the file for the MSSM and perform an one-loop running of the gauge couplings from 1 TeV to 10<sup>16</sup> GeV via
|
|
|
|
|
|
<< "[\$SARAH-Directory]/Output/MSSM/RGEs/RunRGEs.m";
|
|
|
solution = RunRGEs[{g1->0.45, g2->0.63, g3->1.04},3,16, TwoLoop->False];
|
|
|
|
|
|
Here, we saved the output of <span>RunRGEs</span> as variable, which can be used as follows
|
|
|
|
|
|
{g1[16], g2[16], g3[16]} /. solution[[/1|1]];
|
|
|
Plot[g1[x]/. solution[[/1|1]], {x,3,16}];
|
|
|
|
|
|
The result is shown here: [gerahmt|zentriert](/Datei:RGEsMSSM.png "wikilink")
|
|
|
|
|
|
See also
|
|
|
--------
|
|
|
|
|
|
- [Generic RGE calculation](/Generic_RGE_calculation "wikilink") |