SKIRT / SKIRT9

SKIRT version 9 -- advanced radiative transfer in dusty systems
http://www.skirt.ugent.be
GNU Affero General Public License v3.0
38 stars 31 forks source link

Introduce the line emission spectra for normalized sources #88

Closed petercamps closed 3 years ago

petercamps commented 3 years ago

Description This update introduces the concept of line emission spectra for point sources and sources using a built-in geometry. The emission spectrum of these sources is constant across space and is defined through a single SED instance. This update does not affect imported sources where the emission spectrum can vary across space.

In this implementation, line spectra are defined by a list of wavelengths and the corresponding luminosities (integrated over the line profile). All photon packets are emitted exactly at the line wavelengths, i.e. the lines have zero width. As a result, the specific luminosity of the line spectrum is numerically undefined. This means that the line spectrum must be normalized by defining the integrated luminosity over a wavelength range that includes one or more of the line wavelengths. Also, line spectra are always fully distinct from regular "continuum" spectra. To combine line and continuum spectra, one can configure two sources with the same spatial geometry.

Specifically:

Motivation We implemented these features as a way to experiment with line emission spectra in preparation for a (much) more involved implementation of line emission for imported sources. Also, having these features available for built-in geometries allows configuring toy models for testing and benchmarking.

Tests Test cases were added for the new features. The existing test cases still run with the following exception. The LyaSEDDecorator and LyaSEDFamilyDecorator classes have spectra properties that now should be of type ContSED rather than the more general SED. This requires a small change in the ski file. The PTS ski file upgrade procedure handles this change (see the corresponding PTS9 pull request).