From 89ea3cafc79fc3dfb490249c52f145dd9d329045 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Beno=C3=AEt=20Sierro?= <benoit.sierro@bluewin.ch>
Date: Wed, 22 Feb 2023 11:42:06 +0100
Subject: [PATCH] fixed loss factor

---
 src/scgenerator/physics/fiber.py | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/src/scgenerator/physics/fiber.py b/src/scgenerator/physics/fiber.py
index 68efb12..c77eaf8 100644
--- a/src/scgenerator/physics/fiber.py
+++ b/src/scgenerator/physics/fiber.py
@@ -903,7 +903,8 @@ def capillary_loss(wl: np.ndarray, he_mode: tuple[int, int], core_radius: float)
         loss in 1/m
     """
     chi_silica = abs(mat.Sellmeier.load("silica").chi(wl))
-    nu_n = 0.5 * (chi_silica + 2) / np.sqrt(chi_silica)
+    # the real loss alpha is 2*Im(n_eff), which differs from the notation of the paper
+    nu_n = (chi_silica + 2) / np.sqrt(chi_silica)
     return nu_n * (u_nm(*he_mode) * wl / pipi) ** 2 * core_radius**-3