Figure1Data.csv contains x and y values as well as y half-width ranges to reproduce
figure 1.

Figure2Data.csv contains mode eigenvalues in each column.  To reproduce figure
2, plot histograms of the (negative) of the real and imaginary components of
each column.

Figure3Data.csv contains x and y data to reproduce the line in figure 3 - all
three plots are contained in this single line, the plots differ only in the axis limits.

Figure4aData.csv contains the data to reproduce the lines in Fig 4a, and
Figure4b_and_cData.csv that for lines in figures 4b and 4c (which have the same x
values).

For each of the various plots of figure 5, there are separate files for the full quantum
field-theoretical results and the stochastic results, giving x and y
coordinates.

For each plot of Figure 6 there is a data file with the dimensionless
polarization density as a function of x.

For each of the various plots of figure 7, there are separate files for the full quantum
field-theoretical results and the stochastic results, giving x and y
coordinates.

For figure 8, each subplot has a data file for the stochastic results, and a
separate file containing the full quantum field-theoretical results and the
single atom effective coupling strength approximation.

For figure 9, date is provided in two ways.  Data for subplots a,c,e, and g are
provided with the z data in the comma delimited file of the corresponding name. 
In these files, rows correspond to different x values, and columns to different
x prime values.  The values for x and x prime are tabulated in Figure9xdata.csv
(common to all 4 subplots).  Note that the plots show only the absolute value, while
the data gives the complex value of the appropriate correlation function.  For
subplots b,d,f and h, the plots are formed of the product of two one-body
correlation functions - those one body correlation functions are tabulated (with
the relevant x variable) in the corresponding csv files.  Two dimensional plots
can be formed from this data by plotting the matrix formed by the outer product
of the two vectors formed by the one-body correlation functions.

The eigenvalues used to generate figure 10 are in Figure10data.csv.  Each column
represents an eigenvalue, with the last column containing the superradiant
eigenvalues.  To obtain the figures in the paper, plot histograms of (negative
of) the real and imaginary components of the super and sub-radiant columns. 
Note that this data is in units where $g_0^2 = 9.4799e+04$, and $\bar{kappa} =
$3.4211e+02 - 3.4211e+04i$.

The data the subplots in figure 11 are found in the correspondingly named files,
showing the Mott insulating (MI) and superfluid (SF) results as a function of
the dimensionless cavity detuning. For subplots c and f, there are separate
files for the Mott insulating (MI) and superfluid (SF) results, due to different
x ranges.