README: Moiré-Induced Optical Nonlinearities: Single and Multiphoton Resonances 

Folder containing the data for the manuscript Moiré-Induced Optical Nonlinearities: Single and Multiphoton Resonances by A. Camacho-Guardian and N. R. Cooper.

METHODS 
As described in following, the data was obtained either by exact diagonalization of the Lindblad operator or the solution of the extended Gross-Pitaevskii as defined in the main text.

FIGURES
Figure 2
a) Phase diagram obtained by exact diagonalization of the Lindblad operator. Regimes of high and low density are separated by a bistability that can be accessed via hysteresis. 
      Fig2a.dat
(b) and (c) Exciton number obtained by exact diagonalization of the Lindblad operator. Fig. 2(a) corresponds to the low branch where the pump strength is varied from below. Fig. 2(b) illustrates the hysteresis as the pump strength is tuned from an initial large value. 
    Fig2b.dat
    Fig2c.dat 
(d) Purity for the low and hysteresis branch for a cross-section of the set of parameters in Fig. 2(a).
   Fig2d_1.dat
   Fig2d_2.dat
(e) Phase diagram including saturation effects obtained by exact diagonalization of the Lindblad operator.
    Fig2e.dat
Figure 3
Exciton amplitude obtained for the extended Gross-Pitaevskii for U=0
   Fig3_PRL_Non_Interacting.dat
 Three data sets for the numerical solution of extended Gross-Pitaevskii (eGP) corresponding to the multi-stability of the eGP in the presence of interactions. 
  Fig3_PRL_GP1.dat
   Fig3_PRL_GP2.dat
  Fig3_PRL_GP3.dat
Full quantum treatments of the excitons
  Fig3_PRL_Non_Quantum. dat

Figure 4
Full quantum treatment of the excitons in the presence of exciton gain obtained by exact diagonalization of the Lindblad operator including the term in Eq. 7 described in the main text.
Principal figure following the colour code of the figure.
  Fig4_Black.dat
  Fig4_Green.dat
  Fig4_Red.dat

Inset
  Fig4_Inset_Black.dat
  Fig4_Inset_Blue.dat
  Fig4_Inset_Red.dat