1-dimensional hard rods

From SklogWiki
Revision as of 14:06, 19 February 2008 by Nice and Tidy (talk | contribs)
Jump to navigation Jump to search

A 1-dimensional system having hard sphere interactions. The statistical mechanics of this system can be solved exactly (see Ref. 1).

Canonical Ensemble: Configuration Integral

Consider a system of length defined in the range Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \left[ 0, L \right] } .

Our aim is to compute the partition function of a system of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \left. N \right. } hard rods of length Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \left. \sigma \right. } .

Model:

  • External Potential; the whole length of the rod must be inside the range:
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_{0}(x_i) = \left\{ \begin{array}{lll} 0 & ; & \sigma/2 < x < L - \sigma/2 \\ \infty &; & {\rm elsewhere}. \end{array} \right. }

where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \left. x_k \right. } is the position of the center of the k-th rod.

Consider that the particles are ordered according to their label: Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle x_0 < x_1 < x_2 < \cdots < x_{N-1} } ; taking into account the pair potential we can write the canonical partition function (configuration integral) of a system of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle N } particles as:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{ Z \left( N,L \right)}{N!} = \int_{\sigma/2}^{L+\sigma/2-N\sigma} d x_0 \int_{x_0+\sigma}^{L+\sigma/2-N\sigma+\sigma} d x_1 \cdots \int_{x_{i-1}+\sigma}^{L+\sigma/2-N\sigma+i \sigma} d x_i \cdots \int_{x_{N-2}+\sigma}^{L+\sigma/2-N\sigma+(N-1)\sigma} d x_{N-1}. }

Variable change: Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \left. \omega_k = x_k - \left(k+\frac{1}{2}\right) \sigma \right. }  ; we get:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{ Z \left( N,L \right)}{N!} = \int_{0}^{L-N\sigma} d \omega_0 \int_{\omega_0}^{L-N\sigma} d \omega_1 \cdots \int_{\omega_{i-1}}^{L-N\sigma} d \omega_i \cdots \int_{\omega_{N-2}}^{L-N\sigma} d \omega_{N-1}. }

Therefore:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{ Z \left( N,L \right)}{N!} = \frac{ (L-N\sigma )^{N} }{N!}. }
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Q(N,L) = \frac{ (L-N \sigma )^N}{\Lambda^N N!}. }

Thermodynamics

Helmholtz energy function

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \left. A(N,L,T) = - k_B T \log Q \right. }

In the thermodynamic limit (i.e. Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle N \rightarrow \infty; L \rightarrow \infty} with Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho = \frac{N}{L} } , remaining finite):

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle A \left( N,L,T \right) = N k_B T \left[ \log \left( \frac{ N \Lambda} { L - N \sigma }\right) - 1 \right]. }

Equation of state

From the basic thermodynamics, the pressure [linear tension in this case] Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \left. p \right. } can be written as:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle p = - \left( \frac{ \partial A}{\partial L} \right)_{N,T} = \frac{ N k_B T}{L - N \sigma}; }
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Z = \frac{p L}{N k_B T} = \frac{1}{ 1 - \eta}, }

where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \eta \equiv \frac{ N \sigma}{L} } ; is the fraction of volume (length) occupied by the rods.

References

  1. Lewi Tonks "The Complete Equation of State of One, Two and Three-Dimensional Gases of Hard Elastic Spheres", Physical Review 50 pp. 955- (1936)
  2. 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)
  3. L. van Hove, "Sur L'intégrale de Configuration Pour Les Systèmes De Particules À Une Dimension", Physica, 16 pp. 137-143 (1950)