plasma_properties package

Submodules

plasma_properties.error module

exception plasma_properties.error.DimensionMismatchError

Bases: plasma_properties.error.Error

Raised when arrays don’t match size

exception plasma_properties.error.Error

Bases: Exception

Base class for other exceptions

exception plasma_properties.error.TypeError

Bases: plasma_properties.error.Error

Raised when type is incorrect

exception plasma_properties.error.ValueError

Bases: plasma_properties.error.Error

Raised when value is incorrect

plasma_properties.parameters module

class plasma_properties.parameters.Params(Am, mass_density, T, Z)

Bases: object

Compute common plasma parameters.

Parameters

• Am (float or arrary_like) – Atomic mass of element (or isotope) in units of grams [g].

• mass_density (float or array_like) – Range of mass densities in units of grams per cubic centimeter [g/cc].

• T (float or array_like) – Temperature range in units of electron-volts [eV]

• Z (int or arrray_like) – Atomic number for each element

• units_out (str) – Unit system for resulting transport coefficient. Default: dimensionless “star” units.

aws()

Wigner-Seitz radius

\[a_{ws}= (4 \pi n_i/3)^{-1/3}\]

degen()

degeneracy parameter

\[\theta = E_f/T\]

where \(E_F\) is the Fermi energy \(E_F = \frac{\hbar^2}{2m_e}(3 \pi^2 n_e)^{2/3}\).

gamma()

Coulomb coupling parameter

\[\Gamma = \dfrac{\langle Z \rangle^2 e^2}{a_i T}\]

kappa()

Inverse electron screening length

\[\kappa = \dfrac{a_i}{\lambda_{TF}},\]

where \(\lambda_{TF}\) is the Thomas-Fermi screening length.

wp()

Ion plasma frequency

\[\omega_p = \left(\dfrac{4 n_i \langle Z \rangle^2 e^2}{m_i}\right)^{1/2}\]

plasma_properties.transport module

class plasma_properties.transport.SM(Am, mass_density, T, Z, units_out='star')

Bases: object

Generate the Stanton and Murillo transport coefficients. test

Parameters

• Am (float or arrary_like) – Atomic mass of element (or isotope) in units of grams [g].

• mass_density (float or array_like) – Range of mass densities in units of grams per cubic centimeter [g/cc].

• T (float or array_like) – Temperature range in units of electron-volts [eV]

• Z (int or arrray_like) – Atomic number for each element

• units_out (str) – Unit system for resulting transport coefficient. Default: dimensionless “star” units.

References

1

Stanton, Liam G., and Michael S. Murillo. “Ionic transport in high-energy-density matter.” Physical Review E 93.4 (2016): 043203.

knm(g, n, m)

Computes the plasma parameters (e.g. ion plasma frequency, ion-sphere radius, coupling parameter, etc.).

Parameters

• g (float or array_like) – Plasma parameter (eq. 54 from [1])

• n (int) – Subscript for collision intergral Knm (eq. C22 from [1])

• m (int) – Subscript for collision integral Knm (eq. C22 from [1])

Returns

knm – Fit to collision integral (eqs. C22-C24 from [1])

Return type

array_like

plasma_params(Z, ni, mi)

Computes the plasma parameters (e.g. ion plasma frequency, ion-sphere radius, coupling parameter, etc.).

Parameters

• Z (int) – Atomic number

• ni (float) – Number density [1/cc]

• mi (float) – Atomic mass [g]

plot(X, Y, xaxis='temeprature')

Very rough preliminary plotting method…. in development.

self_diffusion()

Computes the self-diffusion coefficient using eq. 56 from [1].

Returns

D – Self-diffusion coefficients for system parameters.

Return type

float or array_like

Note

The structure in which D is returned is in the form of a 3D array (1D if single value inputs). The first axis corresponds to density, the second is temperature, and the third is atomic number.

thermal_conductivity()

Computes the ion thermal conductivity coefficient using eq. 82 from [1].

Returns

K – Ion thermal conductivity coefficients for system parameters.

Return type

float or array_like

Note

The structure in which K is returned is in the form of a 3D array (1D if single value inputs). The first axis corresponds to density, the second is temperature, and the third is atomic number.

viscosity()

Computes the viscosity coefficient using eq. 75 from [1].

Returns

eta – Viscosity coefficients for system parameters.

Return type

float or array_like

Note

The structure in which eta is returned is in the form of a 3D array (1D if single value inputs). The first axis corresponds to density, the second is temperature, and the third is atomic number.

class plasma_properties.transport.YVM(Am, mass_density, T, Z, units_out='star')

Bases: object

Compute viscosity from the Murillo viscosity model.

Parameters

• Am (float or arrary_like) – Atomic mass of element (or isotope) in units of grams [g].

• mass_density (float or array_like) – Range of mass densities in units of grams per cubic centimeter [g/cc].

• T (float or array_like) – Temperature range in units of electron-volts [eV]

• Z (int or arrray_like) – Atomic number for each element

• units_out (str) – Unit system for resulting transport coefficient. Default: dimensionless “star” units.

References

2

Murillo, M. S. “Viscosity estimates for strongly coupled Yukawa systems.” Physical Review E 62.3 (2000): 4115.

viscosity()

plasma_properties.zbar module

class plasma_properties.zbar.MeanIonization(Am, mass_density, T, Z)

Bases: object

Compute the mean ionization state (<Z> or Zbar) for a given system.

Parameters

• Am (float or arrary_like) – Atomic mass of element (or isotope) in units of grams [g].

• mass_density (float or array_like) – Range of mass densities in units of grams per cubic centimeter [g/cc].

• T (float or array_like) – Temperature range in units of electron-volts [eV]

• Z (int or arrray_like) – Atomic number for each element

tf_zbar()

Thomas Fermi Zbar model.

References

Finite Temperature Thomas Fermi Charge State using R.M. More, “Pressure Ionization, Resonances, and the Continuity of Bound and Free States”, Adv. in Atomic Mol. Phys., Vol. 21, p. 332 (Table IV).

Module contents