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_70.inp',skiprows=2)
data2=np.loadtxt('dustkappa_1_100.inp',skiprows=2)
data3=np.loadtxt('dustkappa_osshenn_thick5.inp',skiprows=2)
data4=np.loadtxt('dustkappa_osshenn_thin5.inp',skiprows=2)
cagurto
commented
7 years ago 
Data: Using pyrmg100.lnk (Magnesium-iron silicates (glassy)/Mg SiO(3) [2.71 g/ccm]) pyrmg70.lnk (Magnesium-iron silicates (glassy)/Mg(0.7) Fe(0.3) SiO(3) [3.01 g/ccm]) from http://www.astro.uni-jena.de/Laboratory/OCDB/amsilicates.html
and
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_70.inp',skiprows=2) data2=np.loadtxt('dustkappa_1_100.inp',skiprows=2) data3=np.loadtxt('dustkappa_osshenn_thick5.inp',skiprows=2) data4=np.loadtxt('dustkappa_osshenn_thin5.inp',skiprows=2)
plot wavelenght v/s kappa_abs in cm^2/gram
plt.plot(data1[:,0],data1[:,1],color='b',label ='Amorphous Silicates (3.01 gr/cm) $q=3.5$') plt.plot(data2[:,0],data2[:,1],color='g',label ='Amorphous Silicates (2.71 gr/cm) $q=3.5$') plt.plot(data3[:,0],data3[:,1],color='r',label ='THIN ice mantles $q=3.5$') plt.plot(data4[:,0],data4[:,1],color='c',label ='THICK ice mantles $q=3.5$') plt.legend() plt.xscale('log') plt.yscale('log') axes = plt.gca() axes.set_xlim([1,1000]) axes.set_ylim([1e-2,1e10]) plt.xlabel("$\lambda\; [\mu\mathrm{m}]$") plt.ylabel("$\kappa\; [\mathrm{cm}^2/\mathrm{g}]$") plt.show()
import numpy as np import matplotlib.pyplot as plt data1=np.loadtxt('dustkappa_5_70.inp',skiprows=2) data2=np.loadtxt('dustkappa_5_100.inp',skiprows=2) data3=np.loadtxt('dustkappa_osshenn_thick5.inp',skiprows=2) data4=np.loadtxt('dustkappa_osshenn_thin5.inp',skiprows=2)
plot wavelenght v/s kappa_abs in cm^2/gram
plt.plot(data1[:,0],data1[:,1],color='b',label ='Amorphous Silicates (3.01 gr/cm) $q=3.5$') plt.plot(data2[:,0],data2[:,1],color='g',label ='Amorphous Silicates (2.71 gr/cm) $q=3.5$') plt.plot(data3[:,0],data3[:,1],color='r',label ='THIN ice mantles $q=3.5$') plt.plot(data4[:,0],data4[:,1],color='c',label ='THICK ice mantles $q=3.5$') plt.legend() plt.xscale('log') plt.yscale('log') axes = plt.gca() axes.set_xlim([1,1000]) axes.set_ylim([1e-2,1e10]) plt.xlabel("$\lambda\; [\mu\mathrm{m}]$") plt.ylabel("$\kappa\; [\mathrm{cm}^2/\mathrm{g}]$") plt.show()
import numpy as np import matplotlib.pyplot as plt data1=np.loadtxt('dustkappa_3_70.inp',skiprows=2) data2=np.loadtxt('dustkappa_3_100.inp',skiprows=2) data3=np.loadtxt('dustkappa_osshenn_thick5.inp',skiprows=2) data4=np.loadtxt('dustkappa_osshenn_thin5.inp',skiprows=2)
plot wavelenght v/s kappa_abs in cm^2/gram
plt.plot(data1[:,0],data1[:,1],color='b',label ='Amorphous Silicates (3.01 gr/cm) $q=3.5$') plt.plot(data2[:,0],data2[:,1],color='g',label ='Amorphous Silicates (2.71 gr/cm) $q=3.5$') plt.plot(data3[:,0],data3[:,1],color='r',label ='THIN ice mantles $q=3.5$') plt.plot(data4[:,0],data4[:,1],color='c',label ='THICK ice mantles $q=3.5$') plt.legend() plt.xscale('log') plt.yscale('log') axes = plt.gca() axes.set_xlim([1,1000]) axes.set_ylim([1e-2,1e10]) plt.xlabel("$\lambda\; [\mu\mathrm{m}]$") plt.ylabel("$\kappa\; [\mathrm{cm}^2/\mathrm{g}]$") plt.show()
import numpy as np import matplotlib.pyplot as plt data1=np.loadtxt('dustkappa_2_70.inp',skiprows=2) data2=np.loadtxt('dustkappa_2_100.inp',skiprows=2) data3=np.loadtxt('dustkappa_osshenn_thick5.inp',skiprows=2) data4=np.loadtxt('dustkappa_osshenn_thin5.inp',skiprows=2)
plot wavelenght v/s kappa_abs in cm^2/gram
plt.plot(data1[:,0],data1[:,1],color='b',label ='Amorphous Silicates (3.01 gr/cm) $q=3.5$') plt.plot(data2[:,0],data2[:,1],color='g',label ='Amorphous Silicates (2.71 gr/cm) $q=3.5$') plt.plot(data3[:,0],data3[:,1],color='r',label ='THIN ice mantles $q=3.5$') plt.plot(data4[:,0],data4[:,1],color='c',label ='THICK ice mantles $q=3.5$') plt.legend() plt.xscale('log') plt.yscale('log') axes = plt.gca() axes.set_xlim([1,1000]) axes.set_ylim([1e-2,1e10]) plt.xlabel("$\lambda\; [\mu\mathrm{m}]$") plt.ylabel("$\kappa\; [\mathrm{cm}^2/\mathrm{g}]$") plt.show()
import numpy as np import matplotlib.pyplot as plt data1=np.loadtxt('dustkappa_4_70.inp',skiprows=2) data2=np.loadtxt('dustkappa_4_100.inp',skiprows=2) data3=np.loadtxt('dustkappa_osshenn_thick5.inp',skiprows=2) data4=np.loadtxt('dustkappa_osshenn_thin5.inp',skiprows=2)
plot wavelenght v/s kappa_abs in cm^2/gram
plt.plot(data1[:,0],data1[:,1],color='b',label ='Amorphous Silicates (3.01 gr/cm) $q=3.5$') plt.plot(data2[:,0],data2[:,1],color='g',label ='Amorphous Silicates (2.71 gr/cm) $q=3.5$') plt.plot(data3[:,0],data3[:,1],color='r',label ='THIN ice mantles $q=3.5$') plt.plot(data4[:,0],data4[:,1],color='c',label ='THICK ice mantles $q=3.5$') plt.legend() plt.xscale('log') plt.yscale('log') axes = plt.gca() axes.set_xlim([1,1000]) axes.set_ylim([1e-2,1e10]) plt.xlabel("$\lambda\; [\mu\mathrm{m}]$") plt.ylabel("$\kappa\; [\mathrm{cm}^2/\mathrm{g}]$") plt.show()