 |
The SKIRT project
advanced radiative transfer for astrophysics
|
|
The SKIRT project
advanced radiative transfer for astrophysics
|
#include <MeanInterstellarDustMix.hpp>
Protected Member Functions | |
| MeanInterstellarDustMix () | |
| string | resourceNameForOpticalProps () const override |
 Protected Member Functions inherited from SingleGrainDustMix | |
| SingleGrainDustMix () | |
| double | getOpticalProperties (const Array &lambdav, const Array &thetav, Array &sigmaabsv, Array &sigmascav, Array &asymmparv, Table< 2 > &S11vv, Table< 2 > &S12vv, Table< 2 > &S33vv, Table< 2 > &S34vv, ArrayTable< 2 > &sigmaabsvv, ArrayTable< 2 > &sigmaabspolvv) override |
| virtual string | resourceNameForMuellerMatrix () const |
| virtual string | resourceNameForOpticalProps () const =0 |
| Protected Member Functions inherited from DustMix | |
| DustMix () | |
| virtual double | getOpticalProperties (const Array &lambdav, const Array &thetav, Array &sigmaabsv, Array &sigmascav, Array &asymmparv, Table< 2 > &S11vv, Table< 2 > &S12vv, Table< 2 > &S33vv, Table< 2 > &S34vv, ArrayTable< 2 > &sigmaabsvv, ArrayTable< 2 > &sigmaabspolvv)=0 |
| void | informAvailableWavelengthRange (Range available) |
| virtual size_t | initializeExtraProperties (const Array &lambdav) |
| void | setupSelfAfter () override |
| Protected Member Functions inherited from MaterialMix | |
| MaterialMix () | |
| Configuration * | config () const |
| Random * | random () const |
| void | setupSelfBefore () override |
| Protected Member Functions inherited from SimulationItem | |
| SimulationItem () | |
| virtual bool | offersInterface (const std::type_info &interfaceTypeInfo) const |
| virtual void | setupSelfAfter () |
| virtual void | setupSelfBefore () |
| Protected Member Functions inherited from Item | |
| Item () | |
Private Types | |
| using | BaseType = SingleGrainDustMix |
| using | ItemType = MeanInterstellarDustMix |
Friends | |
| class | ItemRegistry |
Additional Inherited Members | |
| Public Types inherited from DustMix | |
| enum class | ScatteringMode { HenyeyGreenstein , MaterialPhaseFunction , SphericalPolarization , SpheroidalPolarization } |
| Public Types inherited from MaterialMix | |
| enum class | DynamicStateType { None , Primary , Secondary , PrimaryIfMergedIterations } |
| enum class | MaterialType { Dust , Electrons , Gas } |
| Public Member Functions inherited from DustMix | |
| double | asymmpar (double lambda) const override |
| Array | emissionSpectrum (const MaterialState *state, const Array &Jv) const override |
| DisjointWavelengthGrid * | emissionWavelengthGrid () const override |
| Array | emissivity (const Array &Jv) const override |
| bool | hasContinuumEmission () const override |
| double | indicativeTemperature (const MaterialState *state, const Array &Jv) const override |
| double | mass () const override |
| MaterialType | materialType () const override |
| double | opacityAbs (double lambda, const MaterialState *state, const PhotonPacket *pp) const override |
| double | opacityExt (double lambda, const MaterialState *state, const PhotonPacket *pp) const override |
| double | opacitySca (double lambda, const MaterialState *state, const PhotonPacket *pp) const override |
| void | peeloffScattering (double &I, double &Q, double &U, double &V, double &lambda, Direction bfkobs, Direction bfky, const MaterialState *state, const PhotonPacket *pp) const override |
| void | performScattering (double lambda, const MaterialState *state, PhotonPacket *pp) const override |
| virtual ScatteringMode | scatteringMode () const |
| double | sectionAbs (double lambda) const override |
| double | sectionExt (double lambda) const override |
| const Array & | sectionsAbs (double lambda) const override |
| const Array & | sectionsAbspol (double lambda) const override |
| double | sectionSca (double lambda) const override |
| vector< StateVariable > | specificStateVariableInfo () const override |
| const Array & | thetaGrid () const override |
| Public Member Functions inherited from MaterialMix | |
| virtual double | asymmpar (double lambda) const |
| virtual Array | emissionSpectrum (const MaterialState *state, const Array &Jv) const |
| virtual DisjointWavelengthGrid * | emissionWavelengthGrid () const |
| virtual Array | emissivity (const Array &Jv) const |
| virtual bool | hasContinuumEmission () const |
| virtual DynamicStateType | hasDynamicMediumState () const |
| virtual bool | hasExtraSpecificState () const |
| virtual bool | hasLineEmission () const |
| virtual bool | hasNegativeExtinction () const |
| virtual bool | hasPolarizedAbsorption () const |
| virtual bool | hasPolarizedEmission () const |
| virtual bool | hasPolarizedScattering () const |
| virtual bool | hasResonantScattering () const |
| virtual bool | hasScatteringDispersion () const |
| virtual bool | hasStochasticDustEmission () const |
| virtual double | indicativeTemperature (const MaterialState *state, const Array &Jv) const |
| virtual void | initializeSpecificState (MaterialState *state, double metallicity, double temperature, const Array ¶ms) const |
| bool | isDust () const |
| bool | isElectrons () const |
| bool | isGas () const |
| virtual bool | isSpecificStateConverged (int numCells, int numUpdated, int numNotConverged, MaterialState *currentAggregate, MaterialState *previousAggregate) const |
| virtual Array | lineEmissionCenters () const |
| virtual Array | lineEmissionMasses () const |
| virtual Array | lineEmissionSpectrum (const MaterialState *state, const Array &Jv) const |
| virtual double | mass () const =0 |
| virtual MaterialType | materialType () const =0 |
| virtual double | opacityAbs (double lambda, const MaterialState *state, const PhotonPacket *pp) const =0 |
| virtual double | opacityExt (double lambda, const MaterialState *state, const PhotonPacket *pp) const =0 |
| virtual double | opacitySca (double lambda, const MaterialState *state, const PhotonPacket *pp) const =0 |
| virtual vector< SnapshotParameter > | parameterInfo () const |
| virtual void | peeloffScattering (double &I, double &Q, double &U, double &V, double &lambda, Direction bfkobs, Direction bfky, const MaterialState *state, const PhotonPacket *pp) const =0 |
| virtual void | performScattering (double lambda, const MaterialState *state, PhotonPacket *pp) const =0 |
| virtual double | sectionAbs (double lambda) const =0 |
| virtual double | sectionExt (double lambda) const =0 |
| virtual const Array & | sectionsAbs (double lambda) const |
| virtual const Array & | sectionsAbspol (double lambda) const |
| virtual double | sectionSca (double lambda) const =0 |
| virtual vector< StateVariable > | specificStateVariableInfo () const =0 |
| virtual const Array & | thetaGrid () const |
| virtual UpdateStatus | updateSpecificState (MaterialState *state, const Array &Jv) const |
| 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 |
The MeanInterstellarDustMix class represents a population of identical dust grains with properties approximating those of a mixture that is appropriate for the typical interstellar dust medium. The model consists of a mixture of carbonaceous grains and amorphous silicate grains. Carbonaceous grains are PAH-like when small, and graphite-like when large. (see Li & Draine 2001, ApJ, 554, 778). Size distributions are taken from Weingartner & Draine (2001, ApJ, 548, 296), more specifically, they refer to the case A model for \(R_V=3.1\), renormalized following Draine (2003, ApJ, 598, 1017). Specifically, the grain abundances relative to H have been reduced by a factor 0.93 relative to the \(R_V=3.1\) size distribution with \(\text{[C/H]}_{\text{PAH}} = 60~{\text{ppm}}\) in Weingartner & Draine (2001, ApJ, 548, 296). The PAH C abundance relative to H is assumed to be \(\text{[C/H]}_{\text{PAH}} = 0.93 \times 60~{\text{ppm}} = 55.8~{\text{ppm}}\).
The dust mass per hydrogen nucleon is obtained from information in the header of the data file. The optical data were downloaded from Bruce Draine's home page https://www.astro.princeton.edu/~draine/ and more specifically from ftp://ftp.astro.princeton.edu/draine/dust/mix/kext_albedo_WD_MW_3.1_60_D03.all
Extreme forward scattering for X-ray wavelengths
The properties for this dust mix as discussed above are given down to wavelengths as short as \(10^{-4}\) micron, i.e. well into the soft X-ray wavelength range. The tabulated cross sections seem to be very accurate, however the average scattering angle cosine values in this regime are unrealistic (e.g., equal to one for the shortest wavelengths).
Draine 2003c (ApJ, 98, 1026–1037) presents experimental data for the phase function in the soft X-ray wavelength range and proposes an analytical phase function to model the observed extreme forward scattering. Rather than attempting to implement this phase function, we approximate it through the Henyey-Greenstein phase function just as we do for other wavelengths. For wavelengths shorter than 0.01 micron, the value of the HG asymmetry parameter \(g\) that causes the HG function to best approximate the Draine 2003c function for extreme forward angles can be written as \(g = 1 - \lambda/(\mu\mathrm{m})\), i.e. one minus the wavelength in micron. We thus replace the asymmetry parameter values in that wavelength range by this improved appoximation before storing the result in the SKIRT resource table.
|
inlineprotected |
Default constructor for concrete Item subclass MeanInterstellarDustMix : "a typical interstellar dust mix (mean properties)" .
|
overrideprotectedvirtual |
This function returns the name of the stored table resource tabulating the basic optical properties for this dust mix.
Implements SingleGrainDustMix.