Brillouin and Langevin functions
The function is usually applied (see below) in the context where x is a real variable and J is a positive integer or half-integer. In this case, the function varies from -1 to 1, approaching +1 as and -1 as .
The function is best known for arising in the calculation of the magnetization of an ideal paramagnet. In particular, it describes the dependency of the magnetization on the applied magnetic field and the total angular momentum quantum number J of the microscopic magnetic moments of the material. The magnetization is given by:
- 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.
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 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
With N magnetic moments per unit volume, the magnetization density is
where the constants and are defined to be
In the classical limit, the moments can be continuously aligned in the field and can assume all values (). The Brillouin function is then simplified into the Langevin function, named after Paul Langevin:
For small values of x, the Langevin function can be approximated by a truncation of its Taylor series:
An alternative better behaved approximation can be derived from the Lambert's continued fraction expansion of tanh(x):
For small enough x, both approximations are numerically better than a direct evaluation of the actual analytical expression, since the latter suffers from Loss of significance.
The inverse Langevin function L−1(x) is defined on the open interval (−1, 1). For small values of x, 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 x. One popular approximation, valid on the whole range (−1, 1), has been published by A. Cohen:
|Complexity||Optimal approximation||Maximum relative error [%]|
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.
When i.e. when is small, the expression of the magnetization can be approximated by the Curie's law:
where is a constant. One can note that is the effective number of Bohr magnetons.
When , the Brillouin function goes to 1. The magnetization saturates with the magnetic moments completely aligned with the applied field:
- C. Kittel, Introduction to Solid State Physics (8th ed.), pages 303-4 ISBN 978-0-471-41526-8
- Darby, M.I. (1967). "Tables of the Brillouin function and of the related function for the spontaneous magnetization". Br. J. Appl. Phys. 18 (10): 1415–1417. Bibcode:1967BJAP...18.1415D. doi:10.1088/0508-3443/18/10/307.
- Takacs, Jeno (2016). "Approximations for Brillouin and its reverse function". COMPEL. 35 (6): 2095.
- Johal, A. S.; Dunstan, D. J. (2007). "Energy functions for rubber from microscopic potentials" (PDF). Journal of Applied Physics. 101 (8): 084917. Bibcode:2007JAP...101h4917J. doi:10.1063/1.2723870.
- Cohen, A. (1991). "A Padé approximant to the inverse Langevin function". Rheologica Acta. 30 (3): 270–273. doi:10.1007/BF00366640.
- Jedynak, R. (2015). "Approximation of the inverse Langevin function revisited". Rheologica Acta. 54 (1): 29–39. doi:10.1007/s00397-014-0802-2.
- Kröger, M. (2015). "Simple, admissible, and accurate approximants of the inverse Langevin and Brillouin functions, relevant for strong polymer deformations and flows". J Non-Newton Fluid Mech. 223: 77–87. doi:10.1016/j.jnnfm.2015.05.007.
- Petrosyan, R. (2016). "Improved approximations for some polymer extension models". Rheologica Acta. arXiv:1606.02519. doi:10.1007/s00397-016-0977-9.
- Jedynak, R. (2017). "New facts concerning the approximation of the inverse Langevin function". Journal of Non-Newtonian Fluid Mechanics. 249: 8–25. doi:10.1016/j.jnnfm.2017.09.003.
- Jedynak, R. (2018). "A comprehensive study of the mathematical methods used to approximate the inverse Langevin function". Mathematics and Mechanics of Solids: 1–25. doi:10.1177/1081286518811395.
- Benítez, J.M.; Montáns, F.J. (2018). "A simple and efficient numerical procedure to compute the inverse Langevin function with high accuracy". Journal of Non-Newtonian Fluid Mechanics: 153–163. arXiv:1806.08068. doi:10.1016/j.jnnfm.2018.08.011.