Molecular Simulation/Dilute gases
Ideal Gas Law edit
Gasses are often described using ideal gas law, which relates the pressure of a gas to its density with a simple expression.
Where is the density of the gas, is the pressure of the gas, is the Boltzmann constant, equal to , and is the temperature in Kelvin. The expression is derived from approximating a gas as point masses that have only kinetic energy and experience perfectly elastic collisions. Unfortunately, this theory breaks down at densities where intermolecular forces become significant, that is, when the potential energy is non-zero. Ideal gas law is suitable for only very dilute gasses.
Virial Theorem edit
To more accurately describe the properties of dilute gasses the Virial equation of state is used. Virial theorem accounts for the effects of intermolecular forces, through an expansion into higher order functions of the density. Mathematically, this is described using infinite power series, where and are the second and third Virial coefficients.
At low densities the deviations from ideal gas behavior can be sufficiently described in the second Virial coefficient, .
Where is the volume of a gas and is the configurational integral. The configurational integral for the second Virial coefficient is the contribution of every possible pair of positions weighted over its Boltzmann distribution. For the third Virial coefficient the configurational integral would be the contribution of three interacting particles. By extension, the configurational integral for any nth Virial coefficient would be the contribution of n particles interacting. For this reason higher Virial coefficients are quite complicated to derive, fortunately they are only necessary for describing gasses at pressures above 10atm. The configurational integral for two interacting particles is as follows:
Where is the potential energy of the interaction of a single pair of particles, and and are the positions of particles 1 and 2. The second virial coefficient can be written in terms of pairwise intermolecular interaction potential, , if the position of particle 2 is defined relative to the position of particle 1. The distance, , would then be the distance between the two interacting particles. The equation derived from this modification is as follows:
A different is derived for the hard sphere potential and the Lennard-Jones potential.
Hard Sphere Potential edit
The hard sphere model approximates particles as hard spheres that cannot overlap, if the spheres are not overlapping then the potential energy is zero and if they are overlapping then the potential energy is infinitely high. This approximation represents the very strong short range Pauli repulsion forces. The equation for the potential energy is as follows:
Where is the potential energy and is the radius of the hard sphere. Integrating the configurational integral for the hard sphere potential gives, , as the second Virial coefficient. This model is crude and only accounts for repulsive forces, a slightly more accurate model would be the Lennard-Jones potential model.
Lennard-Jones Potential edit
The Lennard-Jones potential is a combination of a polynomial repulsion term, , and a London dispersion attractive term, .
The and terms can be expanded, and internal energy, , can be expressed as:
Where is the potential well depth, is the intercept, and is the distance between the particles. The second Virial coefficient derived from the Lennard-Jones potential has no analytical solution and must be solved numerically.
The Lennard-Jones model is a more accurate than the hard sphere model as it accounts for attractive interactions and the repulsive term is more realistic than hard sphere repulsion. That being said, it is still limited in the fact that only London dispersion attractive interactions are considered, making the model applicable to only nobel gasses.
- McQuarrie, D. A. Statistical thermodynamics; University Science Books: Mill Valley, CA, 1997.