This repository contains the data used to make the figures in
"Controlling lasing around exceptional points in coupled nanolasers"
(A. Fischer et al.; https://doi.org/10.1038/s44310-024-00006-9).

The full repository should be available under https://doi.org/10.5281/zenodo.10683232

The files are named according to which figure they correspond to, and are as follows:

# Figure 1
* fig1c_meta.csv - The pump power ratio, and the pump powers and fluences on each disk for 
* fig1c_counts.csv - The counts to plot the three spectra
* fig1c_wavelength.cvs - Corresponding wavelength (nm) to the spectra

Figure 1d was obtained by solving the coupled mode theory matrix for its eigenvalues and plotting the imaginary and real part

H(\tau) = \begin{pmatrix} \omega_0 + \delta_1 + iP_1 & \kappa \\ 
    \kappa & \omega_0+\delta_2+iP_2(\tau) \end{pmatrix}

* fig1d_meta.csv - Gives the input data for the CMT 
* fig1d_nodetuning1.csv - Real and Imaginary part of mode 1 for no detuning (dashed line in figure)  
* fig1d_nodetuning2.csv - Real and Imaginary part of mode 2 for no detuning (dashed line in figure)  
* fig1d_detuning1.csv - Real and Imaginary part of mode 1 with detuning and corresponding pump ration P1/P2 (coloured scatter plot in figure)  
* fig1d_detuning2.csv - Real and Imaginary part of mode 2 with detuning and corresponding pump ration P1/P2 (coloured scatter plot in figure)  



# Figure 2
For the theory plots in Fig. 2 (an and d) non-linear CMT described in the SI of the paper was used.

* fig2a_meta.csv - Input values for CMT
* fig2a_counts_mode1 - 2d array of mode intensity for mode 1
* fig2a_counts_mode2 - 2d array of mode intensity for mode 2
* fig2a_P1grid.csv - 2d array of P1 values for colormap of modes
* fig2a_P2grid.csv - 2d array of P2 values for colormap of modes
* fig2a_threshold_mode1 - P1 and P2 array for threshold line of mode1
* fig2a_threshold_mode2 - P1 and P2 array for threshold line of mode2

* fig2b_counts.csv - 2d array of lasing intensity (counts)
* fig2b_F1.csv - 2d array of fluence on disk 1 in uJ cm^(-2) pulse^(-1)
* fig2b_F2.csv - 2d array of fluence on disk 2 in uJ cm^(-2) pulse^(-1)
* fig2b_thresholds.csv - Threshold of mode 1 and 2 to plot the threshold lines and also to distinguish the two modes 

* fig2c_meta.csv - Input values for CMT
* fig2c_counts_mode1.csv - 2d array of mode intensity for mode 1
* fig2c_counts_mode2.csv - 2d array of mode intensity for mode 2
* fig2c_P1P2_colormap.csv - P1 and P2 array for axes in colormap of mode 1 and 2
* fig2c_threshold_mode1.csv - P1 and P2 array to plot threshold line of mode 1 (simple CMT -> dashed line)
* fig2c_threshold_mode2.csv - P1 and P2 array to plot threshold line of mode 2 (simple CMT -> dashed line)
* fig2c_threshold_stability_mode1.csv - P1 and P2 array to plot stability threshold line of mode 1 (non-linear CMT -> solid line)
* fig2c_threshold_stability_mode2.csv - P1 and P2 array to plot stability threshold line of mode 2 (non-linear CMT -> solid line)

* fig2d_F1.csv - 2d array of fluence on disk 1 in uJ cm^(-2) pulse^(-1)
* fig2d_F2.csv - 2d array of fluence on disk 2 in uJ cm^(-2) pulse^(-1)
* fig2d_mode1.csv - 2d array of intensity (counts) of mode 1 (blue mode in figure)
* fig2d_mode2.csv - 2d array of intensity (counts) of mode 2 (blue mode in figure)
* fig2d_overlap1.csv - 2d array of intensity (counts) of overlapping/ middle section where both modes are present (brown).
* fig2d_overlap2.csv - 2d array of intensity (counts) of mode 1 (blue mode in figure) 
* fig2d_thresholds_mode1 - Data for threshold points of mode 1 in uJ cm^(-2) pulse^(-1) (1. column = P2 values, 2. & 3. = P1 values of start & end lasing gap)
* fig2d_thresholds_mode2 - Data for threshold points of mode 2 in uJ cm^(-2) pulse^(-1) (1. column = P1 values, 2. & 3. = P2 values of start & end lasing gap)



# Figure 3
* fig3a_counts.csv - 2d array of lasing intensity (counts)
* fig3a_F1.csv - Fluence ơn disk 1 values (uJ cm^(-2) pulse^(-1)) used for y-axis
* fig3a_wavelength.csv - Wavelength values used for x-axis

* fig3b.csv - First three columns give the Power on disk 2, lasing gap start and lasing gap end in uJ cm^(-2) pulse^(-1). Last three columns give the same quantities in W.

* fig3c.csv - Three columns giving Lasing gap size (uJ cm^(-2) pulse^(-1)), Wavelength shift (nm), and Error (nm)



# Figure 4
* fig4a_counts.csv - 2d array of emission counts of the laser used for the upper plot
* fig4a_wavelength.csv - Wavelength values used for x-axis
* fig4a_P1.csv - Powers on disk 1 (uJ cm^(-2) pulse^(-1)) corresponding to the data in fig4a_counts.csv 
* fig4a_spectra.csv - Two spectra (rows 18 and 26 (starting from 0) in fig4a_counts.csv) plotted in the lower part of the figure
* fig4a_spectraPowers.csv - The two powers on disk 1 (uJ cm^(-2) pulse^(-1)) corresponding to the two spectra

* fig4b_total.csv - Peak wavelengths for each measurement number
* fig4b_lower.csv - Peak wavelengths and measurement number for lower P1 power (longer wavelength) 
* fig4b_higher.csv - Peak wavelengths and measurement number for higher P1 power (shorter wavelength) 


If you would like the raw data and how to process them to yield the above data do not hesitate to contact me under: anna.fischer@ibm.com



