8.4. Observation constraints on masses and radii (NICER observatory)#
In this tutorial, you will learn how to extract the observational constraints on the masses adn radii obtained by NICER observatory.
Import the libraries that will be employed in this tutorial.
# Import numpy
import numpy as np
# Import matplotlib
import matplotlib.pyplot as plt
# Import nucleardatapy package
import nucleardatapy as nuda
You can simply print out the properties of the nuda’s function that we will use:
# Explore the nucleardatapy module to find the correct attribute
print(dir(nuda.astro.setupMR))
['__class__', '__delattr__', '__dict__', '__dir__', '__doc__', '__eq__', '__format__', '__ge__', '__getattribute__', '__getstate__', '__gt__', '__hash__', '__init__', '__init_subclass__', '__le__', '__lt__', '__module__', '__ne__', '__new__', '__reduce__', '__reduce_ex__', '__repr__', '__setattr__', '__sizeof__', '__str__', '__subclasshook__', '__weakref__', 'print_latex', 'print_output']
Complete list of available sources:
sources, sources_lower = nuda.astro.mr_sources( )
print('sources:',sources)
sources: ['J0030+0451', 'J0740+6620', 'J0437-4715']
source = 'J0030+0451'
obss = nuda.astro.mr_obss( source = source )
print(f'source: {source}, obss: {obss}')
source: J0030+0451, obss: [1, 2, 3, 4]
for source in sources:
obss = nuda.astro.mr_obss( source = source )
print(f'source: {source}, obss: {obss}')
for obs in obss:
mr = nuda.astro.setupMR( source = source, obs = obs )
mr.print_output()
mrav = nuda.astro.setupMRAverage( source = source, obss = obss )
mrav.print_output()
source: J0030+0451, obss: [1, 2, 3, 4]
- Print output:
source: J0030+0451
obs: 1
mass: 1.44 in Mo
sigma(mass): 0.15 0.14 in Mo
rad: 13.02 in km
sigma(mass): 1.24 1.06 in km
compactness: 0.163
sigma(comp): 0.008 0.009
latexCite:
ref: M.C. Miller, F.K. Lamb, A.J. Dittmann, aet al., ApJL 887, L24 (2019).
label: J0030 Miller 2019
note: write notes about this observation.
- Print output:
source: J0030+0451
obs: 2
mass: 1.34 in Mo
sigma(mass): 0.15 0.16 in Mo
rad: 12.71 in km
sigma(mass): 1.14 1.19 in km
compactness: 0.156
sigma(comp): 0.008 0.01
latexCite:
ref: T.E. Riley, A.L. Watts, S. Bogdanov, P.S. Ray, et al., ApJ 887, L21 (2019).
label: J0030 Riley 2019
note: write notes about this observation.
- Print output:
source: J0030+0451
obs: 3
mass: 1.7 in Mo
sigma(mass): 0.18 0.19 in Mo
rad: 14.44 in km
sigma(mass): 0.88 1.05 in km
compactness: 0.179
sigma(comp): 0.011 0.022
latexCite:
ref: S. Vinciguerra, T. Salmi, A.L. Watts, D. Choudhury, et al., ApJ 961, 62 (2024).
label: J0030 Vinciguerra 2024a
note: write notes about this observation.
- Print output:
source: J0030+0451
obs: 4
mass: 1.4 in Mo
sigma(mass): 0.13 0.12 in Mo
rad: 11.71 in km
sigma(mass): 0.88 0.83 in km
compactness: 0.1773
sigma(comp): 0.0056 0.0074
latexCite:
ref: S. Vinciguerra, T. Salmi, A.L. Watts, D. Choudhury, et al., ApJ 961, 62 (2024).
label: J0030 Vinciguerra 2024b
note: write notes about this observation.
cmax: 0.212
cmin: 0.11299999999999999
- No output for source J0030+0451 (average). To get output, write 'verb_output = True' in env.py.
source: J0740+6620, obss: [1, 2, 3]
- Print output:
source: J0740+6620
obs: 1
mass: 2.08 in Mo
sigma(mass): 0.07 0.07 in Mo
rad: 13.7 in km
sigma(mass): 2.6 1.5 in km
compactness: 0.0
sigma(comp): 1.0 1.0
latexCite: MCMiller:2021
ref: M.C. Miller, F.K. Lamb, A.J. Dittmann, S. Bogdanov, et al., ApJL 918, L28 (2021).
label: J0740 Miller 2021
note: write notes about this observation.
- Print output:
source: J0740+6620
obs: 2
mass: 2.072 in Mo
sigma(mass): 0.067 0.066 in Mo
rad: 12.39 in km
sigma(mass): 1.3 0.98 in km
compactness: 0.0
sigma(comp): 1.0 1.0
latexCite: TERiley:2021
ref: T.E. Riley, A.L. Watts, P.S. Ray, S. Bogdanov, et al., ApJL 918, L27 (2021).
label: J0740 Riley 2021
note: write notes about this observation.
- Print output:
source: J0740+6620
obs: 3
mass: 2.073 in Mo
sigma(mass): 0.069 0.069 in Mo
rad: 12.49 in km
sigma(mass): 1.28 0.88 in km
compactness: 0.0
sigma(comp): 1.0 1.0
latexCite: TSalmi:2024
ref: T. Salmi, D. Choudhury, Y. Kini, T.E. Riley et al., ApJ 974, 294 (2024).
label: J0740 Salmi 2024
note: write notes about this observation.
cmax: 3.0
cmin: -3.0
- No output for source J0740+6620 (average). To get output, write 'verb_output = True' in env.py.
source: J0437-4715, obss: [1]
- Print output:
source: J0437-4715
obs: 1
mass: 1.418 in Mo
sigma(mass): 0.037 0.037 in Mo
rad: 11.36 in km
sigma(mass): 0.945 0.629 in km
compactness: 0.1847
sigma(comp): 0.0097 0.0143
latexCite: DChoudhury:2024
ref: D. Choudhury, T. Salmi, S. Vinciguerra, T.E. Riley, et al., ApJL 971, L20 (2024).
label: J0437 Choudhury 2024
note: write notes about this observation.
cmax: 0.2138
cmin: 0.1418
- No output for source J0437-4715 (average). To get output, write 'verb_output = True' in env.py.
Figure:
sources = [ 'J0030+0451', 'J0740+6620', 'J0437-4715' ]
sources_av = [ 'J0030+0451', 'J0740+6620' ]
nuda.fig.astro_setupMR_fig( None, sources, sources_av )
Plot name: None
source: J0030+0451, obss: [1, 2, 3, 4]
- table: J0030+0451 & 1 & ${13.02}^{1.24}_{-1.06}$ & $1.440^{0.15}_{-0.14}$ & $0.163^{0.008}_{-0.009}$ & \cite{} \\\\
- table: J0030+0451 & 2 & ${12.71}^{1.14}_{-1.19}$ & $1.340^{0.15}_{-0.16}$ & $0.156^{0.008}_{-0.01}$ & \cite{} \\\\
- table: J0030+0451 & 3 & ${14.44}^{0.88}_{-1.05}$ & $1.700^{0.18}_{-0.19}$ & $0.179^{0.011}_{-0.022}$ & \cite{} \\\\
- table: J0030+0451 & 4 & ${11.71}^{0.88}_{-0.83}$ & $1.400^{0.13}_{-0.12}$ & $0.1773^{0.0056}_{-0.0074}$ & \cite{} \\\\
source: J0740+6620, obss: [1, 2, 3]
- table: J0740+6620 & 1 & ${13.70}^{2.6}_{-1.5}$ & $2.080^{0.07}_{-0.07}$ & $0.0^{1.0}_{-1.0}$ & \cite{MCMiller:2021} \\\\
- table: J0740+6620 & 2 & ${12.39}^{1.3}_{-0.98}$ & $2.072^{0.067}_{-0.066}$ & $0.0^{1.0}_{-1.0}$ & \cite{TERiley:2021} \\\\
- table: J0740+6620 & 3 & ${12.49}^{1.28}_{-0.88}$ & $2.073^{0.069}_{-0.069}$ & $0.0^{1.0}_{-1.0}$ & \cite{TSalmi:2024} \\\\
source: J0437-4715, obss: [1]
- table: J0437-4715 & 1 & ${11.36}^{0.945}_{-0.629}$ & $1.418^{0.037}_{-0.037}$ & $0.1847^{0.0097}_{-0.0143}$ & \cite{DChoudhury:2024} \\\\
cmax: 0.187
cmin: 0.126
- table: J0030+0451 & av & $12.89\pm1.166$ & $1.390\pm0.158$ & $0.159\pm0.009$ & \\\\
cmax: 3.0
cmin: -3.0
- table: J0740+6620 & av & $13.10\pm1.663$ & $2.075\pm0.068$ & $0.0\pm0.987$ & \\\\
source: J0030+0451
source: J0740+6620
source: J0437-4715