Estimate Jones Polynomials for Ising anyons on IBM Quantum devices
This is code to accompany "Estimating the Jones polynomial for Ising anyons on noisy quantum computers" by Chris N. Self, Sofyan Iblisdir, Gavin K. Brennen, and Konstantinos Meichanetzidis https://arxiv.org/abs/2210.11127
Combine this repo with the experimental data in zenodo https://doi.org/10.5281/zenodo.7194971 to explore other data fittings or plot the results.
Have any questions? Contact me here: christopher.self@quantinuum.com.
Requirements
The main requirements are:
pandas ~= 1.4
qiskit ~= 0.25
jupyterlab
luigiEverything will definitely work with these versions. Newer versions of pandas may not be able to load the data files. Newer versions of qiskit should be compatible, maybe with small fixes. luigi is only used by the script generate_fits.py.
How to use
The notebook example-knots-circuits.ipynb shows a walkthrough of how to generate the Hadamard test circuits and the stretched zero-noise extrapolation copies. These can be submitted to quantum backends, such as IBMQ, to generate new experimental data.
Once the experimental data is downloaded from zenodo and the 'analysis' folder is placed inside this repo the view... notebooks can be used to plot the results. The download contains the ZNE fits, but generate_fits.py can be re-run or extended to compute new zero-noise extrapolation fits. The script generate_fits.py uses the pipeline tool luigi, see the luigi docs for how to run it.
Files and folders
example-knots-circuits.ipynbwalkthrough for how to generate the Hadamard test circuitsgenerate_fits.pyscript to perform ZNE fits to the raw dataknots_circuits.pyfunctions to generate Hadamard test circuitsview-hadamard-test-results.ipynbplot the results for the complex partition functionview-jones-poly-estimates.ipynbplot the final results for the Jones polynomial