The LineCollection Object¶
Emission lines are bright features in a spectrum caused by specific atomic transitions. In astronomy, these lines are often used to identify the presence of certain elements in celestial objects, such as galaxies or nebulae. The luminosity of these lines can vary significantly, and they can be influenced by various factors such as redshift, dust extinction, and the physical conditions of the emitting gas.
In Synthesizer emission lines are encapsulated by the LineCollection object. This object holds an arbitrary collection of emission lines, with each luminosity, continuum, line ID, and wavelength stored in arrays. However, as well as providing an interface to these arrays, the LineCollection can also be treated as a dictionary.
In this section, we will demonstrate how to create a LineCollection object, how to access its properties, and how to manipulate it.
Creating a LineCollection¶
Before demonstrating anything, we need a LineCollection to demonstrate with. A LineCollection can be created from scratch, taking user specified luminosities, continuum, line IDs, and wavelengths.
[2]:
from synthesizer import LineCollection
lines = LineCollection(
"ArbitraryLine",
lam=912 * angstrom,
lum=1e28 * erg / s,
cont=1e28 * erg / s / Hz,
)
However, in practice, it will be far more common to use LineCollection objects produced by or derived from other Synthesizer objects (e.g. from combining a Galaxy and an EmissionModel).
For the rest of this section, we will take lines directly from a Grid object (more information here). For more details on line generation see the grid lines and galaxy lines notebooks for thorough demonstrations of how to generate lines from each.
[3]:
from synthesizer.grid import Grid
grid = Grid("test_grid")
# Extract the lines for a specific point in the grid
lines = grid.get_lines(
grid.get_grid_point(log10ages=7.0 * Myr, metallicity=0.01)
)
Printing a summary of the LineCollection¶
If we want a summary of what a line contains we can simply print it. This will give us a table showing all the different attributes attached to a LineCollection object and their values (or a summary of their values for large arrays).
[4]:
print(lines)
print(lines.line_ids)
+-------------------------------------------------------------------------------------------------+
| LINECOLLECTION |
+---------------------------+---------------------------------------------------------------------+
| Attribute | Value |
+---------------------------+---------------------------------------------------------------------+
| nlines | 215 |
+---------------------------+---------------------------------------------------------------------+
| ndim | 1 |
+---------------------------+---------------------------------------------------------------------+
| nlam | 215 |
+---------------------------+---------------------------------------------------------------------+
| shape | (215,) |
+---------------------------+---------------------------------------------------------------------+
| available_ratios | [BalmerDecrement, N2, S2, O1, R2, R3, R23, |
| | O32, Ne3O2] |
+---------------------------+---------------------------------------------------------------------+
| available_diagrams | [OHNO, BPT-NII, VO78-SII, VO78-OI] |
+---------------------------+---------------------------------------------------------------------+
| elements | [He, O, O, He, Si, Si, He, H, O, N, N, Si, |
| | Si, Si, O, O, O, C, C, C, Si, O, O, Si, O, |
| | O, N, Si, Si, C, C, Ne, He, O, C, O, O, Al, |
| | N, Mg, Si, Si, C, C, C, Si, O, C, C, Fe, Fe, |
| | Fe, Fe, C, Fe, Fe, Fe, Fe, Fe, Fe, Mg, Fe, |
| | Fe, Mg, Mg, Mg, Fe, Fe, Ar, Fe, He, Ne, Ne, |
| | Fe, Fe, O, O, H, H, Fe, H, H, H, Ne, He, H, |
| | Fe, Ne, H, S, S, H, Fe, Fe, Fe, Fe, H, Fe, |
| | Fe, Fe, Fe, O, Fe, Fe, Fe, Fe, Fe, He, Fe, |
| | H, Fe, Fe, Fe, Fe, Fe, O, Fe, Fe, Fe, O, Fe, |
| | Fe, Fe, Fe, Fe, Fe, Fe, Fe, Fe, Fe, Fe, Fe, |
| | Fe, Fe, Fe, Fe, Fe, Fe, Fe, Fe, He, He, Fe, |
| | O, O, Fe, N, H, N, Ni, He, S, S, Ar, Fe, Fe, |
| | Ca, Ca, Ni, Fe, Ni, Fe, Ar, O, O, O, Cl, Fe, |
| | Fe, Fe, S, Fe, H, Fe, Fe, Fe, S, H, H, S, |
| | S, S, He, He, H, O, O, O, O, O, Ni, Fe, Fe, |
| | H, Fe, Fe, Fe, Fe, Fe, Fe, Fe, Fe, Fe, Fe, |
| | Fe, Fe, Fe, Fe, H, Fe, Ni, Fe, Si, H, Si,] |
+---------------------------+---------------------------------------------------------------------+
| line_ids (215,) | [He 2 1025.27A, O 6 1031.91A, O 6 1037.61A, ...] |
+---------------------------+---------------------------------------------------------------------+
| lam (215,) | 1.03e+03 Å -> 2.48e+04 Å (Mean: 5.98e+03 Å) |
+---------------------------+---------------------------------------------------------------------+
| luminosity (215,) | 0.00e+00 erg/s -> 3.34e+34 erg/s (Mean: 3.52e+32 erg/s) |
+---------------------------+---------------------------------------------------------------------+
| continuum (215,) | 2.32e+17 erg/(Hz*s) -> 6.24e+19 erg/(Hz*s) (Mean: 1.62e+19 erg) |
+---------------------------+---------------------------------------------------------------------+
| cont (215,) | 2.32e+17 erg/(Hz*s) -> 6.24e+19 erg/(Hz*s) (Mean: 1.62e+19 erg) |
+---------------------------+---------------------------------------------------------------------+
| continuum_llam (215,) | 7.40e+28 erg/(s*Å) -> 2.52e+31 erg/(s*Å) (Mean: 5.76e+30 erg/(s*Å)) |
+---------------------------+---------------------------------------------------------------------+
| energy (215,) | 5.00e-01 eV -> 1.21e+01 eV (Mean: 3.76e+00 eV) |
+---------------------------+---------------------------------------------------------------------+
| equivalent_width (215,) | 0.00e+00 Å -> 7.09e+04 Å (Mean: 5.37e+02 Å) |
+---------------------------+---------------------------------------------------------------------+
| lum (215,) | 0.00e+00 erg/s -> 3.34e+34 erg/s (Mean: 3.52e+32 erg/s) |
+---------------------------+---------------------------------------------------------------------+
| nu (215,) | 1.21e+14 Hz -> 2.92e+15 Hz (Mean: 9.09e+14 Hz) |
+---------------------------+---------------------------------------------------------------------+
| vacuum_wavelengths (215,) | 1.03e+03 Å -> 2.48e+04 Å (Mean: 5.98e+03 Å) |
+---------------------------+---------------------------------------------------------------------+
['He 2 1025.27A' 'O 6 1031.91A' 'O 6 1037.61A' 'He 2 1084.94A'
'Si 2 1179.59A' 'Si 3 1206.50A' 'He 2 1215.13A' 'H 1 1215.67A'
'O 5 1218.34A' 'N 5 1238.82A' 'N 5 1242.80A' 'Si 2 1260.42A'
'Si 2 1264.74A' 'Si 2 1265.00A' 'O 1 1302.17A' 'O 1 1304.86A'
'O 1 1306.03A' 'C 2 1334.53A' 'C 2 1335.66A' 'C 2 1335.71A'
'Si 4 1393.75A' 'O 4 1399.78A' 'O 4 1401.16A' 'Si 4 1402.77A'
'O 4 1404.81A' 'O 4 1407.38A' 'N 4 1486.50A' 'Si 2 1526.71A'
'Si 2 1533.43A' 'C 4 1548.19A' 'C 4 1550.77A' 'Ne 4 1601.45A'
'He 2 1640.41A' 'O 1 1641.31A' 'C 1 1657.91A' 'O 3 1660.81A'
'O 3 1666.15A' 'Al 2 1670.79A' 'N 3 1749.67A' 'Mg 6 1806.00A'
'Si 3 1882.71A' 'Si 3 1892.03A' 'C 3 1906.68A' 'C 3 1908.73A'
'C 1 1992.01A' 'Si 7 2146.64A' 'O 3 2320.95A' 'C 2 2325.40A'
'C 2 2326.93A' 'Fe 2 2395.63A' 'Fe 2 2399.24A' 'Fe 2 2406.66A'
'Fe 2 2410.52A' 'C 1 2582.90A' 'Fe 2 2598.37A' 'Fe 2 2607.09A'
'Fe 2 2611.87A' 'Fe 2 2613.82A' 'Fe 2 2625.67A' 'Fe 2 2628.29A'
'Mg 7 2628.89A' 'Fe 2 2631.05A' 'Fe 2 2631.32A' 'Mg 5 2782.76A'
'Mg 2 2795.53A' 'Mg 2 2802.71A' 'Fe 4 2829.36A' 'Fe 4 2835.74A'
'Ar 4 2853.66A' 'Fe 4 3094.96A' 'He 1 3187.74A' 'Ne 5 3345.82A'
'Ne 5 3425.88A' 'Fe 7 3586.32A' 'Fe 6 3662.50A' 'O 2 3726.03A'
'O 2 3728.81A' 'H 1 3734.37A' 'H 1 3750.15A' 'Fe 7 3758.92A'
'H 1 3770.63A' 'H 1 3797.90A' 'H 1 3835.38A' 'Ne 3 3868.76A'
'He 1 3888.64A' 'H 1 3889.05A' 'Fe 5 3891.28A' 'Ne 3 3967.47A'
'H 1 3970.07A' 'S 2 4068.60A' 'S 2 4076.35A' 'H 1 4101.73A'
'Fe 2 4243.97A' 'Fe 2 4276.84A' 'Fe 2 4287.39A' 'Fe 2 4319.62A'
'H 1 4340.46A' 'Fe 2 4346.86A' 'Fe 2 4352.79A' 'Fe 2 4358.37A'
'Fe 2 4359.33A' 'O 3 4363.21A' 'Fe 2 4413.78A' 'Fe 2 4416.27A'
'Fe 2 4452.10A' 'Fe 2 4474.90A' 'Fe 3 4658.05A' 'He 2 4685.68A'
'Fe 2 4814.54A' 'H 1 4861.32A' 'Fe 2 4874.50A' 'Fe 2 4889.62A'
'Fe 2 4905.35A' 'Fe 2 4923.92A' 'Fe 2 4947.39A' 'O 3 4958.91A'
'Fe 2 4973.40A' 'Fe 3 4985.87A' 'Fe 2 5005.52A' 'O 3 5006.84A'
'Fe 2 5018.44A' 'Fe 2 5020.25A' 'Fe 2 5049.30A' 'Fe 2 5072.41A'
'Fe 2 5111.64A' 'Fe 2 5158.01A' 'Fe 2 5158.79A' 'Fe 2 5169.03A'
'Fe 6 5176.04A' 'Fe 2 5184.80A' 'Fe 2 5261.63A' 'Fe 3 5270.40A'
'Fe 2 5273.36A' 'Fe 2 5284.10A' 'Fe 2 5333.66A' 'Fe 2 5376.47A'
'Fe 2 5412.67A' 'Fe 2 5433.15A' 'Fe 2 5527.36A' 'Fe 7 5720.71A'
'He 1 5875.61A' 'He 1 5875.64A' 'Fe 7 6086.97A' 'O 1 6300.30A'
'O 1 6363.78A' 'Fe 2 6516.08A' 'N 2 6548.05A' 'H 1 6562.80A'
'N 2 6583.45A' 'Ni 2 6666.80A' 'He 1 6678.15A' 'S 2 6716.44A'
'S 2 6730.82A' 'Ar 3 7135.79A' 'Fe 2 7155.17A' 'Fe 2 7172.00A'
'Ca 2 7291.47A' 'Ca 2 7323.89A' 'Ni 2 7377.83A' 'Fe 2 7388.17A'
'Ni 2 7411.61A' 'Fe 2 7452.56A' 'Ar 3 7751.11A' 'O 1 8446.25A'
'O 1 8446.36A' 'O 1 8446.76A' 'Cl 2 8578.70A' 'Fe 2 8616.95A'
'Fe 2 8891.93A' 'Fe 2 9051.95A' 'S 3 9068.62A' 'Fe 2 9226.63A'
'H 1 9229.02A' 'Fe 2 9267.56A' 'Fe 2 9399.04A' 'Fe 2 9470.94A'
'S 3 9530.62A' 'H 1 9545.97A' 'H 1 1.00494m' 'S 2 1.02867m'
'S 2 1.03205m' 'S 2 1.03364m' 'He 1 1.08291m' 'He 1 1.08303m'
'H 1 1.09381m' 'O 1 1.12863m' 'O 1 1.12864m' 'O 1 1.12869m'
'O 1 1.12870m' 'O 1 1.12873m' 'Ni 2 1.19102m' 'Fe 2 1.25668m'
'Fe 2 1.27877m' 'H 1 1.28181m' 'Fe 2 1.29427m' 'Fe 2 1.32055m'
'Fe 2 1.32777m' 'Fe 2 1.37181m' 'Fe 2 1.53348m' 'Fe 2 1.59948m'
'Fe 2 1.64355m' 'Fe 2 1.66377m' 'Fe 2 1.67688m' 'Fe 2 1.71113m'
'Fe 2 1.74494m' 'Fe 2 1.79711m' 'Fe 2 1.80002m' 'Fe 2 1.80940m'
'H 1 1.87510m' 'Fe 2 1.89541m' 'Ni 2 1.93877m' 'Fe 2 1.95361m'
'Si 6 1.96247m' 'H 1 2.16553m' 'Si 7 2.48071m']
Accessing specific lines¶
A LineCollection behaves somewhat like a dictionary, we can access line data by passing strings when we index the LineCollection. Unlike a dictionary, however, we are not limited to singular keys, we can pass multiple at once in a list or tuple to instead return a subset of lines. Not only that, we can also defined composite lines by passing a string of comma separated lines. We demonstrate each of these methods to axis data below.
Note that in reality the data is actually stored in contiguous arrays for efficiency.
Line IDs follow the same convention as in Cloudy, whereby a line is usually represented as {atomic or molecular notation} {ionisation state} {wavelength}. The ionisation state is in the usual astronomical notation, e.g. H 1 for atomic hydrogen, but different for molecules, with the state denoted by ‘+’ or ‘-’ (e.g. HCO+).
Extracting a single line¶
To extract a single line we simply pass the name of the line we want to extract.
[5]:
print(lines["H 1 1215.67A"])
+-------------------------------------------------+
| LINECOLLECTION |
+--------------------+----------------------------+
| Attribute | Value |
+--------------------+----------------------------+
| nlines | 1 |
+--------------------+----------------------------+
| id | 'H 1 1215.67A' |
+--------------------+----------------------------+
| ndim | 1 |
+--------------------+----------------------------+
| nlam | 1 |
+--------------------+----------------------------+
| shape | (1,) |
+--------------------+----------------------------+
| elements | [H, ] |
+--------------------+----------------------------+
| line_ids | ['H 1 1215.67A'] |
+--------------------+----------------------------+
| lam | [1215.67] Å |
+--------------------+----------------------------+
| luminosity | [3.33710923e+34] erg/s |
+--------------------+----------------------------+
| continuum | [2.3214967e+17] erg/(Hz*s) |
+--------------------+----------------------------+
| cont | [2.3214967e+17] erg/(Hz*s) |
+--------------------+----------------------------+
| continuum_llam | [4.70931104e+29] erg/(s*Å) |
+--------------------+----------------------------+
| energy | [10.19883629] eV |
+--------------------+----------------------------+
| equivalent_width | [70861.94139551] Å |
+--------------------+----------------------------+
| lum | [3.33710923e+34] erg/s |
+--------------------+----------------------------+
| nu | [2.46606775e+15] Hz |
+--------------------+----------------------------+
| vacuum_wavelengths | [1215.67] Å |
+--------------------+----------------------------+
Extracting composite lines¶
To extract a composite line (e.g. a doublet or triplet) we pass a string of comma separated lines. This will return a new LineCollection object with the composite line as the only line in the collection.
[6]:
print(lines["Mg 7 2628.89A, Fe 2 2631.05A"]) # doublet
print(lines["Mg 7 2628.89A, Fe 2 2631.05A, Fe 2 2631.32A"]) # triplet
+-------------------------------------------------------+
| LINECOLLECTION |
+--------------------+----------------------------------+
| Attribute | Value |
+--------------------+----------------------------------+
| nlines | 1 |
+--------------------+----------------------------------+
| id | 'Mg 7 2628.89A, Fe 2 2631.05A' |
+--------------------+----------------------------------+
| ndim | 1 |
+--------------------+----------------------------------+
| nlam | 1 |
+--------------------+----------------------------------+
| shape | (1,) |
+--------------------+----------------------------------+
| elements | [Mg, Fe] |
+--------------------+----------------------------------+
| line_ids | ['Mg 7 2628.89A, Fe 2 2631.05A'] |
+--------------------+----------------------------------+
| lam | [2629.97] Å |
+--------------------+----------------------------------+
| luminosity | [6.7217364e+30] erg/s |
+--------------------+----------------------------------+
| continuum | [4.70482496e+19] erg/(Hz*s) |
+--------------------+----------------------------------+
| cont | [4.70482496e+19] erg/(Hz*s) |
+--------------------+----------------------------------+
| continuum_llam | [2.03921297e+31] erg/(s*Å) |
+--------------------+----------------------------------+
| energy | [4.71428165] eV |
+--------------------+----------------------------------+
| equivalent_width | [0.32962405] Å |
+--------------------+----------------------------------+
| lum | [6.7217364e+30] erg/s |
+--------------------+----------------------------------+
| nu | [1.13990828e+15] Hz |
+--------------------+----------------------------------+
| vacuum_wavelengths | [2630.75425269] Å |
+--------------------+----------------------------------+
+----------------------------------------------------------------------+
| LINECOLLECTION |
+--------------------+-------------------------------------------------+
| Attribute | Value |
+--------------------+-------------------------------------------------+
| nlines | 1 |
+--------------------+-------------------------------------------------+
| id | 'Mg 7 2628.89A, Fe 2 2631.05A, Fe 2 2631.32A' |
+--------------------+-------------------------------------------------+
| ndim | 1 |
+--------------------+-------------------------------------------------+
| nlam | 1 |
+--------------------+-------------------------------------------------+
| shape | (1,) |
+--------------------+-------------------------------------------------+
| elements | [Mg, Fe, Fe] |
+--------------------+-------------------------------------------------+
| line_ids | ['Mg 7 2628.89A, Fe 2 2631.05A, Fe 2 2631.32A'] |
+--------------------+-------------------------------------------------+
| lam | [2630.42] Å |
+--------------------+-------------------------------------------------+
| luminosity | [3.64386398e+31] erg/s |
+--------------------+-------------------------------------------------+
| continuum | [7.06006223e+19] erg/(Hz*s) |
+--------------------+-------------------------------------------------+
| cont | [7.06006223e+19] erg/(Hz*s) |
+--------------------+-------------------------------------------------+
| continuum_llam | [3.05899689e+31] erg/(s*Å) |
+--------------------+-------------------------------------------------+
| energy | [4.71347515] eV |
+--------------------+-------------------------------------------------+
| equivalent_width | [1.19119571] Å |
+--------------------+-------------------------------------------------+
| lum | [3.64386398e+31] erg/s |
+--------------------+-------------------------------------------------+
| nu | [1.13971327e+15] Hz |
+--------------------+-------------------------------------------------+
| vacuum_wavelengths | [2631.20435998] Å |
+--------------------+-------------------------------------------------+
Extracting a subset of lines¶
To extract a subset of lines we pass a list or tuple of the names of the lines we want to extract.
[7]:
print(
lines[("Mg 7 2628.89A", "Fe 2 2631.05A", "Fe 2 2631.32A", "Mg 5 2782.76A")]
)
+-----------------------------------------------------------------------------------------------+
| LINECOLLECTION |
+-------------------------+---------------------------------------------------------------------+
| Attribute | Value |
+-------------------------+---------------------------------------------------------------------+
| nlines | 4 |
+-------------------------+---------------------------------------------------------------------+
| ndim | 1 |
+-------------------------+---------------------------------------------------------------------+
| nlam | 4 |
+-------------------------+---------------------------------------------------------------------+
| shape | (4,) |
+-------------------------+---------------------------------------------------------------------+
| elements | [Mg, Fe, Fe, Mg] |
+-------------------------+---------------------------------------------------------------------+
| line_ids (4,) | [Mg 7 2628.89A, Fe 2 2631.05A, Fe 2 2631.32A, ...] |
+-------------------------+---------------------------------------------------------------------+
| lam (4,) | 2.63e+03 Å -> 2.78e+03 Å (Mean: 2.67e+03 Å) |
+-------------------------+---------------------------------------------------------------------+
| luminosity (4,) | 0.00e+00 erg/s -> 2.97e+31 erg/s (Mean: 9.11e+30 erg/s) |
+-------------------------+---------------------------------------------------------------------+
| continuum (4,) | 2.35e+19 erg/(Hz*s) -> 2.73e+19 erg/(Hz*s) (Mean: 2.45e+19 erg) |
+-------------------------+---------------------------------------------------------------------+
| cont (4,) | 2.35e+19 erg/(Hz*s) -> 2.73e+19 erg/(Hz*s) (Mean: 2.45e+19 erg) |
+-------------------------+---------------------------------------------------------------------+
| continuum_llam (4,) | 1.02e+31 erg/(s*Å) -> 1.06e+31 erg/(s*Å) (Mean: 1.03e+31 erg/(s*Å)) |
+-------------------------+---------------------------------------------------------------------+
| energy (4,) | 4.46e+00 eV -> 4.72e+00 eV (Mean: 4.65e+00 eV) |
+-------------------------+---------------------------------------------------------------------+
| equivalent_width (4,) | 0.00e+00 Å -> 2.91e+00 Å (Mean: 8.93e-01 Å) |
+-------------------------+---------------------------------------------------------------------+
| lum (4,) | 0.00e+00 erg/s -> 2.97e+31 erg/s (Mean: 9.11e+30 erg/s) |
+-------------------------+---------------------------------------------------------------------+
| nu (4,) | 1.08e+15 Hz -> 1.14e+15 Hz (Mean: 1.12e+15 Hz) |
+-------------------------+---------------------------------------------------------------------+
| vacuum_wavelengths (4,) | 2.63e+03 Å -> 2.78e+03 Å (Mean: 2.67e+03 Å) |
+-------------------------+---------------------------------------------------------------------+
Note that this list can include composite lines as well as single lines.
[8]:
print(
lines[("Mg 7 2628.89A, Fe 2 2631.05A", "Fe 2 2631.32A", "Mg 5 2782.76A")]
)
+---------------------------------------------------------------------------------------+
| LINECOLLECTION |
+--------------------+------------------------------------------------------------------+
| Attribute | Value |
+--------------------+------------------------------------------------------------------+
| nlines | 3 |
+--------------------+------------------------------------------------------------------+
| ndim | 1 |
+--------------------+------------------------------------------------------------------+
| nlam | 3 |
+--------------------+------------------------------------------------------------------+
| shape | (3,) |
+--------------------+------------------------------------------------------------------+
| elements | [Mg, Fe, Fe, Mg] |
+--------------------+------------------------------------------------------------------+
| line_ids | ['Mg 7 2628.89A, Fe 2 2631.05A' 'Fe 2 2631.32A' 'Mg 5 2782.76A'] |
+--------------------+------------------------------------------------------------------+
| lam | [2629.97 2631.32 2782.76] Å |
+--------------------+------------------------------------------------------------------+
| luminosity | [6.72173640e+30 2.97169034e+31 0.00000000e+00] erg/s |
+--------------------+------------------------------------------------------------------+
| continuum | [4.70482496e+19 2.35523727e+19 2.72580856e+19] erg/(Hz*s) |
+--------------------+------------------------------------------------------------------+
| cont | [4.70482496e+19 2.35523727e+19 2.72580856e+19] erg/(Hz*s) |
+--------------------+------------------------------------------------------------------+
| continuum_llam | [2.03921297e+31 1.01978363e+31 1.05527233e+31] erg/(s*Å) |
+--------------------+------------------------------------------------------------------+
| energy | [4.71428165 4.71186299 4.45543968] eV |
+--------------------+------------------------------------------------------------------+
| equivalent_width | [0.32962405 2.91404006 0. ] Å |
+--------------------+------------------------------------------------------------------+
| lum | [6.72173640e+30 2.97169034e+31 0.00000000e+00] erg/s |
+--------------------+------------------------------------------------------------------+
| nu | [1.13990828e+15 1.13932345e+15 1.07732057e+15] Hz |
+--------------------+------------------------------------------------------------------+
| vacuum_wavelengths | [2630.75425269 2632.10457459 2783.58110253] Å |
+--------------------+------------------------------------------------------------------+
Line aliases¶
Any of these operations can also be performed by passing some common line aliases defined in the emissions submodule.
[9]:
from synthesizer.emissions import line_aliases
for alias in line_aliases:
print(f"{alias} -> {line_aliases[alias]}")
# We can use any of these aliases in place of the long form ID
print(lines["Hb"])
Hb -> H 1 4861.32A
Ha -> H 1 6562.80A
Hg -> H 1 4340.46A
O1 -> O 1 6300.30A
O2b -> O 2 3726.03A
O2r -> O 2 3728.81A
O2 -> O 2 3726.03A, O 2 3728.81A
O3b -> O 3 4958.91A
O3r -> O 3 5006.84A
O3 -> O 3 4958.91A, O 3 5006.84A
Ne3 -> Ne 3 3868.76A
N2 -> N 2 6583.45A
S2 -> S 2 6730.82A, S 2 6716.44A
+--------------------------------------------------+
| LINECOLLECTION |
+--------------------+-----------------------------+
| Attribute | Value |
+--------------------+-----------------------------+
| nlines | 1 |
+--------------------+-----------------------------+
| id | 'H 1 4861.32A' |
+--------------------+-----------------------------+
| ndim | 1 |
+--------------------+-----------------------------+
| nlam | 1 |
+--------------------+-----------------------------+
| shape | (1,) |
+--------------------+-----------------------------+
| elements | [H, ] |
+--------------------+-----------------------------+
| line_ids | ['H 1 4861.32A'] |
+--------------------+-----------------------------+
| lam | [4861.32] Å |
+--------------------+-----------------------------+
| luminosity | [2.43960665e+33] erg/s |
+--------------------+-----------------------------+
| continuum | [1.31331716e+19] erg/(Hz*s) |
+--------------------+-----------------------------+
| cont | [1.31331716e+19] erg/(Hz*s) |
+--------------------+-----------------------------+
| continuum_llam | [1.66602648e+30] erg/(s*Å) |
+--------------------+-----------------------------+
| energy | [2.55042238] eV |
+--------------------+-----------------------------+
| equivalent_width | [1464.32645177] Å |
+--------------------+-----------------------------+
| lum | [2.43960665e+33] erg/s |
+--------------------+-----------------------------+
| nu | [6.16689414e+14] Hz |
+--------------------+-----------------------------+
| vacuum_wavelengths | [4862.6779924] Å |
+--------------------+-----------------------------+
Getting line ratios¶
In addition to the line extraction methods above, we can also get known ratios by passing the ratio name as an index. These ratios are defined in the line_ratios submodule (more details below).
[10]:
ratio = lines["BalmerDecrement"]
print(ratio)
2.9180493435290717
Alternatively, we could use the get_ratio method. This can also take a ratio name but also enables the passing of a list of two lines to get the ratio of.
[11]:
ratio2 = lines.get_ratio("BalmerDecrement")
ratio3 = lines.get_ratio(["H 1 4861.32A", "H 1 6562.80A"])
print(ratio2)
print(ratio3)
2.9180493435290717
0.34269468479604454
Getting diagnostic diagrams¶
We can also get diagnostic diagrams by passing the name of the diagram as an index. These diagrams are are also defined in the line_ratios submodule.
[12]:
dia = lines["BPT-NII"]
print(dia)
(0.11011948908117017, 3.2777857932958203)
Or we can use the get_diagram method to get a diagnostic diagram by passing the name of the diagram or a list of lists of lines.
[13]:
dia2 = lines.get_diagram("BPT-NII")
print(dia2)
dia3 = lines.get_diagram(
[["O 3 5006.84A", "H 1 4861.32A"], ["N 2 6583.45A", "H 1 6562.80A"]]
)
print(dia3)
(0.11011948908117017, 3.2777857932958203)
(3.2777857932958203, 0.11011948908117017)
Arithmetic operations¶
We can perform various arithmetic operations on a LineCollection object. These operations are performed element-wise on the underlying line data.
[14]:
# Get a subset for the demonstration
lines_subset = lines[("Fe 2 2631.05A", "Fe 2 2631.32A")]
Adding LineCollections¶
We can add two LineCollection objects together. This will return a new LineCollection object with the same lines as the two input LineCollection objects, but with the line data added together.
[15]:
new_line = lines_subset + lines_subset
print(new_line.luminosity / lines_subset.luminosity)
[2. 2.] dimensionless
Scaling a LineCollection¶
We can scale (multiply) a LineCollection object by a scalar. This will return a new LineCollection object with the same lines as the input LineCollection object, but with the line luminosity scaled accordingly.
[16]:
new_lines = lines_subset * 4
print(new_lines.luminosity / lines_subset.luminosity)
[4. 4.] dimensionless
Iterating over a LineCollection¶
We can loop over the individual lines just as we would loop over an array or list. Here we will demonstrate this on a subset.
[17]:
lines_subset = lines[
("Mg 7 2628.89A", "Fe 2 2631.05A", "Fe 2 2631.32A", "Mg 5 2782.76A")
]
for line in lines_subset:
print(
f"ID: {line.id}, Wavelength: {line.lam}, "
f"Luminosity: {line.luminosity}, Continuum: {line.continuum}"
)
ID: Mg 7 2628.89A, Wavelength: [2628.89] Å, Luminosity: [0.] erg/s, Continuum: [2.35015265e+19] erg/(Hz*s)
ID: Fe 2 2631.05A, Wavelength: [2631.05] Å, Luminosity: [6.7217364e+30] erg/s, Continuum: [2.35467231e+19] erg/(Hz*s)
ID: Fe 2 2631.32A, Wavelength: [2631.32] Å, Luminosity: [2.97169034e+31] erg/s, Continuum: [2.35523727e+19] erg/(Hz*s)
ID: Mg 5 2782.76A, Wavelength: [2782.76] Å, Luminosity: [0.] erg/s, Continuum: [2.72580856e+19] erg/(Hz*s)
Computing fluxes¶
By default a line contains the rest frame luminosity and continuum luminosity. If we instead want the flux we can use the get_flux or get_flux0 methods which each populate the flux and continuum_flux attributes on the line. The latter of these will compute the rest frame flux at a distance of 10pc.
[18]:
lines.get_flux0()
print(lines["H 1 4861.32A"].flux, lines["H 1 4861.32A"].continuum_flux)
[2.0389607e-07] erg/(cm**2*s) [1.09763681e+11] nJy
While the former requires an astropy.cosmology object and a redshift to compute the observer frame flux.
[19]:
from astropy.cosmology import Planck18 as cosmo
lines.get_flux(cosmo, 8.0)
print(lines["H 1 4861.32A"].flux, lines["H 1 4861.32A"].continuum_flux)
[3.01857317e-27] erg/(cm**2*s) [1.62499308e-09] nJy
Concatenating LineCollections¶
If rather than adding the values of two LineCollection objects together we want to instead concatenate them along the first axis, we can use the concatenate method. This will return a new LineCollection object with all the lines from the two input LineCollection objects.
[20]:
import numpy as np
from unyt import Hz, angstrom, erg, s
from synthesizer.emissions import LineCollection
# Simulating having two sets of ndim=2 lines we need to combine
lines1 = LineCollection(
line_ids=["O 3 5006.84A", "H 1 4861.32A"],
lam=[5006.84 * angstrom, 4861.32 * angstrom],
lum=np.array([[1.0, 2.0], [3.0, 4.0]]) * erg / s,
cont=np.array([[0.1, 0.2], [0.3, 0.4]]) * erg / s / Hz,
)
lines2 = LineCollection(
line_ids=["O 3 5006.84A", "H 1 4861.32A"],
lam=[5006.84 * angstrom, 4861.32 * angstrom],
lum=np.array([[5.0, 6.0], [7.0, 8.0]]) * erg / s,
cont=np.array([[0.5, 0.6], [0.7, 0.8]]) * erg / s / Hz,
)
print(lines1.shape, lines2.shape)
lines3 = lines1.concat(lines2)
print(lines3.shape)
(2, 2) (2, 2)
(4, 2)
Extending a LineCollection with new lines¶
If we want to add new lines to a LineCollection object we can use the extend method. This will return a new LineCollection object with all the lines from the original LineCollection object and the new lines. Note that the shape of the new lines must match the shape of the existing lines.
[21]:
# Set up two line collections with different lines
lines1 = LineCollection(
line_ids=["O 3 5006.84A", "H 1 4861.32A"],
lam=[5006.84 * angstrom, 4861.32 * angstrom],
lum=np.array([[1.0, 2.0], [3.0, 4.0]]) * erg / s,
cont=np.array([[0.1, 0.2], [0.3, 0.4]]) * erg / s / Hz,
)
lines2 = LineCollection(
line_ids=["N 2 6583.45A", "H 1 6562.80A"],
lam=[6583.45 * angstrom, 6562.80 * angstrom],
lum=np.array([[5.0, 6.0], [7.0, 8.0]]) * erg / s,
cont=np.array([[0.5, 0.6], [0.7, 0.8]]) * erg / s / Hz,
)
print(lines1.shape, lines2.shape)
lines3 = lines1.extend(lines2)
print(lines3.shape)
(2, 2) (2, 2)
(2, 4)
Blending lines¶
Lines in a LineCollection can be blended based on a given wavelength resolution using the get_blended_lines method. This method takes a set of wavelength bins, and returns a new LineCollection containing lines blended within each bin.
[22]:
import numpy as np
from unyt import angstrom
# Get a subset of lines to blend
lines_subset = lines["H 1 4861.32A", "O 3 5006.84A", "O 3 4958.91A"]
print("Before blending:")
print(lines_subset)
# Blend the lines onto an arbitrary wavelength grid
lam_bins = np.arange(4000, 7000, 1000) * angstrom
blended_lines = lines_subset.get_blended_lines(lam_bins)
print("After blending:")
print(blended_lines)
Before blending:
+-----------------------------------------------------------------------------------+
| LINECOLLECTION |
+--------------------+--------------------------------------------------------------+
| Attribute | Value |
+--------------------+--------------------------------------------------------------+
| nlines | 3 |
+--------------------+--------------------------------------------------------------+
| ndim | 1 |
+--------------------+--------------------------------------------------------------+
| nlam | 3 |
+--------------------+--------------------------------------------------------------+
| shape | (3,) |
+--------------------+--------------------------------------------------------------+
| available_ratios | [R3, ] |
+--------------------+--------------------------------------------------------------+
| elements | [H, O, O] |
+--------------------+--------------------------------------------------------------+
| line_ids | ['H 1 4861.32A' 'O 3 5006.84A' 'O 3 4958.91A'] |
+--------------------+--------------------------------------------------------------+
| lam | [4861.32 5006.84 4958.91] Å |
+--------------------+--------------------------------------------------------------+
| luminosity | [2.43960665e+33 7.99650802e+33 2.67784170e+33] erg/s |
+--------------------+--------------------------------------------------------------+
| continuum | [1.31331716e+19 1.34528228e+19 1.33876459e+19] erg/(Hz*s) |
+--------------------+--------------------------------------------------------------+
| obslam | [43751.88 45061.56 44630.19] Å |
+--------------------+--------------------------------------------------------------+
| flux | [3.01857317e-27 9.89423625e-27 3.31334607e-27] erg/(cm**2*s) |
+--------------------+--------------------------------------------------------------+
| continuum_flux | [1.62499308e-09 1.66454415e-09 1.65647970e-09] nJy |
+--------------------+--------------------------------------------------------------+
| cont | [1.31331716e+19 1.34528228e+19 1.33876459e+19] erg/(Hz*s) |
+--------------------+--------------------------------------------------------------+
| continuum_llam | [1.66602648e+30 1.60881719e+30 1.63212143e+30] erg/(s*Å) |
+--------------------+--------------------------------------------------------------+
| energy | [2.55042238 2.47629629 2.50023076] eV |
+--------------------+--------------------------------------------------------------+
| equivalent_width | [1464.32645177 4970.42677272 1640.71228946] Å |
+--------------------+--------------------------------------------------------------+
| lum | [2.43960665e+33 7.99650802e+33 2.67784170e+33] erg/s |
+--------------------+--------------------------------------------------------------+
| nu | [6.16689414e+14 5.98765804e+14 6.04553134e+14] Hz |
+--------------------+--------------------------------------------------------------+
| obsnu | [6.85210460e+13 6.65295338e+13 6.71725704e+13] Hz |
+--------------------+--------------------------------------------------------------+
| vacuum_wavelengths | [4862.6779924 5008.23663693 4960.29390118] Å |
+--------------------+--------------------------------------------------------------+
After blending:
+--------------------------------------------------------------------+
| LINECOLLECTION |
+--------------------+-----------------------------------------------+
| Attribute | Value |
+--------------------+-----------------------------------------------+
| nlines | 2 |
+--------------------+-----------------------------------------------+
| ndim | 1 |
+--------------------+-----------------------------------------------+
| nlam | 2 |
+--------------------+-----------------------------------------------+
| shape | (2,) |
+--------------------+-----------------------------------------------+
| available_ratios | [R3, ] |
+--------------------+-----------------------------------------------+
| elements | [H, O, O] |
+--------------------+-----------------------------------------------+
| line_ids | ['H 1 4861.32A, O 3 4958.91A' 'O 3 5006.84A'] |
+--------------------+-----------------------------------------------+
| lam | [4910.115 5006.84 ] Å |
+--------------------+-----------------------------------------------+
| luminosity | [5.11744835e+33 7.99650802e+33] erg/s |
+--------------------+-----------------------------------------------+
| continuum | [2.65208175e+19 1.34528228e+19] erg/(Hz*s) |
+--------------------+-----------------------------------------------+
| cont | [2.65208175e+19 1.34528228e+19] erg/(Hz*s) |
+--------------------+-----------------------------------------------+
| continuum_llam | [3.29779976e+30 1.60881719e+30] erg/(s*Å) |
+--------------------+-----------------------------------------------+
| energy | [2.52507718 2.47629629] eV |
+--------------------+-----------------------------------------------+
| equivalent_width | [1551.77655325 4970.42677272] Å |
+--------------------+-----------------------------------------------+
| lum | [5.11744835e+33 7.99650802e+33] erg/s |
+--------------------+-----------------------------------------------+
| nu | [6.10560971e+14 5.98765804e+14] Hz |
+--------------------+-----------------------------------------------+
| vacuum_wavelengths | [4911.48594295 5008.23663693] Å |
+--------------------+-----------------------------------------------+
Plotting Lines¶
To plot lines you can use a LineCollection instance’s plot_lines method. This can be incomprehensible if you plot all lines so we’ll pass a subset to the subset argument.
[23]:
_ = lines.plot_lines(
figsize=(12, 5),
subset=[
"He 2 1025.27A",
"He 2 1084.94A",
"Si 2 1179.59A",
"Si 3 1206.50A",
"He 2 1215.13A",
"H 1 1215.67A",
"Si 2 1264.74A",
"O 1 1302.17A",
"O 1 1304.86A",
],
xlimits=(None, 1350),
ylimits=(10**28.0, 10**36.0),
)
Line Ratios¶
To help with the analysis of emission lines, Synthesizer provides a submodule called line_ratios. This submodule contains definitions for commonly used line ratios and diagnostic diagrams that can be derived from the emission lines in a LineCollection.
We can print out the available line ratios with the available_ratios constant.
[24]:
from synthesizer.emissions import line_ratios
print(line_ratios.available_ratios)
# Print the ratios and their definitions
for ratio in line_ratios.available_ratios:
print(
f"{ratio} = {line_ratios.ratios[ratio][0]} "
f"/ {line_ratios.ratios[ratio][1]}"
)
('BalmerDecrement', 'N2', 'S2', 'O1', 'R2', 'R3', 'R23', 'O32', 'Ne3O2')
BalmerDecrement = H 1 6562.80A / H 1 4861.32A
N2 = N 2 6583.45A / H 1 6562.80A
S2 = S 2 6730.82A, S 2 6716.44A / H 1 6562.80A
O1 = O 1 6300.30A / H 1 6562.80A
R2 = O 2 3726.03A / H 1 4861.32A
R3 = O 3 5006.84A / H 1 4861.32A
R23 = O 3 4958.91A, O 3 5006.84A, O 2 3726.03A, O 2 3728.81A / H 1 4861.32A
O32 = O 3 5006.84A / O 2 3726.03A
Ne3O2 = Ne 3 3868.76A / O 2 3726.03A
Diagnostic Diagrams¶
In addition to line ratios, the submodule also provides a set of common diagnostic diagrams that can be used, for example, to classify the ionizing source.
[25]:
print(line_ratios.available_diagrams)
for diagram in line_ratios.available_diagrams:
print(
f"{diagram} = {line_ratios.diagrams[diagram][0]} "
f"/ {line_ratios.diagrams[diagram][1]}"
)
('OHNO', 'BPT-NII', 'VO78-SII', 'VO78-OI')
OHNO = ['O 3 5006.84A', 'H 1 4861.32A'] / ['Ne 3 3868.76A', 'O 2 3726.03A, O 2 3728.81A']
BPT-NII = ['N 2 6583.45A', 'H 1 6562.80A'] / ['O 3 5006.84A', 'H 1 4861.32A']
VO78-SII = ['S 2 6730.82A, S 2 6716.44A', 'H 1 6562.80A'] / ['O 3 5006.84A', 'H 1 4861.32A']
VO78-OI = ['O 1 6300.30A', 'H 1 6562.80A'] / ['O 3 5006.84A', 'H 1 4861.32A']
Combining line_ratios with a LineCollection¶
[26]:
from synthesizer.grid import Grid
grid = Grid("test_grid")
lines = grid.get_lines((1, 8))
As shown in the line collection docs, we can measure some predefined line ratios by passing the name of a ratio. We can also import all line ratios and loop over all pre-defined ratios.
[27]:
for ratio_id in lines.available_ratios:
ratio = lines.get_ratio(ratio_id)
print(f"{ratio_id}: {ratio:.2f}")
BalmerDecrement: 2.93
N2: 0.07
S2: 0.07
O1: 0.00
R2: 0.76
R3: 3.82
R23: 6.48
O32: 5.00
Ne3O2: 0.23
Similarly, we can also do same for all the diagnostic diagrams.
[28]:
for diagram_id in lines.available_diagrams:
diagram0, diagram1 = lines.get_diagram(diagram_id)
print(f"{diagram_id}: {diagram0:.2f} / {diagram1:.2f}")
OHNO: 3.82 / 0.13
BPT-NII: 0.07 / 3.82
VO78-SII: 0.07 / 3.82
VO78-OI: 0.00 / 3.82
BPT Classification Regions¶
The lines_ratios module also provides functions defining BPT classification regions from Kewley+13 and Kauffmann+03. These functions take a set of x values and return the y values that define the classification regions.
[29]:
import numpy as np
logNII_Ha = np.arange(-2.0, 1.0, 0.01)
logOIII_Hb_kewley = line_ratios.get_bpt_kewley01(logNII_Ha)
logOIII_Hb_kauffman = line_ratios.get_bpt_kauffman03(logNII_Ha)
We also provide functions to plot these classifications either on an existing set of axes or on a new figure.
[30]:
fig, ax = line_ratios.plot_bpt_kewley01(
logNII_Ha, show=False, fig=None, ax=None, label="Kewley+01"
)
_, _ = line_ratios.plot_bpt_kauffman03(
logNII_Ha, fig=fig, ax=ax, label="Kauffman+03", show=True
)
/home/runner/work/synthesizer/synthesizer/src/synthesizer/emissions/line_ratios.py:174: RuntimeWarning: overflow encountered in power
ax.loglog(10**logNII_Ha, 10**logOIII_Hb, **kwargs)
/home/runner/work/synthesizer/synthesizer/src/synthesizer/emissions/line_ratios.py:263: RuntimeWarning: overflow encountered in power
ax.loglog(10**logNII_Ha, 10**logOIII_Hb, **kwargs)
For a more interesting example, see the metallicity dependence plots in the grid lines example.
Observer frame lines¶
Just like an Sed lines can be shifted into the observer frame by calling the get_flux method with a cosmology and redshift.
[31]:
from astropy.cosmology import Planck18 as cosmo
lines.get_flux(cosmo, z=8)
for line in lines:
print(f"{line.id}: {line.flux} @ {line.obslam}")
He 2 1025.27A: [5.65058266e-31] erg/(cm**2*s) @ [9227.43] Å
O 6 1031.91A: [0.] erg/(cm**2*s) @ [9287.19] Å
O 6 1037.61A: [0.] erg/(cm**2*s) @ [9338.49] Å
He 2 1084.94A: [8.42308419e-31] erg/(cm**2*s) @ [9764.46] Å
Si 2 1179.59A: [3.50224506e-28] erg/(cm**2*s) @ [10616.31] Å
Si 3 1206.50A: [6.92626953e-28] erg/(cm**2*s) @ [10858.5] Å
He 2 1215.13A: [2.8904509e-30] erg/(cm**2*s) @ [10936.17] Å
H 1 1215.67A: [2.39667888e-25] erg/(cm**2*s) @ [10941.03] Å
O 5 1218.34A: [3.89705894e-36] erg/(cm**2*s) @ [10965.06] Å
N 5 1238.82A: [1.07296821e-31] erg/(cm**2*s) @ [11149.38] Å
N 5 1242.80A: [9.93728929e-32] erg/(cm**2*s) @ [11185.2] Å
Si 2 1260.42A: [1.15402657e-28] erg/(cm**2*s) @ [11343.78] Å
Si 2 1264.74A: [1.82934312e-28] erg/(cm**2*s) @ [11382.66] Å
Si 2 1265.00A: [2.48203028e-28] erg/(cm**2*s) @ [11385.] Å
O 1 1302.17A: [2.72957116e-29] erg/(cm**2*s) @ [11719.53] Å
O 1 1304.86A: [2.16645567e-28] erg/(cm**2*s) @ [11743.74] Å
O 1 1306.03A: [7.45460414e-29] erg/(cm**2*s) @ [11754.27] Å
C 2 1334.53A: [1.11022304e-28] erg/(cm**2*s) @ [12010.77] Å
C 2 1335.66A: [9.41121034e-29] erg/(cm**2*s) @ [12020.94] Å
C 2 1335.71A: [1.23974178e-28] erg/(cm**2*s) @ [12021.39] Å
Si 4 1393.75A: [2.01583556e-28] erg/(cm**2*s) @ [12543.75] Å
O 4 1399.78A: [6.45345335e-33] erg/(cm**2*s) @ [12598.02] Å
O 4 1401.16A: [3.6478445e-33] erg/(cm**2*s) @ [12610.44] Å
Si 4 1402.77A: [1.82395863e-28] erg/(cm**2*s) @ [12624.93] Å
O 4 1404.81A: [1.59664596e-33] erg/(cm**2*s) @ [12643.29] Å
O 4 1407.38A: [6.30487061e-33] erg/(cm**2*s) @ [12666.42] Å
N 4 1486.50A: [3.39833901e-31] erg/(cm**2*s) @ [13378.5] Å
Si 2 1526.71A: [3.33845478e-29] erg/(cm**2*s) @ [13740.39] Å
Si 2 1533.43A: [1.10665763e-28] erg/(cm**2*s) @ [13800.87] Å
C 4 1548.19A: [9.11506355e-29] erg/(cm**2*s) @ [13933.71] Å
C 4 1550.77A: [2.77126455e-28] erg/(cm**2*s) @ [13956.93] Å
Ne 4 1601.45A: [2.29531955e-34] erg/(cm**2*s) @ [14413.05] Å
He 2 1640.41A: [1.40680641e-29] erg/(cm**2*s) @ [14763.69] Å
O 1 1641.31A: [7.82569717e-33] erg/(cm**2*s) @ [14771.79] Å
C 1 1657.91A: [2.2315061e-30] erg/(cm**2*s) @ [14921.19] Å
O 3 1660.81A: [5.01434002e-29] erg/(cm**2*s) @ [14947.29] Å
O 3 1666.15A: [1.47222065e-28] erg/(cm**2*s) @ [14995.35] Å
Al 2 1670.79A: [0.] erg/(cm**2*s) @ [15037.11] Å
N 3 1749.67A: [2.66495731e-29] erg/(cm**2*s) @ [15747.03] Å
Mg 6 1806.00A: [0.] erg/(cm**2*s) @ [16254.] Å
Si 3 1882.71A: [1.11400673e-27] erg/(cm**2*s) @ [16944.39] Å
Si 3 1892.03A: [8.5387779e-28] erg/(cm**2*s) @ [17028.27] Å
C 3 1906.68A: [1.18254979e-27] erg/(cm**2*s) @ [17160.12] Å
C 3 1908.73A: [8.09855969e-28] erg/(cm**2*s) @ [17178.57] Å
C 1 1992.01A: [2.5339833e-35] erg/(cm**2*s) @ [17928.09] Å
Si 7 2146.64A: [0.] erg/(cm**2*s) @ [19319.76] Å
O 3 2320.95A: [1.05266723e-28] erg/(cm**2*s) @ [20888.55] Å
C 2 2325.40A: [7.91883425e-28] erg/(cm**2*s) @ [20928.6] Å
C 2 2326.93A: [4.53133697e-28] erg/(cm**2*s) @ [20942.37] Å
Fe 2 2395.63A: [1.29649718e-28] erg/(cm**2*s) @ [21560.67] Å
Fe 2 2399.24A: [2.52801671e-29] erg/(cm**2*s) @ [21593.16] Å
Fe 2 2406.66A: [2.42803988e-29] erg/(cm**2*s) @ [21659.94] Å
Fe 2 2410.52A: [2.19388745e-29] erg/(cm**2*s) @ [21694.68] Å
C 1 2582.90A: [9.19219269e-35] erg/(cm**2*s) @ [23246.1] Å
Fe 2 2598.37A: [2.54417017e-29] erg/(cm**2*s) @ [23385.33] Å
Fe 2 2607.09A: [3.52567486e-29] erg/(cm**2*s) @ [23463.81] Å
Fe 2 2611.87A: [1.1613029e-28] erg/(cm**2*s) @ [23506.83] Å
Fe 2 2613.82A: [2.17729741e-29] erg/(cm**2*s) @ [23524.38] Å
Fe 2 2625.67A: [7.21957855e-28] erg/(cm**2*s) @ [23631.03] Å
Fe 2 2628.29A: [1.61156244e-29] erg/(cm**2*s) @ [23654.61] Å
Mg 7 2628.89A: [0.] erg/(cm**2*s) @ [23660.01] Å
Fe 2 2631.05A: [3.19540206e-29] erg/(cm**2*s) @ [23679.45] Å
Fe 2 2631.32A: [1.73715197e-28] erg/(cm**2*s) @ [23681.88] Å
Mg 5 2782.76A: [0.] erg/(cm**2*s) @ [25044.84] Å
Mg 2 2795.53A: [8.02434109e-27] erg/(cm**2*s) @ [25159.77] Å
Mg 2 2802.71A: [3.98757654e-27] erg/(cm**2*s) @ [25224.39] Å
Fe 4 2829.36A: [1.15162538e-27] erg/(cm**2*s) @ [25464.24] Å
Fe 4 2835.74A: [2.08677982e-27] erg/(cm**2*s) @ [25521.66] Å
Ar 4 2853.66A: [3.18725256e-30] erg/(cm**2*s) @ [25682.94] Å
Fe 4 3094.96A: [5.09940036e-28] erg/(cm**2*s) @ [27854.64] Å
He 1 3187.74A: [1.13394388e-27] erg/(cm**2*s) @ [28689.66] Å
Ne 5 3345.82A: [0.] erg/(cm**2*s) @ [30112.38] Å
Ne 5 3425.88A: [0.] erg/(cm**2*s) @ [30832.92] Å
Fe 7 3586.32A: [0.] erg/(cm**2*s) @ [32276.88] Å
Fe 6 3662.50A: [7.40075927e-36] erg/(cm**2*s) @ [32962.5] Å
O 2 3726.03A: [2.9692536e-26] erg/(cm**2*s) @ [33534.27] Å
O 2 3728.81A: [2.42047249e-26] erg/(cm**2*s) @ [33559.29] Å
H 1 3734.37A: [1.38716031e-27] erg/(cm**2*s) @ [33609.33] Å
H 1 3750.15A: [1.68832777e-27] erg/(cm**2*s) @ [33751.35] Å
Fe 7 3758.92A: [0.] erg/(cm**2*s) @ [33830.28] Å
H 1 3770.63A: [2.13538573e-27] erg/(cm**2*s) @ [33935.67] Å
H 1 3797.90A: [2.03111586e-27] erg/(cm**2*s) @ [34181.1] Å
H 1 3835.38A: [2.78894604e-27] erg/(cm**2*s) @ [34518.42] Å
Ne 3 3868.76A: [6.86260161e-27] erg/(cm**2*s) @ [34818.84] Å
He 1 3888.64A: [2.44783151e-27] erg/(cm**2*s) @ [34997.76] Å
H 1 3889.05A: [4.08478836e-27] erg/(cm**2*s) @ [35001.45] Å
Fe 5 3891.28A: [1.51100351e-32] erg/(cm**2*s) @ [35021.52] Å
Ne 3 3967.47A: [2.09041798e-27] erg/(cm**2*s) @ [35707.23] Å
H 1 3970.07A: [6.10739091e-27] erg/(cm**2*s) @ [35730.63] Å
S 2 4068.60A: [4.65576228e-28] erg/(cm**2*s) @ [36617.4] Å
S 2 4076.35A: [1.48742819e-28] erg/(cm**2*s) @ [36687.15] Å
H 1 4101.73A: [9.84415221e-27] erg/(cm**2*s) @ [36915.57] Å
Fe 2 4243.97A: [2.16565528e-28] erg/(cm**2*s) @ [38195.73] Å
Fe 2 4276.84A: [6.18895861e-29] erg/(cm**2*s) @ [38491.56] Å
Fe 2 4287.39A: [2.06262239e-28] erg/(cm**2*s) @ [38586.51] Å
Fe 2 4319.62A: [1.96380977e-29] erg/(cm**2*s) @ [38876.58] Å
H 1 4340.46A: [1.80242068e-26] erg/(cm**2*s) @ [39064.14] Å
Fe 2 4346.86A: [5.12283016e-29] erg/(cm**2*s) @ [39121.74] Å
Fe 2 4352.79A: [2.99879551e-29] erg/(cm**2*s) @ [39175.11] Å
Fe 2 4358.37A: [9.69498736e-30] erg/(cm**2*s) @ [39225.33] Å
Fe 2 4359.33A: [1.52897605e-28] erg/(cm**2*s) @ [39233.97] Å
O 3 4363.21A: [4.51598391e-28] erg/(cm**2*s) @ [39268.89] Å
Fe 2 4413.78A: [1.08133599e-28] erg/(cm**2*s) @ [39724.02] Å
Fe 2 4416.27A: [3.2207237e-29] erg/(cm**2*s) @ [39746.43] Å
Fe 2 4452.10A: [6.94911722e-29] erg/(cm**2*s) @ [40068.9] Å
Fe 2 4474.90A: [3.38786109e-29] erg/(cm**2*s) @ [40274.1] Å
Fe 3 4658.05A: [4.78702735e-27] erg/(cm**2*s) @ [41922.45] Å
He 2 4685.68A: [1.45352047e-30] erg/(cm**2*s) @ [42171.12] Å
Fe 2 4814.54A: [5.89812354e-29] erg/(cm**2*s) @ [43330.86] Å
H 1 4861.32A: [3.88498023e-26] erg/(cm**2*s) @ [43751.88] Å
Fe 2 4874.50A: [2.44164662e-29] erg/(cm**2*s) @ [43870.5] Å
Fe 2 4889.62A: [4.10035971e-29] erg/(cm**2*s) @ [44006.58] Å
Fe 2 4905.35A: [2.87648034e-29] erg/(cm**2*s) @ [44148.15] Å
Fe 2 4923.92A: [5.7102486e-31] erg/(cm**2*s) @ [44315.28] Å
Fe 2 4947.39A: [6.73708485e-29] erg/(cm**2*s) @ [44526.51] Å
O 3 4958.91A: [4.96478819e-26] erg/(cm**2*s) @ [44630.19] Å
Fe 2 4973.40A: [1.94059827e-29] erg/(cm**2*s) @ [44760.6] Å
Fe 3 4985.87A: [5.67401246e-28] erg/(cm**2*s) @ [44872.83] Å
Fe 2 5005.52A: [3.83593774e-29] erg/(cm**2*s) @ [45049.68] Å
O 3 5006.84A: [1.48364449e-25] erg/(cm**2*s) @ [45061.56] Å
Fe 2 5018.44A: [8.54823713e-31] erg/(cm**2*s) @ [45165.96] Å
Fe 2 5020.25A: [1.11648068e-29] erg/(cm**2*s) @ [45182.25] Å
Fe 2 5049.30A: [4.53221013e-29] erg/(cm**2*s) @ [45443.7] Å
Fe 2 5072.41A: [1.09297812e-28] erg/(cm**2*s) @ [45651.69] Å
Fe 2 5111.64A: [6.37239499e-29] erg/(cm**2*s) @ [46004.76] Å
Fe 2 5158.01A: [1.81464493e-29] erg/(cm**2*s) @ [46422.09] Å
Fe 2 5158.79A: [5.55802769e-28] erg/(cm**2*s) @ [46429.11] Å
Fe 2 5169.03A: [4.91094331e-30] erg/(cm**2*s) @ [46521.27] Å
Fe 6 5176.04A: [5.2539497e-35] erg/(cm**2*s) @ [46584.36] Å
Fe 2 5184.80A: [1.57117879e-29] erg/(cm**2*s) @ [46663.2] Å
Fe 2 5261.63A: [2.15175748e-28] erg/(cm**2*s) @ [47354.67] Å
Fe 3 5270.40A: [2.83827959e-27] erg/(cm**2*s) @ [47433.6] Å
Fe 2 5273.36A: [1.26775567e-28] erg/(cm**2*s) @ [47460.24] Å
Fe 2 5284.10A: [2.40912141e-31] erg/(cm**2*s) @ [47556.9] Å
Fe 2 5333.66A: [1.06343621e-28] erg/(cm**2*s) @ [48002.94] Å
Fe 2 5376.47A: [4.87427059e-29] erg/(cm**2*s) @ [48388.23] Å
Fe 2 5412.67A: [1.03302113e-29] erg/(cm**2*s) @ [48714.03] Å
Fe 2 5433.15A: [3.89851421e-29] erg/(cm**2*s) @ [48898.35] Å
Fe 2 5527.36A: [7.65688622e-29] erg/(cm**2*s) @ [49746.24] Å
Fe 7 5720.71A: [0.] erg/(cm**2*s) @ [51486.39] Å
He 1 5875.61A: [2.65542531e-27] erg/(cm**2*s) @ [52880.49] Å
He 1 5875.64A: [1.59322608e-27] erg/(cm**2*s) @ [52880.76] Å
Fe 7 6086.97A: [0.] erg/(cm**2*s) @ [54782.73] Å
O 1 6300.30A: [5.44262504e-28] erg/(cm**2*s) @ [56702.7] Å
O 1 6363.78A: [1.7408629e-28] erg/(cm**2*s) @ [57274.02] Å
Fe 2 6516.08A: [2.15379485e-32] erg/(cm**2*s) @ [58644.72] Å
N 2 6548.05A: [2.86571137e-27] erg/(cm**2*s) @ [58932.45] Å
H 1 6562.80A: [1.14020153e-25] erg/(cm**2*s) @ [59065.2] Å
N 2 6583.45A: [8.45510006e-27] erg/(cm**2*s) @ [59251.05] Å
Ni 2 6666.80A: [1.707319e-29] erg/(cm**2*s) @ [60001.2] Å
He 1 6678.15A: [1.31643434e-27] erg/(cm**2*s) @ [60103.35] Å
S 2 6716.44A: [3.63562026e-27] erg/(cm**2*s) @ [60447.96] Å
S 2 6730.82A: [4.47494538e-27] erg/(cm**2*s) @ [60577.38] Å
Ar 3 7135.79A: [4.62192733e-27] erg/(cm**2*s) @ [64222.11] Å
Fe 2 7155.17A: [6.16181788e-28] erg/(cm**2*s) @ [64396.53] Å
Fe 2 7172.00A: [1.15577289e-28] erg/(cm**2*s) @ [64548.] Å
Ca 2 7291.47A: [8.68430458e-28] erg/(cm**2*s) @ [65623.23] Å
Ca 2 7323.89A: [5.72443745e-28] erg/(cm**2*s) @ [65915.01] Å
Ni 2 7377.83A: [4.5503282e-28] erg/(cm**2*s) @ [66400.47] Å
Fe 2 7388.17A: [8.79636012e-29] erg/(cm**2*s) @ [66493.53] Å
Ni 2 7411.61A: [2.86134557e-29] erg/(cm**2*s) @ [66704.49] Å
Fe 2 7452.56A: [1.99400656e-28] erg/(cm**2*s) @ [67073.04] Å
Ar 3 7751.11A: [1.10840395e-27] erg/(cm**2*s) @ [69759.99] Å
O 1 8446.25A: [3.15516393e-30] erg/(cm**2*s) @ [76016.25] Å
O 1 8446.36A: [7.46260811e-30] erg/(cm**2*s) @ [76017.24] Å
O 1 8446.76A: [6.67509048e-30] erg/(cm**2*s) @ [76020.84] Å
Cl 2 8578.70A: [1.28281768e-28] erg/(cm**2*s) @ [77208.3] Å
Fe 2 8616.95A: [2.87931811e-28] erg/(cm**2*s) @ [77552.55] Å
Fe 2 8891.93A: [1.33644426e-28] erg/(cm**2*s) @ [80027.37] Å
Fe 2 9051.95A: [4.16926659e-28] erg/(cm**2*s) @ [81467.55] Å
S 3 9068.62A: [1.33324267e-26] erg/(cm**2*s) @ [81617.58] Å
Fe 2 9226.63A: [1.52803013e-28] erg/(cm**2*s) @ [83039.67] Å
H 1 9229.02A: [1.09239602e-27] erg/(cm**2*s) @ [83061.18] Å
Fe 2 9267.56A: [6.15250418e-29] erg/(cm**2*s) @ [83408.04] Å
Fe 2 9399.04A: [1.58391238e-28] erg/(cm**2*s) @ [84591.36] Å
Fe 2 9470.94A: [6.23407188e-29] erg/(cm**2*s) @ [85238.46] Å
S 3 9530.62A: [3.36748736e-26] erg/(cm**2*s) @ [85775.58] Å
H 1 9545.97A: [1.58464001e-27] erg/(cm**2*s) @ [85913.73] Å
H 1 1.00494m: [2.39304071e-27] erg/(cm**2*s) @ [90444.6] Å
S 2 1.02867m: [1.23624914e-28] erg/(cm**2*s) @ [92580.3] Å
S 2 1.03205m: [1.68243394e-28] erg/(cm**2*s) @ [92884.5] Å
S 2 1.03364m: [1.21471119e-28] erg/(cm**2*s) @ [93027.6] Å
He 1 1.08291m: [1.26098867e-27] erg/(cm**2*s) @ [97461.9] Å
He 1 1.08303m: [7.88026968e-27] erg/(cm**2*s) @ [97472.7] Å
H 1 1.09381m: [3.88905497e-27] erg/(cm**2*s) @ [98442.9] Å
O 1 1.12863m: [1.29758863e-30] erg/(cm**2*s) @ [101576.7] Å
O 1 1.12864m: [1.14842379e-30] erg/(cm**2*s) @ [101577.6] Å
O 1 1.12869m: [7.47788841e-31] erg/(cm**2*s) @ [101582.1] Å
O 1 1.12870m: [4.3288366e-31] erg/(cm**2*s) @ [101583.] Å
O 1 1.12873m: [1.53108619e-30] erg/(cm**2*s) @ [101585.7] Å
Ni 2 1.19102m: [9.91764318e-29] erg/(cm**2*s) @ [107191.8] Å
Fe 2 1.25668m: [3.08254612e-28] erg/(cm**2*s) @ [113101.2] Å
Fe 2 1.27877m: [1.31563394e-29] erg/(cm**2*s) @ [115089.3] Å
H 1 1.28181m: [7.04574693e-27] erg/(cm**2*s) @ [115362.9] Å
Fe 2 1.29427m: [2.6608098e-29] erg/(cm**2*s) @ [116484.3] Å
Fe 2 1.32055m: [9.86452595e-29] erg/(cm**2*s) @ [118849.5] Å
Fe 2 1.32777m: [1.7466112e-29] erg/(cm**2*s) @ [119499.3] Å
Fe 2 1.37181m: [5.64607133e-29] erg/(cm**2*s) @ [123462.9] Å
Fe 2 1.53348m: [6.29330124e-28] erg/(cm**2*s) @ [138013.2] Å
Fe 2 1.59948m: [2.52685249e-28] erg/(cm**2*s) @ [143953.2] Å
Fe 2 1.64355m: [4.73085405e-27] erg/(cm**2*s) @ [147919.5] Å
Fe 2 1.66377m: [1.30755721e-28] erg/(cm**2*s) @ [149739.3] Å
Fe 2 1.67688m: [2.37739659e-28] erg/(cm**2*s) @ [150919.2] Å
Fe 2 1.71113m: [5.32991462e-29] erg/(cm**2*s) @ [154001.7] Å
Fe 2 1.74494m: [6.24025676e-29] erg/(cm**2*s) @ [157044.6] Å
Fe 2 1.79711m: [1.18786152e-28] erg/(cm**2*s) @ [161739.9] Å
Fe 2 1.80002m: [2.30361457e-28] erg/(cm**2*s) @ [162001.8] Å
Fe 2 1.80940m: [1.19368259e-27] erg/(cm**2*s) @ [162846.] Å
H 1 1.87510m: [1.43569345e-26] erg/(cm**2*s) @ [168759.] Å
Fe 2 1.89541m: [4.24530428e-29] erg/(cm**2*s) @ [170586.9] Å
Ni 2 1.93877m: [9.673886e-29] erg/(cm**2*s) @ [174489.3] Å
Fe 2 1.95361m: [1.47163854e-28] erg/(cm**2*s) @ [175824.9] Å
Si 6 1.96247m: [0.] erg/(cm**2*s) @ [176622.3] Å
H 1 2.16553m: [1.23661296e-27] erg/(cm**2*s) @ [194897.7] Å
Si 7 2.48071m: [0.] erg/(cm**2*s) @ [223263.9] Å
Additionally, just like an Sed object, you can also calculate the rest frame flux at 10 parsecs using the get_flux0 method (again both on a Line and a LineCollection).
[32]:
lines.get_flux0()
for line in lines:
print(f"{line.id}: {line.flux} @ {line.obslam}")
lines[
"He 2 1025.27A",
"He 2 1084.94A",
"Si 2 1179.59A",
"Si 3 1206.50A",
"He 2 1215.13A",
"H 1 1215.67A",
"Si 2 1264.74A",
"O 1 1302.17A",
"O 1 1304.86A",
].plot_lines(ylimits=(None, 10**38.0), xlimits=(None, 1500))
He 2 1025.27A: [3.81680857e-11] erg/(cm**2*s) @ [1025.27] Å
O 6 1031.91A: [0.] erg/(cm**2*s) @ [1031.91] Å
O 6 1037.61A: [0.] erg/(cm**2*s) @ [1037.61] Å
He 2 1084.94A: [5.68955484e-11] erg/(cm**2*s) @ [1084.94] Å
Si 2 1179.59A: [2.36566736e-08] erg/(cm**2*s) @ [1179.59] Å
Si 3 1206.50A: [4.67849893e-08] erg/(cm**2*s) @ [1206.5] Å
He 2 1215.13A: [1.95241773e-10] erg/(cm**2*s) @ [1215.13] Å
H 1 1215.67A: [1.61888871e-05] erg/(cm**2*s) @ [1215.67] Å
O 5 1218.34A: [2.63235296e-16] erg/(cm**2*s) @ [1218.34] Å
N 5 1238.82A: [7.24759638e-12] erg/(cm**2*s) @ [1238.82] Å
N 5 1242.80A: [6.71235751e-12] erg/(cm**2*s) @ [1242.8] Å
Si 2 1260.42A: [7.79512265e-09] erg/(cm**2*s) @ [1260.42] Å
Si 2 1264.74A: [1.23566947e-08] erg/(cm**2*s) @ [1264.74] Å
Si 2 1265.00A: [1.67654117e-08] erg/(cm**2*s) @ [1265.] Å
O 1 1302.17A: [1.84374802e-09] erg/(cm**2*s) @ [1302.17] Å
O 1 1304.86A: [1.46337946e-08] erg/(cm**2*s) @ [1304.86] Å
O 1 1306.03A: [5.035374e-09] erg/(cm**2*s) @ [1306.03] Å
C 2 1334.53A: [7.49924221e-09] erg/(cm**2*s) @ [1334.53] Å
C 2 1335.66A: [6.35700608e-09] erg/(cm**2*s) @ [1335.66] Å
C 2 1335.71A: [8.37410465e-09] erg/(cm**2*s) @ [1335.71] Å
Si 4 1393.75A: [1.36163984e-08] erg/(cm**2*s) @ [1393.75] Å
O 4 1399.78A: [4.35912499e-13] erg/(cm**2*s) @ [1399.78] Å
O 4 1401.16A: [2.46401566e-13] erg/(cm**2*s) @ [1401.16] Å
Si 4 1402.77A: [1.23203241e-08] erg/(cm**2*s) @ [1402.77] Å
O 4 1404.81A: [1.07848913e-13] erg/(cm**2*s) @ [1404.81] Å
O 4 1407.38A: [4.25876156e-13] erg/(cm**2*s) @ [1407.38] Å
N 4 1486.50A: [2.29548177e-11] erg/(cm**2*s) @ [1486.5] Å
Si 2 1526.71A: [2.25503167e-09] erg/(cm**2*s) @ [1526.71] Å
Si 2 1533.43A: [7.47515892e-09] erg/(cm**2*s) @ [1533.43] Å
C 4 1548.19A: [6.15696731e-09] erg/(cm**2*s) @ [1548.19] Å
C 4 1550.77A: [1.87191073e-08] erg/(cm**2*s) @ [1550.77] Å
Ne 4 1601.45A: [1.55042335e-14] erg/(cm**2*s) @ [1601.45] Å
He 2 1640.41A: [9.5025789e-10] erg/(cm**2*s) @ [1640.41] Å
O 1 1641.31A: [5.28603683e-13] erg/(cm**2*s) @ [1641.31] Å
C 1 1657.91A: [1.50731918e-10] erg/(cm**2*s) @ [1657.91] Å
O 3 1660.81A: [3.38704469e-09] erg/(cm**2*s) @ [1660.81] Å
O 3 1666.15A: [9.94443358e-09] erg/(cm**2*s) @ [1666.15] Å
Al 2 1670.79A: [0.] erg/(cm**2*s) @ [1670.79] Å
N 3 1749.67A: [1.8001032e-09] erg/(cm**2*s) @ [1749.67] Å
Mg 6 1806.00A: [0.] erg/(cm**2*s) @ [1806.] Å
Si 3 1882.71A: [7.52479999e-08] erg/(cm**2*s) @ [1882.71] Å
Si 3 1892.03A: [5.76770267e-08] erg/(cm**2*s) @ [1892.03] Å
C 3 1906.68A: [7.98778899e-08] erg/(cm**2*s) @ [1906.68] Å
C 3 1908.73A: [5.47034774e-08] erg/(cm**2*s) @ [1908.73] Å
C 1 1992.01A: [1.71163396e-15] erg/(cm**2*s) @ [1992.01] Å
Si 7 2146.64A: [0.] erg/(cm**2*s) @ [2146.64] Å
O 3 2320.95A: [7.11046907e-09] erg/(cm**2*s) @ [2320.95] Å
C 2 2325.40A: [5.34894829e-08] erg/(cm**2*s) @ [2325.4] Å
C 2 2326.93A: [3.06078981e-08] erg/(cm**2*s) @ [2326.93] Å
Fe 2 2395.63A: [8.75747134e-09] erg/(cm**2*s) @ [2395.63] Å
Fe 2 2399.24A: [1.7076037e-09] erg/(cm**2*s) @ [2399.24] Å
Fe 2 2406.66A: [1.64007218e-09] erg/(cm**2*s) @ [2406.66] Å
Fe 2 2410.52A: [1.48190885e-09] erg/(cm**2*s) @ [2410.52] Å
C 1 2582.90A: [6.20906586e-15] erg/(cm**2*s) @ [2582.9] Å
Fe 2 2598.37A: [1.7185149e-09] erg/(cm**2*s) @ [2598.37] Å
Fe 2 2607.09A: [2.38149352e-09] erg/(cm**2*s) @ [2607.09] Å
Fe 2 2611.87A: [7.84427223e-09] erg/(cm**2*s) @ [2611.87] Å
Fe 2 2613.82A: [1.47070274e-09] erg/(cm**2*s) @ [2613.82] Å
Fe 2 2625.67A: [4.87662087e-08] erg/(cm**2*s) @ [2625.67] Å
Fe 2 2628.29A: [1.0885648e-09] erg/(cm**2*s) @ [2628.29] Å
Mg 7 2628.89A: [0.] erg/(cm**2*s) @ [2628.89] Å
Fe 2 2631.05A: [2.1584036e-09] erg/(cm**2*s) @ [2631.05] Å
Fe 2 2631.32A: [1.17339696e-08] erg/(cm**2*s) @ [2631.32] Å
Mg 5 2782.76A: [0.] erg/(cm**2*s) @ [2782.76] Å
Mg 2 2795.53A: [5.42021517e-07] erg/(cm**2*s) @ [2795.53] Å
Mg 2 2802.71A: [2.69349503e-07] erg/(cm**2*s) @ [2802.71] Å
Fe 4 2829.36A: [7.77890329e-08] erg/(cm**2*s) @ [2829.36] Å
Fe 4 2835.74A: [1.40956067e-07] erg/(cm**2*s) @ [2835.74] Å
Ar 4 2853.66A: [2.15289885e-10] erg/(cm**2*s) @ [2853.66] Å
Fe 4 3094.96A: [3.44450054e-08] erg/(cm**2*s) @ [3094.96] Å
He 1 3187.74A: [7.65946983e-08] erg/(cm**2*s) @ [3187.74] Å
Ne 5 3345.82A: [0.] erg/(cm**2*s) @ [3345.82] Å
Ne 5 3425.88A: [0.] erg/(cm**2*s) @ [3425.88] Å
Fe 7 3586.32A: [0.] erg/(cm**2*s) @ [3586.32] Å
Fe 6 3662.50A: [4.99900331e-16] erg/(cm**2*s) @ [3662.5] Å
O 2 3726.03A: [2.00564672e-06] erg/(cm**2*s) @ [3726.03] Å
O 2 3728.81A: [1.63496062e-06] erg/(cm**2*s) @ [3728.81] Å
H 1 3734.37A: [9.36987505e-08] erg/(cm**2*s) @ [3734.37] Å
H 1 3750.15A: [1.1404176e-07] erg/(cm**2*s) @ [3750.15] Å
Fe 7 3758.92A: [0.] erg/(cm**2*s) @ [3758.92] Å
H 1 3770.63A: [1.44239259e-07] erg/(cm**2*s) @ [3770.63] Å
H 1 3797.90A: [1.37196125e-07] erg/(cm**2*s) @ [3797.9] Å
H 1 3835.38A: [1.88385407e-07] erg/(cm**2*s) @ [3835.38] Å
Ne 3 3868.76A: [4.63549305e-07] erg/(cm**2*s) @ [3868.76] Å
He 1 3888.64A: [1.65344086e-07] erg/(cm**2*s) @ [3888.64] Å
H 1 3889.05A: [2.75915886e-07] erg/(cm**2*s) @ [3889.05] Å
Fe 5 3891.28A: [1.02064008e-12] erg/(cm**2*s) @ [3891.28] Å
Ne 3 3967.47A: [1.41201815e-07] erg/(cm**2*s) @ [3967.47] Å
H 1 3970.07A: [4.12536961e-07] erg/(cm**2*s) @ [3970.07] Å
S 2 4068.60A: [3.14483558e-08] erg/(cm**2*s) @ [4068.6] Å
S 2 4076.35A: [1.00471562e-08] erg/(cm**2*s) @ [4076.35] Å
H 1 4101.73A: [6.64944605e-07] erg/(cm**2*s) @ [4101.73] Å
Fe 2 4243.97A: [1.46283881e-08] erg/(cm**2*s) @ [4243.97] Å
Fe 2 4276.84A: [4.18046628e-09] erg/(cm**2*s) @ [4276.84] Å
Fe 2 4287.39A: [1.39324301e-08] erg/(cm**2*s) @ [4287.39] Å
Fe 2 4319.62A: [1.32649789e-09] erg/(cm**2*s) @ [4319.62] Å
H 1 4340.46A: [1.21748413e-06] erg/(cm**2*s) @ [4340.46] Å
Fe 2 4346.86A: [3.4603267e-09] erg/(cm**2*s) @ [4346.86] Å
Fe 2 4352.79A: [2.02560145e-09] erg/(cm**2*s) @ [4352.79] Å
Fe 2 4358.37A: [6.54868942e-10] erg/(cm**2*s) @ [4358.37] Å
Fe 2 4359.33A: [1.03278003e-08] erg/(cm**2*s) @ [4359.33] Å
O 3 4363.21A: [3.05041925e-08] erg/(cm**2*s) @ [4363.21] Å
Fe 2 4413.78A: [7.3041184e-09] erg/(cm**2*s) @ [4413.78] Å
Fe 2 4416.27A: [2.17550765e-09] erg/(cm**2*s) @ [4416.27] Å
Fe 2 4452.10A: [4.6939319e-09] erg/(cm**2*s) @ [4452.1] Å
Fe 2 4474.90A: [2.28840423e-09] erg/(cm**2*s) @ [4474.9] Å
Fe 3 4658.05A: [3.23350141e-07] erg/(cm**2*s) @ [4658.05] Å
He 2 4685.68A: [9.81811918e-11] erg/(cm**2*s) @ [4685.68] Å
Fe 2 4814.54A: [3.98401543e-09] erg/(cm**2*s) @ [4814.54] Å
H 1 4861.32A: [2.62419413e-06] erg/(cm**2*s) @ [4861.32] Å
Fe 2 4874.50A: [1.64926315e-09] erg/(cm**2*s) @ [4874.5] Å
Fe 2 4889.62A: [2.76967687e-09] erg/(cm**2*s) @ [4889.62] Å
Fe 2 4905.35A: [1.94298101e-09] erg/(cm**2*s) @ [4905.35] Å
Fe 2 4923.92A: [3.85711123e-11] erg/(cm**2*s) @ [4923.92] Å
Fe 2 4947.39A: [4.55071003e-09] erg/(cm**2*s) @ [4947.39] Å
O 3 4958.91A: [3.35357383e-06] erg/(cm**2*s) @ [4958.91] Å
Fe 2 4973.40A: [1.31081918e-09] erg/(cm**2*s) @ [4973.4] Å
Fe 3 4985.87A: [3.83263474e-08] erg/(cm**2*s) @ [4985.87] Å
Fe 2 5005.52A: [2.59106731e-09] erg/(cm**2*s) @ [5005.52] Å
O 3 5006.84A: [1.00215984e-05] erg/(cm**2*s) @ [5006.84] Å
Fe 2 5018.44A: [5.77409211e-11] erg/(cm**2*s) @ [5018.44] Å
Fe 2 5020.25A: [7.54151084e-10] erg/(cm**2*s) @ [5020.25] Å
Fe 2 5049.30A: [3.0613796e-09] erg/(cm**2*s) @ [5049.3] Å
Fe 2 5072.41A: [7.38275772e-09] erg/(cm**2*s) @ [5072.41] Å
Fe 2 5111.64A: [4.30437236e-09] erg/(cm**2*s) @ [5111.64] Å
Fe 2 5158.01A: [1.22574126e-09] erg/(cm**2*s) @ [5158.01] Å
Fe 2 5158.79A: [3.75429031e-08] erg/(cm**2*s) @ [5158.79] Å
Fe 2 5169.03A: [3.31720314e-10] erg/(cm**2*s) @ [5169.03] Å
Fe 6 5176.04A: [3.54889424e-15] erg/(cm**2*s) @ [5176.04] Å
Fe 2 5184.80A: [1.06128678e-09] erg/(cm**2*s) @ [5184.8] Å
Fe 2 5261.63A: [1.45345124e-08] erg/(cm**2*s) @ [5261.63] Å
Fe 3 5270.40A: [1.91717749e-07] erg/(cm**2*s) @ [5270.4] Å
Fe 2 5273.36A: [8.56333052e-09] erg/(cm**2*s) @ [5273.36] Å
Fe 2 5284.10A: [1.62729329e-11] erg/(cm**2*s) @ [5284.1] Å
Fe 2 5333.66A: [7.18321044e-09] erg/(cm**2*s) @ [5333.66] Å
Fe 2 5376.47A: [3.29243175e-09] erg/(cm**2*s) @ [5376.47] Å
Fe 2 5412.67A: [6.97776522e-10] erg/(cm**2*s) @ [5412.67] Å
Fe 2 5433.15A: [2.63333595e-09] erg/(cm**2*s) @ [5433.15] Å
Fe 2 5527.36A: [5.17200982e-09] erg/(cm**2*s) @ [5527.36] Å
Fe 7 5720.71A: [0.] erg/(cm**2*s) @ [5720.71] Å
He 1 5875.61A: [1.7936646e-07] erg/(cm**2*s) @ [5875.61] Å
He 1 5875.64A: [1.0761791e-07] erg/(cm**2*s) @ [5875.64] Å
Fe 7 6086.97A: [0.] erg/(cm**2*s) @ [6086.97] Å
O 1 6300.30A: [3.67633909e-08] erg/(cm**2*s) @ [6300.3] Å
O 1 6363.78A: [1.17590359e-08] erg/(cm**2*s) @ [6363.78] Å
Fe 2 6516.08A: [1.45482743e-12] erg/(cm**2*s) @ [6516.08] Å
N 2 6548.05A: [1.93570688e-07] erg/(cm**2*s) @ [6548.05] Å
H 1 6562.80A: [7.70173846e-06] erg/(cm**2*s) @ [6562.8] Å
N 2 6583.45A: [5.71118066e-07] erg/(cm**2*s) @ [6583.45] Å
Ni 2 6666.80A: [1.15324564e-09] erg/(cm**2*s) @ [6666.8] Å
He 1 6678.15A: [8.89214118e-08] erg/(cm**2*s) @ [6678.15] Å
S 2 6716.44A: [2.45575853e-07] erg/(cm**2*s) @ [6716.44] Å
S 2 6730.82A: [3.02269889e-07] erg/(cm**2*s) @ [6730.82] Å
Ar 3 7135.79A: [3.12198103e-07] erg/(cm**2*s) @ [7135.79] Å
Fe 2 7155.17A: [4.16213349e-08] erg/(cm**2*s) @ [7155.17] Å
Fe 2 7172.00A: [7.80691855e-09] erg/(cm**2*s) @ [7172.] Å
Ca 2 7291.47A: [5.86600182e-08] erg/(cm**2*s) @ [7291.47] Å
Ca 2 7323.89A: [3.86669539e-08] erg/(cm**2*s) @ [7323.89] Å
Ni 2 7377.83A: [3.07361785e-08] erg/(cm**2*s) @ [7377.83] Å
Fe 2 7388.17A: [5.94169217e-09] erg/(cm**2*s) @ [7388.17] Å
Ni 2 7411.61A: [1.9327579e-09] erg/(cm**2*s) @ [7411.61] Å
Fe 2 7452.56A: [1.34689496e-08] erg/(cm**2*s) @ [7452.56] Å
Ar 3 7751.11A: [7.48695482e-08] erg/(cm**2*s) @ [7751.11] Å
O 1 8446.25A: [2.13122389e-10] erg/(cm**2*s) @ [8446.25] Å
O 1 8446.36A: [5.04078045e-10] erg/(cm**2*s) @ [8446.36] Å
O 1 8446.76A: [4.5088346e-10] erg/(cm**2*s) @ [8446.76] Å
Cl 2 8578.70A: [8.66507014e-09] erg/(cm**2*s) @ [8578.7] Å
Fe 2 8616.95A: [1.94489785e-08] erg/(cm**2*s) @ [8616.95] Å
Fe 2 8891.93A: [9.02730251e-09] erg/(cm**2*s) @ [8891.93] Å
Fe 2 9051.95A: [2.81622152e-08] erg/(cm**2*s) @ [9051.95] Å
S 3 9068.62A: [9.00567669e-07] erg/(cm**2*s) @ [9068.62] Å
Fe 2 9226.63A: [1.03214108e-08] erg/(cm**2*s) @ [9226.63] Å
H 1 9229.02A: [7.37882575e-08] erg/(cm**2*s) @ [9229.02] Å
Fe 2 9267.56A: [4.15584235e-09] erg/(cm**2*s) @ [9267.56] Å
Fe 2 9399.04A: [1.06988796e-08] erg/(cm**2*s) @ [9399.04] Å
Fe 2 9470.94A: [4.21093902e-09] erg/(cm**2*s) @ [9470.94] Å
S 3 9530.62A: [2.27464235e-06] erg/(cm**2*s) @ [9530.62] Å
H 1 9545.97A: [1.07037945e-07] erg/(cm**2*s) @ [9545.97] Å
H 1 1.00494m: [1.61643123e-07] erg/(cm**2*s) @ [10049.4] Å
S 2 1.02867m: [8.35051286e-09] erg/(cm**2*s) @ [10286.7] Å
S 2 1.03205m: [1.13643648e-08] erg/(cm**2*s) @ [10320.5] Å
S 2 1.03364m: [8.20503011e-09] erg/(cm**2*s) @ [10336.4] Å
He 1 1.08291m: [8.51762141e-08] erg/(cm**2*s) @ [10829.1] Å
He 1 1.08303m: [5.32289901e-07] erg/(cm**2*s) @ [10830.3] Å
H 1 1.09381m: [2.6269465e-07] erg/(cm**2*s) @ [10938.1] Å
O 1 1.12863m: [8.76484378e-11] erg/(cm**2*s) @ [11286.3] Å
O 1 1.12864m: [7.75727748e-11] erg/(cm**2*s) @ [11286.4] Å
O 1 1.12869m: [5.05110186e-11] erg/(cm**2*s) @ [11286.9] Å
O 1 1.12870m: [2.92400654e-11] erg/(cm**2*s) @ [11287.] Å
O 1 1.12873m: [1.03420536e-10] erg/(cm**2*s) @ [11287.3] Å
Ni 2 1.19102m: [6.69908712e-09] erg/(cm**2*s) @ [11910.2] Å
Fe 2 1.25668m: [2.08217261e-08] erg/(cm**2*s) @ [12566.8] Å
Fe 2 1.27877m: [8.88673472e-10] erg/(cm**2*s) @ [12787.7] Å
H 1 1.28181m: [4.75920253e-07] erg/(cm**2*s) @ [12818.1] Å
Fe 2 1.29427m: [1.79730167e-09] erg/(cm**2*s) @ [12942.7] Å
Fe 2 1.32055m: [6.66320793e-09] erg/(cm**2*s) @ [13205.5] Å
Fe 2 1.32777m: [1.17978641e-09] erg/(cm**2*s) @ [13277.7] Å
Fe 2 1.37181m: [3.8137613e-09] erg/(cm**2*s) @ [13718.1] Å
Fe 2 1.53348m: [4.25094677e-08] erg/(cm**2*s) @ [15334.8] Å
Fe 2 1.59948m: [1.7068173e-08] erg/(cm**2*s) @ [15994.8] Å
Fe 2 1.64355m: [3.19555794e-07] erg/(cm**2*s) @ [16435.5] Å
Fe 2 1.66377m: [8.83217869e-09] erg/(cm**2*s) @ [16637.7] Å
Fe 2 1.67688m: [1.60586408e-08] erg/(cm**2*s) @ [16768.8] Å
Fe 2 1.71113m: [3.6002064e-09] erg/(cm**2*s) @ [17111.3] Å
Fe 2 1.74494m: [4.21511673e-09] erg/(cm**2*s) @ [17449.4] Å
Fe 2 1.79711m: [8.02366818e-09] erg/(cm**2*s) @ [17971.1] Å
Fe 2 1.80002m: [1.55602641e-08] erg/(cm**2*s) @ [18000.2] Å
Fe 2 1.80940m: [8.06298784e-08] erg/(cm**2*s) @ [18094.] Å
H 1 1.87510m: [9.69770272e-07] erg/(cm**2*s) @ [18751.] Å
Fe 2 1.89541m: [2.86758283e-09] erg/(cm**2*s) @ [18954.1] Å
Ni 2 1.93877m: [6.53443605e-09] erg/(cm**2*s) @ [19387.7] Å
Fe 2 1.95361m: [9.94050162e-09] erg/(cm**2*s) @ [19536.1] Å
Si 6 1.96247m: [0.] erg/(cm**2*s) @ [19624.7] Å
H 1 2.16553m: [8.35297033e-08] erg/(cm**2*s) @ [21655.3] Å
Si 7 2.48071m: [0.] erg/(cm**2*s) @ [24807.1] Å
[32]:
(<Figure size 800x600 with 1 Axes>,
<Axes: xlabel='$ \\lambda / \\AA$', ylabel='$L / $erg s$^{-1}$'>)
