1- Grains with THIN ice mantles. 0.5 times the volume of the refractory components / File --> thin5 (coagulated at 105 cm-3)
2- Grains with THICK ice mantles. 4.5 times the volume of the refractory components/ File --> thick5 (coagulated at 105 cm-3)
import numpy as np
import matplotlib.pyplot as plt
data1=np.loadtxt('dustkappa_1.inp',skiprows=2)
data2=np.loadtxt('dustkappa_osshenn_thick5.inp',skiprows=2)
data3=np.loadtxt('dustkappa_osshenn_thin5.inp',skiprows=2)
data4=np.loadtxt('dustkappa_5.inp',skiprows=2)
Data:
Jena Magnesium-iron silicates (glassy) (Jäger et al. 1994, Dorschner et al. 1995)) Mg(0.7) Fe(0.3) SiO(3) [3.01 g/ccm] / File --> pyrmg70.lnk 5 different grain sizes pyrmg70.lnk "MIE" 1 0.0 -0.99999999 -0.99999999 1 -3.5 3.01000 -2.0 pyrmg70.lnk "MIE" 1 0.0 6.4714920e-09 6.4714920e-09 1 -3.5 3.01000 -2.0 pyrmg70.lnk "MIE" 1 0.0 1.0000000 1.0000000 1 -3.5 3.01000 -2.0 pyrmg70.lnk "MIE" 1 0.0 2.0000000 2.0000000 1 -3.5 3.01000 -2.0 pyrmg70.lnk "MIE" 1 0.0 3.0000000 3.0000000 1 -3.5 3.01000 -2.0
Ossenkopf-Henning grains
1- Grains with THIN ice mantles. 0.5 times the volume of the refractory components / File --> thin5 (coagulated at 105 cm-3) 2- Grains with THICK ice mantles. 4.5 times the volume of the refractory components/ File --> thick5 (coagulated at 105 cm-3)
import numpy as np import matplotlib.pyplot as plt data1=np.loadtxt('dustkappa_1.inp',skiprows=2) data2=np.loadtxt('dustkappa_osshenn_thick5.inp',skiprows=2) data3=np.loadtxt('dustkappa_osshenn_thin5.inp',skiprows=2) data4=np.loadtxt('dustkappa_5.inp',skiprows=2)
plot wavelenght v/s kappa_abs in cm^2/gram
plt.plot(data1[:,0],data1[:,1],color='b',label ='$q=-0.991$') plt.plot(data2[:,0],data2[:,1],color='g',label ='$Thin-IM$') plt.plot(data3[:,0],data3[:,1],color='r',label ='$Thick-IM$') plt.plot(data4[:,0],data4[:,1],color='m',label ='$q=3.000$') plt.legend() plt.xscale('log') plt.yscale('log') plt.xlabel("$\lambda\; [\mu\mathrm{m}]$") plt.ylabel("$\kappa\; [\mathrm{cm}^2/\mathrm{g}]$") plt.show()
Paper: http://articles.adsabs.harvard.edu/cgi-bin/nph-iarticle_query?1994A%26A...291..943O&data_type=PDF_HIGH&whole_paper=YES&type=PRINTER&filetype=.pdf