This readme is for the data submitted in conjunction with the paper:
Title: "Enhanced optical cross section via collective coupling of atomic dipoles in a 2D array"
Authors: Robert J. Bettles, Simon A. Gardiner, and Charles S. Adams

The data was produced using Matlab and the figures produced in Inkscape. The data in the files described below is the raw data used in the paper figures; the data included is sufficient to reproduce all figures in the paper. 




========================================
Explanations of structure for data files
========================================


-------------------------
Files included in /Fig1/:
-------------------------

==================
(1.1) Fig1_N2D.csv
==================
2D number density of the uniform random distribution in the main plot of Fig. 1, in units of 1/lambda_0^2 (lambda_0 is the transition wavelength).

====================
(1.2) Fig1_rdips.csv
====================
Dipole positions in the N=100 random atomic monolayer in Fig 1, in units of lambda_0. Each set of three columns represents the (x,y,z) position components for a different number density. Each row is a different dipole in the lattice.

================
(1.3) Fig1_T.csv
================
Transmission in the main panel of Fig. 1. Each row is a different number density, corresponding to Fig1_N2D.csv. The columns correspond to multiple random configurations. For each cell in Fig1_T.csv, a different random configuration of positions was realized -- blue solid line.

=================
(1.4) Fig1_BL.csv
=================
Same as in (1.3) except assuming the dipoles are non-interacting (2D Beer-Lambert model) -- black dotted line.

=========================
(1.5) Fig1inset_Emag.csv
=========================
Electric field magnitude plotted in the inset of Fig. 1. Row i, collumn j corresponds to the ith x-position and jth z-position.

==========================
(1.6) Fig1inset_rdips.csv
==========================
Dipole position vectors of the inset. The columns are the (x,y,z) components and the rows are the different dipole positions. In units of lambda_0.

======================
(1.7) Fig1inset_x.csv
======================
The x values of the inset, in units of lambda_0.

======================
(1.8) Fig1inset_z.csv
======================
The z values of the inset, in units of lambda_0.




-------------------------
Files included in /Fig2/:
-------------------------

==================
(2.1) Fig2_N2D.csv
==================
2D number density of the triangular lattice in the main plot of Fig. 1, in units of 1/lambda_0^2.

==================
(2.2) Fig2_sep.csv
==================
Nearest neighbor spacing of the triangular lattice in the main plot of Fig. 1, in units of lambda_0.

====================
(2.3) Fig2_rdips.csv
====================
Dipole positions in the N=102 triangular lattice in Fig 2, in units of lambda_0. Each set of three columns represents the (x,y,z) position components for a different nearest neighbor spacing. Each row is a different dipole in the lattice.

================
(2.4) Fig2_T.csv
================
Transmission as a function of nearest neighbor spacing through the interacting triangular lattice, main panel of Fig. 2 -- blue solid line.

=================
(2.5) Fig2_BL.csv
=================
Same as in (2.4) except assuming the dipoles are non-interacting (2D Beer-Lambert model) -- black dotted line.

=========================
(2.6) Fig2inset_Emag.csv
=========================
Electric field magnitude plotted in the inset of Fig. 2. Row i, collumn j corresponds to the ith x-position and jth z-position.

==========================
(2.7) Fig2inset_rdips.csv
==========================
Dipole position vectors of the inset. The columns are the (x,y,z) components and the rows are the different dipoles positions. In units of lambda_0, the transition wavelength.

======================
(2.8) Fig2inset_x.csv
======================
The x values of the inset, in units of lambda_0.

======================
(2.9) Fig2inset_z.csv
======================
The z values of the inset, in units of lambda_0.




-------------------------
Files included in /Fig3/:
-------------------------

=================
(3.1) Fig3a_T.csv
(3.2) Fig3b_T.csv
(3.3) Fig3c_T.csv
=================
Transmission as a function of detuning in Fig. 3(a-c) -- blue solid lines.

====================
(3.4) Fig3a_T_BL.csv
(3.5) Fig3b_T_BL.csv
(3.6) Fig3c_T_BL.csv
====================
Transmission through a non-interacting array (2D Beer-Lambert model) as a function of detuning in Fig. 3(a-c) -- black dotted lines.

=====================
(3.7) Fig3a_rdips.csv
(3.8) Fig3b_rdips.csv
(3.9) Fig3c_rdips.csv
=====================
Dipole positions in Fig. 3(a-c). The three columns are the (x,y,z) components of the dipole position vectors. The rows are the different dipoles. 

========================
(3.10) Fig3_Detuning.csv
========================
Detuning of the driving field with respect to the bare dipole transition frequency in the main plot of Fig. 3. Normalized by the line halfwidth, gamma_0.

===========================
(3.11) Fig3inset_0.67_T.csv
(3.12) Fig3inset_0.92_T.csv
===========================
Transmission as a function of detuning in the inset of Fig. 3. 0.67 (purple solid line) and 0.92 (red solid line) refer to the nearest neighbor lattice spacing in units of lambda_0. 

===============================
(3.13) Fig3inset_0.67_rdips.csv
(3.14) Fig3inset_0.92_rdips.csv
===============================
Dipole positions in Fig. 3 inset. The three columns are the (x,y,z) components of the dipole position vectors. The rows are the different dipoles. 0.67 (purple solid line) and 0.92 (red solid line) refer to the nearest neighbor lattice spacing in units of lambda_0.

=============================
(3.15) Fig3inset_Detuning.csv
=============================
Detuning of the driving field with respect to the bare dipole transition frequency in Fig.3 inset. Normalized by the line halfwidth, gamma_0.



-------------------------
Files included in /Fig4/:
-------------------------

=======================
(4.1) Fig4_sep.csv
(4.2) Fig4inset_sep.csv
=======================
Lattice spacing a as a function of lambda_0, in the main plot (4.1) and inset (4.2). 

==================================
(4.3) Fig4a_Ndef=0_V0=inf_T.csv
(4.4) Fig4a_Ndef=0_V0=5000ER_T.csv
(4.5) Fig4a_Ndef=0_V0=500ER_T.csv
(4.6) Fig4a_Ndef=0_V0=50ER_T.csv
(4.7) Fig4a_Ndef=0_V0=5ER_T.csv
==================================
Transmission in Fig.4a as a function of lattice spacing through a square lattice with N=100 lattice sites, each occupied by a single dipole. (4.3) is for infinite trap depth. (4.4-4.7) are for trap depths of 5000ER, 500ER, 50ER and 5ER respectively, where ER is the recoil energy. Each column is a separate run, each row is a different nearest neighbour lattice spacing. For each cell in the .csv files (not just for each column), a different random realization of the dipole positions is taken, as described in the main text and Supplemental Material. 

==================================
(4.8) Fig4b_Ndef=0_V0=inf_T.csv
(4.9) Fig4b_Ndef=10_V0=inf_T.csv
(4.10) Fig4b_Ndef=20_V0=inf_T.csv
(4.11) Fig4b_Ndef=30_V0=inf_T.csv
(4.12) Fig4b_Ndef=40_V0=inf_T.csv
(4.13) Fig4b_Ndef=50_V0=inf_T.csv
==================================
Transmission in Fig.4b as a function of lattice spacing through a square lattice with N=100 lattice sites. As in (4.3-4.7), the files (4.8-4.13) have columns which are different realizations and rows which are different lattice spacings (except for (4.8) which has just one realization). For each realization, Ndef atoms are removed at random from the lattice. (4.8-4.13) correspond to 0, 10, 20, 30, 40 and 50 defects respectively. Each plot has infinite trap depth. 

=======================================
(4.14) Fig4a_Ndef=0_V0=inf_rdips.csv
(4.15) Fig4a_Ndef=0_V0=5000ER_rdips.csv
(4.16) Fig4a_Ndef=0_V0=500ER_rdips.csv
(4.17) Fig4a_Ndef=0_V0=50ER_rdips.csv
(4.18) Fig4a_Ndef=0_V0=5ER_rdips.csv
=======================================
Dipole positions for the plots in (4.3-4.7). In (4.4-4.7), we average over 200 different random configurations. We however only include the first as an example. Each set of three columns corresponds to the (x,y,z) vector components for a different nearest neighbor lattice spacing. Each row corresponds to a different atom in the lattice.

=====================================
(4.19) Fig4b_Ndef=0_V0=inf_rdips.csv
(4.20) Fig4b_Ndef=10_V0=inf_rdips.csv
(4.21) Fig4b_Ndef=20_V0=inf_rdips.csv
(4.22) Fig4b_Ndef=30_V0=inf_rdips.csv
(4.23) Fig4b_Ndef=40_V0=inf_rdips.csv
(4.24) Fig4b_Ndef=50_V0=inf_rdips.csv
=====================================
Same as (4.14-4.18) but for the lines in Fig.4b, this time averaged over 400 different random configurations.

======================================
(4.25) Fig4inset_Ndef=0_V0=inf_T.csv
(4.26) Fig4inset_Ndef=10_V0=50ER_T.csv
======================================
Transmission of the two lines in Fig.4b inset. 

==========================================
(4.27) Fig4inset_Ndef=0_V0=inf_rdips.csv
(4.28) Fig4inset_Ndef=10_V0=50ER_rdips.csv
==========================================
Dipole positions of the two lines in Fig.4b inset, as in (4.14-24).




===================================================================
For queries, contact Robert Bettles at rjbettles89@gmail.com

