Ionization Map of 4U0114+65

Ionization Map and Area¶

In this example, we will create an ionization map using real data from the HMXRB 4U 0114+65. During observations across four spectra, differences in luminosity and the Fe Kalpha line were noted.

The Fe Kalpha line manifests when the ionization parameter is below 100. We will investigate the relationship between the emitting area, characterized by an ionization parameter between 25 and 100, and the intensity (flux) of the Fe Kalpha line.

In [22]:

Copied!

import xraybinaryorbit
from xraybinaryorbit import *
import xraybinaryorbit
from xraybinaryorbit import *

In [23]:

Copied!





luminosity = np.array([1.30628470e+36, 6.24733539e+35, 9.38737128e+34, 4.12532189e+35])

#Lower bound, value, higher bound
feka_data = [
    [0.000054, 0.000105, 0.000157],
    [0.000004, 0.000030, 0.000056],
    [0.000000, 0.000007, 0.000015],
    [0.000012, 0.000033, 0.000055]
]

feka = pd.DataFrame(feka_data)
luminosity = np.array([1.30628470e+36, 6.24733539e+35, 9.38737128e+34, 4.12532189e+35])

#Lower bound, value, higher bound
feka_data = [
    [0.000054, 0.000105, 0.000157],
    [0.000004, 0.000030, 0.000056],
    [0.000000, 0.000007, 0.000015],
    [0.000012, 0.000033, 0.000055]
]

feka = pd.DataFrame(feka_data)

In [24]:

Copied!

params=[0.1,1.53,0.,100.0,37.0,1200.0,1.507962e-06,0.8,1.30628470e+36,50,50000]

chi1,rho1,area1 = ionization_map_phase(size_in_Rstar=3, min_color=None,max_color=None, save_plot=False, name="ionization_map"
                                  ,load_directly=False,parameter_list=params)
params=[0.1,1.53,0.,100.0,37.0,1200.0,1.507962e-06,0.8,1.30628470e+36,50,50000]

chi1,rho1,area1 = ionization_map_phase(size_in_Rstar=3, min_color=None,max_color=None, save_plot=False, name="ionization_map"
                                  ,load_directly=False,parameter_list=params)

max color coefficient is 101.08
min color coefficient is 2.93

No description has been provided for this image

In [25]:

Copied!

params=[0.1,1.53,0,100.0,37.0,1200.0,1.507962e-06,0.8,6.24733539e+35,25.0,100.0]

chi2,rho2,area2 = ionization_map_phase(size_in_Rstar=3, min_color=5.86, max_color=209.98, 
                                  save_plot=False, name="ionization_map",parameter_list=params)
params=[0.1,1.53,0,100.0,37.0,1200.0,1.507962e-06,0.8,6.24733539e+35,25.0,100.0]

chi2,rho2,area2 = ionization_map_phase(size_in_Rstar=3, min_color=5.86, max_color=209.98, 
                                  save_plot=False, name="ionization_map",parameter_list=params)

In [26]:

Copied!

params=[0.1,1.53,0.,100.0,37.0,1200.0,1.507962e-06,0.8,9.38737128e+34,25.0,100.0]

chi3,rho3,area3= ionization_map_phase(size_in_Rstar=3, min_color=5.86, max_color=109.98, 
                                  save_plot=False, name="ionization_map",parameter_list=params)
params=[0.1,1.53,0.,100.0,37.0,1200.0,1.507962e-06,0.8,9.38737128e+34,25.0,100.0]

chi3,rho3,area3= ionization_map_phase(size_in_Rstar=3, min_color=5.86, max_color=109.98, 
                                  save_plot=False, name="ionization_map",parameter_list=params)

In [27]:

Copied!

params=[0.1,1.53,0,100.0,37.0,1200.0,1.507962e-06,0.8,4.12532189e+35,25.0,100.0]

chi4,rho4,area4 = ionization_map_phase(size_in_Rstar=20, min_color=5.86, max_color=109.98, 
                                  save_plot=False, name="ionization_map",parameter_list=params)
params=[0.1,1.53,0,100.0,37.0,1200.0,1.507962e-06,0.8,4.12532189e+35,25.0,100.0]

chi4,rho4,area4 = ionization_map_phase(size_in_Rstar=20, min_color=5.86, max_color=109.98, 
                                  save_plot=False, name="ionization_map",parameter_list=params)

In [28]:

Copied!





total_area= np.array([area1,area2,area3,area4])

yerr = np.transpose(np.diff(feka))
plt.figure(figsize=(5, 8))
plt.errorbar(total_area,feka[1], yerr=yerr, fmt=".")
plt.xlabel(r"Area where $\xi$ 25-100 (R$_{*}^{2}$)")
plt.ylabel(r"Fe k$_{\alpha}$ Area (photons s$^{-1}$ cm$^{-2}$)")
total_area= np.array([area1,area2,area3,area4])

yerr = np.transpose(np.diff(feka))
plt.figure(figsize=(5, 8))
plt.errorbar(total_area,feka[1], yerr=yerr, fmt=".")
plt.xlabel(r"Area where $\xi$ 25-100 (R$_{*}^{2}$)")
plt.ylabel(r"Fe k$_{\alpha}$ Area (photons s$^{-1}$ cm$^{-2}$)")

Out[28]:

Text(0, 0.5, 'Fe k$_{\\alpha}$ Area (photons s$^{-1}$ cm$^{-2}$)')

In [ ]: