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.