1-dimensional hard rods

1-dimensional hard rods (sometimes known as a Tonks gas ) consist of non-overlapping line segments of length $\sigma$ who all occupy the same line which has length $L$. One could also think of this model as being a string of hard spheres confined to 1 dimension (not to be confused with 3-dimensional hard rods). The model is given by the intermolecular pair potential:

${\mathrm{\Phi }}_{12}(x_i,x_j)=\{\begin{array}{ccc}0& ;& |x_i-x_j|>\sigma \\ \mathrm{\infty }& ;& |x_i-x_j|<\sigma \end{array}$

where $x_k$ is the position of the center of the k-th rod, along with an external potential. Thus, the Boltzmann factor is

$e_{ij}:=e^{-\beta {\mathrm{\Phi }}_{12}(x_i,x_j)}=\mathrm{\Theta }(|x_i-x_j|-\sigma )=\{\begin{array}{ccc}1& ;& |x_i-x_j|>\sigma \\ 0& ;& |x_i-x_j|<\sigma \end{array}$

The whole length of the rod must be inside the range:

$V_0(x_i)=\{\begin{array}{ccc}0& ;& \sigma /2<x_i<L-\sigma /2\\ \mathrm{\infty }& ;& .\end{array}$

Canonical Ensemble: Configuration Integral

The statistical mechanics of this system can be solved exactly. Consider a system of length $L$ defined in the range $[0,L]$. The aim is to compute the partition function of a system of $N$ hard rods of length $\sigma$. Consider that the particles are ordered according to their label: $x_0<x_1<x_2<\cdots <x_{N-1}$; taking into account the pair potential we can write the canonical partition function of a system of $N$ particles as:

$\begin{array}{rr}\frac{Z(N,L)}{N!}& =\underset{\sigma /2}{\overset{L-\sigma /2}{\int }}d{x}_0\underset{\sigma /2}{\overset{L-\sigma /2}{\int }}d{x}_1\cdots \underset{\sigma /2}{\overset{L-\sigma /2}{\int }}d{x}_{N-1}\prod _{i=1}^{N-1}{e}_{i-1,i}\\ & =\underset{\sigma /2}{\overset{L+\sigma /2-N\sigma }{\int }}d{x}_0\underset{x_0+\sigma }{\overset{L+\sigma /2-N\sigma +\sigma }{\int }}d{x}_1\cdots \underset{x_{i-1}+\sigma }{\overset{L+\sigma /2-N\sigma +i\sigma }{\int }}d{x}_i\cdots \underset{x_{N-2}+\sigma }{\overset{L+\sigma /2-N\sigma +(N-1)\sigma }{\int }}d{x}_{N-1}.\end{array}$

Variable change: ${\omega }_k=x_k-(k+\frac{1}{2})\sigma$ ; we get:

$\begin{array}{rr}\frac{Z(N,L)}{N!}& =\underset{0}{\overset{L-N\sigma }{\int }}d{\omega }_0\underset{{\omega }_0}{\overset{L-N\sigma }{\int }}d{\omega }_1\cdots \underset{{\omega }_{i-1}}{\overset{L-N\sigma }{\int }}d{\omega }_i\cdots \underset{{\omega }_{N-2}}{\overset{L-N\sigma }{\int }}d{\omega }_{N-1}\\ & =\underset{0}{\overset{L-N\sigma }{\int }}d{\omega }_0\cdots \underset{{\omega }_{i-1}}{\overset{L-N\sigma }{\int }}d{\omega }_i\frac{(L-N\sigma -{\omega }_i{)}^{N-1-i}}{(N-1-i)!}=\underset{0}{\overset{L-N\sigma }{\int }}d{\omega }_0\frac{(L-N\sigma -{\omega }_0{)}^{N-1}}{(N-1)!}\end{array}$

Therefore:

$\frac{Z(N,L)}{N!}=\frac{(L-N\sigma {)}^N}{N!}.$

$Q(N,L)=\frac{(L-N\sigma {)}^N}{{\mathrm{\Lambda }}^{N}N!}.$

Thermodynamics

Helmholtz energy function

$A(N,L,T)=-k_{B}T\log Q$

In the thermodynamic limit (i.e. $N\to \mathrm{\infty };L\to \mathrm{\infty }$ with $\rho =\frac{N}{L}$, remaining finite):

$A(N,L,T)=N{k}_{B}T[\log (\frac{N\mathrm{\Lambda }}{L-N\sigma })-1].$

Equation of state

Using the thermodynamic relations, the pressure (linear tension in this case) $p$ can be written as:

$p=-{(\frac{\partial A}{\partial L})}_{N,T}=\frac{N{k}_{B}T}{L-N\sigma };$

The compressibility factor is

$Z=\frac{pL}{N{k}_{B}T}=\frac{1}{1-\eta }=\underset{Z_{\mathrm{i}\mathrm{d}}}{\underbrace{1}}+\underset{Z_{\mathrm{e}\mathrm{x}}}{\underbrace{\frac{\eta }{1-\eta }}},$

where $\eta \equiv \frac{N\sigma }{L}$; is the fraction of volume (i.e. length) occupied by the rods. 'id' labels the ideal and 'ex' the excess part.

It was shown by van Hove that there is no fluid-solid phase transition for this system (hence the designation Tonks gas).

Chemical potential

The chemical potential is given by

$\mu ={(\frac{\partial A}{\partial N})}_{L,T}=k_{B}T(\ln \frac{\rho \mathrm{\Lambda }}{1-\rho \sigma }+\frac{\rho \sigma }{1-\rho \sigma })=k_{B}T(\ln \frac{\rho \mathrm{\Lambda }}{1-\eta }+\frac{\eta }{1-\eta })$

with ideal and excess part separated:

$\beta \mu =\underset{\beta {\mu }_{\mathrm{i}\mathrm{d}}}{\underbrace{\ln (\rho \mathrm{\Lambda })}}+\underset{\beta {\mu }_{\mathrm{e}\mathrm{x}}}{\underbrace{\ln \frac{1}{1-\eta }+\frac{\eta }{1-\eta }}}$

Isobaric ensemble: an alternative derivation

Adapted from Reference . If the rods are ordered according to their label: $x_0<x_1<x_2<\cdots <x_{N-1}$ the canonical partition function can also be written as:

$Z=\underset{0}{\overset{x_1}{\int }}d{x}_0\underset{0}{\overset{x_2}{\int }}d{x}_1\cdots \underset{0}{\overset{L}{\int }}d{x}_{N-1}f(x_1-x_0)f(x_2-x_1)\cdots f(x_0+L-x_{N-1}),$

where $N!$ does not appear one would have $N!$ analogous expressions by permuting the label of the (distinguishable) rods. $f(x)$ is the Boltzmann factor of the hard rods, which is $0$ if $x<\sigma$ and $1$ otherwise.

A variable change to the distances between rods: $y_k=x_k-x_{k-1}$ results in

$Z=\underset{0}{\overset{\mathrm{\infty }}{\int }}d{y}_0\underset{0}{\overset{\mathrm{\infty }}{\int }}d{y}_1\cdots \underset{0}{\overset{\mathrm{\infty }}{\int }}d{y}_{N-1}f(y_0)f(y_1)\cdots f(y_{N-1})\delta (\sum _{i=0}^{N-1}{y}_i-L):$

the distances can take any value as long as they are not below $\sigma$ (as enforced by $f(y)$) and as long as they add up to $L$ (as enforced by the Dirac delta). Writing the later as the inverse Laplace transform of an exponential:

$Z=\underset{0}{\overset{\mathrm{\infty }}{\int }}d{y}_0\underset{0}{\overset{\mathrm{\infty }}{\int }}d{y}_1\cdots \underset{0}{\overset{\mathrm{\infty }}{\int }}d{y}_{N-1}f(y_0)f(y_1)\cdots f(y_{N-1})\frac{1}{2\pi i}\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}ds\exp [-s(\sum _{i=0}^{N-1}{y}_i-L)].$

Exchanging integrals and expanding the exponential the $N$ integrals decouple:

$Z=\frac{1}{2\pi i}\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}ds{e}^{Ls}{\{\underset{0}{\overset{\mathrm{\infty }}{\int }}dyf(y)e^{-sy}\}}^N.$

We may proceed to invert the Laplace transform (e.g. by means of the residues theorem), but this is not needed: we see our configuration integral is the inverse Laplace transform of another one,

$Z^{\prime }(s)={\{\underset{0}{\overset{\mathrm{\infty }}{\int }}dyf(y)e^{-sy}\}}^N,$

so that

$Z^{\prime }(s)=\underset{0}{\overset{\mathrm{\infty }}{\int }}ds{e}^{Ls}Z(L).$

This is precisely the transformation from the configuration integral in the canonical ($N,T,L$) ensemble to the isobaric ($N,T,p$) one, if one identifies $s=p/kT$. Therefore, the Gibbs energy function is simply $G=-kT\log Z^{\prime }(p/kT)$, which easily evaluated to be $G=kTN\log (p/kT)+p\sigma N$. The chemical potential is $\mu =G/N$, and by means of thermodynamic identities such as $\rho =\partial p/\partial \mu$ one arrives at the same equation of state as the one given above.

References

Related reading

• L. van Hove "Quelques Propriétés Générales De L'intégrale De Configuration D'un Système De Particules Avec Interaction", Physica, 15 pp. 951-961 (1949)
• Zevi W. Salsburg, Robert W. Zwanzig, and John G. Kirkwood "Molecular Distribution Functions in a One-Dimensional Fluid", Journal of Chemical Physics 21 pp. 1098-1107 (1953)
• Robert L. Sells, C. W. Harris, and Eugene Guth "The Pair Distribution Function for a One-Dimensional Gas", Journal of Chemical Physics 21 pp. 1422-1423 (1953)
• Donald Koppel "Partition Function for a Generalized Tonks' Gas", Physics of Fluids 6 609 (1963)
• J. L. Lebowitz, J. K. Percus and J. Sykes "Time Evolution of the Total Distribution Function of a One-Dimensional System of Hard Rods", Physical Review 171 pp. 224-235 (1968)
• Gerardo Soto-Campos, David S. Corti, and Howard Reiss "A small system grand ensemble method for the study of hard-particle systems", Journal of Chemical Physics 108 pp. 2563-2570 (1998)
• Paolo V. Giaquinta "Entropy and Ordering of Hard Rods in One Dimension", Entropy 10 pp. 248-260 (2008)