Working with the PHOENIX models¶

The PHOENIX models are under active development. Check back soon!

[1]:

from gollum.phoenix import PHOENIXSpectrum
import numpy as np
%config InlineBackend.figure_format='retina'

[2]:

spec = PHOENIXSpectrum(teff=2800, logg=5.0)

[3]:

spec.wavelength # "Angstroms"

[3]:

$[8038.01,~8038.02,~8038.03,~\dots,~12848.94,~12848.96,~12848.98] \; \mathrm{\mathring{A}}$

[4]:

spec.flux.unit # "egs / s / cm^2 / cm"

[4]:

$\mathrm{\frac{erg}{s\,cm^{3}}}$

We can normalize the spectrum, which divides the spectrum by the median value, rendering the flux units dimensionless

[5]:

spectrum = spec.normalize()
np.median(spectrum.flux)

[5]:

$1 \; \mathrm{}$

[6]:

ax = spectrum.plot()
ax.set_xlim(8160, 8220)
ax.set_ylim(0, 1.5)

[6]:

(0.0, 1.5)

Neat! There are many great uses for the PHOENIX models. We will be adding support for other model grids in the future.