I How Do Symmetry Properties Affect Boltzmann Distributions for Bosonic Particles?

  • I
  • Thread starter Thread starter jostpuur
  • Start date Start date
  • Tags Tags
    Boltzmann Bosons
jostpuur
Messages
2,112
Reaction score
19
Introduction and warm up

Suppose H_0 is some Hamilton's operator that has eigenvectors |\psi_n\rangle for n\in\mathbb{N} with some eigenvalues E_n so that
<br /> H_0|\psi_n\rangle = E_n|\psi_n\rangle,\quad \forall n\in\mathbb{N}.<br />

Suppose we define a new Hamilton's operator by setting H=H_0\otimes \textrm{id}+\textrm{id}\otimes H_0. Now one possible set of eigenvectors is given by a product |\psi_n\rangle\otimes |\psi_{n&#039;}\rangle for all (n,n&#039;)\in\mathbb{N}^2, and these have eigenvalues E_n+E_{n&#039;}.

Suppose we insist that only those vectors |\psi\rangle, which satisfy the symmetry property
<br /> (\langle x|\otimes\langle y|)|\psi\rangle = (\langle y|\otimes \langle x|)|\psi\rangle<br />
are allowed. If the original H_0 was interpreted as acting on a state of some particle, now H will be acting an a pair of bosonic particles.

Now the acceptable eigenvectors will be (up to a factor) |\psi_n\rangle\otimes |\psi_{n&#039;}\rangle + |\psi_{n&#039;}\rangle\otimes |\psi_n\rangle.

Suppose we would like to write down a formula for the set of pairs (E,|\psi\rangle) which would contain all the eigenvalues and eigenvectors. There is two obvious ways to write down the set. The simplest formula for the set of pairs is
<br /> \Big\{\big(E_n+E_{n&#039;}, |\psi_n\rangle\otimes |\psi_{n&#039;}\rangle + |\psi_{n&#039;}\rangle\otimes |\psi_n\rangle\big)\;\Big|\; n,n&#039;\in\mathbb{N}\Big\}<br />
However, this formula contains some pairs in redundant way, and we might also write down the set of pairs as
<br /> \Big\{\big(E_n+E_{n&#039;}, |\psi_n\rangle\otimes |\psi_{n&#039;}\rangle + |\psi_{n&#039;}\rangle\otimes |\psi_n\rangle\big)\;\Big|\; n,n&#039;\in\mathbb{N},\; n\leq n&#039;\Big\}<br />
These two formulas define the same set, because the elements of these two sets are the same. As long as we are only interested in the set itself, there is no trouble in sight.

The problems start here

Suppose we would like to apply Boltzmann's distribution to this pair of bosons under the assumption that they interact with warm surroundings of some temperature T. There are two obvious options available.

The first option is that we begin by allowing the domain of the index pair (n,n&#039;) to be the full set \mathbb{N}^2. There does not seem to be anything wrong with this choice alone. Once this decision has been made, it seems reasonable to define, perhaps axiomatically, the probability distribution to be proportional to the function
<br /> (n,n&#039;)\mapsto e^{-\frac{E_n+E_{n&#039;}}{T}},\quad\forall n,n&#039;\in\mathbb{N}<br />
Suppose we regret the choice of domain, and would like to shrink it with the additional costraint n\leq n&#039;. In the shrunk domain the new probability distribution, which would be equivalent to the just axiomatically defined one, would be proportional to
<br /> (n,n&#039;)\mapsto \left\{\begin{array}{ll}<br /> e^{-\frac{2E_n}{T}},\quad &amp; n=n&#039; \\<br /> 2e^{-\frac{E_n+E_{n&#039;}}{T}},\quad &amp; n&lt;n&#039; \\<br /> \end{array}\right.<br />

The second option is that we begin by allowing the domain of the index pair (n,n&#039;) to be the set \{(n,n&#039;)\;|\;n,n&#039;\in\mathbb{N},\; n\leq n&#039;\}. There does not yet seem to be anything wrong with this choice either. Once this decision has been made, it seems reasonable to define, perhaps axiomatically, the probability distribution to be proportional to the function
<br /> (n,n&#039;)\mapsto e^{-\frac{E_n+E_{n&#039;}}{T}},\quad\forall n,n&#039;\in\mathbb{N},\; n\leq n&#039;<br />
Suppose we regret the choice of domain, and would like to extend it to the full set \mathbb{N}^2. In the extended domain the new probability distribution, which would be equivalent to the just axiomatically defined one, would be proportional to
<br /> (n,n&#039;)\mapsto \left\{\begin{array}{ll}<br /> e^{-\frac{2E_n}{T}},\quad &amp; n=n&#039; \\<br /> \frac{1}{2}e^{-\frac{E_n+E_{n&#039;}}{T}},\quad &amp; n\neq n&#039; \\<br /> \end{array}\right.<br />

The question: Now we have two different Boltzmann's distributions, which both look reasonable, but since they are different, they cannot be both right. So is one of these right, and the other one wrong? Which way around would be the correct answer?
 
Last edited:
Physics news on Phys.org
You have to use the proper partition function. In the first case, the partition function will have factor 1/2 to account for the over-counting of indistinguishable cases.
 
While using the domain \mathbb{N}^2, in the first case the partition function will be
<br /> Z(T) = \sum_{n,n&#039;=0}^{\infty} e^{-\frac{E_n+E_{n&#039;}}{T}}<br />
and the probability distribution will be
<br /> p(n,n&#039;) = \frac{1}{Z(T)}e^{-\frac{E_n+E_{n&#039;}}{T}},\quad n,n&#039;\in\mathbb{N}<br />
Did you mean to claim that the partition function should be something different, with some factor \frac{1}{2} somewhere? I don't believe that that would be right, because bringing such factors in would imply that
<br /> \sum_{n,n&#039;=0}^{\infty}p(n,n&#039;)=1<br />
would no longer hold.

Equivalently, while using the domain \{(n,n&#039;)\;|\;n,n&#039;\in\mathbb{N},n\leq n&#039;\}, in the first case the partition function will be
<br /> Z(T) = \sum_{n=0}^{\infty} e^{-\frac{2E_n}{T}} + 2\underset{n&lt;n&#039;}{\sum_{n,n&#039;=0}^{\infty}} e^{-\frac{E_n+E_{n&#039;}}{T}}<br />

In the second case, which will not be equivalent to the first one, while using the domain \{(n,n&#039;)\;|\;n,n&#039;\in\mathbb{N},n\leq n&#039;\}, the partition function will be
<br /> Z(T) = \underset{n\leq n&#039;}{\sum_{n,n&#039;=0}^{\infty}} e^{-\frac{E_n+E_{n&#039;}}{T}}<br />

Did you mean that this second option looks more right than the first one? The second option has some kind of over-counting removed, and it looks nice, but how do you know that the over-counting was wrong in the first place?
 
Last edited:
jostpuur said:
While using the domain \mathbb{N}^2, in the first case the partition function will be
<br /> Z(T) = \sum_{n,n&#039;=0}^{\infty} e^{-\frac{E_n+E_{n&#039;}}{T}}<br />
The sum in the partition should run over all states. If you can't distinguish between ##(n,n')## and ##(n',n)##, then these two cases correspond to the same state, hence the sum you have there is not correct: you are over-counting.

See, e.g., http://theory.physics.manchester.ac.uk/~judith/stat_therm/node84.html
 
  • Like
Likes Jilang
I explained two different ways to define the Boltzmann's distribution, and put forward a question that which one of them is right. Your answer is the same as that the first option was wrong, and the second option was right.
 
Insights auto threads is broken atm, so I'm manually creating these for new Insight articles. Towards the end of the first lecture for the Qiskit Global Summer School 2025, Foundations of Quantum Mechanics, Olivia Lanes (Global Lead, Content and Education IBM) stated... Source: https://www.physicsforums.com/insights/quantum-entanglement-is-a-kinematic-fact-not-a-dynamical-effect/ by @RUTA
If we release an electron around a positively charged sphere, the initial state of electron is a linear combination of Hydrogen-like states. According to quantum mechanics, evolution of time would not change this initial state because the potential is time independent. However, classically we expect the electron to collide with the sphere. So, it seems that the quantum and classics predict different behaviours!

Similar threads

Back
Top