Chemical pressure and vacancies in crystals under nonhydrostatic stress
Abstract
Inhomogeneous stress is a driving force for diffusion of vacancies in crystals. The other way around a non-uniform distribution of vacancies induces stress. The accepted theory of composition-stress coupling in crystals is Larché-Cahn (LC) theory. Composition is defined in terms of the population of lattice sites with a limit of one-particle per site. Coupling is modelled by adding a compositional strain term to the elastic strain in the constitutive relation for stress. An alternative mechanism, proposed here, is letting the site binding energy vary with spatial (deformed) density. This generates a chemical pressure adding to the elastic stress. The governing equations for this model are derived using non-equilibrium continuum thermodynamic methods. The gradient of the chemical pressure has a dual function acting both as the drift force for migration and as an effective internal one body force in the Cauchy equation for the elastic stress. The practical evaluation presented in the paper is restricted to equilibrium properties. We examine the effect of externally applied non-hydrostatic stress, either in the form of surface tractions or a one-body force density (gravitation). The model system is a one component crystal with a fixed number of lattice sites. The number of particles can be variable but is always smaller than the number of lattice sites. The linear elastic response is modelled by the standard Lamé stress tensor. The results for systems deformed by surface tractions are in qualitative agreement with LC theory allowing for differences in the expression for the elastic moduli. Deviations are more serious for the crystal deformed by gravitation.
Related
- Large-Area SnS Crystals by Controlled Sn-S Chemical Vapor Deposition: Growth Optimization, Morphology, and Raman Characterization
- Enabling the Ambient Pressure Growth of ScB2 Crystals for AlGaN Power Electronics
- The electromagnetic stress tensor in cubic crystals and amorphous solids
- Stress-enhanced diffusion in Yukawa crystals
- Enhanced Long-Wavelength Fluctuations and Interaction-Stress Correlations in Active Crystals
BibTeX
@misc{sprik2026chemical,
title = {{Chemical pressure and vacancies in crystals under nonhydrostatic stress}},
author = {Michiel Sprik},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.02435},
eprint = {2609.02435},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.02435},
note = {CC BY-NC-ND 4.0},
}
RIS
TY - UNPB
TI - Chemical pressure and vacancies in crystals under nonhydrostatic stress
AU - Michiel Sprik
PY - 2026
DA - 2026/09/02
PB - arXiv
DO - 10.48550/arxiv.2609.02435
UR - https://arxiv.org/abs/2609.02435
AB - Inhomogeneous stress is a driving force for diffusion of vacancies in crystals. The other way around a non-uniform distribution of vacancies induces stress. The accepted theory of composition-stress coupling in crystals is Larché-Cahn (LC) theory. Composition is defined in terms of the population of lattice sites with a limit of one-particle per site. Coupling is modelled by adding a compositional strain term to the elastic strain in the constitutive relation for stress. An alternative mechanism, proposed here, is letting the site binding energy vary with spatial (deformed) density. This generates a chemical pressure adding to the elastic stress. The governing equations for this model are derived using non-equilibrium continuum thermodynamic methods. The gradient of the chemical pressure has a dual function acting both as the drift force for migration and as an effective internal one body force in the Cauchy equation for the elastic stress. The practical evaluation presented in the paper is restricted to equilibrium properties. We examine the effect of externally applied non-hydrostatic stress, either in the form of surface tractions or a one-body force density (gravitation). The model system is a one component crystal with a fixed number of lattice sites. The number of particles can be variable but is always smaller than the number of lattice sites. The linear elastic response is modelled by the standard Lamé stress tensor. The results for systems deformed by surface tractions are in qualitative agreement with LC theory allowing for differences in the expression for the elastic moduli. Deviations are more serious for the crystal deformed by gravitation.
ER -