CC BY
35
ТЕХНОЛОГИЧЕСКИЕ ОСНОВЫ ПОВЫШЕНИЯ НАДЕЖНОСТИ И КАЧЕСТВА ИЗДЕЛИЙ
УДК 538. 6
ОЦЕНКА АНТИРЕЗОНАНСА ПРИ ВОЗДЕЙСТВИИ ЛАЗЕРНОГО ИЗЛУЧЕНИЯ НА НАНОКОМПОЗИТНУЮ СТРУКТУРУ
Г. С. Макеева, О. А. Голованов, В. Я. Савицкий
Введение
Воспроизводимость результатов исследования новых наномасштабных свойств и функциональных возможностей фотонно-кристаллических устройств в терагерцевом диапазоне обусловливает необходимость единого подхода к их математическому моделированию, который позволит преодолеть ограничения, поставленные традиционными в оптике методами проектирования. Этот подход должен быть основан на идеях вычислительной электродинамики, успешно реализованных в системах автоматизированного моделирования и проектированиятехники СВЧ.
В этой связи актуальной является задача теоретического исследования дифракции электромагнитных волн и возникающих геометрических, размерных эффектов в анизотропных магнитных наноструктурах на основе периодических решеток ферромагнитных нанопроволок и магнитных фотонных кристаллов в терагерцевом диапазоне на математических моделях электродинамического уровня строгости.
Строгая постановка краевой задачи дифракции электромагнитных волн на решетке ферромагнитных нанопроволок заключается в следующем. Необходимо решить систему уравнений Максвелла
с электродинамическими граничными условиями совместно с уравнением Ландау-Лифшица движения вектора намагниченности в ферромагнетике с учетом обменного взаимодействия
1. Решение задачи дифракции методом автономных блоков с магнитными нановключениями и каналами Флоке
тогН(Г) = е0 е(Е + аЕ(О;
(1)
™ ч (В(^ гс* Е ^ ) =--— ;
сМ (г) Ж
-(( (г )Й эф (г)) + юг (( (г)-М (г)), (2)
где Е(г), Й(г) - векторы напряженности электрического и магнитного полей; М(/) - вектор намагниченности среды; В (г)- вектор магнитной индукции; Йэф(г) = Й(?) + Й(}(г) - суммарное
эффективное поле, включающее поле обменного взаимодействия Й (г)=q У2М(г), q - константа
обменного взаимодействия; £ - относительная диэлектрическая проницаемость среды; а - электропроводность среды; £0 - электрическая постоянная; у - гиромагнитное отношение; юг - частота релаксации; %0 - статическая восприимчивость.
Рассмотрим задачу дифракции плоской однородной электромагнитной волны на двумерной периодической решетке ферромагнитных нанопроволок, расположенных в среде с относительными диэлектрической и магнитной проницаемостями £у, цу (рис. 1). Падающая ТЕМ-волна (с волновым вектором к , амплитудой С1+1)(ю), частотой ю) распространяется поперечно (вдоль оси г) по отношению к направлению постоянного поля подмагничивания Й00 = Йоёу, приложенного
вдоль оси ферромагнитных нанопроволок.
Элементарную ячейку периодической решетки с геометрическими размерами а, Ь, с, содержащую ферромагнитную нанопроволоку радиуса г и длиной I, представим в виде АБ с виртуальными каналами Флоке и магнитным нановключением [1]. Ферромагнитные нанопроволоки находятся между входными сечениями 51, 52 автономного блока (АБ), рассматриваемого в качестве волноводного трансформатора.
а) б) в)
Рис. 1. Дифракция электромагнитной волны на периодической решетке магнитных нанопроволок: а - ориентация падающей ТЕМ-волны с волновым вектором к; б - двумерная решетка продольно намагниченных ферромагнитных нанопроволок; в - автономный блок (АБ) с каналами Флоке, содержащий магнитную нанопроволоку: 52 - входные сечения АБ, V - внутренняя область АБ, V - область магнитного нановключения
Трехмерная краевая задача дифракции для уравнений (1), (2) с электродинамическими граничными условиями и условиями неасимптотического излучения решена с помощью вычислительного алгоритма на основе метода автономных блоков с магнитными нановключениями и каналами Флоке.
Предлагаемым методом проведен электродинамический расчет ^-параметров - элементов ЛЦ' (а, Р - номера входных сечений 5а, а = 1,2; Р = 1,2) многомодовой многоканальной
матрицы рассеяния двумерной периодической решетки ферромагнитных нанопроволок в зависимости от их геометрических размеров и периода решетки, а также величины постоянного поля
подмагничивания в терагерцевом диапазоне частот. Результаты электродинамического расчета
модуля коэффициентов отражения
Я
и прохождения
Я
ТЕМ-волны через решетки ферро-
магнитных нанопроволок в зависимости от геометрических параметров нанопроволокии периода нанорешетки при постоянном поле подмагничивания Н0 = 75 Э на частоте / = 30 ТГц приведены на рис. 2.
Я„
я;
0,6
0,4
0,2
2 >
У ^ 4 \ V.
^ -----
10
20
30
40
50
60 2г, нм
Рис. 2. Коэффициенты отражения Я^1 и прохождения Я^2 ТЕМ-волны через магнитные нанорешетки в зависимости от диаметра нанопроволок 2г: Н0 = Н0у0 ; период решетки а = 3,5г, Ь = 1,251, с = 2г ;
длина нанопроволоки I = 300 нм ; кривые 1, 3 - модуль коэффициента отражения Я
кривые 2, 4 - модуль коэффициента прохождения Я
12
■ - Н0 = 75 Э ;--------Н0 = 0
Результаты математического моделирования дифракции ТЕМ-волны на периодической решетке магнитных нанопроволок (кривые 1, 2 рис. 2) получены при следующих параметрах ферромагнетика: намагниченность насыщения 4лМ0 = 21580 Гс, проводимость а = 1,03-105 Ом-1 -см"1, константа обменного взаимодействия д = 2,2 -10"9 Э - см2, параметр диссипации юг = 0,0023уН0 а = 0,0023. Ферромагнитные нанопроволоки расположены в немагнитной матрице - диэлектрической среде с относительной диэлектрической проницаемостью = 5 и относительной магнитной
проницаемостью среды ц = 1.
Для сравнения на рис. 2 приведены результаты электродинамического расчета модуля ко-
эффициентов отражения
Я
и прохождения
Я
(кривые 3,4 рис. 2) ТЕМ-волны через решетку
электрических нанопроволок (металлических нанопроволок в отсутствие поля подмагничивания
Н 0 = 0).
Как следует из результатов математического моделирования зависимости коэффициентов
отражения
Я
и прохождения
Я
ТЕМ-волны через периодическую решетку магнитных
нанопроволок имеют экстремумы (кривые 1, 2 рис. 2) при изменении геометрических размеров
магнитной нанорешетки. Для сравнения отметим, что коэффициент прохождения
Я
ТЕМ-волны через решетку электрических нанопроволок плавно уменьшается и его значение
Я
плавно стремится
стремится к нулю (кривая 4 рис. 2), а значение коэффициента отражения к единице (кривая 3 рис. 2) при увеличении диаметра нанопроволоки. Это объясняется тем, что
если диаметр 2г нанопроволок становится больше, чем толщина скин-слоя 5 = -
ю
в металле
(а = 1,03 -105 Ом 1 -см 1, ц = ц0 на частоте / = 30ТГц, § =30 нм), то при 2г >§ электромагнитная волна отражается от решетки электрических нанопроволок с коэффициентом отражения 1 (кривая 3 рис. 2 при 2г > 30 нм) как от эффективной «металлической пленки».
»11 »11
2. Зависимость S-параметров двумерной решетки магнитных нанопроволок от поля намагничивания вблизи точки антирезонанса
На электродинамическом уровне строгости с помощью вычислительного алгоритма на основе метода автономных блоков проведено численное исследование управляемости «^-параметрами матрицы рассеяния двумерной периодической решетки ферромагнитных нанопроволок с помощью внешнего магнитного поля в терагерцовом диапазоне частот.
12
Рассчитанные зависимости модуля коэффициента прохождения Лп ТЕМ-волны через периодическую решетку ферромагнитных (железо) нанопроволок от величины поля намагничивания Н0 при изменении геометрических размеров нанорешетки приведены на рис. 3. Как следует
из результатов математического моделирования, зависимости коэффициентов прохождения от приложенного внешнего постоянного магнитного поля Н0 имеют максимумы (кривые 1-4
»12 »11
рис. 3), при этом максимальное значение
»12 »11
и определяющее этот максимум значение поля
намагничивания Н 0 зависят также и от геометрических размеров магнитной нанорешетки (диа-
»12 »11
имеет мак-
метра 2г нанопроволоки периода решетки а = 3,5г). Коэффициент прохождения
симальное значение при Н0 = 800 Э для 2г = 10 нм. Однако оптимальное значение внешнего постоянного поля намагничивания Н 0, реализующее этот максимум, существенно уменьшается при увеличении диаметра нанопроволок и становится равным Н0 = 75 Э при 2г = 25 нм.
Рис. 3. Зависимость »12 от Н0 (/ = 30 ТГц; а = 3,5г, Ь = 1,25/, с = 2г ; / = 300 нм): 1 - 2г = 10 нм; 2 - 2г = 15 нм; 3 - 2г = 20 нм; 4 - 2г = 25 нм
Дадим следующую трактовку полученным результатам математического моделирования. Толщина скин-слоя в гиромагнитной проводящей среде отличается от толщины скин-слоя § в металле с М = 1 и может быть записана в виде
§1 = § , (3)
где М/ = |це/| + ц/ , М/ = Ме/ _ ^М"/ - эффективная скалярная магнитная проницаемость, которая зависит от компонент тензора магнитной проницаемости м .
Коэффициент затухания однородной плоской волны в гиромагнитной проводящей среде определяется
k' = ^^ (^Ф+ ^^ 1 ' ^ / 8 = 1/ §!• (4)
Толщина скин-слоя (3) зависит от действительной ц'^ и мнимой ц"е^ частей эффективной скалярной магнитной проницаемости цеф . Следовательно, значение 81 максимально при минимальных значениях ц'^ и , т.е. в точке антирезонанса эффективной магнитной проницаемости цеф . Явление антирезонансазаключается в обращении в нуль действительной части эффективной магнитной проницаемости ц'^ = 0. При этом в отличие от резонанса мнимая часть ц"^ (при наличии диссипациив ферромагнетике) имеет минимальное значение ц"^ ^ 0 . В неограниченном ферромагнетике условие антирезонанса имеет вид
mares = ЮИ + ЮМ , (5)
где юае - частота антирезонанса, юИ = уИ0, юм = у4яМ0.
Толщина скин-слоя в ферромагнитном металле 81 =8 / ^цф зависит от магнитной проницаемости, а следовательно, от ширины резонансной кривой, в которую вносит вклад и обменное взаимодействие.
Глубина проникновения поля в ферромагнитный металл 81 =8 / ^цф увеличивается вблизи точки антирезонанса (5) при цф = |цеф | + цф ^ 0, что приводит к уменьшению коэффициента затухания kэлектромагнитной волны. На частоте юае в пределах ширины линии ферромаг-81
нитного антирезонанса имеется минимум поглощения электромагнитной волны и резкое возрастание амплитуды волны, проходящей через эффективную «пленку» ферромагнитного металла
(толщиной 2г < 8—). При 8— =8/ ^цф > 2г коэффициент прохождения увеличивается.
Для решетки с толщиной 2г, сравнимой с глубиной проникновения поля 81 в ферромагнитный металл, расчет просачивающейся мощности требует учета электродинамических граничных условий, и его следует проводить обязательно с учетом обменного взаимодействия.
Для металлических образцов с ограниченными размерами в форме эллипсоида (предельными случаями являются бесконечно тонкая пластинка и бесконечно тонкий цилиндр) условие антирезонанса запишем в виде
Юае =Юо + ЮМ , (6)
где ю0 - собственная частота ферромагнитного резонанса в эллипсоиде.
Внутреннее магнитное поле в ферромагнитном эллипсоиде, помещенном во внешнее магнитное поле И 0 , равно
И0вН = И0 -НМ0, (7)
где N - тензор размагничивания.
Как следует из результатов математического моделирования (рис. 3), максимум коэффици-
управляется полем намагничивания И0 и зависит от геометрии и соотношения размеров магнитной нанорешетки (диаметра 2г нанопроволоки периода решетки а = 3,5г). При расстояниях между ферромагнитными нанопроволоками диаметром 2г = 25 нм больших, чем а > 100 нм [2, 3], моделью слоя невзаимодействующих магнитных нанопроволок является бесконечно тонкий цилиндр с продольным подмагничиванием внешним магнитным полем И0 = И0 у0,
ента прохождения
„12 К11
размагничивающие факторы N = N = 2п, Nz = 0. В этом случае собственная частота ферромагнитного резонанса однородного типа колебаний равна
» = Н0 + 2яМ0. (8)
У
При дальнейшем сокращении диаметра и расстояний между ферромагнитными нанопрово-локами а < 60 нм (переходе в масштаб длины обменного взаимодействия) обменное взаимодействие в системе сильно взаимодействующих нанопроволок, связанных силами магнитного ди-поль-дипольного взаимодействия, играет доминирующую роль. При малом диаметре 2г = 10 нм и расстояниях а = 17,5 нм двумерная магнитная нанорешетка приближается по свойствам к квазисплошной ферромагнитной тонкой пленке с касательным направлением намагниченности, при этом размагничивающие факторы Nx = 0, N = 4, Nz = 0.
В случае магнитной нанорешетки с высокой плотностью собственная частота ферромагнитного резонанса однородного типа прецессии намагниченностиравна
( \2
= Н0 (И0 + 4кИ0). (9)
Из (8), (9) следует, что для того, чтобы реализовать антирезонансный минимум поглощения в точке антирезонанса (6), требуется приложить большее значение внешнего магнитного поля Н 0 в случае магнитной нанорешетки с высокой плотностью (кривая 1 рис. 3), чем в случае магнитной нанорешетки с малой плотностью (кривая 4 рис. 3).
3. Результаты математического моделирования размерных эффектов, возникающих при дифракции в терагерцовом диапазоне
Анализ полученных результатов математического моделирования дифракции электромагнитного излучения на исследуемой магнитной нанокомпозитной структуре в виде решеток из нанопро-волоки показывает, что при фиксированной длине проволочных элементов снижение относительного коэффициента затухания ТЕМ-волны происходит как с увеличением диаметра нанопроволоки, так и с увеличением напряженности поля намагничивания (рис. 4). Наибольшие значения коэффициента отражения возникают при длине нанопроволоки менее 200 нм (рис. 5). Эффект антирезонанса проявляется при использовании ячеек нанорешеток размером более 500 нм.
0 200 400 600 800 1000
1 2
4 3
1Д'?1 Г"
-10
-20
-30
Рис. 4. Относительный коэффициент затухания ТЕМ-волны, проходящей через двумерную магнитную нанорешетку в зависимости от длины / нанопроволок при изменении величины поля намагничивания Н0 :
|^12| - коэффициент прохождения при Н0 Ф 0; ^^ - при Н0 = 0; Н0 = Н0у0 ; а = 3,5г, Ь = 1,25/; с = 2г; / = 30 ТГц; кривые: 1 - 2г = 10 нм, Н0 = 450 Э; 2 - 2г = 15 нм, Н0 = 105 Э ; 3 - 2г = 20 нм, Н0 = 75 Э ;
4 - 2г = 25 нм, Н0 = 50 Э
R
0,8 0,6 0,4 0,2
\a
ЛЧ J
0 2 DO 4 00 600 8 00 1000 /, нм
Рис. 5. Модуль коэффициента отражения Я^1 ТЕМ-волны от двумерной магнитной нанорешетки в зависимости от длины нанопроволок l при изменении величины поля намагничивания H0 :
H0 = H0y0; f = 30 ТГц; a = 3,5r, b = 1,25l; с = 2r ; кривые: 1 - 2r = 10 нм, H0 = 450 Э; 2 - 2r = 15 нм, H0 = 105 Э; 3 - 2r = 20 нм, H0 = 75 Э ; 4 - 2r = 25 нм, H0 = 50 Э
Разработанный подход к математическому моделированию размерных эффектов на нано-решетках из магнитных нанопроволок, базирующийся на методах вычислительной электродинамики, позволяетпреодолеть ограничения известных методов анализа в магнитостатическом приближении.
Список литературы
1. Исследование эффекта антирезонанса при воздействии лазерного излучения на магнитную нанокомпо-зитную структуру / Г. С. Макеева, О. А Голованов, В. Я. Савицкий, О. В. Филиппов // Труды Междунар. симп. Надежность и качество. - 2015. - Т. 2. - С. 142 - 145.
2. M. Pardavi-Horvath, P. E. Si, M. Vazquez, W. O. Rosa, and G. Badini // J. Appl. Phys. - 2008. - 103, 07D517.
3. Голованов, О. А. Оптические фильтры на основе опаловых матриц из диэлектрических наносфер / О. А. Голованов, В. Я. Савицкий, О. В. Филиппов // Надежность и качество сложных систем. - 2015. -№ 1. - С. 17-23.
Макеева Галина Степановна
доктор физико-математических наук, профессор, кафедра радиотехники и радиотехнических систем, Пензенский государственный университет (440026, Россия, г. Пенза, ул. Красная, 40) E-mail: radiotech@pnzgu.ru
Голованов Олег Александрович
доктор физико-математических наук, профессор,
кафедра общепрофессиональных дисциплин,
Военный учебно--научный центр Сухопутных войск
«ОВА ВС РФ» (филиал г. Пенза)
(440005, Россия, г. Пенза-5)
E-mail: golovanovol@mail.ru
Савицкий Владимир Яковлевич
доктор технических наук, профессор, кафедра общепрофессиональных дисциплин, Военный учебно-научный центр Сухопутных войск «ОВА ВС РФ» (филиал г. Пенза) (440005, Россия, г. Пенза-5) E-mail: W.savis@gmail.com
Makeeva Galina Stepanovna
doctor of physical and mathematical sciences,
professor, sub-department of radiotechnics
and radiotechnic systems,
Penza State University,
(440026, 40 Krasnaya street, Penza, Russia)
Golovanov Oleg Aleksandrovich
doctor of physical and mathematical sciences,
professor, sub-department of general trade science,
Militarian Educational Scientific Center
of the Land Forces «OVA VS RF» (branch Penza)
(440005, Penza-5 city, Russia)
Savitskiy Vladimir Yakovlevich doctor of technical sciences, professor, sub-department of general trade science, Militarian Educational Scientific Center of the Land Forces «OVA VS RF» (branch Penza) (440005, Penza-5 city, Russia)
Аннотация. Проведено математическое моделирование дифракции ТЕМ-волны на периодической решетке магнитных нанопроволок, базирующееся на решении трехмерной краевой задача дифракции для уравнений Максвелла с электродинамическими граничными условиями совместно с уравнением Ландау-Лифшица с учетом обменного взаимодействия. С помощью вычислительного алгоритма на основе метода автономных блоков с виртуальными каналами Флоке и магнитными нановключениями проведен электродинамический расчет ^-параметров многомодовой матрицы рассеяния двумерной решетки продольно намагниченных ферромагнитных нанопроволок. Сделан вывод, что максимум коэффициента прохождения ТЕМ-волны через решетку ферромагнитных нанопроволок управляется полем намагничивания вблизи точки антирезонанса в те-рагерцовом диапазоне частот и зависит от геометрии и соотношения размеров магнитной наноре-шетки.
Ключевые слова: дифракция ТЕМ-волны, двумерная решетка продольно намагниченных ферромагнитных нанопроволок, электродинамический расчет ^-параметров, управление полем намагничивания, антирезонанс, терагерцовый диапазон частот.
Abstract. Mathematical modelling of diffraction of TEM-waves on a periodic lattice magnetic nanowire, based on the decision three-dimensional regional a problem of diffraction for equations Maxwell with electrodynamic boundary conditions together with equation Landau-Lifshitsa taking into account exchange interaction is spent. By means of computing algorithm on the basis of a method of independent blocks with virtual channels the Flock and magnetic nanoim-pregnations electrodynamic calculation of S-parametres mods matrixes of dispersion of a two-dimensional lattice is carried out is longitudinal magnetised ferromagnetic nanowire. The conclusion is drawn, that the maximum of factor of passage of TEM-waves through a lattice ferromagnetic nanowire copes a field of magnetisation near to an antiresonance point in terahertz a range of frequencies and depends on geometry and a parity of the sizes magnetic nanolattice.
Key words: diffraction of TEM-waves, a two-dimensional lattice it is longitudinal magnetised ferromagnetic nanowire, electrodynamic calculation of S-parametres, management of a magnetisation field, an antiresonance, terahertz a range of frequencies.
УДК 538. 6
Макеева, Г. С.
Оценка антирезонанса при воздействии лазерного излучения на нанокомпозитную структуру /
Г. С. Макеева, О. А. Голованов, В. Я. Савицкий // Надежность и качество сложных систем. - 2015. -№ 3 (11). - С. 53-60.
