Symmetry, X-Ray Diffraction (XRD), Band Structures, and Elastic Tensors
symmetry and asymmetry; reflections on X-rays; structure, both electronic and crystalline; stressing crystals
Molecules and materials often contain atoms that are identical except for rotations, translations, or reflections. Identifying these symmetry elements is useful for understanding the properties of the structure; entire classes are taught in graduate school to identify point groups and space groups and to help determine the corresponding properties (many a computational chemist has a cherished copy of Cotton’s Chemical Applications of Group Theory lying around somewhere).
Symmetry is often used to:
simplify electronic structure calculations: many matrix elements and integrals must be zero due to symmetry considerations, thus simplifying calculations;
inform spectral interpretation: IR-active vibrations transform as one of the Cartesian coordinates (x, y, or z) under the dipole selection rule;
determine rotational spectral line spacing and relative intensities; and
correct overcounting of entropy from the ideal-gas, rigid-rotor harmonic-oscillator (RRHO) approximation when a molecule has proper rotational symmetry elements (Cₙ axes, etc.).
Previously, the Rowan platform used the pymsymm package for computing point-group symmetry corrections to RRHO. Since it was written in C and didn’t support periodic systems, it was going to be difficult to extend with various features we were interested in. We therefore rewrote our point-group code in Python (achieving nearly the same speed as the original C code, without the risk of segfaults). For the 230 space groups of periodic systems, we added spglib. Now, molecules and crystals have their symmetry attached, and it can be viewed in the basic calculation workflow.
Molecules are often imported and drawn with imperfect symmetry, and now they can be symmetrized in the structure builder, allowing users to quickly fix nearly symmetric geometries.
Symmetry labels have also been added to orbitals in the electronic properties workflow for easier identification and use in chemistry classes.
XRD
For crystal structures, the X-ray diffraction pattern of the material is a standard fingerprint, as the height and positioning of the peaks indicate which of the 230 space groups it could belong to. When X-rays hit a crystal, the X-rays are absorbed and then re-emitted by the electrons surrounding each atom (Thomson scattering). If X-rays are shone on the crystal from various angles, the constructive and destructive interference patterns generated in the scattering can be used to infer the symmetry of the crystal (Bragg’s law: nλ = 2d sin θ).
The X-ray diffraction pattern can be directly predicted from the material’s structure and is useful for identifying materials, especially polymorphs. Real samples will have defects, mixtures of phases, preferred orientation, temperature effects, and broadened peaks, making the observed pattern slightly different, but they often align well with the predicted XRD pattern.
Carbon has numerous allotropes (single-element polymorphs) and the major ones have strikingly different XRD patterns. Changing from cubic diamond to lonsdaleite (hexagonal diamond) markedly lowers the crystal symmetry from the high-symmetry cubic space group Fd-3m to the hexagonal P6₃/mmc, which breaks equivalences among lattice planes and results in many more distinct XRD peaks.
Band Structures and Density of States
While XRD probes the positions of the nuclei and the resulting symmetry, we can also compute the electronic structure of a crystal. The periodic analogs of molecular orbitals are the band structure and density of states, which govern conductivity, optical absorption, reactivity, and more. These are shaped by the crystal’s lattice, bonding, and the way the electronic wavefunction repeats through space.
The band gap is the separation between the valence band maximum (VBM) and conduction band minimum (CBM); the periodic counterpart of the HOMO–LUMO gap. Its size and location determine the material’s electronic and optical behavior. A direct gap occurs when the VBM and CBM sit at the same point in the band structure, and an indirect gap occurs when they sit at different points. Direct-gap materials conserve momentum during electron–hole recombination and emit light efficiently, making them ideal for optoelectronics; indirect-gap materials are poor emitters, so charge carriers travel farther before recombining, which is useful in semiconductors.
Carbon nanotubes illustrate how geometry alone can flip this behavior. A nanotube is a sheet of graphene rolled into a cylinder, and the rolling direction changes which electronic states are available around the circumference, even though the atoms and bonding are identical. Two types of nanotubes are armchair and zigzag nanotubes.
Zigzag nanotubes are rolled along one of graphene primitive lattice vector direction (n,0), while armchair nanotubes are rolled along the sum of the two primitive lattice vectors (n,n). This zone-folding changes how graphene’s electronic bands align after rolling: all armchair nanotubes are metallic (zero or very small band gap), while zigzag nanotubes are metallic when n is a multiple of 3 and semiconducting otherwise.
Subtle geometric differences influence the electronic properties not only in carbon nanotubes but across a wide range of materials. For instance, VO₂ undergoes a dramatic metal–insulator transition between its low-temperature monoclinic (insulating) and high-temperature rutile (metallic) structures, accompanied by significant changes in its electronic structure. More broadly, entire fields such as twistronics, straintronics, and the study of topological semimetals exploit how small changes in geometry can dramatically alter band structures, density of states, and conductivity.
Elastic Tensor
While XRD characterizes a material’s atomic arrangement and band structures characterize a material’s electronic structure, the elastic tensor characterizes a material’s response to mechanical stress. The full tensor is a 6×6 matrix, and each entry Ci,j relates an applied stress to the resulting strain (deformation) in the elastic regime, together describing the stiffness and compressibility of the material.
A few scalar properties distill the full tensor into intuitive quantities: the bulk modulus (K) measures resistance to compression, the shear modulus (G) measures resistance to shearing, and their ratio, Pugh’s ratio K/G, is a reliable predictor of ductility. Materials with K/G > 1.75 tend to deform plastically before fracturing where those with lower values tend to cleave brittlely.
Silicon is a textbook brittle material: high stiffness and a Pugh’s ratio of 1.40, well below the ductility threshold. By contrast, gold has a Pugh’s ratio above 1.75, reflecting the ability of its metallic bonds to redistribute stress before fracture.
The elastic tensor also encodes directional anisotropy; a cubic crystal like silicon is nearly isotropic (C₁₁ ≈ C₂₂ ≈ C₃₃ and C₄₄ ≈ C₅₅ ≈ C₆₆), while lower-symmetry materials can be orders of magnitude stiffer along one axis than another. This directional information, often invisible to scalar moduli alone, matters for designing layered coatings, piezoelectrics, and materials under directional loads.
We are excited to see how these new workflows will be used in materials discovery pipelines, and have several more materials workflows in progress. If you work in the materials space, we’d be happy to chat about what you would like to see next on our platform, or how you can integrate our workflows to your existing discovery pipelines via rowan-python.
 | A guest post by
|
Rowan is getting more cooler by the day!