 |
The SKIRT project
advanced radiative transfer for astrophysics
|
|
The SKIRT project
advanced radiative transfer for astrophysics
|
#include <CastelliKuruczSEDFamily.hpp>
Public Member Functions | |
| CastelliKuruczSEDFamily (SimulationItem *parent) | |
| double | cdf (Array &lambdav, Array &pv, Array &Pv, const Range &wavelengthRange, const Array ¶meters) const override |
| Range | intrinsicWavelengthRange () const override |
| vector< SnapshotParameter > | parameterInfo () const override |
| double | specificLuminosity (double wavelength, const Array ¶meters) const override |
| Public Member Functions inherited from SimulationItem | |
| template<class T> | |
| T * | find (bool setup=true) const |
| template<class T> | |
| T * | interface (int levels=-999999, bool setup=true) const |
| virtual string | itemName () const |
| void | setup () |
| string | typeAndName () const |
| Public Member Functions inherited from Item | |
| Item (const Item &)=delete | |
| virtual | ~Item () |
| void | addChild (Item *child) |
| const vector< Item * > & | children () const |
| virtual void | clearItemListProperty (const PropertyDef *property) |
| void | destroyChild (Item *child) |
| virtual bool | getBoolProperty (const PropertyDef *property) const |
| virtual vector< double > | getDoubleListProperty (const PropertyDef *property) const |
| virtual double | getDoubleProperty (const PropertyDef *property) const |
| virtual string | getEnumProperty (const PropertyDef *property) const |
| virtual int | getIntProperty (const PropertyDef *property) const |
| virtual vector< Item * > | getItemListProperty (const PropertyDef *property) const |
| virtual Item * | getItemProperty (const PropertyDef *property) const |
| virtual string | getStringProperty (const PropertyDef *property) const |
| int | getUtilityProperty (string name) const |
| virtual void | insertIntoItemListProperty (const PropertyDef *property, int index, Item *item) |
| Item & | operator= (const Item &)=delete |
| Item * | parent () const |
| virtual void | removeFromItemListProperty (const PropertyDef *property, int index) |
| virtual void | setBoolProperty (const PropertyDef *property, bool value) |
| virtual void | setDoubleListProperty (const PropertyDef *property, vector< double > value) |
| virtual void | setDoubleProperty (const PropertyDef *property, double value) |
| virtual void | setEnumProperty (const PropertyDef *property, string value) |
| virtual void | setIntProperty (const PropertyDef *property, int value) |
| virtual void | setItemProperty (const PropertyDef *property, Item *item) |
| virtual void | setStringProperty (const PropertyDef *property, string value) |
| void | setUtilityProperty (string name, int value) |
| virtual string | type () const |
Protected Member Functions | |
| CastelliKuruczSEDFamily () | |
| void | setupSelfBefore () override |
| Protected Member Functions inherited from SEDFamily | |
| SEDFamily () | |
| Protected Member Functions inherited from SimulationItem | |
| SimulationItem () | |
| virtual bool | offersInterface (const std::type_info &interfaceTypeInfo) const |
| virtual void | setupSelfAfter () |
| Protected Member Functions inherited from Item | |
| Item () | |
Private Types | |
| using | BaseType |
| using | ItemType |
Private Attributes | |
| StoredTable< 4 > | _table |
Friends | |
| class | ItemRegistry |
An instance of the CastelliKuruczSEDFamily class represents the Castelli-Kurucz atlas of stellar atmosphere models for a wide range of metallicities, effective temperatures, and gravities. These models were computed by Fiorella Castelli (Castelli & Kurucz 2003, IAU Symposium 210, Modelling of Stellar Atmospheres, Uppsala, Sweden, eds. N.E. Piskunov, W.W. Weiss. and D.F. Gray, 2003, ASP-S210) and include improvements upon those previously provided by by Kurucz (1990).
The luminosity values given by the models are normalized for a stellar surface equal to the unit sphere, so these values must be multipled by the stellar surface, i.e. \(4\pi\,R^2\).
The parameter space is discretized on 8 metallicities, 76 temperatures, and 11 gravity values for a total of 6688 possible models. From these possibilities, 3808 models are actually available and 2880 have not been calculated. The parameter ranges for available models are summarized in the table below, with in all cases \( 6.325 \times 10^{-5} \le Z \le 6.325 \times 10^{-2} \).
| Effective temperature (K) | Gravity ( \(\log_{10} (g/\mathrm{cm}\,\mathrm{s}^{-2})\)) | |
|---|---|---|
| 3500 <= Teff <= 6000 | --> | 0.0 <= logg <= 5.0 |
| 6000 < Teff <= 7500 | --> | 0.5 <= logg <= 5.0 |
| 7500 < Teff <= 8250 | --> | 1.0 <= logg <= 5.0 |
| 8250 < Teff <= 9000 | --> | 1.5 <= logg <= 5.0 |
| 9000 < Teff <= 11750 | --> | 2.0 <= logg <= 5.0 |
| 11750 < Teff <= 19000 | --> | 2.5 <= logg <= 5.0 |
| 19000 < Teff <= 26000 | --> | 3.0 <= logg <= 5.0 |
| 26000 < Teff <= 31000 | --> | 3.5 <= logg <= 5.0 |
| 31000 < Teff <= 39000 | --> | 4.0 <= logg <= 5.0 |
| 39000 < Teff <= 49000 | --> | 4.5 <= logg <= 5.0 |
| 49000 < Teff <= 50000 | --> | logg = 5.0 |
The following table lists the properties of solar metallicity stars ( \(Z=0.02\)) of different spectral types and luminosity classes as suggested by Castelli and Kurucz based on other sources.
| Stellar type | Temperature (K) | Gravity ( \(\log_{10} (g/\mathrm{cm}\,\mathrm{s}^{-2})\)) |
|---|---|---|
| O3V | 44852 | +3.92 |
| O5V | 40862 | +3.92 |
| O5.5V | 39865 | +3.92 |
| O6V | 38867 | +3.92 |
| O6.5V | 37870 | +3.92 |
| O7V | 36872 | +3.92 |
| O7.5V | 35874 | +3.92 |
| O8V | 34877 | +3.92 |
| O8.5 | 33879 | +3.92 |
| O9V | 32882 | +3.92 |
| O9.5 | 31884 | +3.92 |
| B0V | 30000 | +3.90 |
| B1V | 25400 | +3.90 |
| B3V | 18700 | +3.94 |
| B5V | 15400 | +4.04 |
| B8V | 11900 | +4.04 |
| A0V | 9520 | +4.14 |
| A1V | 9230 | +4.10 |
| A3V | 8270 | +4.20 |
| A5V | 8200 | +4.29 |
| F0V | 7200 | +4.34 |
| F2V | 6890 | +4.34 |
| F5V | 6440 | +4.34 |
| F8V | 6200 | +4.40 |
| G0V | 6030 | +4.39 |
| G2V | 5860 | +4.40 |
| G8V | 5570 | +4.50 |
| K0V | 5250 | +4.49 |
| K2V | 4780 | +4.5 |
| K4V | 4560 | +4.5 |
| K5V | 4350 | +4.54 |
| K7V | 4060 | +4.5 |
| M0V | 3850 | +4.59 |
| M2V | 3580 | +4.64 |
| M6V | 3050 | +5.00 |
| B0III | 29000 | +3.34 |
| B5III | 15000 | +3.49 |
| G0III | 5850 | +2.94 |
| G5III | 5150 | +2.54 |
| K0III | 4750 | +2.14 |
| K5III | 3950 | +1.74 |
| M0III | 3800 | +1.34 |
| BOI | 26000 | +2.84 |
| B5I | 13600 | +2.44 |
| AOI | 9730 | +2.14 |
| A5I | 8510 | +2.04 |
| F0I | 7700 | +1.74 |
| F5I | 6900 | +1.44 |
| G0I | 5550 | +1.34 |
| G5I | 4850 | +1.14 |
| K0I | 4420 | +0.94 |
| K5I | 3850 | +0.00 |
| M0I | 3650 | -0.10 |
| M2I | 3600 | -0.10 |
The SEDs are tabulated over a wavelength range from 0.009 \(\mu\mathrm{m}\) to 160 \(\mu\mathrm{m}\) with the spectral resolution shown in the figure below.
The data for the models were downloaded from http://www.stsci.edu/hst/observatory/crds/castelli_kurucz_atlas.html and converted to SKIRT stored table format for inclusion as a resource file. The stored table is opened during setup, and it is subsequently interpolated to the desired parameters and wavelength grid when needed.
When imported from a text column file, the properties of the star-forming region represented by this SED family must appear in the following order, and have the specified default units unless these units are overridden by column header info:
\[ R\,(\mathrm{km}) \quad Z\,(\mathrm{dimless}) \quad T_\mathrm{eff}\,(\mathrm{K}) \quad g\,(\mathrm{m}\,\mathrm{s}^{-2}) \]
where \(R\) is the stellar radius, \(T_\mathrm{eff}\) is the effective stellar surface temperature, \(Z\) is the stellar metallicity, and \(g\) is the gravity at the stellar surface.
|
inlineprotected |
Default constructor for concrete Item subclass CastelliKuruczSEDFamily: "a Castelli-Kurucz SED family for stellar atmospheres".
|
explicit |
This constructor can be invoked programmatically by classes that use a hard-coded SED family (as opposed to selected through the ski file). Before the constructor returns, the newly created object is hooked up as a child to the specified parent in the simulation hierarchy (so it will automatically be deleted), and its setup() function has been called.
|
overridevirtual |
This function constructs both the normalized probability density function (pdf) and the corresponding normalized cumulative distribution function (cdf) for the SED with the specified parameters over the specified wavelength range. The function returns the normalization factor. The number and type of parameters must match the information returned by the parameterInfo() function; if not the behavior is undefined.
Implements SEDFamily.
|
overridevirtual |
This function returns the intrinsic wavelength range of the SED family. It retrieves this range from the underlying stored table.
Implements SEDFamily.
|
overridevirtual |
This function returns the number and type of parameters used by this particular SED family as a list of SnapshotParameter objects. Each of these objects specifies unit information and a human-readable descripton for the parameter.
Implements SEDFamily.
|
overrideprotectedvirtual |
This function opens the appropriate resource file (in SKIRT stored table format).
Reimplemented from SimulationItem.
|
overridevirtual |
This function returns the specific luminosity \(L_\lambda\) (i.e. radiative power per unit of wavelength) for the SED with the specified parameters at the specified wavelength, or zero if the wavelength is outside of the SED's intrinsic wavelength range. The number and type of parameters must match the information returned by the parameterInfo() function; if not the behavior is undefined.
Implements SEDFamily.