УДК 537.9
ОБРАТНЫЙ МАГНИТОЭЛЕКТРИЧЕСКИЙ ЭФФЕКТ В СЛОИСТЫХ СТРУКТУРАХ ПЬЕЗОЭЛЕКТРИК—ФЕРРОМАГНЕТИК И ПЬЕЗОЭЛЕКТРИК—ФЕРРОМАГНЕТИК— ПЬЕЗОЭЛЕКТРИК В ФОРМЕ ПРЯМОУГОЛЬНОЙ ПЛАСТИНКИ
Д.А.Филиппов, Г.С.Радченко*, В.М.Лалетин**, Т.О.Фирсова, Т.А.Галкина
INVERSE MAGNETOELECTRIC EFFECT IN LAYERED STRUCTURES PIEZOELECTRIC— FERROMAGNETIC AND PIEZOELECTRIC—FERROMAGNETIC—PIEZOELECTRIC IN THE SHAPE OF A RECTANGULAR PLATE
D.A.Filippov, G.S.Radchenko*, V.M.Laletin**, T.O.Firsova, T.A.Galkina
Политехнический институт НовГУ, Dmitry.Filippov@novsu.ru *Научно-исследовательский институт физики Южного федерального университета, Ростов-на-Дону **Институт технической акустики, Беларусь, laletin57@rambler.ru
Представлена теория обратного магнитоэлектрического эффекта в слоистых структурах пьезоэлектрик— ферромагнетик и пьезоэлектрик—ферромагнетик—пьезоэлектрик в форме прямоугольной пластинки, основанная на совместном решении уравнений эластодинамики и электростатики отдельно для магнитострикционной и пьезоэлектрической фаз с учетом условий на границе раздела между фазами. Получены выражения для коэффициента обратного магнитоэлектрического преобразования через параметры, характеризующие магнитострикционную и пьезоэлектрическую фазы. Рассчитаны теоретические зависимости коэффициента обратного магнитоэлектрического преобразования от частоты переменного электрического. Результаты расчетов хорошо согласуются с экспериментальными данными. Ключевые слова: обратный магнитоэлектрический эффект, магнитострикция, пьезоэлектричество, электромеханический резонанс
The theory of inverse magnetoelectric effect in layered structures, based on a joint solution of equations of electrostatics and elastodynamics for magnetostrictive and piezoelectric phases taking into account the boundary condition at the interface is presented. The equations for frequency dependence of the inverse ME coefficient in the region of electromechanical resonance was obtained. The dependence of the inverse ME coefficient on the frequency of the alternating electric field were calculated. The calculation results are in good agreement with the experimental data.
Keywords: inverse magnetoelectric effect, magnetostriction, piezoelectricity, electromechanical resonance
1. Введение
Обратный, или инверсный, магнитоэлектрический эффект интересен тем, что на его основе можно создавать трансформаторы, имеющие от классических всего одну обмотку [1,2]. В отличие от автотрансформатора, который также имеет одну обмотку, в магнитоэлектрическом трансформаторе входная и выходная цепи гальванически развязаны, что создает большие преимущества перед автотрансформатором, который не имеет гальванической развязки между входной и выходной цепью.
Магнитоэлектрический (МЭ) эффект, лежащий в основе рабаты МЭ трансформатора, теоретически был предсказан [3,4] и экспериментально обнаружен [5,6] более полувека назад. Уникальность этого эффекта заключается том, что он относится к так называемым перекрестным эффектам и связывает между собой величины, имеющие разные тензорные размерности — аксиальный вектор с полярным и, наоборот, полярный с аксиальным. Прямой МЭ эффект заключается в индуцировании электрического напряжения на обкладках конденсатора, диэлектриком которого является МЭ материал. При обратном эффекте происходит изменение намагниченности МЭ материала при помещении его в электрическое поле. Обратный МЭ эффект интересен тем, что позволяет управлять магнитными свойствами вещества электрическим полем, что значительно упрощает конструкцию управляющих устройств, заменяя катушки индуктивности обкладками конденсаторов. Например, в СВЧ области спектра в МЭ кристаллах можно при помощи электрического поля управлять сдвигом линии магнитного резонанса [7-9], что позволяет создавать на этой основе такие приборы, как СВЧ- и фазовращатели, фильтры, аттенюаторы и т.д. Несмотря на многочисленные исследования, выполненные на монокристаллах, широкого применения на практике МЭ эффект не нашел. Это объясняется малостью величины эффекта и тем, что в большинстве материалов МЭ эффект наблюдается при температурах значительно ниже комнатной. Малость МЭ эффекта в монокристаллах обусловлена тем, что механизмом его возникновения является совместное действие нечетной части потенциала внутрикри-сталлического поля, внешнего электрического поля и спин-орбитального взаимодействия [10]. Спин-орбитальное взаимодействие является релятивистским эффектом, чем и объясняется малость МЭ эффекта в монокристаллах.
В композиционных МЭ материалах, состоящих из механически связанных между собой магнитост-
рикционной и пьезоэлектрической фаз, величина эффекта на несколько порядков больше и наблюдается он при комнатной температуре. Механизмом возникновения МЭ эффекта в композиционных материалах является механическое взаимодействие магнитост-рикционной и пьезоэлектрической подсистем. При прямом эффекте в переменном магнитном поле в магнитострикционной компоненте возникают деформации, которые посредством механического взаимодействия передаются в пьезоэлектрическую фазу, в которой вследствие пьезоэффекта происходит изменение поляризации, что приводит к возникновению напряжения на обкладках конденсатора. При обратном МЭ эффекте в результате подачи электрического напряжения на обкладки конденсатора в пьезоэлектрической компоненте возникают деформации, которые передаются в магнитострикционную фазу, где вследствие магнитострикции происходит изменение намагниченности, что приводит к возникновению ЭДС индукции на катушке, внутри которой находится образец.
При теоретическом описании МЭ эффекта используются либо метод эффективных параметров [11,12], либо метод, основанный на решении уравнений эластодинамики и электростатики отдельно для пьезоэлектрической и магнитострикционной фаз при учете условий на границе раздела между фазами [13,14]. Метод эффективных параметров применим тогда, когда характерные размеры композита много меньше длины звуковой волны. Он хорошо подходит для описания МЭ эффекта в объемных композитах. Описание обратного МЭ эффекта на основе этого метода выполнено в работах [15,16]. Однако для слоистых магнитострикционно-пьезоэлектрических структур, когда длина звуковой волны соизмерима с толщиной слоев, этот метод неприменим. Ранее в работах [17-19] была представлена теория обратного МЭ эффекта, основанная на решении уравнений эла-стодинамики и электростатики отдельно для магни-тострикционной и пьезоэлектрической фаз с учетом условий на границе раздела. Однако в данных работах предполагалось, что смещения и деформации магнитострикционного и пьезоэлектрического слоев одинаковы и граница раздела между слоями никак не учитывалась. При обратном МЭ эффекте под действием электрического поля колебания среды сначала возбуждаются в пьезоэлектрическом слое. Затем посредством тангенциальных напряжений эти колебания передаются через границу раздела в магнетик, где вследствие магнитострикции происходит изменение намагниченности. При распространении колебаний по толщине образца, вследствие того, что
магнитная фаза является пассивной по отношению к электрическому полю, которое возбуждает колебания в пьезоэлектрике, происходит изменение амплитуды колебаний среды, в результате чего деформации и напряжения в магнетике и пьезоэлектрике будут неоднородными по толщине образца. В данной работе этот механизм возбуждения колебаний учтен в явном виде путем введения граничных условий на интерфейсе между фазами, что делает модель более точной.
2. Обратный эффект в двухслойной пластине
При обратном эффекте в качестве объекта исследований используют образец из магнитоэлектрического материала, помещенного внутрь катушки с числом витков N. Схематичное изображение такой структуры приведено на рис. 1.
Рис.1. Схематичное изображение структуры: 1 — магнетик, 2 — пьезоэлектрик, 3 — электроды, 4 — непроводящий каркас, 5 — обмотка катушки
Для теоретического описания эффекта в качестве модели рассмотрим образец в виде прямоугольной пластинки длиной L и шириной W, представляющий собой двухслойную магнитострикци-онно-пьезо-электрическую структуру с толщиной магнетика tm и пьезоэлектрика ^ (рис.2). На верхней и нижней поверхности пластинки нанесены тонкие металлические электроды. Образец предварительно поляризован в направлении, перпендикулярно слоям пластинки. Будем считать пластинку узкой, т.е. L>>W. Это позволяет ограничиться рассмотрением колебаний, распространяющихся только по длине пластинки, вследствие того, что колебания, распространяющиеся по ширине пластинки, будут более высокочастотными. Начало системы координат совместим с границей раздела слоев, а ось 2 направим вертикально вверх, перпендикулярно слоям. Ограничимся рассмотрением поперечного эффекта вследствие того, что он гораздо легче реализуется в эксперименте, чем продольный. В этом случае поле подмагничивания Ны^, направлено по оси X перпендикулярно вектору поляризации Р. Рассмотрим пла-нарные колебания, которые распространяются вдоль оси X. Переменное электрическое поле с частотой ю возбуждает в пьезоэлектрической подсистеме упругие колебания, которые посредством касательных напряжений передаются через границу раздела в
магнитострикционный слой, в результате чего происходит изменение намагниченности.
Рис.2. Модель двухслойной структуры: 1 — магнетик, 2 — пьезоэлектрик, 3 — электроды
Поскольку пластинка узкая, то отличными от нуля компонентами тензора напряжений будут только 1хх и 1хг , где индекс а = р для пьезоэлектрика и а = т для магнетика. При данной ориентации полей уравнения тензора деформаций для пьезоэлектрической и магнитострикционной фаз аSij и х — проекция вектора магнитной индукции тВХ имеют следующий вид:
=—РТ + И Е
°хх ¿хх ' "-ХХ^^П
YР
1
Р V = РТ
/-1 1х1 °р
(1)
(2)
где рТхх и рТх2 — компоненты тензора напряжений в пьезоэлектрической фазе, Ур,Ор — модуль Юнга и модуль сдвига пьезоэлектрика, Ихх,1 — пьезоэлектрический тензор, Ег — г-проекция вектора напряженности электрического поля.
Для магнитострикционной фазы аналогичные уравнения запишутся в следующей форме:
т V =1 тт + пв в
I т хх,х
тV = 1 тт ¿хг—^ хг,
т
вх = Цхх{Нх + ^ х ттхх)
(3)
(4)
(5)
тт
где 1хх и 1хг — компоненты тензора напряжений в магнитострикционной фазе, 1т, От — модуль Юнга и
в дА
модуль сдвига магнетика, д =дв — пьезомагнит-
ный коэффициент, Нх — напряженность магнитного поля, А — магнитострикция магнетика.
Уравнение движения для смещений среды имеет вид:
даТг:1
Ра
",2 а ^агр
д их д Тхх
дt2
дх
дг
(6)
где ра — плотность магнетика или пьезоэлектрика,
агт
1ц — тензор напряжений.
+
Для расчета обратного МЭ эффекта воспользуемся методом, разработанным в [14]. Решение уравнения для вектора смещения среды (6) представим в виде плоских волн, амплитуда которых изменяется по толщине образца:
а и(х, z)=аg(z)(A1 ^(и^ - Ъх)+А2 sm(юt - Ъх)), (7) где А1 и А2 — постоянные интегрирования, аg(z) — некоторая функция. Для нахождения этой функции подставим выражение (7) в уравнение (6), в результате чего получим два дифференциальных уравнения для функций pg(z) и mg(z) с четырьмя неизвестными постоянными интегрирования, которые определяются с помощью граничных условий. На границе раздела в точке г = 0 равны смещения сред и также равны тангенциальные составляющие тензора напряжений, т.е. имеют место следующие равенства:
ри(х,0) =ти(х,0); (х,0)=тТхг(х,0). (8) На свободных верхней и нижней поверхностях образца тангенциальные компоненты тензора напряжений равны нулю, т.е. имеют место следующие граничные условия:
рГХ2 (х, г = ^) = 0 и тТХ2 (х, г = -т) = 0.
(9)
Используя граничные условия (8,9), получим систему уравнений, условие совместности которой приводит к следующему уравнению:
Ш(Кт ) = Ур1 р tg(к р ), (10)
где кт =xmtm и кр =хр1р — безразмерные перемен-
ные. Здесь введены обозначения %т =-2(1+у)
"4 - ъ2
V
Гр = 2(1+ У)
4 - ъ 2
Vm = №, Vp = ^ — скоро-V Рт \ Рр
сти продольных волн в магнетике и пьезоэлектрике соответственно.
Выполняя вычисления, аналогичные проведенным в работе [14], для компонент тензора напряжений в магнитострикционной фазе получим уравнения в виде:
тТхх = ^(Ъх)(ехр(-2Кт) X
X ехр(ХтГ) + ехр(-ХтГ)) - g1¡L,xBx 1 (11)
У^В + dxx,xEz
Ъ COS(к)(1 + ехр(-2Кт)) У :1Ь(Кт) +tg(Кp)
где А2 =
введены безразмерные параметры к=ЪЬ/2, у=У^т/Ур?р.
Коэффициент обратного МЭ преобразования структуры определим как отношение среднего значения индукции магнитного поля в образце к величине приложенного внешнего электрического поля, ее вызвавшей, т.е.:
т ь/2
< Вх > =_
< Ег > Ь • Гт
^ ^Вх (х, z)dxdz
1 0 -Ь/2
Ег
(12)
Подставляя выражение для тензора напряжений (13) в уравнение для магнитной индукции (5) и выполняя интегрирование с учетом закона полного тока и условия разомкнутой цепи, для коэффициента обратного МЭ преобразования получим выражение в виде:
аВ = УтУ-хх^^хх^ -д—. (13)
Здесь введены следующие обозначения:
д= 1и(к) шсо/лскт)+гёС^р) ^ (14)
к кт V Кт К р )
Да =1 + К2т(1 -УД), К2т = УтУхх(ч1,х )2. (15)
Как следует из (13), при значении параметра Да = 0 имеет место резонансное увеличение эффекта.
В области низких частот обратный МЭ коэффициент практически не зависит от частоты. В предельном случае низких частот, разлагая в ряд уравнение (13) по малым параметрам к , кт, кр и ограничиваясь первыми членами разложения, для его низкочастотного значения получим выражение в следующем виде:
аВ = УтУххЧВхх/хх,2 /(1 +У+Кт)- (16)
Выходное напряжение, индуцируемое на катушке, найдем, используя закон Фарадея для электромагнитной индукции. Полагая, что все магнитное поле сосредоточено в магнетике, для индуцированного напряжения получим выражение в виде:
>т
'Тр
где Vп — входное напряжение, подаваемое на обкладки структуры, Vinv — выходное напряжение, снимаемое с катушки индуктивности.
Как следует из (17), коэффициент передачи напряжения определяется выражением
=юа в—ФШ*
(17)
ЪТ = юа В*тЖЫ. 'р
(18)
Его зависимость от частоты также имеет резонансный характер, а численное значение определяется не только характеристиками образца, но и числом витков катушки.
3. Сравнение с экспериментом
В области низких частот индуцируемое вследствие обратного МЭ эффекта напряжение на катушке и коэффициент передачи напряжения линейно
возрастает с частотой, но по мере приближения к области электромеханического резонанса зависимость становится нелинейной. На рис.3 представлена рассчитанная по формуле (17) зависимость индуцированного напряжения Vinv от частоты приложенного электрического поля для структуры Р2Т—№—Р2Т в области электромеханического резонанса, а на рис.4 — в низкочастотной области спектра. При расчетах использовались следующие параметры: длина Ь = 9,5 мм, ширина Ж = 4,5 мм, толщина каждого слоя пьезо-
и
т
р
В
электрика tp = 0,5 мм, толщина никеля tm = 0,2 мм, модули Юнга Ym = 210 ГПа и Yp = 6 ГПа для Ni и PZT соответственно, pm = 8900 кг/м3 и pm = 7000 кг/м3 их плотности. Значение пьезомодуля принималось равным p2XX = -150 м/В, величина пьезомагнитного ко-
эффициента q [20].
определялась из результатов работы
Рис.3. Зависимость индуцированного напряжения и^ от частоты приложенного электрического поля для структуры PZT—N—PZT в области электромеханического резонанса
2.0-г
>
н
.с О
Рис.4. Зависимость индуцированного напряжения Unv от частоты приложенного электрического поля для структуры PZT—Ni—PZT в низкочастотной области
Как видно из рисунка, в области низких частот наблюдается линейная зависимость Vinv от частоты, а
на частоте порядка 180 кГц наблюдается резонансное увеличение эффекта. Рассчитанные теоретические частотные зависимости Vinv хорошо согласуются с экспериментальными данными, представленными в работе [17].
На рис.5 представлена рассчитанная частотная зависимость коэффициента обратного МЭ преобразования ав от частоты для структуры D-терфенол—
ЦТС. При расчетах использовались параметры структуры, представленные в работе [21]. Параметры D-терфенола, используемые в расчетах, брались из работы [22].
Рис.5. Частотная зависимость коэффициента обратного МЭ преобразования ав от частоты для структуры D-терфенол— ЦТС
Как видно из рис.5, теоретическая зависимость находится в хорошем согласовании с экспериментальной, данной в работе [21].
Заключение
Таким образом, как следует из представленных выше выражений для величин, характеризующих обратный МЭ эффект, в области низких частот коэффициент обратного МЭ преобразования не зависит от частоты, в то время как индуцированное на катушке индуктивности напряжение и коэффициент передачи напряжения линейно возрастают с увеличением частоты. По мере приближения к области электромеханического резонанса зависимость указанных величин становится нелинейной и на частоте электромеханического резонанса наблюдается пиковое увеличение эффекта. Частота электромеханического резонанса зависит не только от длины образца, но также от толщины магнитострикционного и пьезоэлектрического слоев. Величина коэффициента обратного МЭ преобразования нелинейным образом зависит от толщины магнитострикционного и пьезоэлектрического слоев, причем кривая имеет максимум, положение которого зависит от характеристик магнитострикционного и пьезоэлектрического слоев. Представленные в работе теоретические зависимости хорошо согласуются с экспериментальными данными.
Работа выполнена при частичной поддержке Министерства образования и науки РФ в рамках проектной части государственного задания проект №я11.177.2014/К и гранта конкурса РФФИ — БРФФИ проект №216-52-00184 Бел_а.
1. Филиппов Д.А., Галкина Т.А., Лалетин В.М., Srinivasan G. Преобразователь напряжения на основе инверсного магнитоэлектрического эффекта // Письма в ЖТФ. 2012. Т.38. №2. С.82-86.
2. Филиппов Д.А., Антоненков О.В., Радченко Г.С., Галкина Т.А. Кольцевой преобразователь напряжения на основе инверсного магнитоэлектрического эффекта // Вестник НовГУ. 2014. №81. С.61-63.
3. Ландау Л.Д., Лифшиц Е.М. Электродинамика сплошных сред. М.: ГИФМЛ, 1959. 532 с.
4. Дзялошинский И.Б. К вопросу о магнитоэлектрическом эффекте в антиферромагнетиках // ЖЭТФ. 1960. Т.37. С.881-882.
5. Астров Д.Н. Магнитоэлектрический эффект в окиси хрома // ЖЭТФ. 1961. Т.40. С.1035-1041.
6. Folen V.J., Rado G.T., Stalder E.W. Anysotropy of the mag-netoelectric effect in &2O3 // Phys. Rev. Lett. 1961. V.6. №11. P.607-608.
7. Kita E., Siratori K., Tasaki A.J. Electronic shift in the anti-ferromagnetic resonance and the mechanism of the parallel magnetoelectric effect of Cr2O3 // J. Phys. Soc. Jap. 1979. Vol.46. №3. P.1033-1034.
8. Бичурин М.И., Петров В.М. Влияние электрического поля на спектр антиферромагнитного резонанса в борате железа // ФТТ. 1987. Т.29. №8. С.2509-2510.
9. Zhou H.M., Chen O., Deng J.-H. Ferromagnetic resonance frequency shift model of laminated magnetoelectric structure tuned by electric field // Chinese Physics B. 2014. V.23. №04752.
10. Bichurin M.I., Filippov D.A. The microscopic mechanism of the magnetoelectric effect in the microwave range // Ferroelectric. 1997. V.204. №1-4. P.225-232.
11. Harshe G., Dougherty J.O., Newnham R.E. Theoretical modelling of multilayer magnetoelectric composites // Int. J. Appl. Electromagn. Mater. 1993. V.4. P.145-159.
12. Bichurin M.I., Petrov V.M., Srinivasan G. Theory of low-frequency magnetoelectric coupling in magnetostrictive-piezoelectric bilayers // Phys. Rev. B. 2003. V.68. P.054402 (113).
13. Филиппов Д.А. Теория магнитоэлектрического эффекта в гетерогенных структурах на основе ферромагнетик-пьезоэлектрик // ФТТ. 2005. Т.47. Вып.6. С.1082-1084.
14. Филиппов Д.А., Лалетин В.М., Galichyan T.A. Магнитоэлектрический эффект в двухслойной магнитострикци-онно-пьезоэлектрической структуре // ФТТ. 2013. Т.55. Вып.9. С.1728-1733.
15. Филиппов Д.А., Галкина Т.А., Srinivasan G. Инверсный магнитоэлектрический эффект в феррит-пьезоэлектрических структурах // Письма в ЖТФ. 2010. Т.36. Вып. 21. С.23-28.
16. Филиппов Д.А., Галкина Т.А., Лалетин В.М., Srinivasan G. Инверсный магнитоэлектрический эффект в дискообразных образцах из феррит-пьезоэлектрических композитов // ФТТ. 2011. Т.53. №9. С.1737-1742.
17. Fetisov Y.K., Petrov V.M., Srinivasan G. Inverse Magneto-electric Effects in a Ferromagnetic-Piezoelectric Layered Structure // J. Mater. Res. 2007. Vol.22. №8. P.2074-2080.
18. Tong B., Yang X.F., Ouyang J. et al. Giant converse magnetoelectric effect of AlN-(Fe90Co 10)78Si12B10 thin film composites // Journal of Alloys and Compounds. 2013. V.563. P.51-54.
19. Bichurin M.I., Petrov V.M., Petrov R.V. Direct and inverse magnetoelectric effect in layered composites in electromechanical resonance range: A review // JMMM. 2012. V.324. P.3548-3549.
20. Burdin D.A., Chashin D. V., Ekonomov N. A. et al. Nonlinear magnetoelectric effects in planar ferromagnetic-piezoelectric structures // JMMM. 2014. V.358-359. P.98-104.
21. Калгин А.В., Гриднев С.А., Gribe Z.H. Особенности обратного магнитоэлектрического эффекта в двухслойных композитах Tb0.12Dy0.2Fe0.68-PbZr0.53Ti0.4703 // ФТТ. 2014. Т.56. Вып.7. С.1278-1283.
22. Lim S.H., Kim S.R., Kang S.Y. et al. Magnetostrictive properties of polymer-bonded Terfenol-D composites // JMM. 1999. V.191. P. 113-121.
References
1. Filippov D.A., Galkina T.A.,. Laletin V.M, Srinivasan G. Preobrazovatel' napriazheniia na osnove inversnogo magnitoelektricheskogo effekta [Voltage transformer based on inverse magnetoelectric effect]. Pis'ma v ZhTF -Technical Physics Letters, 2012, no. 1, vol. 38, pp. 93-95.
2. Filippov D.A., Antonenkov O.V., Radchenko G.S., Galkina T.A. Kol'tsevoi preobrazovatel' napriazheniia na osnove inversnogo magnitoelektricheskogo effekta [Ring voltage converter based on inverse magnetoelectric effect]. Vestnik NovGU. Ser. Tekhnicheskie nauki - Vestnik NovSU. Issue: Engineering Sciences, 2014, no. 81, pp. 61-63.
3. Landau L.D., Lifshitz E.M. Elektrodinamika sploshnykh sred [Electrodynamics of continuous media]. Moscow, GIFML Publ., 1959. 532 p.
4. Dzialoshinskii I.B. K voprosu o magnitoelektricheskom effekte v antiferromagnetikakh [On the magneto-electrical effect in antiferromagnets]. Zhurnal eksperimental'noi i teo-reticheskoi fiziki (ZhETF) - Journal of Experimental and Theoretical Physics (JETP), 1960, vol. 10, no. 3, p. 628629.
5. Astrov D.N. Magnitoelektricheskii effekt v okisi khroma [Magnetoelectric effect in chromium oxide]. Zhurnal eksperimental'noi i teoreticheskoi fiziki (ZhETF) - Journal of Experimental and Theoretical Physics (JETP), 1961, vol. 13, vo. 4, pp. 729-733.
6. Folen V.J., Rado G.T., Stalder E.W. Anysotropy of the magnetoelectric effect in Cr2O3. Physical Review Letters, 1961, vol. 6, no. 11, pp. 607-608.
7. Kita E., Siratori K., Tasaki A.J. Electronic shift in the antiferromagnetic resonance and the mechanism of the parallel magnetoelectric effect of Cr2O3. Journal of the Physical Society of Japan, 1979, vol. 46, no. 3, pp. 1033-1034.
8. Bichurin M.I., Petrov V.M. Vliianie elektricheskogo polia na spektr antiferromagnitnogo rezonansa v borate zheleza [Effect of an electric field on the spectrum of antiferromagnetic resonance in ferrum borate]. Fizika tverdogo tela - Physics of the Solid State, 1987, vol. 29, no. 8, pp. 2509-2510.
9. Zhou H.-M., Chen Q., Deng J.-H. Ferromagnetic resonance frequency shift model of laminated magnetoelectric structure tuned by electric field. Chinese Physics B, 2014, vol. 23, no. 04752.
10. Bichurin M.I., Filippov D.A. The microscopic mechanism of the magnetoelectric effect in the microwave range. Ferroelectric, 1997, vol. 204, no. 1-4, pp. 225-232.
11. Harshe G., Dougherty J.P., Newnham R.E. Theoretical modelling of multilayer magnetoelectric composites. Int. J. Appl. Electromagn. Mater., 1993, vol. 4, pp. 145-159.
12. Bichurin M.I., Petrov V.M., Srinivasan G. Theory of low-frequency magnetoelectric coupling in magnetostrictive-piezoelectric bilayers. Physical Review B, 2003, vol. 68, p. 054402 (1-13).
13. Filippov D.A. Teoriia magnitoelektricheskogo effekta v geterogennykh strukturakh na osnove ferromagnetik-p'ezoelektrik [Theory of magnetoelectric effect in ferromagnetic-piezoelectric heterostructures]. FTT - Phys. Solid State, 2005, vol. 47, no. 6, pp. 1118-1121.
14. Filippov D.A., Laletin V.M., Galichyan T.A. Magnitoelektricheskii effect v dvukhsloinoi magni-tostriktsionno-p'ezoelektricheskoi strukture [Magnetoelectric effect in a magnetostrictive-piezoelectric bilayer structure]. Fizika tverdogo tela - Physics of the Solid State, 2013, vol. 55, no. 9, pp. 1840-1845.
15. Filippov D.A., Galkina T.A., Srinivasan G. Inversnyi magnitoelektricheskii effekt v ferrit-p'ezoelektricheskikh strukturakh [Inverse magnetoelectric effect in serrite-piezoelectric structures]. Pis'ma v ZhTF - Technical Physics Letters, 2010, vol. 36, no. 11, pp. 984-986.
16. Filippov D.A., Galkina T.A., Laletin V.M., Srinivasan G. Inversnyi magnitoelektricheskii effekt v diskoobraznykh obraztsakh iz ferrit-p'ezoelektricheskikh kompozitov [Inverse magnetoelectric effect in disk-shaped samples of
ferrite piezoelectric composites]. Fizika tverdogo tela - 20.
Physics of the Solid State, 2011, vol. 53, no. 9, pp. 18321838.
17. Fetisov Y.K., Petrov V.M., Srinivasan G. Inverse magnetoelectric effects in a ferromagnetic-piezoelectric layered structure. Journal of Materials Research, 2007, vol. 21. 22, no. 8, pp. 2074-2080.
18. Tong B., Yang X.F., Ouyang J., Lin G.Q., Chen S. Giant converse magnetoelectric effect of AlN-(Fe90Co 10)78Si12B10 thin film composites. Journal of Alloys and Compounds, 2013, vol. 563, pp. 51-54.
19. Bichurin M.I., Petrov V.M., Petrov R.V. Direct and inverse magnetoelectric effect in layered composites in 22 electromechanical resonance range: A review. Journal of Magnetism and Magnetic Materials, 2012, vol. 324, pp.35483549.
Burdin D.A., Chashin D.V., Ekonomov N.A., Fetisov L.Y., Fetisov Y.K., Srinivasan G. Sreenivasulu G. Nonlinear magnetoelectric effects in planar ferromagnetic-piezoelectric structures. Journal of Magnetism and Magnetic Materials, 2014, vol. 358-359, pp. 98-104.
Kalgin A.V., Gridnev S.A., Gribe Z.H. Osobennosti obratnogo magnitoelektricheskogo effekta v dvukhsloinykh kompozitakh Tb0.12Dy0.2Fe0.68-PbZr0.53Ti0.4703
[Specific features of the inverse magnetoelectric effect in two-layered Tb0.12Dy0.2Fe0.68-PbZr0.53Ti0.4703
composites]. Fizika tverdogo tela - Physics of the Solid State, 2014, vol. 56, no. 7, pp. 1327-1333. Lim S.H., Kim S.R., Kang S.Y., Park J.K., Nam J.T., Son D. Magnetostrictive properties of polymer-bonded Terfenol-D composites. Journal of Micromechanics and Microengineering, 1999, vol. 191, pp. 113-121.