#include <ThemisDustMix.hpp>

Inheritance diagram for ThemisDustMix:

[legend]

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 GrainPopulation *	population (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 GrainSizeDistribution *	populationSizeDistribution (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

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 &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 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

Public Member Functions inherited from MultiGrainPopulationInterface
virtual	~MultiGrainPopulationInterface ()

virtual int	numPopulations () const =0

virtual double	populationBulkDensity (int c) const =0

virtual string	populationGrainType (int c) const =0

virtual double	populationMass (int c) const =0

virtual const GrainSizeDistribution *	populationSizeDistribution (int c) const =0

virtual Range	populationSizeRange (int c) const =0

virtual double	totalMass () const =0

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

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

Protected Member Functions inherited from MultiGrainPopulationInterface
	MultiGrainPopulationInterface ()

Private Types
using	BaseType = MultiGrainDustMix

using	ItemType = ThemisDustMix

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

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

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:

ThemisDustMix.hpp

Public Member Functions

Protected Member Functions

Private Types

Private Attributes

Friends