*d*P = -P *d**ξ* Σ n_{s}σ_{s} = -P
*d**ξ*/*l*
(2)

where the sum over s accounts for all different types of
centers, and *l* = 1/(Σn_{s}σ_{s})
is the mean free path.

The probability that the last collision of the particle, before reaching a cross section plane at some
*ξ* , was along the element *d**ξ'*, at
a position *ξ'*. (
*ξ'* < *ξ*, figure (1)), is:

P(*ξ'*) *d**ξ'*/*l* =
*e*^{(ξ' - ξ)/l}*d**ξ'*/*l*
(3)

*
**l *is determined by the properties of the medium
and the name "mean free path" is justified, since

_{-∞}∫^{ξ}(*ξ' - ξ*) P(*ξ'*)
*d**ξ'*/*l* = l, that is, *l*
is the probable distance where the particle comes from.

The velocity distribution __v__ of a particle, at the cross section plane
*ξ*, is obtained by adding the local velocity distributions at
all the distance elements along the line up to *ξ*. Each
distribution is multiplied by the probability that the particle comes from that distance,
namely, that it has collided for the last time with the medium at that place :

f(*ξ*,v) = _{-∞}∫^{ξ}f_{0}(*ξ'*,v) P(*ξ'*)
*d**ξ'*/*l*
(4)

where v = |__v|__ and f_{0}( *ξ'*,
v) is the local equilibrium distribution at *ξ'*. If the equilibrium distribution
does not change significantly along the mean free path *l*, it will be possible to use the
first terms of Taylor series:

f_{0}(*ξ'*,v) = f_{0}(*ξ*,v) +
(∂f_{0}/∂*ξ*) (*ξ' - ξ*) +
(∂f_{0}/∂v) (v(*ξ'*) - v(*ξ*))
(5)

The change in velocity is a result of force acting on the particle. For
example, a gravitational field, or an electric field for a charged particle.
Assuming that the force acts in the x-direction, the velocity change will
be:

v(*ξ'*)^{2} - v(*ξ*)^{2}
≈ 2v (v' - v) = 2a (*ξ'* - *ξ*) *Cos*(θ)
(6)

where a is the acceleration, equal to eE/m for an electric field E.

Substituting (6) and (5) into (4), integrating over
*ξ'*, and finally substituting x
and θ for *ξ*
according to (1):

f(x, v, θ) = f_{0}(x, v) - *l* *Cos*(θ) [( ∂f_{0} / ∂x) + (a/v)
(∂f_{0}/∂v)]
(7)

The non-equilibrium distribution (7) is the linear Boltzmann transport
equation in a somewhat non-traditional form.

f_{0}(x, v) is the local equilibrium distribution given by Maxwell-Boltzmann or
Fermi-Dirac distributions. The medium properties enter through the mean free
path *l*, and external forces
through the acceleration a. The standard form of the linear
Boltzmann transport equation includes a constant relaxation time τ,

∂f/∂t + v grad_{x}(f) + a grad_{v}(f) = -(f - f_{0})/τ
, where τ =
*l*/v .

**3. The currents **

A particle moving at a velocity __v__ will cross a plane section at x during a
time Δt if its distance from
the plane is less than

|__v__| *Cos*(θ)
Δt.

The particle current through a unit cross section at x is obtained
by summing all the velocities of all directions in space (spatial angle of 4π),
each velocity weighed by the distribution probability:

J_{n} =
∫f(x, v, θ) |__v__| *Cos*(θ)*d*^{3}__v__
(8)

where *d*^{3}__v__ = 2π v^{2 }*d*v
*Sin*(θ) *d*θ. Substitution of the distribution function
(7) into (8) yields the current:

J_{n} = -(4π*l*/3) [(eE/m) ∫(∂f_{0}/∂v)
v^{2} *d*v + ∫(∂f_{0}/∂x) v^{3} *d*v]
(9)

where the f_{0} contribution in (7) is zero, and
the trigonometric integral is

∫*Cos*^{2}(θ)* Sin*(θ) *d*θ = 2/3.

Since the energy associated with each particle is (1/2)mv^{2} the
energy current will be:

J_{u} = (m/2)
∫f(x, v, θ) |__v__|^{3} *Cos*(θ) *d*^{3}__v__
(10)

and by similar calculation:

J_{u} = -(4π*l*/3) [(eE/m) ∫(∂f_{0}/∂v)
v^{4} *d*v + ∫(∂f_{0}/∂x) v^{5} *d*v]
(11)

The currents will now be calculated by applying Maxwell-Boltzmann
distribution:

f_{0}(x, v) = n(x) (m/(2πkT))^{3/2}*e*^{-m v^2/(2kT)}
(12)

where n(x) is the local particle density. The
pre-exponential factor in (12) is determined by the normalization condition:

_{0}∫^{∞}4π v^{2} f_{0} *d*v = n.

Since, (by (12)):

∂f_{0}/∂v = -(mv/kT) f_{0}
(13)

then:

J_{n} = (4π*l*/3) (eE/kT - *d* /*d*x) ∫f_{0} v^{3} *d*v
(14)

Transforming to a dimensionless variable x = (m/2kT)^{1/2} v,
the integral will be:

∫f_{0} v^{3} *d*v = n (m/(2πkT))^{3/2} (2kT/m)^{2}
∫x^{3} *e*^{-x^2} *d*x =
(nπ) (kT/2πm)^{1/2},

since
∫x^{3} *e*^{-x^2} *d*x = Γ(2)/2 = 1/2 (Appendix A).

The particle current is then:

J_{n} = [(4*l*/3)/(2πm)^{1/2}]
[eEn/(kT)^{1/2} - *d*(n(kT)^{1/2})/*d*x]
(15)

A similar calculation yields the energy current:

J_{u} = (4π*l*/3) (m/2) (eE/kT - *d* /*d*x) ∫f_{0} v^{5} *d*v
(16)

The integral is:

∫f_{0} v^{5} *d*v =
n (m/(2πkT))^{3/2} (2kT/m)^{3}
∫x^{5} *e*^{-x^2} *d*x =
(nπ) (2kT/m) (kT/2πm)^{1/2},

since
∫x^{5} *e*^{-x^2} *d*x = Γ(3)/2 = 1 (Appendix A).

The energy current is then:

J_{u} = 2[(4*l*/3)/(2πm)^{1/2}] [(kT)^{1/2} eEn -
*d*(n(kT)^{3/2})/*d*x]
(17)

or, by combining (15) and (17):

J_{u} = 2kT J_{n} - 2[(4*l*/3)/(2πm)^{1/2}]
n (kT)^{1/2} *d*(kT)/*d*x
(18)

This equation expresses the energy current as a sum of two terms. The first
is convection, the energy associated with the particle flow; and the second is
heat conduction, which is proportional to the temperature gradient and
independent of the particle current.

In electrical conductors (but not in semiconductors) the charge density does
not depend on the temperature. The electrical current J_{q} will then be (by (15):

J_{q} = eJ_{n} = e [(4*l*/3)/(2πm)^{1/2}]
(1/kT)^{1/2 } [eEn - (n/2) *d*(kT)/*d*x]
(19)

**4. The transport coefficients**

The transport coefficients are
determined directly from the currents in equations (15)-(19). The diffusivity
D, defined by

J_{n}(*d*T/*d*x = 0, E = 0) = -D *d*n/*d*x, is:

D = [(4*l*/3)/(2πm)^{1/2}] (kT)^{1/2}
(20)

The electrical conductivity σ,

defined by J_{q} =
eJ_{n}(*d*n/*d*x = 0, *d*T/*d*x = 0) = σE, is:

σ = [(4*l*/3)/(2πm)^{1/2}]e^{2} n/(kT)^{1/2}
(21)

The thermal conductivity, defined by J_{u}(J_{n} = 0) = -K
*d*T/*d*x, is:

K = 2[(4*l*/3)/(2πm)^{1/2}] (kT)^{1/2} k n
(22)

By comparing these coefficients, the Wiedemann-Franz law:

K/σ = 2kT k/e^{2}
(23)

and Einstein relation:

D/μ = D/(σ/ne) = kT/e
(24)

follow directly. μ is the mobility

The transport coefficients, and the relationships between them, are
determined from the currents (equations (15)-(19)) calculated from Boltzmann
transport equation (7).

**5. Particles flow within a medium with a temperature gradient**

The
particle flow may be that of a gas within a narrow pipe, a porous medium, or
that of atoms diffusing within a solid. The current of uncharged particles is,
(by (15)):

J_{n} = -[(4*l*/3)/(2πm)^{1/2}]
[*d*(n(kT)^{1/2})/*d*x]
(25)

If there is no time dependent accumulation or depletion of particles the
current will not depend on x and then the equation can be
integrated to yield:

L J_{n} = [(4*l*/3)/(2πm)^{1/2}]
[n_{H}(kT_{H})^{1/2} - n_{C}(kT_{C})^{1/2}]
(26)

where L is the length of the medium, and the subscripts
C and H denote the cold and hot ends. If
the ideal gas law, p = n kT, is obeyed, then n (kT)^{1/2} can be replaced in equations
(25)-(26) by p/(kT)^{1/2}. Another important
result is that the current is determined by the thermodynamic parameters at the
ends of the medium and not by their variations along it. In the special case

J = 0 the following expression becomes a system
invariant:

n (kT)^{1/2} = p/(kT)^{1/2} = Const
(27)

In particular n_{C} (kT_{C})^{1/2} =
n_{H} (kT_{H })^{1/2}, and there is no dependence on the
medium properties. This is Knudsen law (Knudsen (1952)).

**6. Irreversible thermodynamics**

Appendix-B, equations (64)-(83), summarizes equilibrium thermodynamics and some results that are applied here.
Consider the non-equilibrium system of section-2, figure-1, where the temperature and particle
density vary along the x-direction.
The system can be divided into sub-volumes ΔV_{i} =
A Δx_{i}
where A is the area cross-section.
If Δx_{i} is equal to the mean free path then
local equilibrium thermodynamic functions may be used within each
volume ΔV_{i}.
Assume that S(U, N) is the local entropy defined in a sub-volume.
The change of entropy ΔS within two neighboring sub-volumes,

i and i + 1, will be:

ΔS_{i} =
(∂S/∂U)_{i} ΔU_{i} + (∂S/∂N)_{i} ΔN_{i}
(28)

ΔS_{i+ 1} =
(∂S/∂U)_{i + 1} ΔU_{i+ 1} + (∂S/∂N)_{i+ 1} ΔN_{i+ 1}
(29)

Since particles and energy flow from one sub volume to the other:

ΔU_{i+ 1} = -ΔU_{i} = ΔU
(30)

ΔN_{i+ 1} = -ΔN_{i} = ΔN
(31)

The entropy of the two sub volumes will be:

ΔS = ΔS_{i} + ΔS_{i+ 1} =
Δ(∂S/∂U) ΔU + Δ(∂S/∂N) ΔN
(32)

Assuming that the entropy is locally not far from equilibrium, the equilibrium values conditions
∂S/∂U = 1/T and ∂S/∂N = -μ/T
(equations (70)-(71)) may be applied.

ΔS = ΔU Δ(1/T) - ΔN Δ(μ/T)
(33)

Dividing the local thermodynamic parameters and functions by the local volume ΔV yield their densities,

s = S/ΔV, u = U/ΔV, n = N/ΔV.
The parameter differences are replaced by their gradients,
in this example by the derivative along the x direction,

Δ(1/T) = *l* *d*(1/T)/*d*x, and
Δ(μ/T) = *l* *d*(μ/T)/*d*x. Finally, the time derivative of (33) yields:

*d*(Δs)/*d*t = *l* *d*(1/T)/*d*x *d*(Δu)/*d*t –
*l* *d*(μ/T)/dx *d*(Δn)/*d*t
(34)

If the local entropy is not far from equilibrium, the rate of entropy production may be expressed in terms of flows
and forces:

*d*(Δs)/*d*t = X_{u} J_{u} + X_{n} J_{n}
(35)

where, by comparison to (34), J_{u} = *d*(Δu)/*d*t,
and J_{n} = *d*(Δn)/*d*t.
The forces are:

X_{u} = *d*(1/T)/*d*x
(36)

X_{n} = -*d*(μ/T)/*d*x
(37)

The flows are related to the forces by linear phenomenological equations:

J_{n} = -L_{11} *d*(μ/T)/*d*x + L_{12} *d*(1/T)/*d*x
(38 )

J_{u} = -L_{21} *d*(μ/T)/*d*x + L_{22} *d*(1/T)/*d*x
(39 )

Onsager theorem postulates that the linear coefficients, L_{ij}, are symmetrical:
L_{12} = L_{21}.

Inserting the chemical potential of mono-atomic ideal gas,

μ = -T ∂S/∂N =
kT *ln*[(T^{3/2}/n)/(T_{0}^{3/2}/n_{0})],

calculated in appendix-B equation (83), the derivative will be:

*d*(μ/T)/*d*x =
(3k/2) *d*(*ln*(T))/*d*x - k *d*(*ln*(n))/*d*x

= (3k/2T) *d*T/*d*x - (k/n) *d*n/*d*x
(40)

The flows will be:

J_{n} = L_{11} (k/n) *d*n/*d*x
- [L_{11} (3k/2T) + L_{12} (1/T)^{2}] (*d*T/*d*x)
(41 )

J_{u} = L_{21} (k/n) *d*n/*d*x
- [L_{21} (3k/2T) + L_{22} (1/T)^{2}] (*d*T/*d*x)
(42 )

The particle flow J_{n} and the energy flow
J_{u} are coupled by cross terms with
symmetrical linear coefficients. This result is how far irreversible thermodynamics takes us. It gives little information,
if any, since the physics of the system is included in the explicit expressions of the L_{ij}s
that are not calculated by the theory.

However, the particle and energy currents (15)-(18), without external forces, given by Boltzmann's transport calculation,
with

A = [(4*l*/3)/(2πm)^{1/2}], are:

J_{n} = -A
[(kT)^{1/2} (*d*n/*d*x) + ((n/2)/(kT)^{1/2}) *d*(kT)/*d*x]
(43)

J_{u} = -2A
[(kT)^{3/2} (*d*n/*d*x) + (3n/2) (kT)^{1/2} *d*(kT)/*d*x]
(44)

These equations (43)-(44) do not contain unknown coefficients. Explicit values of the linear coefficients of irreversible thermodynamics
may be calculated by comparing them to the similar flow equations (41)-(42):

L_{11} = -A (n/k) (kT)^{1/2}
(45)

L_{12} = -2A (n/k) (kT)^{3/2}
(46 )

L_{21} = -2A (n/k) (kT)^{3/2}
(47 )

L_{22} = -6A (n/k) (kT)^{5/2}
(48 )

The cross terms, L_{12} and
L_{21}, are equal, so that they verify Onsager's theorem.

**7. Charge Transport through a junction**

Electrons in solids occupy energy levels in pairs of opposite spins.
Lower levels are fully occupied and upper levels are empty, and the level population is
determined by Fermi-Dirac statistics. An electron may be scattered by a center within a medium only to
an empty level. For this reason low energy electrons are not scattered by the medium, since all the
nearby levels are occupied. Maxwell-Boltzmann velocity distribution is not applicable in this case.

Semiconductors are a special case. In these materials a fully occupied valence band is separated
from an empty conduction band by a forbidden gap without energy levels.
Doping the material with donors will introduce electrons to the conduction band and make it n-type.
Acceptors will introduce holes to the valence band and make it p-type. A donor atom has an
extra electron at a level just below the bottom of the conduction band. Thermal excitation
will transform it to the conduction band where it is free to move. An acceptor atom misses an electron
just above the valence band. Electron caught by the acceptor from the valence band will leave a hole
behind it that is free to move within the valence band. Practically, all donors are ionized so that the
electron-density is equal to the donor density. Similarly, the hole-density is equal to the acceptor density.

Maxwell-Boltzmann velocity distribution may apply to semiconductors where the charge density,
either electrons or holes, is lower than the energy level density.

Consider two pieces of a semiconductor with higher n_{L} and lower
n_{R} electron densities.
The left is n-type where n_{L} is equal to the donor density.
The right may be n-type or p-type. For p-type
material n_{R} is determined by the condition

p_{R} n_{R} = p n = Const,
that is general to semiconductors. p_{R} is equal to the acceptor density.

If the two pieces are brought into contact to form a junction, electrons will start to flow from
the left n_{L} side to the right
n_{R} side.
This flow generates space charge region, electric field,
and potential difference V_{c} between the
two pieces that stop further electron flow. V_{c} is
the contact potential. It is not directly measureable by a voltmeter, for example, because the voltmeter
probes make their own contact potentials with the two ends of the junction. The sum of all junction potentials is then zero.

If the two pieces are separated after being in contact, their capacity will drop, the voltage on
them will grow and it will be possible to measure it. Methods that measure the contact potential are based on this phenomenon.

It is assumed that the junction is perpendicular to the x-direction, and that the depth of the space charge layer in it is
small compared to the electron's mean free path. The thermal velocity of electrons at the left side must
have a minimal component v_{x} in
order to overcome the potential barrier and cross to the right side:
(1/2) m v_{x}^{2} > eV_{c}.
Electrons at the right side do not have any barrier for crossing back to the left.

An electron that crosses the junction from left to right, will be slowed down by the electric field
and lose the energy eV_{c}. This energy will be thermally
regained from its vicinity. An electron moving the opposite way will be accelerated by the field
and thermally deliver the field-gained energy to its vicinity.

For this case Boltzmann transport equation f
is approximated by the first non zero term in equation (7), which is
the equilibrium distribution f_{0}.
The current, equation (8),through the contact will be:

J_{n} = ∫_{(vx > v0)} f_{0}(v) v_{x} *d*^{3}__v__ +
∫_{(vx < 0)} f_{0}(v) v_{x} *d*^{3}__v__ =
J_{n}^{1} –J_{n}^{2}
(49)

where: v_{0} =(2 eV_{c} / m)^{1/2}.

It is convenient to use here Cartesian coordinates
*d*^{3}__v__ = *d*v_{x}
*d*v_{y} *d*v_{z}.
By inserting Maxwell-Boltzmann distribution, equation-(12) , for f_{0} :

J_{n}^{1} = n_{L} (m /2πkT)^{3/2}
_{0}∫^{∞} *e*^{-αz} *d*v_{z}
_{0}∫^{∞} *e*^{-αy} *d*v_{y}
_{v0}∫^{∞} *e*^{-αx} v_{x} *d*v_{x}
(50)

J_{n}^{2} = n_{R} (m /2πkT)^{3/2}
_{-∞}∫^{0} *e*^{-αz} *d*v_{z}
_{-∞}∫^{0} *e*^{-αy} *d*v_{y}
_{-∞}∫^{0} *e*^{-αx} v_{x} *d*v_{x}
(51)

where: α_{i} = mv_{i}^{2} / (2kT), i = x, y, z.

Using the formulas, appendix-A, (a > 0):

_{0}∫^{∞} *e*^{-(a x^2)} *d*x = (π/a)^{1/2}/2
(52)

and

_{0}∫^{∞} *e*^{-(a x^2)} x *d*x = 1/(2a)
(53)

the integrals will be:

_{0}∫^{∞} *e*^{-αi} *d*v_{i} = (πkT/m)^{1/2}
(54)

_{v0}∫^{∞} *e*^{-αx} v_{x}
*d*v_{x} = (kT/m) *e*^{-eVc/ kT}
(55)

_{-∞}∫^{0} *e*^{-αx} v_{x}
*d*v_{x} = (kT/m)
(56)

And the overall current:

J_{n} = J_{n}^{1} – J_{n}^{2} =
(kT/8πm)^{1/2} n_{L} (*e*^{-eVc/kT} –
n_{R}/n_{L})
(57)

Since the system is at equilibrium, J_{n} = 0, the potential difference at the contact will be:

*e*^{-eVc/kT} = n_{R}/n_{L}
(58)

or:

V_{c} = (kT/e) *ln*(n_{L}/n_{R})
(59)

If an external voltage source V is connected to the system, and the voltage falls on the
junction, the current will be:

J_{n} = (kT/8πm)^{1/2} n_{L} (*e*^{e(V – Vc)/kT} –
n_{R}/n_{L})
(60)

Or, by (58):

J_{n} = (kT/8πm)^{1/2} n_{R} (*e*^{eV/kT} – 1)
(61)

This formula is the I–V diode equation.

**8. Thermoelectric effects**

**Thomson effect:** Generation or absorption of heat in a current carrying conductor.

Generation or absorption of heat along a conductor is given by the derivative of the energy current,

J_{u} = 2kT J_{n} - 2((4*l*/3)/(2πm)^{1/2})
n (kT)^{1/2} *d*(kT)/*d*x
,
equation-(18), with respect to
x.
Assuming that the particle flow does not depend on x,
*d*J_{n}/*d*x = 0, and assuming also
that the temperature change is linear with the distance, *d*(kT)/*d*x = Const,
the heat generation or absorption will be:

*d*J_{u}/*d*x =
[2 J_{n} -((4*l*/3)/(2πm)^{1/2})
n / (kT)^{1/2}] *d*(kT)/*d*x
(62)

**Peltier effect:** Heating or cooling a junction by electric current.

An electron that crosses a junction will be either slowed down, or accelerated by a contact potentiall difference
V_{c}, depending on its polarity.
Thus, it will absorb or deliver an amount of heat e V_{c} to the junction.
V_{c} is given by
V_{c} = (kT/e) *ln*(n_{L}/n_{R}),
equation (59).

**Seebeck effect:**
Direct power generation from temperature difference.

Direct conversion of thermal energy to electrical energy by a junction is discussed
by applying the original Carnot method to a junction device operating in a reversible closed cycle.
The "Carnot machine" shown in figure-2 has some resemblance to devices
that are used to measure contact potentials (Klein et al (2003)).

Figure-2: A "Carnot machine" comprises of two junctions connected in parallel. One with two semiconductors in permanent
contact, and one with the two parts separable. The distance between them is externally controlled.

Two junctions are connected in parallel. In the left junction two semiconductors are in permanent contact while in the
right junction they are separable and the distance between them is externally controlled.
The connecting wires are made of the same semiconductor material so that there is only one contact junction in the device.
The contact potential difference within the junction is given by
V_{c} = (kT/e) *ln*(n_{L}/n_{R}),
equation-(59).
However, there is no external voltage difference between the junction ends.
Since the junction parts are electrically charged, there will be attraction force between them. Distance change
between the separated parts requires mechanical work.

The machine operates in a four step cycle between a hotter heat-bath T_{H}
and a colder heat-bath T_{C}.

In step-1 the machine is attached to the hot bath T_{H}
and the distance between the two parts of the right junction
is reduced by moving the upper part down. The junction potential V_{H}
does not change during this isothermal movement and the increasing capacity of the right junction
drives charge flow Δq in the two junction circuit.
The flow through the left junction absorbs heat

Δq V_{H}
from the hot bath (Peltier effect), and the right junction performs
an equivalent work on the external force that holds it in place.

In step-2 the upper part of the right junction continues to move down while
the machine is separated from the hot bath T_{H}.
During this adiabatic step the machine cools down and the movement continues
until it reaches the temperature T_{C} of the cold bath,
while the junction potential decreases from V_{H}
to V_{C}.

In step-3 the machine is attached to the cold bath T_{C} and
the distance between the parts of the right junction
is increased by moving the upper part back up. The junction potential V_{C}
does not change during this isothermal movement and the decreasing capacity of the right junction
drives charge flow Δq in the opposite direction to step-1.
The flow through the left junction delivers heat Δq V_{C}
to the cold bath (Peltier effect the other way) and the external force performs
an equivalent work on the right junction.

In step-4 the upper part of the right junction continues to move up while
the machine is separated from the cold bath T_{C}.
During this step the machine heats up and the adiabatic movement continues
until it reaches the temperature T_{H} of the hot bath.
The junction potential increases from V_{C}
to V_{H} and the machine completes the cycle and returns to its original state.

The adiabatic work done by the machine on the external force in step-2 is equal to the adiabatic work
done by the external force on the machine in step-4. The electro-motive force V of
the machine is the work done by moving a unit charge in a complete closed cycle. This work is equal to the sum of
of step-1 and step-3 for a unit charge e.

Using
V_{c} = (kT/e) *ln*(n_{L}/n_{R}),
equation (59), this work is equal to:

eV =
k(T_{H} - T_{C}) *ln*(n_{L}/n_{R})
(63)

Combining equations (59) and (63) yields Carnot's efficiency
1 - T_{C} / T_{H}.

**9. Summary and conclusions**

Boltzmann transport equation relates
the properties of a non-equilibrium system, expressed by a non-equilibrium
distribution, in terms of local equilibrium distributions. The equation is
derived by assuming random movement of particles within a medium.

Boltzmann transport equation is applied to calculation of the general
currents in a medium with particle and temperature gradients. The currents
determine the transport coefficients of the medium.

The calculated diffusivity, electrical and thermal conductivities verify
Wiedemann-Franz and Einstein Laws.

The calculated current of uncharged particles verifies Knudsen law.

See:

Expansion of an Ideal Gas

Osmosis and Thermodynamics

van't Hoff's Evidence

**Appendix A: integrals**

_{0}∫^{ ∞}x^{n
}*exp*(-a^{2}x^{2})
*d*x = Γ((n+1)/2) / 2a^{n+1}

Γ is the gamma function:

Γ(n) = (n – 1)!

Γ(*x* + 1) = *x* Γ(*x*)

Γ(*x*)
Γ(*x *+ 1/2) = Γ(2*x*)
π^{1/2} / 2^{2x - 1}

Γ(1/2) = π^{1/2}

Γ(3/2) = (1/2) π^{1/2}

Γ(5/2) = (3/4) π^{1/2 }

**Appendix B: Summary of equilibrium thermodynamics**

Each state of a quantum system is assumed to have the same probability. Ω(U)
is the system's density of states as a function
of the energy U, and the entropy is:

S(U) = k ln(Ω(U))
(64)

where k is Boltzmann's constant. The entropy is maximal when the system is
at thermodynamic equilibrium.

The temperature is:

1/T = *d*S(U)/*d*U
(65)

x is a parameter of the system, volume, number of particles, etc.,
and its conjugate force X is:

δU = -X δx
(66)

The first law of thermodynamics, for a system in thermal contact with a heat bath, of a temperature
T, is:

δU = δQ – Σ_{i}
X_{i }δx_{i} = δQ – pδV + μδN + …
(67)

where the volume V, the number of particles N, ... are
parameters of the system, the pressure p, the chemical potential
μ, ... are their conjugate forces, and δQ is
the heat flow from the heat bath. The first law is conservation of energy.

The overall entropy S*, for a system at equilibrium with a heat bath,
is the sum of the system entropy S(U, x_{i})
and the bath entropy S_{T}(U_{T}).

S* = S(U, x_{i}) + S_{T}(U_{T})
(68)

The maximum condition for the entropy is:

0 = *d*S* = (∂S/∂U) *d*U + Σ_{i }(∂S/∂x_{i}) *d*x_{i} +
(∂S_{T}/∂U_{T}) *d*U_{T}

= (∂S/∂U) *d*U + Σ_{i} (∂S/∂x_{i}) *d*x_{i} +
(∂S_{T}/∂U_{T}) (-*d*U – Σ_{i}
X_{i} *d*x_{i})

= [(∂S/∂U) – (∂S_{T}/∂U_{T})] *d*U + Σ_{i}
[(∂S/∂x_{i}) – (∂S_{T}/∂U_{T}) X_{i}] *d*x_{i}
(69)

Therefore, the conditions for thermodynamic equilibrium are:

(∂S/∂U) = (∂S_{T}/∂U_{T}) = 1/T
(70)

(∂S/∂x_{i}) = X_{i}/T
(71)

**Reversible process**

When a parameter change δx is slow enough the system will smoothly
adjust its equilibrium state to it. If the change δx is reversed
the system will reverse its response and return to its original state.
A process that at any moment can reverse its direction, in response to reverse of a parameter
change, is a reversible process.

The energy conservation for an equilibrium system in contact with a heat bath is, by (67):

δQ = δU + Σ_{i} X_{i} δx_{i}
(72)

δQ is the heat supplied from the bath.
Using the equilibrium conditions (70), (71):

δS = (∂S/∂U) δU + Σ_{i} (∂S/∂x_{i}) δx_{i}
= (1/T) δU + Σ_{i} (X_{i}/T) δx_{i} = δQ/T
(73)

or

δS = (1/T) δU + (p/T) δV - (μ/T)δN + ... = δQ/T
(74)

The equation δS = δQ/T is the classical definition of entropy.

The overall entropy change in a reversible process is zero,
ΔS^{overall} = 0. If a system is thermally isolated from its vicinity,
then the process is adiabatic and the entropy change of the system is zero,
ΔS^{system} = 0. If a system is in thermal contact with a heat bath,
then the process is isothermal and the sum entropy change of the system
and the bath is zero, ΔS^{system} + ΔS^{bath} = 0.
Entropy change between two states of a system may be calculated
by following a reversible path between them.

**Entropy change in ideal gas**

A system is defined by three parameters out of four, the pressure p,
number of particles N,
volume V, and the temperature T.
These parameters are connected by the ideal gas
equation: p V = N kT.
A standard system is similarly defined by three parameters out
of p_{0}, N_{0},
V_{0},
and T_{0}, where p_{0} V_{0}
= N_{0} kT_{0}.

n = N/V is the particle density.

The entropy change is calculated for moving ΔN particles, by a reversible process,

from a standard system (p_{0}, N_{0},
V_{0}, T_{0})
to a specific system (p, N,
V, T)

Consider the standard system (p_{0}, N_{0},
V_{0}, T_{0}),
in contact with a heat bath. One wall of the system has a small hole in it, and there is a moving partition
that encloses an additional volume ΔV_{0}.

Figure 3: Entropy of an ideal gas.

The partition is first attached to the wall with the hole, so that the initial
volume ΔV_{0} is zero.
The partition then separates from the wall, moves and increases
the volume ΔV_{0} to some small value,
and removes ΔN particles from the system.
If ΔV_{0} is
much smaller than the volume V_{0},
ΔV_{0} << V_{0}, then the gas pressure
p_{0} will almost not change during the movement.
The isothermal work of pushing the partition by the gas is
W_{1} = p_{0} ΔV_{0}.
The energy for this process, ΔQ_{1}, comes from the heat bath:

ΔQ_{1} = W_{1} =
p_{0} ΔV_{0} = ΔN kT_{0}
(75)

The entropy change for this process is:

ΔS_{1} = ΔQ_{1}/T_{0} = ΔN k
(76)

The hole is now closed. ΔV_{0} is separated
from the volume V_{0}
and is isolated from its vicinity. Now the partition reverses direction,
compresses the gas and heats it adiabatically to the temperature T and
pressure p’. During this process there is no entropy change.

The pressure increase in
adiabatic compression of mono-atomic ideal gas, is (Appendix C):

p’ = p_{0} (T/T_{0})^{5/2}
(77)

The volume ΔV’ is now connected to a heat
bath of a temperature T
and the partition continues to compress it to a
volume ΔV until the pressure
reaches p.

The isothermal work done during this process is:

W_{2} = _{p’}∫^{p}*p* *d*V
= ΔN kT _{p’}∫^{p}*d**p*/*p* = ΔN kT *ln*(p/p’)
(78)

And the entropy change is (by (77)):

ΔS_{2} = ΔQ_{2}/T =
ΔN k *ln*[(p/p_{0}) (T_{0}/T)^{5/2}]
(79)

The volume ΔV is connected now through a hole to the system
(p, N, V, T)
and the partition pushes the ΔN particles into it by reducing
ΔV to zero .

The isothermal work for this last process is:

W_{3} = p ΔV = ΔN kT
(80)

and the entropy change is:

ΔS_{3} = ΔQ_{3}/T = -ΔN k
(81)

The overall entropy change, during these three processes, is:

ΔS = ΔS_{1} + ΔS_{2}
+ ΔS_{3} = ΔS_{2},
since ΔS_{3} = - ΔS_{1}.

Therefore:

ΔS = ΔN k *ln*[(p/p_{0})
(T_{0}/T)^{5/2}]

= ΔN k *ln*((n/n_{0}) (T_{0}/T)^{3/2})

= ΔN k *ln*[(n/T^{3/2})/(n_{0}/T_{0}^{3/2})]
(82)

The chemical potential is (by (71), (74), (82)):

μ = -T ∂S/∂N =
kT *ln*[(T^{3/2}/n)/(T_{0}^{3/2}/n_{0})]
(83)

**Appedix-C: Adiabatic expansion of a mono-atomic ideal gas**

**The gas pressure**

A particle in a box of volume V moves with a velocity
__v__ at an angle φ to the
perpendicular of a wall as shown in figure-4.

Figure-4: Elastic collision of a particle with a wall.

It will collide with the wall during a time Δt, if its distance
from the wall is less than |__v__| Cosφ
Δt. The number
ΔN_{col} of particle collisions with a unit area of
the wall during Δt is:

ΔN_{col} = [n ∫f_{0}(v) |__v__| Cosφ
d^{3}__v__] Δt
(84)

n is the particle density, f_{0}(v)
is the probability that the particle velocity is
v, normalized by

∫4π
v^{2 }f_{0}(v) dv = 1, and
d^{3}__v__ = 2π
v^{2 }dv Sinφ dφ
. Substituting in equation (84):

ΔN_{col} = n (∫2π
v^{2 }f_{0}(v) v dv) (∫Cosφ
Sinφ dφ)
Δt
(85)

The three-dimensional average velocity of particles is
<v> =∫4π
v^{2 }f_{0}(v) v dv. <v> is not zero since the integration is
only over positive values. The trigonometric integral over half the sphere is

∫Cosφ
* *Sinφ* *dφ
= 1, ( 0 < φ < π/2),
therefore:

ΔN_{col} = n (<v>/2)
Δt
(86)

Upon a collision with the wall the particle will transfer a momentum 2mv Cosφ
to the wall. The overall momentum ΔM, delivered
to a unit area of the wall during Δt,
is the sum of momenta delivered by all the colliding particles that come from all
directions of the half sphere. Similarly to equations (84) - (85):

ΔM = n 2m (∫2π
v^{2 }f_{0}(v) v^{2} dv) (∫Cos^{2}φ
Sinφ dφ) Δt
(87)

The three-dimensional average is
<v^{2}> =∫4π
v^{2 }f_{0}(v) v^{2} dv,

∫Cos^{2}φ
Sinφ dφ = 1/3, and the pressure
p =ΔM /
Δt is therefore:

p = n m <v^{2}>/3
(88)

The thermal average is calculated by applying Maxwell-Boltzmann distribution:

f_{0}(v) = (m / 2πkT)^{3 / 2}
exp(-mv^{2} / 2kT) (89)

to the calculation of
<v> and <v^{2}> (appendix A):

<v> = 4π ∫f_{0}(v) v^{3}dv = 2 (2kT /
πm)^{1/2}
(90)

<v^{2}> = 4π
∫f_{0}(v) v^{4} dv = 3kT / m (91)

Substituting <v^{2}> in equation (88)
yields the pressure:

p = n kT (92)

Equation (92) is the state equation of an ideal gas.

**Adiabatic expansion of a mono-atomic ideal gas**

In adiabatic system the particles do not exchange energy with the walls by inelastic collisions and the only energy transfer
mechanism is via elastic collisions with a moving partition.

A single particle in a box of volume V and of unit cross section,
with one wall a moving partition, collides

ΔN_{col} / N
= (v/2) Δt / V times with it
during a time Δt (equation (86)).
Since the relative velocity between the particle and the moving partition is conserved, during
each collision the particle will lose (or gain) a velocity 2v_{M},
twice the partition velocity. A partition
that moves in the x-direction will change the v_{x} component
of the particle velocity. However, the velocity change will eventually spread to the
v_{y} and v_{z} components.
Therefore, the change of a particle velocity during
Δt will be:

dv = -(1/3) 2v_{M} (v/2)
Δt / V = -(v/3) dV/V
(93)

since v_{M}Δt = dV.
The equation of adiabatic expansion of a mono-atomic gas, without internal degrees of freedom of gas particles, is calculated by integration of
equation (93):

v/v_{0} = (V_{0}/V)^{1/3}
(94)

and by equation (88):

p/p_{0} = (V_{0}/V)^{5/3}

or, by the state equation (92), p = n kT:

T/T_{0} = (V_{0}/V)^{2/3}
or
p/p_{0} = (T/T_{0})^{5/2}
(95)

Since:

d(T/T_{0}) =
(-2/3) (V/V_{0})^{-5/3} d(V/V_{0}))
(96)

the work of adiabatic compression or expansion of the gas will be:

W = _{0}∫^{V} pdV =
p_{0}V_{0} _{0}∫^{V}(V/V_{0})^{-(5/3)}
d(V/V_{0})

= (3/2) p_{0}V_{0} _{0}∫^{T} d(T/T_{0})

= (3/2) (p_{0}V_{0}/T_{0})
(T - T_{0}) = (3/2) N k (T - T_{0})
(97)

**References**

C. Kittel, Elementary Statistical Physics (Wiley, New
York, 1958) pp. 159-165, 192-201.

M. Knudsen, Kinetic Theory of Gases,
(Methuen, London, 1952) pp.33-34.

F.W. Sears, G.L. Salinger,
Thermodynamics, Kinetic Theory, and Statistical Thermodynamics (Addison-Wesley,
3rd ed 1975) pp 331-336, 355-361.

J.M. Ziman, Principles of the Theory of
Solids (Cambridge University Press, 1969) pp 116-123, 179-200.

U Klein, W Vollmann, P Jose Abatti, Contact Potential Differences Measurement:
Short History and Experimental Setup for Classroom Demonstration,
IEEE Trans. Education, Vol. 46, No. 3 (2003) pp 338-344.

On the net: May, 2002. Irreversible thermodynamics added: October 2011.
Charge Transport through a junction and thermoelectric effects added: September 2012.

By the author:

- "Thermoelectric Effect Peltier Seebeck and Thomson",

Abstract: http://urila.tripod.com/Thermoelectric_abstract.htm

Full page: http://urila.tripod.com/Thermoelectric.pdf, February 2014.
- "Osmosis Desalination and Carnot",
http://urila.tripod.com/Osmosis_Carnot.htm, December 2012.
- "Light Scattering",
http://urila.tripod.com/scatter.htm, August (2011).
- "The Sun and the Moon a Riddle in the Sky",
http://urila.tripod.com/moon.htm, July (2011).
- "Osmosis and thermodynamics",
American Journal of Physics, Vol 75 (11), pp. 997-998, November (2007).
- "van't Hoff's Evidence", http://urila.tripod.com/evidence.htm, October (2007).
- "Osmosis and Thermodynamics", http://urila.tripod.com/osmotic.htm,
January (2007).
- "Expansion of an ideal gas", http://urila.tripod.com/expand.htm, December (2002).
- "Optimizing the Efficiency of Reverse Osmosis Seawater Desalination", http://urila.tripod.com/Seawater.htm, May (2002).
- "Boltzmann Transport Equation", http://urila.tripod.com/Boltzmann.htm, May (2002).
- "Energy of Seawater Desalination", http://urila.tripod.com/desalination.htm, April (2000).
- "Avogadro's number atomic and molecular weight",
http://urila.tripod.com/mole.htm, April (2000).
- "Vapor Pressure, Boiling and Freezing Temperatures of a Solution", http://urila.tripod.com/colligative.htm, December (1998).
- "Osmosis Reverse Osmosis and Osmotic Pressure what they are", http://urila.tripod.com/index.htm, February (1998).
- "Calculation of linear coefficients in irreversible processes by kinetic arguments",
American Journal of Physics, Vol 46 (11), pp. 1163-1164, November (1978).
- "Derivation of some basic properties of ideal gases and solutions from processes of elastic collisions",
Journal of Chemical Education, Vol 55 (6), pp. 369-371, June (1978).

Links:

- Thermodynamics Research Laboratory,
http://www.uic.edu/~mansoori/Thermodynamics.Educational.Sites_html
- Thermodynamik - Warmelehre,
http://www.schulphysik.de/thermodyn.html
- The Blind Men and the Elephant
- My Spin on Lunacy
- Five Weeks in a Balloon
- The first man I saw
- "Faster, Faster!"
- Perfection can't be rushed
- The man higher up
- Brains
- The First-Class Passenger
- other