Optical constants:

Si  : Draine2003

aC  : Zubko  ACH2 mix

Ice : Hudgins+93 of ice at T=10K +extend wavelength range using Warren+08  


; Volume fraction of AGN dust model using fluffy grains: with and
; without ice mas ratio: M(aC+Si)/M(ice) =1.
 
; Porosity for aC and Si grains is Vv=10% and no extra porosity is
; assumed for ice.
 
Vv      = 0.1
rhice   = 1.0
rhc     = 1.6
rhsi    = 2.7
wmolc   = 12.
wmolsi  = 134.5
 
 
; Table 2 in LIPSIII for wmolsi = 135 the mass of the MRN grains aC
; and aSi grains has a ratio m(Si):m(C) = 10:20, hence :
msc  = 0.5  ;   =: (Vsi* rhsi) / (Vc * rhc)
 
;
Vc  = (1.-Vv) / ( 1 + msc * rhc/rhsi)
Vsi =  1. - Vv - Vc
 
if Vc+Vsi+Vv ne 1. then print, ' check total volume Vc+Vsi+Vv ', 100*(Vc+Vsi+Vv)
 
 
 
print, ' Volume ratios (%) of fluffy grains without ice:
print, 'Vc, Vsi, Vc, Vv, Vtot'
print, [Vsi,  VaC, Vv , Vc+Vsi+Vv  ] * 100
;       .21   .69  .10  = 100
 
; In this volume the dust mass is:
 
Mc  = Vc * rhc
Msi = Vsi * rhsi
Mv  = 0.
 
 
 
print, ' Mc, Msi, Mc + Msi'
print,   Mc, Msi, Mc + Msi
 
 
; the volumne ratios in 1 g of dust are therefore:
Vc  = Mc/rhc
Vsi = Msi/rhsi
 
 
 
print, ' Check shall be as above:
print, ' Volume ratios (%) of fluffy grains without ice:
print, 'Vc, Vsi, Vv, Vtot'
print, [Vsi,  VaC, Vv , Vc+Vsi+Vv  ] * 100
;       .21   .69  .10  = 100
 
 
; For euqal amount of mass in eis m_eis as for the fluffy grains:
mice = Mc+Msi
 
Vice = mice/rhice
 
; renormalize to totla volume of 100%
Vnorm = Vice + Vv + Vc + Vsi
 
print, ' Volume ratios (%) of fluffy grains with ice: '
print, [Vsi,  VaC, Vv , Vice, Vc+Vsi+Vv+Vice  ] * 100
;       .08  .26   .4  .62  = 100

 
 
end