scientific commentaries
On symmetries of higher-order elastic constants
aTerminal Effects Division, DEVCOM ARL, Aberdeen Proving Ground, MD, 21005-5066, USA
*Correspondence e-mail: john.d.clayton1.civ@army.mil
In elastic crystals, a hyperelastic description is conventionally assumed, and the strain energy potential is idealized as a Taylor-series expansion in strain about an unstrained reference state. Coefficients of quadratic terms are second-order or linear elastic constants. Coefficients of higher-order terms are elastic constants of third order, fourth order, and so on. Recently published work by Telyatnik [Acta Cryst. (2024), A80, 394–404] extends prior knowledge of symmetry properties for anisotropic elastic constants of single crystals, as well as transversely isotropic and isotropic solids, to terms up to sixth order. Effective elastic constants for polycrystalline aggregates, with possible anisotropy, were reported by Telyatnik, in the same article, to the same order. A terse summary of nonlinear crystal elasticity and independent elastic constants of orders two and three are given in this commentary for context. Methods and results of Telyatnik, anticipated to be of great utility to crystal elasticity research, are then highlighted.
Keywords: nonlinear elasticity; crystal symmetry; continuum mechanics.
1. Nonlinear elasticity of crystals
Linear elasticity accurately describes the mechanical response of elastic solids when deformations are small and the stress–strain response is linear, meaning Hooke's law applies. Linear elasticity theory (Voigt, 1910; Love, 1927; Hearmon, 1946) can be successfully applied to many, if not most, problems in structural mechanics. A nonlinear theory, on the other hand, is needed for accurate descriptions of mechanics of solids when deformations are large or when linearity breaks down.
All known solids are ultimately nonlinear. No finite volume of any real material can be compressed to a point of infinitesimal size: its bulk modulus must eventually increase with decreasing volume to prevent this. Nonlinear elasticity can describe wave propagation in pre-stressed crystals (Thurston & Brugger, 1964; Thurston et al., 1966; Thurston, 1974) and short-range core effects from lattice defects (Teodosiu, 1982). At regimes departing more from linearity (Chang & Barsch, 1967), or in shock compression, elastic constants up to order four have been measured (Fowles, 1967; Graham, 1972). In strong crystals such as quartz, sapphire and diamond, large uniaxial compressive deformations can be reached before inelasticity from dislocation motion, deformation or fracture ensues (Clayton, 2019). Chen et al. (2020) showed the importance of anisotropic constants up to fifth order on the stability of silicon.
Symmetries of second-order constants for anisotropic elasticity of crystals have been known for over a century (Voigt, 1910; Hearmon, 1946). Independent third-order constants have been known for all crystal classes since works by Fumi (1951) and Hearmon (1953), and those of fourth-order since the work of Brendel (1979). Measurements of third-order constants often rely on sound speeds in pre-stressed crystals (Thurston & Brugger, 1964; Brugger, 1965; Thurston, 1974). Contemporary first-principles atomic simulations involving density functional theory (DFT) have been used to predict constants up to orders four and five (Chen et al., 2020; Pandit & Bongiorno, 2023). Since elastic constants are ultimately related to interatomic forces (Born & Huang, 1954; Wallace, 1972), characteristics of elastic constants give insight into atomic-scale physics, and vice versa. Higher-order constants are associated with anharmonicity (Hiki & Granato, 1966; Hiki, 1981).
A brief primer is given here; theoretical presentations of nonlinear elasticity for anisotropic crystals are available in books on the subject (Thurston, 1974; Teodosiu, 1982; Clayton, 2011). Let be the spatial position vector of a material particle that occupied a reference position at some initial time . Cartesian coordinates (xk,XK) with are used for Euclidean 3-space, and repeated indices are summed. The deformation gradient is the two-point tensor
Per the Cauchy–Born rule, primitive ; Clayton, 2011). The classical strain measure for nonlinear crystal elasticity is the Green–Lagrange strain (Wallace, 1972; Thurston, 1974; Clayton, 2011):
vectors of a crystalline material at point deform affinely with (Born & Huang, 1954Strain energy per unit initial volume is . Isentropic or isothermal conditions are implied; elastic constants are isentropic or isothermal values. Cauchy (true) stress is
Using Greek indices for Voigt notation (Thurston, 1974; Clayton, 2011), series expansion of W gives
The constant datum energy is W0, first-order for a stress-free reference state, second-order constants are , third-order are , and so on. Constants have implied symmetries:
For materials (i.e. crystal classes) of lowest (i.e. triclinic) symmetry, and have 21 and 56 independent components, respectively.
Denote by an orthogonal matrix () that belongs to the symmetry group of transformations for a given material (e.g. its Laue group or crystal class). Then yields constraints among the elastic constants of each order dictated by that symmetry group. The greater the intrinsic symmetry, the more expansive the symmetry group, and heuristically the fewer independent elastic constants of a given order. For example, if a material is isotropic, can be any rotation, limiting the number of independent to two and to three.
For crystal structures, standard conventions are used to relate coordinate axes to directions in the lattice (Brainerd et al., 1949); an example is shown in Fig. 1 for a triclinic crystal. Independent second- and third-order elastic constants (Fumi, 1951; Hearmon, 1953) for all 11 Laue groups and isotropic solids are reported for easy reference in Tables 1 and 2, following Brugger (1965), Thurston (1974) and Clayton (2011).
|
|
Also of theoretical interest are Cauchy symmetries arising if all interatomic forces are central, as from a pair potential (Love, 1927; Born & Huang, 1954). In tensor form, second-order constants then obey ; similar constraints arise at higher orders. In cubic classes of greatest symmetry, , and in Voigt notation. In a Cauchy isotropic solid, and each contain but one independent constant, and Poisson's ratio .
Symmetric second-order strain tensors differing from of (2) have been used in nonlinear elastic potentials akin to (4), with noted advantages for describing DFT (Nielsen, 1986) and shock (Clayton, 2019) data. Symmetries of elastic constant tensors of all orders are unchanged (e.g. Tables 1 and 2 remain valid) so long as the strain tensor has components referred to {XK}; any such strain transforms the same under the action of . Values of are identical for all such strain measures, but values of and constants of successively higher orders generally differ, and (3) is transformed for strain different from . Another example is linear elasticity, for which the small strain tensor is with displacement . To first order in , . Informally, when , and W is truncated at order two, (3) and (4) give Hooke's law:
The effective elastic constants of heterogeneous solids (e.g. polycrystalline aggregates) depend on the properties and orientations of constituents (Gnaupel-Herold, 2023). If an aggregate has a certain target symmetry, estimates or bounds on its effective `averaged' constants can be obtained from theoretical averaging schemes. For second-order constants, well known estimates include upper and lower bounds of Voigt and Reuss, respectively, Hill's proposition (Hill, 1952) that interpolates between the two, and `self-consistent' models (Kroner, 1958; De Wit, 1997). Such methods were extended to third-order constants for isotropic (Barsch, 1968; Lubarda, 1997) and textured (Johnson, 1985; Kube & Turner, 2016) polycrystals. Isotropic Voigt averages of fourth-order constants for aggregates of cubic crystals were derived by Krasilnikov & Vekilov (2019).
2. Nonlinear crystal elasticity to sixth order
Telyatnik (2024) recently developed numerically efficient algorithms for symbolic computations of effective elastic constants of orders two through six for polycrystalline aggregates having overall target symmetries of any transverse isotropy or full isotropy. Constituent crystallites can have any anisotropy. Effective constants are defined as arithmetic averages over the minimal set of symmetry operations generating a given target symmetry. In cases of targeted transverse or full isotropy, continuous integrals replace discrete averages. For efficiency, nested calculations for higher symmetries apply precomputed averages from lower symmetries. Gauss–Jordan elimination is used to find algebraic relationships among all elastic constants of a given order, for each target symmetry. Computations exceed capabilities of existing tools supplementing Vol. D of the International Tables for Crystallography (Authier & Zarembowitch, 2003).
Independent elastic constants and symmetry relationships for all crystal classes, transverse isotropy and full isotropy, for orders two through six, are provided in Appendix A of Telyatnik (2024) and Telyatnik (2021). Previously, such information (e.g. as in Tables 1 and 2) was available only for constants of all crystal classes and isotropy, to order four (Brendel, 1979). Averages for independent elastic constants, again up to sixth order, are given in Appendix B of Telyatnik (2024) and Telyatnik (2021). These include the target symmetries of all crystal classes, transverse isotropy and isotropy. Other appendices (Telyatnik, 2021; Telyatnik, 2024) list all components and rotation matrices. For anisotropic target symmetries, averages are defined as described in the preceding paragraph. These averages do not incorporate data on local crystal orientations (i.e. distribution functions) included in some prior definitions of effective third-order constants (Johnson, 1985; Kube & Turner, 2016). Therefore, anisotropic aggregate constants of Telyatnik (2024) cannot be expected to reproduce effective constants of textured polycrystalline metals, for example. However, these anisotropic aggregate values can serve as higher-symmetry approximations for constants of classes of lower true symmetry (Telyatnik, 2021).
For isotropic target symmetry, averaged constants (Telyatnik, 2024) are consistent with Voigt's postulate. Previous isotropic averages were usually limited to second and third orders, with the latter for constituents having cubic or hexagonal symmetry (Barsch, 1968; Lubarda, 1997; Kube & Turner, 2016). Previously, the highest known order of derived, isotropic Voigt-averaged constants was four, as reported in Appendix A of Krasilnikov & Vekilov (2019) for crystallites of general anisotropy (e.g. triclinic symmetry) and verified independently by Telyatnik (2024). Therefore, the isotropic Voigt-type averages for five independent fifth-order constants and seven independent sixth-order constants derived by Telyatnik (2024) appear to be a new and valuable contribution to nonlinear elasticity theory for solid crystals.
Finally, note that for some solids under extreme strain [e.g. rubbery polymers (Ogden, 1984) and biological tissues (Fung, 1993)] Taylor polynomials like (4) can be cumbersome so are often replaced with other functional forms (e.g. exponentials) needing fewer constants to fit data. These solids are often idealized as incompressible, for which (3) and (4) are inappropriate.
References
Authier, A. & Zarembowitch, A. (2003). International Tables for Crystallography Vol. D, Physical Properties of Crystals, edited by A. Authier, pp. 72–98. Dordrecht: Springer. Google Scholar
Barsch, G. (1968). J. Appl. Phys. 39, 3780–3793. CrossRef CAS Web of Science Google Scholar
Born, M. & Huang, K. (1954). Dynamical Theory of Crystal Lattices. London: Oxford University Press. Google Scholar
Brainerd, J. et al. (1949). Proc. Inst. Radio Eng. 37, 1378–1395. Google Scholar
Brendel, R. (1979). Acta Cryst. A35, 525–533. CrossRef CAS IUCr Journals Web of Science Google Scholar
Brugger, K. (1965). J. Appl. Phys. 36, 759–768. CrossRef Web of Science Google Scholar
Chang, Z. & Barsch, G. (1967). Phys. Rev. Lett. 19, 1381–1382. CrossRef CAS Web of Science Google Scholar
Chen, H., Zarkevich, N., Levitas, V., Johnson, D. & Zhang, X. (2020). npj Comput. Mater. 6, 115. Web of Science CrossRef Google Scholar
Clayton, J. (2011). Nonlinear Mechanics of Crystals. Dordrecht: Springer. Google Scholar
Clayton, J. (2019). Nonlinear Elastic and Inelastic Models for Shock Compression of Crystalline Solids. Cham: Springer. Google Scholar
De Wit, R. (1997). J. Appl. Cryst. 30, 510–511. CrossRef CAS Web of Science IUCr Journals Google Scholar
Fowles, R. (1967). J. Geophys. Res. 72, 5729–5742. CrossRef CAS Web of Science Google Scholar
Fumi, F. (1951). Phys. Rev. 83, 1274–1275. CrossRef CAS Web of Science Google Scholar
Fung, Y.-C. (1993). Biomechanics: Mechanical Properties of Living Tissues, 2nd ed. New York: Springer. Google Scholar
Gnäupel-Herold, T. (2023). J. Appl. Cryst. 56, 1658–1673. Web of Science CrossRef IUCr Journals Google Scholar
Graham, R. (1972). J. Acoust. Soc. Am. 51, 1576–1581. CrossRef CAS Web of Science Google Scholar
Hearmon, R. (1946). Rev. Mod. Phys. 18, 409–440. CrossRef CAS Web of Science Google Scholar
Hearmon, R. F. S. (1953). Acta Cryst. 6, 331–340. CrossRef IUCr Journals Web of Science Google Scholar
Hiki, Y. (1981). Annu. Rev. Mater. Sci. 11, 51–73. CrossRef CAS Web of Science Google Scholar
Hiki, Y. & Granato, A. (1966). Phys. Rev. 144, 411–419. CrossRef CAS Web of Science Google Scholar
Hill, R. (1952). Proc. Phys. Soc. A, 65, 349–354. CrossRef Web of Science Google Scholar
Johnson, G. (1985). J. Appl. Mech. 52, 659–663. CrossRef Web of Science Google Scholar
Krasilnikov, O. & Vekilov, Y. (2019). Phys. Rev. B, 100, 134107. Web of Science CrossRef Google Scholar
Kroner, E. (1958). Z. Phys. 151, 504–518. CAS Google Scholar
Kube, C. & Turner, J. (2016). J. Elast. 122, 157–177. Web of Science CrossRef Google Scholar
Love, A. (1927). A Treatise on the Mathematical Theory of Elasticity, 4th ed. Cambridge University Press. Google Scholar
Lubarda, V. (1997). J. Mech. Phys. Solids, 45, 471–490. CrossRef CAS Web of Science Google Scholar
Nielsen, O. (1986). Phys. Rev. B, 34, 5808–5819. CrossRef CAS Web of Science Google Scholar
Ogden, R. (1984). Non-Linear Elastic Deformations. Chichester: Ellis-Horwood. Google Scholar
Pandit, A. & Bongiorno, A. (2023). Comput. Phys. Commun. 288, 108751. Web of Science CrossRef Google Scholar
Telyatnik, R. (2021). Mendeley Data, V2, https://data.mendeley.com/datasets/mf8rbjzwmw/2. Google Scholar
Telyatnik, R. S. (2024). Acta Cryst. A80, 394–404. Web of Science CrossRef IUCr Journals Google Scholar
Teodosiu, C. (1982). Elastic Models of Crystal Defects. Berlin: Springer. Google Scholar
Thurston, R. (1974). Handbuch der Physik, edited by C. Truesdell, vol. VI, pp. 109–308. Berlin: Springer. Google Scholar
Thurston, R. & Brugger, K. (1964). Phys. Rev. 133, A1604–A1610. CrossRef Google Scholar
Thurston, R., McSkimin, H. & Andreatch, P. (1966). J. Appl. Phys. 37, 267–275. CrossRef CAS Web of Science Google Scholar
Voigt, W. (1910). Lehrbuch der Kristallphysik. Leipzig: Teubner. Google Scholar
Wallace, D. (1972). Thermodynamics of Crystals. New York: John Wiley and Sons. Google Scholar
This article is published by the International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.