Intensity vs distance for grid on phiH

[SUVENDU RAKSHIT]

Hello,

I have a question. It may be very stupid one.

I will be very happy if someone help me.

My model is:

table power law

hden=10  

stop column density 23

phi(H)=17 to 25 in steps of 0.5(say) # A small grid on ionized photon flux

..........

I want to have a plot of intensity of a line vs distance(r) from the central source for a object of known Luminosity(L)? How I can get that?

r=sqrt(Q/4*pi*phi(H)). Is it possible to have relation between Q and L from pycloudy or do I need to calculate it from the cloudy manual definition( considering the continuum shape in my case table power law)?

Thank you very much.

Regards,

Suvendu

[Christophe Morisset]

Ops, sorry for the delay, I guess you have found yourself a solution... The relation between Q and L is not that obvious, given that both are integrals of the flux, but not on the same boundaries, so it depends on the shape of the SED, and thus on the T* or the power index. You can always run a 1-zone Cloudy model and look at the first informations from the output file. With a few values of the index or T* you can make a table and then interpolate to have the relation between Q and L for any value.

Hope it helps,

Christophe

[SUVENDU RAKSHIT]

Hello, 

Thank you very much. I run a 1 zone model as you suggested 

table power law

phiH=17;

nH=10;

stop zone 1;

but I am unable to figure out where is the L  and Q in the output file!!

(The question I asked you before to have a plot of radial emissivity with radius for a grid over phiH and nH.

It would be very easy if I understand this sentence "For the adopted mean ionizing continuum source luminosity 

L_ion=10^44.14 erg/s, log(phiH)=20 corresponds to a source distance of R=15(H0/70km/s) light-days." 

I am missing a factor 2. 

r=sqrt[ L/(4pi*phiH*h*nu ] 

lam=91.2 nm)

Thank you.

Regards,

Suvendu

[Christophe Morisset]

Hi,

At the very beginning of the output file, you will find something like this (see below), where L, Q and a lot of other informations are given. Part of this information is read by pyCloudy and you can also have it from the CloudyModel object. Obviously, you cannot have Q and Phi at the same time, it depends on the geometry (spherical or plan parallel).

           4171CellPeak1.00E+00   Lo 1.00e-08=910.79cm   Hi-Con:1.14E+01 Ryd   E(hi):7.35E+06Ryd     E(hi):     100.01 MeV

           L(nu>1ryd):  39.3414   Average nu:1.489E+00   L( X-ray):   0.0000   L(BalC):  39.1931     Q(Balmer C):  50.0693

           Q(1.0-1.8):  49.7350   Q(1.8-4.0):  49.1213   Q(4.0-20):  46.7256   Q(20--):   0.0000     Ion pht flx:5.379E+14

           L(gam ray):   0.0000   Q(gam ray):   0.0000   L(Infred):  37.8505   Alf(ox):   0.0000     Total lumin:  39.5827

           log L/Lsun:   5.9999   Abs bol mg: -10.2497   Abs V mag:  -6.0267   Bol cor:  -4.2230     nuFnu(Bbet):  37.9523

           U(1.0----):1.794E+04   U(4.0----):1.411E+01   T(En-Den):1.076E+02   T(Comp):4.787E+04     nuJnu(912A):2.070E+04

           Occ(FarIR):6.777E-04   Occ(H n=6):2.344E-10   Occ(1Ryd):9.501E-13   Occ(4R):6.963E-17     Occ (Nu-hi):5.695E-27

           Tbr(FarIR):1.074E-06   Tbr(H n=6):1.027E-06   Tbr(1Ryd):1.502E-07   Tbr(4R):4.400E-11     Tbr (Nu-hi):1.021E-20