From af5cefc8649e0ebce4e0a3b6ebc79edd3e94c505 Mon Sep 17 00:00:00 2001
From: aibishop <125414009+aibishop@users.noreply.github.com>
Date: Thu, 16 Feb 2023 14:18:14 +1100
Subject: [PATCH] Update optical_tweezers.py

---
 monashspa/PHS3000/optical_tweezers.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/monashspa/PHS3000/optical_tweezers.py b/monashspa/PHS3000/optical_tweezers.py
index ecc6ef2..dd5be24 100644
--- a/monashspa/PHS3000/optical_tweezers.py
+++ b/monashspa/PHS3000/optical_tweezers.py
@@ -66,7 +66,7 @@ def trap_k_theory(r, w, alpha, eccentricity, I):
     c = 299792458
 
     # calculate k
-    k = np.sqrt(alpha)/4*alpha/c*I*w*2*np.pi*eccentricity**3/xi**3*(
+    k = alpha/c*I*w*2*np.pi*eccentricity**3/xi**3*(
             np.sqrt(np.pi/2)*((xi*a/eccentricity)**2-1)*np.exp(-a**2/2)*scipy.special.erf(xi*a/(np.sqrt(2)*eccentricity)) 
             + (xi*a/eccentricity)*np.exp(-a**2/(2*eccentricity**2))
         )