Changeset 288:63ccb2315004 for source/class/pf

Timestamp:

Sep 14, 2011, 12:59:42 PM (13 years ago)

Author:

Evgeny Stambulchik <Evgeny.Stambulchik@…>

Branch:

default

Phase:

public

Message:

Added oscillator strength (for absorption) using Kramers quasi-classical
formula + Menzel&Pekeris(1935) for the Gaunt factor.

Location:

source/class/pf

Files:

• Plasma.js (modified) (1 diff)
• base/Bohr.js (modified) (1 diff)

source/class/pf/Plasma.js

@@ -55 +55 @@
             return Math.pow(Z_core, 4)*
               (pf.base.Bohr.LymanA(n_u) + pf.base.Bohr.LymanA(n_l));
+        },
+        
+        getOscillatorStrength: function()
+        {
+            var n_u = this.n_u;
+            var n_l = this.n_l;
+            
+            return pf.base.Bohr.getOscillatorStrength(n_u, n_l);
         },
         

source/class/pf/base/Bohr.js

@@ -80 +80 @@
             }
         },
+
+        /* Quasi-classical absorption oscillator strength */
+        getKramersF: function(n_u, n_l)
+        {
+            return 32/(3*Math.PI*Math.sqrt(3))*
+                n_u*n_u*n_u*n_l/Math.pow(n_u*n_u - n_l*n_l, 3);
+        },
+
+        getFGauntFactor : function(n_u, n_l)
+        {
+            var g = 1;
+            
+            if (n_u == 2 && n_l == 1) {
+                g = 0.717;
+            } else
+            if (n_u == 3 && n_l == 1) {
+                g = 0.765;
+            } else
+            if (n_u == 3 && n_l == 2) {
+                g = 0.757;
+            } else
+            if (n_u - n_l == 1) {
+                g = 0.782;
+            } else {
+                var r2 = (n_l/n_u)*(n_l/n_u);
+                var t = Math.pow((1 - r2)*n_l, 2/3);
+                g = 1 - 0.1728*(1 + r2)/t;
+                     // - 0.0496*(1 + 484/15*r2 + r2*r2)/(t*t);
+            }
+            
+            return g;
+        },
+
+        getOscillatorStrength : function(n_u, n_l)
+        {
+            var g = this.getFGauntFactor(n_u, n_l);
+            
+            return g*this.getKramersF(n_u, n_l);
+        },
         
         classicalIonizationField : function(Z_core, n)