Peeter Joot's (OLD) Blog.

Math, physics, perl, and programming obscurity.

Posts Tagged ‘reservoir’

An updated compilation of notes, for ‘PHY452H1S Basic Statistical Mechanics’, Taught by Prof. Arun Paramekanti

Posted by peeterjoot on March 27, 2013

Here’s my second update of my notes compilation for this course, including all of the following:

March 27, 2013 Fermi gas

March 26, 2013 Fermi gas thermodynamics

March 26, 2013 Fermi gas thermodynamics

March 23, 2013 Relativisitic generalization of statistical mechanics

March 21, 2013 Kittel Zipper problem

March 18, 2013 Pathria chapter 4 diatomic molecule problem

March 17, 2013 Gibbs sum for a two level system

March 16, 2013 open system variance of N

March 16, 2013 probability forms of entropy

March 14, 2013 Grand Canonical/Fermion-Bosons

March 13, 2013 Quantum anharmonic oscillator

March 12, 2013 Grand canonical ensemble

March 11, 2013 Heat capacity of perturbed harmonic oscillator

March 10, 2013 Langevin small approximation

March 10, 2013 Addition of two one half spins

March 10, 2013 Midterm II reflection

March 07, 2013 Thermodynamic identities

March 06, 2013 Temperature

March 05, 2013 Interacting spin

plus everything detailed in the description of my first update and before.

Posted in Math and Physics Learning. | Tagged: , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , | 1 Comment »

Relativistic generalization of statistical mechanics

Posted by peeterjoot on March 22, 2013

[Click here for a PDF of this post with nicer formatting (especially if my latex to wordpress script has left FORMULA DOES NOT PARSE errors.)]


I was wondering how to generalize the arguments of [1] to relativistic systems. Here’s a bit of blundering through the non-relativistic arguments of that text, tweaking them slightly.

I’m sure this has all been done before, but was a useful exercise to understand the non-relativistic arguments of Pathria better.

Generalizing from energy to four momentum

Generalizing the arguments of section 1.1.

Instead of considering that the total energy of the system is fixed, it makes sense that we’d have to instead consider the total four-momentum of the system fixed, so if we have N particles, we have a total four momentum

\begin{aligned}P = \sum_i n_i P_i = \sum n_i \left( \epsilon_i/c, \mathbf{p}_i \right),\end{aligned} \hspace{\stretch{1}}(1.2.1)

where n_i is the total number of particles with four momentum P_i. We can probably expect that the n_i‘s in this relativistic system will be smaller than those in a non-relativistic system since we have many more states when considering that we can have both specific energies and specific momentum, and the combinatorics of those extra degrees of freedom. However, we’ll still have

\begin{aligned}N = \sum_i n_i.\end{aligned} \hspace{\stretch{1}}(1.2.2)

Only given a specific observer frame can these these four-momentum components \left( \epsilon_i/c, \mathbf{p}_i \right) be expressed explicitly, as in

\begin{aligned}\epsilon_i = \gamma_i m_i c^2\end{aligned} \hspace{\stretch{1}}(1.0.3a)

\begin{aligned}\mathbf{p}_i = \gamma_i m \mathbf{v}_i\end{aligned} \hspace{\stretch{1}}(1.0.3b)

\begin{aligned}\gamma_i = \frac{1}{{\sqrt{1 - \mathbf{v}_i^2/c^2}}},\end{aligned} \hspace{\stretch{1}}(1.0.3c)

where \mathbf{v}_i is the velocity of the particle in that observer frame.

Generalizing the number if microstates, and notion of thermodynamic equilibrium

Generalizing the arguments of section 1.2.

We can still count the number of all possible microstates, but that number, denoted \Omega(N, V, E), for a given total energy needs to be parameterized differently. First off, any given volume is observer dependent, so we likely need to map

\begin{aligned}V \rightarrow \int d^4 x = \int dx^0 \wedge dx^1 \wedge dx^2 \wedge dx^3.\end{aligned} \hspace{\stretch{1}}(1.0.4)

Let’s still call this V, but know that we mean this to be four volume element, bounded in both space and time, referred to a fixed observer’s frame. So, lets write the total number of microstates as

\begin{aligned}\Omega(N, V, P) = \Omega \left( N, \int d^4 x, E/c, P^1, P^2, P^3 \right),\end{aligned} \hspace{\stretch{1}}(1.0.5)

where P = ( E/c, \mathbf{P} ) is the total four momentum of the system. If we have a system subdivided into to two systems in contact as in fig. 1.1, where the two systems have total four momentum P_1 and P_2 respectively.

Fig 1.1: Two physical systems in thermal contact


In the text the total energy of both systems was written

\begin{aligned}E^{(0)} = E_1 + E_2,\end{aligned} \hspace{\stretch{1}}(1.0.6)

so we’ll write

\begin{aligned}{P^{(0)}}^\mu = P_1^\mu + P_2^\mu = \text{constant},\end{aligned} \hspace{\stretch{1}}(1.0.7)

so that the total number of microstates of the combined system is now

\begin{aligned}\Omega^{(0)}(P_1, P_2) = \Omega_1(P_1) \Omega_2(P_2).\end{aligned} \hspace{\stretch{1}}(1.0.8)

As before, if \bar{{P}}^\mu_i denotes an equilibrium value of P_i^\mu, then maximizing eq. 1.0.8 requires all the derivatives (no sum over \mu here)

\begin{aligned}\left({\partial {\Omega_1(P_1)}}/{\partial {P^\mu_1}}\right)_{{P_1 = \bar{{P_1}}}}\Omega_2(\bar{{P}}_2)+\Omega_1(\bar{{P}}_1)\left({\partial {\Omega_2(P_2)}}/{\partial {P^\mu}}\right)_{{P_2 = \bar{{P_2}}}}\times\frac{\partial {P_2^\mu}}{\partial {P_1^\mu}}= 0.\end{aligned} \hspace{\stretch{1}}(1.0.9)

With each of the components of the total four-momentum P^\mu_1 + P^\mu_2 separately constant, we have {\partial {P_2^\mu}}/{\partial {P_1^\mu}} = -1, so that we have

\begin{aligned}\left({\partial {\ln \Omega_1(P_1)}}/{\partial {P^\mu_1}}\right)_{{P_1 = \bar{{P_1}}}}=\left({\partial {\ln \Omega_2(P_2)}}/{\partial {P^\mu}}\right)_{{P_2 = \bar{{P_2}}}},\end{aligned} \hspace{\stretch{1}}(1.0.10)

as before. However, we now have one such identity for each component of the total four momentum P which has been held constant. Let’s now define

\begin{aligned}\beta_\mu \equiv \left({\partial {\ln \Omega(N, V, P)}}/{\partial {P^\mu}}\right)_{{N, V, P = \bar{{P}}}},\end{aligned} \hspace{\stretch{1}}(1.0.11)

Our old scalar temperature is then

\begin{aligned}\beta_0 = c \left({\partial {\ln \Omega(N, V, P)}}/{\partial {E}}\right)_{{N, V, P = \bar{{P}}}} = c \beta = \frac{c}{k_{\mathrm{B}} T},\end{aligned} \hspace{\stretch{1}}(1.0.12)

but now we have three additional such constants to figure out what to do with. A first start would be figuring out how the Boltzmann probabilities should be generalized.

Equilibrium between a system and a heat reservoir

Generalizing the arguments of section 3.1.

As in the text, let’s consider a very large heat reservoir A' and a subsystem A as in fig. 1.2 that has come to a state of mutual equilibrium. This likely needs to be defined as a state in which the four vector \beta_\mu is common, as opposed to just \beta_0 the temperature field being common.

Fig 1.2: A system A immersed in heat reservoir A’


If the four momentum of the heat reservoir is P_r' with P_r for the subsystem, and

\begin{aligned}P_r + P_r' = P^{(0)} = \text{constant}.\end{aligned} \hspace{\stretch{1}}(1.0.13)


\begin{aligned}\Omega'({P^\mu_r}') = \Omega'(P^{(0)} - {P^\mu_r}) \propto P_r,\end{aligned} \hspace{\stretch{1}}(1.0.14)

for the number of microstates in the reservoir, so that a Taylor expansion of the logarithm around P_r' = P^{(0)} (with sums implied) is

\begin{aligned}\ln \Omega'({P^\mu_r}') = \ln \Omega'({P^{(0)}}) +\left({\partial {\ln \Omega'}}/{\partial {{P^\mu}'}}\right)_{{P' = P^{(0)}}} \left( P^{(0)} - P^\mu \right)\approx\text{constant} - \beta_\mu' P^\mu.\end{aligned} \hspace{\stretch{1}}(1.0.15)

Here we’ve inserted the definition of \beta^\mu from eq. 1.0.11, so that at equilibrium, with \beta_\mu' = \beta_\mu, we obtain

\begin{aligned}\Omega'({P^\mu_r}') = \exp\left( - \beta_\mu P^\mu \right)=\exp\left( - \beta E \right)\exp\left( - \beta_1 P^1 \right)\exp\left( - \beta_2 P^3 \right)\exp\left( - \beta_3 P^3 \right).\end{aligned} \hspace{\stretch{1}}(1.0.16)

Next steps

This looks consistent with the outline provided in by Lubos to the stackexchange “is there a relativistic quantum thermodynamics” question. I’m sure it wouldn’t be too hard to find references that explore this, as well as explain why non-relativistic stat mech can be used for photon problems. Further exploration of this should wait until after the studies for this course are done.


[1] RK Pathria. Statistical mechanics. Butterworth Heinemann, Oxford, UK, 1996.

Posted in Math and Physics Learning. | Tagged: , , , , , , , , , , , , , | Leave a Comment »