The code Figure14.py needs the following input files:
phasediagram_gsoverlap_OBC_Npart=2_hermTrue_D0=1.6000_B0=12.0000_rtrunc=-1.0000.txt
phasediagram_gsoverlap_OBC_Npart=2_hermTrue_D0=2.1000_B0=12.0000_rtrunc=-1.0000.txt
phasediagram_gsoverlap_OBC_Npart=3_hermTrue_D0=1.6000_B0=12.0000_rtrunc=-1.0000.txt
phasediagram_gsoverlap_OBC_Npart=3_hermTrue_D0=2.1000_B0=12.0000_rtrunc=-1.0000.txt
containing the overlaps between the model Laughlin ground state and the true ground state of respective angular momentum sector.

They can be created by running:
python phasediagram_overlap.py

The Npart, D0, as well as the ranges of d0 and Vharm0 can be changed by editing the code (see comments inside)

Note that the code produces also files phasediagram_gsnorm..., which contain the projection of the true ground state on the subspace of configurations which make up the Laughlin ground state. This provides some information on what is the reason that the overlap is not perfect: is it because of admixture of other configurations, or the change of weights within the same configurations? These files are not used by the plotting script.

Also, note that the calculation of the phase diagram for 3 particles at D0=2.1 takes a long time (more than 1 day) on a laptop.

With the input files present, the figure is created by running:
python Figure14.py
