The SKIRT project
advanced radiative transfer for astrophysics
ThemisDustMix Class Reference

#include <ThemisDustMix.hpp>

Inheritance diagram for ThemisDustMix:

Public Member Functions

int numHydrocarbonSizes () const
int numSilicateSizes () const
Public Member Functions inherited from MultiGrainDustMix
Array emissivity (const Array &Jv) const override
bool hasStochasticDustEmission () const override
int numPopulations () const override
const GrainPopulationpopulation (int c) const
double populationBulkDensity (int c) const override
string populationGrainType (int c) const override
double populationMass (int c) const override
double populationNormalization (int c) const
const GrainSizeDistributionpopulationSizeDistribution (int c) const override
Range populationSizeRange (int c) const override
double totalMass () const override
Public Member Functions inherited from DustMix
double asymmpar (double lambda) const override
Array emissionSpectrum (const MaterialState *state, const Array &Jv) const override
DisjointWavelengthGridemissionWavelengthGrid () 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
bool 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 ArraysectionsAbs (double lambda) const override
const ArraysectionsAbspol (double lambda) const override
double sectionSca (double lambda) const override
vector< StateVariablespecificStateVariableInfo () const override
const ArraythetaGrid () const override
Public Member Functions inherited from MaterialMix
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 void initializeSpecificState (MaterialState *state, double metallicity, double temperature, const Array &params) 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 vector< SnapshotParameterparameterInfo () const
virtual bool scatteringEmulatesSecondaryEmission () 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 ItemgetItemProperty (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)
Itemoperator= (const Item &)=delete
Itemparent () 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
Public Member Functions inherited from MultiGrainPopulationInterface
virtual ~MultiGrainPopulationInterface ()

Protected Member Functions

 ThemisDustMix ()
void setupSelfBefore () override
Protected Member Functions inherited from MultiGrainDustMix
 MultiGrainDustMix ()
void addPopulation (const GrainPopulation *population)
void addPopulation (GrainComposition *composition, GrainSizeDistribution *sizeDistribution, int numSizes, GrainPopulation::NormalizationType normType, double normValue)
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
size_t initializeExtraProperties (const Array &lambdav) override
Protected Member Functions inherited from DustMix
 DustMix ()
void informAvailableWavelengthRange (Range available)
void setupSelfAfter () override
Protected Member Functions inherited from MaterialMix
 MaterialMix ()
Configurationconfig () const
Randomrandom () const
Protected Member Functions inherited from SimulationItem
 SimulationItem ()
virtual bool offersInterface (const std::type_info &interfaceTypeInfo) const
Protected Member Functions inherited from Item
 Item ()
Protected Member Functions inherited from MultiGrainPopulationInterface
 MultiGrainPopulationInterface ()

Private Types

using BaseType
using ItemType

Private Attributes

int _numHydrocarbonSizes
int _numSilicateSizes

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 }

Detailed Description

The ThemisDustMix class represents the THEMIS model for dust in the diffuse interstellar medium described by Jones et al. 2017 (A&A, 602, A46) and the references therein. In this model, there are two families of dust particles: amorphous silicate and amorphous hydrocarbon. For the silicate, it is assumed that 50% of the mass is amorphous enstatite, and that the remaining half is amorphous forsterite, and the size distribution is a lognormal distribution (the same for both populations). For the amorphous hydrocarbon, the size distribution is a combination of a lognormal and a power-law distribution.

For each of the populations, the user can configure the number of grain size bins used to discretize the thermal emission calculations for the dust mix.

The following is an extract from the GRAIN_J17.DAT file specifying the THEMIS dust model for the DustEM code, described by Compiègne et al. 2011 (AA, 525, A103):

# grain type, nsize, type keywords, Mdust/MH, rho, amin, amax, alpha/a0 [, at, ac, gamma (ED)] [, au, zeta, eta (CV)]
# cgs units
aOLM5  100 logn      0.255E-02  2.190E+00  1.00E-07   4900.0E-07     8.00E-07    1.00E+00
aPyM5  100 logn      0.255E-02  2.190E+00  1.00E-07   4900.0E-07     8.00E-07    1.00E+00
CM20   100 logn      0.600E-03  1.510E+00  0.50E-07   4900.0E-07     7.00E-07    1.00E+00
CM20   100 plaw-ed   0.170E-02  1.600E+00  0.40E-07   4900.0E-07    -5.00E-00   10.00E-07   50.0E-07  1.000E+00

The aOlM5 (olivine) grain type represents the amorphous silicates with forsterite-normative composition, aPyM5 (pyroxene) represents the amorphous silicates with enstatite-normative composition, and CM20 represents the amorphous carbonaceous dust grains in the THEMIS model.

This class uses grain compositions obtained through the DustEmGrainComposition class, respectively selecting the grain types aOLM5, aPyM5, and CM20. It further uses the size distributions defined in the ModifiedPowerLawGrainSizeDistribution and LogNormalGrainSizeDistribution classes.

Constructor & Destructor Documentation

◆ ThemisDustMix()

ThemisDustMix::ThemisDustMix ( )
inlineprotected

Default constructor for concrete Item subclass ThemisDustMix: "a THEMIS (Jones et al. 2017) dust mix".

Member Function Documentation

◆ numHydrocarbonSizes()

int ThemisDustMix::numHydrocarbonSizes ( ) const
inline

This function returns the value of the discoverable integer property numHydrocarbonSizes: "the number of grain size bins for each of the hydrocarbon populations".

The minimum value for this property is "1".

The default value for this property is given by the conditional value expression "5".

◆ numSilicateSizes()

int ThemisDustMix::numSilicateSizes ( ) const
inline

This function returns the value of the discoverable integer property numSilicateSizes: "the number of grain size bins for each of the silicate populations".

The minimum value for this property is "1".

The default value for this property is given by the conditional value expression "5".

◆ setupSelfBefore()

void ThemisDustMix::setupSelfBefore ( )
overrideprotectedvirtual

This function adds the relevant grain populations to the dust mix

Reimplemented from MaterialMix.


The documentation for this class was generated from the following file: