Dynamic polarizability: make sure states with pqn < ground state are …

…not counted twice (see #72)
nikolasibalic · Jun 15, 2021 · 202447b · 202447b
1 parent f40bc02
commit 202447b
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/arc/calculations_atom_single.py b/arc/calculations_atom_single.py
@@ -2388,7 +2388,8 @@ def defineBasis(self, nMin, nMax):
         self.basis = []
         self.lifetimes = []
 
-        for n1 in np.arange(self.nMin, self.nMax + 1):
+        for n1 in np.arange(max(self.nMin, self.atom.groundStateN),
+                            self.nMax + 1):
             lmin = self.l - 1
             if (lmin < - 0.1):
                 lmin = self.l + 1