previous up next contents Previous: 3.2 Maximum Entropy Closure Up: 3.2 Maximum Entropy Closure Next: 3.2.2 Critique and Modifications

3.2.1 Distribution Families

Depending on the set of constraints different classes of distribution functions are obtained. For our applications the constraints can be divided in two classes: the first consists of a set of determined moments, the second consists of additional constraints (related to the isotropicity condition) on the distribution functions.

Here we will only study the closure from the reduction conditions 2.40. It is also possibly to study the closures from other constraints, e.g., isotropic symmetric part or the shifting (linear-isotropic) assumption, but the corresponding families cannot be described using elementary functions.

The maximum entropy closure from the reduction conditions is easily achieved, by noting that they determine additional moments, which have to be included in the ansatz.

Hence we get the following family of distributions

$\displaystyle f(r,k_x) = \exp(a_0 + a_1 k_x + a_2 r^2 + b_2 k_x^2 + a_3 k_x r^2 + a_4 r^4 + b_4 k_x^2 r^2 + a_5 k_x r^5)   .$ (3.8)

Then the 7 free parameters have to be determined from the 5 known moments plus the two moments which are additionally known from the reduction conditions.

If a function $ f$ is written as

$\displaystyle f = f(k_x,r)$ (3.9)

with $ r$ the spherical (not cylindrical, though we assume cylindrical symmetry) radius, then the integration takes the form

$\displaystyle \int_{\mathbb{R}^3} f d^3x = \int_0^{\infty} dr 2 r \pi \int_{-r}^{r} dk_x f(k_x,r)   .$ (3.10)

previous up next contents Previous: 3.2 Maximum Entropy Closure Up: 3.2 Maximum Entropy Closure Next: 3.2.2 Critique and Modifications


R. Kosik: Numerical Challenges on the Road to NanoTCAD