Voigt notation

In fields like physics and materials science, describing how a material deforms under stress is fundamental. This relationship, however, is not a simple one-to-one correspondence; it is governed by tensors, complex mathematical objects with multiple directional components. Representing the full stress-strain relationship involves a fourth-rank stiffness tensor with a staggering 81 components, creating a significant challenge for calculation and analysis. This article addresses this complexity by introducing Voigt notation, a brilliantly practical scheme for simplifying these tensor equations. Across the following chapters, we will explore the core principles of this powerful method and its wide-ranging applications. The first chapter, "Principles and Mechanisms," delves into how Voigt notation works, the physical trade-offs required to use it, and how it elegantly captures a material's essential properties. Following that, "Applications and Interdisciplinary Connections" will demonstrate how this notation is an indispensable tool in mechanics, optics, and crystal physics, bridging abstract theory with tangible engineering design.

In many areas of science and engineering, it is often necessary to track multiple interacting physical quantities simultaneously, similar to an accounting process. For example, when an external force is applied to an object like a block of rubber, it does not simply compress in the direction of the force; it also deforms in other directions, such as bulging out to the sides. This relationship between a push (​​stress​​) and the resulting deformation (​​strain​​) is not a simple one-to-one affair. Both stress and strain are ​​tensors​​—objects that describe properties that have a magnitude and multiple directions.

Let's write down the full relationship, called Hooke's Law for materials: $\sigma_{ij} = C_{ijkl} \varepsilon_{kl}$. The stress tensor $\sigma_{ij}$ and strain tensor $\varepsilon_{kl}$ each have 9 components in 3D space, describing forces and deformations on every face in every direction. The object connecting them, $C_{ijkl}$, is the fourth-rank ​​stiffness tensor​​. In principle, it has $3 \times 3 \times 3 \times 3 = 81$ components! Writing this out is a chore, and programming it into a computer is a recipe for bugs and headaches. This is a classic bookkeeping problem. How can we simplify our accounts without losing the physics?

Around the turn of the 20th century, Woldemar Voigt had a wonderfully practical idea. He noticed that the stress and strain tensors are symmetric (e.g., the shear stress on the top face in the x-direction is the same as the shear stress on the side face in the y-direction, $\sigma_{xy} = \sigma_{yx}$). So, instead of 9 components for each, there are really only 6 independent ones. Why not, he thought, just list them in a column? We can map the two indices of the tensor to a single index, like this:

$11 \to 1$ $22 \to 2$ $33 \to 3$ $23 \text{ and } 32 \to 4$ $13 \text{ and } 31 \to 5$ $12 \text{ and } 21 \to 6$

Suddenly, our $3 \times 3$ symmetric tensors $\boldsymbol{\sigma}$ and $\boldsymbol{\varepsilon}$ become $6 \times 1$ vectors (which we can call $\mathbf{T}$ and $\mathbf{S}$). What was a monstrous rank-4 tensor $C_{ijkl}$ with 81 components now becomes a much friendlier $6 \times 6$ matrix, $\mathbf{C}$. The complicated tensor equation becomes a familiar matrix equation: $\mathbf{T} = \mathbf{C} \mathbf{S}$. This is the essence of ​​Voigt notation​​. It's a re-indexing scheme, a piece of clever bookkeeping.

But wait. Physics is a jealous guardian of its laws. We can't just change our notation for free. There must be a catch, a price to pay for this convenience. The price is revealed when we consider the physics of ​​work and energy​​. The work done per unit volume to deform a material is given by the elegant tensor contraction $\dot{W} = \sigma_{ij} \dot{\varepsilon}_{ij}$ (where the dot means "rate of change"). Let's write it out:

$\dot{W} = \sigma_{11}\dot{\varepsilon}_{11} + \sigma_{22}\dot{\varepsilon}_{22} + \sigma_{33}\dot{\varepsilon}_{33} + \sigma_{12}\dot{\varepsilon}_{12} + \sigma_{21}\dot{\varepsilon}_{21} + \dots$

Because of symmetry ($\sigma_{12} = \sigma_{21}$), this becomes:

$\dot{W} = \sigma_{11}\dot{\varepsilon}_{11} + \sigma_{22}\dot{\varepsilon}_{22} + \sigma_{33}\dot{\varepsilon}_{33} + 2\sigma_{12}\dot{\varepsilon}_{12} + 2\sigma_{13}\dot{\varepsilon}_{13} + 2\sigma_{23}\dot{\varepsilon}_{23}$

Look at that! A factor of 2 has appeared in front of the shear terms. If we want our new Voigt notation to preserve this physical law—if we want the work to be a simple dot product $\dot{W} = \mathbf{T}^{\mathsf{T}}\dot{\mathbf{S}}$—we have to make a choice. We define our Voigt stress vector simply as $T_1 = \sigma_{11}, \dots, T_6 = \sigma_{12}$. But to make the math work out, we must define the Voigt strain vector by absorbing those factors of 2. We define the ​​engineering shear strain​​:

$S_1 = \varepsilon_{11}$, $S_2 = \varepsilon_{22}$, $S_3 = \varepsilon_{33}$ $S_4 = 2\varepsilon_{23}$, $S_5 = 2\varepsilon_{13}$, $S_6 = 2\varepsilon_{12}$

This is the physicist's bargain. We get a simple matrix form, but we must remember that the shear components of our strain "vector" $\mathbf{S}$ are twice the tensor shear strains. This isn't an arbitrary choice; it's a requirement to ensure our notation is ​​work-conjugate​​, a fancy way of saying it correctly represents the flow of energy.

With this bargain struck, we can now look at the grand prize: the $6 \times 6$ stiffness matrix $\mathbf{C}$. It turns out that this matrix isn't just an arbitrary collection of 36 numbers. A deep principle of thermodynamics requires that for an elastic material, the work you do on it is stored as potential energy (strain energy, $W$). This means the order in which you apply strains doesn't matter for the final energy stored. A mathematical consequence of this physical fact is that the stiffness tensor must have a "major symmetry," $C_{ijkl} = C_{klij}$. In our Voigt notation, this translates to a beautifully simple condition: the matrix $\mathbf{C}$ must be symmetric!

$\mathbf{C} = \mathbf{C}^{\mathsf{T}}$

A candidate matrix that is not symmetric cannot represent a real elastic material, because it would violate the principle of energy conservation. Because of this symmetry, the 36 components are reduced to just $\frac{6 \times (6+1)}{2} = 21$ independent constants for the most general anisotropic material (a triclinic crystal). These 21 numbers are the material's elastic "soul."

Twenty-one constants are better than 81, but it's still a lot. Fortunately, most materials have internal symmetries. The atoms in a crystal are arranged in a repeating, orderly lattice. This internal order imposes strict rules on the stiffness matrix. If a crystal looks the same after you rotate it, its elastic response must also be the same. This powerful idea makes many of the 21 constants zero and forces others to be equal.

• For a ​​cubic crystal​​ (like salt or iron), the matrix simplifies dramatically. There are only 3 independent constants: $C_{11}$, $C_{12}$, and $C_{44}$.
  $$\mathbf{C}_{\text{cubic}} = \begin{pmatrix} C_{11} C_{12} C_{12} 0 0 0 \\ C_{12} C_{11} C_{12} 0 0 0 \\ C_{12} C_{12} C_{11} 0 0 0 \\ 0 0 0 C_{44} 0 0 \\ 0 0 0 0 C_{44} 0 \\ 0 0 0 0 0 C_{44} \end{pmatrix}$$
• For a ​​transversely isotropic​​ material (like a hexagonal crystal or a fiber-reinforced composite), which is symmetric around one axis, we have 5 independent constants.

Each crystal symmetry class has its own characteristic stiffness matrix "fingerprint." We can read the underlying symmetry of a material by looking at the structure of its Voigt matrix.

The numbers in the matrix are not just abstract entries; they tell a physical story.

First, they tell us about ​​anisotropy​​—the property of having different characteristics in different directions. The constants in the matrix are only "constant" if you measure them along the crystal's natural axes. If you rotate your coordinate system, the components of the stiffness tensor transform, and the numbers in your Voigt matrix change! For a cubic crystal rotated by 45 degrees around the z-axis, the new $C'_{11}$ component becomes a mixture of the old constants: $C'_{11} = \frac{1}{2}(C_{11} + C_{12}) + C_{44}$. The resistance to stretching in this new direction is different. We can even quantify the degree of anisotropy with numbers like the Zener ratio, $A = 2C_{44} / (C_{11}-C_{12})$, which is 1 for a perfectly isotropic material and different from 1 for an anisotropic cubic crystal.

Second, they tell us about ​​stability​​. A physical material must be stable. If you poke it, it should resist and store energy; it shouldn't spontaneously explode or collapse. This means the strain energy $W = \frac{1}{2} \mathbf{S}^{\mathsf{T}}\mathbf{C}\mathbf{S}$ must always be positive for any deformation. Mathematically, this means the matrix $\mathbf{C}$ must be ​​positive definite​​. This imposes a set of inequalities that the elastic constants must obey, for instance, that all the diagonal elements must be positive ($C_{11} > 0, C_{44} > 0, \dots$) and that all eigenvalues of the matrix must be positive. These are not mathematical curiosities; they are fundamental conditions for a material's existence.

The true beauty of Voigt notation is that it provides a universal language that extends far beyond simple elasticity. Consider ​​piezoelectricity​​, the remarkable property of some crystals (like quartz) to generate a voltage when squeezed. Here, the mechanical world of stress and strain is coupled to the electrical world of electric fields and charge displacement.

Using Voigt notation, we can write the coupled constitutive equations with breathtaking clarity and elegance:

$\mathbf{S} = \mathbf{s}^{E} \mathbf{T} + \mathbf{d}^{\mathsf{t}} \mathbf{E}$ $\mathbf{D} = \mathbf{d} \mathbf{T} + \boldsymbol{\epsilon}^{T} \mathbf{E}$

Here, $\mathbf{S}$ (strain) and $\mathbf{D}$ (electric displacement) are the "outputs," while $\mathbf{T}$ (stress) and $\mathbf{E}$ (electric field) are the "inputs." The matrices tell the whole story: $\mathbf{s}^E$ is the material's compliance (the inverse of stiffness) measured at a constant electric field, $\boldsymbol{\epsilon}^T$ is its dielectric permittivity measured at constant stress, and the matrix $\mathbf{d}$ is the piezoelectric coefficient that orchestrates the dance between the two worlds. The superscripts $E$ and $T$ are part of this precise language, telling us the conditions of the measurement.

This formalism allows us to uncover profound physics. For instance, we can rearrange these equations to find the effective stiffness of a piezoelectric material under different electrical conditions. The stiffness at constant electric displacement, $\mathbf{c}^D$, turns out to be different from the stiffness at constant electric field, $\mathbf{c}^E$. The relationship is $\mathbf{c}^{D} = \mathbf{c}^{E} + \mathbf{e}^{\mathsf{T}} (\boldsymbol{\epsilon}^{S})^{-1} \mathbf{e}$, where $\mathbf{e}$ and $\boldsymbol{\epsilon}^S$ are other material property matrices. The term $\mathbf{e}^{\mathsf{T}} (\boldsymbol{\epsilon}^{S})^{-1} \mathbf{e}$ represents "piezoelectric stiffening." The material actually becomes mechanically stiffer just by changing the electrical boundary conditions! This is a real, measurable effect, and it falls naturally out of the clean, powerful bookkeeping of Voigt notation.

From a messy accounting problem to a universal language describing the interplay of forces and fields, Voigt notation is a testament to how the right mathematical tool can not only simplify our work but also deepen our understanding of the inherent beauty and unity of the physical world.

Now that we have acquainted ourselves with the rules and structure of Voigt notation, you might be tempted to think of it as just a clever bit of bookkeeping, a convenient shorthand to tidy up our equations. And it is certainly that! But to leave it there would be like describing a grand cathedral as merely a pile of stones. The true power and beauty of a good notation lie not in what it abbreviates, but in what it reveals. Voigt notation transforms the unwieldy mathematics of tensors into a practical and intuitive language, allowing us to converse with materials, understand their personalities, and even predict their behavior across a breathtaking range of disciplines. It is a bridge from abstract principles to tangible reality.

Let's begin on the home turf of Voigt notation: the world of mechanics and materials science. Imagine you are holding a piece of rubber. You pull on it, and it stretches in the direction you pull. The relationship is simple. Now, imagine you are holding a crystal. You pull on it along one of its axes. It stretches, of course, but it might also shrink or shear in other directions! This is the essence of anisotropy—the properties of the material depend on the direction you are looking.

Generalized Hooke's law, $\sigma_{ij} = C_{ijkl} \epsilon_{kl}$, is the full, formal description of this behavior, but it's a beast to work with. Enter Voigt notation. The relationship becomes a crisp matrix equation, $\mathbf{T} = \mathbf{C} \mathbf{S}$, or its inverse, $\mathbf{S} = \mathbf{s} \mathbf{T}$. The stiffness matrix $\mathbf{C}$ or compliance matrix $\mathbf{s}$ is no longer just a collection of numbers; it becomes a fingerprint of the crystal's mechanical "personality".

For a cubic crystal, the compliance matrix has a beautifully simple form with only three independent constants: $s_{11}$, $s_{12}$, and $s_{44}$. Suppose we apply a simple uniaxial stress $T_0$ along one axis, say the $z$-axis, so that the only non-zero stress component is $T_3 = T_0$. The matrix equation immediately tells us the full story of the strain:

The result is $S_3 = s_{11} T_0$, which is the direct strain, and $S_1 = S_2 = s_{12} T_0$, which is the transverse strain (the Poisson effect). The off-diagonal term $s_{12}$ is the link, the crosstalk in the conversation between stress and strain. But what if we are interested in a direction that is not along the primary axes, like the main diagonal of the cube? The formalism handles this with ease. By transforming the strain tensor back to its full form, we can find the strain along any arbitrary direction, revealing how the crystal's internal structure dictates its response to external forces.

This matrix representation is not just descriptive; it is predictive. If an engineer has the stiffness matrix for an orthorhombic crystal, they can calculate crucial design parameters like the directional Young's modulus for any orientation within the material. It even allows us to ask deep questions about the nature of materials themselves. For instance, a cubic crystal is anisotropic, while glass is isotropic. What is the bridge between them? We can define a quantity called the Zener anisotropy factor, $A = \frac{2C_{44}}{C_{11} - C_{12}}$. If we demand that our cubic crystal behaves isotropically, its three elastic constants must conspire in such a way that this factor becomes exactly 1. Voigt notation thus provides a quantitative measure of how "crystal-like" a material is.

The framework is also wonderfully extensible. What happens when we heat a crystal? It wants to expand. This thermal strain can be represented by a simple vector, $\boldsymbol{\alpha} \Delta T$. The total strain vector $\mathbf{S}$ is then just the sum of the mechanical and thermal parts: $\mathbf{S} = \mathbf{s} \mathbf{T} + \boldsymbol{\alpha} \Delta T$. This unified equation lets us solve fascinating coupled-physics problems. For example, if we want to heat a tetragonal crystal but prevent it from getting any longer along its main axis, how much stress must we apply? The Voigt equation gives a direct and elegant answer: we need a compressive stress that precisely counteracts the thermal expansion, a value determined by the thermal expansion coefficient and the compliance constant along that axis.

The power of Voigt notation extends far beyond the realm of pushes and pulls into the subtle and beautiful world of light. The way light propagates through a crystal is determined by the material's dielectric impermeability tensor, $\boldsymbol{\eta}$ (the inverse of the dielectric tensor), which defines an "index ellipsoid." The shape and orientation of this ellipsoid tell an incoming light wave what refractive indices it will experience.

In many advanced optical devices, we want to control this ellipsoid on demand. This is the magic of the electro-optic or Pockels effect. By applying an electric field $\mathbf{E}$, we can induce a small change in the impermeability tensor: $\Delta \eta_{i} = \sum_{k} r_{ik} E_k$, where $r_{ik}$ is the electro-optic tensor in Voigt form.

Consider a Pockels cell, a cornerstone of modern optics used for high-speed switching and modulation. A crystal like KDP (potassium dihydrogen phosphate) is placed in an electric field. Let's say we apply the field $E_z$ along the optic axis. The Voigt notation tells us which components of the impermeability tensor will change. For KDP's $\bar{4}2m$ symmetry, the relevant coefficient is $r_{63}$. This means the applied field induces a $\Delta \eta_6$ component, which corresponds to an $xy$ term in the index ellipsoid equation. A cross-section of the ellipsoid that was a perfect circle now becomes an ellipse, rotated by 45 degrees. This creates two new principal axes with different refractive indices. A light wave traveling along the z-axis will be split into two polarization components that travel at different speeds, accumulating a phase difference. By applying a specific voltage, the "half-wave voltage" $V_{\pi}$, we can make this phase difference exactly $\pi$, effectively rotating the polarization of the light. The expression for this critical voltage falls directly out of the Voigt analysis. The same principles apply to different geometries, such as a transverse Pockels cell where the field is applied perpendicular to the light's path.

The story gets even more exciting in the realm of nonlinear optics, where intense laser light itself can alter the optical properties of a material. Phenomena like sum-frequency generation (SFG), where two photons of different frequencies combine to create a new photon at the sum frequency, are governed by the second-order nonlinear susceptibility tensor, $d_{ijk}$. Again, Voigt notation comes to our rescue, contracting this to a $3 \times 6$ matrix $d_{il}$. This matrix is nothing less than a recipe book for the experimentalist. Suppose the table of coefficients for a KDP crystal shows a non-zero $d_{36}$ component. What does this mean? The indices tell you everything! The first index, $3$, tells you the generated nonlinear polarization (and thus the output light) will be along the $z$-axis. The second index, $6$, corresponds to the $xy$ product of input fields. To exploit this coefficient, you must therefore send in one input beam polarized along the $x$-axis and another polarized along the $y$-axis. The crystal will dutifully mix them and produce new light polarized along the $z$-axis. The notation is a direct map from a tensor table to a laboratory setup. For quantitative design, one can calculate the "effective" nonlinear coefficient, $d_{\text{eff}}$, for any interaction geometry, which is a scalar value that determines the efficiency of the process. This calculation, a projection of the $d$ tensor onto the field polarizations, is made immeasurably clearer using the Voigt matrix.

We have seen Voigt notation describe elasticity and optics, but its deepest power comes from its intimate relationship with the most fundamental concept in crystal physics: symmetry. Neumann's Principle is the grand decree: any physical property of a crystal must be invariant under all the symmetry operations of the crystal's point group. This means the tensor describing the property must have at least the symmetry of the crystal itself.

This principle acts as a powerful "veto." For example, any crystal that possesses a center of inversion (is centrosymmetric) cannot be piezoelectric. Piezoelectricity, the coupling of polarization and strain, is described by a third-rank tensor like $d_{ijk}$ or $e_{ijk}$. An inversion operation flips the sign of all coordinates, and for a third-rank tensor, this means the tensor component itself flips sign. But invariance requires the component to remain unchanged. The only way for a number to be equal to its negative is if it is zero. Thus, all piezoelectric coefficients in a centrosymmetric crystal must be zero.

What happens when a crystal is not centrosymmetric? Symmetry still imposes powerful constraints. Consider a crystal with $C_{6v}$ symmetry. It has a six-fold rotation axis and mirror planes. Instead of 18 independent piezoelectric coefficients in the $3 \times 6$ Voigt matrix $e_{iJ}$, we can apply the symmetry operations one by one to see which components are forced to be zero and which must be equal to others. The rigorous application of these symmetry constraints reveals that the entire complex tensor is determined by only three independent, non-zero values! This is a colossal simplification, a gift from nature's love of symmetry, made manifest through our notational system.

Perhaps the most profound application lies in understanding phase transitions. Imagine a material in a high-temperature phase with a high-symmetry, centrosymmetric structure like $D_{2h}$. As it cools, it might undergo a phase transition, where the atoms shift slightly into a new arrangement with a lower, non-centrosymmetric symmetry, like $C_s$. In the high-symmetry phase, the material cannot be piezoelectric; the $d_{i\lambda}$ matrix is entirely zero, forbidden by symmetry. But in the low-symmetry phase, the inversion center is gone! The veto is lifted. Suddenly, components of the piezoelectric tensor are allowed to be non-zero. The once-empty $d_{i\lambda}$ matrix springs to life with a pattern of 10 independent non-zero coefficients, a pattern dictated entirely by the new, lower symmetry. This is not just a mathematical curiosity; it is the birth of a physical property, a direct and observable consequence of a change in symmetry, all described perfectly by the population and structure of a Voigt matrix.

In the end, Voigt notation is far more than a compact way to write things down. It is a language that unifies disparate fields, from mechanics and thermodynamics to optics and condensed matter theory. It provides a computational engine that connects abstract symmetry principles to concrete engineering designs. And in doing so, it reveals a profound and elegant unity in the way materials respond to the forces and fields of the universe.