cool_magn.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 */
/*CoolMagn compute magnesium cooling */
#include "cddefines.h"
#include "coolheavy.h"
#include "taulines.h"
#include "phycon.h"
#include "dense.h"
#include "ligbar.h"
#include "lines_service.h"
#include "atoms.h"
#include "cooling.h"
 
void CoolMagn(void)
{
        realnum cs, 
          csdum, 
          csoi;
        double cs2s2p, 
          cs2s3p;
 
        DEBUG_ENTRY( "CoolMagn()" );
 
        /* Mg I 2853 
         * fit to Dima's  integration of
         * >>refer      mg1     cs      Leep, D., & Gallagher, A. 1976, Phys Rev A, 13, 148 */
        cs = (realnum)(5.21e-4*phycon.te/phycon.te10);
        PutCS(cs,TauLines[ipMgI2853]);
        atom_level2(TauLines[ipMgI2853]);
 
        /* Mg I 2026,  */
        MakeCS(TauLines[ipMgI2026]);
        atom_level2(TauLines[ipMgI2026]);
 
        /* Mg 1 4571, data from Mendoza
         * >>refer      mg1     as      Mendoza, C. 1982, in Planetary Nebulae, IAU Symp No. 103,
         * >>refercon ed by D.R. Flower, (D. Reidel: Holland), 143
         * cs set to O I */
        csoi = (realnum)(2.68e-5*phycon.te*(1. + 1.67e-6*phycon.te - 2.95e-10*phycon.te*
          phycon.te));
        csoi = (realnum)MAX2(0.1,csoi);
        PutCS(csoi/9.,TauLines[ipT4561]);
        atom_level2(TauLines[ipT4561]);
 
        /* Mg II 2798
         * cs from 
         * >>refer      mg2     cs      Sigut, A., & Pradhan, A.K., 1994, J Phys B sub
         * refer not publ as of '97, cs agrees fairly well with Harrington et al.
         * previous reference (~'82) */
        cs2s2p = 4.50*phycon.te10;
        PutCS(cs2s2p,TauLines[ipT2796]);
        PutCS(cs2s2p*0.5,TauLines[ipT2804]);
        PutCS(1.0,*TauDummy);
        atom_level3(TauLines[ipT2804],*TauDummy,TauLines[ipT2796]);
        /* call PutCS( cs , t2800 )
         * call atom_level2( t2800 )
         *
         * following used in MAGNES for photo destruction rate */
        if( atoms.PopLevels[0] > 0. )
        {
                atoms.popmg2 = (realnum)((atoms.PopLevels[2] + atoms.PopLevels[1])/
                  atoms.PopLevels[0]);
        }
        else
        {
                atoms.popmg2 = 0.;
        }
 
        /* MG IV 4.487 MIC 
         * cs 
         * >>referold   mg4     cs      Saraph, H.E. & Tully, J.A. 1994, A&AS, 107, 29 */
        /*cs = (realnum)MIN2(0.425,0.180*phycon.te05*phycon.te02);
        cs = (realnum)MAX2(0.356,cs);*/
        /* >>chng 06 jul 06-Humeshkar Nemala*/
        /*>> refer Mg IV cs Berrington,K.A., Saraph,H. E., & Tully, J.A. 1998,A&AS,129,161 */
        /*This is the cs of the transition between the levels of 2P^o term(J=1/2 - J=3/2)*/
        if(phycon.te < 4E5)
        {
                cs = (realnum)(0.155*phycon.te07*phycon.te01*
                        phycon.te003*phycon.te0005*phycon.te0001);
        }
        else
        {
                cs = (realnum)(5.124/((phycon.te20/phycon.te02)*
                        phycon.te007*phycon.te0005*phycon.te0001));
        }
        PutCS(cs,TauLines[ipTMg4]);
 
        atom_level2(TauLines[ipTMg4]);
 
        /* MG V 5.61, 13.54 micron, 
         * >>refer      mg5     cs      Butler, K., & Zeippen, C.J. 1994, A&AS, 108, 1
         * >>chng 96 jul 16 had been constant 0.3 */
        cs = (realnum)MIN2(0.311,0.11*phycon.te10);
        PutCS(cs,TauLines[ipTMg14]);
 
        cs = (realnum)MIN2(1.06,0.339*phycon.te10);
        PutCS(cs,TauLines[ipTMg6]);
        /* >>chng 96 jul 16 had been constant 0.3 */
        cs = (realnum)MIN2(0.297,0.0745*phycon.te10*phycon.te02);
        PutCS(cs,*TauDummy);
        atom_level3(TauLines[ipTMg6],TauLines[ipTMg14],*TauDummy);
 
        /* [Mg V] 2751+2893- Ne III-like, cs 
         * >>refer      mg5     as      Mendoza, C., & Zeippen, C.J. 1987, MNRAS, 224, 7p
         * >>chng 96 aug 5 to three level atom
         * c2751 = atom_pop2(1.33,9.,5.,2.4,5.14e4,xmg(5))*7.11e-12
         * dCooldT = dCooldT + c2751*5.14e4*tsq1
         * call CoolAdd( 'Mg 5' , 2751 , C2751 )
         *
         * following is 2-1 transition, both 2928 and 2783, 
         * >>refer      mg5     cs      Butler, K., & Zeippen, C.J. 1994, A&AS, 108, 1 */
        PutCS(1.187,TauLines[ipxMg52855]);
 
        cs = (realnum)MIN2(0.278,0.0171*phycon.te20*phycon.te05/
          phycon.te005/phycon.te003);
        cs = (realnum)MAX2(0.182,cs);
 
        /* 3-2 transition, 2417.5 */
        PutCS(cs,TauLines[ipxMg52417]);
 
        /* 3-1 transition, 1324.58 */
        PutCS(0.153,TauLines[ipxMg51325]);
 
        atom_level3(TauLines[ipxMg52855],TauLines[ipxMg52417],TauLines[ipxMg51325]);
 
        /* Mg VI, 1806- OII like, data 
         * >>refer      mg6     all     Kafatos, M., & Lynch, J.P. 1980, ApJS, 42, 611 */
        /* >>refer      mg6     as      Becker, Butler, Zeippen, 1989, A&A 221, 375 
         * >>refer      mg6     cs      Ramsbottom & Bell 1997, A&AS 125, 543 */
        CoolHeavy.c1806 = atom_pop2(0.6,4.,10.,0.1,7.974e4,dense.xIonDense[ipMAGNESIUM][5])*
          1.11e-11;
        CoolAdd("Mg 6",1806,CoolHeavy.c1806);
 
        /* [Mg VII] IR lines at 5.517 and 9.03 microns, 
         * carbon-like, 
         * >>refer      mg7     cs      Lennon, D.J. Burke, V.M. 1994, A&AS, 103, 273 */
        if( phycon.alogte < 4.4 )
        {
                cs = (realnum)(0.027*phycon.te30/phycon.te03*phycon.te003*phycon.te001);
        }
        else
        {
                cs = 0.44f;
        }
        PutCS(cs,TauLines[ipfsMg790]);
        if( phycon.alogte < 4.6 )
        {
                cs = (realnum)(MIN2(1.456,0.0577*phycon.te30*phycon.te02/phycon.te001/
                  phycon.te001));
                csdum = (realnum)(8.275e-3*phycon.sqrte/phycon.te10/phycon.te001);
        }
        else
        {
                cs = (realnum)(3.257/(phycon.te05*phycon.te02*phycon.te003*phycon.te003));
                csdum = (realnum)(1.456/(phycon.te10*phycon.te01*phycon.te005));
        }
        PutCS(cs,TauLines[ipfsMg755]);
        PutCS(csdum,*TauDummy);
        /* atom_level3(  t10,t21,t20) */
        atom_level3(TauLines[ipfsMg790],TauLines[ipfsMg755],*TauDummy);
 
        /* [mg vii] 2510, 2629
         * c2596 = atom_pop2( 1.7,9.,5.,10.,5.76e4,dense.xIonDense(12,7))*7.96e-12
         * dCooldT = dCooldT + c2596 * 5.76e4*tsq1
         * call CoolAdd( 'Mg 7' , 2596 , C2596 )
         *
         * >>chng 96 aug 5, converted to 3 level atom */
        cs = (realnum)MIN2(0.22,0.3622/(phycon.te05*phycon.te02*phycon.te003));
        PutCS(cs,TauLines[ipxMg71190]);
 
        /* 2-1 transitions, 2509.2A+2629.1A together */
        cs = (realnum)MIN2(1.067,0.247*phycon.te10*phycon.te03*phycon.te005);
        PutCS(cs,TauLines[ipxMg72569]);
 
        /* 3-2 transition, 2261.5A */
        cs = (realnum)MIN2(0.542,3.863/(phycon.te20*phycon.te03*
          phycon.te01/phycon.te003));
        cs = (realnum)MAX2(0.3735,cs);
        PutCS(cs,TauLines[ipxMg72261]);
 
        atom_level3(TauLines[ipxMg72569],TauLines[ipxMg72261],TauLines[ipxMg71190]);
        /* atom_level3(  t10,t21,t20)
         *
         * Mg VIII 3.03 micron, data from 
         * >>refer      mg8     as      Chandra, S. 1982, SoPh, 75, 133
         * cs from 
         * >>refer      mg8     cs      Zhang, H.L., Graziani, M., Pradhan, A.K. 1994, A&A, 283, 319 */
        PutCS(1.0,TauLines[ipxMg08303]);
        atom_level2(TauLines[ipxMg08303]);
        /* fs303 = atom_pop2(0.26,2.,4.,0.324,4752.,xmg(8))*6.58E-13
         * call CoolAdd( 'Mg 8' , 303 , FS303 )
         *
         * Mg IX 704.5, 1909-like, A from
         * >>refer      mg9     as      Muhlethaler, H.P., & Nussbaumer, H. 1976, A&A 48, 109
         * AtomSeqBeryllium line, cs data from 
         * >>refer      mg9     cs      Keenan, F.P. Berrington, K.A., Burke, P.G., Dufton, P.L.,
         * >>refercon Kingston, A.E. 1986, PhyS 34, 216
         * A's 
         * >>refer      mg9     as      Fleming, J., Bell, K.L, Hibbert, A., Vaeck, N., Godefroid, M.R.
         * >>refercon 1996, MNRAS, 279, 1289 */
        cs = (realnum)(0.98288 - 0.23766*phycon.alogte + 0.014334*POW2(phycon.alogte));
        cs = (realnum)MAX2(0.01,cs);
        PutCS(cs,TauLines[ipT705]);
        atom_level2(TauLines[ipT705]);
 
        /* Mg X 610
         * >>refer      mg10    cs      Cochrane, D.M., & McWhirter, R.W.P. 1983, PhyS, 28, 25 */
        ligbar(12,TauLines[ipTMg610],TauLines[ipT58],&cs2s2p,&cs2s3p);
        PutCS(cs2s2p,TauLines[ipTMg610]);
        PutCS(cs2s2p*0.5,TauLines[ipTMg625]);
        PutCS(1.0,*TauDummy);
        atom_level3(TauLines[ipTMg625],*TauDummy,TauLines[ipTMg610]);
 
        /* Mg X 58A */
        PutCS(cs2s3p,TauLines[ipT58]);
        atom_level2(TauLines[ipT58]);
        return;
}