No abstract

The vibrational energy of crystals is known to propagate in a fixed number of phonon branches allowed by symmetry. In the realm of nonlinear dynamics, however, additional nonlinear propagating modes are possible. Nonlinear propagating modes have unique properties that are important in many disciplines including optical communications, conducting polymers, biology, magnetism, and nuclear physics. Yet, despite the crucial importance of crystal lattice vibrations in fundamental and applied science, such additional propagating modes have not been observed in ordinary crystals. Here, we show that propagating modes exist beyond the phonons in fluorite-structured thoria, urania, and natural calcium fluoride using neutron scattering and first-principles calculations. These modes are observed at temperatures ranging from 5 K up to 1200 K, extend to frequencies 30–40% higher than the maximum phonon frequency, and travel at velocities comparable to or higher than the fastest phonon. The nonlinear origin of the modes is explained in part via perturbation theory, which approximately accounts for nonlinearity. Given that these modes are still clearly observed at 5 K, they are likely an inherent feature of the quantum ground state. The existence of these waves in three-dimensional crystals may have ramifications for material properties.

Phonons and their interactions are necessary for determining a wide range of materials properties. Here we present four independent advances which facilitate the computation of phonons and their interactions from first-principles. First, we implement a group-theoretical approach to construct the order N Taylor series of a d-dimensional crystal purely in terms of space group irreducible derivatives (ID), which guarantees symmetry by construction and allows for a practical means of communicating and storing phonons and their interactions. Second, we prove that the smallest possible supercell which accommodates N given wavevectors in a d-dimensional crystal is determined using the Smith Normal Form of the matrix formed from the corresponding wavevectors; resulting in negligible computational cost to find said supercell, in addition to providing the maximum required multiplicity for uniform supercells at arbitrary N and d. Third, we develop a series of finite displacement methodologies to compute phonons and their interactions which exploit the first two developments: lone and bundled irreducible derivative (LID and BID) approaches. LID computes a single ID, or as few as possible, at a time in the smallest supercell possible, while BID exploits perturbative derivatives for some order less than N (e.g. Hellman-Feynman forces) in order to extract all ID in the smallest possible supercells using the fewest possible computations. Finally, we derive an equation for the order N volume derivatives of the phonons in terms of the order N = N + 2 ID. Given that the former are easily computed, they can be used as a stringent, infinite ranged test of the ID. Our general framework is illustrated on graphene, yielding irreducible phonon interactions to fifth order. Additionally, we provide a cost analysis for the rock-salt structure at N = 3, demonstrating a massive speedup compared to popular finite displacement methods in the literature.