#include <MeanInterstellarDustMix.hpp>

Inheritance diagram for MeanInterstellarDustMix:

[legend]

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

Detailed Description

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 ${[C/H]}_{PAH} = 60 ppm$ in Weingartner & Draine (2001, ApJ, 548, 296). The PAH C abundance relative to H is assumed to be ${[C/H]}_{PAH} = 0.93 \times 60 ppm = 55.8 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 - λ / (μ 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.

Constructor & Destructor Documentation

◆ MeanInterstellarDustMix()

MeanInterstellarDustMix::MeanInterstellarDustMix ( )

inlineprotected

Default constructor for concrete Item subclass MeanInterstellarDustMix : "a typical interstellar dust mix (mean properties)" .

Member Function Documentation

◆ resourceNameForOpticalProps()

string MeanInterstellarDustMix::resourceNameForOpticalProps ( ) const

overrideprotectedvirtual

This function returns the name of the stored table resource tabulating the basic optical properties for this dust mix.

Implements SingleGrainDustMix.

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

MeanInterstellarDustMix.hpp

Protected Member Functions

Private Types

Friends

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MeanInterstellarDustMix()

Member Function Documentation

◆ resourceNameForOpticalProps()