Brillouin and Langevin functions. The Brillouin and Langevin functions are a pair of special functions that appear when studying an idealized paramagnetic material in statistical mechanics. These functions are named after French physicists Paul Langevin and Léon Brillouin who contributed to the microscopic understanding of magnetic properties of matter.

The Langevin function is derived using statistical mechanics, and describes how magnetic dipoles are alignment by an applied field. The Brillouin function was developed later to give an explanation that considers quantum physics. The Langevin function could then be a seen as a special case of the more general Brillouin function if the quantum number  would be infinite ().

Brillouin function for paramagnetism  

The Brillouin function arises when studying magnetization of an ideal paramagnet. In particular, it describes the dependency of the magnetization  on the applied magnetic field , defined by the following equation:

The function  is usually applied in the context where  is a real variable and a function of the applied field . In this case, the function varies from -1 to 1, approaching +1 as  and -1 as .

The total angular momentum quantum number  is a positive integer or half-integer. Considering the microscopic magnetic moments of the material. The magnetization is given by:

where
 is the number of atoms per unit volume,
 the g-factor,
 the Bohr magneton,
 is the ratio of the Zeeman energy of the magnetic moment in the external field  to the thermal energy :

 is the Boltzmann constant and  the temperature.

Note that in the SI system of units  given in Tesla stands for the magnetic field, , where  is the auxiliary magnetic field given in A/m and  is the permeability of vacuum.

{| class="toccolours collapsible collapsed" width="80%" style="text-align:left"
!Click "show" to see a derivation of this law:
|-
|A derivation of this law describing the magnetization of an ideal paramagnet is as follows. Let z be the direction of the magnetic field. The z-component of the angular momentum of each magnetic moment (a.k.a. the azimuthal quantum number) can take on one of the 2J+1 possible values -J,-J+1,...,+J. Each of these has a different energy, due to the external field B: The energy associated with quantum number m is

(where g is the g-factor, μB is the Bohr magneton, and x is as defined in the text above). The relative probability of each of these is given by the Boltzmann factor:

where Z (the partition function) is a normalization constant such that the probabilities sum to unity. Calculating Z, the result is:
.
All told, the expectation value of the azimuthal quantum number m is
.
The denominator is a geometric series and the numerator is a type of arithmetico–geometric series, so the series can be explicitly summed. After some algebra, the result turns out to be

With N magnetic moments per unit volume, the magnetization density is
.
|}

High-field saturation limit  

When , the Brillouin function goes to 1. The magnetization saturates with the magnetic moments completely aligned with the applied field:

Low field/high-temperature limit  

For low fields the curve appears almost linear, and could be replaced by a linear slope as in Curie's law of paramagnetism.  When  (i.e. when  is small) the expression of the magnetization can be approximated by:

and equivalent to Curie's law with the constant given by 

Using  as the effective number of Bohr magnetons.

Note that this is only valid for low fields in paramagnetism. Ferromagnetic materials still has a spontaneous magnetization at low fields (below the Curie-temperature), and the susceptibility must then instead be explained by Curie–Weiss law.

Two-state case (spin-1/2)  

The most simple case of the Brillouin function would be the case of , when the function simplifies to the shape of a tanh-function. Then written as

 

This could be linked to Ising's model, for a case with two possible spins: either up or down. Directed in parallel or antiparallel to the applied field.

This is then equivalent to a 2-state particle: it may either align its magnetic moment with the magnetic field or against it. So the only possible values of magnetic moment are then  and . If so, then such a particle has only two possible energies,  when it is aligned with the field and  when it is oriented opposite to the field.

Langevin function for classical paramagnetism  

The Langevin function () was named after Paul Langevin who published two papers with this function in 1905 to describe paramagnetism by statistical mechanics. Written as:

It could be derivated by describing how magnetic moments are aligned by a magnetic field, considering the statistical thermodynamics. One derivation could be seen here:

{| class="toccolours collapsible collapsed" width="80%" style="text-align:left"
!Click "show" to see a derivation of this function:
|-
| The energy for each magnetic moment will be

where  is the angle between the magnetic moment and
the magnetic field (which we take to be pointing in the 
coordinate.) The corresponding partition function is

Here we could simplify to replace  
We see there is no dependence on the  angle, and also we can
change variables to  to obtain

Now, the expected value of the  component of the magnetization (the other two are seen to be null (due to integration over ), as they should) will be given by

To simplify the calculation, we see this can be written as a differentiation of :

Carrying out the derivation we find

where  is the Langevin function:

|}

The Langevin function can also be derived as the  classical limit of the Brillouin function, if the magnetic moments can be continuously aligned in the field and the quantum number  would be able to assume all values (). The Brillouin function is then simplified into the langevin function.

Classical or quantum approach? 

Langevin function is often seen as the classical theory of paramagnetism, while the Brillouin function is the quantum theory of paramagnetism. When Langevin published the theory paramagnetism in 1905 it was before the adoption of quantum physics. Meaning that Langevin only used concepts of classical physics.

Niels Bohr showed in his thesis that classical statistical mechanics can not be used to explain paramagnetism, and that quantum theory has to be used. This would later be known as the Bohr–Van Leeuwen theorem. The magnetic moment would later be explained in quantum theory by the Bohr magneton (), which is used in the Brillouin function.

It could be noted that there is a difference in the approaches of Langevin and Bohr, since Langevin assumes a magnetic polarization  as the basis for the derivation, while Bohr start the derivation from motions of electrons and a model of an atom. Langevin is still assuming a fix magnetic dipole. This could be expressed as by J. H. Van Vleck:
"When Langevin assumed that the magnetic moment of the atom or molecule had a fixed value , he was quantizing the system without realizing it ". This makes the Langevin function to be in the borderland between classical statisitcal mechanics and quantum theory (as either semi-classical or semi-quantum).

Langevin function for electric polarization

The Langevin function could also be used to describe electric polarization, in the specific case when the polarization is explained by orientation of (electrically polarized) dipoles. So that the electric polarization is given by:

but here for an electric dipole moment  and an electric field  (insead of the magntic equivalents), that is

Simplified functions 
For small values of , the Langevin function can be approximated by a truncation of its Taylor series:

The first term of this series expansion is equivalent to Curie's law, when writing it as

 

An alternative, better behaved approximation can be derived from the
Lambert's continued fraction expansion of :

For small enough , both approximations are numerically better than a direct evaluation of the actual analytical expression, since the latter suffers from catastrophic cancellation for  where .

Inverse Langevin function 
The inverse Langevin function () is without an explicit analytical form, but there exist several approximations.

The inverse Langevin function  is defined on the open interval (−1, 1). For small values of , it can be approximated by a truncation of its Taylor series

and by the Padé approximant

Since this function has no closed form, it is useful to have approximations valid for arbitrary values of . One popular approximation, valid on the whole range (−1, 1), has been published by A. Cohen:

This has a maximum relative error of 4.9% at the vicinity of . Greater accuracy can be achieved by using the formula given by R. Jedynak:

valid for . The maximal relative error for this approximation is 1.5% at the vicinity of x = 0.85. Even greater accuracy can be achieved by using the formula given by M. Kröger:

The maximal relative error for this approximation is less than 0.28%. More accurate approximation was reported by R. Petrosyan:

valid for . The maximal relative error for the above formula is less than 0.18%.

New approximation given by R. Jedynak, is the best reported approximant at complexity 11:

valid for . Its maximum relative error is less than 0.076%.

Current state-of-the-art diagram of the approximants to the inverse Langevin function 
presents the figure below. It is valid for the rational/Padé approximants,

A recently published paper by R. Jedynak, provides a series of the optimal approximants to the inverse Langevin function. The table below reports the results with correct asymptotic behaviors,.

Comparison of relative errors for the different optimal rational approximations, which were computed with constraints (Appendix 8 Table 1)

Also recently, an efficient near-machine precision approximant, based on spline interpolations, has been proposed by Benítez and Montáns, where Matlab code is also given to generate the spline-based approximant and to compare many of the previously proposed approximants in all the function domain.

Inverse Brillouin function 

Approximations could also be used to express the inverse Brillouin function () . Takacs proposed the following approximation to the inverse of the Brillouin function:

where the constants  and  are defined to be

See also 

Paramagnetism
Bohr magneton
Paul Langevin

References 

Eponymous laws of physics
Magnetism