This folder contains all the code needed to perform the optical binding calculations using an AOM.
There are two types of Julia codes to perform the simulation:

- calculate_force_x.jl / calculate_force_y.jl: Code to compute the optical binding force on the two particles for several Gaussian beam positions, as described in the SI. 'x' is for parallel and 'y' is for perpendicular polarization.
     The output is a file named a_260_f1_xdata.dat/a_260_f1_ydata.dat (or the same with f2 for particle 2) that contains the optical binding force. Lines are distances and columns are beam positions.
- get_potential_x.jl / get_potential_y.jl: Code to do the weighted average and perform the numerical integration. 
     The output is a file where the first column is the distance and the second column is the optical binding potential.

This folder also contains the beam intensity profiles used to perform the integration. These are file starting with 'um_intensity'
