ФИЗИКА PHYSICS
УДК 004.031.4; 025.4.036
DOI: 10.18413/2075-4639-2018-50-1-47-54
ПОВЕРХНОСТНЫЕ СОСТОЯНИЯ В НЕЛИНЕЙНОМ ПОЛУПРОСТРАНСТВЕ
SURFACE STATES IN NONLINEAR HALF-SPACE
С.Е. Савотченко S.E. Savotchenko
Белгородский государственный технологический университет имени В.Г. Шухова, ул. Костюкова, 46, 308012, Белгород, Россия
Belgorod State Technological University named after V.G. Shukhov, Kostyukova Str., 46, 308012, Belgorod, Russia
E-mail: [email protected]
Аннотация
Рассмотрены вопросы существования поверхностных колебательных состояний в полуограниченных ангармонических кристаллах с различными знаками нелинейности. Предложена модель, математическая формулировка которой представляет собой одномерную краевую задачу для нелинейного уравнения Шредингера на полуоси. В рассматриваемой системе в зависимости от значения частоты получены несколько типов стационарных колебательных состояний, описывающих локальные возбуждения вблизи поверхности. Получены дисперсионные соотношения, определяющие значения частот таких состояний. В предельных случаях получены выражения для частоты стационарных колебательных состояний в явном виде.
Abstract
The problems of the existence of surface vibrational states in semibounded anharmonic crystals with different signs of nonlinearity are considered. A model is proposed whose mathematical formulation is a one-dimensional boundary-value problem for the nonlinear Schrödinger equation on the half-axis. In the system under consideration, depending on the frequency value, several types of stationary vi-brational states are obtained that describe local excitations near the surface. Dispersion relations that determine the frequencies of such states are obtained. In limiting cases, expressions for the frequency of stationary vibrational states are obtained in explicit form.
Ключевые слова: нелинейное уравнение Шредингера, поверхность, солитон, локализованные состояния.
Keywords: nonlinear Schrödinger equation, surface, soliton, localized states.
Введение
Изучение различных поверхностных состояний в кристаллах, в том числе в ангармонических, проводится достаточно давно [1-8]. Нелинейные свойства, которые демонстрируют поверхностные волны, наблюдались во многих экспериментах [9, 10]. Интерес к исследованию нелинейных поверхностных состояний возрастает в связи с интенсивным развитием физики нелинейных явлений и их применением в технических приложениях. Для
технических приложении нелинейных поверхностных состоянии важна возможность распространения устойчивых слабозатухающих импульсов или волновых пакетов [11].
При теоретическом изучении нелинейных поверхностных волн звуковой природы в упругих средах, как правило, используются уравнения Буссинеска и различные его обобщения [4-6]. В тоже время, нелинейные волны электромагнитной природы в оптических средах с показателем преломления, зависящим квадратично от модуля амплитуды напряженности электрического поля, описываются нелинейным уравнением Шредингера (НУШ) [2, 11-13].
НУШ широко используется при формулировке математических моделей описания локализованных вблизи дефектов малоамплитудных нелинейных колебаний в ангармонических кристаллах [14]. Решения НУШ описывают локализованные колебательные состояния с частотами вне зоны сплошного спектра, зависящими от амплитуды колебаний, в ангармоническом кристалле с точечным дефектом [15].
Для описания возбуждений в средах с пространственной дисперсией и коллективных возбуждений молекулярных цепочках используются обобщения НУШ с производными четвертого порядка [16-22]. В [23-25] было предложено модифицировать потенциал взаимодействия возбуждения в нелинейной среде с дефектом, обладающего внутренней структурой, при учете дальнодействующих сил межатомного взаимодействия. В [26] на основе НУШ были проанализированные особенности взаимодействия локализованных состояний вблизи границы раздела нелинейных сред.
В данной работе предлагается описать различные виды стационарных колебательных состояний, которые возникают вблизи поверхности нелинейной среды. Для этого предлагается простая модель, использующая НУШ. В рамках данной модели предполагается, что возбуждение свободно распространяется вдоль поверхности кристалла, а неличные свойства среды существенны в направлении, перпендикулярном поверхности. Тогда задача сводится к одномерной. В результате появляется возможность аналитически описать локализацию возбуждений вблизи поверхности солитонными решениями НУШ, форма которых зависит от знака параметра нелинейности, а амплитуда и частота определяются характеристиками среды и поверхности как плоского дефекта.
1. Формулировка математической модели
Рассмотрим ангармонический полубесконечный кристалл в дальнейшем который будем называть нелинейной средой. Плоская поверхность кристалла лежит в плоскости уОх и проходит через начало координат перпендикулярно оси Ох. Кристалл занимает полупространство в области х > 0.
Будем предполагать, что возмущение параметров кристалла, создаваемое его поверхностью как плоским дефектом, сосредоточено на расстояниях, существенно меньших размеров изучаемых возбуждений, и поэтому считается локальным. Будем рассматривать процессы локализации возбуждений вблизи поверхности на основе НУШ (принято й=1):
г а^-УМ2 ^, (1)
дг 2т
где А = + + - трехмерный оператор Лапласа, т - эффективная масса возбужде-
дх ду дх
ний, у - параметр нелинейности среды (может иметь различные знаки), характеризующий ангармонизм взаимодействия элементарных возбуждений в кристалле.
Если считать, что возбуждение в форме волны свободно может распространяться вдоль поверхности, то задача нахождения стационарных состояний с энергией Е сводится к одномерной с помощью подстановки в НУШ
х, у, г, г) = и(х) ехр(к у + к^ - гЕг).
Поскольку нас будут интересовать стационарные колебательные состояния, то удобно ввести обозначения: ra2 = E - квадрат частоты колебаний, s2 = 1I2m - скорость распространения возбуждения. В результате НУШ (1) становится одномерным и принимает вид:
s2u" (x) + (ra2 - ra¡2)u(x) + у|u(x)|2u(x) = О, (2)
где ra 2 = (к2 + к 2)s2 - значение уровня дна энергетической зоны в сечении изоэнергетиче-ской поверхности.
Поскольку НУШ (2) рассматривается на полуоси x > 0, то для определения поверхностных стационарных состояний следует задать граничные условия, которые целесообразно выбрать в виде:
u(+0) = uo; (3)
u " (+0) = Q, (4)
где параметры u0 и Q описывают возмущение характеристик кристалла вблизи поверхности. В частности, они характеризуют изменение массы атомов поверхностного слоя и силового взаимодействия между атомами поверхностного и приповерхностного слоев. В силу локальности такого возмущения характеристик кристалла вблизи поверхности данные величины можно считать постоянными.
Таким образом, математическая формулировка нахождения стационарных колебательных состояний в рамках предложенной модели сводится к решению краевой задачи на полуоси для НУШ (2) с граничными условиями (3), (4). Искомое решение данной краевой задачи должно быть ограничено: |u(x) < const V xe(0; œ).
В линейном бесконечном одномерном кристалле без дефектов распространяются свободные волны с квадратичным законом дисперсии. В ангармоническом кристалле без дефектов могут распространяться локализованные возбуждения - солитоны, и нелинейные волны, описываемые периодическими решениями НУШ, которые выражаются через эллиптические функции, называемые кноидальными волнами. В зависимости от знака параметра нелинейности существуют различные типы решений. В данной работе рассматриваются только локализованные состояния НУШ.
2. Поверхностные состояния в кристалле с отрицательным ангармонизмом
Сначала рассмотрим случай отрицательного ангармонизма, когда у < 0. Для удобства обозначим параметр нелинейности g = -у > 0. В рассматриваемой модели существует несколько типов локализованных колебательных состояний в различных частотных диапазонах.
2.1. Поверхностное состояние первого типа
В области частот ra > ra0 уравнение (2), имеет решение вида:
u(x) = Atthqt(x - x0) , (5)
параметры которого получаются после подстановки (5) в уравнение (2). Волновое число имеет вид:
q2 = (ra2 -ra2)I2s2, (6)
а амплитуда решения (5):
A2 = 2s VI g . (7)
Такое решение справедливо и ограничено всюду на всей числовой оси при -œ < x < œ и описывает свободно распространяющееся возбуждение в форме ступени. Для такого случая бесконечного кристалла частота ra и положение «центра» возбуждения Х0 являются свободными параметрами. Если At>0, то решение (5) называют кинком, а если At<0, то - анти-кинком. Именно антикинк удовлетворяет требованию |u| ^ 0 при x ^ +œ.
Для случая полуограниченного кристалла функция (5) является решением уравнения
(2) на полуоси х > 0, но ее параметры теперь должны определяться из граничных условий
(3) и (4). Волновое число и амплитуда определяются такими же выражениями (6) и (7).
Подстановка (5) в граничные условия (3) и (4) позволяет получить в явном виде частоту колебаний данного стационарного состояния:
ю2 = ®2 +®2, (8)
-"о + юь
где ю2 = 4Я + ), а также выражение, определяющее параметр хо:
5'
42
х0 =--— ЛпЬ —
ю
ь
Я ио
Л
V 5 юь у
(9)
Подстановка (8) в (6) и (7) позволяет выразить волновое число и амплитуду через параметры полуограниченного кристалла:
=ю2/252, (10)
А2 = ю\/ я . (11)
Из (8) и (9) следует, что поверхностное состояние рассматриваемого вида, которое описывается функцией (5), существует при Q > 0 .
Из (9) следует, что при ио < 0 положение центра возбуждения хо > 0, то есть будет находиться в области, занимаемой полуограниченным кристаллом. В противоположной ситуации при ио > 0, приводящей к хо < 0, решение (5) будет удовлетворять требованию исчезновения возбуждения при удалении от поверхности кристалла на бесконечности.
Для ненагруженной, то есть свободной от примесей поверхности, параметр ио = о, тогда из (9) следует, что хо = о. В этом случае из (8) получается частота локализованного вблизи свободной поверхности состояния:
ю2 =ю0 + Qs(Ьg )1/2. (12)
При этом волновое число и амплитуда выражаются через характеристики среды и дефекта и примут вид:
42 = (Я/2)1/2 Q/5 , (13)
А2 = 5Q(2/я)1/2 . (14)
Таким образом, решение вида (5), удовлетворяющее граничным условиям (3) и (4), является поверхностным состоянием с полностью определяемыми через характеристики кристалла параметрами.
2.2. Поверхностное состояние второго типа
В области частот ю < юо уравнение (2) имеет решение вида:
и(х) = А5/$&ч(х - х0), (15)
Параметры решения (15) имеют вид:
42 = (ю2-ю2)/52, (16)
А2 = 25 2д2/ я . (17)
Для бесконечного кристалла решение (15) будет всюду ограничено при выполнении требования: хо ^ 0. Как и для решения (5) на всей числовой оси, частота ю и положение «центра» возбуждения хо являются свободными параметрами решения (15). В силу такого требования обычно при построении физических моделей в безграничных кристаллах решение вида (15) не используется. Однако оно находит применение для описания эффектов локализации возбуждений вблизи дефектов [15].
В случае полубесконечного кристалла для ограниченности решения (15) должно выполняться требование: хо < 0. Тогда решение (15) ограничено на полуоси х > 0 и удовлетворяет условию на бесконечности: \и\ ^ 0 при х ^
Подстановка (15) в граничные условия (3) и (4) приводит к соотношениям:
(ш02-®2)2(u0V 2 + ш02— ю2)1/2 = —«0Qs( g / 2)2, (18)
^ -ю2)1/2 f 2 f2
shqx0 =
g
(19)
Выражение (18) представляет собой дисперсионное соотношение, определяющее частоту поверхностных колебательных состояний с функцией вида (15), а выражение (19) определяет параметр хо.
В предельном случае малых значений qxo << 1 из (19) можно получить параметр хо в явном виде: х0 = (2/g)12s /u0 . Поверхностное состояние данного типа реализуется при uo < 0.
Таким образом, решение вида (15), удовлетворяющее граничным условиям (3) и (4), является поверхностным состоянием с полностью определяемыми через характеристики кристалла параметрами в отличие от свободно распространяющих возбуждений в бесконечном кристалле.
3. Поверхностные состояния в кристалле с положительным ангармонизмом
Рассмотрим теперь случай положительного ангармонизма, когда у > 0. Для такого знака параметра нелинейности НУШ (2) имеет решение солитонного типа в диапазоне частот ю < юо:
«( х) = A / A q( х — х0). (20)
Волновое число q в (20) определяется выражением (16), а амплитуда имеет вид:
A2= 2s 2q2/у. (21)
Решение (20) без наложения дополнительных требований всюду ограничено и удовлетворяет условию на бесконечности: |u| ^ 0 при х ^ да.
Для случая полуограниченного кристалла функция (20) является решением уравнения (2) на полуоси х > 0, но ее параметры теперь должны определяться из граничных условий (3) и (4).
Подстановка (20) в граничные условия (3) и (4) приводит к соотношениям:
(ю2 — ю2)2^ —ю2 — «02у /2)"° = u0Qs(y / 2)2, (22)
chqx = (ю2 — ю2)1'2 f2Г2. (23)
«0 Ш
Выражение (22) представляет собой дисперсионное соотношение, определяющее частоту поверхностных колебательных состояний с функцией вида (20), а выражение (23) определяет параметр хо.
Из (22) следует, что максимальная частота такого колебательного состояния ниже значения юо, и она определяется выражением ю^ =ю2 — «2у /2. Другими словами, локализованные поверхностные состояния в среде с положительной нелинейностью существуют в диапазоне частот ю < ютах. Получается, что в полуограниченном кристалле допустимый диапазон существования стационарных состояний НУШ с положительной нелинейностью ограничивается более низкой частотой, чем в бесконечном кристалле.
В предельном случае малых значений qxo << 1 из (23) вытекает, что частота поверхностного колебательного состояния (20) приближается к наибольшему возможному значению ю = ютах.
u
0
Заключение
На основе модели полуограниченного ангармонического кристалла, возбуждения в котором описываются НУШ, рассмотрены вопросы существования различных видов стационарных локализованных вдоль поверхности состояний в широком интервале частот. Рассмотрены случаи как положительной, так и отрицательной нелинейности сред.
Математическая формулировка модели представляет собой краевую задачу для НУШ на полуоси. Получены решения сформулированной краевой задачи, описывающие локализацию возбуждений вблизи поверхности.
Параметры модели можно разделить на две группы. Первая - это параметры, характеризующие свойства среды, в которой распространяется возбуждение. К ним относятся скорость распространения возбуждения s, параметр нелинейности у, граничная частота Ш0. Вторая группа - это характеристики поверхности как плоского дефекта, описывающие изменение вблизи поверхности таких свойств кристаллической решетки как массы атомов поверхностного слоя U0 и константы силового взаимодействия между атомами поверхностного и приповерхностного слоев Q.
В случае бесконечного кристалла без дефектов в нелинейных средах, возбуждения в которых описываются НУШ, могут распространяться свободные солитоны, форма которых определяется знаком нелинейности и частотой. Такие солитоны являются двухпараметри-ческими в том смысле, что соответствующее решение НУШ содержит два свободных параметра, как правило, в качестве которых выступают частота ш и положение максимума Х0. Остальные параметры солитоны выражаются через них, а также через характеристики среды, то есть параметры первой группы. Если же рассматривать полуограниченный кристалл, то за счет наличия поверхности как плоского дефекта параметры поверхностных состояний будут определяться полностью все через параметры первой и второй группы. Другими словами, поверхностные состояния описываются солитонными решениями НУШ, не имеющими свободных параметров.
Показано, что в зависимости от частоты в среде с отрицательным ангармонизмом могут существовать два вида стационарных колебательных поверхностных состояний, а в среде с положительным ангармонизмом - одного вида. При частотах, ниже граничного значения Ш0 могут существовать поверхностные колебательные состояния как в средах с отрицательным ангармонизмом, так и в средах с положительным ангармонизмом. При частотах, выше граничного значения Ш0 поверхностные колебательные состояния могут существовать только в средах с отрицательным ангармонизмом.
В среде с положительным ангармонизмом меняется диапазон существования стационарных состояний НУШ. В полубесконечном кристалле он ограничивается сверху более низкой частотой, чем в бесконечном кристалле.
Полученные в данной работе теоретические результаты дополняют исследования поверхностных упругих нелинейных волн, распространяющихся вблизи свободной поверхности [3], а также поверхностных магнитных поляритонов в ферромагнитной полуограниченной среде [27], поскольку в них приведены решения НУШ на полуоси одного типа только для самофокусирующией среды с кубической нелинейностью для нулевых граничных условий на поверхности.
Список литературы References
1. Павло В.И., Солодов И.Ю. 1977. Нелинейные свойства поверхностных волн в твердых телах. ФТТ, 19, 10: 2948-2954.
Pavlo V.I., Solodov I.Ju. 1977. Nelinejnye svojstva poverhnostnyh voln v tverdyh telah. FTT, 19, 10: 2948-2954
2. Михалаке Д., Назмитдинов Р.Г., Федянин В.К. 1989. Нелинейные оптические волны в слоистых структурах. Физика элементарных частиц и атомного ядра, 20, 1: 198-253.
Mihalake D., Nazmitdinov R.G., Fedjanin V.K. 1989. Nelinejnye opticheskie volny v sloistyh strukturah. Fizika jelementarnyh chastic i atomnogo jadra, 20, 1: 198-253.
3. Горенцвейг В.И., Кившарь Ю.С., Косевич А.М., Сыркин Е.С. 1990. Самоиндуцированные нелинейные сдвиговые поверхностные акусктические волны в кристаллах. ФНТ, 16, 11: 1472-1482.
Gorencvejg V.I., Kivshar' Ju.S., Kosevich A.M., Syrkin E.S. 1990. Samoinducirovannye nelinejnye sdvigovye poverhnostnye akuskticheskie volny v kristallah. FNT, 16, 11: 1472-1482
4. Mayer A.P. 1995. Surface acoustic waves in nonlinear elastic media. Physics Reports, 256, 45: 237-366.
5. Ковалев А.С., Сыркин Е.С., Можен Ж.А. 2002. Многомерные и поверхностные солитоны в нелинейной упругой среде. ФНТ, 28, 6: 635-647.
Kovalev A.S., Syrkin E.S., Mozhen Zh.A. 2002. Mnogomernye i poverhnostnye solitony v nelinejnoj uprugoj srede. FNT, 28, 6: 635-647
6. Ковалев А.С., Майер А.П., Соколова Е.С., Экль К. 2002. Солитоны в упругих пластинах. ФНТ, 28, 10: 1092-1102.
Kovalev A.S., Majer A.P., Sokolova E.S., Jekl' K. 2002. Solitony v uprugih plastinah. FNT, 28, 10: 1092-1102
7. Федоров Л.В., Ляхомская К.Д. 1997. Нелинейные поверхностные волны при учете эффекта насыщения. Письма в ЖЭТФ, 23, 23: 36-39.
Fedorov L.V., Ljahomskaja K.D. 1997. Nelinejnye poverhnostnye volny pri uchete jeffekta nasysh-henija. Pis'ma v ZhJeTF, 23, 23: 36-39
8. Коровай О.В., Хаджи П.И. 2003.Нелинейные поверхностные волны в симметричной трехслойной структуре, обусловленные генерацией экситонов и биэкситонов в полупроводниках. ФТТ,45, 2: 364-368.
Korovaj O.V., Hadzhi P.I. 2003.Nelinejnye poverhnostnye volny v simmetrichnoj trehslojnoj strukture, obuslovlennye generaciej jeksitonov i bijeksitonov v poluprovodnikah. FTT,45, 2: 364-368
9. Наянов В.И. 1986. Поверхностные акустические кноидальные волны и солитоны в структуре LiNbO3-пленка SiO. Письма в ЖЭТФ, 44: 245-249.
Najanov V.I. 1986. Poverhnostnye akusticheskie knoidal'nye volny i solitony v strukture LiNbO3-plenka SiO. Pis'ma v ZhJeTF, 44: 245-249
10. Усиевич Б. А., Нурлигареев Д.Х., Сычугов В.А., Ивлева Л.И., Лыков П.А., Богодаев Н.В. 2010. Нелинейные поверхностные волны на границе фоторефрактивного кристалла. Квантовая электроника, 40, 5: 437-440.
Usievich B. A., Nurligareev D.H., Sychugov V.A., Ivleva L.I., Lykov P.A., Bogodaev N.V. 2010. Nelinejnye poverhnostnye volny na granice fotorefraktivnogo kristalla. Kvantovaja jelektronika, 40, 5: 437-440
11. Kivshar U.S., Kosevich A.M., Chubykalo O.A. 1990. Radiative effects in the theory of beam propagation at nonlinear interfaces. Phys. Rev. A. :41: 3: 1677-1688.
12.Ахмедиев Н.Н., Корнеев В.И., Кузьменко Ю.В. 1985. Возбуждение нелинейных поверхностных волн гауссовыми световыми пучками. ЖЭТФ, 88, 1: 107-115.
Ahmediev N.N., Korneev V.I., Kuz'menko Ju.V. 1985. Vozbuzhdenie nelinejnyh poverhnostnyh voln gaussovymi svetovymi puchkami. ZhJeTF, 88, 1: 107-115
13. Abdullaev F.Kh., Baizakov B.B., Umarov B.A. 1998. Resonance phenomena in interaction of a spatial soliton with the modulated interface of two nonlinear media. Optics Communications, 156: 341-346.
14. Герасимчук И.В., Ковалев А.С. 2000. Локализация нелинейных волн в слоистых средах. ФНТ, 26, 8: 799-809.
Gerasimchuk I.V., Kovalev A.S. 2000. Lokalizacija nelinejnyh voln v sloistyh sredah FNT, 26, 8: 799-809
15. Богдан М.М., Герасимчук И.В., Ковалев А.С. 1997. Динамика и устойчивость локализованных мод в нелинейных средах с точечными дефектами. ФНТ, 23, 2: 197-207.
Bogdan M.M., Gerasimchuk I.V., Kovalev A.S. 1997. Dinamika i ustojchivost' lokalizovannyh mod v nelinejnyh sredah s tochechnymi defektami. FNT, 23, 2: 197-207.
16. Косевич А.М., Савотченко С.Е. 1999. Особенности динамики одномерных дискретных систем с взаимодействием не только ближайших соседей и роль высшей дисперсии в солитонной динамике. ФНТ, 25, 7: 737-747.
Kosevich A.M., Savotchenko S.E. 1999. Osobennosti dinamiki odnomernyh diskretnyh sistem s vzaimodejstviem ne tol'ko blizhajshih sosedej i rol' vysshej dispersii v solitonnoj dinamike. FNT, 25, 7: 737-747.
17. Kosevich A.M., Savotchenko S.E. 2000.Forced vibrations and resonance wave scattering on impurity in 1D discrete lattice with nearest- and next-nearest neighbors interaction. Physica B, 284-288: 1551-1552.
18. Cавотченко C.E. 2000. Рассеяние волн дефектами в средах с пространственной дисперсией и безызлучательные динамические солитоны. Известия высших учебных заведений. Физика, 43, 10: 876-881.
Savotchenko S.E. 2000. Rassejanie voln defektami v sredah s prostranstvennoj dispersiej i bezyzlu-chatel'nye dinamicheskie solitony. Izvestija vysshih uchebnyh zavedenij. Fizika, 43, 10: 876-881
19. Cавотченко C.E. 2004. Локализация волн вблизи интерфейса нелинейных сред с пространственной дисперсией. Известия высших учебных заведений. Физика, 47, 5: 79-84.
Savotchenko S.E. 2004. Lokalizacija voln vblizi interfejsa nelinejnyh sred s prostranstvennoj dispersiej. Izvestija vysshih uchebnyh zavedenij. Fizika, 47, 5: 79-84.
20. Krasilnikov V.V., Savotchenko S.E. 2004. Peculiarities of solito n motion in molecular systems with high dispersion. Phys. Stat. Sol. (c), 1, 11: 2757-2760.
21. Cавотченко C.E. 2005. Нелинейные коллективные возбуждения в квазиодномерных структурах при наличии пространственной дисперсии. Известия вузов. Физика, 48, 9: 24-27.
Savotchenko S.E. 2005. Nelinejnye kollektivnye vozbuzhdenija v kvaziodnomernyh strukturah pri nalichii prostranstvennoj dispersii. Izvestija vuzov. Fizika, 48, 9: 24-27.
22. Cавотченко C.E. 2006. Особенности нелинейной динамики квазичастиц в молекулярных структурах с водородными связями при наличии взаимодействия не только ближайших соседей. Известия вузов. Физика, 49, 2: 52-56.
Savotchenko S.E. 2006. Osobennosti nelinejnoj dinamiki kvazichastic v molekuljarnyh strukturah s vodorodnymi svjazjami pri nalichii vzaimodejstvija ne tol'ko blizhajshih sosedej. Izvestija vuzov. Fizika, 49, 2: 52-56.
23. Красильников В.В., Cавотченко C.E. 2015. Особенности рассеяния частиц точечным дефектом с внутренней структурой. Известия тульского государственного университета. Eстествен-ные науки, 4: 178-183.
Krasil'nikov V.V., Savotchenko S.E. 2015. Osobennosti rassejanija chastic tochechnym defektom s vnutrennej strukturoj. Izvestija tul'skogo gosudarstvennogo universiteta. Estestvennye nauki, 4: 178-183.
24. Cавотченко C.E. 2016. Особенности взаимодействия волны с точечным дефектом в среде с пространственной дисперсией. Известия Воронежского государственного педагогического университета, 1(270): 196-199.
Savotchenko S.E. 2016. Osobennosti vzaimodejstvija volny s tochechnym defektom v srede s prostranstvennoj dispersiej. Izvestija Voronezhskogo gosudarstvennogo pedagogicheskogo universiteta, 1(270): 196-199.
25. Cавотченко C.E. 2016. Особенности локализации нелинейных возбуждений вблизи дефекта с внутренней структурой. Вестник Воронежского государственного университета. Cерия: Физика. Математика, 4: 51-59.
Savotchenko S.E2016.. Osobennosti lokalizacii nelinejnyh vozbuzhdenij vblizi defekta s vnutrennej strukturoj. Vestnik Voronezhskogo gosudarstvennogo universiteta. Serija: Fizika. Matematika, 4: 51-59.
26. Cавотченко C.E. 2017. Взаимодействие локализованных состояний вблизи границы раздела нелинейных сред. Конденсированные среды и межфазные границы, 2: 291-295.
Savotchenko, S.E. 2017. Vzaimodejstvie lokalizovannyh sostojanij vblizi granicy razdela nelinejnyh sred. Kondensirovannye sredy i mezhfaznye granicy, 2: 291-295.
27. Дикштейн И^. Никитов Д.С, Никитов C.A. 1998. Нелинейные самолокализованные поверхностные магнитные поляритоны в ферромагнитной среде. ФТТ, 40, 10: 1885-1889.
Dikshtejn I.E. Nikitov D.S., Nikitov S.A. 1998. Nelinejnye samolokalizovannye poverhnostnye magnitnye poljaritony v ferromagnitnoj srede. FTT, 40, 10: 1885-1889.