3

In classic thermodynamics one can derive the Maxwell Boltzmann statistics by solving a Lagrange multipliers equation. In this process a new parameter $\beta$ is introduced to take account of the total energy constraint,

$$L = \beta(E-\sum_i p_i \epsilon_i)+...$$

and it is later connected to $\frac{1}{k_B T}$ by inserting our knowledge about perfect gases. On a side note, about every book around seems to skip the passage with a frustrating "it can be shown that".

My problem is that, given similar Lagrange multipliers equation, I'm not sure whether one can in general relate $E$ and $\beta$: it is evident that $\beta$ has no physical significance by itself, and it's always bounded to some parameter in the constraints that it is multiplying.

For example, say that I have an array of cells that can be either occupied by 0, 1 or 2 dots. I have some constraints, one of which is the density of the dots $\rho$. The Lagrange multipliers equation has the form

$$L = \gamma(\rho-p_1-2p_2)+...$$

where $p_i$ is the probability of having a cell with $i$ dots. Is there some way to eliminate the $\gamma$ parameter in the final explicit formula for $p_1$, and use $\rho$ instead?

Qmechanic
  • 201,751
  • I more or less explained it all in this post: http://physics.stackexchange.com/questions/231017/is-the-boltzmann-constant-really-that-important/231065#231065 – DanielSank Feb 18 '17 at 04:49

2 Answers2

3

You probably realize that the Lagrange multiplier $\gamma$ that you introduced is nothing but the chemical potential, right?

Coming back to the Lagrange multiplier $\beta$, you have many ways to convince yourself that it is inversely proportional to the temperature. I will just focus on some of them:

  • The first one consists of devising a simple way to measure experimentally $\beta$. The simplest model you can think of is evidently the ideal gas. You can work out easily that the pressure in this case will be inversely proportional to $\beta$. Now, it is logically very important to remember that the canonical ensemble characterized by the set $(\beta, N, V)$ is model independent in the same way that setting the temperature in usual thermodynamics has nothing to do with the system under study. Therefore, if you find it to be an inverse temperature in one case, it has to be it for every case you can imagine.

  • The second one might actually be easier. The Lagrange multiplier method arises in a more global strategy called a statistical inference that has been much popularized by E.T. Jaynes in the second half of the last century. The idea is to start from Shannon's entropy and to maximize it with respect to some constraints that characterize the statistical ensemble. If you replace the probability within Shannon's entropy by the canonical one, you get the expression for the entropy in the canonical ensemble. It reads:

    $S(\beta,N,V) = -k_B \sum_{m} \frac{e^{-\beta \epsilon_m}}{Q(\beta, N,V)}\ln \left( \frac{e^{-\beta \epsilon_m}}{Q(\beta, N,V)} \right)$

    where the sum runs over all the microstates $m$ and $Q$ is the canonical partition function.

    Now you can compute what is the value of $\frac{\partial S}{\partial \epsilon}$ and you will find out that it is $k_B \beta$...that is $1/T$ if you trust thermodynamics.

  • There is actually another easy way that consists in starting from the microcanonical ensemble and define the temperature from there. If you then subdivide a system into a big part and a small part, you will find directly the Boltzmann weight with the temperature of the big system for any microstate of the small part. Now, if you separately model this small system with your Lagrange multiplier you will find the same law but with $\beta$ instead of $1/k_B T$. Now, this is true whatever the system and if we moreover assume that the partition function is a monotonic function of $\beta$ (unless we are at a phase transition), we have to conclude that $\beta = 1/k_B T$.

gatsu
  • 7,202
1

Check this link http://www.irc.leeds.ac.uk/~phy6pdo/teaching/thermo/hw/part5.pdf

The Langrange multiplier is simply d(function to be maximized)/d(constrain function), and it is related to physical quantities in this case. You can read through the document.