Ultrasonic elastic responses in monopole lattice
Abstract
The latest experimental advances have extended the scenario of coupling mechanical degrees of freedom in chiral magnets (MnSi/MnGe) to the topologically nontrivial skyrmion crystal and even monopole lattices. Equipped with a spinwave theory highlighting the topological features, we devise an interacting model for acoustic phonons and magnons to explain the experimental findings in a monopole lattice with a topological phase transition, i.e., annihilation of monopoleantimonopole pairs. We reproduce the anisotropic magnetoelastic modulations of elastic moduli: drastic ultrasonic softening around the phase transition and a multipeakandtrench fine structure for sound waves parallel and orthogonal to the magnetic field, respectively. Comparison with experiments indicates that the magnetoelastic coupling induced by DzyaloshinskiiMoriya interaction is comparable to that induced by the exchange interaction. Other possibilities such as elastic hardening are also predicted. The study implies that the monopole defects and their motion in MnGe play a crucial role.

March 2017
1 Introduction
Topology is now playing a more and more significant role in condensed matter physics. Apart from an upsurge in the focus on topological classification of quantum phases of matter, another field bearing ideas of topology, magnetism inhabited by stripes, vortices, domain walls, etc., has been under experimental and theoretical investigations for a long history[1]. In the recent decade, the realization of topologically nontrivial Skyrmion or Skyrmion crystal (SkX) in chiral magnets[2, 3, 4, 5, 6, 7, 8, 9, 10, 11] revives the old idea originally proposed as a hadron model[12], serving as a new scenario of the interplay between the orbital and spin of electrons and ions, and offering plenty of brand new phenomena[13, 14, 15, 16, 17, 18, 19, 20, 21, 22] together with the potential for application in magnetic storage[23, 24, 25, 26, 27, 28, 29]. Symmetry breaking of spins in noncentrosymmetric chiral magnets, which bears both the Heisenberg exchange interaction (EXI) and the DzyaloshinskiiMoriya interaction (DMI)[30, 31, 32] due to spinorbit coupling, can embody the Skyrmion texture[2, 3]. A minimal Hamiltonian for isolated Skyrmions or a Skyrmion lattice in spatial dimensions includes the EXI, the Blochtype DMI, and the Zeeman energy[33, 34] (setting throughout this paper)
(1) 
in which is the microscopic lattice constant. In twodimensional (2D) thin films, usually magnetic anisotropies can also play a key role in stabilizing the Skyrmion phase. Dissimilar to other modulated magnetic structures like the conical and helical phases that can be realized therein, a Skyrmion winds a sphere certain times and is characterized by the topological Skyrmion number[35]. In such magnetic systems, the involvement of topology is often enhanced or even induced by strong correlation effect. Indeed, based on an adiabatic approximation for the realspace Berry phase produced by the fixedlength spin texture of Skyrmion, the constraint drawn by the strong coupling with itinerant electrons can be described by the emergent electromagnetic fields (EEMF)[36, 37, 38, 15, 39].
Apart from the most common Skyrmion lattice as vertical magnetic field induced triangular lattices of Skyrmions observed in chiral magnets, the coalescence or bisection of columnar Skyrmion tubes has been observed experimentally in a bulk material[40]. Such merging points are in fact singularities or defects (the Bloch points[41, 42]) in the spin texture
where the spin moment . Specifically, they can be regarded as monopoles or antimonopoles[40, 39, 43] in terms of the EEMF and can create or annihilate Skyrmions, which are the tubes in three dimensions (3D). In a Skyrmion lattice with these defects, the Skyrmion tubes may not penetrate the sample and can go to an end at the monopole defects at a depth. For simplicity, we henceforth may use ’monopole’ to refer to both monopoles and antimonopoles. A twoSkyrmionmerging model based on the nonlinear sigma model has been used to study its effect on tranport[44]. Another regular solution for such ended Skyrmions has been investigated in the conical phase[45]. Based on the LandauLifshitzGilbert equation, people also studied the dynamics and energetics of the monopole defects[46, 47].
Partly owing to the intrinsic 3D nature of the monopole defects, there turns out to be no experimental or theoretical evidence for the formation of a lattice of the monopoles in thin films[2, 3, 48, 49, 50]. Nevertheless, it should be more feasible to be realized in 3D bulk chiral magnets. First being theoretically suggested as a superposition of multiple noncoplanar spin spirals[51], a 3D lattice of monopole defects found corroboration from the thermodynamic calculation followed[52]. Experimentally confirming these expectations, from the analysis of the electric and thermal Hall effects and the Nernst effect in contrast to those of conventional MnSilike ones, a 3D bulk MnGe under pressure was strongly inferred to possess Skyrmion tubes and a simultaneous lattice of the monopoleantimonopole defects[53, 54, 55]. Notably, with a lattice variant of the minimal Hamiltonian Eq. (1), latest Monte Carlo simulation confirmed the formation of a monopoleantimonopole lattice in 3D as long as the magnetic period becomes short[56], which is just the situation of MnGe observed in those experiments. This endorses our inclusion of the EXI and DMI in Sec. 2 among other possible magnetic interactions. Further, a recent study using Lorentz transmission electron microscopy clearly revealed the magnetic structure of MnGe in concordance with the threeorthogonalspiral model anticipated[57].
It is necessary to pay attention to the distinction between two points of view, i.e., the spin texture and the original spin moment field , of the magnetic ordering background, upon which we shall mainly study the effect due to fluctuations. The latter, , which is for constructing the multispiral spin density waves, is nonsingular and topologically trivial because any configuration mapped onto a 3ball can be smoothly connected to . However, this work aims at highlighting the effect characteristic for the topologically nontrivial defects and the orientational field , which is inherently described by the second homotopy group of a 2sphere that cannot shrink to a point, naturally appears to be physically more relevant. On the other hand, for MnGe, in the strong correlation regime, localized spins’ moments can hardly vary in magnitude and the saturated magnetization is large and only vanishes at the singular defects. Indeed, this choice of the unitlength field already proved to be appropriate for the description of the influence on conduction electron in the strongly correlated MnGe material, where the 3D monopole lattice even exists in the absence of external magnetic field[43, 39].
Expecting new effects due to the monopole defects, we initiated our theoretical investigation to study the phenomena emerging from the coupling between other components and the topologically nontrivial monopole lattice in the hosting material MnGe by calculating the influence from spinwave excitations therein. Our first study focuses on the longitudinal electric transport in MnGe massively influenced by the quantum fluctuations of the EEMF due to spin waves and identified a topological phase transition of strong correlation genesis[39]. One more intriguing possibility is the coupling to phonons in the hosting solids. Indeed, unconventional ultrasonic responses of elastic stiffness have been observed recently. Experimentalists employed longitudinal sound waves to the Skyrmion phases in both MnSi[58, 59] and MnGe[43] with the sound propagating direction parallel (mode) or orthogonal (mode) to the applied external magnetic field along the axis as shown in Fig. 1. They recognized the SkX phase by a distinct elastic stiffness with anisotropy for the two modes in MnSi. Dissimilarly, varying with external magnetic field, MnGe showed not only a drastic softening of mode stiffness but also a multipeakandtrench fine structure of mode stiffness, which are much stronger than those in the SkX phase of MnSi. The range of the magnetic field for the softening considerably coincides with the topological phase transition aforementioned, which is the destruction of the monopole lattice. Fluctuations associated with such qualitative structural transition is therefore not surprising to produce dramatic modification in mechanical properties. The EEMF not only helps describe itinerant electrons coupled with localized spins but also captures the essential feature of the nontrivial spin texture, hence we find it rather informative even for the magnetoelastic response to be discussed. And it is reasonable to owe the qualitatively different responses of MnGe compared with MnSi to the influence from the monopoles/antimonopoles on the magnon fluctuations therein.
Towards this end, we reuse our spinwave theory and propose in this paper an interaction theory for magnons and phonons induced by both EXI and DMI, which reproduces and explains the magnetoelastic experimental features. We refer the reader to the separate papers[39, 43] for the detailed description of the monopole lattice consisting of three orthogonal spin spirals, the spinwave theory of a Skyrmion lattice or monopole lattice, and the experimental confirmation. In Sec. 2, we first devise the magnetoelastic interactions and derive the effective phonon theory by integrating out the magnon degrees of freedom. The influence of the magnons on the renormalized phonon excitations is discussed in Sec. 3.1 with an emphasis on the renormalized phonon linear mode. In Sec. 3.2 and Sec. 3.3, we compare the experimental observations with our theoretical predictions and explain the physical origin of the magneticfielddependent evolution of the hybridized excitation spectra. We emphasize the distinction between the monopole lattice in MnGe and the conventional SkX in MnSi in the rest of Sec. 3 and conclude in Sec. 4.
2 Theoretical Models
2.1 Magnetoelastic coupling
As aforesaid, it is the unitnorm constrained spin texture that produces the EEMF and yields the nontrivial topology of the Skyrmion and the monopole rather than the bare spin moment itself, which manifests the strong electron correlation therein[39, 43]. If were used, the form factors to be calculated in the following would assume too simple forms to give rise to the complex enough magneticfield dependence and the correct theory to reproduce experiments. We thus make use of instead of to derive the appropriate magnetoelastic interactions from the Hamiltonian Eq. (1). This is the minimal model capable of producing a stable monopole lattice as discussed previously. And this study shows that EXI and DMI are adequate and essential to capture the correct physics herein.
When longitudinal sound waves are artificially generated in the material, the lattice structure will be perturbed, which is described by the longitudinal phonon. However, the strengths of EXI and DMI depend on lattice bond lengths hence the magnetoelastic coupling (Fig. 2). Such a perturbation to the magnetic interactions by mechanical degrees of freedom can be accounted for by expanding the strengths of EXI and DMI up to the linear order of phonon degrees of freedom
(2a)  
(2b) 
wherein is the lattice displacement from equilibrium in the continuum limit, is the coefficient of this expansion, and summation over repeated indices is understood henceforth. We emphasize four features of this expansion. 1) In conformity with translational symmetry, either or depends on rather than lattice displacement itself. 2) Alongside, obviously deals with the longitudinal phonon. 3) All spatial derivatives match in the direction since any , reflecting magnetic interactions between two adjacent spins sitting at the endpoints of a lattice bond along the direction in a lattice model, should be affected by the longitudinal phonon propagating along the same direction through . 4) As stated in 3), these magnetoelastic couplings take anisotropy into account so as to produce the delicate experimental observations.
The next step is to expand as in Eqs. (2a)(2b) to incorporate the spinwave degrees of freedom and the superscript means the ground state value or setting after taking the derivative. For a spin spiral indexed by , this is the phase of the constituent spin density wave, which also indicates the spatial shift of the Skyrmion lattice or the position of the monopole defects. And is the magnetization along the propagation axis of spin spiral , generating rotation around this axis. Note that we consider the situation where there are three orthogonal spin spirals and it is those quantities’ deviations away from their static mean field values that constitute the spinwave fields in the expansion[39]. Without loss of generality, we thus set the static value of any to and denote the fluctuation in by .
Starting from the EXI and DMI induced magnetoelastic interaction parts in Eqs. (2a)(2b), up to terms bilinear in spinwave and phonon fields, one arrives at
(2ca)  
(2cb) 
wherein we temporarily omit all the prefactors for simplicity. One crucial criterion is the translational invariance in the continuum model to be derived, i.e., for the spinwave field , only its spatial derivative can enter simply because is the (phase) displacement field of spin spiral . Therefore, when , i.e., for the field, applied upon the parenthesis in the first term in either Eq. (2ca) or Eq. (2cb) is equivalent to a spatial derivative. And we notice that the typical wave vector of the periodic function in the parenthesis is much larger than that of the slowly varying spinwave field . Since Eqs. (2ca)(2cb) will be part of a Lagrangian density, a spatial integration can be carried assuming is uniform, resulting in zero in the end for the total derivative. This procedure exactly eliminates the term linear in field, otherwise the translational invariance would be violated. However, this is unnecessary for , i.e., for the fluctuating field . Therefore, in either Eq. (2ca) or Eq. (2cb), for , only the second term remains.
Thus for either EXI or DMI case, we attain two terms
(2cd) 
in total, wherein but whenever . We call the functions preceding the bilinear fields form factors, which render the coupling between magnon and phonon nonuniform in space. They are defined as and actually record the information of the complicated magnetic structure affecting the EXI or DMI induced magnetoelastic interactions. After Fourier transform to the momentum space, we take the spatial average over a magnetic unit cell of the form factors, making the couplings dependent only on the variable uniform magnetization along the axis. Therefore we have an exactly solvable theory without couplings between unequal momenta. This uniform magnetization directly relates to the external magnetic field applied on the system along the axis. The magnetic field dependence of the ultrasonic responses is due to the fact that the magnetic texture of the monopole lattice varies with respect to [39].
2.2 Effective theory of phonon
A standard free theory of the longitudinal phonon reads
(2ce) 
wherein is an elastic constant, is the mass density and we work in imaginary time henceforth. For MnGe, spinwave takes the form[39]
(2cf) 
wherein . Note that the dependent is the emergent magnetic field that characterizes the nontrivial topology of the Skyrmions or monopoles in the system. And we get the magnetoelastic interactions
(2cg) 
from Sec. 2.1. Because the three parts, Eqs. (2ce)(2cf)(2cg), comprising the full theory are all bilinaear, we can diagonalize the actions in the energymomentum space. Below, is used as a shorthand for the argument , wherein is a generic complex frequency that can equal for instance. The longitudinal phonon action is transformed to
in which matrix and . The spinwave action is transformed to
in which matrix can be solved and . The action of magnetoelastic interactions becomes
wherein vector , matrices .
Now the full theory is constituted of three parts , which are the free phonon theory, the free magnon theory, and their coupling due to the magnetoelastic interactions in turn. Next, we can integrate out the harmonic magnon bath in the pathintegral formalism, leading to the partition function expressed as
Because of the bilinearity, after a Gaussian functional integration, we get the effective action for the phonons
(2ch) 
The expression in the bracket is nothing but multiplying the inverse of the Matsubara Green’s function matrix of the effective phonon theory in the energymomentum space.
3 Methods and results
3.1 Renormalized phonon spectrum
Instead of the hardly accessible analytic dispersion relations, we resort to inspecting the spectral function of this effective theory. In conformity with experimental investigations, we focus on the effects due to sound waves propagating along and directions by looking at the diagonal as a function of momentum along direction and frequency , wherein retarded Green’s function comes from the analytic continuation . First of all, we confirm the property of the bosonic spectral function[60] that it is always positive for and negative for . Secondly, the renormalized phonon spectrum is reflected in plot by functiontype ridge structures. Thirdly, setting a realistic but small enough , we can extract the information of the phonon excitations from the plot by identifying the ridge structures.
Equation (2ce) itself can only give banal longitudinal phonon excitations of linear dispersion with the acoustic velocity . On the other hand, the magnon theory possesses an excitation spectrum composed of three distinct modes[39]. For MnGe, there are two kinds of gapless excitation, one acoustic mode and one quadratic mode when is small, and another excitation with an energy gap , as shown in Fig. 3. Because of the spin Berry phase of the topologically nontrivial spin textures of Skyrmions/monopoles, any two unparallel spin spirals are possible to mingle with each other. This gives rise to quadratic couplings like in Eq. (2cf), mixing up conjugate fields and hence the quadratic mode[39, 15, 61]. An exception is that these three modes degenerate into the linear gapless mode when . These gapless excitations are the NambuGoldstone bosons due to the broken translation symmetry. In the case of electric transport[39], the softest one affects the electrons’ motion the most. However, all the gapless modes will play a crucial role in the lowenergy physics describing the interplay with the linear acoustic phonons.
The salient point is that because of the magnetoelastic coupling, the three magnon modes will hybridize with the phonon mode, giving rise to rich possibilities of renormalized excitation spectra. One interesting instance occurs when the phonon mode intersects with the magnon modes in the dispersion plot, resulting in mutual repulsion and reconstruction of the involved dispersion curves. In the gross, despite possible reconstructions, one is in general still able to relate the new modes to their respective precedents before hybridization, which will henceforth be used as convenient tags of the new modes in the effective phonon theory. Even without any intersection in the dispersion curves, the intensity of a new magnon mode (sharpness of the ridge structure) in a spectral function plot directly reflects the degree of hybridization that influences the (approximately) linear renormalized phonon mode. We have the relation (seen from Eq. (2ce)) between the velocity of the renormalized phonon linear mode and the new stiffness . This phonon linear mode after hybridization is of the major importance because it is this stiffness that is measured experimentally at different external magnetic fields as the ultrasonic responses. Its velocity is the slope of the corresponding dispersion ridge extracted from scanning the spectral function, e.g., in the shaded region in Fig. 4.
Because of the experimental difficulty in determining the accurate value of , and the magnetoelastic couplings , one has to search for the correct control parameters corresponding to the real materials. We consider two, and , for the effective theory. One is the ratio between strengths of different magnetoelastic couplings , denoted as . The other is the ratio of the velocity of the magnon linear mode, , to the original unperturbed acoustic velocity , denoted as . We also give an estimate of the characteristic energy scales in the MnGe experiment. The sound wave propagating in MnGe crystal has frequency and wavelength [43]. The original acoustic velocity is hence and the phonon energy is . As for the material MnGe[54, 57], is about 4Å and the melting temperature () of the magnetic order can be used to estimate the strength of EXI . We set for MnGe. Thus, it is straightforward to obtain and , using the typical value when . On the other hand, in our theoretical study, we set , , and hence . Consequently, or is fixed by . In addition to the foregoing variable , we set since a too large magnetoelastic coupling unrealistically alters the phonon spectrum while a tiny one renders the effect feeble to detect.
The first message from the above is that is of the same order as . This implies that we had better tune not far from unity if we were to explain the experiment. The second message is , which means the lowenergy phonon excitations and hence gapless magnons around the long wavelength limit play a major role. Because the matrix in Eq. (2ch) contains two individual parts due to EXI and DMI, the combination becomes not simply a summation of the separate effects, but one with inevitable interference between the two types of magnetoelastic interactions. Indeed, we have seen distinct profiles when changing within a typical range by step . Note that the sign difference between two magnetoelastic interactions is possible as implied by the sign change in DMI[62, 63, 64].
3.2 Rich possiblities of magnetoelastic responses
The value of strongly affects the excitation spectrum of the effective phonon theory in a clearer manner as compared with . Imagine drawing the original unperturbed phonon dispersion line in Fig. 3, whose slope might be smaller or larger than that of the blue magnon linear mode, providing a natural bipartite classification. 1) . The phonon mode lies below the blue magnon linear mode and intersects with the green magnon quadratic mode. On one hand, if both the magnon quadratic mode itself and the magnetoelastic interactions are strong enough in intensity, depending on the details of the coupling, the intersection becomes an anticrossing of two reconstructed modes repelling each other in the effective phonon spectrum. On the other hand, near the structural phase transition where fluctuations become immense, the blue magnon mode strongly repels the phonon mode downward as a result of the magnonphonon interaction, serving as the major cause of the softening effect in Fig. 5 in agreement with the experimental data shown Fig. 6. 2) . The phonon mode lies above the blue magnon linear mode and intersects with the orange magnon gapped mode. In spite of this, the high energy scale of the gap renders itself irrelevant for the long wavelength phonons. Similar to 1), near the phase transition, the blue mode repels the phonon mode upward, which is an elastic hardening prediction from our theory. As seen in Fig. 7, the renormalized phonon velocity is in general larger than its original unity value. It is expected to be realized by changing either or . Recent studies of controllable DMI in Skyrmion lattices by exerting strains[62, 65] or varying compositions[64] can be candidates for realizing this case. Last but not least, when , the foregoing distinction becomes vague. As a minor reassurance, we indeed saw evident, albeit complex transitions from softening to hardening in when traverses the range .
We discuss more details about the comparison between theory and experiment. Figure 5 shows the theoretical result that reproduces the experimentally observed signals in Fig. 6, including the drastic decrease in the velocity of the mode and the multipeakandtrench fine structure for the mode. In this model calculation, we cannot produce a temperature dependence simply because the spectral function in this case doesn’t depend on temperature. However, as suggested by basically the same feature at different temperatures marked by the gray dots in Fig. 6, the current theory is able to illustrate the essential physics therein. The external magnetic field in Fig. 6 is in general not simply proportional to the parameter in Fig. 5 since should be regarded as an approximation in modeling the effect on the deformation of the magnetic structure due to the possibly complex magnetization process. Nevertheless, this comparison suffices to highlight the key magnetoelastic responses. The reliability of this result is supported by the fact that within the range of the control parameters, the basic characters hold all along. After all, we can notice that in Fig. 5 the most evident changes in the velocity (stiffness) always reside around the destruction of the monopole lattice, i.e., the monopoleantimonopole pair annihilation at [39], and especially this topological phase transition clearly manifests itself by the drastic softening in Fig. 5(b). As for the magnitude of the softening with respect to the original situation, our theory produces , which is a bit smaller than the experimental value between 6 and 10. Despite this discrepancy, we highlight the excellent consistency between the theoretical and experimental fine structures. Besides the clear match for singletrench mode case of Fig. 5(b), all three trenches and two peaks in Fig. 5(a) find their counterparts in the experiment. On the other hand, if we set or to zero, it becomes impossible to reproduce the experimental signals. Taken as a whole, these strongly endorse our theory and can be regarded as a new way to find out some quantities temporarily outside the reach of experimental detection.
Another aspect worth discussing is why the experimental features are reproduced when is not far away from unity although we have . The resolution consists in DMI’s unique sensitivity to minute strains[62, 63, 65]. Some anticrossing points mixing up spinup and spindown bands in the band structure contribute massively to the spinorbit interaction hence the DMI. The phonon induced strain modifies the band structure slightly, which, however, may considerably change the DMI because the position of the Fermi level in regard to the nearby anticrossing points can drastically change the contribution from these points. Experimentally, typically ten times larger relative modulation of DMI than that of EXI is observed[65], which is consistent with our theoretical finding.
3.3 Magneticfielddependent evolution of the hybridized excitation spectra
Obviously, noticing the range of , Fig. 5 belongs to the foregoing case. A careful inspection of the dispersion profiles, i.e., spectral function plots Fig. 8 (mode) and Fig. 9 (mode), leads us to the following explanation.
For the mode, the dominant modes appear to be the linear excitations of magnons (upper branch) and phonons (lower branch). Increasing , the slope of the magnon mode almost remains the same while its intensity gradually increases () to some maximum value around the phase transition and drops down afterward (). The larger the intensity becomes, the stronger repulsion is exerted to the phonon mode underneath, explaining the drastic softening. The reason why this case lacks the participation of the magnon quadratic mode lies in the magnetoelastic interaction in Eq. (2cd). The magnon quadratic mode originates from the quadratic term in our spinwave theory in Eq. (2cf). Thus, if it were to largely affect the phonon linear mode in the effective phonon theory, field must have an adequate coupling with phonons. The relevant form factors hereof, when , is diagonal as seen from the calculation, making the mode (mode) phonon primarily couples with (). But only and in the quadratic term are important since and are vanishing.
For the mode, the interplay takes place mainly between the original phonon linear mode and the magnon linear and quadratic modes as seen from their high intensities in Fig. 9. We can first faintly recognize the three magnon modes (gapped, linear and quadratic, from left to right) together with the always dominant phonon linear mode when , which is in contrast to the case with only one magnon linear dispersion present. This is because of the degeneracy of magnon excitations when mentioned in Sec. 3.1. Note that the gapped magnon mode disappears in most plots since increasing enlarges the gap beyond the plot range. For and , a typical (anti)crossing or reconstruction of the magnon and phonon linear modes occurs at some lowenergy scale inside the plot range. The mutual repulsion between them makes the reconstructed phonon linear mode, especially the part right to the crossing point, move downward, explaining the first small trench in Fig. 5(a). Once the crossing point gets higher (), the lowenergy part of the phonon mode simply bounces back and gives rise to the first small peak near in Fig. 5(a). Thereafter, as increases till the vicinity of , the magnon linear mode above the phonon linear mode becomes stronger and stronger while the magnon quadratic mode keeps moving down till , which constructively pushes the phonon mode downward, yielding the deepest trench. And then, abruptly, a transient reverse procedure is observed approximately from to plots, hence the dramatic upsurge in Fig. 5(a). Next, the magnon quadratic mode recovers all the way back and crosses over the phonon mode at lower and lower energy scales while the magnon linear mode becomes stronger till , which is again a constructive effect of dragging downward the phonon mode. Finally, a fading reunion of the three magnon modes are observed from to , which is natural for an induced ferromagnetism of too large .
In a nutshell, the experimentally observed magnetoelastic phenomena are the consequences of two aspects that vary with or the external magnetic field. One is the magnon quadratic mode generated at nonzero (together with the herein insignificant gapped mode), whose reciprocating shift in the spectral function plot in Fig. 9 is clearly controlled by the emergent magnetic field as explained in Sec. 3.1 in terms of Eq. (2cf). The other is the family of form factors directly affecting the intensities of and hybridization between various modes. They show nonmonotonous behavior upon increasing and vary strongly near the phase transition. An intricate integration of these two aspects leads to the rich experimental features. A detailed inspection of the underlying ground state spin configuration provides more insights[39]. Here we only recapitulate two key aspects. One is the nonmonotonic profile of , whose maximum near and fast dip around make the form factors’ variation more perceivable. Note that the realistic lattice cutoff introduced to the monopole defects in the calculation will also postpone the complete destruction of the monopole lattice to some value larger than the ideal value . The other, from a more intriguing viewpoint, is the nontrivial contribution from the (anti)monopoles and their characteristic collisionandannihilation motion. In contrast to the SkX in MnSi discussed below, we have argued for the crucial role of the topological defects, i.e, the (anti)monopoles, in the magnetoresistivity and especially the topological phase transition of the destruction of the monopole lattice. Here in the magnetoelastic phenomena, not only but also the form factors are reflecting the rich facets of the spin textures. It is the coupling between phonons and the (anti)monopoles that causes all the complexities.
3.4 Relation to triangular lattice of Skyrmion tubes in MnSi
Now we briefly discuss the magnetoelastic couplings in the SkX of MnSi. Except from the modification in the spinwave theory, we can apply the similar form of couplings to MnSi. Nevertheless, we observe vanishing couplings when the sound wave propagates along the cylindrical symmetry axis of Skyrmion tubes in MnSi. For the EXI induced coupling, the form factors are simply zero due to the translational symmetry along the direction. Although this is not the case for the coupling induced by DMI, the spatial averages of corresponding form factors turn out to be zero in the end. As for the perpendicular propagating case, we observe nonvanishing magnetoelastic effects, albeit negligibly smaller as compared to MnGe. On the other hand, the experimental signals () of the SkX in MnSi are quite smaller than those () detected in MnGe[43]. Therefore, our form of magnetoelastic couplings turns out to be the leading order effect in the MnGe monopole lattice. Based on the solid state mechanics and macroscopic thermodynamics, many more complex higher order terms and fitting parameters are incorporated into a GinzburgLandau free energy calculation[66] for the ultrasonic signals in MnSi[58, 59]. This comparison actually lends credence to the aforementioned essential role played by the monopole defects that are absent in MnSi. In other words, the contribution from monopole defects to the magnetoelastic effects, if present at all, dominates and is captured by the formalism developed in this work.
4 Concluding remark
The study of ultrasonic elastic responses possesses strong motivation from the experimental findings. We not only explain the observed softening effect but also predict new issues, e.g., hardening and other patterns of the dependence of the stiffness on the magnetic field, which in return provides a way to determine some experimentally inaccessible physical quantities. Based on a wellestablished spinwave theory from the previous magnetoresistivity study, we are able to identify once again the nontrivial features of the monopole lattice, especially the topological phase transition signifying strong correlations. Thanks to the agreement with the experimental observations, this magnetoelasticity study, together with the magnetoresistivity one, can establish the importance of the topological nature of the spin configuration in strongly correlated electronic systems. As a whole, they speak for a crucial role played by the monopole defects in chiral magnet MnGe. In particular, these studies pave the way for even more intriguing scenarios of coupling topologically nontrivial objects with other systems. They show the rich physics therein and help us gain insights for further investigations towards plenty of manipulation methods for Skyrmionics applications.
Acknowledgments
We thank Yoichi Nii for useful discussions and the indispensable experimental results. X.X.Z was partially supported by the Panasonic Scholarship and by Japan Society for the Promotion of Science through Program for Leading Graduate Schools (ALPS) and GrantinAid for JSPS Fellows (No. 16J07545). This work was supported by JSPS GrantinAid for Scientific Research (No. 24224009) and JSPS GrantinAid for Scientific Research on Innovative Areas (No. 26103006) from MEXT, Japan, and ImPACT Program of Council for Science, Technology and Innovation (Cabinet office, Government of Japan). This work was also supported by CREST, Japan Science and Technology Agency.
References
 A. P. Malozemoff. Magnetic Domain Walls in Bubble Materials: Advances in Materials and Device Research. Academic Press, 1 1979.
 A. N. Bogdanov and D. A. Yablonskii. Thermodynamically stable ”vortices” in magnetically ordered crystals. the mixed state of magnets. JETP, 95:182, 1989.
 U. K. Rößler, A. N. Bogdanov, and C. Pfleiderer. Spontaneous skyrmion ground states in magnetic metals. Nature, 442(7104):797–801, Aug. 2006.
 S. Mühlbauer, B. Binz, F. Jonietz, C. Pfleiderer, A. Rosch, A. Neubauer, R. Georgii, and P. Böni. Skyrmion lattice in a chiral magnet. Science, 323(5916):915–919, 2009.
 Akira Tonomura, Xiuzhen Yu, Keiichi Yanagisawa, Tsuyoshi Matsuda, Yoshinori Onose, Naoya Kanazawa, Hyun Soon Park, and Yoshinori Tokura. Realspace observation of skyrmion lattice in helimagnet thin samples. Nano Letters, 12(3):1673–1677, 2012. PMID: 22360155.
 X. Z. Yu, N. Kanazawa, Y. Onose, K. Kimoto, W. Z. Zhang, S. Ishiwata, Y. Matsui, and Y. Tokura. Near roomtemperature formation of a skyrmion crystal in thinfilms of the helimagnet FeGe. Nature Materials, 10(2):106–109, Dec. 2010.
 X. Z. Yu, Y. Onose, N. Kanazawa, J. H. Park, J. H. Han, Y. Matsui, N. Nagaosa, and Y. Tokura. Realspace observation of a twodimensional skyrmion crystal. Nature, 465(7300):901–904, Jun. 2010.
 W. Münzer, A. Neubauer, T. Adams, S. Mühlbauer, C. Franz, F. Jonietz, R. Georgii, P. Böni, B. Pedersen, M. Schmidt, A. Rosch, and C. Pfleiderer. Skyrmion lattice in the doped semiconductor . Phys. Rev. B, 81:041203, Jan. 2010.
 S. Seki, X. Z. Yu, S. Ishiwata, and Y. Tokura. Observation of skyrmions in a multiferroic material. Science, 336(6078):198–201, 2012.
 Stefan Heinze, Kirsten von Bergmann, Matthias Menzel, Jens Brede, André Kubetzka, Roland Wiesendanger, Gustav Bihlmayer, and Stefan Blügel. Spontaneous atomicscale magnetic skyrmion lattice in two dimensions. Nature Physics, 7(9):713–718, 2011.
 Stefan Heinze, Kirsten von Bergmann, Matthias Menzel, Jens Brede, André Kubetzka, Roland Wiesendanger, Gustav Bihlmayer, and Stefan Blügel. Spontaneous atomicscale magnetic skyrmion lattice in two dimensions. Nature Physics, 7(9):713–718, Jul. 2011.
 T.H.R. Skyrme. A unified field theory of mesons and baryons. Nuclear Physics, 31:556 – 569, 1962.
 Minhyea Lee, W. Kang, Y. Onose, Y. Tokura, and N. P. Ong. Unusual hall effect anomaly in under pressure. Phys. Rev. Lett., 102:186601, May 2009.
 A. Neubauer, C. Pfleiderer, B. Binz, A. Rosch, R. Ritz, P. G. Niklowitz, and P. Böni. Topological hall effect in the phase of . Phys. Rev. Lett., 102:186602, May 2009.
 Jiadong Zang, Maxim Mostovoy, Jung Hoon Han, and Naoto Nagaosa. Dynamics of skyrmion crystals in metallic thin films. Phys. Rev. Lett., 107:136804, Sep. 2011.
 Karin Everschor, Markus Garst, Benedikt Binz, Florian Jonietz, Sebastian Mühlbauer, Christian Pfleiderer, and Achim Rosch. Rotating skyrmion lattices by spin torques and field or temperature gradients. Phys. Rev. B, 86:054432, Aug. 2012.
 R. Ritz, M. Halder, M. Wagner, C. Franz, A. Bauer, and C. Pfleiderer. Formation of a topological nonfermi liquid in . Nature, 497(7448):231–234, May 2013.
 Wei Li, Chiming Jin, Renchao Che, Wensen Wei, Langsheng Lin, Lei Zhang, Haifeng Du, Mingliang Tian, and Jiadong Zang. Emergence of skyrmions from rich parent phases in the molybdenum nitrides. Phys. Rev. B, 93:060409, Feb 2016.
 H Velkov, O Gomonay, M Beens, G Schwiete, A Brataas, J Sinova, and R A Duine. Phenomenology of currentinduced skyrmion motion in antiferromagnets. New Journal of Physics, 18(7):075016, 2016.
 F. Jonietz, S. Mühlbauer, C. Pfleiderer, A. Neubauer, W. Münzer, A. Bauer, T. Adams, R. Georgii, P. Böni, R. A. Duine, K. Everschor, M. Garst, and A. Rosch. Spin transfer torques in at ultralow current densities. Science, 330(6011):1648–1651, 2010.
 X.Z. Yu, N. Kanazawa, W.Z. Zhang, T. Nagai, T. Hara, K. Kimoto, Y. Matsui, Y. Onose, and Y. Tokura. Skyrmion flow near room temperature in an ultralow current density. Nature Communications, 3:988, Aug. 2012.
 Lingyao Kong and Jiadong Zang. Dynamics of an insulating skyrmion under a temperature gradient. Phys. Rev. Lett., 111:067203, Aug 2013.
 N S Kiselev, A N Bogdanov, R SchÃ¤fer, and U K Rößler. Chiral skyrmions in thin magnetic films: new objects for magnetic storage technologies? Journal of Physics D: Applied Physics, 44(39):392001, 2011.
 Albert Fert, Vincent Cros, and João Sampaio. Skyrmions on the track. Nature Nanotechnology, 8(3):152–156, Mar. 2013.
 Niklas Romming, Christian Hanneken, Matthias Menzel, Jessica E. Bickel, Boris Wolter, Kirsten von Bergmann, André Kubetzka, and Roland Wiesendanger. Writing and deleting single magnetic skyrmions. Science, 341(6146):636–639, 2013.
 R. Tomasello, E. Martinez, R. Zivieri, L. Torres, M. Carpentieri, and G. Finocchio. A strategy for the design of skyrmion racetrack memories. Sci. Rep., 4:6784, Oct. 2014.
 Wataru Koshibae, Yoshio Kaneko, Junichi Iwasaki, Masashi Kawasaki, Yoshinori Tokura, and Naoto Nagaosa. Memory functions of magnetic skyrmions. Japanese Journal of Applied Physics, 54(5):053001, 2015.
 Xichao Zhang, Motohiko Ezawa, and Yan Zhou. Magnetic skyrmion logic gates: conversion, duplication and merging of skyrmions. Sci. Rep., 5:9400, Mar. 2015.
 Christian Hanneken, AndrÃ© Kubetzka, Kirsten von Bergmann, and Roland Wiesendanger. Pinning and movement of individual nanoscale magnetic skyrmions via defects. New Journal of Physics, 18(5):055009, 2016.
 I. Dzyaloshinskii. A thermodynamic theory of âweakâ ferromagnetism of antiferromagnetics. Journal of Physics and Chemistry of Solids, 4(4):241 – 255, 1958.
 Tôru Moriya. Anisotropic superexchange interaction and weak ferromagnetism. Phys. Rev., 120:91–98, Oct 1960.
 A. Fert and Peter M. Levy. Role of anisotropic exchange interactions in determining the properties of spinglasses. Phys. Rev. Lett., 44:1538–1541, Jun 1980.
 P Bak and M H Jensen. Theory of helical magnetic structures and phase transitions in mnsi and fege. Journal of Physics C: Solid State Physics, 13(31):L881, 1980.
 S. V. Maleyev. Cubic magnets with dzyaloshinskiimoriya interaction at low temperature. Phys. Rev. B, 73:174402, May 2006.
 R. Rajaraman. Solitons and Instantons, Volume 15: An Introduction to Solitons and Instantons in Quantum Field Theory (NorthHolland Personal Library). North Holland, 1 edition, 4 1987.
 Naoto Nagaosa and Yoshinori Tokura. Emergent electromagnetism in solids. Physica Scripta, 2012(T146):014020, Jan. 2012.
 T. Schulz, R. Ritz, A. Bauer, M. Halder, M. Wagner, C. Franz, C. Pfleiderer, K. Everschor, M. Garst, and A. Rosch. Emergent electrodynamics of skyrmions in a chiral magnet. Nature Physics, 8(4):301–304, Feb. 2012.
 N. Nagaosa, X. Z. Yu, and Y. Tokura. Gauge fields in real and momentum spaces in magnets: monopoles and skyrmions. Philosophical Transactions of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, 370(1981):5806–5819, Nov. 2012.
 XiaoXiao Zhang, Andrey S. Mishchenko, Giulio De Filippis, and Naoto Nagaosa. Electric transport in threedimensional skyrmion/monopole crystal. Phys. Rev. B, 94:174428, Nov 2016.
 P. Milde, D. Köhler, J. Seidel, L. M. Eng, A. Bauer, A. Chacon, J. Kindervater, S. Mühlbauer, C. Pfleiderer, S. Buhrandt, C. Schütte, and A. Rosch. Unwinding of a skyrmion lattice by magnetic monopoles. Science, 340(6136):1076–1080, 2013.
 G E Volovik. Linear momentum in ferromagnets. Journal of Physics C: Solid State Physics, 20(7):L83, 1987.
 P. R. Kotiuga. The algebraic topology of bloch points. IEEE Transactions on Magnetics, 25(5):3476–3478, Sep 1989.
 N. Kanazawa, Y. Nii, X.X. Zhang, A. S. Mishchenko, G. De Filippis, F. Kagawa, Y. Iwasa, N. Nagaosa, and Y. Tokura. Critical phenomena of emergent magnetic monopoles in a chiral magnet. Nature Communications, 7:11622, may 2016.
 Rina Takashima and Satoshi Fujimoto. Electrodynamics in skyrmions merging. Journal of the Physical Society of Japan, 83(5):054717, 2014.
 Filipp N. Rybakov, Aleksandr B. Borisov, Stefan Blügel, and Nikolai S. Kiselev. New type of stable particlelike states in chiral magnets. Phys. Rev. Lett., 115:117201, Sep 2015.
 Christoph Schütte and Achim Rosch. Dynamics and energetics of emergent magnetic monopoles in chiral magnets. Phys. Rev. B, 90:174432, Nov. 2014.
 ShiZeng Lin and Avadh Saxena. Dynamics of dirac strings and monopolelike excitations in chiral magnets under a current drive. Phys. Rev. B, 93:060401, Feb 2016.
 M. N. Wilson, A. B. Butenko, A. N. Bogdanov, and T. L. Monchesky. Chiral skyrmions in cubic helimagnet films: The role of uniaxial anisotropy. Phys. Rev. B, 89:094411, Mar 2014.
 A. B. Butenko, A. A. Leonov, U. K. Rößler, and A. N. Bogdanov. Stabilization of skyrmion textures by uniaxial distortions in noncentrosymmetric cubic helimagnets. Phys. Rev. B, 82:052403, Aug 2010.
 H Wilhelm, M Baenitz, M Schmidt, C Naylor, R Lortz, U K Rößler, A A Leonov, and A N Bogdanov. Confinement of chiral magnetic modulations in the precursor region of fege. Journal of Physics: Condensed Matter, 24(29):294204, 2012.
 B. Binz and A. Vishwanath. Theory of helical spin crystals: Phases, textures, and properties. Phys. Rev. B, 74:214408, Dec. 2006.
 JinHong Park and Jung Hoon Han. Zerotemperature phases for chiral magnets in three dimensions. Phys. Rev. B, 83:184406, May 2011.
 N. Kanazawa, Y. Onose, T. Arima, D. Okuyama, K. Ohoyama, S. Wakimoto, K. Kakurai, S. Ishiwata, and Y. Tokura. Large topological hall effect in a shortperiod helimagnet . Phys. Rev. Lett., 106:156603, Apr. 2011.
 N. Kanazawa, J.H. Kim, D. S. Inosov, J. S. White, N. Egetenmeyer, J. L. Gavilano, S. Ishiwata, Y. Onose, T. Arima, B. Keimer, and Y. Tokura. Possible skyrmionlattice ground state in the B chirallattice magnet as seen via smallangle neutron scattering. Phys. Rev. B, 86:134425, Oct. 2012.
 Y. Shiomi, N. Kanazawa, K. Shibata, Y. Onose, and Y. Tokura. Topological nernst effect in a threedimensional skyrmionlattice phase. Phys. Rev. B, 88:064409, Aug 2013.
 SeongGyu Yang, YeHua Liu, and Jung Hoon Han. Formation of a topological monopole lattice and its dynamics in threedimensional chiral magnets. Phys. Rev. B, 94:054420, Aug 2016.
 Toshiaki Tanigaki, Kiyou Shibata, Naoya Kanazawa, Xiuzhen Yu, Yoshinori Onose, Hyun Soon Park, Daisuke Shindo, and Yoshinori Tokura. Realspace observation of shortperiod cubic lattice of skyrmions in . Nano Letters, 15(8):5438–5442, 2015. PMID: 26237493.
 Y. Nii, A. Kikkawa, Y. Taguchi, Y. Tokura, and Y. Iwasa. Elastic stiffness of a skyrmion crystal. Phys. Rev. Lett., 113:267203, Dec. 2014.
 A. E. Petrova and S. M. Stishov. Field evolution of the magnetic phase transition in the helical magnet mnsi inferred from ultrasound studies. Phys. Rev. B, 91:214402, Jun 2015.
 Gerald D. Mahan. ManyParticle Physics (Physics of Solids and Liquids). Springer, 3rd ed. 2000 edition, 10 2000.
 Olga Petrova and Oleg Tchernyshyov. Spin waves in a skyrmion crystal. Phys. Rev. B, 84:214433, Dec 2011.
 Takashi Koretsune, Naoto Nagaosa, and Ryotaro Arita. Control of dzyaloshinskiimoriya interaction in : a firstprinciples study. Sci. Rep., 5:13302, Aug. 2015.
 J. Gayles, F. Freimuth, T. Schena, G. Lani, P. Mavropoulos, R. A. Duine, S. Blügel, J. Sinova, and Y. Mokrousov. Dzyaloshinskiimoriya interaction and hall effects in the skyrmion phase of . Phys. Rev. Lett., 115:036602, Jul. 2015.
 K. Shibata, X. Z. Yu, T. Hara, D. Morikawa, N. Kanazawa, K. Kimoto, S. Ishiwata, Y. Matsui, and Y. Tokura. Towards control of the size and helicity of skyrmions in helimagnetic alloys by spin–orbit coupling. Nature Nanotechnology, 8(10):723–728, Sep. 2013.
 K. Shibata, J. Iwasaki, N. Kanazawa, S. Aizawa, T. Tanigaki, M. Shirai, T. Nakajima, M. Kubota, M. Kawasaki, H. S. Park, D. Shindo, N. Nagaosa, and Y. Tokura. Large anisotropic deformation of skyrmions in strained crystal. Nature Nanotechnology, 10(7):589–592, Jun. 2015.
 Y. Hu and B. Wang. Unified Theory of Magnetoelastic Effects in B20 compounds. ArXiv eprints, 1604.02766, April 2016.