cool_chlo.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 */
/*CoolChlo compute chlorine cooling */
#include "cddefines.h"
#include "taulines.h"
#include "coolheavy.h"
#include "dense.h"
#include "phycon.h"
#include "thermal.h"
#include "lines_service.h"
#include "atoms.h"
#include "cooling.h"
 
void CoolChlo(void)
{
        double a21, 
          a31, 
          a32, 
          a41, 
          a42, 
          a43, 
          a51, 
          a52, 
          a53, 
          a54, 
          cs, 
          cs01,
          cs02,
          cs12, 
          cs13, 
          cs14, 
          cs15, 
          cs23, 
          cs24, 
          cs25, 
          cs34, 
          cs35, 
          cs45, 
          p[5], 
          p3,
          tused;
        realnum
          p2, 
          pop2,
          rate;
        static double gcl3[5]={4.,4.,6.,2.,4.};
        static double excl3[4]={18053.,66.,11693.,95.};
 
        DEBUG_ENTRY( "CoolChlo()" );
 
        /* >>chng >>03 nov 09, add this line, from
         * >>refer      Cl1     cs      Hollenbach, D. & McKee, C.F. 1989, ApJ, 342, 306 */
        /* rates are said to be ok over range 30 - 3000K */
        tused = MAX2( 30. , phycon.te );
        tused = MIN2( 3000. , phycon.te  );
        tused /= 100.;
 
        /* HM89 give deexcitation rates, must convert into electron collision strength, as expected
         * by the code's infrastructure */
        /* electron collision strength */
        rate = (realnum)( 4.7e-8 * dense.eden + 
                /* >>chng 05 jul 05, eden to cdsqte */
                /*8.3e-10*pow( tused , 0.17 ) * dense.xIonDense[ipHYDROGEN][0]) / dense.eden);*/
                8.3e-10*pow( tused , 0.17 ) * dense.xIonDense[ipHYDROGEN][0] );
        /* possible for atomic hydrogen density to be vary small, causing zero rate coef,
         * which triggers thrown assert - guard against this */
        LineConvRate2CS( TauLines[ipCl1_11m]  , SDIV(rate) );
        atom_level2(TauLines[ipCl1_11m]);
 
        /* [Cl II] 14.3678 mic, 33.281 mic*/
        /* >>refer      cl2     as      Mendoza, C., & Zeippen, C.J. 1983, MNRAS, 202, 981*/
        /* the following cs were about 2x smaller */
        /* >>referold   cl2     cs      Krueger, T.K., & Czyzak, S.J. 1970, Pro Roy Soc Lond, 318, 531 */
        /* >>chgn 03 feb 24, change to following collision strengths */
        /* >>refer      cl2     cs      Wilson, N.J., & Bell, K.L. 2002, MNRAS, 331, 389 */
        /* order of 3P ground term, 2, 1, 0 from lowest to highest */
 
        /* this is the 3P J=1 -> 2, 14.4 micron line */
        cs12 = 17.5 / (phycon.te10*phycon.te04 );
        /*PutCS(2.17,TauLines[ipfsCl214]);*/
        PutCS( cs12 , TauLines[ipfsCl214] );
 
        /* this is the 3P J=0 -> 1, 33.3 micron line */
        cs01 = 4.85 / (phycon.te10*phycon.te02 );
        PutCS( cs01 , TauLines[ipfsCl233] );
        /*PutCS(0.93,TauLines[ipfsCl233]);*/
 
        /* the 0 - 2 transition */
        cs02 = 4.51 / (phycon.te10*phycon.te04 );
        PutCS( cs02 , *TauDummy );
        /*PutCS(1.00,TauDummy);*/
 
        /* atom_level3(  t10,t21,t20) */
        atom_level3(TauLines[ipfsCl214],TauLines[ipfsCl233],*TauDummy);
 
        /* [Cl II] 8578.7, 9123.6, 6161.8, 3677.9 */
        /* >>refer      Cl2     As      Mendoza, C., & Zeippen, C.J. 1983, MNRAS, 202, 981 */
        /* following numbering in terms of level position, 1, 2, then 3 */
        a21 = 0.133;
        a31 = 1.33;
        a32 = 2.06;
        cs13 = 1.01;
        /* this is 10x what is in the paper, as per comment made in 
         * >>refer      Cl2     cs      Keenan, F.P., Aller, L.H., Exter, K.M., Hyung, S., & 
         * >>refercon   Pollacco, D.L. 2003, ApJ, 584, 385 */
        cs23 = 1.49;
        cs12 = 8.389;
        /* POP3(G1,G2,G3,O12,O13,O23,A21,A31,A32,E12,E23,P2,ABUND,GAM2) */
        p3 = atom_pop3(9.,5.,1.,cs12,cs13,cs23,a21,a31,a32,1.576e4,23344.,
          &p2,dense.xIonDense[ipCHLORINE][1],0.,0.,0.);
        /*p3 = atom_pop3(9.,5.,1.,3.86,0.456,1.15,a21,a31,a32,1.576e4,23344.,
          &p2,dense.xIonDense[ipCHLORINE][1],0.);*/
 
        /* [Cl II] 8578.7, 9123.6 doublet, both together */
        CoolHeavy.c8579 = p2*a21*2.32e-12;
        CoolAdd("Cl 2",8579,CoolHeavy.c8579);
 
        /* [Cl II] 6161.8 auroral line */
        CoolHeavy.c6164 = p3*a32*3.23e-12;
        CoolAdd("Cl 2",6164,CoolHeavy.c6164);
 
        /* [Cl II] 3677.9 */
        CoolHeavy.c3679 = p3*a31*5.41e-12;
        CoolAdd("Cl 2",3679,CoolHeavy.c3679);
 
        /* [Cl III] this is a [SII] - like doublet, vac lam=5519, 5539
         * all data from 
         * >>refer      cl3     all     Mendoza, C. 1982, in Planetary Nebulae, IAU Symp No. 103,
         * >>refercon ed by D.R. Flower, (D. Reidel: Holland), 143 */
        cs12 = 1.26;
        a21 = 4.83e-3;
 
        cs13 = 1.88;
        a31 = 7.04e-4;
 
        cs14 = 0.627;
        a41 = 0.305;
 
        cs15 = 1.26;
        a51 = 0.754;
 
        cs23 = 3.19;
        a32 = 3.22e-6;
 
        cs24 = 1.24;
        a42 = 0.303;
 
        cs25 = 1.91;
        a52 = 0.323;
 
        cs34 = 1.38;
        a43 = 0.100;
 
        cs35 = 3.33;
        a53 = 0.316;
 
        cs45 = 1.34;
        a54 = 7.65e-6;
 
        double Cooling , CoolingDeriv;
        atom_pop5(gcl3,excl3,cs12,cs13,cs14,cs15,cs23,cs24,cs25,cs34,cs35,
                cs45,a21,a31,a41,a51,a32,a42,a52,a43,a53,a54,p,
                dense.xIonDense[ipCHLORINE][2]  , &Cooling , &CoolingDeriv, 0.,0.,0.,0.);
 
        CoolHeavy.Cl5539 = p[1]*a21*3.59e-12;
        CoolHeavy.Cl5519 = p[2]*a31*3.61e-12;
        CoolHeavy.Cl3354 = p[3]*a41*5.93e-12;
        CoolHeavy.Cl3344 = p[4]*a51*5.95e-12;
        CoolHeavy.Cl8504 = p[3]*a42*2.34e-12;
        CoolHeavy.Cl8436 = p[4]*a42*2.36e-12;
        CoolHeavy.Cl8552 = p[3]*a43*2.33e-12;
        CoolHeavy.Cl8483 = p[4]*a53*2.35e-12;
 
        /* following are whole multiplets */
        CoolHeavy.c5525 = CoolHeavy.Cl5539 + CoolHeavy.Cl5519;
        CoolHeavy.c3350 = CoolHeavy.Cl3354 + CoolHeavy.Cl3344;
        CoolHeavy.c8494 = CoolHeavy.Cl8504 + CoolHeavy.Cl8436 + CoolHeavy.Cl8552 + 
          CoolHeavy.Cl8483;
 
        // total cooling from 5-level atom
        CoolAdd("Cl 3",5525,Cooling);
        thermal.dCooldT += CoolingDeriv;
 
        /* [CL IV], like [OIII]
         * cs from 
         * >>refer      cl4     cs      Galavis, M.E., Mendoza, C., & Zeippen, C.J. 1995, A&AS, 111, 347 */
        a21 = 0.251;
        cs12 = 6.437;
        a32 = 2.80;
        cs23 = MIN2(2.1,0.0450*phycon.te30*phycon.te03*phycon.te03);
        a31 = 2.50;
        cs13 = 1.922;
        /* POP3(G1,G2,G3,O12,O13,O23,A21,A31,A32,E12,E23,P2,ABUND,GAM2) */
        p3 = atom_pop3(9.,5.,1.,cs12,cs13,cs23,a21,a31,a32,2.24e4,3.11e4,&pop2,
          dense.xIonDense[ipCHLORINE][2],0.,0.,0.);
        /* whole 2-1 transition */
        CoolHeavy.c8047 = pop2*a21*2.48e-12;
        CoolAdd("Cl 4",8047,CoolHeavy.c8047);
        thermal.dCooldT += CoolHeavy.c8047*(2.24e4*thermal.tsq1 - thermal.halfte);
        /* 3-1 transition */
        CoolHeavy.c3119 = p3*a31*6.38e-12;
        CoolAdd("Cl 4",3119,CoolHeavy.c3119);
        /* 3-2 transition */
        CoolHeavy.c5324 = p3*a32*3.74e-12;
        CoolAdd("Cl 4",5324,CoolHeavy.c5324);
 
        /* [Cl IV] fine structure lines, 20.354, 11.741 microns */
        cs = MIN2(2.7,6.637/(phycon.te10*phycon.te03*phycon.te01));
        cs = MAX2(1.6,cs);
        PutCS(cs,TauLines[ipCl04203]);
 
        cs = MIN2(8.0,15.65/phycon.te10);
        PutCS(cs,TauLines[ipCl04117]);
 
        cs = MIN2(2.0,5.805/(phycon.te10*phycon.te03));
        PutCS(cs,*TauDummy);
 
        /* atom_level3(  t10,t21,t20) */
        atom_level3(TauLines[ipCl04203],TauLines[ipCl04117],*TauDummy);
 
        /* fixit - add Cl V 6.71 micron using cs from
         * >>refer      cl5     cs      Saraph, H.E., & Storey, P.J., 1999, A&AS, 134, 369 */
 
        /* [Cl IX] 7335 A, 
         * >>referold   cl9     as      Saraph, H.E. & Tully, J.A. 1994, A&AS, 107, 29 */
        /* >>refer      cl9     as      Berrington,K.A.,Saraph, H.E. & Tully, J.A. 1998, A&AS, 129, 161 */
        /*>>chng 06 jul 18 Changes made-Humeshkar Nemala*/
        if(phycon.te < 2.03E5)
        {
                cs = (realnum)(0.1175*phycon.te07*phycon.te004*phycon.te0001);
        }
        else if(phycon.te < 5.11E5)
        {
                cs = (realnum)((60.7989E-02)/(phycon.te05*phycon.te01*phycon.te0004));
        }
        else if(phycon.te < 1.284E6)
        {
                cs = (realnum)(0.274857); 
        }
        else
        {
                cs = (realnum)(27.327963/(phycon.te30*phycon.te02*phycon.te007));
        }
        /*PutCS(0.28,TauLines[ipCl973]);*/
        PutCS(cs,TauLines[ipCl973]);
        atom_level2(TauLines[ipCl973]);
 
        return;
}