Galactic object models#

Tow models have been added to handle galactic objects: 1) GalacticBeamed to hadle galatcit objects that are not beamed 2) GalacticUnbeamed to hadle galatcit objects that are beamed

The redshift will be set 0, and not used, and a new parameter, DL_cm will be used in place.

The only difference, is that the GalacticBeamed has beaming set as default to 10 and set free, whilst GalacticUnbeamed has beaming set as default to 1 and frozen

from jetset.jet_model import GalacticUnbeamed,GalacticBeamed
from astropy import constants as const
from astropy import units as u
import matplotlib.pyplot as plt
import numpy as np

import jetset
print('tested with',jetset.__version__)

tested with 1.4.0rc0

Galactic Unbeamed leptonic example#

This is an example for a PWN.

Please notice, that for a realistic model you should build a composite model with Radio and Wind region, and define a custom emitting region. A proper thermal external field will be provided in the next version. The ``CMB`` is the same used for jets

pwn=GalacticUnbeamed(emitters_distribution='bkn',verbose=False,emitters_type='electrons',distance=2*u.kpc,name='pwn')
pwn.add_EC_component(['EC_CMB'])
pwn.parameters.R.val=(2*u.pc).to('cm').value
pwn.parameters.B.val=1E-4
pwn.parameters.gamma_break.val=1E8
pwn.parameters.gmax.val=8E9
pwn.parameters.gmin.val=5E5
pwn.parameters.p.val=3.2
pwn.parameters.p_1.val=3.8
pwn.parameters.N.val=1E-8
pwn.show_model()

--------------------------------------------------------------------------------
model description:
--------------------------------------------------------------------------------
type: GalacticUnbeamed
name: pwn
geometry: spherical

electrons distribution:
 type: bkn
 gamma energy grid size:  201
 gmin grid : 2.000000e+00
 gmax grid : 1.000000e+06
 normalization:  True
 log-values:  False
 ratio of cold protons to relativistic electrons: 1.000000e+00

radiative fields:
 seed photons grid size:  100
 IC emission grid size:  100
 source emissivity lower bound :  1.000000e-120
 spectral components:
   name:Sum, state: on
   name:Sum, hidden: False
   name:Sync, state: self-abs
   name:Sync, hidden: False
   name:SSC, state: on
   name:SSC, hidden: False
   name:EC_CMB, state: on
   name:EC_CMB, hidden: False
external fields transformation method: disk

SED info:
 nu grid size jetkernel: 1000
 nu size: 500
 nu mix (Hz): 1.000000e+06
 nu max (Hz): 1.000000e+30

flux plot lower bound   :  1.000000e-30

--------------------------------------------------------------------------------

WARNING: AstropyDeprecationWarning: 'classic' backend for show_in_notebook() is deprecated as of 6.1. Instead, use the supported backend 'ipydatagrid'. [astropy.table.table]

Table length=10

model name	name	par type	units	val	phys. bound. min	phys. bound. max	log	frozen
pwn	R	region_size	cm	6.171355e+18	1.000000e+03	1.000000e+30	False	False
pwn	B	magnetic_field	gauss	1.000000e-04	0.000000e+00	--	False	False
pwn	NH_cold_to_rel_e	cold_p_to_rel_e_ratio		1.000000e+00	0.000000e+00	--	False	True
pwn	gmin	low-energy-cut-off	lorentz-factor*	5.000000e+05	1.000000e+00	1.000000e+09	False	False
pwn	gmax	high-energy-cut-off	lorentz-factor*	8.000000e+09	1.000000e+00	1.000000e+15	False	False
pwn	N	emitters_density	1 / cm3	1.000000e-08	0.000000e+00	--	False	False
pwn	gamma_break	turn-over-energy	lorentz-factor*	1.000000e+08	1.000000e+00	1.000000e+09	False	False
pwn	p	LE_spectral_slope		3.200000e+00	-1.000000e+01	1.000000e+01	False	False
pwn	p_1	HE_spectral_slope		3.800000e+00	-1.000000e+01	1.000000e+01	False	False
pwn	DL_cm	distance	cm	6.171355e+21	0.000000e+00	3.085678e+24	False	False

--------------------------------------------------------------------------------

pwn.eval()
p=pwn.plot_model(frame='src')

pwn.save_model('pwn.pkl')

GalacticUnbeamed.load_model("pwn.pkl")

--------------------------------------------------------------------------------
model description:
--------------------------------------------------------------------------------
type: GalacticUnbeamed
name: pwn
geometry: spherical

electrons distribution:
 type: bkn
 gamma energy grid size:  201
 gmin grid : 5.000000e+05
 gmax grid : 8.000000e+09
 normalization:  True
 log-values:  False
 ratio of cold protons to relativistic electrons: 1.000000e+00

radiative fields:
 seed photons grid size:  100
 IC emission grid size:  100
 source emissivity lower bound :  1.000000e-120
 spectral components:
   name:Sum, state: on
   name:Sum, hidden: False
   name:Sync, state: self-abs
   name:Sync, hidden: False
   name:SSC, state: on
   name:SSC, hidden: False
   name:EC_CMB, state: on
   name:EC_CMB, hidden: False
external fields transformation method: disk

SED info:
 nu grid size jetkernel: 1000
 nu size: 500
 nu mix (Hz): 1.000000e+06
 nu max (Hz): 1.000000e+30

flux plot lower bound   :  1.000000e-30

--------------------------------------------------------------------------------

WARNING: AstropyDeprecationWarning: 'classic' backend for show_in_notebook() is deprecated as of 6.1. Instead, use the supported backend 'ipydatagrid'. [astropy.table.table]

Table length=12

model name	name	par type	units	val	phys. bound. min	phys. bound. max	log	frozen
pwn	gmin	low-energy-cut-off	lorentz-factor*	5.000000e+05	1.000000e+00	1.000000e+09	False	False
pwn	gmax	high-energy-cut-off	lorentz-factor*	8.000000e+09	1.000000e+00	1.000000e+15	False	False
pwn	N	emitters_density	1 / cm3	1.000000e-08	0.000000e+00	--	False	False
pwn	gamma_break	turn-over-energy	lorentz-factor*	1.000000e+08	1.000000e+00	1.000000e+09	False	False
pwn	p	LE_spectral_slope		3.200000e+00	-1.000000e+01	1.000000e+01	False	False
pwn	p_1	HE_spectral_slope		3.800000e+00	-1.000000e+01	1.000000e+01	False	False
pwn	R	region_size	cm	6.171355e+18	1.000000e+03	1.000000e+30	False	False
pwn	R_H	region_position	cm	1.000000e+17	0.000000e+00	--	False	True
pwn	B	magnetic_field	gauss	1.000000e-04	0.000000e+00	--	False	False
pwn	NH_cold_to_rel_e	cold_p_to_rel_e_ratio		1.000000e+00	0.000000e+00	--	False	True
pwn	beam_obj	beaming		1.000000e+00	1.000000e-04	--	False	False
pwn	DL_cm	distance	cm	6.171355e+21	0.000000e+00	3.085678e+24	False	False

--------------------------------------------------------------------------------

None

pwn.energetic_report()

WARNING: AstropyDeprecationWarning: 'classic' backend for show_in_notebook() is deprecated as of 6.1. Instead, use the supported backend 'ipydatagrid'. [astropy.table.table]

Table length=19

name	type	units	val
BulkLorentzFactor	jet-bulk-factor		1.000000e+00
U_e	Energy dens. rest. frame	erg / cm3	7.480812e-09
U_p_cold	Energy dens. rest. frame	erg / cm3	1.503276e-11
U_B	Energy dens. rest. frame	erg / cm3	3.978874e-10
U_Synch	Energy dens. rest. frame	erg / cm3	1.512730e-11
U_Disk	Energy dens. rest. frame	erg / cm3	0.000000e+00
U_BLR	Energy dens. rest. frame	erg / cm3	0.000000e+00
U_DT	Energy dens. rest. frame	erg / cm3	0.000000e+00
U_CMB	Energy dens. rest. frame	erg / cm3	4.185297e-13
U_Star	Energy dens. rest. frame	erg / cm3	0.000000e+00
U_seed_tot	Energy dens. rest. frame	erg / cm3	1.512730e-11
L_Sync	Lum. rest. frame.	erg / s	2.170466e+38
L_SSC	Lum. rest. frame.	erg / s	2.312959e+35
L_EC_Disk	Lum. rest. frame.	erg / s	0.000000e+00
L_EC_BLR	Lum. rest. frame.	erg / s	0.000000e+00
L_EC_DT	Lum. rest. frame.	erg / s	0.000000e+00
L_EC_CMB	Lum. rest. frame.	erg / s	2.139695e+35
L_EC_Star	Lum. rest. frame.	erg / s	0.000000e+00
NH_cold_to_rel_e	cold_p_to_rel_e_ratio		1.000000e+00

Galactic Unbeamed hadronic pp example#

This is an example for a galactic unbeamed hadronic (pp) object

gal_hadronic=GalacticUnbeamed(emitters_distribution='plc',verbose=False,emitters_type='protons',distance=2*u.kpc,name='gal_hadronic_unbeamed')
gal_hadronic.parameters.R.val=1E18
gal_hadronic.parameters.N.val=1000
gal_hadronic.parameters.B.val=1E-3

gal_hadronic.show_model()

--------------------------------------------------------------------------------
model description:
--------------------------------------------------------------------------------
type: GalacticUnbeamed
name: gal_hadronic_unbeamed
geometry: spherical

protons distribution:
 type: plc
 gamma energy grid size:  201
 gmin grid : 2.000000e+00
 gmax grid : 1.000000e+06
 normalization:  True
 log-values:  False

radiative fields:
 seed photons grid size:  100
 IC emission grid size:  100
 source emissivity lower bound :  1.000000e-120
 spectral components:
   name:Sum, state: on
   name:Sum, hidden: False
   name:Sync, state: self-abs
   name:Sync, hidden: False
   name:SSC, state: on
   name:SSC, hidden: False
   name:PP_gamma, state: on
   name:PP_gamma, hidden: False
   name:PP_neutrino_tot, state: on
   name:PP_neutrino_tot, hidden: False
   name:PP_neutrino_mu, state: on
   name:PP_neutrino_mu, hidden: False
   name:PP_neutrino_e, state: on
   name:PP_neutrino_e, hidden: False
   name:Bremss_ep, state: on
   name:Bremss_ep, hidden: False
external fields transformation method: disk

SED info:
 nu grid size jetkernel: 1000
 nu size: 500
 nu mix (Hz): 1.000000e+06
 nu max (Hz): 1.000000e+30

flux plot lower bound   :  1.000000e-30

--------------------------------------------------------------------------------

WARNING: AstropyDeprecationWarning: 'classic' backend for show_in_notebook() is deprecated as of 6.1. Instead, use the supported backend 'ipydatagrid'. [astropy.table.table]

Table length=10

model name	name	par type	units	val	phys. bound. min	phys. bound. max	log	frozen
gal_hadronic_unbeamed	R	region_size	cm	1.000000e+18	1.000000e+03	1.000000e+30	False	False
gal_hadronic_unbeamed	B	magnetic_field	gauss	1.000000e-03	0.000000e+00	--	False	False
gal_hadronic_unbeamed	T_esc_e_secondaries	escape_time	R / c	1.000000e+00	1.000000e+00	--	False	False
gal_hadronic_unbeamed	gmin	low-energy-cut-off	lorentz-factor*	2.000000e+00	1.000000e+00	1.000000e+09	False	False
gal_hadronic_unbeamed	gmax	high-energy-cut-off	lorentz-factor*	1.000000e+06	1.000000e+00	1.000000e+15	False	False
gal_hadronic_unbeamed	N	emitters_density	1 / cm3	1.000000e+03	0.000000e+00	--	False	False
gal_hadronic_unbeamed	NH_pp	target_density	1 / cm3	1.000000e+00	0.000000e+00	--	False	False
gal_hadronic_unbeamed	gamma_cut	turn-over-energy	lorentz-factor*	1.000000e+04	1.000000e+00	1.000000e+09	False	False
gal_hadronic_unbeamed	p	LE_spectral_slope		2.000000e+00	-1.000000e+01	1.000000e+01	False	False
gal_hadronic_unbeamed	DL_cm	distance	cm	6.171355e+21	0.000000e+00	3.085678e+24	False	False

--------------------------------------------------------------------------------

gal_hadronic.eval()
p=gal_hadronic.plot_model(frame='src')
p.setlim(y_min=1E33)

gal_hadronic.energetic_report()

WARNING: AstropyDeprecationWarning: 'classic' backend for show_in_notebook() is deprecated as of 6.1. Instead, use the supported backend 'ipydatagrid'. [astropy.table.table]

Table length=20

name	type	units	val
BulkLorentzFactor	jet-bulk-factor		1.000000e+00
U_e	Energy dens. rest. frame	erg / cm3	3.408013e-08
U_B	Energy dens. rest. frame	erg / cm3	3.978874e-08
U_p	Energy dens. rest. frame	erg / cm3	2.386056e+01
U_p_target	Energy dens. rest. frame	erg / cm3	1.503276e-03
U_Synch	Energy dens. rest. frame	erg / cm3	1.108748e-10
U_Disk	Energy dens. rest. frame	erg / cm3	0.000000e+00
U_BLR	Energy dens. rest. frame	erg / cm3	0.000000e+00
U_DT	Energy dens. rest. frame	erg / cm3	0.000000e+00
U_CMB	Energy dens. rest. frame	erg / cm3	0.000000e+00
U_Star	Energy dens. rest. frame	erg / cm3	0.000000e+00
U_seed_tot	Energy dens. rest. frame	erg / cm3	1.108748e-10
L_Sync	Lum. rest. frame.	erg / s	4.176989e+37
L_SSC	Lum. rest. frame.	erg / s	2.253964e+34
L_EC_Disk	Lum. rest. frame.	erg / s	0.000000e+00
L_EC_BLR	Lum. rest. frame.	erg / s	0.000000e+00
L_EC_DT	Lum. rest. frame.	erg / s	0.000000e+00
L_EC_CMB	Lum. rest. frame.	erg / s	0.000000e+00
L_EC_Star	Lum. rest. frame.	erg / s	0.000000e+00
L_pp_gamma	Lum. rest. frame.	erg / s	1.220405e+40

Galactic Beamed hadronic pp example#

from jetset.jet_model import GalacticUnbeamed,GalacticBeamed
from astropy import constants as const
from astropy import units as u
import matplotlib.pyplot as plt
import numpy as np

gal_hadronic=GalacticBeamed(emitters_distribution='plc',verbose=False,emitters_type='protons',distance=2*u.kpc,name='gal_hadronic_beamed')
gal_hadronic.parameters.R.val=1E18
gal_hadronic.parameters.N.val=1000
gal_hadronic.parameters.B.val=1E-3

gal_hadronic.show_model()

--------------------------------------------------------------------------------
model description:
--------------------------------------------------------------------------------
type: GalacticBeamed
name: gal_hadronic_beamed
geometry: spherical

protons distribution:
 type: plc
 gamma energy grid size:  201
 gmin grid : 2.000000e+00
 gmax grid : 1.000000e+06
 normalization:  True
 log-values:  False

radiative fields:
 seed photons grid size:  100
 IC emission grid size:  100
 source emissivity lower bound :  1.000000e-120
 spectral components:
   name:Sum, state: on
   name:Sum, hidden: False
   name:Sync, state: self-abs
   name:Sync, hidden: False
   name:SSC, state: on
   name:SSC, hidden: False
   name:PP_gamma, state: on
   name:PP_gamma, hidden: False
   name:PP_neutrino_tot, state: on
   name:PP_neutrino_tot, hidden: False
   name:PP_neutrino_mu, state: on
   name:PP_neutrino_mu, hidden: False
   name:PP_neutrino_e, state: on
   name:PP_neutrino_e, hidden: False
   name:Bremss_ep, state: on
   name:Bremss_ep, hidden: False
external fields transformation method: blob

SED info:
 nu grid size jetkernel: 1000
 nu size: 500
 nu mix (Hz): 1.000000e+06
 nu max (Hz): 1.000000e+30

flux plot lower bound   :  1.000000e-30

--------------------------------------------------------------------------------

WARNING: AstropyDeprecationWarning: 'classic' backend for show_in_notebook() is deprecated as of 6.1. Instead, use the supported backend 'ipydatagrid'. [astropy.table.table]

Table length=12

model name	name	par type	units	val	phys. bound. min	phys. bound. max	log	frozen
gal_hadronic_beamed	R	region_size	cm	1.000000e+18	1.000000e+03	1.000000e+30	False	False
gal_hadronic_beamed	R_H	region_position	cm	1.000000e+17	0.000000e+00	--	False	True
gal_hadronic_beamed	B	magnetic_field	gauss	1.000000e-03	0.000000e+00	--	False	False
gal_hadronic_beamed	T_esc_e_secondaries	escape_time	R / c	1.000000e+00	1.000000e+00	--	False	False
gal_hadronic_beamed	beam_obj	beaming		1.000000e+01	1.000000e-04	--	False	False
gal_hadronic_beamed	gmin	low-energy-cut-off	lorentz-factor*	2.000000e+00	1.000000e+00	1.000000e+09	False	False
gal_hadronic_beamed	gmax	high-energy-cut-off	lorentz-factor*	1.000000e+06	1.000000e+00	1.000000e+15	False	False
gal_hadronic_beamed	N	emitters_density	1 / cm3	1.000000e+03	0.000000e+00	--	False	False
gal_hadronic_beamed	NH_pp	target_density	1 / cm3	1.000000e+00	0.000000e+00	--	False	False
gal_hadronic_beamed	gamma_cut	turn-over-energy	lorentz-factor*	1.000000e+04	1.000000e+00	1.000000e+09	False	False
gal_hadronic_beamed	p	LE_spectral_slope		2.000000e+00	-1.000000e+01	1.000000e+01	False	False
gal_hadronic_beamed	DL_cm	distance	cm	6.171355e+21	0.000000e+00	3.085678e+24	False	False

--------------------------------------------------------------------------------

gal_hadronic.eval()
p=gal_hadronic.plot_model(frame='src')
p.setlim(y_min=1E36)

gal_hadronic.energetic_report()

WARNING: AstropyDeprecationWarning: 'classic' backend for show_in_notebook() is deprecated as of 6.1. Instead, use the supported backend 'ipydatagrid'. [astropy.table.table]

Table length=40

name	type	units	val
BulkLorentzFactor	jet-bulk-factor		1.000000e+01
U_e	Energy dens. blob rest. frame	erg / cm3	3.408013e-08
U_B	Energy dens. blob rest. frame	erg / cm3	3.978874e-08
U_p	Energy dens. blob rest. frame	erg / cm3	2.386056e+01
U_p_target	Energy dens. blob rest. frame	erg / cm3	1.503276e-03
U_Synch	Energy dens. blob rest. frame	erg / cm3	1.108748e-10
U_Synch_DRF	Energy dens. disk rest. frame	erg / cm3	1.108748e-06
U_Disk	Energy dens. blob rest. frame	erg / cm3	0.000000e+00
U_BLR	Energy dens. blob rest. frame	erg / cm3	0.000000e+00
U_DT	Energy dens. blob rest. frame	erg / cm3	0.000000e+00
U_CMB	Energy dens. blob rest. frame	erg / cm3	0.000000e+00
U_Star	Energy dens. blob rest. frame	erg / cm3	0.000000e+00
U_Disk_DRF	Energy dens. disk rest. frame	erg / cm3	0.000000e+00
U_BLR_DRF	Energy dens. disk rest. frame	erg / cm3	0.000000e+00
U_DT_DRF	Energy dens. disk rest. frame	erg / cm3	0.000000e+00
U_CMB_DRF	Energy dens. disk rest. frame	erg / cm3	0.000000e+00
U_Star_DRF	Energy dens. disk rest. frame	erg / cm3	0.000000e+00
U_seed_tot	Energy dens. blob rest. frame	erg / cm3	1.108748e-10
L_Sync_rf	Lum. blob rest. frame.	erg / s	4.176989e+37
L_SSC_rf	Lum. blob rest. frame.	erg / s	2.253964e+34
L_EC_Disk_rf	Lum. blob rest. frame.	erg / s	0.000000e+00
L_EC_BLR_rf	Lum. blob rest. frame.	erg / s	0.000000e+00
L_EC_DT_rf	Lum. blob rest. frame.	erg / s	0.000000e+00
L_EC_CMB_rf	Lum. blob rest. frame.	erg / s	0.000000e+00
L_EC_Star_rf	Lum. blob rest. frame.	erg / s	0.000000e+00
L_pp_gamma_rf	Lum. blob rest. frame.	erg / s	1.220405e+40
jet_L_Sync	jet Lum.	erg / s	1.039013e+39
jet_L_SSC	jet Lum.	erg / s	5.606665e+35
jet_L_EC_Disk	jet Lum.	erg / s	0.000000e+00
jet_L_EC_BLR	jet Lum.	erg / s	0.000000e+00
jet_L_EC_Star	jet Lum.	erg / s	0.000000e+00
jet_L_EC_DT	jet Lum.	erg / s	0.000000e+00
jet_L_EC_CMB	jet Lum.	erg / s	0.000000e+00
jet_L_pp_gamma	jet Lum.	erg / s	3.035720e+41
jet_L_rad	jet Lum.	erg / s	3.046115e+41
jet_L_kin	jet Lum.	erg / s	2.235985e+50
jet_L_tot	jet Lum.	erg / s	2.235985e+50
jet_L_e	jet Lum.	erg / s	3.193665e+41
jet_L_B	jet Lum.	erg / s	3.728622e+41
jet_L_p	jet Lum.	erg / s	2.235985e+50

gal_hadronic.save_model('gal_hadronic_beamed.pkl')

gal_hadronic_new=GalacticBeamed.load_model('gal_hadronic_beamed.pkl')

Galactic object models

Contents

Galactic object models#

Galactic Unbeamed leptonic example#

Galactic Unbeamed hadronic pp example#

Galactic Beamed hadronic pp example#