Changeset 289:a3d60ac74412 for source/class

Timestamp:

Sep 14, 2011, 6:21:42 PM (13 years ago)

Author:

Evgeny Stambulchik <Evgeny.Stambulchik@…>

Branch:

default

Phase:

public

Message:

Einstein coefficient.

Location:

source/class/pf

Files:

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

source/class/pf/Plasma.js

@@ -62 +62 @@
             var n_l = this.n_l;
             
-            return pf.base.Bohr.getOscillatorStrength(n_u, n_l);
+            return pf.base.Bohr.oscillatorStrength(n_u, n_l);
+        },
+        
+        getEinsteinA: function()
+        {
+            var n_u = this.n_u;
+            var n_l = this.n_l;
+            var Z_core = this.r.getQ() + 1;
+            var m_core = this.r.getM();
+            
+            return pf.base.Bohr.EinsteinA(Z_core, m_core, n_u, n_l);
         },
         

source/class/pf/base/Bohr.js

@@ -8 +8 @@
         // alpha^2
         alpha2 : 5.325135e-05,
+        
+        // alpha^3
+        alpha3 : 3.885939e-07,
         
         // Speed of light (1/alpha)
@@ -82 +85 @@
 
         /* Quasi-classical absorption oscillator strength */
-        getKramersF: function(n_u, n_l)
+        KramersF: function(n_u, n_l)
         {
             return 32/(3*Math.PI*Math.sqrt(3))*
@@ -88 +91 @@
         },
 
-        getFGauntFactor : function(n_u, n_l)
+        fGauntFactor : function(n_u, n_l)
         {
             var g = 1;
@@ -113 +116 @@
         },
 
-        getOscillatorStrength : function(n_u, n_l)
+        oscillatorStrength : function(n_u, n_l)
         {
-            var g = this.getFGauntFactor(n_u, n_l);
+            var g = this.fGauntFactor(n_u, n_l);
             
-            return g*this.getKramersF(n_u, n_l);
+            return g*this.KramersF(n_u, n_l);
+        },
+        
+        EinsteinA : function(Z_core, m_core, n_u, n_l)
+        {
+            var f  = this.oscillatorStrength(n_u, n_l);
+            var dE = this.transitionEnergy(Z_core, m_core, n_u, n_l);
+            
+            return 2*this.alpha3*dE*dE*f*(n_l/n_u)*(n_l/n_u);
         },