Solves the Richardson-Gaudin equations for integrable SO(5) Hamiltonians.
Along with standard packages (numpy, scipy, pandas,
matplotlib) included in the Anaconda distribution,
this code requires the quspin package for
exact-diagonalization.
To install quspin (if you have Anaconda installed),
run
conda install -c weinbe58 quspin
in your terminal.
You should also have the multiprocessing and concurrent
packages installed.
Some of these scripts require a file named context.json
containing a filepath for saving CSVs of results. Format is as follows:
{"results_filepath": "/path/to/your/folder"}
To solve the rational SO(5) model, run
python solve_rg_eqs.py
in your terminal and input the requested parameters. Recommended: L = 4, Ne = Nw = 2, G = 1, and dg = 0.01.
solve_rg_eqs.py solves the Richardson-Gaudin
equations and derive observables for the SO(5) RG models.
exact_diag.py performs exact diagonalization
on the same Hamiltonians using QuSpin.
spectral_fun.py computes spectral functions
using exact diagonalization.
utils.py contains small functions used in the other files.