A General Algorithm for Minimal Uncoupled Stress--Strain Calculations of Elastic Constant
Abstract
We have developed a general algorithm for constructing the minimum number of strain configurations required to extract second-order elastic constants from first-principles stress-strain calculations without algebraic coupling between the constants targeted for extraction. The method operates on the sparsity structure of the symmetry-reduced elastic stiffness matrix, represented by an auxiliary binary matrix, and generates admissible strain patterns algorithmically rather than prescribing separate strain sets for individual crystal symmetries. Each independent elastic constant is obtained from a single stress component, while additional stress relations generated by the same deformation provide internal consistency checks. Because the six-dimensional strain space contains only 2^6-1=63 non-empty strain patterns, the minimality of the resulting sets can be verified exactly for every elastic-symmetry type considered. The construction reproduces established high-efficiency strain sets in cases where they are minimal and yields improved sets for symmetry classes in which previously tabulated schemes require additional configurations. Because the formulation is based on stress-strain derivatives, it can also be readily extended to finite temperatures using consistently sampled ensemble-averaged thermodynamic stresses. Explicit strain sets are derived for all elastic-symmetry types spanning the 230 crystallographic space groups, and their accuracy and computational efficiency are benchmarked against established approaches.
Related
- Scalable construction of force-constant basis sets for large-scale anharmonic lattice-dynamical calculations
- Study of gravitational radiation at 𝒥 for the general type-D Einstein-Maxwell solution with cosmological constant
- Thermal equilibrium in general relativity and the special nature of the cosmological constant
- GæniceN: a minimal model to relax magnetization states self-consistently in general artificial spin ices
- General circuit mapping algorithm for neutral atom quantum computers
BibTeX
@misc{simak2026general,
title = {{A General Algorithm for Minimal Uncoupled Stress--Strain Calculations of Elastic Constant}},
author = {S. I. Simak and Utkarsh Singh},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.11733},
eprint = {2610.11733},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.11733},
}
RIS
TY - UNPB
TI - A General Algorithm for Minimal Uncoupled Stress--Strain Calculations of Elastic Constant
AU - S. I. Simak
AU - Utkarsh Singh
PY - 2026
DA - 2026/10/08
PB - arXiv
DO - 10.48550/arxiv.2610.11733
UR - https://arxiv.org/abs/2610.11733
AB - We have developed a general algorithm for constructing the minimum number of strain configurations required to extract second-order elastic constants from first-principles stress-strain calculations without algebraic coupling between the constants targeted for extraction. The method operates on the sparsity structure of the symmetry-reduced elastic stiffness matrix, represented by an auxiliary binary matrix, and generates admissible strain patterns algorithmically rather than prescribing separate strain sets for individual crystal symmetries. Each independent elastic constant is obtained from a single stress component, while additional stress relations generated by the same deformation provide internal consistency checks. Because the six-dimensional strain space contains only 2^6-1=63 non-empty strain patterns, the minimality of the resulting sets can be verified exactly for every elastic-symmetry type considered. The construction reproduces established high-efficiency strain sets in cases where they are minimal and yields improved sets for symmetry classes in which previously tabulated schemes require additional configurations. Because the formulation is based on stress-strain derivatives, it can also be readily extended to finite temperatures using consistently sampled ensemble-averaged thermodynamic stresses. Explicit strain sets are derived for all elastic-symmetry types spanning the 230 crystallographic space groups, and their accuracy and computational efficiency are benchmarked against established approaches.
ER -