An astrophysics based Python code to simulate flaring activity of a magnetar using the Monte Carlo statistical technique. The MAGNETO code generates a sequence of flares from a single magnetar over its lifetime to determine how the time between two flares and flare energies evolve over the magnetar's age.
To install the MAGNETO code, run the following commands in a terminal:
> git clone https://github.com/Kushaalkumar-pothula/MAGNETO-Code.git
> cd MAGNETO-Code
Now you should have the MAGNETO code successfully installed on your computer. Now you can proceed to run the code.
Dependencies: Python 3, NumPy and Matplotlib. There are two ways to run MAGNETO: By following the default method or the custom method. The custom method gives more freedom over physical parameters.
- Open MAGNETO-Code directory in your terminal.
- Run the following command:
> bash run.sh
This will run the simulation with some default values of the magnetar's magntic field energy and power-law index.
In this method of running MAGNETO, you will have control over the physical aspects of the magnetar. You can specify the values of the respective parameters using options while running scripts. To get help with them, you can use the -h
flag while running the script:
> python energy_evol_modular.py -h
usage: energy_evol_modular.py [-h] E [E ...] a [a ...]
Simulate flaring activity of a magnetar using the Monte Carlo method.
positional arguments:
E Total magetic energy of the magnetar
a Power of time in energy-time relation, known as alpha
optional arguments:
-h, --help show this help message and exit
Similarly, you can run the interval evolution simulation:
> python interval_evol_modular.py
Kushaal Kumar Pothula (https://kushaalkumarpothula.wordpress.com/). Initially written in 2020, when I was 15.