La fuerza en la dirección $j$ ejercida sobre los átomos debido al par de haces en esa dirección está dada por
$$
F_j = F_j^+ + F_j^-,
$$
dónde
$$
F_j^\pm = \pm \hslash k \dfrac{\gamma_{22}}{2} \times
$$
$$
\times \dfrac{ I/I_0 }{ 1 + I/I_0 + 4(\delta \mp kv_j \mp \beta_j r_j)^2/\gamma_{22}^2 }.
$$
$\gamma_{22}$ es la tasa de desintegración espontánea de los átomos del estado |2⟩ al estado |1⟩,
$k$ es el número de onda, $I$ es la intensidad del rayo láser,
$I_0$ es la intensidad de saturación de la transición |1⟩ - |2⟩,
$\delta$ es la desafinación del láser con respecto a la resonancia de transición,
$v$ es la velocidad del átomo y $\beta_j$ está definida por
$$
\beta_j = \dfrac{g\mu_B}{\hslash} \dfrac{\partial B_j}{\partial j}.
$$
$g$ es el factor giromagnético, $\mu_B$ es el magnetón de Bohr y $B_j$ el componente $j$
del campo magnético generado por el par de bobinas de Helmholtz.
La dinámica de los átomos se basa en la segunda ley de Newton y una variable aleatoria $\epsilon = 0,5$
Se agregó m/s para simular las fluctuaciones aleatorias en las velocidades atómicas debido a emisiones espontáneas:
$$
v(t) \rightarrow v(t) + \epsilon.
$$
Lea mas: