Аннотация:
Рассматривается одномерный оператор Шрёдингера на единичном отрезке с условием Дирихле, возмущенный оператором сдвига. Основной результат посвящен асимптотике собственных значений этого оператора по номеру, которая является равномерной по сдвигу. В ней явно выделены члены, порождаемые оператором сдвига. Установлено, что система собственных и присоединенных функций для рассматриваемого оператора образует базис Бари в пространстве функций, интегрируемых с квадратом на единичном отрезке.
Библиография: 29 наименований.
Дифференциально-разностные уравнения являются важным и интересным объектом для изучения. Сдвиги в данных уравнениях часто интерпретируют как запаздывание (или опережение), суть которого состоит в том, что протекание физического процесса определяется его характеристиками не только в данный момент времени, но и в предыдущие моменты. В этом случае независимой переменной придается смысл времени. Задачи, в которых сдвиг присутствует в пространственных переменных, возникают в теории многослойных пластин и оболочек, теории многомерных диффузионных процессов, а также теории нелинейных оптических систем с двумерной обратной связью. Исследованию таких уравнений посвящена весьма обширная литература (см. [1]–[4] и приведенные там ссылки). Также имеется большое число статей, касающихся изучения краевых задач для эллиптических функционально-дифференциальных уравнений. Здесь прежде всего выделим книгу [5] и недавние работы [6]–[11]. В цитированных источниках исследовались вопросы устойчивости, существования, единственности и гладкости решений. Стоит также отметить статьи [12]–[14], посвященные эволюционным дифференциально-разностным уравнениям. Помимо них, изучались уравнения с нелокальностями вида растяжений или сжатий (см. [15]–[18]).
Несмотря на то, что исследование дифференциально-разностных уравнений ведется достаточно активно, вопросы спектральных свойств соответствующих операторов мало затрагивались. Отметим обзорную статью [6], в которой рассматривались вопросы характера спектра, а также полноты и минимальности собственных и присоединенных (корневых) функций сильно эллиптических дифференциально-разностных операторов. Кроме того, в работе [19] был получен критерий сохранения гладкости обобщенных собственных функций для такого класса операторов.
Отдельно подчеркнем, что дифференциально-разностные операторы, в которых сдвиг рассматривался бы как параметр, ранее детально не исследовались и подобное возмущение не является классическим примером в современной теории возмущений. Насколько нам известно, была изучена лишь одна подобная задача [16], [18; гл. 5, § 5.3, § 5.4] о поведении собственных значений дифференциального оператора с растяжениями и сжатиями при малых вариациях области. В приведенных работах было показано, что собственные значения такого оператора сходятся к собственным значениям оператора на невозмущенной области и были получены оценки скорости сходимости.
Вопросы резольвентной сходимости общих операторов с малыми переменными сдвигами рассматривались в нашей недавней работе [20]. Здесь изучались эллиптические операторы четного порядка $2m$ в произвольных областях с общими краевыми условиями, в которых также могут присутствовать малые сдвиги. При этом предполагалось, что сдвиги допускаются только в членах младшего порядка, а главная часть остается классическим дифференциальным оператором. Было показано, что в пределе малых сдвигов возникает дифференциальный оператор, соответствующий нулевым сдвигам. Основной результат указанной работы посвящен доказательству равномерной резольвентной сходимости в норме операторов, действующих из $L_2$ в $W_2^m$. Кроме того, были установлены явные оценки скорости сходимости. На основе этих результатов была показана сходимость спектров и псевдоспектров возмущенных операторов к аналогичным множествам предельных операторов.
В настоящей статье мы продолжаем исследование операторов со сдвигами в младших членах. А именно, мы рассматриваем оператор Шрёдингера на отрезке с краевыми условиями Дирихле, возмущенным оператором исследуемого вида. Величина сдвига может быть произвольной, включая случай сколь угодно малого сдвига. Наш основной результат посвящен асимптотике собственных значений для больших номеров, которая является равномерной по сдвигу. Помимо этого мы также устанавливаем базисность в $L_2(0, 1)$ системы корневых функций для рассматриваемого оператора. Интерес к изучению указанных вопросов связывается с двумя основными причинами: нестандартное возмущение сдвигом и нетривиальные высокочастотные эффекты, которые оно порождает. Хотя задача о спектральных асимптотиках для задачи Штурма–Лиувилля является классической и хорошо исследованной (см., например, монографии [21], [22]), наличие возмущения, зависящего от параметра, естественным образом влечет зависимость остаточного члена в спектральных асимптотиках не только от номера, но и от такого параметра. При одновременном росте номера собственного значения и изменении параметра оценка остатка может разрушаться. В качестве простейшего примера можно привести функцию вида $\varepsilon n/(\varepsilon^6 n^3+1)$. Для каждого фиксированного $\varepsilon$ при больших $n$ она стремится к нулю как $O(n^{-2})$. Однако если при этом уменьшать $\varepsilon$, скажем, выбирая $\varepsilon=n^{-1}$, то эта функция перестает быть малой, а при $\varepsilon=n^{-1/2}$ даже начинает расти!
Применение классической теории возмущений также не дает подходящего решения данной проблемы, так как даже наличие равномерной резольвентной сходимости означает сходимость спектров лишь в каждом компакте на комплексной плоскости. При этом скорость сходимости, вообще говоря, зависит от выбора компакта. Таким образом, наиболее разумный подход – это вывести асимптотику по большим $n$ с оценкой остатка, равномерной по параметру, описывающему возмущение. Для получения такой асимптотики требуется модификация имеющихся подходов и техник. Решение данной проблемы будет осуществляться с помощью метода подобных операторов, предложенного А. Г. Баскаковым. Основная идея этого метода состоит в построении преобразования подобия, позволяющего подходящим образом представить оператор в виде блочно-диагональной матрицы, блоки которой удается проанализировать. Это делает метод подобных операторов весьма эффективным средством спектрального анализа различных классов дифференциальных и разностных операторов (см. [23]–[27] и используемую там литературу).
§ 2. Постановка задачи и формулировка основных результатов
В пространстве $L_2(0, 1)$ рассмотрим самосопряженный оператор $\mathcal{A}$ вида $\mathcal{A}y=-y''$ на области определения $\mathring{W}_2^2(0, 1)$. Здесь и далее через $\mathring{W}_2^j(0, 1)$,$j\geqslant 1$, обозначается подпространство функций из пространства Соболева $W_2^j(0, 1)$, обращающихся в нуль на концах отрезка.
Основным объектом исследования настоящей работы является возмущение оператора $\mathcal{A}$, использующее оператор сдвига. Для описания такого оператора нам понадобятся оператор продолжения нулем вне интервала $(0, 1)$, действующий из $L_2(0, 1)$ в $L_2(\mathbb{R})$ и оператор сужения на $(0, 1)$, действующий из $L_2(\mathbb{R})$ в $L_2(0, 1)$. Они задаются формулами
Оператор сдвига $\mathcal{T}^\alpha$ в $L_2(\mathbb{R})$ определим равенством $(\mathcal{T}^\alpha y)(x)=y(x+\alpha)$,$x\in \mathbb{R}$, где $\alpha\in[0,1]$ – величина сдвига.
Пусть $V$ и $S$ – комплекснозначные функции из пространства $C^1[0, 1]$. Возмущающий оператор $\mathcal{B}^\alpha$ в пространстве $L_2(0, 1)$ введем следующим образом:
где второе слагаемое в правой части понимается как обычное умножение функции $S$ на выражение в скобках. На произвольную функцию $y\in L_2(0, 1)$ этот оператор действует по правилу
где функция $y$ считается продолженной нулем вне отрезка $[0, 1]$, а результат действия сужается на данный отрезок. При $\alpha=0$ оператор $\mathcal{B}^\alpha$ превращается в обычный оператор умножения на потенциал $-V$, а при $\alpha=1$ – в оператор умножения на потенциал $-V+S$.
В настоящей работе изучается поведение спектра оператора $\mathcal{H}^\alpha=\mathcal{A}-\mathcal{B}^\alpha$ в пространстве $L_2(0, 1)$ на области определения $\mathfrak{D}(\mathcal{H}^\alpha):=\mathring{W}_2^2(0,1)$. Наш первый предварительный результат описывает основные свойства оператора $\mathcal{H}^\alpha$.
Теорема 1. Оператор $\mathcal{H}^\alpha$ является $m$-секториальным, и соответствующая ему секториальная замкнутая форма в пространстве $L_2(0, 1)$ задается равенством
на области определения $\mathfrak{D}(\mathfrak{h}^\alpha)=\mathring{W}_2^1(0, 1)$. Существует число $\Lambda$, не зависящее от $\alpha$, такое, что полуплоскость $\operatorname{Re}\lambda\leqslant \Lambda$ попадает в резольвентное множество оператора $\mathcal{H}^\alpha$ для каждого $\alpha\in[0,1]$. При этом оператор $\mathcal{H}^\alpha$ имеет компактную резольвенту, и его спектр состоит из счетного числа собственных значений с единственной точкой накопления в бесконечности.
В силу теоремы 1 собственные значения оператора $\mathcal{H}^\alpha$ расположены в полуплоскости $\operatorname{Re}\lambda\geqslant \Lambda$. Пронумеруем их в порядке возрастания их модулей и обозначим через $\lambda_n$,$n\in\mathbb{N}$. Через $\chi_I=\chi_I(x)$ определим характеристическую функцию отрезка $I$ на вещественной прямой.
Наш первый основной результат посвящен спектральным асимптотикам оператора $\mathcal{H}^\alpha$ при больших номерах $n$.
Теорема 2. Для собственных значений $\lambda_n$ оператора $\mathcal{H}^\alpha$ справедлива асимптотика
при $n\to +\infty$. Оценка остаточного члена равномерна по $\alpha\in[0,1]$.
Второй основной результат статьи описывает базисность корневых функций оператора $\mathcal{H}^\alpha$. Напомним, что базисы, которые порождаются системами проекторов, квадратично близкими к полным и минимальным системам ортогональных проекторов, называются базисами Бари.
Теорема 3. Система корневых функций оператора $\mathcal{H}^\alpha$ образует базис Бари в пространстве $L_2(0, 1)$.
Обсудим кратко полученные результаты. В теореме 2 установлена равномерная спектральная асимптотика (2.3), которую можно переформулировать следующим образом: существует фиксированное $n_0>0$ такое, что при $n>n_0$ выполнено равенство (2.3), где остаток оценивается величиной $Cn^{-2}$ с постоянной $C$, не зависящей от $\alpha$ и $n$. Формула (2.3) описывает поведение всего ансамбля собственных значений, за исключением их конечного числа, с номерами $n\geqslant n_0$ равномерно по $\alpha$.
При $\alpha=0$ в (2.3) пропадают все слагаемые, содержащие функцию $S$, и мы получаем классическую асимптотику собственных значений оператора Шрёдингера на отрезке (см. [21; гл. 1, § 5]):
При $\alpha=1$ исчезают слагаемые, содержащие $\sin (\pi n\alpha)$ и $\sin (2\pi n\alpha)$, и член с характеристической функцией, и асимптотика вновь оказывается классической:
Если $0<\alpha<1$, то в асимптотике собственных значений присутствуют три неклассических дополнительных слагаемых, которые порождаются исключительно оператором сдвига. Очевидные оценки
позволяют достаточно точно описать характер зависимости четвертого и пятого слагаемых в (2.3) от большого номера $n$ и параметра $\alpha$.
Третий член в асимптотике (2.3) устроен качественно иным образом. Переписав в нем произведение синусов как разность косинусов и проинтегрировав по частям, легко получаем, что
где оценка остатка равномерна по $\alpha$. Аналогично можно преобразовать и второе слагаемое в асимптотике (2.3). В результате эта асимптотика принимает вид
где оценка остатка равномерна по $\alpha$. Таким образом, при достаточно больших $n$ собственные значения оператора $\mathcal{H}^\alpha$ локализуются вдоль следующей кривой в комплексной плоскости:
Локализация происходит вблизи точек на этой кривой, которые соответствуют значениям параметра $t=\pi n$. Кривая напоминает график периодической функции, причем при $\alpha=0$ никаких осцилляций нет, а при возрастании параметра $\alpha$ от нуля период осцилляций кривой уменьшается и осцилляции становятся более заметными. При $\alpha\to 1-0$ коэффициент $a_\alpha$ стремится к нулю и осциллирующий косинус в определении кривой $K_\alpha$ пропадает. Асимптотика (2.7) отличается от классической асимптотики третьим слагаемым с $\cos(\pi n \alpha)$, и его наличие описывает высокочастотный эффект, порождаемый сдвигом. Одновременно с этим подчеркнем, что основной результат статьи [20] в применении к оператору $\mathcal{H}^\alpha$ утверждает сходимость резольвенты этого оператора к резольвенте предельного оператора $\mathcal{H}^0$ по норме операторов, действующих из $L_2(0, 1)$ в $W_2^1(0, 1)$, и при $\alpha\to+0$ каждое собственное значение оператора $\mathcal{H}^\alpha$ сходится к соответствующему собственному значению оператора $\mathcal{H}^0y=-y'' + Vy$ на интервале $(0,1)$ с условиями Дирихле. Скорость сходимости неравномерна по номеру $n$, и возмущение оператором сдвига порождает описанный высокочастотный эффект.
Кривая $K_\alpha$ и примерное положение собственных значений показано на рис. 1 для потенциалов $V(x)=4x+6\mathrm{i} x^2$,$S(x)=\frac{1}{2}(1+\mathrm{i} x)$. В этом случае
Рис. 1, (a) и (b) приведены для значений $\alpha=0.2$ и $\alpha=0.7$. Помимо кривой, здесь также изображены три серии точек. Две серии, отмеченные темно-серыми квадратами и кругами, показывают значения суммы первых двух членов классических асимптотик (2.4) и (2.5). Третья серия точек, отмеченная черными ромбами, соответствует сумме первых членов асимптотики (2.7). Таким образом, из рис. 1 видно, что при указанных значениях $\alpha$ имеется существенный разрыв между членами классической асимптотики (2.4) и членами асимптотики (2.7), что является проявлением описанного выше высокочастотного эффекта. Аналогичный разрыв с асимптотикой (2.5) заметно уменьшается при $\alpha\to 1-0$.
Еще одна наша теорема посвящена базисности Бари корневых функций для такого класса операторов. Известные до настоящего момента результаты затрагивали только вопросы полноты и минимальности системы корневых функций для дифференциально-разностных операторов. Более того, отсутствуют даже результаты по классической базисности Рисса. Используемая в статье техника позволяет получить необходимые оценки на спектральные проекторы, и из теоремы Бари–Маркуса (см. [28; гл. VI, теорема 5.2]) непосредственно следует базисность Бари системы корневых функций в $L_2(0, 1)$. Отметим, что согласно той же теореме Бари–Маркуса из базисности Бари вытекает базисность Рисса, и потому полученный в теореме 3 результат более общий, чем классическая базисность Рисса корневых функций.
В заключение опишем структуру статьи. В § 3 устанавливается $m$-секториальность оператора $\mathcal{H}^\alpha$, а также доказывается теорема 1. В § 4 мы кратко опишем идею используемого метода исследования, а также проведем необходимые построения операторов подобия. В § 5 приведены дополнительные технические вычисления, которые вместе с результатами § 4 используются для доказательства теорем 2 и 3 в § 6.
Всюду далее символом $C$ с индексами внизу обозначаем различные константы, не зависящие от $\alpha$,$k$,$j$,$n$,$m$.
§ 3. $m$-секториальность и компактность резольвенты
Доказательство теоремы 1. Из теоремы 3 работы [20; § 5] следует, что форма $\mathfrak{h}^\alpha$ секториальна и замкнута, а ее числовая область расположена в секторе вида $\{z\in\mathbb{C}\colon |{\operatorname{Im}z}|\leqslant C_0(\operatorname{Re}z-C_1)\}$, где $C_0$ и $C_1$ – некоторые фиксированные постоянные, не зависящие от $\alpha$, причем $C_0>0$. Поэтому в силу первой теоремы о представлении [29; гл. VI, § 2.1, теорема 2.1] данной форме соответствует некоторый $m$-секториальный оператор, который обозначим через $\widetilde{\mathcal{H}}^\alpha$. Покажем, что $\widetilde{\mathcal{H}}^\alpha=\mathcal{H}^\alpha$. Пусть $u\in\mathfrak{D}(\mathcal{H}^\alpha)$,$v\in\mathfrak{D}(\mathfrak{h}^\alpha)$. Тогда интегрированием по частям легко установить, что
$$
\begin{equation*}
(\mathcal{H}^\alpha u, v)_{L_2(0, 1)}=\mathfrak{h}^\alpha(u, v).
\end{equation*}
\notag
$$
Отсюда следует, что $\mathfrak{D}(\mathcal{H}^\alpha)\subseteq \mathfrak{D}(\widetilde{\mathcal{H}}^\alpha)$ и оператор $\widetilde{\mathcal{H}}^\alpha$ является продолжением оператора $\mathcal{H}^\alpha$. Поэтому для доказательства совпадения операторов $\mathcal{H}^\alpha$ и $\widetilde{\mathcal{H}}^\alpha$ достаточно проверить, что $\mathfrak{D}(\mathcal{H}^\alpha)\supseteq \mathfrak{D}(\widetilde{\mathcal{H}}^\alpha)$.
Согласно первой теореме о представлении [29; гл. VI, § 2.1, теорема 2.1] область определения оператора $\widetilde{\mathcal{H}}^\alpha$ состоит из функций $u\in\mathfrak{D}(\mathfrak{h}^\alpha)$, для каждой из которых существует соответствующая функция $h\in L_2(0, 1)$, удовлетворяющая равенству
и в силу стандартных теорем о повышении гладкости решений эллиптических краевых задач сразу заключаем, что $u\in\mathring{W}_2^2(0, 1)= \mathfrak{D}(\mathcal{H}^\alpha)$. Поэтому $\mathfrak{D}(\widetilde{\mathcal{H}}^\alpha)\subseteq \mathfrak{D}(\mathcal{H}^\alpha)$, откуда уже следуют равенства $\mathfrak{D}(\widetilde{\mathcal{H}}^\alpha) =\mathfrak{D}(\mathcal{H}^\alpha)$ и $\widetilde{\mathcal{H}}^\alpha=\mathcal{H}^\alpha$. Учитывая описанные выше свойства формы $\mathfrak{h}^\alpha$, видим, что полуплоскость $\operatorname{Re}\lambda\leqslant C_1-1$ лежит в резольвентном множестве оператора $\mathcal{H}^\alpha$ для каждого $\alpha\in [0, 1]$.
Резольвента оператора $\mathcal{H}^\alpha$ является ограниченным оператором в $L_2(0, 1)$, а по теореме Банаха об обратном операторе – ограниченным оператором из $L_2(0, 1)$ в $W_2^2(0, 1)$. Из компактности вложения пространства $W_2^2(0, 1)$ в $L_2(0, 1)$ следует компактность резольвенты оператора $\mathcal{H}^\alpha$. Согласно теореме об операторе с компактной резольвентой [29; гл. 3, § 6.8, теорема 6.29] спектр оператора $\mathcal{H}^\alpha$ состоит из счетного числа собственных значений с единственной точкой накопления в бесконечности. Теорема 1 доказана.
Хорошо известно, что подобные операторы имеют одинаковые спектры (см. более подробно в [24; лемма 2.10] и [26; лемма 1]) и основная идея метода подобных операторов заключается в построении такого преобразования подобия, которое позволит свести задачу об изучении собственных значений и соответствующих собственных подпространств оператора $\mathcal{H}^\alpha$ к исследованию этих же характеристик для некоторого оператора блочно-диагонального вида.
Необходимое преобразование строится в два этапа. На первом исходный оператор $\mathcal{H}^\alpha$ преобразуется в оператор $\mathcal{A}-\mathcal{B}_0-\widetilde{\mathcal{B}}^\alpha$, где $\mathcal{B}_0$ – некоторый оператор с блочно-диагональной матрицей, а $\widetilde{\mathcal{B}}^\alpha$ принадлежит идеалу операторов Гильберта–Шмидта в $L_2(0, 1)$, обозначаемому через $\mathfrak{S}_2(L_2(0, 1))$ с нормой $\|\,{\cdot}\,\|_2$ (см. [28; гл. 3, § 9]). Конкретный вид оператора $\mathcal{B}_0$ будет определен ниже. На втором этапе осуществляется преобразование $\mathcal{A}-\mathcal{B}_0-\widetilde{\mathcal{B}}^\alpha$ в оператор $\mathcal{A}-\mathcal{B}_0-\mathcal{Q}$, где $\mathcal{Q}$ не только принадлежит $\mathfrak{S}_2(L_2(0, 1))$, но и имеет блочно-диагональную матрицу. Общую схему, идеи и основные определения данного метода можно найти в [24].
Начнем с очевидных свойств невозмущенного оператора $\mathcal{A}$. Он самосопряжен, имеет компактную резольвенту, его собственные значения и соответствующие ортонормированные в $L_2(0, 1)$ собственные функции и проекторы Рисса задаются в виде
Все собственные значения простые, а система собственных функций образует ортонормированный базис в $L_2(0, 1)$.
Применение леммы 1 из [20] ко второму слагаемому в операторе $\mathcal{B}^\alpha$ позволяет утверждать, что $\mathcal{B}^\alpha$ является подчиненным оператору $\mathcal{A}$. Банахово пространство операторов в $L_2(0,1)$, подчиненных $\mathcal{A}$, обозначим через $\mathfrak{L}_{\mathcal{A}}(L_2(0, 1))$. Норма на этом пространстве вводится равенством
Из [20] следует, что норма $\|\mathcal{B}^\alpha\|_{\mathcal{A}}$ ограничена равномерно по $\alpha$.
Через $\sigma(\mathcal{A})$ и $\rho(\mathcal{A})$ обозначим спектр и резольвентное множество оператора $\mathcal{A}$. Оператор, осуществляющий первое преобразование подобия, будет задаваться с помощью операторов $J$,$\Gamma\colon\mathfrak{S}_2(L_2(0, 1))\to\mathfrak{S}_2(L_2(0, 1))$, определяемых формулами
Данные операторы корректно определены и ограничены (см. [24; лемма 3.4]). Продолжения $J$ и $\Gamma$ на пространство $\mathfrak{L}_{\mathcal{A}}(L_2(0, 1))$ (которые будут обозначаться теми же символами) задаются соотношениями
где $\mathcal{X}\in\mathfrak{L}_{\mathcal{A}}(L_2(0, 1))$,$\mathcal{P}_{(m)}= \sum_{k=1}^m\mathcal{P}_k$. Отметим, что $J_1=J$,$\Gamma_1=\Gamma$ и каждый оператор $J_m$ является ортогональным проектором. При этом $J_m$ и $\Gamma_m$ отличаются от соответствующих операторов $J$ и $\Gamma$ на оператор конечного ранга. Поскольку $\mathcal{B}^\alpha\in\mathfrak{L}_{\mathcal{A}}(L_2(0, 1))$, то для него корректно определены операторы $J$ и $\Gamma$ вида (4.2) и операторы $J_m$ и $\Gamma_m$ вида (4.3).
Через $B(L_2(0, 1))$ обозначается пространство линейных ограниченных операторов, действующих в $L_2(0, 1)$. Согласно [24; предположение 3.2], [26; теорема 9] возможность первого преобразования подобия гарантируется следующей леммой.
Лемма 1. Существует достаточно большое не зависящее от $\alpha$ число $m\in \mathbb{N}$ такое, что операторы $\mathcal{B}^\alpha$,$J_m\mathcal{B}^\alpha$ и $\Gamma_m\mathcal{B}^\alpha$ обладают свойствами:
(I) $\Gamma_m\mathcal{B}^\alpha\in B(L_2(0, 1))$ и $\|\Gamma_m\mathcal{B}^\alpha\|_2\leqslant 1/2$;
(II) $(\Gamma_m\mathcal{B}^\alpha)\mathfrak{D}(\mathcal{A})\subset \mathfrak{D}(\mathcal{A})$;
(V) для любого $\delta>0$ существует $\lambda_\delta\in\rho(\mathcal{A})$ такое, что $\|\mathcal{B}^\alpha(\mathcal{A}-\lambda_\delta I)^{-1}\|<\delta$ для всех $\alpha$.
Для доказательства леммы предварительно докажем пару вспомогательных технических утверждений и установим различные оценки для операторов $\Gamma_m\mathcal{B}^\alpha$ и $\mathcal{B}^\alpha\Gamma_m\mathcal{B}^\alpha$.
Оператор $\mathcal{B}^\alpha$ представим в виде $\mathcal{B}^\alpha=\mathcal{B}^\alpha_{1}+\mathcal{B}^\alpha_{2}+\mathcal{B}^\alpha_{3}$, где
Матрицы операторов $\mathcal{B}^\alpha_{\varkappa}$,$\varkappa=1, 2, 3$, обозначим через $(b_{kj}^\varkappa)_{k, j=1}^\infty$, где $b_{kj}^\varkappa=(\mathcal{B}^\alpha_{\varkappa} e_j, e_k)_{L_2(0,1)}$,$\varkappa=1, 2, 3$. Для них верны соотношения
В этих формулах произведения синусов представим в виде суммы произведений синусов и косинусов от аргументов $(j\pm k)x$, а затем проинтегрируем по частям, перенося производную на $S$ и $V$. Это сразу дает оценки
Лемма 2. Операторы $\Gamma_m\mathcal{B}^\alpha$ являются операторами Гильберта–Шмидта, и существует не зависящее от $\alpha$ число $m_0\in\mathbb{N}$ такое, что $\|\Gamma_m\mathcal{B}^\alpha\|_2\leqslant 1/2$ при $m>m_0$. Кроме того, справедливы оценки
Доказательство. Покажем, что $\Gamma \mathcal{B}^\alpha$ является оператором Гильберта–Шмидта. Согласно [28; гл. III, § 9] для этого необходимо доказать сходимость ряда $\sum_{\substack{k, j=1 \\ k\ne j}}^\infty|(\Gamma \mathcal{B}^\alpha e_j, e_k)_{L_2(0, 1)}|^2$. Используя второе представление из (4.2), неравенства (4.6) и определения (4.1), получаем
Следовательно, оператор $\Gamma \mathcal{B}^\alpha$ является оператором Гильберта–Шмидта. Согласно (4.3) оператор $\Gamma_m\mathcal{B}^\alpha$ отличается от $\Gamma\mathcal{B}^\alpha$ на оператор конечного ранга, и потому $\Gamma_m\mathcal{B}^\alpha\in\mathfrak{S}_2(L_2(0, 1))$.
Теперь покажем, что существует число $m\in\mathbb{N}$, не зависящее от $\alpha$, такое, что $\|\Gamma_m\mathcal{B}^\alpha\|_2\leqslant 1/2$. Используя оценки (4.6) и определение (4.3) оператора $\Gamma_m\mathcal{B}^\alpha$, аналогично вычислениям в (4.8) получаем
В правой части полученного неравенства стоит остаток фиксированного ряда $\sum_{\substack{k, j=1\\ k\ne j}}^\infty 1/((k-j)^4(k+j)^2)$, сходимость которого была доказана в (4.8). Поэтому остаток стремится к нулю при $m\to+\infty$, и за счет подходящего выбора числа $m\in\mathbb{N}$ всегда можно добиться выполнения соотношения $\|\Gamma_m\mathcal{B}^\alpha\|_2\leqslant 1/2$ сразу для всех $\alpha$.
Осталось установить неравенства (4.7). Используя определения (4.1) и (4.2), свойство дизъюнктности системы проекторов $\{\mathcal{P}_n\}_{n=1}^\infty$,$n\in\mathbb{N}$, и оценки (4.6), имеем
Полностью повторив приведенные выше рассуждения для $\|(\Gamma\mathcal{B}^\alpha)\mathcal{P}_n\|_2^2$, мы увидим, что эта норма оценивается той же величиной. Тогда из (4.3) непосредственно следует, что
Откуда получаем неравенства (4.7). Лемма доказана.
Лемма 3. Операторы $\mathcal{B}^\alpha\Gamma_m\mathcal{B}^\alpha$ и $(\Gamma_m\mathcal{B}^\alpha)J_m\mathcal{B}^\alpha$ являются операторами Гильберта–Шмидта. Кроме того, имеет место оценка
Доказательство. Операторы $\mathcal{B}^\alpha$ и $J_m\mathcal{B}^\alpha$ очевидно ограничены, и потому с учетом леммы 2 выполнено $\mathcal{B}^\alpha\Gamma_m\mathcal{B}^\alpha, (\Gamma_m\mathcal{B}^\alpha)J_m\mathcal{B}^\alpha \in \mathfrak{S}_2(L_2(0, 1))$.
Из последней оценки и (4.10) следует (4.9). Лемма доказана.
Доказательство леммы 1. Свойства (I) и (III) выполнены в силу лемм 2 и 3. Свойства (II) и (IV) проверяются непосредственным вычислением, повторяя рассуждения из [24; лемма 3.4].
Пусть $\lambda_\delta$ – достаточно большое по модулю отрицательное число. Тогда
Теперь мы можем воспользоваться известными результатами метода подобных операторов.
Теорема 4. Оператор $\mathcal{H}^\alpha$ подобен оператору $\mathcal{A}-J_m\mathcal{B}^\alpha-\widetilde{\mathcal{B}}^\alpha$ с $m>m_0$, где $m_0$ из леммы 2, и
Используя условия (I) и (III) леммы 1, мы можем разложить оператор $(I+\Gamma_m\mathcal{B}^\alpha)^{-1}$ в ряд Неймана, что вместе с формулой (4.13) дает представление (4.11).
Осталось установить оценки из (4.12). Применим к (4.11) проектор $\mathcal{P}_n$,$n\in \mathbb{N}$, справа и с обеих сторон и вычислим разность полученных равенств. Тогда
Подставляя в последнюю формулу неравенства (4.7), (4.9) и $\|\Gamma_m\mathcal{B}^\alpha\|_2\leqslant 1/2$, мы получаем первую оценку из (4.12). Вторая оценка следует непосредственно из (4.8), (4.13) и равномерной по $\alpha$ ограниченности оператора $\mathcal{B}^\alpha$. Теорема 4 доказана.
4.2. Второе преобразование подобия
Согласно теореме 4 оператор $\mathcal{H}^\alpha$ подобен оператору $\mathcal{A}-J_m\mathcal{B}^\alpha-\widetilde{\mathcal{B}}^\alpha$, где $\widetilde{\mathcal{B}}^\alpha$ принадлежит $\mathfrak{S}_2(L_2(0, 1))$. В то же время оператор $J_m\mathcal{B}^\alpha$ только ограничен. Поэтому для корректного применения основных теорем метода подобных операторов (см. [24; § 3]) в качестве невозмущенного оператора мы будем рассматривать оператор $\mathcal{A}-J_m\mathcal{B}^\alpha$, а в качестве возмущения – оператор $\widetilde{\mathcal{B}}^\alpha$.
Теорема 5. Существует число $m_1\in\mathbb{N}$, не зависящее от $\alpha$, такое, что
Доказательство. По построению оператор $J_m$ является проектором, поэтому $\|J_m\|_2=1$. Исходя из [25; лемма 3], леммы 2, а также второй оценки из (4.12), получаем
при $m>m_1$ сразу для всех значений $\alpha$ при подходящем выборе числа $m_1\in\mathbb{N}$. Применяя теперь [26; теорема 7], получаем оставшиеся утверждения теоремы. Теорема доказана.
§ 5. Диагональные элементы матрицы оператора $\mathcal{B}^\alpha\Gamma \mathcal{B}^\alpha$
Для доказательства основных теорем 2 и 3 необходимо провести анализ диагональных элементов матрицы $(c_{nn})_{n\in\mathbb{N}}$ оператора $\mathcal{B}^\alpha\Gamma \mathcal{B}^\alpha$. Они представимы в виде
Стандартным образом решение (5.4) задачи (5.3) можно найти в интегральном виде и учесть условие ортогональности из (5.4). После несложных технических преобразований это решение можно представить в виде
Отсюда, с учетом (5.8), вытекают соотношения (5.7). Лемма доказана.
Данное утверждение позволяет получить асимптотику функции $U$ по $n$ равномерно по $\alpha$ и $x\in[0,1]$. А именно, разбивая интервал интегрирования в формуле для $U_1$ на $(0,x)$ и $(x,1)$ и применяя доказанную лемму поточечно для всех $x\in[0,1]$, в результате рутинных вычислений получаем
Это представление, (5.6) и установленные выше асимптотики для $U_1$,$U_2$ подставим в скалярное произведение $(U, (\mathcal{B}^\alpha)^*e_n)_{L_2(0,1)}$ и воспользуемся затем леммой 4 для вычисления асимптотики полученных интегралов. В полученных соотношениях затем сделаем замену переменной $x-\alpha\mapsto x$ и поменяем местами $x$ и $t$. Тогда
где оценка остатка равномерна по $\alpha$. Замена местами переменных $x$ и $t$ в последнем интеграле в этой формуле позволяет его упростить, что вместе с (5.5) доказывает асимптотику (5.2).
§ 6. Доказательство основных теорем
Доказательство теоремы 2. Из теоремы 5 и [26; лемма 1] следует, что спектр оператора $\mathcal{H}^\alpha$ совпадает со спектром оператора $\mathcal{A}-J_m\mathcal{B}^\alpha-J_m\mathcal{X}_*$ при $m>m_1$ и, кроме того, $\mathcal{A}-J_m\mathcal{B}^\alpha-J_m\mathcal{X}_*$ перестановочен со всеми проекторами $\mathcal{P}_{(m)}$,$\mathcal{P}_n$,$n\geqslant m+1$. Следовательно, подпространства $\mathfrak{H}_{(m)}=\operatorname{Im}\mathcal{P}_{(m)}$,$\mathfrak{H}_n=\operatorname{Im}\mathcal{P}_n$,$n\geqslant m+1$, инвариантны для оператора $\mathcal{A}-J_m\mathcal{B}^\alpha-J_m\mathcal{X}_*$. Имеют место равенства
где $\mathfrak{A}_{(m)}=(\mathcal{A}-J_m\mathcal{B}^\alpha-J_m\mathcal{X}_*)|_{\mathfrak{H}_{(m)}}$ и $\mathfrak{A}_n=(\mathcal{A}-J_m\mathcal{B}^\alpha-J_m\mathcal{X}_*)|_{\mathfrak{H}_n}$ – сужения оператора $\mathcal{A}-J_m\mathcal{B}^\alpha-J_m\mathcal{X}_*$ на подпространства $\mathfrak{H}_{(m)}$ и $\mathfrak{H}_n$,$n\geqslant m+1$, соответственно. Для доказательства равенств (6.1) необходимо в точности повторить рассуждения из [23; замечание 2]. Ввиду конечномерности подпространства $\mathfrak{H}_{(m)}$ из представления (6.1) следует, что множество $\sigma(\mathfrak{A}_{(m)})$ содержит не более чем $m$ элементов. При этом подпространства $\mathfrak{H}_n$ являются одномерными.
Тождества (6.1) показывают, что для вычисления асимптотики собственных значений оператора $\mathcal{H}^\alpha$ для больших $n$ необходимо описать множество $\sigma(\mathfrak{A}_n)$ для больших $n$. Имеет место соотношение
Рассмотрим сужение этого оператора на подпространство $\mathfrak{H}_n$. Тогда $\mathfrak{A}_n=\mathcal{A}_n-(J\mathcal{B}^\alpha)_n-(J\widetilde{\mathcal{B}}^\alpha)_n -(J(\mathcal{X}_*-\widetilde{\mathcal{B}}^\alpha))_n$,$n\geqslant m+1$. Для вычисления всех слагаемых будем использовать матричное представление этих операторов в базисе из $\mathfrak{H}_n$.
Матрица оператора $\mathcal{B}^\alpha$ имеет вид $(b_{kj})_{k, j\in\mathbb{N}}$ с элементами, определенными в (5.1). Из определения (4.11) оператора $\widetilde{\mathcal{B}}^\alpha$ видно, что основной вклад в асимптотику будет вносить оператор $\mathcal{B}^\alpha\Gamma \mathcal{B}^\alpha$, для диагональных элементов $(c_{nn})_{n\in\mathbb{N}}$ которого верна асимптотика (5.2). Остальные слагаемые в (4.11) имеют порядок $O(n^{-2})$ и попадут в остаточный член. Это явно видно, если применить оператор $J_m$ к (4.11), а также использовать неравенства (4.7), (4.9) и очевидное тождество $J_m((\Gamma_m\mathcal{B}^\alpha)J_m\mathcal{B}^\alpha)=0$.
Подробнее остановимся на оценке для слагаемого $(J(\mathcal{X}_*-\widetilde{\mathcal{B}}^\alpha))_n$ для больших $n$. В [25; теорема 8] для такого класса операторов были установлены следующие неравенства:
Таким образом, мы установили все необходимые соотношения для $\mathfrak{A}_n$ при $n\geqslant m+1$. Тогда из приведенных выше рассуждений заключаем, что в представлении (6.1) множество $\sigma(\mathfrak{A}_{(m)})$ является конечным и содержащим не более $m$ элементов, а множества $\sigma(\mathfrak{A}_n)$,$n\geqslant m+1$, являются одноточечными, $\sigma(\mathfrak{A}_n)=\{\lambda_n\}$,$n\geqslant m+1$, причем $\lambda_n$ имеет вид (2.3). Теорема доказана.
Доказательство теоремы 3. В условиях теорем 4 и 5 при $m>m_1$ имеет место равенство
Через $\widetilde{\mathcal{P}}_n$,$n\geqslant m+1$, обозначим проектор Рисса, построенный по множеству $\{\lambda_n\}_{n\geqslant m+1}$, где $\lambda_n$ определяется в теореме 2. Согласно [26; лемма 1] и (6.2) имеет место равенство
Оценим обе части последнего тождества. Из теорем 4 и 5 непосредственно следуют неравенства $\|\Gamma_m\mathcal{B}^\alpha\|_2\leqslant 1/2$,$\|\Gamma_m\mathcal{X}_*\|_2\leqslant 1/2$. Учитывая их, мы выводим
Далее, теорема 5, [24; теорема 3.1] и третье неравенство из (4.12) дают $\|\mathcal{X}_*\|_2\leqslant 4\|\widetilde{\mathcal{B}}^\alpha\|_2\leqslant 4C_{20}$. Аналогичным образом и по той же схеме оценивается $\|\mathcal{P}_n(\Gamma_m\mathcal{X}_*)\|_2$. Таким образом,
Так как система проекторов $\{\mathcal{P}_n\}_{n=1}^\infty$ является ортогональной системой проекторов на одномерное подпространство $\operatorname{Span}\{e_n\}$, то согласно теореме Бари–Маркуса [28; гл. 6, теорема 5.2] система корневых функций оператора $\mathcal{H}^\alpha$ образует базис Бари в пространстве $L_2(0, 1)$. Теорема доказана.
Авторы выражают благодарность рецензенту за многочисленные замечания, которые позволили улучшить первоначальную версию статьи.
Список литературы
1.
Г. Г. Онанов, А. Л. Скубачевский, “Дифференциальные уравнения с отклоняющимися аргументами в стационарных задачах механики деформируемого тела”, Прикл. мех., 15:5 (1979), 30–47; англ. пер.: G. G. Onanov, A. L. Skubachevskii, “Differential equations with displaced arguments in stationary problems in the mechanics of a deformable body”, Soviet Appl. Mech., 15:5 (1979), 391–397
2.
B. Y. Rubinstein, L. M. Pismen, “Resonant two-dimensional patterns in optical cavities with a rotated beam”, Phys. Rev. A, 56:5 (1997), 4264–4272
3.
Self-organization in optical systems and applications in information technology, Springer Ser. Synergetics, 66, eds. M. A. Vorontsov, W. B. Miller, Springer-Verlag, Berlin, 1995, xvi+247 pp.
4.
А. В. Разгулин, “Об автоколебаниях в нелинейной параболической задаче с преобразованным аргументом”, Ж. вычисл. матем. и матем. физ., 33:1 (1993), 69–80; англ. пер.: A. V. Razgulin, “Self-excited oscillations in the nonlinear parabolic problem with transformed argument”, Comput. Math. Math. Phys., 33:1 (1993), 61–70
5.
A. L. Skubachevskii, Elliptic functional differential equations and applications, Oper. Theory Adv. Appl., 91, Birkhäuser Verlag, Basel, 1997, x+293 pp.
6.
А. Л. Скубачевский, “Краевые задачи для эллиптических функционально-дифференциальных уравнений и их приложения”, УМН, 71:5(431) (2016), 3–112; англ. пер.: A. L. Skubachevskii, “Boundary-value problems for elliptic functional-differential equations and their applications”, Russian Math. Surveys, 71:5 (2016), 801–906
7.
D. Neverova, A. Skubachevskii, “On the smoothness of generalized solutions to boundary value problems for strongly elliptic differential-difference equations on a boundary of neighboring subdomains”, Russ. J. Math. Phys., 22:4 (2015), 504–517
8.
A. L. Skubachevskii, “Elliptic differential-difference operators with degeneration and the Kato square root problem”, Math. Nachr., 291:17-18 (2018), 2660–2692
9.
В. В. Лийко, А. Л. Скубачевский, “Смешанные задачи для сильно эллиптических дифференциально-разностных уравнений в цилиндре”, Матем. заметки, 107:5 (2020), 693–716; англ. пер.: V. V. Liiko, A. L. Skubachevskii, “Mixed problems for strongly elliptic differential-difference equations in a cylinder”, Math. Notes, 107:5 (2020), 770–790
10.
G. A. Kamenskii, Extrema of nonlocal functionals and boundary value problems for functional differential equations, Nova Sci. Publ., New York, 2007, x+225 pp.
11.
А. В. Разгулин, “Устойчивость бифуркационных автоколебаний в нелинейной параболической задаче с преобразованным аргументом”, Ж. вычисл. матем. и матем. физ., 33:10 (1993), 1499–1508; англ. пер.: A. V. Razgulin, “The stability of self-excited bifurcation oscillations in a nonlinear parabolic problem with transformed argument”, Comput. Math. Math. Phys., 33:10 (1993), 1323–1330
12.
A. B. Muravnik, “On the Cauchy problem for differential-difference parabolic equations with high-order nonlocal terms of general kind”, Discrete Contin. Dyn. Syst., 16:3 (2006), 541–561
13.
A. Muravnik, “On the half-plane Dirichlet problem for differential-difference elliptic equations with several nonlocal terms”, Math. Model. Nat. Phenom., 12:6 (2017), 130–143
14.
А. Б. Муравник, “Функционально-дифференциальные параболические уравнения: интегральные представления и качественные свойства решений задачи Коши”, Уравнения