Changeset 367:e1a3480589f2

Timestamp:

Apr 14, 2016, 2:41:03 PM (8 years ago)

Author:

Evgeny Stambulchik <Evgeny.Stambulchik@…>

Branch:

default

Phase:

public

Message:

Implemented rotational Stark broadening for Delta-n=1 transitions.

File:

• source/class/pf/Plasma.js (modified) (4 diffs)

source/class/pf/Plasma.js

@@ -83 +83 @@
         getQsHwhm: function(s)
         {
-            var n_u = this.r.getN_u();
-            var n_l = this.r.getN_l();
-
             var ppi = this.getPpiFactor(s);
             var rpi = this.getRpiFactor(s);
@@ -91 +88 @@
             // 1.4385 is S_1 HWHM
             var hwhm = 1.4385*this.getDelta0(s)*ppi*rpi;
-            
-            // rude correction for alpha lines
-            if (n_u - n_l == 1) {
-                hwhm /= 2;
-            }
             
             return hwhm;
@@ -134 +126 @@
         getStarkHwhm: function(s)
         {
+            var n_u = this.r.getN_u();
+            var n_l = this.r.getN_l();
+
             var w_qs = this.getQsHwhm(s);
-            var f  = this.getStarkQuasistaticity(s);
-            
-            return f*w_qs;
+            
+            // "rotational" model for delta-n = 1 and QC for anything else
+            if (n_u - n_l == 1) {
+                var R = this.getStarkRatio(s);
+                if (R != 0) {
+                    return w_qs/(2.0/R + R/1.2);
+                } else {
+                    return 0;
+                }
+            } else {
+                var f = this.getStarkQuasistaticity(s);
+
+                return f*w_qs;
+            }
+            
         },
         
         getTotalStarkHwhm: function()
         {
+            var n_u = this.r.getN_u();
+            var n_l = this.r.getN_l();
+
             var w = 0;
             
@@ -149 +159 @@
                 var w_s = this.getStarkHwhm(s);
 
-                w += Math.pow(w_s, 3/2);
-            }
-            
-            return Math.pow(w, 2/3);
+                if (n_u - n_l == 1) {
+                    w += w_s;
+                } else {
+                    w += Math.pow(w_s, 3/2);
+                }
+            }
+            
+            if (n_u - n_l == 1) {
+                return w;
+            } else {
+                return Math.pow(w, 2/3);
+            }
         },