The team, which included Assistant Professor Tatsuya Mori of the University of Tsukuba, Assistant Professor Hideyuki Mizuno of the University of Tokyo, Specially Appointed Associate Professor Yasuhiro Fujii of the University of Osaka, and Professor Akitoshi Koreeda of Ritsumeikan University, developed a continuum model designed to describe how terahertz electromagnetic waves interact with the vibrational dynamics of glass. The crucial innovation lies in what the model chooses to include. Rather than treating glass as a perfectly uniform, featureless medium, the framework incorporates elastic heterogeneity, the well-established fact that the stiffness of a glass varies from place to place at the nanometer scale. On top of that mechanical disorder, the model adds microscopic charge fluctuations occurring at atomic and molecular scales, which are what actually allow an electromagnetic wave to couple to the vibrations of the material.

This combination is physically motivated. When a terahertz wave travels through a solid, its oscillating electric field pushes on the charges inside the material. In a perfectly ordered crystal, those charges respond in sharply defined, symmetry-selected ways. In a glass, where the atomic arrangement is disordered and the mechanical properties fluctuate in space, the vibrational modes are a complicated mixture, and the way they couple to light is correspondingly messy. The new model captures this messiness in a controlled continuum description, translating the elastic heterogeneity and charge fluctuations of the microscopic world into quantities that can be computed and compared directly with optical measurements.

To test the model, the researchers applied it to glycerol glass, a representative and well-studied molecular glass. The comparison was demanding: the model had to reproduce both the real and imaginary parts of the complex dielectric function, the quantity that determines how a material absorbs, stores, and transmits electromagnetic energy. Across the frequency range from 0.3 to 2.5 terahertz, the calculated dielectric response matched the experimentally measured values closely. That agreement across both components of a complex response function, rather than just a single absorption curve, is a strong indication that the underlying physics has been captured rather than merely fitted.

Perhaps the most striking success of the model concerns a long-standing mystery in glass physics known as the boson peak. Glasses exhibit an anomalous excess of vibrational states at a characteristic frequency, and the dynamics of the material change qualitatively around this frequency. The new calculations captured a clear transition in the terahertz dielectric response: below the boson-peak frequency, the response resembles a resonance-like behavior, while above it the response shifts into a broad relaxation-like form. Being able to reproduce this crossover quantitatively links the optical behavior of glass directly to the crossover in its vibrational dynamics, providing the kind of microscopic-to-macroscopic connection that has been lacking.

The analysis also revealed which part of the vibrational dynamics dominates the terahertz response. Around the boson peak, the researchers found that the transverse contribution, associated with shear-like atomic motions perpendicular to the wave direction, predominantly shapes both the real and imaginary parts of the dielectric function. The longitudinal contribution, tied to compressive motions along the wave direction, is comparatively small. This finding points to transverse shear dynamics as the pivotal player in how glass interacts with terahertz light, a conclusion with real physical weight because it identifies the specific class of disordered vibrations that engineers would need to control in order to tune terahertz absorption.

Another demanding test involved the infrared light-vibration coupling coefficient, a quantity that measures how effectively vibrational motion converts into optical absorption. Near the boson peak, experiments show that this coefficient depends nearly linearly on frequency. The model replicated this nearly linear frequency dependence, adding an independent piece of evidence that the framework is describing the correct physics rather than reproducing one dataset by coincidence. Matching the dielectric function, the resonance-to-relaxation crossover, and the coupling coefficient within a single model is a notable consistency achievement for a material as structurally complex as glass.

The practical implications extend into one of the most anticipated technology areas of the coming decade. Terahertz frequencies promise wireless links with enormous bandwidth, advanced imaging, and new photonic devices, but the technology has been held back in part by the difficulty of finding materials that transmit terahertz waves with low loss. Because glass is a foundational material for optics, packaging, and fiber technology, understanding precisely why it becomes lossy at terahertz frequencies is a prerequisite for designing around the problem. The new framework provides a way to quantitatively link terahertz absorption to charge fluctuations using material-specific mechanical properties, which means that measured or simulated mechanical characteristics of a glass can be translated into predictions of its terahertz dielectric behavior.

That predictive capability points toward a concrete engineering goal: the design and evaluation of glass materials with low permittivity and low dielectric loss, characteristics that are increasingly important for terahertz communications and photonic technologies. Instead of screening candidate glasses through laborious terahertz measurements alone, researchers could use the continuum model to estimate how changes in composition, and therefore in elastic heterogeneity and charge fluctuations, would alter terahertz transmission. The work, published in Physical Review B under the title Continuum model for the terahertz dielectric response of glasses, was supported by JSPS KAKENHI grants and by the Asahi Glass Foundation, reflecting both national and industry interest in the problem.

For a phenomenon that has been observed experimentally for some time but never quantitatively explained, the arrival of a working model marks a genuine turning point. It reframes the terahertz opacity of glass not as an arbitrary material limitation but as a direct consequence of elastic disorder and charge fluctuations acting together at the nanoscale. As terahertz technology moves from laboratory demonstrations toward real-world deployment in communications and sensing, the ability to predict and ultimately engineer the dielectric response of disordered solids could shape which materials make it into the devices of the future. What began as an effort to explain why glass dims at terahertz frequencies may end up helping to brighten the entire terahertz technology landscape.

Subject of Research: A continuum model linking the terahertz dielectric response of glass to elastic heterogeneity and microscopic charge fluctuations

Article Title: New model explains why glass becomes less transparent to terahertz light

Article References: New model explains why glass becomes less transparent to terahertz light. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: glass, terahertz, dielectric function, boson peak, elastic heterogeneity, charge fluctuations, vibrational dynamics, molecular glass, glycerol glass, terahertz communications, photonic technologies, continuum model

At ambient conditions, zirconia exists in a monoclinic crystal structure belonging to the space group P21/c, while at elevated temperatures approaching 2650 kelvin it transforms into a highly symmetric cubic fluorite structure with space group Fm3̄m. These are not merely cosmetic differences. In the cubic phase, every zirconium atom is surrounded by eight equivalent oxygen atoms at identical bond lengths of 2.21 angstroms and ideal tetrahedral angles of 109.47 degrees, producing a perfectly isotropic material whose optical properties are the same in every direction. In the monoclinic phase, that coordination collapses to a distorted sevenfold environment, with seven distinct zirconium-oxygen bond lengths ranging from 2.02 to 2.17 angstroms and bond angles scattered wildly between 52.95 and 107.07 degrees. The new research demonstrates that this symmetry breaking is precisely what transforms zirconia from an optically uniform solid into a directionally sensitive, birefringent crystal.

The team carried out their investigation using density functional theory with the plane-wave pseudopotential method implemented in the CASTEP code, treating exchange and correlation through the generalized gradient approximation of Perdew, Burke, and Ernzerhof. Wavefunctions were expanded to a kinetic energy cutoff of 550 electronvolts, and Brillouin-zone integration employed a dense 12 by 12 by 12 Monkhorst-Pack grid with total energy convergence tighter than 10 microelectronvolts per atom. Because standard semi-local functionals systematically underestimate band gaps, the researchers applied an empirical scissors operator of 1.5 electronvolts to shift conduction bands upward into alignment with experimental benchmarks, a computationally efficient strategy that preserves the quantum mechanical description of symmetry breaking while avoiding the prohibitive cost of hybrid functionals or GW many-body calculations across the dense k-point meshes needed for high-resolution directional optics.

The structural results anchor the study in experimental reality. The optimized cubic lattice constant came out at 5.113 angstroms, within about 0.45 percent of the experimental value of 5.09 angstroms, and the monoclinic parameters of a equals 5.19, b equals 5.24, and c equals 5.38 angstroms with a beta angle of 99.23 degrees matched measured data closely. Energy-volume curves fitted with the Murnaghan equation of state confirmed that the monoclinic phase, with a minimum ground-state energy of minus 8715.4763 electronvolts, is thermodynamically more stable than the cubic phase at minus 8642.1649 electronvolts, consistent with decades of experimental observation. Interestingly, the cubic lattice is mechanically stiffer, with a bulk modulus of 210 gigapascals against 185 gigapascals for the monoclinic form, reflecting the stronger and more uniform bonding network of the high-symmetry structure.

The electronic structure calculations reveal how symmetry reshapes the fundamental band gap itself. The cubic phase exhibits a direct gap of 4.83 electronvolts, with both the valence band maximum and conduction band minimum sitting at the Gamma point of the Brillouin zone. The monoclinic polymorph, by contrast, shows an indirect gap of 5.17 electronvolts, with the conduction band minimum displaced to the B point. The researchers trace this shift to crystal-field effects: in the cubic lattice, the Oh-like crystal field preserves orbital degeneracy so that the highest occupied oxygen 2p states and the lowest unoccupied zirconium 4d states converge at the same momentum vector, enabling direct optical transitions. In the monoclinic lattice, the low-symmetry C2h crystal field splits the zirconium 4d t2g and eg manifolds, and the highly non-uniform bond lengths force the p-d orbital overlaps to become directional and asymmetric. By Bloch’s theorem, this real-space asymmetry displaces the conduction band minimum in reciprocal space, making the fundamental excitation indirect and dependent on phonon coupling to conserve crystal momentum.

The density of states analysis adds a clear chemical picture. The upper valence band from minus 6 electronvolts to the Fermi level is dominated by oxygen 2p states hybridized with zirconium 4d orbitals, while the conduction band edge is governed by unoccupied zirconium 4d states, confirming that the fundamental optical excitations are oxygen 2p to zirconium 4d charge-transfer transitions. The flatter, less dispersive bands of the monoclinic phase signal a larger effective mass and tighter spatial localization of charge carriers, a direct electronic fingerprint of the reduced orbital overlap in the distorted sevenfold coordination environment.

The optical analysis is where the study delivers its most striking quantitative insight. In the cubic phase, the imaginary part of the dielectric function shows a single sharp resonance at 9.52 electronvolts with an intensity of 10.85, the signature of degenerate electronic transitions responding identically in all directions. In the monoclinic phase, that degeneracy is lifted: the main absorption peak splits into three distinct responses along the principal crystallographic axes, peaking at 10.02 electronvolts with intensity 9.15 along the [001] direction, 9.82 electronvolts with intensity 8.32 along [010], and 9.95 electronvolts with intensity 7.81 along [100], which also carries a distinct shoulder at 7.90 electronvolts. Even the absorption edge itself splits directionally, beginning at 4.00 electronvolts along [010], 4.10 along [100], and 4.35 along [001]. The static dielectric constant tells the same story, falling from an isotropic 4.91 in the cubic phase to direction-dependent values of 4.63, 4.61, and 4.32 in the monoclinic phase, yielding a static birefringence of 0.31, roughly 7.2 percent, between the [100] and [001] axes.

Beyond the visible and ultraviolet regime, the calculations resolve a robust collective electronic resonance at 37 to 38 electronvolts in both polymorphs. Unlike ordinary valence plasmons in the 15 to 25 electronvolt range, this high-energy feature originates from deep semi-core transitions of zirconium 4p electrons into unoccupied 4d states, corresponding to the atomic zirconium N2,3 edge. The sharpness of this peak in the cubic phase, contrasted with its broader, attenuated profile in the monoclinic lattice, mirrors the structural homogeneity of the high-symmetry crystal, and the calculated plasma frequency aligns closely with experimental reflection electron energy loss and transmission electron energy loss spectroscopy measurements in the 36 to 40 electronvolt window, providing an independent validation of the computational framework.

The practical implications are considerable. Because both phases absorb almost nothing in the visible spectrum while rising sharply in the deep ultraviolet between 4.00 and 4.35 electronvolts, zirconia is a natural fit for solar-blind ultraviolet photodetectors that ignore visible solar noise, and for deep-ultraviolet transparent conducting oxides. The pronounced birefringence of the monoclinic phase opens the door to polarization filters, wave plates, phase-matching layers in integrated photonic circuits, and directional anti-reflective coatings whose performance can be tuned simply by choosing the crystallographic orientation. Meanwhile, the intense high-energy plasmonic screening and reflectivity behavior support the use of zirconia films as radiation-resistant shielding layers and thermal barrier coatings for aerospace components exposed to extreme electromagnetic and thermal loads.

The authors frame their work as a quantitative design framework rather than a purely descriptive exercise. By explicitly mapping the microscopic symmetry breaking of the cubic-to-monoclinic transition onto macroscopic directional optical tensors across an energy continuum extending to 40 electronvolts, the study moves beyond the isotropic averages that dominated earlier modeling and validates its predictions against experimental benchmarks at every stage. Future extensions, the researchers suggest, could incorporate finite-temperature lattice dynamics to track how thermal expansion modulates anisotropy in operating devices, simulate oxygen vacancies and non-stoichiometry relevant to resistive switching memories and catalysis, and model doped systems such as yttria-stabilized zirconia. For now, the message is clear: in zirconia, the geometry of a handful of atoms around each zirconium site is not a crystallographic footnote but the master switch controlling how the material interacts with light, and engineers who learn to exploit it gain a powerful new lever for the photonics and extreme-environment technologies of the next decade.

Subject of Research: First-principles investigation of crystal symmetry effects on the electronic structure and optical anisotropy of cubic and monoclinic zirconium dioxide polymorphs

Article Title: Unveiling the role of crystal symmetry in the electronic structure and optical anisotropy of ZrO 2 polymorphs: a first-principles study

Article References: Alsubaie, A. S., & Mahmoud, K. H. (2026). Unveiling the role of crystal symmetry in the electronic structure and optical anisotropy of ZrO2 polymorphs: a first-principles study. Results in Physics, Article 108718. https://doi.org/10.1016/j.rinp.2026.108718

Image Credits: AI Generated

DOI: 10.1016/j.rinp.2026.108718

Keywords: zirconium dioxide, zirconia polymorphs, crystal symmetry, density functional theory, optical anisotropy, band gap, dielectric function, birefringence, UV photodetectors, thermal barrier coatings, first-principles calculations, optoelectronics