jaxrts.plasma_physics

This submodule contains basic formulas used in plasma physics.

Functions

Debye_Hueckel_screening_length(n, T[, Z])

Calculate the Debye-Hückel screening length.

calculate_mean_free_charge_more(rho, m_A, ...)

Finite Temperature Thomas Fermi Charge State using an analytical fit provided by [More, 1985] p.

chem_pot_interpolationIchimaru(T, n_e)

Interpolation function for the chemical potential between the classical and quantum region, originally from [Ichimaru, 2018], eqn.(3.147).

chem_pot_sommerfeld_fermi_interpolation(T, n_e)

Interpolation function for the chemical potential of a non-interacting (ideal) fermi gas given in the paper of [Cowan, 2019].

compton_energy(probe_energy, scattering_angle)

coulomb_potential_fourier(Z1, Z2, k)

The Fourier transform of the Coloumb potential.

coupling_param(Z1, Z2, n_e, T_e)

Returns the degree of interparticle coupling with corresponding charge numbers Z1 and Z2 at temperature T_e and density n_e.

degeneracy_param(n_e, T_e)

Calculate the plasma degeneracy parameter.

epsilon_from_susceptibility(xi, k)

Calculate the dielectric function from a full susceptibility xi

fermi_dirac(k, chem_pot, T)

Return the Fermi-Dirac distribution.

fermi_energy(n_e)

Calculate the Fermi energy of an ideal fermi gas from a given electron density.

fermi_temperature(n_e)

Calculates the Fermi temperature of an ideal fermi gas from a given electron density.

fermi_wavenumber(n_e)

interparticle_spacing(Z1, Z2, n_e)

kin_energy(k)

Kinetic energy of a free electron with wavevector k.

noninteracting_susceptibility_from_eps_RPA(...)

Calculates the non-interacting susceptilibily from a given dielectric function epsilon in RPA.

plasma_frequency(electron_density)

Calculate the plasma frequency \(\omega_\text{pe}\)

susceptibility_from_epsilon(epsilon, k)

Calculate the full susceptilibily from a given dielectric function epsilon by inverting

temperature_interpolation(n_e, T_e[, power])

Interpolate between electron temperature and Fermi temperature.

therm_de_broglie_wl(T)

thomson_momentum_transfer(energy, angle)

Momentum transfer \(k = \mid\vec{k}\mid\), assuming that the absolute value of the momentum for incoming and scattered light is only slightly changed.

wiegner_seitz_radius(n_e)

Return the Wiegner Seitz Radius \(r_s\).