%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Data statement  
%
% for "Using entanglement to discern phases in the disordered one-
% dimensional Bose-Hubbard model" 
% by Andrew M. Goldsborough and Rudolf A. R\"{o}mer
%
% 09/08/2015
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Contents:

README.txt
one_bdbcorr.zip
one_ee.zip
one_eespec.z01
one_eespec.z02
one_eespec.zip
one_energy_nminus1.zip
one_energy_nplus1.zip
one_nexpect.zip
one_output.zip
half_bdbcorr.zip
half_ee.zip
half_eespec.z01
half_eespec.z02
half_eespec.zip
half_nexpect.zip
half_output_apbc.zip
half_output_pbc.zip
two_bdbcorr.zip
two_ee.zip
two_eespec.z01
two_eespec.z02
two_eespec.z03
two_eespec.z04
two_eespec.z05
two_eespec.z06
two_eespec.z07
two_eespec.zip
two_energy.zip
two_energy_nminus1.zip
two_energy_nplus1.zip
two_nexpect.zip

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

The eespec folders are too large to be uploaded into the WRAP system and are 
therefore split into sections. To extract type the following into the linux 
terminal (e.g. for filling one):

cat one_eespec.z01 one_eespec.z02 one_eespec.zip > one_eespec_full.zip
zip -FF eespec_full.zip --out eespec_fixed.zip
unzip eespec_fixed.zip

For the fixing stage (zip -FF) select skip when it asks for additional files.
The zipping was performed in SUSE Linux Enterprise Desktop 11. The split
folders should have sizes 4.7Gb, 5.2Gb and 13Gb for fillings one, half and 
two respectively.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

The folder names are labelled by the filling fraction (one, half and two) 
then the observable type:
bdbcorr - b, b^{\dagger} correlation function
ee - entanglement entropy
eespec - entanglement spectrum, the singular values for each SVD
energy - graound state energy
energy_nminus1 - energy for filling -1
energy_nplus1 - energy for filling +1
nexpect - number operator expectation value
output - terminal output containing energy
output_apbc - output for anti-periodic boundary conditions
output_pbc - output for periodic boundary conditions

The ground state energy was originally scraped from the output file, but 
newer versions of the DMRG program outputted it as a separate file for 
convenience.

The energy_nminus1 and energy_nplus1 are for use when calculating the energy
gap in the filling one and two phase diagrams.

The output_apbc and output_pbc are for calculating the superfluid stiffness
for half filling in figure 4a.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Within each folder each file is named according to the variables used to 
create it. Usually this takes the form:

L_tstr_tdis_Ustr_Udis_mustr_mudis_seed_chi_type_5.txt

where:
L = length
tstr = t strength, which is set to 1 throughout
tdis = diorder on t, unused,set to zero
Ustr = U strength
Udis = disorder on U, unused, set to zero
mustr = strength of chemical potential, set to zero
mudis = disorder on mu
seed = seed for random number generator
chi = number of states kept in DMRG
type = data type (bdbcorr, ee, etc.)
5 = number of bosons per site in DMRG, set to 5

for the energy_nminus1 and energy n_plus1 there is an extra variable for 
number of bosons (npart):

L_npart_tstr_tdis_Ustr_Udis_mustr_mudis_seed_chi_energy_5.txt

the output file names are of the form:

disorderedBH5_100_100_0_200_0_0_100_1-100_200.pbs.o2642477.29
disorderedBH5_L_tstr_tdis_Ustr_Udis_mustr_mudis_seed_chi.pbs.oxxxxxxx

where the variables are as before but are multiplied by 100 so they can be 
generated in a bash script without floating point numbers. The "xxxxxxx"
is a job number from the queueing system. Later executions batch executed 
and are of the form:

disorderedBH5_L_tstr_tdis_Ustr_Udis_mustr_mudis_1-100_chi.pbs.oxxxxxxx.seed

Many of apbc/pbc runs are called "simple" and only produced energy outputs
to reduce execution time as the periodic boundaries require higher chi and
more sweeps.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Within each file the data in each column is as follows:
bdbcorr - site 1, site 2, correlation
ee - split position, entanglement entropy
eespec - split position, singular values from SVD at that point
energy - energy
energy_nminus1 - energy
energy_nplus1 - energy
nexpect - site, n expectation value
output - terminal output
output_apbc - terminal output 
output_pbc - output

The seed used in filling two has the formula:

1234*i + 5678

to make sure that the choice of seed is not related to the similar shape
of the phase diagram.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%