cool_argo.cpp

Go to the documentation of this file.

/* This file is part of Cloudy and is copyright (C)1978-2013 by Gary J. Ferland and
 * others.  For conditions of distribution and use see copyright notice in license.txt */
/*CoolArgo compute argon cooling */
#include "cddefines.h"
#include "coolheavy.h"
#include "phycon.h"
#include "ligbar.h"
#include "taulines.h"
#include "dense.h"
#include "thermal.h"
#include "lines_service.h"
#include "atoms.h"
#include "cooling.h"
 
void CoolArgo(void)
{
        realnum a12, 
          a13, 
          a14, 
          a15, 
          a23, 
          a24, 
          a25, 
          a34, 
          a35, 
          a45, 
          cs12, 
          cs13, 
          cs14, 
          cs15, 
          cs23, 
          cs24, 
          cs25, 
          cs34, 
          cs35, 
          cs45, 
          pop2, 
          popn3;
        double cs, p[5];
        static double gAr4[5]={4.,4.,6.,2.,4.};
        static double exAr4[4]={21090.4,128.9,13636.2,177.1};
 
        DEBUG_ENTRY( "CoolArgo()" );
 
        /* Argon II 6.98 micron
         * trans prob from 
         * >>refer      Ar2     AS      Nussbaumer, H., & Storey, P. J. 1988, A&A, 200, L25
         * Y(ik) from
         * >>refer      Ar2     CS      Pelan, J., & Berrington, K. A. 1995, A&AS, 110, 209 */
        PutCS(3.1,TauLines[ipTAr7]);
        atom_level2(TauLines[ipTAr7]);
 
        /* A III 7136, 7751, 3109, 5192 CS
         * >>refer      Ar3     CS      Galavis, M. E., Mendoza, C., & Zeippen, C. J. 1995, A&AS, 111, 347
         * >>chng 97 mar 19, break up into three level atom */
        cs12 = 4.825f;
        cs13 = 0.841f;
        cs23 = (realnum)MIN2(1.30,3.296/(phycon.te10*phycon.te01/phycon.te001/
          phycon.te001));
        a12 = 0.3963f;
        a23 = 2.59f;
        a13 = 3.952f;
        /* POP3(G1,G2,G3,O12,O13,O23,A21,A31,A32,E12,E23,P2,ABUND,GAM2) */
        popn3 = (realnum)(atom_pop3(9.,5.,1.,cs12,cs13,cs23,a12,a13,a23,1.955e4,2.770e4,
          &pop2,dense.xIonDense[ipARGON][2],0.,0.,0.));
        CoolHeavy.c7136 = pop2*a12*2.7e-12;
        CoolHeavy.c5192 = popn3*a23*3.83e-12;
        CoolHeavy.c3109 = popn3*a13*6.40e-12;
        CoolAdd("Ar 3",7136,CoolHeavy.c7136);
        CoolAdd("Ar 3",5192,CoolHeavy.c5192);
        CoolAdd("Ar 3",3109,CoolHeavy.c3109);
        /* add to deriv */
        /* >>chng 01 mar 10, did not have second pair of lines included */
        thermal.dCooldT += CoolHeavy.c7136*1.955e4*thermal.tsq1 + 
                (CoolHeavy.c5192+CoolHeavy.c3109)*4.73e4*thermal.tsq1;
 
        /* Ar III 21.8(J=0,1), 9.0 (J=1,2) mircon lines,
         * >>refer      Ar3     CS      Galavis, M. E., Mendoza, C., & Zeippen, C. J. 1995, A&AS, 111, 347 */
        PutCS(3.1,TauLines[ipTAr9]);
        cs = MIN2(1.384,3.110/(phycon.te10/phycon.te001*phycon.te001));
        PutCS(cs,TauLines[ipTAr22]);
        if( phycon.te < 1e4 )
        {
                cs = 0.671;
        }
        else
        {
                cs = MIN2(0.906,0.150*phycon.te20/phycon.te02/phycon.te02*
                  phycon.te003);
        }
        PutCS(cs,*TauDummy);
        atom_level3(TauLines[ipTAr9],TauLines[ipTAr22],*TauDummy);
 
        /* Argon IV 4711+4740, 7335 lines (O II like)
         * CS from
         * >>refer      Ar4     CS      Zeippen, C. J., Le Bourlot, J., & Butler, K. 1987, A&A, 188, 251
         * >>chng 97 jan 31, increase to full 5 level atom
         * Ar IV, cs data from 
         * >>refer      Ar4     CS      Ramsbottom, C. A., Bell., K. L., & Keenan, F. P. 1997, MNRAS 284, 754
         * differs by 2-3x from older values
         * temp dependence form 
         *>>refer       Ar4     CS      Ramsbottom, C. A., & Bell, K. L. 1997, At. Data Nucl. Data Tables, 66, 65 */
        cs12 = (realnum)MAX2(0.761,0.481*phycon.te05);
        cs12 = (realnum)MIN2(0.853,cs12);
        a12 = 2.23e-2f;
 
        cs13 = (realnum)MAX2(1.14,0.719*phycon.te05);
        cs13 = (realnum)MIN2(1.3,cs13);
        a13 = 1.77e-3f;
 
        cs14 = (realnum)MAX2(0.39,0.108*phycon.te10*phycon.te02*phycon.te02);
        cs14 = (realnum)MIN2(0.5,cs14);
        a14 = 0.862f;
 
        cs15 = (realnum)MAX2(0.78,0.216*phycon.te10*phycon.te02*phycon.te02);
        cs15 = (realnum)MIN2(1.0,cs15);
        a15 = 2.11f;
 
        cs23 = 7.06f;
        a23 = 2.30e-5f;
 
        cs24 = (realnum)MAX2(1.53,0.346*phycon.te10*phycon.te05);
        cs24 = (realnum)MIN2(1.96,cs24);
        a24 = 0.603f;
 
        cs25 = (realnum)MAX2(2.18,0.664*phycon.te10*phycon.te02);
        cs25 = (realnum)MIN2(2.65,cs25);
        a25 = 0.789f;
 
        cs34 = (realnum)MAX2(1.56,0.475*phycon.te10*phycon.te02);
        cs34 = (realnum)MIN2(1.89,cs34);
        a34 = 0.119f;
 
        cs35 = (realnum)MAX2(4.01,1.00*phycon.te10*phycon.te02*phycon.te02);
        cs35 = (realnum)MIN2(5.03,cs35);
        a35 = 0.598f;
 
        cs45 = (realnum)(0.0359*phycon.te20*phycon.te20*phycon.te02*phycon.te02);
        a45 = 4.94e-5f;
 
        /* FIVEL( G(1-5) , ex(wn,1-5), cs12,cs13,14,15,23,24,25,34,35,45,
         *  A21,31,41,51,32,42,52,43,53,54, pop(1-5), abund) */
        double Cooling;
        double CoolingDeriv;
        atom_pop5(gAr4,exAr4,cs12,cs13,cs14,cs15,cs23,cs24,cs25,cs34,cs35,
          cs45,a12,a13,a14,a15,a23,a24,a25,a34,a35,a45,p,
          dense.xIonDense[ipARGON][3],&Cooling , &CoolingDeriv, 0.,0.,0.,0.);
        CoolHeavy.Ar4740 = p[1]*a12*4.20e-12;
        CoolHeavy.Ar4711 = p[2]*a13*4.20e-12;
        CoolHeavy.Ar2868 = p[3]*a14*6.94e-12;
        CoolHeavy.Ar2854 = p[4]*a15*6.94e-12;
        CoolHeavy.Ar7263 = p[3]*a24*2.74e-12;
        CoolHeavy.Ar7171 = p[4]*a25*2.74e-12;
        CoolHeavy.Ar7331 = p[3]*a34*2.74e-12;
        CoolHeavy.Ar7237 = p[4]*a35*2.74e-12;
 
        // total cooling from 5-level atom
        CoolAdd("Ar 4",4740,Cooling);
        thermal.dCooldT +=  CoolingDeriv;
 
        /* Argon V 6435+7007,  
         * >>refer      Ar5     AS      Mendoza, C., & Zeippen, C. J. 1982, MNRAS, 199, 1025
         * >>refer      Ar5     CS      Galavis, M. E., Mendoza, C., & Zeippen, C. J. 1995, A&AS, 111, 347
         * POPEXC( O12,g1,g2,A21,excit,abund); result already*a21, excit in Kelvin */
        if( phycon.te < 1e4 )
        {
                cs12 = 3.09f;
        }
        else
        {
                cs12 = (realnum)MIN2(4.454,0.634*phycon.te20/phycon.te03*phycon.te001*
                  phycon.te001);
        }
        cs13 = 0.56f;
        cs23 = 1.65f;
        a12 = 0.68f;
        a13 = 6.55f;
        a23 = 3.35f;
 
        /* >>chng 01 mar 10, convert from 2 to 3 level atom */
        popn3 = (realnum)(atom_pop3(9.,5.,1.,cs12,cs13,cs23,a12,a13,a23,2.055e4,3.110e4,
          &pop2,dense.xIonDense[ipARGON][4],0.,0.,0.));
        CoolHeavy.c7007 = pop2*a12*2.84e-12;
        CoolHeavy.c4626 = popn3*a23*4.30e-12;
        CoolHeavy.c2691 = popn3*a13*7.39e-12;
 
        CoolAdd("Ar 5",7007,CoolHeavy.c7007);
        CoolAdd("Ar 5",4626,CoolHeavy.c4626);
        CoolAdd("Ar 5",2691,CoolHeavy.c2691);
 
        /* add to deriv */
        /* >>chng 01 mar 10, did not have second pair of lines included */
        thermal.dCooldT += CoolHeavy.c7007*2.055e4*thermal.tsq1 + 
                (CoolHeavy.c4626+CoolHeavy.c2691)*5.17e4*thermal.tsq1;
 
        /* Ar V 3P fine structure lines , A from
         * >>refer      Ar5     AS      Mendoza, C. 1982, in IAU Symp. 103, Planetary 
         * >>refercon   Nebulae, ed. D. R. Flower, (Dordrecht, Holland: D. Reidel Publishing Co.), 143
         * >>refer      Ar5     CS      Galavis, M. E., Mendoza, C., & Zeippen, C. J. 1995, A&AS, 111, 347 */
        cs = MIN2(3.26,26.27/(phycon.te20*phycon.te03*phycon.te01*
          phycon.te003));
        PutCS(cs,TauLines[ipTAr13]);
 
        cs = MIN2(8.47,44.31/(phycon.te20/phycon.te003/phycon.te003));
        PutCS(cs,TauLines[ipTAr8]);
 
        cs = MIN2(1.95,7.280/(phycon.te10*phycon.te05*phycon.te005));
        PutCS(cs,*TauDummy);
 
        atom_level3(TauLines[ipTAr13],TauLines[ipTAr8],*TauDummy);
 
        /* [Ar VI] 4.53 micron, cs from 
         * >>refer      Ar6     CS      Saraph, H. E., & Storey, P. J. 1996, A&AS, 115, 151
         * >>chng 96 dec 11, cs should have been 6.33, caught by Simon Casassus */
        PutCS(6.33,TauLines[ipAr06453]);
 
        atom_level2(TauLines[ipAr06453]);
 
        return;
}