The SKIRT project
advanced radiative transfer for astrophysics
WeingartnerDraineDustMix Class Reference

#include <WeingartnerDraineDustMix.hpp>

Inheritance diagram for WeingartnerDraineDustMix:

Public Types

enum class  Environment : int { MilkyWay , LMC , SMC }
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

Environment environment () const
int numGraphiteSizes () const
int numPAHSizes () 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

 WeingartnerDraineDustMix ()
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

Environment _environment
int _numGraphiteSizes
int _numPAHSizes
int _numSilicateSizes

Friends

class ItemRegistry

Detailed Description

The WeingartnerDraineDustMix class represents realistic dust mixtures consisting of populations of silicate, graphite, and PAH dust grains. It is assumed that 50% of the PAHs are neutral and 50% are ionized. The size distributions of each of these grains are fitted in such a way that they can reproduce the extinction curve of the Milky Way, the LMC, or the SMC, respectively. For details we refer to Weingartner & Draine (2001, ApJ, 548, 296) and Li & Draine (2001, ApJ, 554, 778).

As noted at the end of section 3 in Weingartner & Draine (2001), the absence of PAHs in the SMC dust model might be a special attribute of the specific sight line, so it may be unrealistic.

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 class uses the grain compositions defined in the DraineSilicateGrainComposition, DraineGraphiteGrainComposition, DraineNeutralPAHGrainComposition, and DraineIonizedPAHGrainComposition classes.

The grain size distributions are implementated internally in the class because they are very specific to this particular dust mix. They are given as complicated analytical formulas with predefined constant parameters depending on the selected environment (Milky Way, LMC or SMC). For the graphite and silicate populations, the distribution is a power-law function with a curvature and an exponential cutoff. For the PAH populations, it is the sum of two log-normal distributions, cut off at some upper grain size. The exact forms of the grain size distribution functions can be found in Weingartner & Draine (2001, ApJ, 548, 296).

Member Enumeration Documentation

◆ Environment

enum class WeingartnerDraineDustMix::Environment : int
strong

The enumeration type indicating the environment for the Weingartner-Draine dust.

MilkyWay : "the Milky Way".

LMC : "the Large Magellanic Cloud".

SMC : "the Small Magellanic Cloud".

Constructor & Destructor Documentation

◆ WeingartnerDraineDustMix()

WeingartnerDraineDustMix::WeingartnerDraineDustMix ( )
inlineprotected

Default constructor for concrete Item subclass WeingartnerDraineDustMix: "a Weingartner and Draine (2001) dust mix".

Member Function Documentation

◆ environment()

Environment WeingartnerDraineDustMix::environment ( ) const
inline

This function returns the value of the discoverable Environment enumeration property environment: "the environment determining the dust model".

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

◆ numGraphiteSizes()

int WeingartnerDraineDustMix::numGraphiteSizes ( ) const
inline

This function returns the value of the discoverable integer property numGraphiteSizes: "the number of graphite grain size bins".

The minimum value for this property is "1".

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

◆ numPAHSizes()

int WeingartnerDraineDustMix::numPAHSizes ( ) const
inline

This function returns the value of the discoverable integer property numPAHSizes: "the number of neutral and ionized PAH size bins (each)".

The minimum value for this property is "1".

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

◆ numSilicateSizes()

int WeingartnerDraineDustMix::numSilicateSizes ( ) const
inline

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

The minimum value for this property is "1".

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

◆ setupSelfBefore()

void WeingartnerDraineDustMix::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: