h = 6.626E-34; % Js

hbar = h/(2*pi);

w0 = 2*pi*1.06E9; % s^-1

e = 1.6E-19; % C

a0 = 5.3E-11; % m



E = 1/(2*sqrt(3))*hbar*w0/(e*a0);