==== Front Sci Rep Sci Rep Scientific Reports 2045-2322 Nature Publishing Group UK London 37393307 37808 10.1038/s41598-023-37808-4 Article Influence of CoFeB layer thickness on elastic parameters in CoFeB/MgO heterostructures Shekhar S. shashank.shekhar@amu.edu.pl 1 Mielcarek S. 1 Otani Y. 23 Rana B. 1 Trzaskowska A. 1 1 grid.5633.3 0000 0001 2097 3545 Institute of Spintronics and Quantum Information, Faculty of Physics, Adam Mickiewicz University, Uniwersytetu Poznańskiego 2, 61-614 Poznan, Poland 2 grid.7597.c 0000000094465255 Center for Emergent Matter Science, RIKEN, 2-1 Hirosawa, Wako, 351-0198 Japan 3 grid.26999.3d 0000 0001 2151 536X Institute for Solid State Physics, University of Tokyo, Kashiwa, Chiba 277-8581 Japan 1 7 2023 1 7 2023 2023 13 106685 4 2023 28 6 2023 © The Author(s) 2023 https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The surface acoustic waves, i.e., surface phonons may have huge potential for future spintronic devices, if coupled to other waves (e.g., spin waves) or quasiparticles. In order to understand the coupling of acoustic phonons with the spin degree of freedom, especially in magnetic thin film-based heterostructures, one needs to investigate the properties of phonons in those heterostructures. This also allows us to determine the elastic properties of individual magnetic layers and the effective elastic parameters of the whole stacks. Here, we study frequency versus wavevector dispersion of thermally excited SAWs in CoFeB/MgO heterostructures with varying CoFeB thickness by employing Brillouin light spectroscopy. The experimental results are corroborated by finite element method-based simulations. From the best agreement of simulation results with the experiments, we find out the elastic tensor parameters for CoFeB layer. Additionally, we estimate the effective elastic parameters (elastic tensors, Young’s modulus, Poisson’s ratio) of the whole stacks for varying CoFeB thickness. Interestingly, the simulation results, either considering elastic parameters of individual layers or considering effective elastic parameters of whole stacks, show good agreement with the experimental results. These extracted elastic parameters will be very useful to understand the interaction of phonons with other quasiparticles. Subject terms Materials science Magnetic properties and materials Surfaces, interfaces and thin films http://dx.doi.org/10.13039/501100013920 Uniwersytet im. Adama Mickiewicza w Poznaniu 054/13/SNŚ/0031 Shekhar S. http://dx.doi.org/10.13039/501100004281 Narodowe Centrum Nauki 2020/39/D/ST3/02378 Rana B. issue-copyright-statement© Springer Nature Limited 2023 ==== Body pmcIntroduction Surface acoustic waves (SAWs) are elastic waves that propagate along the surface of elastic materials with their displacement amplitude decaying with depth into the materials, so that the energy of acoustic phonons, associated with SAWs, is mostly confined in the vicinity of the surface. There are many types of SAWs such as Rayleigh waves, Sezawa waves, Pseudo-SAWs, Lamb waves, Love waves, and so on. One of the most used SAWs in modern devices is the Rayleigh waves, which is named after Lord Rayleigh who first reported the propagation and properties of SAWs in 18851. Today, SAWs are commonly used in microelectronics devices, sensors and filters because of their low power consumption, high sensitivity and broad tunability of operational frequency. SAW-based sensors can be utilized for characterizing various physical (e.g., density, viscosity) and chemical properties of materials2. Moreover, these sensors can be used to sense micro-pressures3 and to detect bacteria spores such as Bacillus thuringiensis and E. coli4. Apart from these, SAWs also promise to develop future spintronics devices if coupled to other waves, e.g., spin waves (SWs), i.e., the collective precessional motion of ordered magnetic spins in magnetic materials. It has been demonstrated that the SAWs can be a very useful tool for exciting5–9 and manipulating10–12 SWs and vice versa13–15, which can be a compliment to current CMOS based technology. Furthermore, SAWs have also proven to be a potential tool for the nucleation of magnetic skyrmions16, the creation of magnonic crystals, i.e., artificial crystals for tailoring magnonic band structures17,18, domain wall driving19,20, spin current generator21, nano-oscillator based reservoir computation22. Recently, magnetic thin films and their heterostructures are promising to be potential candidates for future spintronics applications because of their fascinating interfacial properties. Therefore, it is quite essential to understand how acoustic phonons couple with SWs, i.e., magnons and spin degree of freedom, especially, in magnetic thin film heterostructures. At the same time, it is also important to increase the coupling efficiency between phonons and spin degree of freedom. A significant amount of research has been dedicated to the investigation of the coupling between magnons and acoustic phonons in magnetic thin films. It was observed that the magnetic anisotropies and spatial profiles of SWs and acoustic phonons play a crucial role in the coupling phenomena23–26. Generally, the coupling strength can be increased by maximizing the overlapping area of magnon dispersion curves with the phonon dispersion curves. This can be achieved by tuning the dispersion curves for phonons and magnons by playing with the elastic and magnetic parameters of the materials, as properties of acoustic phonons and magnons strongly depend upon the elastic and magnetic parameters of the materials, respectively. CoFeB thin films are known to be one of the most promising ferromagnetic materials for future spintronics devices due to compelling features such as low Gilbert damping27, negligible magneto crystalline anisotropy, high tunnel magnetoresistance28 and large spin polarization29. They have shown potential for many applications like magnetic tunnel junctions27, racetrack memory30, magnetic random-access memory31, read head32, and so on. Therefore, it will be quite promising to investigate magnon–phonon coupling in CoFeB thin films and their heterostructures. While doing so, one should separately investigate phonon and magnon dispersion relations in CoFeB thin film-based heterostructures as a first step. Hence, we have adapted CoFeB/MgO heterostructures in this study for the investigation of phonon dispersion. This is stimulated by the fact that the elastic parameter of one thin layer may be drastically different from its bulk values. Moreover, the effective elastic property of the whole stacks could also be very different than the single layer. The most used technique to study SAWs and estimate the elastic properties of thin films is Brillouin light spectroscopy (BLS). This method allows us to measure acoustic phonons in a non-invasive manner and provides information about the frequency and wavevector of SAWs. Several studies have been performed using the BLS technique on the estimation of elastic parameters for different thin films such as ZnO33, SnO2 and SnS234, ITO35, TiN36; as well as for magnetic multilayers such as [Ni80Fe20/Au/Co/Au]1037, [CoFeB/Au]N38 and topological insulators39. Trzaskowska et al. observed that the change in the effective thickness of the layer (by varying the number of the repetitive layer) has an influence on the dispersion relation of the SAWs leaving the profiles of the elastic wave to be the same38. It was also observed that for multilayered thin films, the elastic properties strongly depend on their synthesis condition40–42, materials of the deposited layers, and the type of substrate43,44. Hence, estimating these elastic properties experimentally as well as theoretically is very important from the application point of view. In this study, we investigate the dispersion relation (the frequency versus wavevector) of SAWs in CoFeB/MgO heterostructures with varying CoFeB thickness by employing BLS. Finite element method (FEM) based simulations are performed to corroborate experimental results. From the best agreement of simulation results with the experiments we find out the elastic tensor parameters for CoFeB layer and estimate the effective elastic parameters of the whole stacks for varying CoFeB thickness. Material and methods Sample fabrication The samples for this study are multilayer structures deposited on a Si (001) substrate with 700 nm SiO2 on it. The multilayer stacking is as follows: Ta (10)/Co20Fe60B20 (t = 1, 1.4, 2, 3, 5, 10 and 20)/MgO (2)/Al2O3 (10), where the numbers in parentheses denote the thickness of the layers in nanometers. The Co20Fe60B20 will be denoted as CoFeB in the rest of the manuscript. The multilayers were deposited by radio frequency (rf) sputtering at a base pressure of 10−8 Torr at room temperature. The details of the sample preparation can also be found in Refs.45,46. A schematic view of the sample is presented in Fig. 1a.Figure 1 (a) The schematic view of the studied samples. The green lines represent the geometry of the experiment. (b) The BLS spectra obtained from the Si/SiO2/Ta/CoFeB(t)/MgO/Al2O3 samples at 32° angle (θ) of incidence of light (corresponds to q = 12.5 μm−1) for different thicknesses (t = 1, 5, 10 and 20 nm) of CoFeB layer. Experimental setup The thermally excited SAWs were studied using a six-pass, tandem Brillouin spectrometer (JRS Scientific Instruments), which ensures a contrast of 1015 refs.47,48. The incident light source was a Nd:YAG single-mode diode-pumped laser of 200 mW output power, which emits the second harmonic of wavelength λ0 = 532 nm (Excelsior, Spectra Physics). Measurements were performed in the backscattering geometry with pp polarization of light and the scattered light was collected using f/8 optics with a focal length of 58 mm. The solid angle of the lens was 0.63 steradians and the free spectral range was 20 GHz. A detailed description of the experimental setup can be found in Refs.49,50. BLS employs the inelastic scattering of incident photons from thermally excited phonons. The wavevector and frequency of the phonons in the studied samples can be characterized by measuring the projection of incident light wavevector q and frequency shift Δf of scattered light. Due to momentum conservation in the scattering process, the wavevector of acoustic waves is equal to the projection of the incident light wavevector in the sample plane. Thus, the wavevector q of the acoustic waves can be expressed as:1 q=4πsinθλ0 where θ is the angle of incidence of light with the normal to the sample surface (Fig. 1a). The angle θ was varied in the range of 10°–80° to change q of the measured phonons in the range of 4–23 µm−1. The phase velocity (νSAW) of SAWs can be correlated with θ by the following relation:2 ΔfSAW=2vSAWsinθλ0 where ΔfSAW is the Brillouin frequency shift for SAW51. Results For transparent materials, the incident light is mostly scattered in the bulk of the material via the photoelastic coupling mechanism, while for opaque materials, the light is scattered from the surface of the material by the surface ripple mechanism. This fact provides an opportunity to probe SAWs in opaque materials. One of the criteria to identify the type of waves is that the Rayleigh-type SAWs in a given material must show a linear relationship between frequency and wavevector49–52. Such criteria perfectly work for homogeneous materials. However, in the case of multilayer films on the substrate, linear dispersion relations are not usually observed. These systems can be divided into two categories: slow-on-fast (when the velocity of the transverse bulk wave in the multilayer is smaller than that of the substrate) and fast-on-slow system (vice versa)53. The next criterion for Rayleigh SAW is that νSAW should be lower than that of the velocity (νT) of the slowest transverse bulk waves. The typical spectra for the studied sample are presented in Fig. 1b. For all the studied samples, the spectra show the Rayleigh waves (marked by RW in the figure). In some spectra, the low-intensity Sezawa waves were also observed (marked by arrows) and identified from FEM-based simulations, which will be discussed later. The measured phase velocity dispersion of Rayleigh SAWs as a function of wavevector × total layer thickness (i.e., qh) is shown in Fig. 2. Here, q is the wavevector (according to Eq. 1) and h is the total layer thickness of the multilayer system, i.e., h=tTa+tCoFeB+tMgO+tAl2O3, where tTa, tCoFeB, tMgO, tAl2O3 are the thicknesses of Ta, CoFeB, MgO, Al2O3 layers, respectively.Figure 2 Phase velocity of Rayleigh SAWs obtained from the Si/SiO2/Ta/CoFeB(t)/MgO/Al2O3 samples for different wavenumbers. The solid curve represents the fit to an exponential decay function. The negative slope of phase velocity dispersion with qh suggests that the studied films can be classified as slow-on-fast systems. In these types of systems, the dispersion curves start with the transverse wave velocity (νT) in the layer and asymptotically decrease to the Rayleigh wave velocity (νR) in the layer deposited on the substrate54. By fitting the experimental data points with an exponential decay function, we extract that the νT (i.e., the velocity at qh = 0) in the studied layer is about 5000 m s−1, whereas the surface Rayleigh wave velocity in the overall stacking layer tends to be about 2500 m s−1. In our study, the penetration depth of Rayleigh waves is in the range of 250 nm to a few micrometers. According to the thickness of the SiO2 layer, we consider the substrate as Si/SiO2. To calculate the Rayleigh wave velocity, all the stacking layers should be considered as an effective layer instead of considering individual layers. This is worth mentioning here that the Rayleigh surface wave velocity in the overall stacking layer is quite different than in any individual layer of the stack. As the problem is not so trivial, we perform FEM-based simulations to determine the velocities of SAWs for large qh as described in section "Discussion". FEM simulation To theoretically understand the dispersion character of the acoustic waves in the studied multilayers, FEM-based simulations were performed using COMSOL Multiphysics software55. The unit cell selected for simulations consists of a long rectangular bar with dimensions 100 (x) × 100 (y) × 3742 (z) nm3 made of a multilayer or effective layer (see Fig. 3a). We consider the substrate to be a uniform elastic half-space on which the multilayers with the determined thickness of each layer are present. The model assumes perfectly bonded and ideally flat parallel layers with zero interfacial thickness, uniform respective layer thickness, no interfacial roughness and defects, and uniform elastic properties within a given layer.Figure 3 (a) The schematic diagram of the unit cell used for simulations. It consists of a long rectangular bar with dimensions 100 (x) × 100 (y) × 3742 (z) nm3 made of a multilayer. (b) Dispersion relation of SAWs obtained from the Si/SiO2/Ta/CoFeB (t = 20)/MgO/Al2O3 sample. Blue points represent the frequencies obtained from the FEM simulation, whereas red points are the frequencies obtained from the experiment. (c) The total displacement profile of SAW modes (Rayleigh wave, Sezawa wave) in studied samples is presented. The color scale is also provided. To reflect the conditions of the elastic semi-space, specific boundary conditions were applied. For the walls perpendicular to the free surface the Bloch-Floquet periodic boundary conditions were applied for each component of displacement:3 aexpiqxx+qyy,bexpiqxx+qyy, where a and b denote the components of the displacement in the x and y direction in the cartesian coordinate system and qx and qy are the wavevector components given as:4 qx=2πcosαλSAW,qy=2πcosβλSAW. Here α and β are the angles made by the wavevector with x and y axes, respectively and λSAW is the wavelength of SAW. The exponential decay of the SAW amplitude was taken into account by applying a fixed boundary condition to the wall opposite to the free surface. Figure 3b shows the dispersion relation of SAWs measured (red points) from the Si/SiO2/Ta/CoFeB (t = 20)/MgO/Al2O3 sample. In the studied multilayers both: the Rayleigh and Sezawa surface waves are visible. To calculate the dispersion characters of SAWs and understand the profile of the acoustic modes we first adapted the elastic parameters and densities of each layer as given in Table 1. Here, the parameters for all the materials, excluding CoFeB, are taken from indicated references. The elastic tensor components for CoFeB were calculated from the best fitting of frequency versus wavevector characters for all the multilayers with varying CoFeB thickness and it is described in the Discussion section. The calculated dispersion relation obtained from the simulations is shown in Fig. 3b by solid blue dots. A very good agreement is obtained between the experimental and simulation dispersion characters, especially for Rayleigh waves, using the estimated elastic parameters of CoFeB. We note here that the materials in the studied samples have different crystallographic symmetry, so the parameters have strickle characters to the materials and their symmetry. Two-dimensional maps of deformation for SAWs are presented in Fig. 3c. A difference in the deformation profile obtained for Rayleigh and Sezawa waves is clearly visible.Table 1 Elastic constant and density of materials used for the simulation. Silicon56 Silicon dioxide57 Tantalum58 MgO59 Al2O360 Elastic constants (in GPa)  c11 166 75 260.91 295.9 466  c22 166 75 260.91 295.9 466  c33 166 75 260.91 295.9 506  c12 64 15 157.43 95.4 127  c13 64 15 157.43 95.4 117  c23 64 15 157.43 95.4 117  c44 80 30 81.82 153.9 235  c55 80 30 81.82 153.9 235  c66 80 30 51.74 153.9 169.5  c14=-c24=c56 0 0 0 0 94 Density (in kg m−3) 2332 2650 16,600 3580 3950 In Fig. 4, the phase velocity dispersions of SAWs for Si/SiO2/Ta/CoFeB (t = 20)/MgO/Al2O3 film are plotted. Here, the red points represent the experimental results, whereas the blue points represent the results obtained from FEM simulations. The dispersion plots reconfirm that the studied multilayer is a slow-on-fast system. The main criterion of the above classification is vTlayer<2·vTsubstrate for slow-on-fast systems, while the criterion for the fast-on-slow system is vTlayer>2·vTsubstrate according to Farnell and Adler53. For slow-on-fast systems, the presence of an elastically soft layer on an elastically hard substrate leads to the reduction of Rayleigh SAWs velocity and formation of higher-order modes known as Sezawa waves, whose displacement profile is presented in Fig. 3c. The simulated phase velocity dispersion graphs allow us to estimate the phase velocities of transverse bulk waves in the layer and in the substrate, which are about 5000 m s−1 and 5740 m s−1, respectively. The phase velocity of Rayleigh SAW in the multilayer is estimated to be 2350 m s−1. Here, the substrate is defined as Si/SiO2, and the layer is defined as Ta/CoFeB/MgO/Al2O3. The estimated values of velocities for all studied samples are the basis for calculations of the elastic parameters of the CoFeB layer and effective layers.Figure 4 Phase velocity dispersion of SAWs, obtained from Si/SiO2/Ta/CoFeB (t = 20)/MgO/Al2O3 sample, are plotted as a function of qh. Red dots are experimental results obtained from BLS measurements, whereas blue dots are simulation results performed by FEM. Discussion For an in-depth understanding of the effect of CoFeB layer thickness on the velocity of SAWs in the studied materials, some additional simulations on SAW propagation were performed by FEM. The penetration depths of Rayleigh-type SAWs in the materials partly depend on the wavevector. We have previously mentioned that the substrate must be taken into account while calculating the characters of SAWs, especially for smaller wavevectors. As the penetration depth is way larger than the total thickness of the multilayer (which is a maximum of 42 nm for a multilayer with 20 nm thick CoFeB layer), the studied multilayers Ta/CoFeB/MgO/Al2O3 can be treated as an effective layer on the substrate. In that way, the effective elastic parameter of the multilayer can be estimated. To do so we first calculate the elastic parameters of CoFeB under two conditions. The best agreement between the experimental results and FEM simulations is treated as the first condition. The second condition is to fulfill the Cauchy relations61,62. We started performing simulations with the value of elastic tensor components of CoFeB as proposed in Ref.23. However, these parameters didn’t fulfill our first condition. So, we modified the parameters in such a way that it fulfills both the above-mentioned conditions. The estimated elastic tensor components for CoFeB are (in GPa):210130130000130210130000130130210000000450000004500000045 Based on these parameters, the calculated phase velocities of bulk acoustic modes and Rayleigh SAW, which propagate on the surface (001) are presented in Fig. 5.Figure 5 The phase velocity of bulk acoustic waves (T1, T2—transverse waves) and Rayleigh wave (R SAW) for CoFeB (a); and SAWs (b) obtained for effective layers with different thicknesses of CoFeB. The value of Rayleigh SAW in CoFeB, calculated according to FEM simulations, is around 2350 m s−1 (Fig. 5a). Interestingly, the same value is also obtained by plotting phase velocity in the effective layer as a function of CoFeB thickness and fitting with exponential function as depicted in Fig. 5b. These results suggest that the layers in the studied samples can be treated as an effective layer. The CoFeB layer has cubic symmetry, but when the CoFeB is included in a multilayer system it becomes more convenient to use Young’s modulus which can be calculated from the estimated elastic tensors of CoFeB. In our case, Young’s modulus E, as described in63, is given as:5 E=σ33ε33=1s33′ where s33′ is the component of stiffness in the measurement plane of the system (in Voigt notation). For a material with cubic symmetry, the s33′ component of stiffness is given as63:6 s33′=s3333′=s11-(s11-s12-12s44)[1-cos4φsin4ψ-sin4φsin4ψ-cos4ψ] where s11,s12, and s44 are the components of stiffness in the crystal plane of the system and φ and ψ are the two arbitrary rotation angles to transform the cubic crystal coordinate system into a lattice plane system. The relations among stiffness components s11,s12, and s44 and elastic constants c11,c12 and c44 are given as:7 s11=c11+c12c11-c12c11+2c12s12=-c12c11-c12c11+2c12s44=1/c44. Using Eqs. (5)–(7), Young’s modulus E can be written as:8 E=2c11-c12c11-c12c442c11+c12-c11+2c122c44-c11+c121-cos4φsin4ψ-sin4φsin4ψ-cos4ψ With the estimated elastic components, Young’s modulus E is plotted for CoFeB as shown in Fig. 6.Figure 6 (a) The 3D view and (b) view in (001) plane of calculated Young’s modulus (E) of CoFeB. We calculate the elastic parameters of the effective layers for varying thicknesses of the CoFeB layer. This is done by adapting the proportion method, where elastic parameters of individual layers are multiplied by the volume fraction of the effective layer, added thereafter, and averaged out with the total volume64,65. There is no unique way to do this, and different literatures propose different procedures66,67. We used ELATE68, which is an open-source online tool for the analysis of elastic tensors, to calculate Young’s modulus and Poisson’s ratio of each layer. Then, we calculated their weighted average based on the thickness of each layer. Figure 7 represents the plot of the calculated values of elastic tensors of the effective layer as a function of total layer thickness (h), which shows a significant variation of effective elastic parameters with h. The accuracy of the elastic tensor values is in the range of 5 GPa. Next, we compare the phase velocity dispersion of SAWs calculated by considering the elastic parameters of individual layers and the elastic parameters of the effective layer (Fig. 8). Both results show great agreement with each other.Figure 7 Calculated values of elastic tensor components for effective layer as a function of total or effective layer thickness h. Figure 8 The phase velocity of SAWs obtained from the Si/SiO2/Ta/Co20Fe60B20 (t = 20)/MgO/Al2O3 samples as individual layers (circles) and as effective layer (triangles). The square points represent the phase velocity obtained from the experiment. Finally, we calculate Young’s modulus and Poisson’s ratio of the effective layer. A multilayer system can consist of various layers with different crystallographic symmetries and different numbers of independent elastic tensor components. However, for the sake of simplicity, we can treat the effective layer as isotropic and define elastic parameters and densities for it. Although, different studies propose different procedures69,70, we employ above mentioned weighted average technique and the calculated values of Young’s modulus and Poisson’s ratio are presented in Table 2. These values will be very useful to analyze the coupling between phonons and magnons.Table 2 Young’s modulus and Poisson’s ratio for uniform effective layers (isotropic). CFB layer thickness (nm) Density (kg m−3) Young’s modulus (GPa) Poisson’s ratio 1 9549 249 0.31 1.4 9505 249 0.31 2 9442 243 0.31 3 9343 238 0.31 5 9167 228 0.32 10 8824 210 0.33 20 8382 187 0.34 Conclusions In conclusion, we have studied frequency versus wavevector dispersion of thermally excited SAWs in Si/SiO2/Ta/CoFeB(t)/MgO/Al2O3 heterostructures with varying CoFeB thickness (t) from 1 to 20 nm. The thermally excited acoustic phonons are detected by employing Brillouin light spectroscopy in backscattering geometry. Both Rayleigh-type and Sezawa-type SAWs are detected in the experiment and their frequency versus wavevector dispersions are measured up to a wavevector of 23 μm−1. The experimental results are corroborated by finite element method based COMSOL simulations. The simulated mode profiles of Rayleigh and Sezawa waves show very distinct spatial features. We calculate the elastic tensor parameters for the CoFeB layer from the best agreement of simulation results with the experiment data points. We further estimate the effective elastic parameters (elastic tensors, Young’s modulus, Poisson’s ratio) of the whole stacks for varying CoFeB thickness by considering the whole stack as a uniform layer. It is observed that the effective elastic parameters significantly vary with the CoFeB layer thickness. Interestingly, the simulation results, either considering elastic parameters of individual layers or considering effective elastic parameters of whole stacks, show good agreement with the experimental results. These extracted elastic parameters will be very useful to understand the interaction of phonons with other quasiparticles such as magnons and skyrmions in a similar type of heterostructures. Acknowledgements SS gratefully acknowledges the financial support from Adam Mickiewicz University under the framework of “Initiative of Excellence—Research University” (ID-UB 054/13/SNŚ/0031). BR acknowledges financial support from the NCN SONATA-16 project with grant number 2020/39/D/ST3/02378. The authors sincerely thank Dr. Katsuya Miura and Dr. Hiromasa Takahashi from Hitachi Ltd., Tokyo, Japan, for providing thin films for the study. Author contributions B.R. and A.T. planned the studies. Experiments, simulations and analysis are carried out by S.S. and A.T. First version of manuscript is prepared by S.S. All the authors discussed the results and contributed to write and improve the quality of manuscript. Data availability The dataset used and/or analyzed during the current study will be made available by the corresponding author on reasonable request. Competing interests The authors declare no competing interests. Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. ==== Refs References 1. Rayleigh L On waves propagated along the plane surface of an elastic solid Proc. Lond. Math. Soc. 1885 s1-17 4 11 10.1112/plms/s1-17.1.4 2. Hossenlopp JM Applications of acoustic wave devices for sensing in liquid environments Appl. Spectrosc. Rev. 2006 41 151 164 10.1080/05704920500510745 3. Caliendo C Hamidullah M Pressure sensing with zero group velocity lamb modes in self-supported a-SiC/c-ZnO membranes J. Phys. D Appl. Phys. 2018 51 385102 10.1088/1361-6463/aad6f3 4. Branch DW Brozik SM Low-level detection of a Bacillus anthracis simulant using Love-wave biosensors on 36°YX LiTaO3 Biosens. Bioelectron. 2004 19 849 859 10.1016/j.bios.2003.08.020 15128104 5. Li X Labanowski D Salahuddin S Lynch CS Spin wave generation by surface acoustic waves J. Appl. Phys. 2017 122 043904 10.1063/1.4996102 6. Tateno S Nozaki Y Nozaki Y Highly nonreciprocal spin waves excited by magnetoelastic coupling in a Ni/Si bilayer Phys. Rev. Appl. 2020 13 034074 10.1103/PhysRevApplied.13.034074 7. Xu M Inverse Edelstein effect induced by magnon-phonon coupling Phys. Rev. B 2018 97 180301 10.1103/PhysRevB.97.180301 8. Dreher L Surface acoustic wave driven ferromagnetic resonance in nickel thin films: Theory and experiment Phys. Rev. B Condens. Matter Mater. Phys. 2012 86 134415 10.1103/PhysRevB.86.134415 9. Weiler M Elastically driven ferromagnetic resonance in nickel thin films Phys. Rev. Lett. 2011 106 117601 10.1103/PhysRevLett.106.117601 21469894 10. Froes D Arana M Sampaio LC Sinnecker JP Acoustic wave surfing: spin waves and spin pumping driven by elastic wave J. Phys. D Appl. Phys. 2021 54 255001 10.1088/1361-6463/abed71 11. Küß M Nonreciprocal magnetoacoustic waves in dipolar-coupled ferromagnetic bilayers Phys. Rev. Appl. 2021 15 034060 10.1103/PhysRevApplied.15.034060 12. Hatanaka D On-chip coherent transduction between magnons and acoustic phonons in cavity magnomechanics Phys. Rev. Appl. 2022 17 034024 10.1103/PhysRevApplied.17.034024 13. Gowtham PG Labanowski D Salahuddin S Mechanical back-action of a spin-wave resonance in a magnetoelastic thin film on a surface acoustic wave Phys. Rev. B 2016 94 014436 10.1103/PhysRevB.94.014436 14. Rovillain P Impact of spin-wave dispersion on surface-acoustic-wave velocity Phys. Rev. Appl. 2022 18 064043 10.1103/PhysRevApplied.18.064043 15. Xu M Nonreciprocal surface acoustic wave propagation via magneto-rotation coupling Sci. Adv. 2020 6 eabb1724 10.1126/sciadv.abb1724 32821833 16. Yokouchi T Creation of magnetic skyrmions by surface acoustic waves Nat. Nanotechnol. 2020 15 361 366 10.1038/s41565-020-0661-1 32231267 17. Kryshtal RG Medved AV Surface acoustic wave in yttrium iron garnet as tunable magnonic crystals for sensors and signal processing applications Appl. Phys. Lett. 2012 100 192410 10.1063/1.4714507 18. Chumak AV Dhagat P Jander A Serga AA Hillebrands B Reverse Doppler effect of magnons with negative group velocity scattered from a moving Bragg grating Phys. Rev. B Condens. Matter. Mater. Phys. 2010 81 140404 10.1103/PhysRevB.81.140404 19. Adhikari A Adenwalla S Surface acoustic waves increase magnetic domain wall velocity AIP Adv. 2021 11 015234 10.1063/9.0000159 20. Vilkov E Magnetic domain wall motion driven by an acoustic wave Ultrasonics 2022 119 106588 10.1016/j.ultras.2021.106588 34601287 21. Kobayashi D Spin current generation using a surface acoustic wave generated via spin-rotation coupling Phys. Rev. Lett. 2017 119 077202 10.1103/PhysRevLett.119.077202 28949686 22. Chowdhury MFF Focused surface acoustic wave induced nano-oscillator based reservoir computing Appl. Phys. Lett. 2022 121 102402 10.1063/5.0110769 23. Babu NKP The interaction between surface acoustic waves and spin waves: The role of anisotropy and spatial profiles of the modes Nano. Lett. 2021 21 946 951 10.1021/acs.nanolett.0c03692 33231459 24. Chen J Hybridized propagation of spin waves and surface acoustic waves in a multiferroic-ferromagnetic heterostructure Phys. Rev. Appl. 2023 19 024046 10.1103/PhysRevApplied.19.024046 25. Babu NKP Interaction between thermal magnons and phonons in a CoFeB/Au multilayer IEEE Magn. Lett. 2019 10 1 5 10.1109/LMAG.2019.2950304 26. Yamamoto K Xu M Puebla J Otani Y Maekawa S Interaction between surface acoustic waves and spin waves in a ferromagnetic thin film J. Magn. Magn. Mater. 2022 545 168672 10.1016/j.jmmm.2021.168672 27. Ikeda S A perpendicular-anisotropy CoFeB–MgO magnetic tunnel junction Nat. Mater. 2010 9 721 724 10.1038/nmat2804 20622862 28. Ikeda S Tunnel magnetoresistance of 604% at 300K by suppression of Ta diffusion in CoFeB/MgO/CoFeB pseudo-spin-valves annealed at high temperature Appl. Phys. Lett. 2008 93 082508 10.1063/1.2976435 29. Huang SX Chen TY Chien CL Spin polarization of amorphous CoFeB determined by point-contact Andreev reflection Appl. Phys. Lett. 2008 92 242509 10.1063/1.2949740 30. Zhang Y Perpendicular-magnetic-anisotropy CoFeB racetrack memory J. Appl. Phys. 2012 111 093925 10.1063/1.4716460 31. Iwata-Harms JM Ultrathin perpendicular magnetic anisotropy CoFeB free layers for highly efficient, high speed writing in spin-transfer-torque magnetic random access memory Sci. Rep. 2019 9 19407 10.1038/s41598-019-54466-7 31857596 32. Tsunekawa K Giant tunneling magnetoresistance effect in low-resistance CoFeB/MgO(001)/CoFeB magnetic tunnel junctions for read-head applications Appl. Phys. Lett. 2005 87 072503 10.1063/1.2012525 33. Carlotti, G., Socino, G., Petri, A. & Verona, E. Elastic constants of sputtered ZnO films. Ultrasonics Symposium Proceedings 295–299. 10.1109/ULTSYM.1987.198972 (1987). 34. Mnari M Cros B Amlouk M Belgacem S Barjon D Study of the elastic properties of sprayed SnO2 and SnS2 layers Can. J. Phys. 2011 77 705 715 10.1139/p99-023 35. Wittkowski T Elastic properties of indium tin oxide films Thin Solid Films 2001 398–399 465 470 10.1016/S0040-6090(01)01373-6 36. Elena M Elastic behaviour of TiN thin films Thin Solid Films 1993 236 209 213 10.1016/0040-6090(93)90671-B 37. Trzaskowska A Mielcarek S Graczykowski B Stobiecki F Surface waves investigation in NiFe/Au/Co/Au multilayers by high-resolution Brillouin spectroscopy J. Alloys Compd. 2012 517 132 138 10.1016/j.jallcom.2011.12.059 38. Trzaskowska A The studies on phonons and magnons in [CoFeB/Au]N multilayers of different number of repetitions J. Magn. Magn. Mater. 2022 549 169049 10.1016/j.jmmm.2022.169049 39. Trzaskowska A Mielcarek S Wiesner M Lombardi F Mroz B Dispersion of the surface phonons in semiconductor/topological insulator Si/Bi2Te3 heterostructure studied by high resolution Brillouin spectroscopy Ultrasonics 2021 117 106526 10.1016/j.ultras.2021.106526 34303926 40. Schultrich B Scheibe H-J Grandremy G Schneider D Elastic modulus of amorphous carbon films Phys. Status Solidi A 1994 145 385 392 10.1002/pssa.2211450219 41. Lofaj F Kabátová M Kvetková L Dobrovodský J The effects of deposition conditions on hydrogenation, hardness and elastic modulus of W-C: H coatings J. Eur. Ceram. Soc. 2020 40 2721 2730 10.1016/j.jeurceramsoc.2019.12.062 42. Huang H Effect of deposition conditions on mechanical properties of low-temperature PECVD silicon nitride films Mater. Sci. Eng. A 2006 435–436 453 459 10.1016/j.msea.2006.07.015 43. Liu K Elastic properties of chemical-vapor-deposited monolayer MoS2, WS2, and their bilayer heterostructures Nano Lett. 2014 14 5097 5103 10.1021/nl501793a 25120033 44. Mukhopadhyay T Mahata A Adhikari S Zaeem MA Effective mechanical properties of multilayer nano-heterostructures Sci. Rep. 2017 7 15818 10.1038/s41598-017-15664-3 29150623 45. Rana B Electric field control of spin waves in ultrathin CoFeB films Phys. Rev. B 2019 100 224412 10.1103/PhysRevB.100.224412 46. Rana B Nonlinear control of damping constant by electric field in ultrathin ferromagnetic films Phys. Rev. Appl. 2020 14 014037 10.1103/PhysRevApplied.14.014037 47. Sandercock, J. R. Trends in brillouin scattering: Studies of opaque materials, supported films, and central modes. 173–206. 10.1007/3540115137_6 (1982). 48. Scarponi F High-performance versatile setup for simultaneous Brillouin-Raman microspectroscopy Phys. Rev. X 2017 7 031015 10.1103/PhysRevX.7.031015 49. Trzaskowska A Mielcarek S Sarkar J Band gap in hypersonic surface phononic lattice of nickel pillars J. Appl. Phys. 2013 114 134304 10.1063/1.4824103 50. Mielcarek S Trzaskowska A Mroz B Andrews T High resolution Brillouin scattering studies of β-Gd2(MoO4)3; the bulk and surface phase transitions J. Phys. Condens. Matter 2005 17 587 10.1088/0953-8984/17/4/003 51. Every AG Measurement of the near-surface elastic properties of solids and thin supported films Meas. Sci. Technol. 2002 13 R21 10.1088/0957-0233/13/5/201 52. Mielcarek S Trzaskowska A Graczykowski B Sarkar J Hypersonic surface waves in 2D titanium nanostructure on silicon Phys. Status Solidi (RRL) Rapid Res. Lett. 2012 6 175 177 10.1002/pssr.201206039 53. Farnell GW Adler EL Elastic wave propagation in thin layers Phys. Acoust. 1972 9 35 127 10.1016/B978-0-12-395670-5.50007-6 54. Kundu T Ultrasonic Nondestructive Evaluation: Engineering and Biological Material Characterization 2003 Boca Raton CRC Press 55. COMSOL Multiphysics® software, COMSOL AB, Sweden. 56. Kaajakari, V. Silicon as an anisotropic mechanical material: A tutorial. http://www.kaajakari.net/~ville/research/tutorials/tutorials.shtml (Accessed 2023). 57. Carlotti G Doucet L Dupeux M Elastic properties of silicon dioxide films deposited by chemical vapour deposition from tetraethylorthosilicate Thin Solid Films 1997 296 102 105 10.1016/S0040-6090(96)09346-7 58. Featherston FH Neighbours JR Elastic constants of tantalum, tungsten, and molybdenum Phys. Rev. 1963 130 1324 1333 10.1103/PhysRev.130.1324 59. Marklund K Mahmoud S. A. Elastic constants of magnesium oxide Phys. Scr. 1971 3 75 10.1088/0031-8949/3/2/005 60. Rao RVGS Elastic constants of alumina Proc Indian Acad Sci Sect A 1949 29 352 360 10.1007/BF03172465 61. Shimizu H Ohnishi M Sasaki S Ishibashi Y Cauchy relation in dense H2O ice VII Phys. Rev. Lett. 1995 74 2820 2823 10.1103/PhysRevLett.74.2820 10058026 62. Krüger JK Cauchy-like relation between elastic constants in amorphous materials Phys. Rev. B 2002 66 12206 10.1103/PhysRevB.66.012206 63. Zhang JM Zhang Y Xu KW Ji V Young’s modulus surface and Poisson’s ratio curve for cubic metals J. Phys. Chem. Solids 2007 68 503 510 10.1016/j.jpcs.2007.01.025 64. Hou Z Wu F Fu X Liu Y Effective elastic, parameters of the two-dimensional phononic crystal Phys. Rev. E 2005 71 037604 10.1103/PhysRevE.71.037604 65. Pereira A Costa M Anflor C Pardal J Leiderman R Estimating the effective elastic parameters of nodular cast iron from micro-tomographic imaging and multiscale finite elements: Comparison between numerical and experimental results Metals 2018 8 695 10.3390/met8090695 66. Wright OB Matsuda O Watching surface waves in phononic crystals Philos. Trans. R. Soc. A 2015 373 20140364 10.1098/rsta.2014.0364 67. Głowacki J Tomanik M Pezowicz C Krauss H Mechanical and histomorphometrical evaluation of false and floating ribs of young adults with idiopathic scoliosis Acta Bioeng. Biomech. 2020 22 3 10 10.37190/ABB-01575-2020-01 32868948 68. Gaillac R Pullumbi P Coudert FX ELATE: an open-source online application for analysis and visualization of elastic tensors J. Phys. Condens. Matter 2016 28 275201 10.1088/0953-8984/28/27/275201 27199239 69. Shrikanth S Knowles KM Neelakantan S Prasad R Planes of isotropic Poisson’s ratio in anisotropic crystalline solids Int. J. Solids Struct. 2020 191–192 628 645 10.1016/j.ijsolstr.2019.10.014 70. Norris A. N. Extreme values of Poisson’s ratio and other engineering moduli in anisotropic materials J. Mech. Mater. Struct. 2006 1 793 812 10.2140/jomms.2006.1.793