Darwin–Fowler method explained

In statistical mechanics, the Darwin–Fowler method is used for deriving the distribution functions with mean probability. It was developed by Charles Galton Darwin and Ralph H. Fowler in 1922–1923.^[1]

Distribution functions are used in statistical physics to estimate the mean number of particles occupying an energy level (hence also called occupation numbers). These distributions are mostly derived as those numbers for which the system under consideration is in its state of maximum probability. But one really requires average numbers. These average numbers can be obtained by the Darwin–Fowler method. Of course, for systems in the thermodynamic limit (large number of particles), as in statistical mechanics, the results are the same as with maximization.

Darwin–Fowler method

In most texts on statistical mechanics the statistical distribution functions

f in Maxwell–Boltzmann statistics, Bose–Einstein statistics, Fermi–Dirac statistics) are derived by determining those for which the system is in its state of maximum probability. But one really requires those with average or mean probability, although – of course – the results are usually the same for systems with a huge number of elements, as is the case in statistical mechanics. The method for deriving the distribution functions with mean probability has been developed by C. G. Darwin and Fowler^[2] and is therefore known as the Darwin–Fowler method. This method is the most reliable general procedure for deriving statistical distribution functions. Since the method employs a selector variable (a factor introduced for each element to permit a counting procedure) the method is also known as the Darwin–Fowler method of selector variables. Note that a distribution function is not the same as the probability – cf. Maxwell–Boltzmann distribution, Bose–Einstein distribution, Fermi–Dirac distribution. Also note that the distribution function f _{i} which is a measure of the fraction of those states which are actually occupied by elements, is given by f _{i} = n _{i/g} _{i} or n _{i=} f _{ig} _{i} , where g _{i} is the degeneracy of energy level i of energy \varepsilon _{i} and n _{i} is the number of elements occupying this level (e.g. in Fermi–Dirac statistics 0 or 1). Total energy E and total number of elements N are then given by E = \sum _{i} n _{i\varepsilon} _{i} and N = \sum n _{i} .

The Darwin–Fowler method has been treated in the texts of E. Schrödinger,^[3] Fowler^[4] and Fowler and E. A. Guggenheim,^[5] of K. Huang,^[6] and of H. J. W. Müller–Kirsten.^[7] The method is also discussed and used for the derivation of Bose–Einstein condensation in the book of R. B. Dingle.^[8]

Classical statistics

For

N = \sum _{in} _{i} independent elements with n _{i} on level with energy \varepsilon _{i} and E = \sum _{in} _{i\varepsilon} _{i} for a canonical system in a heat bath with temperature T we setZ = \sum ^{- E/kT} _{arrangements e} = \sum _{arrangements\prod} _{iz} ^{n _{i}} ^{i} , z _{i} = e ^{- \varepsilon _{i/kT}} . The average over all arrangements is the mean occupation number(n _{i)} _{av} = (\sum _{jn} _{jZ})/(Z) = z _{j (\partial)/(\partial z _{j})} ln Z. Insert a selector variable \omega by settingZ _{\omega} = \sum \prod _{i(\omega} z ^{n _{i}} _{i)} . In classical statistics the N elements are (a) distinguishable and can be arranged with packets of n _{i} elements on level \varepsilon _{i} whose number is (N!)/(\prod _{in} _{i!}) , so that in this caseZ _{\omega} = N!\sum _{n _{i}} \prod _{i ((\omega z ^{n _{i}} _{i)})/(n _{i!})} . Allowing for (b) the degeneracy g _{i} of level \varepsilon _{i} this expression becomesZ _{\omega} = N!\prod ^{inf ty} _{i =1} \left(\sum _{n _{i=0,1,2,\ldots}} ((\omega z ^{n _{i}} _{i)})/(n _{i!}) \right) ^{g _{i}} = N!e ^{\omega\sum _{ig} _{iz} _{i}} . The selector variable \omega allows one to pick out the coefficient of \omega ^{N} which is Z . ThusZ = \left(\sum _{ig} _{iz} ^{N,} ^{i\right)} and hence(n _{j)} _{av} = z _{j (\partial)/(\partial z _{j})} ln Z = N (g ^{- \varepsilon _{j/kT}} _{je})/(\sum _{ig} ^{- \varepsilon _{i/kT}} ^{ie}) . This result which agrees with the most probable value obtained by maximization does not involve a single approximation and is therefore exact, and thus demonstrates the power of this Darwin–Fowler method.

Quantum statistics

We have as above

Z _{\omega} = \sum\prod (\omega z ^{n _{i}} _{i)} , z ^{- \varepsilon _{i/kT}} _{i = e} , where n _{i} is the number of elements in energy level \varepsilon _{i} . Since in quantum statistics elements are indistinguishable no preliminary calculation of the number of ways of dividing elements into packets n _{1,} n _{2,} n _{3,} ... is required. Therefore the sum \sum refers only to the sum over possible values of n _{i} .

In the case of Fermi–Dirac statistics we have

n _{i=0} or n _{i=1} per state. There are g _{i} states for energy level \varepsilon _{i} .Hence we haveZ _{\omega=(1+\omega} z ^{g _{1}} _{1 )} ( 1+ \omega z ^{g _{2}} _{2 )} … = \prod( 1+ \omega z ^{g _{i}} _{i)} .

In the case of Bose–Einstein statistics we have

n _{i=0,1,2,3,} \ldots inf ty. By the same procedure as before we obtain in the present caseZ _{\omega} = ( 1+ \omega z _{1+(\omega} z ^{2} _{1 )} + (\omega z ^{3} _{1 )} + … ) ^{g _{1}} ( 1+ \omega z _{2} + (\omega z ^{2} _{2 )} + … ) ^{g _{2}} … . But1 + \omega z _{1} + (\omega z ^{2} _{1 )} + … = (1)/(( 1- \omega z _{1)}) . ThereforeZ _{\omega=\prod} _{i(1-\omega} z ^{- g _{i}} _{i)} .

Summarizing both cases and recalling the definition of

Z , we have that Z is the coefficient of \omega ^{N} inZ _{\omega=\prod} _{i(1\pm} \omega z ^{\pm g _{i}} _{i)} , where the upper signs apply to Fermi–Dirac statistics, and the lower signs to Bose–Einstein statistics.

Next we have to evaluate the coefficient of

\omega ^{N} in Z _{\omega.} In the case of a function \phi(\omega) which can be expanded as\phi(\omega) = a _{0} + a _{1\omega} + a ^{2} _{2 \omega} + … , the coefficient of \omega ^{N} is, with the help of the residue theorem of Cauchy,a _{N} = (1)/(2 \pi i) \oint (\phi(\omega)d\omega)/(\omega ^{N+1}) . We note that similarly the coefficient Z in the above can be obtained asZ = (1)/(\oint 2 \pi i) (Z _{\omega})/(\omega ^{N+1}) d\omega\equiv (1)/(2 \pi i) \int e ^{f(\omega)} d\omega, wheref(\omega) = \pm\sum _{ig} _{iln} ( 1 \pm \omega z _{i)-(N+1)ln\omega.} Differentiating one obtainsf'(\omega) = (1)/(\omega) \left[\sum _{i (g _{i})/((\omega z ^{-1} _{i)} \pm 1)} - (N +1 )\right], andf''(\omega) = (N +1)/(\omega ^{2}) \mp (1)/(\omega ^{2}) \sum _{i (g _{i})/([(\omega z ^{-1} _{i)} \pm 1 ] ^{2})} . One now evaluates the first and second derivatives of f(\omega) at the stationary point \omega _{0} at which f'(\omega _{0)=0.} . This method of evaluation of Z around the saddle point \omega _{0} is known as the method of steepest descent. One then obtainsZ = (e ^{f(\omega _{0)}})/(\sqrt{ 2 \pi f''(\omega _{0)}) }.We have f'(\omega _{0)} = 0 and hence(N +1 ) = \sum _{i (g _{i})/((\omega _{0 z} ^{-1} ^{i)} \pm 1)} (the +1 being negligible since N is large). We shall see in a moment that this last relation is simply the formulaN = \sum _{in} _{i.} We obtain the mean occupation number (n _{i)} _{av} by evaluating(n _{j)} _{av} = z _{j (d)/(dz _{j})} ln Z = (g _{j})/((\omega _{0 z} ^{-1} ^{j)} \pm 1) = (g _{j})/(e ^{(\varepsilon _{j-\mu)/kT}} \pm 1) , e ^{\mu/kT} = \omega _{0.} This expression gives the mean number of elements of the total of N in the volume V which occupy at temperature T the 1-particle level \varepsilon _{j} with degeneracy g _{j} (see e.g. a priori probability). For the relation to be reliable one should check that higher order contributions are initially decreasing in magnitude so that the expansion around the saddle point does indeed yield an asymptotic expansion.

Further reading

• Book: The Historical Development of Quantum Theory. Mehra. Jagdish . Jagdish Mehra . Rechenberg . Helmut . Helmut Rechenberg . 2000-12-28. Springer Science & Business Media. 9780387951805. en.

Notes and References

1. Web site: Darwin–Fowler method. Encyclopedia of Mathematics. en. 2018-09-27.
2. C. G. . Darwin . R. H. . Fowler . On the partition of energy . Phil. Mag. . 44 . 1922 . 450–479, 823–842 . 10.1080/14786440908565189 .
3. Book: Schrödinger, E. . Statistical Thermodynamics . Cambridge University Press . 1952 .
4. Book: Fowler, R. H. . Statistical Mechanics . Cambridge University Press . 1952 .
5. Book: Fowler, R. H. . E. . Guggenheim . Statistical Thermodynamics . Cambridge University Press . 1960 .
6. Book: Huang, K. . Statistical Mechanics . Wiley . 1963 .
7. Book: Müller–Kirsten, H. J. W. . Basics of Statistical Physics . 2nd . World Scientific . 2013 . 978-981-4449-53-3 .
8. Book: Dingle, R. B. . Asymptotic Expansions: Their Derivation and Interpretation . Academic Press . 1973 . 267–271 . 0-12-216550-0 .