Magnetic moments induce strong phonon
renormalization in FeSi

S. Krannich, Yvan Sidis, D. Lamago, R. Heid, J.-M. Mignot, H. Löhneysen, A.
Ivanov, P. Steffens, T. Keller, L. Wang, et al.

► To cite this version:

S. Krannich, Yvan Sidis, D. Lamago, R. Heid, J.-M. Mignot, et al. Magnetic moments in- duce strong phonon renormalization in FeSi. Nature Communications, 2015, 6 (1), pp.8961. 10.1038/ncomms9961. cea-01820491

HAL Id: cea-01820491

https://cea.hal.science/cea-01820491v1

Submitted on 22 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE

Received 20 Feb 2015 | Accepted 21 Oct 2015 | Published 27 Nov 2015

DOI: 10.1038/ncomms9961

OPEN

# Magnetic moments induce strong phonon renormalization in FeSi

S. Krannich<fnref n="1" />, Y. Sidis², D. Lamago<fnref n="1" />,², R. Heid<fnref n="1" />, J.-M. Mignot², H.v. Löhneysen<fnref n="1" />,³, A. Ivanov⁴, P. Steffens⁴, T. Keller⁵,⁶, L. Wang<fnref n="1" />, E. Goering⁷ & F. Weber<fnref n="1" />

The interactions of electronic, spin and lattice degrees of freedom in solids result in complex phase diagrams, new emergent phenomena and technical applications. While electron-phonon coupling is well understood, and interactions between spin and electronic excitations are intensely investigated, only little is known about the dynamic interactions between spin and lattice excitations. Noncentrosymmetric FeSi is known to undergo with increasing temperature a crossover from insulating to metallic behaviour with concomitant magnetic fluctuations, and exhibits strongly temperature-dependent phonon energies. Here we show by detailed inelastic neutron-scattering measurements and *ab initio* calculations that the phonon renormalization in FeSi is linked to its unconventional magnetic properties. Electronic states mediating conventional electron-phonon coupling are only activated in the presence of strong magnetic fluctuations. Furthermore, phonons entailing strongly varying Fe-Fe distances are damped via dynamic coupling to the temperature-induced magnetic moments, highlighting FeSi as a material with direct spin-phonon coupling and multiple interaction paths.

<footnote n="1">¹Institute for Solid State Physics, Karlsruhe Institute of Technology, 76021 Karlsruhe, Germany. ²Laboratoire Léon Brillouin (CEA-CNRS), CEA Saclay, F-91911 Gif-sur-Yvette, France. ³Physics Institute, Karlsruhe Institute of Technology, 76049 Karlsruhe, Germany. ⁴Institute Laue-Langevin, BP 156, 38042 Grenoble, France. ⁵Max Planck Institute for Solid State Research, D-70569 Stuttgart, Germany. ⁶Max Planck Society, Outstation at the FRM II, D-85748 Garching, Germany. ⁷Max Planck Institute for Intelligent Systems, 70569 Stuttgart, Germany. Correspondence and requests for materials should be addressed to F.W. (email: frank.weber@kit.edu).</footnote>
The family of transition metal silicides crystallizing in the noncentrosymmetric cubic B20 structure exhibits a variety of complex phenomena, which are of interest to both basic and applied science¹⁻³. In MnSi, the helical arrangement of Mn spins results in exotic phases such as partial order under hydrostatic pressure⁴ or a skyrmion lattice in an applied magnetic field⁵. The latter is also found in Fe₁₋ₓCoₓSi (ref. 6). Pure FeSi has been under intense research for over 40 years owing to its puzzling non-magnetic ground state⁷,⁸ and recent experiments opened the perspective for potential applications in spintronics⁹ and thermoelectricity¹⁰,¹¹.

FeSi, a semiconductor with a small bandgap $\Delta$ of 60–70 meV (refs 7,8,12), shows similarities to the class of heavy-fermion Kondo insulators. Its foremost peculiarity is a large temperature-induced paramagnetic moment of > 2 $\mu_B$, while the ground state is non-magnetic⁷,¹³. The mechanism responsible for the temperature-activated magnetization and the strong temperature dependence of the bandgap in FeSi has remained controversial ever since its discovery back in the 1960s (refs 7,8). Inelastic polarized neutron scattering showed a sharp increase with increasing temperature of magnetic scattering close to reciprocal lattice points¹⁴,¹⁵, indicating the existence of ferromagnetic correlations¹⁶. Recently, thermal lattice disorder¹⁷, as well as strong electron-electron correlations¹⁸, has been invoked to explain the unusual properties of FeSi. In particular, the strong reduction of the gap $\Delta$ on heating towards room temperature¹² is observed in an unexpectedly low-temperature range with respect to the gap value of 60–70 meV.

Although the average atomic structure of FeSi does not exhibit any notable changes with temperature, anomalous phonon properties have been reported recently from Raman¹² and inelastic neutron-scattering measurements¹⁷. The authors of those publications found a surprisingly strong softening of phonon energies with increasing temperature, which they related to the closing of the charge gap in FeSi. However, the information obtained was limited because of restrictions in the accessible range of temperature and wave vector of Raman measurements. Furthermore, most inelastic neutron-scattering data were obtained on powder samples, not allowing an investigation of the lifetimes of phonon modes, or their directional dependences in reciprocal space.

We have performed inelastic neutron-scattering experiments on a high-quality single crystal of FeSi. This allows us to determine phonon dispersion curves. In addition, we obtain detailed information on the temperature dependence of the phonon linewidths that are directly linked to the lifetimes of lattice vibrations, across the transition from a non-magnetic semiconductor to a metal with a large paramagnetic moment. Large temperature-dependent intrinsic linewidths have been observed for several phonon modes, which simultaneously show decreasing coherent intensities and very strongly temperature-dependent energies. We find that the temperature dependences of linewidths and intensities agree well with that of the temperature-induced magnetic moment. *Ab initio* lattice dynamical calculations demonstrate that reduced phonon lifetimes can be understood only if we include the evolution of magnetism. In conjunction with this behaviour, the reduction of phonon intensities goes far beyond the quasi-harmonic approximation. We argue that the particular ionic motions associated with the investigated phonons can lead to additional renormalization as temperature increases, via a dynamical coupling to the growing magnetic moments.

## Results

To establish an understanding of the lattice dynamical properties of FeSi in detail, we performed extensive measurements of the

phonon modes along various high-symmetry directions in reciprocal space at $T=10$ K. Phonons were primarily investigated using inelastic neutron scattering on the thermal triple-axis spectrometer (TAS) 1T at the ORPHEE reactor at Laboratoire Léon Brillouin at CEA Saclay and results were compared with *ab initio* calculations based on density functional perturbation theory (DFPT). To that end, we used the experimental low-temperature lattice constants¹⁹ and optimized the internal parameter of the unit cell to achieve a force-free equilibrium structure. In the following, all wave vectors are expressed in reciprocal lattice units of $2\pi/a$, where $a=4.48$ Å is the low-temperature lattice constant of the cubic unit cell of FeSi. Overall, we find very good agreement between the observed and calculated energies (Fig. 1a; Supplementary Fig. 1).

We have studied the phonon properties over a wide range of temperature ($10 \text{ K} \le T \le 790 \text{ K}$) at various wave vectors for which phonon selection rules are unambiguous. Note that the rather low symmetry of the B20 crystal structure typically results in unclear phonon selection rules, that is, various phonon modes contribute to the spectral weight distribution in the corresponding energy scans. Accordingly, the linewidths of the different phonon peaks can no longer be accurately determined. Raw data taken at $q = (1/2, 1/2, 1/2)$, that is, at the Brillouin zone boundary **R** point along the [111] direction, are shown in Fig. 1b,c. We find two well-separated phonon peaks, in agreement with the DFPT predictions, at energies of 24.3 and 34.2 meV for $T=10$ K. For simplicity, we call them R1 and R2 modes, respectively, in the following. We observe a very strong softening of the phonon energies, by 14% for both modes, in the investigated temperature range, which is in good agreement with the reported temperature dependence of a peak observed previously in the phonon density of states¹⁷. The detailed study of the present work, however, demonstrates that the softening exhibits quite unusual temperature dependence. Over the whole temperature range between 10 and 790 K, we observe three distinct regions

<figure><img src="image_2.png" /></figure>

**Figure 1 | Representative inelastic neutron-scattering data and dispersion at the R point.** (a) Measured and calculated phonon dispersion in FeSi from the centre of the Brillouin zone towards the **R** point, **Q** = (1.5, 1.5, 1.5), which is the zone boundary in the crystallographic [111] direction. Shown are branches with longitudinal symmetry. (b,c) Measured neutron intensities at the **R** point at temperatures (b) $T=790$ K and (c) $T=10$ K. Solid lines denote fits to the data using a Lorentzian line shape convoluted with the experimental resolution on top of the experimental background (dotted horizontal lines). Dashed (red) vertical lines indicate the observed peak energies and mark the softening of the R1 and R2 modes between $T=10$ and 790 K by 14% of the low-temperature values. Thick horizontal bars denote the calculated experimental resolution. Error bars represent s.d.
<figure><img src="image_3.png" /><figcaption>**Figure 2** | Temperature dependence of the phonon properties at the R point. (a,b) Temperature evolution of the energy of the (a) R1 and (b) R2 phonon modes. The dashed lines denote the calculated energies from the quasi-harmonic approximation using temperature-dependent lattice constants reported in ref. 19 and simulating the effect of thermal expansion (th. exp.). (c) Temperature evolution of intensities of the R1 and R2 modes corrected for the Bose factor and normalized to their values at $T=10$ K. The dashed line indicates the expected evolution because of the temperature-dependent Debye-Waller factor [19] (DW). (d) Temperature dependence of the observed linewidths $\Gamma_{\text{exp}}$ of the R1 and R2 modes. Error bars represent s.d.</figcaption></figure>

(Fig. 2a,b): (1) the temperature dependence is rather weak for $T \le 100$ K. In our calculations, we show that in this temperature range, the softening can be quantitatively attributed to the increase in the lattice constants using values obtained by neutron diffraction¹⁹ (see dashed lines in Fig. 2a,b). (2) At intermediate temperatures ($100 \text{ K} \le T \le 300 \text{ K}$), we observe a strong anomalous phonon softening, with a factor of ~4 stronger decrease of the phonon energy than predicted based on the thermal expansion (see dashed line in Fig. 2a,b). (3) Finally, for $300 \text{ K} \le T \le 790 \text{ K}$, a weaker softening rate is recovered, with a decreasing energy consistent with the measured lattice expansion. We note that many other phonon modes show qualitatively similar temperature dependences, although the strength of the softening varies. There are also some modes for which the energies fully follow the expectations from thermal expansion without an additional softening. Thus, the anomalous softening is related to specific phonon symmetries. A complete report is beyond the scope of this paper and will be published in the future.

We can further analyse the temperature dependence of the phonon intensities and linewidths at wave vectors with unambiguous phonon selection rules. In particular, we can accurately determine the experimental background at the R point, and derive the energy-integrated intensities and the linewidths $\Gamma_{\text{exp}}$ (full width at half maximum) in a first step from Gaussian fits to the raw data of the R1 and R2 modes. Both quantities show anomalous behaviour, again in the similar temperature range 200–400 K. After correcting for the temperature-dependent phonon intensities due to the Bose factor, we find that the energy-integrated intensities of both R1 and R2 show a rapid decrease by 30% in the range $200 \text{ K} \le T \le 400 \text{ K}$, and then level off towards higher $T$ (Fig. 2c). At low temperature, $\Gamma_{\text{exp}}$ agrees well with the calculated experimental resolution of 3 and 3.2 meV for the R1 and R2 modes, respectively. Thus, the intrinsic phonon linewidth $\Gamma_{\text{phon}}$ cannot be resolved in this $T$ range. At higher temperatures, at which $\Gamma_{\text{exp}}$ deviates from the instrumental resolution, we can extract $\Gamma_{\text{phon}}$ by fitting a Lorentzian function, which is typically a good approximation of the spectral function of a damped phonon as long as $\Gamma_{\text{phon}}$ is much smaller than the mode energy $\hbar\omega$, convoluted with the Gaussian experimental resolution to the data (Fig. 1b). For $T \ge 200 \text{ K}$, the linewidth increases strongly and levels off at 500 K (Fig. 2d). The temperature dependences are qualitatively similar for both R1 and R2, with an increase of $\Gamma_{\text{exp}}$ by 30 and 50%, respectively. Generally, phonons exhibit a reduced intensity with increasing temperature due to thermal atomic motions, which reduce the

phonon intensity according to the Debye-Waller factor²⁰. Furthermore, incoherent atomic motions hamper an easy propagation of phonons through the lattice and, hence, lead to smaller phonon lifetimes, that is, larger linewidths. However, the observed strong temperature dependences cannot be explained by these effects. First, the observed temperature dependence of the Debye-Waller factor¹⁹ only accounts for an intensity loss by about 4% (see dashed line in Fig. 2c) in the temperature range investigated. Second, the temperature dependence of the Debye-Waller factor is smooth, in stark contrast to our observations¹⁹.

As already mentioned, we can observe finite values of $\Gamma_{\text{phon}}$ for the two phonon modes R1 and R2 only for $T \ge 300 \text{ K}$ because statistical error bars on $\Gamma_{\text{exp}}$ are about $\pm (0.15-0.25)$ meV due to the relatively low experimental resolution of the TAS. To investigate the evolution of $\Gamma_{\text{phon}}$ at lower temperatures, we performed neutron resonant spin-echo measurements on the polarized thermal neutron TRISP spectrometer at the Heinz-Maier-Leibnitz Neutron Research Facility (FRM-II) in Garching. In this technique, the precession of the neutron spin in magnetic fields before and after being scattered in the sample is used to achieve very high energy resolution²¹⁻²³. The resulting temperature dependence of $\Gamma_{\text{phon}}$ for the R1 mode in the range $15 \text{ K} \le T \le 230 \text{ K}$ is shown in Fig. 3a along with results from the TAS measurements for $T \ge 300 \text{ K}$. In the former, a constant $\Gamma_0 = 0.35$ meV, which we assign to anharmonic interactions and crystal imperfections, was subtracted from the neutron resonant spin-echo results. We find that the increase of $\Gamma_{\text{phon}}$ starts just above $T=100 \text{ K}$, that is, well below room temperature, in agreement with the onset of the strong softening and levels off at high temperatures.

In the following, we will show that the observed phonon anomalies can be related to the magnetism in FeSi, which has long been known to exhibit an unusual temperature dependence on the same $T$ region.

We investigated the magnetic degrees of freedom via polarized inelastic neutron-scattering and magnetization measurements. Results from the latter are shown in Supplementary Fig. 2 and agree well with previous reports⁷,⁸. In Fig. 3a, we plot the (scaled) square of the temperature-induced magnetic moment $M^2(T)$ calculated via $M^2(T) = 3k_B T \times \chi(T)$ and compare it with $\Gamma_{\text{phon}}(T)$ of the R1 and R2 phonon modes. We find good agreement between the increase of $M^2(T)$ and $\Gamma_{\text{phon}}(T)$, indeed indicating a close relationship between magnetism and the lattice dynamical properties in FeSi. Pramagnetic scattering obtained via polarized
<figure><img src="image_4.png" /><figcaption>**Figure 3** | Temperature dependences of lattice and magnetic degrees of freedom in FeSi. (a) Phonon linewidths $\Gamma_{\text{phon}}$ (full width at half maximum) of the R1 and R2 modes obtained from TAS ($T \ge 300$ K) and neutron resonant spin-echo (NRSE; $T \le 230$ K) measurements. Solid lines represent $M^2 = 3k_B T\chi(T)$ based on our magnetization measurements and scaled to fit the high temperature values of $\Gamma_{\text{phon}}$. Dashed arrows indicate the predicted linewidths for the R1 and R2 modes in spin-polarized DFPT calculations. (b) Paramagnetic scattering measured at $15 \text{ K} \le T \le 300 \text{ K}$ employing the differential method. Data are vertically offset for clarity. Colour-coded dashed lines denote the respective zero lines. Solid lines are fits using equation (1). (c) Comparison of the temperature dependence of $\Gamma_{\text{phon}}$ of the R1 and R2 modes with that of the temperature-induced paramagnetic moment $M^2(T)$ extracted from polarized neutron scattering (see b). Values of $M^2(T)$ at $T \ge 300 \text{ K}$ were taken from ref. 15 and scaled using the value at $T = 300 \text{ K}$. All data are normalized to their respective value at $T = 600 \text{ K}$. (d) Comparison of the reduction of the intensities of the R1 and R2 modes with increasing temperature with $M^2(T)$ derived from polarized neutron scattering. Again, all data are normalized to their respective value at $T = 600 \text{ K}$. Error bars represent s.d.</figcaption></figure>

inelastic neutron scattering close to the (1, 1, 0) Bragg position is shown in Fig. 3b (Supplementary Fig. 3; Supplementary Note 1). The data demonstrate that magnetic fluctuations present at room temperature vanish completely at a temperature between 15 and 150 K. In agreement with previous reports$^{14,15}$, the observed scattering can be described by the paramagnetic scattering function

$$S(\mathbf{q}, \omega) = M^2(q) \frac{\hbar\omega/k_B T}{1 - \exp(-\hbar\omega/k_B T) \pi \Gamma^2 + (\hbar\omega)^2} \quad (1)$$

using a single value for the energy width $\Gamma_{\text{magn}} = (10 \pm 1) \text{ meV}$ of the paramagnetic scattering at all temperatures (Fig. 3b). Hence, the increase of magnetic scattering with temperature is qualitatively the same at all energies.

We compare $\Gamma_{\text{phon}}(T)$ of the R1 and R2 modes with $M^2(T)$ extracted from polarized neutron scattering (Fig. 3c) and with the intensity reduction of both phonon modes (Fig. 3d). Here we include results presented in ref. 15 for $M^2(T)$ at $T \ge 300 \text{ K}$. The absolute values were scaled to match our result at room temperature.

We find good agreement in both comparisons demonstrating that also the intensity reduction is linked to the evolution of magnetism in FeSi. At this point, we emphasize that neutron diffraction data do not reveal any significant magneto-elastic coupling, that is, the lattice constant increases smoothly on heating$^{19}$. Hence, any coupling between phonons and magnetism must be dynamic in origin.

We now return to the *ab initio* lattice dynamical calculations that yield a smooth evolution of the phonon energies R1 and R2

when taking only the thermal lattice expansion into account (Supplementary Note 2). Besides phonon energies, we can also calculate the electronic contribution to the phonon linewidth $\gamma$, that is, the coupling strength of a specific phonon mode to the electronic states. Because of the charge gap, there is no Fermi surface and, hence, no electron-phonon coupling (EPC) in FeSi at low temperatures. Indeed, our electronic structure calculations show a strong suppression of the electronic density of states near $E_F$, which is reminiscent of a gap (Supplementary Fig. 4). The charge gap in FeSi decreases strongly on heating towards room temperature$^{12}$ and hence EPC can be present in the metallic regime. However, calculations performed for metallic FeSi (Supplementary Note 2) yield linewidths smaller than 0.1 meV for all phonons, much smaller than the observed values (Fig. 3a). We conclude that a mere metallization of FeSi cannot explain our observations.

Motivated by the very similar temperature evolution of $\Gamma_{\text{phon}}(T)$ and $M^2(T)$, we have performed spin-polarized calculations. Although DFPT, in agreement with experiment, does not yield a stable magnetically ordered ground state for FeSi, we found that such a state can be stabilized by increasing the lattice constant to $a^* = 4.65 \AA$, that is, by nearly 4% with respect to the experimental lattice constant. Spin-polarized calculations with $a^* = 4.65 \AA$ yield a ferromagnetically ordered ground state with an ordered moment of $0.8 \mu_B$ per Fe atom (Supplementary Fig. 4). Although FeSi at elevated temperatures is a paramagnet, the assumption of ferromagnetic short-range correlations is not unreasonable, the scattering due to magnetic fluctuations strongly peaks on approaching reciprocal lattice vectors with large ferromagnetic structure factor$^{15}$, which is a clear signature of
short-range ferromagnetic correlations (Supplementary Fig. 3; Supplementary Note 1)¹⁶. Indeed, properties of FeSi have been interpreted in terms of a nearly ferromagnetic semiconductor²⁴. Finally, closely related compounds such as FeGe show helical magnetic order with a small pitch between neighbouring Fe moments, which on short-length scales resembles ferromagnetic order²⁵ (Supplementary Note 2).

Results for the lattice dynamics for non-magnetic and ferromagnetic FeSi, both with $a^* = 4.65$ Å, are shown for two high-symmetry directions in Fig. 4a,b. Compared with the dispersion presented in Fig. 1a, energies are generally softer as a result of the increased lattice constant. The non-magnetic calculations with different lattice constants show different phonon energies by about 25% (as can be seen by comparing Fig. 4a ($a^* = 4.65$ Å) and Fig. 1a ($a = 4.48$ Å)). Importantly, introducing magnetism in FeSi with the same lattice constant $a^*$ causes additional softening by up to 10%, which is observed at the Brillouin zone centre and over the whole wave vector range, for example, along the [100] direction (Fig. 4b). The magnetism-induced effect on phonons is particularly strong for the R1 mode, with a softening of 35%, and still relatively large for the R2 mode (15%). Furthermore, the spin-polarized calculation predicts strong overall EPC, with an EPC constant of $\lambda \approx 1$, which is similar to values obtained for conventional superconductors with medium to strong EPC²⁶,²⁷.

Accordingly, several phonon branches, in particular the R1 and R2 modes, are strongly damped and exhibit large values of the calculated linewidth $\gamma$ due to EPC (Fig. 4a,b). We were able to investigate two more phonon modes, one at $q = (0.25, 0.25, 0)$

<figure><img src="image_5.png" /><figcaption>**Figure 4** | Results from spin-polarized DFPT calculations. (a,b) Calculated phonon dispersions using $a^* = 4.655$ Å with (solid lines) and without (dashed lines) ferromagnetic spin polarization along the (a) [111] and (b) [100] directions. For clarity, only branches with longitudinal symmetry are shown. The calculated electronic contribution $\gamma$ to the phonon linewidth is represented as error bars in the data of the spin-polarized calculation. In the non-magnetic calculation, results for $\gamma$ are <0.1 meV for all phonon modes. (c,d) Unit cell of FeSi with the calculated eigenvector of the (c) R1 and (d) R2 modes indicated by arrows (qualitative). Coordinates in real space are given in lattice constants units (a, a, a). r.l.u., reciprocal lattice unit.</figcaption></figure>

and the other at the zone centre. Spin-polarized DFPT predicts essentially zero EPC for the former mode but very strong EPC for the latter. Both predictions are verified by our measurements (Supplementary Fig. 5; Supplementary Note 3). Our analysis therefore indicates that temperature-induced magnetism is directly responsible for the increasing phonon linewidths, which can account for the close similarity between $\Gamma_{\text{phon}}(T)$ and $M^2(T)$. We note that our calculations do not include a direct spin-phonon coupling (SPC). Rather, the large linewidths in spin-polarized DFPT originate from a change in the electronic states at the Fermi surface, which are mediating EPC.

## Discussion

On the other hand, our calculations underestimate the absolute value of the linewidths by 30% for the R1 mode and even a factor of three for the R2 mode and do not predict the intensity reduction seen in Fig. 2c. Large phonon linewidths, for which DFPT correctly predicted the wave vector dependence but underestimated the absolute strength, have also been found in some charge density wave compounds²⁸,²⁹ and superconductors²⁶. The likely origin of the discrepancies, in those cases, is that DFPT does not include anharmonic effects, which definitely play an important role in charge density wave materials. For FeSi, an explanation might include a direct SPC. The impact of SPC on the phonon lifetime in a Heisenberg ferromagnet has been calculated, where the lowest-order contribution to the phonon lifetime originates from processes, in which the phonon emits and later reabsorbs two magnons³⁰. In an ordered ferromagnet, the strength of the renormalization depends sensitively on the overlap of the wave vector and energy position of the renormalized phonon with those of a bound state of the two emitted spin waves. In contrast to well-defined spin waves, magnetic intensities in FeSi show qualitatively the same temperature dependences independent of their respective energy (Fig. 3b). Hence, one would expect the same temperature dependences for phonon renormalization effects due to SPC independent of phonons' respective energies. Thus, this scenario of phonons coupling to magnetic excitations naturally explains the concomitant increase of magnetic scattering and phonon renormalization in FeSi (Fig. 3c).

Our results show that there are also phonons, which do not exhibit strong renormalization. Therefore, it is interesting to investigate the phonon displacement patterns to look for common features in modes with strong renormalizations and differences with regard to phonon modes showing no response to the growing magnetism. We have computed the atomic displacement patterns for the R1 and R2 modes (Fig. 4c,d), as well as for the other two modes discussed above (Supplementary Fig. 5) from our DFPT calculations. Indeed, we find that Fe-Fe distances do strongly change for the modes exhibiting large linewidths at high temperatures, whereas the pattern of the mode at $q = (0.25, 0.25, 0)$ with zero EPC features motions of the Fe atoms in the unit cell with much smaller changes in the Fe-Fe distances. This difference directly reflects the strong interaction between magnetic exchange that decisively depends on the Fe-Fe distance and the lattice degrees of freedom. Another interesting scenario considers the coupling between the time-dependent electric polarization induced by vibrating ions and the magnetic moment of the Fe 3d electrons. Although, this situation has—to our knowledge—not yet been considered theoretically, we would also qualitatively expect that phonons, whose eigenvectors exhibit strong changes in the Fe-Fe distances, show the strongest effects. Hence, FeSi appears to be a rare example of SPC (in contrast to the rather common effects on phonons because of magnetoe-elastic coupling). An increasing anharmonicity due to SPC could
explain the reduction of the phonon intensities that is also similar to the temperature dependence of $M^2(T)$ (Fig. 3d).

In conclusion, we have reported temperature-dependent phonon renormalization in FeSi over a wide temperature range. Our analysis using DFPT calculations demonstrates that our observations are directly related to the temperature-induced magnetism. Going beyond the quasi-harmonic approximation, our results suggest that FeSi exhibits a purely dynamic SPC, which enhances the renormalization predicted by DFPT. Therefore, FeSi is not only a very peculiar paramagnet but may serve as a model compound for strongly interacting magnetic and lattice degrees of freedom.

## Methods

**Sample characterization.** Our single crystal of FeSi weighing 35 g was characterized in terms of magnetization and electrical resistivity, as well as Larmor diffraction and crystal mosaicity using neutrons (for results see Supplementary Fig. 2). Magnetization measurements have been performed by a new Quantum Design MPMS3 system, utilizing the 1,000-K oven option above 300 K. The measurements have been performed in vibrating sample magnetometer (VSM) mode with 1-mm peak-to-peak amplitude. The VSM data have been subsequently normalized to room temperature d.c. mode results, to obtain high-quality absolute magnetization values. No further data processing has been performed. The sample was a piece cut from the large single crystal used for neutron scattering with $m=0.033$ g and dimensions $(4 \times 2 \times 1)$ mm$^3$. Standard Four-Point Probe resistivity measurement has been performed using the physical properties measurements system (PPMS) a.c. transport function to control the temperature from 5 to 300 K and excitation current to be 0.1 mA. The sample cut from the large single crystal used for neutron scattering had the dimensions $(4 \times 1 \times 1)$ mm$^3$. The contacts are made by connecting 20-µm Ag wires using conductive silver paint. The distance between the two leads for the voltage measurement was 3.5 mm. Therefore, we calculated the resistivity, $\rho = (\frac{V}{I}) \times (\frac{A}{V})$ (V: voltage; I: excitation current; A: area, that is, $1 \times 1$ mm$^2$; l: distance between leads).

**Inelastic neutron scattering.** Inelastic neutron scattering was used to study the phonons on the thermal TAS 1T at the ORPHEE reactor at Laboratoire Léon Brillouin, CEA Saclay. Double-focusing graphite monochromator and analyser were used for all phonon measurement. The final energy at the analyser was set to 14.7 meV allowing the use of a graphite filter to suppress higher-order scattering. The single crystal was mounted in a closed-cycle refrigerator allowing measurements in the temperature range $5 K \le T \le 790$ K. The paramagnetic scattering in FeSi was studied with polarized thermal neutrons on the TAS IN20 at ILL, Grenoble, using the differential method$^{31}$ with a standard set-up employing Heusler monochromator and analyser and a final energy of 14.7 meV. The intrinsic phonon linewidths $\Gamma_{\text{phon}}$ of the R1 mode at $15 K \le T \le 225$ K were investigated using the neutron resonance spin-echo technique on the polarized thermal TAS TRISP at the neutron research source Hans-Maier Leibnitz (FRM-II) in Garching$^{21,22}$ (for details of the measurement see Supplementary Note 4 and Supplementary Fig. 6). The same instrument was used to perform Larmor diffraction measurements$^2$.

**Density functional theory.** Density functional theory calculations were performed in the framework of the mixed basis pseudopotential method (Meyer, B., Elsässer, C., Lechermann, F. & Fähnle, M. FORTRAN90 Program for Mixed-Basis Pseudopotential Calculations for Crystals (Max-Planck-Institut fur Metall-forschung, Stuttgart)). The exchange-correlation functional was treated in the local density approximation in the Perdew-Wang parameterization. Norm-conserving pseudopotentials for Fe and Si were constructed including the Fe 3s and 3p semicore states in the valence space. The deep potential can be efficiently treated in the mixed-basis scheme, which combines local functions together with plane waves for the representation of the valence states. Local functions of s, p and d symmetry at the Fe sites were combined with plane waves up to 22 Ry. Phonon frequencies and EPC were calculated using the linear response technique or DFPT$^{32}$ in combination with the mixed-basis pseudopotential method$^{33}$. Brillouin zone integrations were performed with a cubic $8 \times 8 \times 8$ k-point mesh (24 points in the irreducible Brillouin zone) in combination with a standard smearing technique using a Gaussian broadening of $\sigma = 50$ meV. Further calculations were performed using $\sigma = 200$ meV and including spin polarization.

## References

1. Manyala, N. et al. Magnetoresistance from quantum interference effects in ferromagnets. *Nature* **404**, 581–584 (2000).

2. Pfleiderer, C., Böni, P., Keller, T., Rössler, U. K. & Rosch, A. Non-Fermi liquid metal without quantum criticality. *Science* **316**, 1871–1874 (2007).

3. Schulz, T. et al. Emergent electrodynamics of skyrmions in a chiral magnet. *Nat. Phys.* **8**, 301–304 (2012).

4. Pfleiderer, C. et al. Partial order in the non-Fermi-liquid phase of MnSi. *Nature* **427**, 227–231 (2004).

5. Mühlbauer, S. et al. Skyrmion lattice in a chiral magnet. *Science* **323**, 915–919 (2009).

6. Münzer, W. et al. Skyrmion lattice in the doped semiconductor $Fe_1-xCo_xSi$. *Phys. Rev. B* **81**, 041203 (2010).

7. Jaccarino, V., Wertheim, G. K., Wernick, J. H., Walker, L. R. & Arajs, S. Paramagnetic excited state of FeSi. *Phys. Rev.* **160**, 476–482 (1967).

8. Wertheim, G. K. et al. Unusual electronic properties of FeSi. *Phys. Lett.* **18**, 89–90 (1965).

9. Manyala, N. et al. Large anomalous Hall effect in a silicon-based magnetic semiconductor. *Nat. Mater.* **3**, 255–262 (2004).

10. Tomczak, J., Haule, K. & Kotliar, G. in *New Materials for Thermoelectric Applications: Theory and Experiment* NATO Science for Peace and Security Series B: Physics and Biophysics (eds Zlatic, V. & Hewson, A.) Ch. 4, 45–57 (Springer, 2013).

11. Sales, B. C., Delaire, O., McGuire, M. A. & May, A. F. Thermoelectric properties of Co-, Ir-, and Os-doped FeSi alloys: Evidence for strong electron-phonon coupling. *Phys. Rev. B* **83**, 125209 (2011).

12. Menzel, D. et al. Electron-phonon interaction and spectral weight transfer in $Fe_1-xCo_xSi$. *Phys. Rev. B* **79**, 165111 (2009).

13. Mandrus, D., Sarrao, J. L., Migliori, A., Thompson, J. D. & Fisk, Z. Thermodynamics of FeSi. *Phys. Rev. B* **51**, 4763–4767 (1995).

14. Shirane, G., Fischer, J. E., Endoh, Y. & Tajima, K. Temperature-induced magnetism in FeSi. *Phys. Rev. Lett.* **59**, 351–354 (1987).

15. Tajima, K., Endoh, Y., Fischer, J. E. & Shirane, G. Spin fluctuations in the temperature-induced paramagnet FeSi. *Phys. Rev. B* **38**, 6954–6960 (1988).

16. Imada, M., Fujimori, A. & Tokura, Y. Metal-insulator transitions. *Rev. Mod. Phys.* **70**, 1039–1263 (1998).

17. Delaire, O. et al. Phonon softening and metallization of a narrow-gap semiconductor by thermal disorder. *Proc. Natl Acad. Sci. USA* **108**, 4725–4730 (2011).

18. Tomczak, J. M., Haule, K. & Kotliar, G. Signatures of electronic correlations in iron silicide. *Proc. Natl Acad. Sci. USA* **109**, 3243–3246 (2012).

19. Vočadlo, L., Knight, K. S., Price, G. D. & Wood, I. G. Thermal expansion and crystal structure of FeSi between 4 and 1173 K determined by time-of-flight neutron powder diffraction. *Phys. Chem. Minerals* **29**, 132–139 (2002).

20. Squires, G. L. Introduction to the Theory of Thermal Neutron Scattering (Dover Publications, 1996).

21. Aynajian, P. et al. Energy gaps and kohn anomalies in elemental superconductors. *Science* **319**, 1509–1512 (2008).

22. Habicht, K., Golub, R., Mezei, F., Keimer, B. & Keller, T. Temperature-dependent phonon lifetimes in lead investigated with neutron-resonance spin-echo spectroscopy. *Phys. Rev. B* **69**, 104301 (2004).

23. Keller, T. et al. Momentum-resolved electron-phonon interaction in lead determined by neutron resonance spin-echo spectroscopy. *Phys. Rev. Lett.* **96**, 225501 (2006).

24. Takahashi, Y. & Moriya, T. A theory of nearly ferromagnetic semiconductors. *J. Phys. Soc. Jpn* **46**, 1451–1459 (1979).

25. Lebech, B., Bernhard, J. & Freltoft, T. Magnetic structures of cubic FeGe studied by small-angle neutron scattering. *J. Phys. Condens. Matter* **1**, 6105–6122 (1989).

26. Weber, F. et al. Phonons and electron-phonon coupling in $YNi_2B_2C$. *Phys. Rev. B* **89**, 104503 (2014).

27. Weber, F. et al. Electron-phonon coupling in the conventional superconductor $YNi_2B_2C$ at high phonon energies studied by time-of-flight neutron spectroscopy. *Phys. Rev. Lett.* **109**, 057001 (2012).

28. Weber, F. et al. Electron-phonon coupling and the soft phonon mode in $TiSe_2$. *Phys. Rev. Lett.* **107**, 266401 (2011).

29. Weber, F. et al. Extended phonon collapse and the origin of the charge-density wave in 2H-Nb$Se_2$. *Phys. Rev. Lett.* **107**, 107403 (2011).

30. Silberglitt, R. Effect of spin waves on the phonon energy spectrum of a heisenberg ferromagnet. *Phys. Rev.* **188**, 786–792 (1969).

31. Shirane, G., Shapiro, S. & Tranquada, J. Neutron Scattering with a Triple-Axis Spectrometer (Cambridge Univ. Press, 2002).

32. Baroni, S., de Gironcoli, S., Dal Corso, A. & Giannozzi, P. Phonons and related crystal properties from density-functional perturbation theory. *Rev. Mod. Phys.* **73**, 515–562 (2001).

33. Heid, R. & Bohnen, K.-P. Linear response in a density-functional mixed-basis approach. *Phys. Rev. B* **60**, R3709–R3712 (1999).

## Acknowledgements

We acknowledge valuable discussion with L. Pintschovius, S.K. and F.W. were supported by the young investigator group VH-NG-840 of the Helmholtz society. T.K. was supported by the DFG through SFB/TRR80. We acknowledge support by Deutsche Forschungsgemeinschaft and Open Access Publishing Fund of Karlsruhe Institute of Technology.
**Author contributions**

Neutron scattering: S.K., Y.S., D.L., J.-M.M., A.I., P.S., T.K. and F.W.; theory(DFPT):
R.H., S.K. and F.W.; resistivity: L.W.; magnetization: E.G.; data analysis: S.K., Y.S., F.W.;
paper writing: F.W., H.v.L., Y.S., S.K., J.-M.M.

**Reprints and permission** information is available online at http://npg.nature.com/
reprintsandpermissions/

**How to cite this article:** Krannich, S. *et al*. Magnetic moments induce strong phonon renormalization in FeSi. *Nat. Commun.* 6:8961 doi: 10.1038/ncomms9961 (2015).

**Additional information**

**Supplementary Information** accompanies this paper at http://www.nature.com/
naturecommunications

**Competing financial interests:** The authors declare no competing financial interests.

<figure><img src="image_6.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0

International License. The images or other third party material in this
article are included in the article's Creative Commons license, unless indicated otherwise
in the credit line; if the material is not included under the Creative Commons license,
users will need to obtain permission from the license holder to reproduce the material.
To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/