Misplaced Pages

Virial stress

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.
Measure of mechanical stress at the atomic scale

In mechanics, virial stress is a measure of stress on an atomic scale for homogeneous systems. The name is derived from Latin vis 'force': "Virial is then derived from Latin as well, stemming from the word virias (plural of vis) meaning forces." The expression of the (local) virial stress can be derived as the functional derivative of the free energy of a molecular system with respect to the deformation tensor.

Volume averaged Definition

The instantaneous volume averaged virial stress is given by

τ i j = 1 Ω ∑ k ∈ Ω [ − m ( k ) ( u i ( k ) − u ¯ i ) ( u j ( k ) − u ¯ j ) + 1 2 ∑ ℓ ∈ Ω ( x i ( ℓ ) − x i ( k ) ) f j ( k ℓ ) ] {\displaystyle \tau _{ij}={\frac {1}{\Omega }}\sum _{k\in \Omega }\left} where

  • k {\displaystyle k} and ℓ {\displaystyle \ell } are atoms in the domain,
  • Ω {\displaystyle \Omega } is the volume of the domain,
  • m ( k ) {\displaystyle m^{(k)}} is the mass of atom k,
  • u i ( k ) {\displaystyle u_{i}^{(k)}} is the i-th component of the velocity of atom k,
  • u ¯ j {\displaystyle {\bar {u}}_{j}} is the j-th component of the average velocity of atoms in the volume,
  • x i ( k ) {\displaystyle x_{i}^{(k)}} is the i-th component of the position of atom k, and
  • f i ( k ℓ ) {\displaystyle f_{i}^{(k\ell )}} is the i-th component of the force applied on atom k by atom ℓ.

At zero kelvin, all velocities are zero so we have τ i j = 1 2 Ω ∑ k , ℓ ∈ Ω ( x i ( ℓ ) − x i ( k ) ) f j ( k ℓ ) {\displaystyle \tau _{ij}={\frac {1}{2\Omega }}\sum _{k,\ell \in \Omega }\left(x_{i}^{(\ell )}-x_{i}^{(k)}\right)f_{j}^{(k\ell )}}

This can be thought of as follows. The τ11 component of stress is the force in the x1-direction divided by the area of a plane perpendicular to that direction. Consider two adjacent volumes separated by such a plane. The 11-component of stress on that interface is the sum of all pairwise forces between atoms on the two sides.

The volume averaged virial stress is then the ensemble average of the instantaneous volume averaged virial stress.

In a three dimensional, isotropic system, at equilibrium the "instantaneous" atomic pressure is usually defined as the average over the diagonals of the negative stress tensor:

P a t = − 1 3 T r ( τ ) . {\displaystyle {\mathcal {P}}_{at}=-{\frac {1}{3}}Tr(\tau ).}

The pressure then is the ensemble average of the instantaneous pressure P a t = ⟨ P a t ⟩ . {\displaystyle P_{at}=\langle {\mathcal {P}}_{at}\rangle .} This pressure is the average pressure in the volume Ω.

Equivalent Definition

It's worth noting that some articles and textbook use a slightly different but equivalent version of the equation

τ i j = 1 Ω ∑ k ∈ Ω [ − m ( k ) ( u i ( k ) − u ¯ i ) ( u j ( k ) − u ¯ j ) − 1 2 ∑ ℓ ∈ Ω x i ( k ℓ ) f j ( k ℓ ) ] {\displaystyle \tau _{ij}={\frac {1}{\Omega }}\sum _{k\in \Omega }\left}

where x i ( k ℓ ) {\displaystyle x_{i}^{(k\ell )}} is the i-th component of the vector oriented from the ℓ-th atoms to the k-th calculated via the difference

x i k ℓ = x i ( k ) − x i ( ℓ ) {\displaystyle x_{i}^{k\ell }=x_{i}^{(k)}-x_{i}^{(\ell )}}

Both equation being strictly equivalent, the definition of the vector can still lead to confusion.

Derivation

The virial pressure can be derived, using the virial theorem and splitting forces between particles and the container or, alternatively, via direct application of the defining equation P = − ∂ F ( N , V , T ) ∂ V {\displaystyle P=-{\tfrac {\partial F(N,V,T)}{\partial V}}} and using scaled coordinates in the calculation.

Inhomogeneous Systems

If the system is not homogeneous in a given volume the above (volume averaged) pressure is not a good measure for the pressure. In inhomogeneous systems the pressure depends on the position and orientation of the surface on which the pressure acts. Therefore, in inhomogeneous systems a definition of a local pressure is needed. As a general example for a system with inhomogeneous pressure you can think of the pressure in the atmosphere of the earth which varies with height.

Instantaneous local virial stress

The (local) instantaneous virial stress is given by:

τ a b ( r → ) = − ∑ i = 1 N δ ( r → − r → ( i ) ) [ m ( i ) u a ( i ) u b ( i ) + 1 2 ∑ j = 1 , j ≠ i N ( r → ( i ) − r → ( j ) ) a f → b ( i j ) ] , {\displaystyle \tau _{ab}({\vec {r}})=-\sum _{i=1}^{N}\delta \left({\vec {r}}-{\vec {r}}^{(i)}\right)\left,}

Measuring the virial pressure in molecular simulations

The virial pressure can be measured via the formulas above or using volume rescaling trial moves.

See also

References

  1. Wasalwar, Yash (May 23, 2023). "Spotlight: the Virial Theorem". Medium. Archived from the original on February 3, 2024.
  2. ^ Morante, S., G. C. Rossi, and M. Testa. "The stress tensor of a molecular system: An exercise in statistical mechanics." The Journal of chemical physics 125.3 (2006): 034101, http://aip.scitation.org/doi/abs/10.1063/1.2214719.
  3. ^ Allen, MP; Tildesley, DJ (1991). Clarendon Press (ed.). Computer Simulations of Liquids. Oxford. pp. 46–50.{{cite book}}: CS1 maint: location missing publisher (link)
  4. Navet, M.; Jamin, E.; Feix, M. R. (1980-02-01). "" Virial " pressure of the classical one-component plasma". Journal de Physique Lettres. 41 (3): 69–73. doi:10.1051/jphyslet:0198000410306900. ISSN 0302-072X. S2CID 122678419.
  5. Loison, Claire (2005). Numerical Simulations of a Smectic Lamellar Phase of Amphiphilic Molecules. Cuvillier Verlag. ISBN 978-3-86537-655-8.
  6. Miguel, Enrique de; Jackson, George (2006-10-30). "The nature of the calculation of the pressure in molecular simulations of continuous models from volume perturbations". The Journal of Chemical Physics. 125 (16): 164109. Bibcode:2006JChPh.125p4109D. doi:10.1063/1.2363381. hdl:10272/9584. ISSN 0021-9606. PMID 17092065.

External links

Category: