Adiabatic Dynamics of Coupled Spins and Phonons in Magnetic Insulators

Shang Ren,John Bonini,Massimiliano Stengel,Cyrus E. Dreyer,David Vanderbilt,and David Vanderbilt

DOI: https://doi.org/10.1103/physrevx.14.011041

2024-03-08

Physical Review X

Abstract:In conventional ab initio methodologies, phonons are calculated by solving equations of motion involving static interatomic force constants and atomic masses. The Born-Oppenheimer approximation, where all electronic degrees of freedom are assumed to adiabatically follow the nuclear dynamics, is also adopted. This approach does not fully account for the effects of broken time-reversal symmetry in systems with magnetic order. Recent attempts to rectify this involve the inclusion of the velocity dependence of the interatomic forces in the equations of motion, which accounts for time-reversal symmetry breaking, and can result in chiral phonon modes with nonzero angular momentum even at the zone center. However, since the energy ranges of phonons and magnons typically overlap, the spins cannot be treated as adiabatically following the lattice degrees of freedom. Instead, phonon and spins must be treated on a similar footing. Focusing on zone-center modes, we propose a method involving Hessian matrices and Berry curvature tensors in terms of both phonon and spin degrees of freedom, and describe a first-principles methodology for calculating these. We then solve Lagrange's equations of motion to determine the energies and characters of the mixed excitations, allowing us to quantify, for example, the energy splittings between chiral pairs of phonons in some cases, and the degree of magnetically induced mixing between infrared and Raman modes in others. The approach is general and can be applied to determine the adiabatic dynamics of any mixed set of slow variables. https://doi.org/10.1103/PhysRevX.14.011041 Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI. Published by the American Physical Society

physics, multidisciplinary

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is: **How to systematically deal with the coupling dynamics between phonons and magnons (spin waves) in magnetic insulators?** Specifically, the authors focus on the dynamic behaviors and interactions of phonons and magnons in magnetic - ordered systems that break time - reversal symmetry (TRS).

### Detailed Explanation

1. **Background and Challenges**
- In traditional ab - initio methods, phonons are calculated by solving the equations of motion involving static inter - atomic force constants and atomic masses. These methods usually adopt the Born - Oppenheimer approximation, assuming that all electronic degrees of freedom can adiabatically follow the motion of nuclei.
- However, this traditional method cannot fully describe the influence of time - reversal symmetry breaking in magnetic - ordered systems. For example, in magnetic materials, the breaking of TRS will cause phonon modes to exhibit chiral characteristics, that is, phonon modes with non - zero angular momentum at the center of the Brillouin zone.
- In addition, since the energy ranges of phonons and magnons usually overlap, the spin cannot be simply regarded as adiabatically following the lattice degrees of freedom. Therefore, it is necessary to deal with phonons and spins simultaneously and consider the coupling between them.

2. **Proposed Method**
- The authors proposed a method based on the Hessian matrix and Berry curvature tensor for simultaneously dealing with phonon and spin degrees of freedom. They determine the energy and characteristics of mixed excitations through Lagrange equations.
- Specifically, they introduced the Hessian matrix and Berry curvature tensor containing phonon and spin degrees of freedom and calculated these matrices by first - principles methods.
- Through this method, the authors were able to quantify, for example, the energy splitting of chiral phonon pairs and the degree of mixing between magnetically - induced infrared and Raman modes.

3. **Significance of the Research**
- This method provides a general theoretical framework, which is not only applicable to coupled phonon and spin systems but can also be extended to the adiabatic dynamics of any mixed slow - variable systems.
- It significantly surpasses previous approximate methods, such as the minimal spin - phonon model, which only considers part of the Berry curvature contributions.

4. **Application and Verification**
- The authors carried out four case studies, covering three - dimensional and two - dimensional ferromagnetic and antiferromagnetic materials. They analyzed the energy splitting and symmetry of chiral phonons in these systems.
- The results show that the energy splitting of chiral phonons depends on the point - group symmetry in the presence of magnetic order, and high - order terms need to be considered to correctly capture the situation where magnetic order breaks spatial symmetry.

### Summary
This paper solves the complex problem of phonon - magnon coupling dynamics in magnetic insulators by introducing a new theoretical framework and calculation method, providing an important tool for understanding the physical properties of these systems.