What it solves

Electrostatic potential for adiabatic electrons, computed from

n_{e} = n_{e 0} exp (\frac{e ϕ}{k _{B} T _{e}}) \approx n_{e 0} (1 + \frac{e ϕ}{k _{B} T _{e 0}}),

where $n_{e 0}$ and $T_{e 0}$ denote electron equilibrium density and temperature, respectively.

This is solved in weak form: find $ϕ \in H^{1}$ such that

\int_{Ω} ψ \frac{n _{e 0} ( x )}{T _{e 0} ( x )} ϕ d x = \int_{Ω} ψ (n_{e} (x) - n_{e 0} (x)) d x \forall ψ \in H^{1} .

The equation is discretized as

σ_{1} M_{n / T}^{0} ϕ = (Λ^{0}, n_{e} - n_{e 0})_{L^{2}},

where $M_{n / T}^{0}$ is a struphy.feec.mass.WeightedMassOperator and $σ_{1}$ is a normalization parameter.

Variables it advances

This propagator still inherits the abstract variable template from Propagator and does not declare its own variables.

This propagator takes no configuration options.

No model in the catalogue currently composes this propagator — it is available for building new models.