On giant amplification of effects of phonon nonspecular reflection at a solitary interface between magnetic and non-magnetic media
- Authors: Gulyaev Y.V.1, Sukhorukova O.S.2,3, Tarasenko A.S.3, Tarasenko S.V.3, Shavrov V.G.1
-
Affiliations:
- Institute of Radio Engineering and Electronics named after Kotelnikov, Russian Academy of Sciences
- Donetsk State University
- Donetsk Institute for Physics and Engineering named after Galkin
- Issue: Vol 514, No 1 (2024)
- Pages: 83-90
- Section: ТЕХНИЧЕСКИЕ НАУКИ
- URL: https://journals.rcsi.science/2686-7400/article/view/261451
- DOI: https://doi.org/10.31857/S2686740024010138
- EDN: https://elibrary.ru/OIBZQK
- ID: 261451
Cite item
Full Text
Abstract
On the secluded interface of semi-infinite magnetic and non-magnetic media, hybridiza-tion of magnetoelastic interaction with inhomogeneous exchange or magneto-dipole interaction can lead to the formation of symmetry-protected bound states in the radiation spectrum of the leaky surface magnon polarons. If a quasi-plane bulk elastic wave falls on the surface of a magnetic medium from outside, and the wave parameters approach to the parameters of the surface “dark” state of the above mentioned type, the nonspecular reflection effects of the first order increase unrestrictedly (in the non-dissipative approximation).
Full Text
В настоящее время практическая реализация идеи использования пространственного перемещения не заряда, а вращательного спинового момента с целью создания на основе магнитных диэлектриков быстродействующих, энергоэффективных устройств функциональной электроники существенно связана с физикой вытекающих магнонных поляронов [1]. В частности, это объясняется тем, что роль подобных гибридных магнон-фононных состояний в акустике слоистых магнитных сред во многом сопоставима с тем вкладом, который в проводимость ионных кристаллов вносят “истинные” носители заряда – поляроны (“электроны в фононной шубе”), на возможность существования которых было впервые указано в [2]. Как следствие, в магнонике, в полной аналогии с акустоэлектроникой, оптимизация работы широкой гаммы устройств, связанных с обработкой сигналов и использующих связь между акустическими волнами и носителями тока, требует детального рассмотрения связи между дисперсионными свойствами распространяющегося вдоль слоистой магнитной гетероструктуры вытекающего магнонного полярона и структурой поля его фононного излучения. В настоящее время особое внимание исследователей открытых волноводных структур в спинтронике привлекает изучение как условий формирования в сплошном спектре излучения распространяющейся вытекающей волны локализованных состояний, обладающих нулевой радиационной шириной (связанных состояний в континууме (ССК)), так и возникающие в этой связи динамические эффекты. Недавно в [3] на примере полностью отражающей, акустически полуоткрытой магнитной сэндвич-структуры (магнитный слой между двумя акустически неэквивалентными упругими полупространствами) в бездиссипативном приближении было показано, что совместный учет магнитоупругого и неоднородного обменного взаимодействий может приводить к формированию интерференционного типа связанных состояний [4] в сплошном спектре фононного излучения вытекающих поверхностных магнонных поляронов, тогда как в окрестности таких “темных” состояний их радиационная ширина может быть сколь угодно малой. В условиях полного внутреннего отражения (ПВО) для падающей извне на поверхность полуоткрытого магнитного слоя плоской однопарциальной объемной упругой волны c частотой ω и волновым вектором k коэффициент отражения имеет вид (один открытый канал излучения)
, (1)
где Zin – входной поверхностный акустический импеданс отражающей структуры, тогда как – поверхностный волновой импеданс полуограниченной акустически менее плотной среды, поддерживающей распространение плоской объемной однопарциальной упругой волны [5]. Пусть q – нормаль к границе раздела сред, a – нормаль к плоскости падения волны. В этом случае, в соответствии с определением “темных” состояний, на плоскости внешних параметров “ω − h” (h ≡ (kb), b = [qa]) структуру входного поверхностного волнового импеданса в (1) в окрестности точки формирования поверхностного ССК (,) и заданном ω, согласно [3], можно представить как
(2)
Расчет показывает, что как при Zin (h, ω) = 0, так и при Zin–1 (h, ω) = 0 мгновенный поток энергии через границу раздела “магнитный слой – акустически менее плотная среда” строго равен нулю. Как известно, в кристаллоакустике существует понятие особой объемной волны (ООВ), характеризующейся аналогичным свойством нормальной компоненты потока энергии на границе раздела сред с акустическим контрастом, но в отражающей среде глубина локализации ООВ отсутствует (см., например, [6]). Учитывая это, в дальнейшем дисперсионные кривые, определяемые на основании (2) как h = hN(ω) (Zin (h, ω) = 0) и h = hD (ω) (Zin–1 (h, ω) = 0), будем называть особой поверхностной волной (см. также [3, 7]) первого (ОПВ-I) и второго (ОПВ-II) типа соответственно. По мере приближения к точке формирования поверхностного ССК (hN(ω) = hD(ω)) в соответствии с (1), (2) и определением ССК
|hD(ω) − hN(ω)| → 0. (3)
Как следует из результатов [3], такое поведение является примером формирования ССК по механизму Фридриха–Винтена [4], поскольку в точке ССК на плоскости внешних параметров “ω − h” имеет место вырождение дисперсионных кривых, отвечающих спектрам вытекающих ОПВ первого и второго типов соответственно. Так как при этом и числитель и знаменатель поверхностного входного волнового импеданса (2) одновременно обращаются в ноль, то в условиях ПВО обращается в ноль также и вклад радиационных полей излучения ОПВ-I и ОПВ-II в сплошной спектр излучательных фононных мод рассматриваемой полуоткрытой структуры (следствие деструктивной интерференции).
Хорошо известно [5], что если из акустически менее плотной среды на поверхность полностью отражающей слоистой структуры падает квазиплоская объемная упругая волна, то после отражения с коэффициентом (1) она будет испытывать лежащее в плоскости падения волны продольное (в направлении b) смещение вдоль границы раздела сред (эффект Шоха):
(4)
В результате при заданном ω и сканировании по углу вследствие (2) ∆(ω, h) будет достигать максимума как на линии ОПВ-I (при h → hN(ω)), так и на линии ОПВ-II (при h → hD(ω)). По мере приближения к точке существования поверхностного ССК (3) величина каждого из этих экстремумов неограниченно (в рамках рассматриваемой модели) возрастает [3].
Если в условиях, подобных выше указанным, на поверхность магнетика падает не квазиплоская, а квазимонохроматическая волна, то вместо (4) будет иметь место запаздывание по времени отраженного импульса с плоским фронтом по отношению к падающему (Вигнеровское запаздывание):
(5)
Используя полученные выше соотношения, несложно убедиться, что в условиях (1) при частотном сканировании и заданном угле падения квазимонохроматической волны (длительность импульса ) на поверхность полуограниченной полностью отражающей структуры в результате существования ОПВ первого и второго типа соотношение (2) можно представить в виде Zin ≈ C(ω(h) – ωN(h))[ω(h) – ωD(h)]–1. Таким образом, для отраженного фононного импульса локальные максимумы запаздывания отраженного сигнала на несущей частоте (эффект Вигнера (5)) c учетом (2) также достигаются как при Zin (h, ω) = 0 (ОПВ-I), так и при Zin–1 (h, ω) = 0 (ОПВ-II). По мере приближения параметров падающего фононного импульса к значениям, отвечающим поверхностному ССК: |ωD(h) − ωN(h)| → 0, максимальное время запаздывания отраженного сигнала на несущей частоте, согласно (5), также будет неограниченно (в рамках рассматриваемой модели) возрастать () [3].
Вместе с тем указанные результаты работы [3] были получены на основе следующих, существенных для целей предлагаемого сообщения, ограничений. Во-первых, осталось неясным, возможно ли такое же резкое усиление незеркальных эффектов первого порядка при полном внутреннем отражении объемных фононов, если имеется только одна граница раздела полуограниченных магнитной и немагнитной сред. Во-вторых, в [3] для отражающей среды плоскость падения упругой с узким угловым или частотным спектром волны отвечала лишь высокосимметричному направлению. Наконец, в-третьих, согласно предложенной в [4] классификации возможных типов ССК, в [3] рассматривались лишь ССК интерференционного типа (механизм Фридриха–Винтгена), однако, помимо ССК, полученных методом “inverse construction”, также в принципе еще возможны и симметрийно защищенные ССК [4]. Однако несмотря на все вышесказанное, возможность существования аналога найденному в [3] сверхусилению незеркальных эффектов фононного отражения первого порядка в области ПВО, но теперь для случая формирования за счет одновременного учета МУ и неоднородного обменного взаимодействий симметрийно защищенных ССК на уединенной границе раздела магнитной и немагнитной сред, до сих пор не изучалась. Ответ на вышеперечисленные вопросы на примере структуры с одним открытым каналом излучения и является целью данного сообщения.
Для простоты и наглядности расчетов предположим, что имеется два полупространства с вектором нормали к плоскости границы раздела q || OY, причем верхнее (y > 0) занято идеальной жидкостью (соответствующие величины будем обозначать знаком тильда), в которой потенциал смещений , плотность и модуль сжатия – и , а вектор смещения и тензор напряжений – ũ и связаны соотношениями [5]:
(6)
Что же касается нижнего полупространства (y < 0), то будем полагать, что оно занято пространственно однородным двухосным антиферромагнетиком (АФМ), легкая магнитная ось которого (OZ) лежит в плоскости границы раздела сред. Чтобы рассмотреть случай симметрийно защищенного поверхностного ССК, являющегося результатом одновременного учета МУ и неоднородного обменного взаимодействий, рассмотрим двухподрешеточную модель (|M1|=|M2|=M0, М0 – намагниченность насыщения подрешеток М1,2) обменно коллинеарного центросимметричного АФМ, обладающего изотропным упругим и МУ взаимодействиями с плотностью термодинамического потенциала в терминах векторов ферро- (m = (М1 + М2) / 2М0) и антиферромагнетизма (l = (M1 − M2) / 2M0) вида [8]:
(7)
Здесь δ, a, b – константы межподрешеточного обмена и двухосной магнитной анизотропии, γ – константа МУ-взаимодействия, λ, μ – коэффициенты Ламе. МУ-динамика рассматриваемой модели АФМ-среды (2) описывается замкнутой системой уравнений, состоящей из основного уравнения механики сплошной среды и уравнений Ландау–Лифшица для векторов m и l [8]:
(8)
Здесь ρ – плотность, u – вектор упругих смещений, g – гиромагнитное отношение, которое будем считать одинаковым для обеих подрешеток, Hr ≡ −δW / δr – эффективное поле, r = m, l. В декартовой системе координат равновесный вектор АФМ l0 направлен вдоль оси OZ, а легкая магнитная плоскость АФМ совпадает с YOZ. Если в модели АФМ (7), (8) пренебречь взаимодействием высокочастотных магнонов (b/a < 1) с фононами, то в этом случае итоговые материальные соотношения будут такие же, как и для упругоизотропной немагнитной среды (см., например, [9]), но с тем исключением, что теперь для упругой волны с частотой ω и произвольно ориентированным волновым вектором k модуль упругости c44 будет эффективным, одновременно обладая как пространственной, так и временной дисперсией:
(9)
где – индуцированная легкоосной анизотропией энергия активации рассматриваемой спиновой волны, – МУ-щель, – скорость спиновых волн [8]. В результате соотношение для спектра нормальных МУ-волн с частотой и произвольно ориентированным волновым вектором k в неограниченной АФМ-среде (7) с учетом (9) можно представить в виде
|Λik (ω, k) – ρω2Iik│ = 0. (10)
Здесь Λik – тензор Кристоффеля, отвечающий за исключением модуля c44 (см. (9)) упругоизотропной немагнитной среде [9], Iik – единичный тензор. В эластостатическом пределе ω << st |k| (10) принимает вид [10]:
(11)
Будем полагать, что на границе раздела между полуограниченными магнитной (y < 0) и немагнитной (y > 0) средами выполнена следующая система граничных условий (Ly ≡ ∂ly / ∂y):
(12)
Пусть q││OY, а b – единичный вектор вдоль линии пересечения плоскости границы раздела сред (y = 0) и плоскости падения объемной упругой волны в идеальной жидкой среде (6). Если – угол между равновесной ориентацией вектора антиферромагнетизма (l0 ⊥ q) и нормалью к плоскости падения волны a (a = [bq]), то в АФМ-среде (7) для краевой задачи (12) с помощью (9), (10) можно получить характеристическое уравнение: полином четвертой степени относительно η2 ≡ −(kq)2 (), причем в условиях ПВО η2 > 0. Если его корни обозначить как , то для АФМ-среды (7) (, ) пространственную структуру вектора упругих смещений u и трех компонент тензора упругих напряжений , определяемых как , при y < 0 можно представить в виде ():
(13)
Здесь Aj и U(j)(σ(j)) – соответственно амплитуда и векторы поляризации, отвечающие парциальной волне ej. Так как мы рассматриваем условия ПВО для падающей из жидкости (6) на поверхность АФМ (7) плоской объемной акустической волны, то, как показывает расчет, при ϑ → 0 в (13)
(14)
По аналогии с методикой расчета из [5], удобно, с помощью граничных условий (12), исключить из дальнейшего рассмотрения три из четырех амплитуд парциальных волн. Если это А2−4, то для формирующейся в АФМ-среде (7), (11)–(14), четырехпарциальной эванесцентной МУ-волны пространственную структуру коллинеарных направлению q компонент векторов u и можно представить в виде
(15)
Таким образом, если в немагнитной среде (6) на границу раздела (12) падает плоская объемная упругая волна, то в условиях ПВО, с учетом введенных выше обозначений, ее коэффициент отражения от поверхности полуограниченной АФМ-среды (7) можно представить в виде (отметим, что, в отличие от [5], здесь мы используем определение поверхностного волнового импеданса, принятое в кристаллоакустике [11]):
(16)
Расчет показывает, что при ϑ → 0 в (15), (16) одновременно (|F31 + F41|) → 0, а соотношение для входного импеданса в (16) принимает вид
(17)
где
(18)
а явный вид соотношений для fA(ω,h), φA(ω,h) не приводится ввиду их громоздкости.
В результате из (16)–(18) следует, что у поверхностного волнового импеданса так же, как и в (2), и числитель и знаменатель могут независимо стремиться к нулю, если одновременно с ϑ → 0, так же и . Таким образом, выполнение условия существования (3) в случае (17), (18) становится возможным, если a || l0 || OZ, b || OX (т.е. ϑ = 0), а . Согласно (18), этот предел – соотношение для спектра, возникающего в этом случае за счет коллективизации вкладов в полуограниченной АФМ-среде (7), (12), (14) магнитоупругого и неоднородного обменного взаимодействий, “темного” состояния: симметрийно защищенного поверхностного ССК (двухпарциальной эласто-обменной поверхностной волны). С учетом (14) в (18) дает
(19)
В эластостатическом пределе (ω << sth) соотношение для спектра симметрийно защищенного поверхностного ССК (19) принимает вид [10]
(20)
Из совместного анализа (19), (20) следует, что на плоскости внешних параметров “ω − h” необходимое условие существования такого поверхностного симметрийно защищенного ССК имеет вид
т.е. возможно, например, в случае низкотемпературных антиферромагнетиков [8], для которых st > c.
При a || l0 || OZ, b || OX и структура V(ϑ = 0) в (16) также значительно упрощается:
(21)
и отвечает коэффициенту отражения плоской объемной волны, падающей из жидкости (6) на поверхность полуограниченного упруго изотропного немагнитного твердого тела. При этом в (21) Zin = 0 () с учетом (14) соответствует спектру поверхностной волны Рэлея в полуограниченной упругоизотропной немагнитной среде [5].
Пусть для падающей извне на поверхность полуограниченного магнетика (7) квазиплоской объемной упругой волны ее частота и угол падения с учетом (17) связаны соотношением Zin (ϑ) = 0. В этом случае качественный вид зависимости максимума продольного смещения отраженного остронаправленного пучка ∆S ≡ ∆(h = hN(ω), ϑ) от управляемого углом ϑ “фактора вытекания” [12] распространяющейся в магнетике несобственной поверхностной МУ-волны в безразмерных переменных графически представлен на рис. 1. Если, следуя [5, 13], считать, что жидкая среда, из которой падает на поверхность магнетика остронаправленный пучок, – ксилол и ширина пучка мм, а частота падающей волны ω ≈ ω0 (в (18)–(20) энергия активации магнона в легкоплоскостном АФМ ω0 ≈ 109 Гц), то при в (4) , а при , .
Рис. 1. Зависимость максимальной величины смещения Шоха ∆S ≡ ∆(h = hN(ω), ϑ) от степени близости угла падения квазиплоской волны к точке формирования поверхностного симметрийно защищенного ССК (3), (19), (20) (ϑ = 0).
Отметим, что для случая падения квазимонохроматической упругой волны, параметры ω и h которой связаны условием Zin(ϑ)=0, вид рис. 1 качественно не изменится, если в нем заменить на ῶ∆t, а на , где ῶ – несущая частота, а – фактор затухания сигнала, распространяющегося в магнетике.
Пусть теперь АФМ (7) находится в той же магнитоакустической конфигурации, но занимает не полупространство (y < 0), а слой (−d < y < d), обе поверхности которого с помощью граничных условий вида
(22)
акустически связаны с полупространствами (d < |y|), заполненными идеальной жидкостью. В (22) верхние индексы “±” относятся к верхнему и нижнему полупространству соответственно. Если эти две полуограниченные среды неэквивалентны между собой по своим физическим свойствам, то, как показывает расчет, для плоской объемной упругой волны, падающей из акустически менее плотной среды на поверхность магнитного слоя, отвечающий ϑ = 0 спектр симметрийно защищенных ССК, формирующийся в этом случае в условиях ПВО (), при одновременном учете МУ и неоднородного обменного взаимодействий с учетом (14) принимает вид (k ∊ XY, q || OY)
(23)
В частном случае η3,4d → ∞ (23) переходит в (19). В окрестности такого темного состояния (при ϑ → 0) остаются в силе соотношения (1)–(5) (включая одновременное стремление к нулю и числителя и знаменателя входного поверхностного волнового импеданса при приближении к точке формирования “темного” состояния). Однако теперь, как следует из (23), в рассматриваемой полуоткрытой магнитной сэндвич-структуре симметрийно защищенные ССК могут формироваться как антиволноводными (если в (23) ), так и волноводными (если в (23) ) модами спектра магнонных поляронов с u || OZ || l0 и k ∊ XY, распространяющихся вдоль слоя легкоплоскостного АФМ (7), у которого при y = ±d σzy = 0, Ly = 0. В этом случае для падающей извне на поверхность магнитного слоя плоской объемной волны структура коэффициента отражения будет совпадать с (21), но теперь Zin = 0, с учетом (14), как при , так и при , будет определять спектр обобщенной волны Лэмба [3, 14], которая сформирована магнонными поляронами с u,k ∊ XY, распространяющимися вдоль рассматриваемого слоя АФМ (7).
Таким образом, уже на уединенной границе раздела немагнитной и магнитной сред вследствие совместного учета магнитоупругого и неоднородного обменного взаимодействий возможно существование симметрийно защищенных поверхностных “темных” состояний в сплошном спектре излучения вытекающих магнонных поляронов. При их формировании, так же как и в случае поверхностных ССК интерференционного типа [3], и числитель и знаменатель входного поверхностного волнового импеданса отражающей среды одновременно стремятся к нулю (усиливается квазистационарность резонансно возбуждаемой вытекающей ОПВ и первого и второго типа). По мере приближения параметров падающей извне на поверхность магнетика квазиплоской (или квазимонохроматической) объемной упругой волны к тем, что отвечают указанному типу “темных” состояний, становится возможным практически неограниченное увеличение (в рамках рассмотренной модели) незеркальных эффектов отражения первого порядка. Следует, однако, учесть, что в данной работе, как и в предыдущей [3], и сам метод расчета, и сделанные на основе полученных результатов выводы относительно возможности реализации сверхусиления незеркальных эффектов отражения основаны на приближении квазиплоской волны, т.е. в первую очередь справедливы для падающих остронаправленных волновых пучков (ширина пучка () много больше длины волны (), отвечающей его осевому вектору) или импульсов с узким частотным спектром (несущая частота импульса много больше обратной длительности импульса). При форма огибающей отраженного пучка крайне незначительно изменяется по сравнению с формой падающего пучка. Однако по мере приближения параметров падающей квазиплоской волны к тем, что характерны для точки ССК (2), (3), по отношению к ширине пучка растет длина пробега резонансно возбуждаемой вытекающей ОПВ (), и она пространственно отделяется от других волновых компонент, формирующих отраженный пучок. В частности, как и в случае боковой волны (см. [13]), по отношению к геометро-акустической части отраженного пучка, резко увеличивается угловой размер области, в котором присутствует поле фононного излучения (область “засветки”) распространяющегося магнонного полярона (вытекающей ОПВ) при одновременном снижении интенсивности излучаемой на бесконечность плоской объемной упругой волны. Аналогично при резонансном импульсном возбуждении рассматриваемой вытекающей ОПВ генерируемый ею сигнал будет отделяться по времени от сигнала других гармоник, составляющих отраженный импульс, а величина этого запаздывания будет неограниченно (в рамках рассматриваемой модели) возрастать по мере приближения на плоскости “ω − h” к точке формирования симметрийно защищенного поверхностного ССК. Наконец, если, в отличие от (5), фронт падающего импульса не является плоским, то в условиях (2), (3) появляется и дополнительный вклад в усиление эффекта Вигнера (5), связанный с одновременным неограниченным усилением также и пространственного эффекта Шоха (4) [15]. Для более корректного учета указанной выше возможности неограниченного роста незеркальных эффектов фононного отражения первого порядка в непосредственной окрестности точки формирования поверхностного ССК используемая в нашей работе схема вычислений должна быть пересмотрена, например, в пользу расчетов на основе конкретизации формы падающего волнового пучка подобно [12, 13]. Однако это предполагается сделать в отдельной работе.
ИСТОЧНИК ФИНАНСИРОВАНИЯ
Работа поддержана Российским научным фондом (проект 20-19-00745-П).
КОНФЛИКТ ИНТЕРЕСОВ
Авторы данной работы заявляют, что у них нет конфликта интересов.
About the authors
Yu. V. Gulyaev
Institute of Radio Engineering and Electronics named after Kotelnikov, Russian Academy of Sciences
Email: shavrov@cplire.ru
Academician of the RAS
Russian Federation, MoscowO. S. Sukhorukova
Donetsk State University; Donetsk Institute for Physics and Engineering named after Galkin
Email: shavrov@cplire.ru
Russian Federation, Donetsk; Donetsk
A. S. Tarasenko
Donetsk Institute for Physics and Engineering named after Galkin
Email: shavrov@cplire.ru
Russian Federation, Donetsk
S. V. Tarasenko
Donetsk Institute for Physics and Engineering named after Galkin
Email: shavrov@cplire.ru
Russian Federation, Donetsk
V. G. Shavrov
Institute of Radio Engineering and Electronics named after Kotelnikov, Russian Academy of Sciences
Author for correspondence.
Email: shavrov@cplire.ru
Russian Federation, Moscow
References
- Maekawa S., Kikkawa T., Chudo H., Ieda J., Saitoh E. Spin and spin current – From fundamentals to recent progress Journal of Applied Physics. 2023. V. 133(2). P. 020902.
- Пекар С.И. Исследования по электронной теории кристаллов. М.–Л.: ГТТЛ, 1951. 256 с.
- Гуляев Ю.В., Сухорукова О.С., Тарасенко А.С., Тарасенко С.В., Шавров В.Г. “Суперрезонансные” состояния в спектре вытекающих поверхностных магнонных поляронов. // ДАН. 2022. Т. 505. № 1. С. 10–15.
- Azzam S.I., Kildishev A.V. Photonic Bound States in the Continuum: from Basics to Applications // Adv. Opt. Mater. 2021. V. 9. P. 2001469.
- Бреховских Л.М. Волны в слоистых средах. М.: Изд-во АН СССР, 1957. 502 c.
- Favretto-Cristini N., Komatitsch D., Carcione J.M., Cavallini F. Elastic surface waves in crystals. Pt 1: Review of the physics // Ultrasonics. 2011. V. 51(6). P. 653–660.
- Гуляев Ю.В., Тарасенко С.В., Шавров В.Г. Электромагнитный аналог вытекающей поверхностной упругой волны первого типа для уединенной границы раздела прозрачных диэлектриков // УФН. 2020. Т. 190. № 9. С. 933–949.
- Ожогин В.И., Преображенский В.Л. Ангармонизм смешанных мод и гигантская акустическая нелинейность антиферромагнетиков // УФН. 1988. Т. 155. № 4. С. 593–621.
- Такер Дж., Рэмптон В. Гиперзвук в физике твердого тела. М.: Мир, 1975. 453 c.
- Гуляев Ю.В., Тарасенко С.В., Шавров В.Г. Спин-волновая акустика антиферромагнитных структур как магнитоакустических метаматериалов // УФН. 2011. Т. 181. № 6. С. 595–626.
- Балакирев М.К., Гилинский И.А. Волны в пьезокристаллах. Новосибирск: Наука, 1982. 240 с.
- Bertoni H.L., Tamir T. Unified theory of Rayleigh-angle phenomena for acoustic beams at liquid-solid interfaces // Appl. Phys. 1973. V. 2. P. 157–172.
- Бреховских Л.М., Годин О.А. Акустика слоистых сред. М.: Наука, 1989. 414 с.
- Auld B.A. Acoustic Fields and Waves in Solids. N.Y.: J. Willey, 1973. 423 p.
- Chauvat D., Emile O., Bretenaker F., A. Le Floch Direct Measurement of the Wigner Delay Associated with the Goos–Hänchen Effect Phys. Rev. Lett. // 2000. V. 84. P. 71–74.
Supplementary files
