Grand Canonical Partition Function for Simple System

PitchAintOne
Messages
1
Reaction score
0

Homework Statement



I would like to calculate the grand canonical partition function (GCPF) for a system in which there are are m lattice sites. A configuration may be specified by the numbers (n_1, n_2, ... , n_m), where n_k = 1 if a particle occupies site k and n_k = 0 if no particle occupies site k. Occupied sites have an associated energy \epsilon (constant) and unoccupied sites have zero associated energy.

Homework Equations



The general form of the GCPF in my book (Chandler) is given like this:

\Xi = \sum_{\nu} e^{ -\beta E_{\nu} + \beta \mu N_{\nu} }

where \nu indicates a summation over all states. (I am confused as to what, exactly, is meant by a "state" in the context of this problem.)

The Attempt at a Solution



For a given state j the number of particles is given by N_j = \sum_{j=1}^{m} n_i (summing over all sites). For the same state j the energy is given by E_j = \epsilon \sum_{i=1}^{m} n_i.

I'm unsure of the correct direction from here. Inserting the expressions for N_j and E_j into \Xi creates a mess of summations. Is that the only way? Is it simplify-able?

I feel that I should be able to calculate the GCPF for just one site and then extend the result to m sites since the sites are independent of one another. Is this possible? If so, how?

Thank you all.
 
Physics news on Phys.org
A state is a set of occupation numbers (n_1,n_2,...,n_m) so when you sum over all possible states, you sum over configurations (0,0,...,0), (1,0,0,...,0), ... , (1,1,...,1). Luckily the system is quite symmetric, so you do not need to sum the states one by one. Instead, you can write the sum over the number of occupied sites like so:

\sum_{\nu} \rightarrow \sum_{n=0}^{m} \omega(n)

where \omega(n) is the density of states, ie. the number of microstates (n1,n2,...,nm) corresponding to the macrostate.
 
Hi, I had an exam and I completely messed up a problem. Especially one part which was necessary for the rest of the problem. Basically, I have a wormhole metric: $$(ds)^2 = -(dt)^2 + (dr)^2 + (r^2 + b^2)( (d\theta)^2 + sin^2 \theta (d\phi)^2 )$$ Where ##b=1## with an orbit only in the equatorial plane. We also know from the question that the orbit must satisfy this relationship: $$\varepsilon = \frac{1}{2} (\frac{dr}{d\tau})^2 + V_{eff}(r)$$ Ultimately, I was tasked to find the initial...
The value of H equals ## 10^{3}## in natural units, According to : https://en.wikipedia.org/wiki/Natural_units, ## t \sim 10^{-21} sec = 10^{21} Hz ##, and since ## \text{GeV} \sim 10^{24} \text{Hz } ##, ## GeV \sim 10^{24} \times 10^{-21} = 10^3 ## in natural units. So is this conversion correct? Also in the above formula, can I convert H to that natural units , since it’s a constant, while keeping k in Hz ?
Back
Top