Теоретическая и математическая физика
RUS  ENG    ЖУРНАЛЫ   ПЕРСОНАЛИИ   ОРГАНИЗАЦИИ   КОНФЕРЕНЦИИ   СЕМИНАРЫ   ВИДЕОТЕКА   ПАКЕТ AMSBIB  
Общая информация
Последний выпуск
Архив
Импакт-фактор
Правила для авторов
Лицензионный договор
Загрузить рукопись

Поиск публикаций
Поиск ссылок

RSS
Последний выпуск
Текущие выпуски
Архивные выпуски
Что такое RSS



ТМФ:
Год:
Том:
Выпуск:
Страница:
Найти






Персональный вход:
Логин:
Пароль:
Запомнить пароль
Войти
Забыли пароль?
Регистрация


Теоретическая и математическая физика, 2023, том 216, номер 1, страницы 148–168
DOI: https://doi.org/10.4213/tmf10421
(Mi tmf10421)
 

Эта публикация цитируется в 1 научной статье (всего в 1 статье)

Асимптотика решений уравнения Гельмгольца в двухслойной среде с локализованной правой частью

А. Ю. Аникин, А. И. Клевин

Институт проблем механики им. А. Ю. Ишлинского Российской академии наук, Москва, Россия
Список литературы:
Аннотация: Построен главный член квазиклассического асимптотического решения уравнения Гельмгольца с локализованной правой частью, содержащего малый параметр. Данное уравнение возникает, например, в задаче океанической акустики, в которой малый параметр является отношением характерного масштаба “вертикальной” координаты к характерным масштабам остальных координат. Уравнение рассматривается в ограниченной по “вертикальной” координате области, которая разделена на два слоя, между которыми коэффициент в уравнении Гельмгольца и производная решения имеют заданный разрыв в виде скачка. Техника построения асимптотики включает в себя операторное разделение переменных (адиабатическое приближение) и применение к полученным после разделения переменных уравнениям недавно разработанного метода построения асимптотик уравнений с локализованной правой частью.
Ключевые слова: уравнение Гельмгольца, уравнение с правой частью, квазиклассическая асимптотика, канонический оператор Маслова.
Финансовая поддержка Номер гранта
Российский научный фонд 21-11-00341
Исследование выполнено за счет гранта Российского научного фонда № 21-11-00341, https://rscf.ru/project/21-11-00341/.
Поступило в редакцию: 13.12.2022
После доработки: 15.02.2023
Англоязычная версия:
Theoretical and Mathematical Physics, 2023, Volume 216, Issue 1, Pages 1036–1054
DOI: https://doi.org/10.1134/S0040577923070103
Реферативные базы данных:
Тип публикации: Статья

1. Введение

В задачах океанической акустики уравнение Гельмгольца

$$ \begin{equation} \biggl(\Delta_x + \frac{\partial^2}{\partial {z}^2} + \frac{\omega^2}{c^2(x,z)} \biggr) u = - f\biggl( \frac{x-x^0}{\mu}, z \biggr) \end{equation} \tag{1} $$
относительно неизвестной функции $u(x, z)$ описывает (см. [1]–[4]) распространение звука с заданной частотой $\omega$ в мелкой воде в трехмерной области с координатами $(x, z) = (x_1, x_2, z)$ с переменной скоростью звука $c(x, z)$ в ней. В океанической акустике под $z$ понимается координата глубины. Переходя к более общему случаю, мы предполагаем, что уравнение рассматривается в области размерности $n+1$ с координатами $(x, z) = (x_1,\ldots, x_n, z)$. В уравнение входит оператор Лапласа $\Delta_x = \partial^2 /\partial x_1^2 + \cdots + \partial^2 / \partial x_n^2$ по переменным $x$. Пусть имеется локализованный в окрестности множества $x=x^0$ источник, заданный правой частью уравнения (1) с помощью функции $f(\xi, z)$, быстро убывающей при $|\xi| \to +\infty$, и малого параметра $\mu$, отвечающего за степень локализации правой части в окрестности множества $x=x^0$. Задача (1) возникает при изучении решений волнового уравнения со скоростью звука $c(x,z)$, описывающих волны, порожденные гармоническим по времени пространственно локализованным источником.

Предположим, что характерный масштаб задачи $d_0$ по переменной $z$ много меньше характерного масштаба задачи $l_0$ по переменной $x$, т. е. $h\equiv d_0/l_0 \ll 1$. Предположим также, что $\mu \sim d_0$, т. е. параметры $\mu$ и $d_0$ являются величинами одного порядка. Замена функции $f(\xi, z)$ на функцию $f(\xi \mu/d_0, z)$ позволяет без ограничения общности считать, что $\mu = d_0$. Перейдем к безразмерным переменным $z=z' d_0$, $x=x'l_0$ и введем величину $k(x,z)= \omega d_0/c(x,z)$. После очевидных переобозначений мы приходим к уравнению

$$ \begin{equation} -h^2 \Delta_x u - \frac{\partial^2 u}{\partial {z}^2} - k^2(x, z) u = f\biggl( \frac{x-x^0}{h}, z \biggr). \end{equation} \tag{2} $$

Пусть задана область $(x, z) \in \mathbb{R}^n\times [z_-, z_+]$, разделенная на два слоя графиком некоторой гладкой функции $D=D(x)\colon \mathbb{R}^n \to (z_-, z_+)$. Предположим, что коэффициент уравнения (2) имеет вид

$$ \begin{equation*} k(x, z) = \begin{cases} k_+(x), & z\in [D(x), z_+], \\ k_-(x), & z\in [z_-, D(x)), \end{cases} \end{equation*} \notag $$
где $k_\pm(x)$ – некоторые гладкие функции. Пусть также задана функция плотности вида
$$ \begin{equation*} \rho(x, z) = \begin{cases} \rho_+(x), & z\in [D(x), z_+], \\ \rho_-(x), & z\in [z_-, D(x)), \end{cases} \end{equation*} \notag $$
где $\rho_\pm(x)$ – некоторые гладкие функции. Таким образом, функции $k(x, z)$, $\rho(x, z)$ могут иметь разрыв на границе раздела двух сред. Для функции $u(x, z)$ потребуем выполнения следующих граничных условий:
$$ \begin{equation} \begin{gathered} \, u|_{z=z_\pm} = 0, \qquad u|_{z=D(x)-0} = u|_{z=D(x)+0}, \\ \frac{1}{\rho_-(x)}\frac{\partial u}{\partial z}\bigg|_{z=D(x)-0} = \frac{1}{\rho_+(x)}\frac{\partial u}{\partial z}\bigg|_{z=D(x)+0}. \end{gathered} \end{equation} \tag{3} $$
Потребуем, чтобы функция $f(\xi, z)$ принадлежала классу Шварца по $x$ и была непрерывна по $z$.

Замечание 1. В океанической акустике данное разделение на два слоя соответствует среде жидкость–дно. При этом обычно рассматривают полупространство $\mathbb{R}^n\times [z_-, +\infty)$. Возникающая в нашей задаче ситуация с областью $\mathbb{R}^n\times [z_-, z_+]$ может рассматриваться как более простое приближение к задаче в полупространстве. Ограничение глубины дна устанавливается с целью устранения сложностей, связанных с возникновением непрерывного спектра соответствующего оператора в полупространстве (см. [5]).

Замечание 2. Мотивировку рассмотрения уравнения с локализованной правой частью можно найти в работе [6]. В реальных приложениях задача поиска асимптотики уравнения с локализованной правой частью может быть более естественной, чем формально обобщающая ее задача поиска асимптотики функции Грина, которой соответствует правая часть в виде дельта-функции (см., например, [7], [8]). При этом сам переход от асимптотики функции Грина к асимптотике задачи с заданной правой частью может потребовать дополнительных вычислений, связанных, как известно, со взятием свертки функции Грина и правой части.

В настоящей работе мы строим главный член квазиклассического асимптотического решения уравнения (2) с граничными условиями (3) относительно малого положительного параметра $h\to +0$. При этом мы доказываем существование самого асимптотического решения – функции $u(x, z, h)$, подстановка которой в уравнение (2) дает асимптотическое тождество $-h^2 \Delta_x u - u_{zz} - k^2 u = f + O(h^\infty)$.

Метод решения задачи включает в себя три основных направления. Во-первых, учитывая наличие в уравнении “медленной” переменной $z$ и “быстрой” переменной $x$, мы будем использовать адиабатическое приближение в виде операторного разделения переменных (см. [9], [10]), основанного на исчислении Фейнмана–Маслова (см. [11], [12]). Во-вторых, мы будем использовать предложенный в работе [6] (см. также [13]) метод построения асимптотического решения уравнений с локализованной правой частью с помощью теории канонического оператора Маслова на лагранжевых многообразиях. В-третьих, при проведении вычислений, связанных с решениями возникающей в задаче гамильтоновой системы, мы будем использовать некоторые методы квантового аналога принципа Мопертюи–Якоби (см. [14], а также, например, [15]–[17]).

Имея в виду неограниченность области по переменной $x$, отметим, цитируя раздел 1.6 работы [6], что в асимптотической теории сложно наложить естественные условия (типа условий излучения на бесконечности или принципа предельного поглощения), которые бы выделяли единственное асимптотическое решение. Там же отмечается, что построенное решение фактически отвечает принципу предельного поглощения.

Мы не изучаем вопрос о существовании точного решения рассматриваемой задачи и о близости к нему построенных нами асимптотических решений. Даже для дифференциальных уравнений этот вопрос представляет значительную трудность и соответствующие рассуждения и доказательства весьма сложны и громоздки (см., например, [18]). В нашем случае уравнение в $\mathbb{R}^2$, получаемое после адиабатической редукции, является псевдодифференциальным, и соответствующий более сложный вопрос должен составить предмет отдельного исследования.

Отметим, что в работе [19] рассматривалась аналогичная задача в однослойной среде.

2. Операторнозначный символ и его спектр

Согласно [10] разделение “медленной” переменной $z$ и “быстрой” переменной $x$ заключается в рассмотрении оператора в левой части уравнения (2) как дифференциального оператора с операторнозначным символом (см. [20], [9], [10]):

$$ \begin{equation*} \mathcal{H}(x, p) = p^2 - \frac{\partial^2}{\partial z^2} - k^2(x, z), \qquad \widehat{\mathcal{H}} \equiv \mathcal{H}\biggl(\overset{2}{x}, -ih\frac{\overset{1}{ \partial}}{\partial x}\biggr) = -h^2\Delta_x - \frac{\partial^2}{\partial z^2} - k^2(x, z), \end{equation*} \notag $$
где $p=(p_1, \ldots, p_n)$ – сопряженная к $x$ переменная. Уравнение (2) кратко может быть записано как $\widehat{\mathcal{H}}u=f$. Символ $\mathcal{H}(x, p)$ мы понимаем как параметрически зависящий от $(x, p)$ оператор в гильбертовом пространстве функций $\chi(z) \in L_2((z_-, z_+), \rho^{-1}(x, z)\, dz)$ с областью определения
$$ \begin{equation} \begin{aligned} \, \{ \chi(z) &\in L_2((z_-, z_+), \rho^{-1}\, dz) \mid \ \chi, \rho^{-1} \chi' \in AC(z_-, z_+), \\ & \mathcal{H}(x, p) \chi \in L_2((z_-, z_+), \rho^{-1}\, dz), \, \chi|_{z=z_-+0} = \chi|_{z=z_+-0} = 0\}. \end{aligned} \end{equation} \tag{4} $$
Скалярное произведение в $L_2((z_-, z_+), \rho^{-1}(x, z)\, dz)$ задается формулой
$$ \begin{equation*} \langle \chi_1(z), \chi_2(z) \rangle = \frac{1}{\rho_-(x)} \int_{z_-}^{D(x)} \chi_1 \overline{\chi_2} \,dz + \frac{1}{\rho_+(x)} \int_{D(x)}^{z_+} \chi_1 \overline{\chi_2}\, dz. \end{equation*} \notag $$

Зададим следующие вспомогательные объекты. Введем аналитические в пространстве $\mathbb{C}$ и вещественные при $\zeta\in \mathbb{R}$ функции

$$ \begin{equation} \mathbf{s}(\zeta) = \frac{\sin \sqrt{\zeta}}{\sqrt{\zeta}} = \sum_{m=0}^\infty \frac{(-1)^m \zeta^m}{(2m+1)!}, \qquad \mathbf{c}(\zeta) = \cos \sqrt{\zeta} = \sum_{m=0}^\infty \frac{(-1)^m \zeta^m}{(2m)!}, \end{equation} \tag{5} $$
а также аналитическую в $\mathbb{C}$ и строго положительную1 при $\zeta \in \mathbb{R}$ функцию
$$ \begin{equation*} \mathbf{b}(\zeta) = \frac{\sqrt{\zeta} - \sin \sqrt{\zeta}}{\zeta\sqrt{\zeta}} = \sum_{m=0}^\infty \frac{(-1)^m \zeta^k}{(2 m+3)!}. \end{equation*} \notag $$
Положим
$$ \begin{equation} \mu_\pm(x, \mathcal{E}) = k_\pm^2(x) -\mathcal{E}, \qquad \chi_\pm(A_\pm, x, \mathcal{E}, z) = A_\pm \mathbf{s}(\mu_\pm (z-z_\pm)^2) (z-z_\pm), \end{equation} \tag{6} $$
$$ \begin{equation} \chi(A_+, A_-, x, \mathcal{E}, z) = \begin{cases} \chi_+(A_+, x, \mathcal{E}, z), & z\in [D(x), z_+], \\ \chi_-(A_-, x, \mathcal{E}, z), & z\in [z_-, D(x)), \end{cases} \end{equation} \tag{7} $$
$$ \begin{equation} v_+(x, \mathcal{E}) = \begin{pmatrix} \mathbf{s}(\mu_+(D-z_+)^2) (D-z_+) \\ \rho_+^{-1} \mathbf{c}(\mu_+ (D-z_+)^2) \end{pmatrix}, \qquad v_-(x, \mathcal{E}) = \begin{pmatrix} \mathbf{s}(\mu_-(D-z_-)^2) (D-z_-) \\ \rho_-^{-1} \mathbf{c}(\mu_- (D-z_-)^2) \end{pmatrix}, \end{equation} \tag{8} $$
$$ \begin{equation} C_\pm(x, \mathcal{E}) = \pm 2 (z_\pm - D)^3 \mathbf{b}(4\mu_\pm (D-z_\pm)^2), \qquad C_0(x, \mathcal{E}) = \frac{v_-^\mathrm{T} v_+ }{v_-^\mathrm{T} v_-}, \end{equation} \tag{9} $$
$$ \begin{equation} a_+(x, \mathcal{E}) = \frac{1}{\sqrt{\rho_+^{-1} C_+ + \rho_-^{-1}C_- C_0^2}}, \qquad a_-(x, \mathcal{E}) = \frac{C_0}{\sqrt{\rho_+^{-1} C_+ + \rho_-^{-1}C_- C_0^2}}, \nonumber \end{equation} \notag $$
$$ \begin{equation} \begin{aligned} \, [b] \mathcal{D}(x, \mathcal{E}) & = -\rho_+(x) \rho_-(x) \det[v_+(x, \mathcal{E}), v_-(x, \mathcal{E})] ={} \\ &= \rho_+ \mathbf{c}(\mu_- (D-z_-)^2) \mathbf{s}(\mu_+(D-z_+)^2) (z_+-D) + {} \\ & \hphantom{={}} + \rho_- \mathbf{c}(\mu_+ (D-z_+)^2) \mathbf{s}(\mu_-(D-z_-)^2) (D-z_-). \end{aligned} \end{equation} \tag{10} $$
Здесь через $v_-^\mathrm{T}$ обозначено транспонирование зависящего от $(x, \mathcal{E})$ столбца $v_-$. В силу положительности $\mathbf{b}(\zeta)$ при вещественных $\zeta$ функции $C_\pm(x, \mathcal{E})$ удовлетворяют неравенствам $C_\pm(x, \mathcal{E}) > 0$. Легко убедиться, что ни одна из двух векторных функций $v_\pm(x, \mathcal{E})$ не обращается в нуль, что делает функцию $C_0(x, \mathcal{E})$ корректно определенной выражением (9).

Лемма 1. Справедливы равенства

$$ \begin{equation} \begin{gathered} \, \int_{z_-}^{D(x)} |\chi_-|^2\, dz = |A_-|^2 C_-, \qquad \int_{D(x)}^{z_+} |\chi_+|^2\, dz = |A_+|^2 C_+, \\ \langle \chi, \chi \rangle = \frac{|A_-|^2 C_-}{\rho_-} + \frac{|A_+|^2 C_+}{\rho_+}. \end{gathered} \end{equation} \tag{11} $$

Доказательство. Последнее равенство следует из первых двух. Непосредственно проверяется, что

$$ \begin{equation*} |A_\pm|^2 \frac{d}{dz} 2 (z-z_\pm)^3 \mathbf{b}(4\mu_\pm (z-z_\pm)^2) = |\chi_\pm|^2. \end{equation*} \notag $$
Таким образом, найдены первообразные подынтегральных функций. Вычисляя интегралы с помощью первообразных и используя определения функций $C_\pm(x, \mathcal{E})$, через которые выражаются значения первообразных в точке $z=D(x)$, мы получаем доказываемые равенства.

Лемма 2. Пусть пара $(\tilde{x}, \widetilde{\mathcal{E}})$ удовлетворяет равенству $\mathcal{D}(\tilde{x}, \widetilde{\mathcal{E}}) = 0$. Тогда справедливо равенство

$$ \begin{equation*} -\frac{\partial \mathcal{D}}{\partial \mathcal{E}}(\tilde{x}, \widetilde{\mathcal{E}}) = \big[C_0^{-1} \rho_- C_+ + C_0 \rho_+ C_-\big]\big|_{(x, \mathcal{E}) = (\tilde{x}, \widetilde{\mathcal{E}})}. \end{equation*} \notag $$
Как следствие, значение $-(\partial \mathcal{D} / \partial \mathcal{E})(\tilde{x}, \widetilde{\mathcal{E}})$ не обращается в нуль и имеет знак, совпадающий со знаком $C_0(\tilde{x}, \widetilde{\mathcal{E}})$.

Доказательство. В произвольной точке $(x, \mathcal{E})$ производная задается выражением

$$ \begin{equation} -\frac{\partial \mathcal{D}}{\partial \mathcal{E}}(x, \mathcal{E}) = \rho_+ \rho_- \det \biggl[\frac{\partial v_+}{\partial \mathcal{E}}, v_-\biggr] + \rho_+ \rho_- \det \biggl[v_+, \frac{\partial v_-}{\partial \mathcal{E}}\biggr]. \end{equation} \tag{12} $$
При $(x, \mathcal{E}) = (\tilde{x}, \widetilde{\mathcal{E}})$ определитель $\det[v_+, v_-]$ обращается в нуль, т. е. ненулевые векторы $v_\pm(\tilde{x}, \widetilde{\mathcal{E}})$ линейно зависимы. Легко доказывается равенство
$$ \begin{equation*} v_+(\tilde{x}, \widetilde{\mathcal{E}}) = C_0(\tilde{x}, \widetilde{\mathcal{E}}) v_-(\tilde{x}, \widetilde{\mathcal{E}}). \end{equation*} \notag $$
Подстановка его в (12) дает
$$ \begin{equation*} \begin{aligned} \, &-\frac{\partial \mathcal{D}}{\partial \mathcal{E}}(\tilde{x}, \widetilde{\mathcal{E}}) = \biggl\{ \rho_+ \rho_- C_0^{-1} \det \biggl[\frac{\partial v_+}{\partial \mathcal{E}}, v_+\biggr] + \rho_+ \rho_- C_0 \det \biggl[v_-, \frac{\partial v_-}{\partial \mathcal{E}}\biggr] \biggr\} \bigg|_{(x, \mathcal{E}) = (\tilde{x}, \widetilde{\mathcal{E}})} ={}\\ & = \biggl\{ \rho_+ \rho_- C_0^{-1} \det \begin{pmatrix} -\mathbf{s}'(\mu_+(D-z_+)^2) (D-z_+)^3 & \mathbf{s}(\mu_+(D-z_+)^2) (D-z_+) \\ -\rho_+^{-1} \mathbf{c}'(\mu_+ (D-z_+)^2) (D-z_+)^2 & \rho_+^{-1} \mathbf{c}(\mu_+ (D-z_+)^2) \end{pmatrix} +{}\\ &\qquad\! + \rho_+ \rho_- C_0 \det \begin{pmatrix} \mathbf{s}(\mu_-(D-z_-)^2) (D-z_-) & \!\!\!-\mathbf{s}'(\mu_-(D-z_-)^2) (D-z_-)^3 \\ \rho_-^{-1} \mathbf{c}(\mu_- (D-z_-)^2) & \!\!\! -\rho_-^{-1} \mathbf{c}'(\mu_- (D-z_-)^2) (D-z_-)^2 \end{pmatrix} \biggr\} \Bigg|_{ \substack{ x=\tilde{x} \\ \mathcal{E}=\widetilde{\mathcal{E}}} }. \end{aligned} \end{equation*} \notag $$
Использование легко получаемого тождества $2\mathbf{b}(4\zeta) = \mathbf{c}(\zeta) \mathbf{s}'(\zeta) - \mathbf{c}'(\zeta) \mathbf{s}(\zeta)$ дает
$$ \begin{equation*} \begin{aligned} \, -\frac{\partial \mathcal{D}}{\partial \mathcal{E}}(\tilde{x}, \widetilde{\mathcal{E}}) &= \bigl[ -\rho_- C_0^{-1} 2\mathbf{b}(4 \mu_+ (D-z_+)^2) (D-z_+)^3 +{}\\ &\qquad\qquad\qquad\qquad+ \rho_+ C_0 2\mathbf{b}(4 \mu_- (D-z_-)^2) (D-z_-)^3\bigr] \bigr|_{ \begin{subarray}{l} x=\tilde{x} \\ \mathcal{E}=\widetilde{\mathcal{E}} \end{subarray} }. \end{aligned} \end{equation*} \notag $$
Остается воспользоваться данными в (9) определениями функций $C_\pm(x, \mathcal{E})$.

Предложение 1. Оператор $\mathcal{H}(x, p)$ с областью определения (4) является самосопряженным и обладает простым дискретным спектром вида $\{p^2 - \mathcal{E}_\varkappa(x)\}_{\varkappa=1}^\infty$, причем

$$ \begin{equation} \cdots < \mathcal{E}_\varkappa(x) < \cdots < \mathcal{E}_2(x) < \mathcal{E}_1(x), \qquad \lim_{\varkappa\to \infty} \mathcal{E}_\varkappa(x) = -\infty. \end{equation} \tag{13} $$
Значения $\mathcal{E}_\varkappa(x)$ гладко зависят от $x$ и определяются как решения относительно $\mathcal{E}$ уравнения
$$ \begin{equation} \mathcal{D}(x, \mathcal{E}) = 0, \end{equation} \tag{14} $$
называемого дисперсионным соотношением. Последовательности собственных значений $\{p^2 - \mathcal{E}_\varkappa(x)\}_{\varkappa=1}^\infty$ соответствует последовательность собственных функций
$$ \begin{equation} \chi_\varkappa(x, z) = \chi(a_+(x, \mathcal{E}_\varkappa(x)), a_-(x, \mathcal{E}_\varkappa(x)), x, \mathcal{E}_\varkappa(x), z), \qquad \varkappa \in \mathbb{N}, \end{equation} \tag{15} $$
которая составляет ортонормированный базис в $L_2((z_-, z_+), \rho^{-1} dz)$.

Доказательство. Спектральные свойства оператора $\mathcal{H}(x, p)$ следуют из теории операторов Штурма–Лиувилля. Утверждение о самосопряженности оператора $\mathcal{H}(x, p)$ с областью определения (4), о простоте и дискретности его спектра, а также о существовании ортонормированного базиса из его собственных функций можно найти, например, в главе 9 книги [21]. Утверждение о выполнении свойства (13) имеется в теореме 4.3.1 в книге [22]. При этом зависимость оператора и его спектра от $p^2$ легко вычленяется.

При фиксированном $x$ пары $(\mathcal{E}_\varkappa, \chi_\varkappa)$ находятся как решения относительно $(\mathcal{E}, \chi)$ следующей спектральной задачи на отрезке $z\in [z_-, z_+]$:

$$ \begin{equation} \chi'' + k^2(x, z) \chi = \mathcal{E} \chi, \qquad \chi|_{z=z_+} = \chi|_{z=z_-} = 0, \end{equation} \tag{16} $$
$$ \begin{equation} \chi|_{D(x) -0} = \chi|_{D(x)+0}, \qquad \frac{1}{\rho_-(x)} \frac{\partial \chi}{\partial z}\Big|_{z=D(x) - 0} = \frac{1}{\rho_+(x)} \frac{\partial \chi}{\partial z}\Big|_{z=D(x) + 0}. \end{equation} \tag{17} $$
Рассмотрение уравнения (16) с первым начальным условием в каждом слое (т. е. по отдельности при $z\in [D(x), z_+]$ и при $z\in [z_-, D(x))$) приводит к уравнениям с граничными условием с не зависящим от $z$ коэффициентом:
$$ \begin{equation*} \chi_\pm'' + k_\pm^2(x) \chi_\pm = \mathcal{E} \chi_\pm, \qquad \chi_\pm|_{z=z_\pm} = 0. \end{equation*} \notag $$
Общие решения этих уравнений с учетом граничных условий имеют следующий вид (где $A_+, A_- \in \mathbb{C}$ – произвольные постоянные): Введенные в (6) выражения для функций $\chi_\pm(A_\pm, x, \mathcal{E}, z)$, которые зависят от параметров $A_\pm$, $x$, $\mathcal{E}$, являются единообразной записью данных решений.

Функции $\chi_\pm(A_\pm, x, \mathcal{E}, z)$ задают согласно соотношениям (7) кусочным образом функцию $\chi(A_+, A_-, x, \mathcal{E}, z)$. При этом на границе $z=D(x)$ должны выполняться условия (17). Заданные выражением (8) векторные функции $v_+(x, \mathcal{E})$, $v_-(x, \mathcal{E})$ удовлетворяют тождествам

$$ \begin{equation} A_\pm v_\pm(x, \mathcal{E}) = \begin{pmatrix} \chi_\pm(A_\pm, x, \mathcal{E}, D(x)) \\ (\rho_\pm(x))^{-1}(\partial \chi_\pm/\partial z)(A_\pm, x, \mathcal{E}, D(x)) \end{pmatrix}. \end{equation} \tag{18} $$
С помощью этих векторных функций условия (17) могут быть записаны в виде равенства
$$ \begin{equation*} A_+ v_+(x, \mathcal{E}) = A_- v_-(x, \mathcal{E}). \end{equation*} \notag $$
Требование отличия от нулевой функции решения спектральной задачи (16), (17) влечет одновременное выполнение неравенств $A_+\neq 0$ и $A_- \neq 0$. Условие $A_+ v_+ = A_- v_-$ для некоторых $A_\pm$ есть требование линейной зависимости $v_+$, $v_-$, которое эквивалентно условию $\det[v_+, v_-] = 0$. Введенное в (14) дисперсионное соотношение $\mathcal{D}(x, \mathcal{E}) = 0$ эквивалентно уравнению $\det[v_+, v_-] = 0$.

Пусть $(\tilde{x}, \widetilde{\mathcal{E}})$ удовлетворяет дисперсионному соотношению (14), т. е. $\mathcal{D}(\tilde{x}, \widetilde{\mathcal{E}}) = 0$. Подстановка этой пары в (18) дает систему уравнений $A_+ v_+(\tilde{x}, \widetilde{\mathcal{E}}) = A_- v_-(\tilde{x}, \widetilde{\mathcal{E}}) $ относительно $A_\pm$. Данной системе удовлетворяют в точности такие пары $(A_+, A_-)$, что $A_- = C_0(\tilde{x}, \widetilde{\mathcal{E}}) A_+$. Равенствами

$$ \begin{equation*} A_+ = a_+(\tilde{x}, \widetilde{\mathcal{E}}), \qquad A_- = a_-(\tilde{x}, \widetilde{\mathcal{E}}) \end{equation*} \notag $$
задается значение пары $(A_+, A_-)$, удовлетворяющее условию $A_- = C_0(\tilde{x}, \widetilde{\mathcal{E}}) A_+$ и условию нормировки $\langle \chi, \chi\rangle = 1$. Это можно увидеть, воспользовавшись последней формулой в (11).

Таким образом, решениям $\{\mathcal{E}_\varkappa(x)\}_{\varkappa = 1}^\infty$ уравнения (14) (дисперсионного соотношения) соответствуют определенные в (15) вещественные нормированные собственные функции $\chi_\varkappa(x, z)$, параметрически зависящие от $x$. При этом из леммы 2 следует, что

$$ \begin{equation} -\frac{\partial \mathcal{D}}{\partial \mathcal{E}}(x, \mathcal{E}_\varkappa(x)) = [C_0^{-1} \rho_- C_+ + C_0 \rho_+ C_-]|_{(x, \mathcal{E}) = (x, \mathcal{E}_\varkappa(x))} \neq 0 \end{equation} \tag{19} $$
при любом $x$. Следовательно, согласно теореме о неявной функции заключаем, что $\mathcal{E}_\varkappa(x)$ гладким образом зависит от $x$.

Замечание 3. Каждая пара $(x, \mathcal{E})$, удовлетворяющая дисперсионному соотношению $\mathcal{D}(x, \mathcal{E}) =0$, должна при этом удовлетворять условию $\mu_+(x, \mathcal{E}) > 0$ или условию $\mu_-(x, \mathcal{E}) > 0$. Действительно, при одновременном выполнении $\mu_-(x, \mathcal{E}) \leqslant 0$ и $\mu_+(x, \mathcal{E}) \leqslant 0$ из разложений (5) видно, что каждый множитель обоих слагаемых в правой части последнего равенства в (10) является строго положительным, откуда следует $\mathcal{D}(x, \mathcal{E}) > 0$.

Замечание 4. Выражения для собственных функций (15) определяются явно заданной функцией $\chi(a_+(x, \mathcal{E}), a_-(x, \mathcal{E}), x, \mathcal{E}, z)$, в аргументы которой подставляются функции $\mathcal{E}_\varkappa(x)$, $\varkappa \in \mathbb{N}$, неявно заданные дисперсионным соотношением (14).

3. Операторное разделение переменных

Воспользуемся адиабатическим приближением в виде операторного разделения переменных (см. [9], [10]). Метод заключается в представлении неизвестного решения уравнения (2) в виде суммы (разложение решения по модам)

$$ \begin{equation} u(x, z, h) = \sum_{\varkappa=1}^\infty \phi_\varkappa\biggl(\overset{2}{x}, -ih \frac{\overset{1}{\partial}}{\partial x},z, h\biggr) u_\varkappa(x, h), \end{equation} \tag{20} $$
члены которой определяются действиями некоторых псевдодифференциальных операторов $\bigl\{\phi_\varkappa(\overset{2}{x}, \overset{1}{-ih \, \partial / \partial x},z, h)\bigr\}_{\varkappa=1}^\infty$ на не зависящие от переменной $z$ новые неизвестные функции $\{u_\varkappa(x, h)\}_{\varkappa=1}^\infty$. Функции $\{u_\varkappa(x, h)\}_{\varkappa=1}^\infty$ должны удовлетворять при этом некоторым псевдодифференциальным уравнениям, в которые переменная $z$ уже не входит. Отметим, что нахождение псевдодифференциальных операторов $\bigl\{\phi_\varkappa(\overset{2}{x}, \overset{1}{-ih \,\partial / \partial x},z, h)\bigr\}_{\varkappa=1}^\infty$ связано (см., например, [10]) с известной подстановкой Пайерлса.

Предложение 2. Найдется счетное семейство пар псевдодифференциальных операторов

$$ \begin{equation*} \widehat{\phi}_\varkappa \equiv \phi_\varkappa\biggl(\overset{2}{x}, -ih \frac{\overset{1}{ \partial}}{ \partial x},z, h\biggr), \qquad \widehat{H}_\varkappa\equiv H_\varkappa\biggl(\overset{2}{x}, -ih\frac{\overset{1}{ \partial}}{ \partial x},h\biggr), \qquad \varkappa \in \mathbb{N}, \end{equation*} \notag $$
с символами в виде скалярных функций2 $\phi_\varkappa(x, p, z, h)$, $H_\varkappa(x, p, h)$, удовлетворяющее равенствам
$$ \begin{equation} \widehat{\mathcal{H}} \widehat{\phi}_\varkappa= \widehat{\phi}_\varkappa \widehat{H}_\varkappa+O(h^\infty), \qquad \varkappa \in \mathbb{N}. \end{equation} \tag{21} $$
При этом последовательность $\{(\widehat{\phi}_\varkappa, \widehat{H}_\varkappa) \}_{\varkappa = 1}^\infty$ можно выбрать такой, что при каждом $\varkappa \in \mathbb{N}$ первые члены разложений
$$ \begin{equation} \phi_\varkappa(x, p, z, h) = \sum_{m=0}^\infty h^m \phi_{\varkappa,m}(x, p, z), \qquad H_\varkappa(x, p, h) = \sum_{m=0}^\infty h^m H_{\varkappa,m}(x, p) \end{equation} \tag{22} $$
по малому параметру $h$ равны
$$ \begin{equation} H_{\varkappa,0}(x, p) =p^2 - \mathcal{E}_\varkappa(x), \qquad \phi_{\varkappa,0}(x, p, z)=\chi_\varkappa(x, z), \end{equation} \tag{23} $$
$$ \begin{equation} H_{\varkappa,1}(x, p) = -2i \sum_{j=1}^n p_j \biggl\langle\chi_\varkappa, \frac{\partial \chi_\varkappa}{\partial x_j} \biggr\rangle. \end{equation} \tag{24} $$

Доказательство. Для каждого $\varkappa \in \mathbb{N}$ подстановка (22) в (21) дает последовательность уравнений

$$ \begin{equation} (\mathcal{H}-H_{\varkappa,0})\phi_{\varkappa,0}=0, \qquad (\mathcal{H}-H_{\varkappa,0})\phi_{\varkappa,m} = H_{\varkappa, m} \phi_{\varkappa,0} + \mathcal{B}_{\varkappa, m}, \qquad m \in \mathbb{N}, \end{equation} \tag{25} $$
где выражение для $\mathcal{B}_{\varkappa, m}(x, p, z)$ задается многочленом от функций $\phi_{\varkappa, m'}(x, p, z)$, $H_{\varkappa, m'}(x, p)$, $0\leqslant m' \leqslant m-1$, от операторнозначной функции $\mathcal{H}(x, p)$ и от производных этих функций по $(x, p)$ различных порядков. Последовательность пар собственных значений и соответствующих нормированных собственных функций $\{(p^2-\mathcal{E}_\varkappa(x), \chi_\varkappa(x, z))\}_{\varkappa=1}^\infty$ оператора $\mathcal{H}(x, p)$ задает решения первого уравнения. Равенствами (23) соответствующим образом заданы главные члены символов $\phi_\varkappa$, $H_\varkappa$. Остальные уравнения решаются последовательно. На каждом шаге неизвестная функция $H_{\varkappa,m}$ может быть выражена с помощью скалярного умножения уравнения на функцию $\phi_{\varkappa,0} = \chi_\varkappa$, принадлежащую ядру самосопряженного оператора $\mathcal{H}(x, p) - (p^2 - \mathcal{E}_\varkappa(x))$, что дает
$$ \begin{equation} H_{\varkappa,m}(x, p) = -\langle \mathcal{B}_{\varkappa, m}(x, p, z), \chi_\varkappa(x, z) \rangle. \end{equation} \tag{26} $$
Функции $\phi_{\varkappa, m}$ определяются своими разложениями по базису $\{\chi_{\varkappa'} \}_{\varkappa' = 1}^\infty$, которые получаются после скалярного умножения (25) на базисные элементы:
$$ \begin{equation*} \langle \phi_{\varkappa, m}, \chi_{\varkappa} \rangle = 0, \qquad \langle \phi_{\varkappa, m}, \chi_{\varkappa'} \rangle = (H_{\varkappa',0}-H_{\varkappa,0} )^{-1} \langle \mathcal{B}_{\varkappa, m}, \chi_{\varkappa'}\rangle, \qquad \varkappa' \neq \varkappa. \end{equation*} \notag $$
При $m=1$ имеем
$$ \begin{equation*} \mathcal{B}_{\varkappa, 1} =i\sum_{j=1}^n\biggl(\frac{\partial \mathcal{H}}{\partial p_j}\frac{\partial \phi_{\varkappa,0}}{\partial x_j}-\frac{\partial\phi_{\varkappa,0}}{\partial p_j}\frac{\partial H_{\varkappa,0}}{\partial x_j}\biggr) = 2i \sum_{j=1}^n p_j \frac{\partial \phi_{\varkappa,0}}{\partial x_j}. \end{equation*} \notag $$
Подстановка этой функции в (26) дает (24).

Пусть символы $\{(\phi_\varkappa, H_\varkappa) \}_{\varkappa = 1}^\infty$ заданы в соответствии с предложением 2. Зададимся вопросом о существовании функций

$$ \begin{equation} f_\varkappa(\xi, h) = \sum_{m=0}^\infty h^m f_{\varkappa, m}(\xi), \qquad \varkappa \in \mathbb{N}, \end{equation} \tag{27} $$
осуществляющих разложение
$$ \begin{equation} f\biggl( \frac{x-x^0}{h}, z \biggr) = \sum_{\varkappa=1}^\infty \phi_\varkappa\biggl(\overset{2}{x}, -ih\frac{\overset{1}{ \partial}}{ \partial x},z, h\biggr) f_\varkappa\biggl( \frac{x-x^0}{h}, h \biggr) + O(h^\infty). \end{equation} \tag{28} $$
Через $A(p, z)$ обозначим следующим образом определенное преобразование Фурье функции $f(\xi, z)$ по переменной $\xi$:
$$ \begin{equation*} A(p, z) = \frac{e^{-i\pi n/4}}{(2\pi)^{n/2}} \int_{\mathbb{R}^n} e^{-i p \cdot \xi} f(\xi, z)\, d\xi, \qquad f( \xi, z) = \frac{e^{i\pi n/4}}{(2\pi)^{n/2}} \int_{\mathbb{R}^n} e^{i p \cdot \xi} A(p, z)\, dp, \end{equation*} \notag $$
где запись $p \cdot \xi$ означает $p_1 \xi_1 + \cdots + p_n \xi_n$. Введем неизвестные асимптотические ряды
$$ \begin{equation*} A_\varkappa(p, h) = \sum_{m=0}^\infty h^m A_{\varkappa, m}(p), \qquad \varkappa \in \mathbb{N}, \end{equation*} \notag $$
заданные фурье-образами членов рядов (27): при $\varkappa \in \mathbb{N}$, $m\in \mathbb{N}\cup\{0\}$ имеем
$$ \begin{equation*} A_{\varkappa, m}(p) = \frac{e^{-i\pi n/4}}{(2\pi)^{n/2}} \int_{\mathbb{R}^n} e^{-i p \cdot \xi} f_{\varkappa, m}(\xi)\, d\xi, \qquad f_{\varkappa, m}( \xi) = \frac{e^{i\pi n/4}}{(2\pi)^{n/2}} \int_{\mathbb{R}^n} e^{i p \cdot \xi} A_{\varkappa, m}(p)\, dp. \end{equation*} \notag $$

Предложение 3. Найдутся функции $\{f_\varkappa(\xi, h)\}_{\varkappa = 1}^\infty$, осуществляющие разложение (28). Фурье-образы их главных членов задаются равенствами

$$ \begin{equation} A_{\varkappa, 0}(p) = \langle A(p, z), \chi_\varkappa(z, x^0)\rangle, \qquad \varkappa \in \mathbb{N}. \end{equation} \tag{29} $$

Доказательство. Равенство (28) может быть представлено в виде

$$ \begin{equation*} \int_{\mathbb{R}^n} e^{i p \cdot (x-x^0)/h} A(p, z)\, dp = \sum_{\varkappa=1}^\infty \int_{\mathbb{R}^n} e^{i p \cdot (x-x^0)/h} \phi_\varkappa(x, p, z, h) A_{\varkappa}(p, h)\, dp + O(h^\infty). \end{equation*} \notag $$
Воспользуемся разложением в ряд Тейлора и, используя равенство
$$ \begin{equation*} (x-x^0)e^{i p \cdot (x-x^0)/h} = -ih\frac{\partial}{\partial p}e^{i p \cdot (x-x^0)/h}, \end{equation*} \notag $$
проинтегрируем по частям достаточное число раз каждый член суммы в правой части:
$$ \begin{equation*} \begin{aligned} \, \int_{\mathbb{R}^n} e^{i p \cdot (x-x^0)/h} &\phi_\varkappa A_{\varkappa}\, dp ={} \\ &= \int_{\mathbb{R}^n} e^{i p \cdot (x-x^0)/h} \sum_{|\alpha|=0}^\infty \frac{(x-x^0)^\alpha}{\alpha!} \frac{\partial^{|\alpha|} \phi_\varkappa}{\partial x^\alpha}(x^0, p, z, h) A_{\varkappa}(p, h)\, dp ={} \\ &= \sum_{|\alpha|=0}^\infty (ih)^{|\alpha|} \int_{\mathbb{R}^n} e^{i p \cdot (x-x^0)/h} \frac{\partial^{|\alpha|}}{\partial p^\alpha} \biggl[ \frac{\partial^{|\alpha|} \phi_\varkappa}{\partial x^\alpha}(x^0, p, z, h) A_{\varkappa}(p, h) \biggr] dp. \end{aligned} \end{equation*} \notag $$
Приравняем сумму по $\varkappa$ данных асимптотических рядов к $A(p, z)$. В главном члене при нулевой степени $h$ получаем уравнение
$$ \begin{equation*} A(p, z) = \sum_{\varkappa=1}^\infty \phi_{\varkappa, 0}(x^0, p, z) A_{\varkappa, 0}(p). \end{equation*} \notag $$
Это уравнение является записью разложения функции $A(p, z)$ по базису, состоящему из элементов $\chi_\varkappa(z, x^0) \equiv \phi_{\varkappa, 0}(x^0, p, z)$, $\varkappa \in \mathbb{N}$. Его решение задается равенствами (29). Действуя по индукции, предположим, что нам известны функции $A_{\varkappa, m-1}(p)$. Приравнивая члены асимптотического ряда при степени $h^m$, $m\in \mathbb{N}$, получаем уравнение вида
$$ \begin{equation*} \mathcal{A}_{\varkappa, m}(p, z) = \sum_{\varkappa=1}^\infty \phi_{\varkappa, 0}(x^0, p, z) A_{\varkappa, m}(p), \end{equation*} \notag $$
где функция $\mathcal{A}_{\varkappa, m}(p, z)$ выражается через производные по $p$ и взятые в точке $x= x^0$ производные по $x$ различных порядков от известных функций $\phi_{\varkappa, m'}$, $A_{\varkappa, m''}$, $0\leqslant m' \leqslant m$, $0\leqslant m'' \leqslant m-1$. Решение этого уравнения задается коэффициентами разложения по базису: $A_{\varkappa, m}(p) = \langle \mathcal{A}_{\varkappa, m}(p, z), \chi(z, x^0)\rangle$. Преобразование Фурье полученного семейства асимптотических рядов $\{A_\varkappa(p, h)\}_{\varkappa=1}^\infty$ дает семейство асимптотических рядов $\{f_\varkappa(\xi, h)\}_{\varkappa=1}^\infty$, обеспечивающее разложение (28).

Подстановка (20) и (28) в уравнение (2) (кратко записываемое как $\widehat{\mathcal{H}}u=f$) сводит задачу об его асимптотическом решении к семейству задач поиска асимптотических решений $u_\varkappa(x, h)$ уравнений

$$ \begin{equation} \widehat H_\varkappa u_\varkappa(x, h)=f_{\varkappa}\biggl( \frac{x-x^0}{h}, h \biggr), \qquad \varkappa \in \mathbb{N}. \end{equation} \tag{30} $$

Замечание 5. Существенное отличие разложений (20) и (28) от обычных разложений по базису состоит в зависимости от операторов $-ih\, \partial / \partial x$ поправок к главным членам $\phi_\varkappa(\overset{2}{x}, \overset{1}{-ih \, \partial / \partial x},z, h)$. В нашем случае главные члены $\phi_{\varkappa,0}(\overset{2}{x}, \overset{1}{-ih \, \partial / \partial x},z)$ не содержат зависимость от $-ih \,\partial / \partial x$, поскольку их символы $\phi_{\varkappa,0}(x, p, z) \equiv \chi_\varkappa(x, z)$ не зависят от $p$. В общей ситуации уже в главные члены могут входить операторы $-ih\, \partial / \partial x$.

4. “Гиперболические” и “эллиптические” моды. Решение редуцированных уравнений с локализованной правой частью

Для каждого фиксированного $\varkappa$ будем строить решения уравнения (30) с помощью теории канонического оператора Маслова на лагранжевых многообразиях в соответствии с работой [6]. Обозначим через

$$ \begin{equation*} \mathbf{H}(x,p)\equiv H_{\varkappa,0}=p^2-\mathcal{E}_\varkappa(x) \end{equation*} \notag $$
главный символ оператора $\widehat H_\varkappa$, а через
$$ \begin{equation*} \mathbf{H}_{\rm sub}(x, p)\equiv H_{\varkappa, 1}+\frac{i}{2}\operatorname{tr} \frac{\partial^2 H_{\varkappa,0}}{\partial p\,\partial x} = H_{\varkappa, 1} \end{equation*} \notag $$
– его субглавный символ.

Потребуем, чтобы $\mathcal{E}_\varkappa(x^0) \neq 0$ при каждом $\varkappa$. В зависимости от знака $\mathcal{E}_\varkappa(x^0)$ возможны следующие ситуации. Мы говорим об “эллиптической” моде с номером $\varkappa$, если $\mathcal{E}_\varkappa(x^0) < 0$. В этом случае $\{p\mid p^2 -\mathcal{E}_\varkappa(x^0)\} = \varnothing$, откуда следует (см. [6]), что уравнение (30) имеет асимптотическое решение, которое на каждом не содержащем точку $x^0$ компакте в $\mathbb{R}^n$ асимптотически мало: $u_\varkappa \in O(h^\infty)$ в данном компакте. В случае $\mathcal{E}_\varkappa(x^0) > 0$ мы говорим о “гиперболической” моде. Из дискретности последовательности $\{\mathcal{E}_\varkappa\}_{\varkappa=1}^\infty$ и ее монотонного убывания на $-\infty$ следует, что количество гиперболических мод конечно. Далее мы предполагаем, что $\mathcal{E}_\varkappa(x^0) > 0$ при данном выбранном $\varkappa$.

Обозначим через $x=X(q,\tau)$, $p=P(q,\tau)$ решение задачи Коши

$$ \begin{equation} \frac{d x}{d\tau}=\frac{\partial \mathbf{H}}{\partial p}, \qquad \frac{d p}{d\tau}=-\frac{\partial \mathbf{H}}{\partial x}, \qquad x|_{\tau=0}=x^0, \qquad p|_{\tau=0}=q \end{equation} \tag{31} $$
для системы Гамильтона в фазовом пространстве $\mathbb{R}^{2n}$ с гамильтонианом $\mathbf{H}$ и с зависящим от параметра $q\in \mathbb{R}^n$ начальным условием.

Введем сферу $\{x=x^0, \ \mathbf{H}(x^0, p)=0 \} \cong \mathbb{S}^{n-1}$ радиуса $\sqrt{\mathcal{E}_\varkappa(x^0)}$ в слое фазового пространства над точкой $x=x^0$. Зададим некоторую ее параметризацию вида3 $(x, p) = (x^0, Q(\psi))$, $\psi \in \mathbb{S}^{n-1}$. Через $\mathcal{X}(\psi,\tau)=X(Q(\psi),\tau)$ и $\mathcal{P}(\psi,\tau)=P(Q(\psi),\tau)$ обозначим компоненты решения системы Гамильтона (31) с начальными условиями на данной сфере. Возникает $n$-мерное лагранжево многообразие

$$ \begin{equation*} \Lambda_+ = \{x=\mathcal{X}(\psi, \tau), \ p=\mathcal{P}(\psi, \tau)\mid \psi \in \mathbb{S}^{n-1}, \ \tau>0\} \cong \mathbb{S}^{n-1} \times \mathbb{R}_{>0} \end{equation*} \notag $$
в фазовом пространстве $\mathbb{R}^{2n}$.

Условие

$$ \begin{equation*} dq_1\wedge \cdots \wedge dq_n = d\sigma_0 \wedge dH(x^0, q) \end{equation*} \notag $$
однозначно определяет некоторую форму объема $d\sigma_0$ на определенной выше сфере $\mathbb{S}^{n-1}$ в слое $(x, p) = (x^0, q)$, $q\in \mathbb{R}^n$, фазового пространства над точкой $x^0$. Гамильтоновой системой форма $d\sigma_0$ степени $n-1$ переносится на лагранжево многообразие $\Lambda_+$: ее координатная запись в координатах $(\psi, \tau)$ совпадает с записью формы объема $d\sigma_0$ на $\mathbb{S}^{n-1}$ в координатах $\psi$. Лагранжево многообразие следует оснастить следующей гладкой инвариантной относительно фазового потока системы Гамильтона мерой (формой объема): $d\mu_+ = -d\sigma_0 \wedge d\tau$.

В слое над $x^0$ имеем $dH(x^0, q) = 2q_1dq_1 + \cdots + 2q_ndq_n$, и форма $d\sigma_0$ оказывается пропорциональной стандартной форме объема на сфере $\mathbb{S}^{n-1}$ в $\mathbb{R}^n$. Она может быть записана как ограничение на $\mathbb{S}^{n-1}$ формы

$$ \begin{equation*} \sum_{j=1}^n \frac{(-1)^{n+j} q_j}{2q^2} dq_1\wedge \cdots \wedge \widehat{dq_j} \wedge \cdots \wedge dq_n, \end{equation*} \notag $$
где запись $\widehat{dq_j}$ обозначает пропуск соответствующего элемента в произведении. При $n=2$, если мы возьмем $Q(\psi)=\sqrt{\mathcal{E}_\varkappa(x^0)}(\cos\psi,\sin\psi)$, получим $d\sigma_0 = (-1/2)d\psi$ и соответственно $d\mu_+ = (1/2)d\psi \wedge d\tau$.

Следуя работе [6], потребуем выполнения следующих условий:

Замечание 6. Входящая в гамильтониан $\mathbf{H}(x,p)=p^2-\mathcal{E}_\varkappa(x)$ функция $\mathcal{E}_\varkappa(x)$ неявно задана дисперсионным соотношением (14). В практических вычислениях задачу Коши (31) удобнее представить в следующем эквивалентном виде:

$$ \begin{equation} \begin{gathered} \, \frac{d x}{d\tau}=2p, \qquad \frac{d p}{d\tau}= \biggl\{\bigl[C_0^{-1} \rho_- C_+ + C_0 \rho_+ C_-\bigr]^{-1} \frac{\partial \mathcal{D}}{\partial x} \biggr\} \Big|_{\mathcal{E} = p^2}, \\ x\big|_{\tau=0}=x^0, \qquad p\big|_{\tau=0}=q. \end{gathered} \end{equation} \tag{32} $$
Для величин $C_0$, $C_+$, $C_-$, $\mathcal{D}$ даны явные выражения от переменных $(x, \mathcal{E})$, функции $\rho_\pm$, $k_\pm$ зависят от переменной $x$. Выражение в правой части равенства для $d p/ d \tau$ получается следующим образом. Оно равно производной от неявной функции, которая равна отношению
$$ \begin{equation*} \frac{\partial \mathcal{E}_\varkappa}{\partial x}(x, \mathcal{E}_\varkappa(x)) = \frac{\partial \mathcal{D} / \partial x}{- \partial \mathcal{D} / \partial \mathcal{E}}(x, \mathcal{E}_\varkappa(x)). \end{equation*} \notag $$
При этом производная в знаменателе не обращается в нуль и вычислена в (19). Наконец, следует воспользоваться справедливым на траекториях гамильтоновой системы равенством $\mathbf{H}= 0$ и подставить $p^2$ вместо $\mathcal{E}_\varkappa(x)$. Таким образом, правые части системы дифференциальных уравнений задаются явными выражениями, зависящими от переменных $x$, $p$. При этом для определения функций $\mathcal{X}(\psi,\tau)$, $\mathcal{P}(\psi,\tau)$ не обязательно решать уравнение (31) во всех точках $x$. Достаточно знать его решение $\mathcal{E}_\varkappa(x^0)$ только в точке $x=x^0$. Данное решение определяет нули гамильтониана $\mathbf{H}(x,p)$ над точкой $x^0$, которые берутся в качестве начальных условий гамильтоновой системы при задании функций $\mathcal{X}(\psi,\tau)$, $\mathcal{P}(\psi,\tau)$. Приведенные рассуждения могут рассматриваться в рамках общего принципа Мопертюи–Якоби (см. [14]–[17]).

Положим

$$ \begin{equation} \begin{gathered} \, \notag S(\psi,\tau)=2\int_0^\tau \mathcal{P}^2(\psi,\tau')\,d\tau', \\ B(\psi,\tau)=A_{\varkappa, 0}(Q(\psi))\exp\biggl(-i\int_0^\tau \mathbf{H}_{\rm sub}(\mathcal{X}(\psi,\tau'), \mathcal{P}(\psi,\tau'))\,d\tau'\biggr). \end{gathered} \end{equation} \tag{33} $$

Пусть $K_{\Lambda_+, d\mu_+}^h$ – канонический оператор на лагранжевом многообразии с мерой $d\mu_+$.

Предложение 4 [6]. Найдется асимптотическое решение уравнения (30) вида

$$ \begin{equation*} u_\varkappa(x, h) = \sqrt{2\pi}e^{\pi i/4} h^{(n-1)/2} K_{\Lambda_+, d\mu_+}^h [B(\psi,\tau)+O(h)]. \end{equation*} \notag $$

Главный член данного асимптотического решения равен

$$ \begin{equation*} u_\varkappa(x, h) \asymp \sqrt{2\pi}e^{\pi i/4} h^{(n-1)/2} K_{\Lambda_+, d\mu_+}^h [B(\psi,\tau)]. \end{equation*} \notag $$

В данные для определения канонического оператора $K_{\Lambda_+, d\mu_+}^h$ входит выбор ветви аргумента не обращающегося в нуль при $\tau\geqslant 0$, $\psi\in \mathbb{S}^1$ якобиана $J^1(\psi, \tau)$, который может быть рассмотрен как элемент $\varepsilon =1$ семейства с параметром $\varepsilon \in [0, 1]$ якобианов4

$$ \begin{equation*} J^\varepsilon(\psi,\tau)=\det\frac{d(\mathcal{X}_1 - i\varepsilon \mathcal{P}_1 )\wedge \cdots \wedge d(\mathcal{X}_n - i\varepsilon \mathcal{P}_n )}{d\mu_+}. \end{equation*} \notag $$
Согласно [6] этот выбор должен быть осуществлен следующим образом. Введем функцию (здесь вектор $q$ рассматривается в виде столбца)
$$ \begin{equation*} \begin{aligned} \, W(q,\tau,\lambda,\nu,\varepsilon)&=\det\frac{\partial (X-i\varepsilon P, \lambda\tau + i\nu(q^2 + \mathcal{E}_\varkappa(x^0))}{\partial (q,\tau)} ={} \\ &= \det \begin{pmatrix} \dfrac{\partial X}{\partial q}- i\varepsilon \dfrac{\partial P}{\partial q} & \dfrac{\partial X}{\partial \tau}- i\varepsilon \dfrac{\partial P}{\partial \tau} \\ 2i\nu q^\mathrm{T} & \lambda \end{pmatrix}. \end{aligned} \end{equation*} \notag $$
Зафиксируем ветвь аргумента функции $W$ такую, что $\operatorname{arg} W(q, 0, 1, 0, 1) = -\pi n/2$. Справедливо равенство $J^\varepsilon(\psi,\tau) = -i W(Q(\psi), \tau, 0, 1, \varepsilon)$, которое может быть использовано в качестве более удобного определения $J^\varepsilon$ без использования формы $d\mu_+$. Функция $J^\varepsilon$ не обращается в нуль при $\varepsilon > 0$ и может обращаться в нуль при $\varepsilon = 0$. Выбор ветви аргумента5 $J^\varepsilon$ осуществляется исходя из условия $\operatorname{arg} J^\varepsilon(\psi, \tau) = \operatorname{arg} W(Q(\psi), \tau, 0, 1, \varepsilon) - \pi/2$. В частности, данное условие определяет выбор аргумента $J^1$.

Предложение 5. Выбранная выше ветвь $\operatorname{arg} J^1(\psi,\tau)$ при $\tau = 0$ удовлетворяет неравенствам $-\pi n/2 -\pi < \operatorname{arg} J^1(\psi, 0) < -\pi n/2$.

Доказательство. Перенесем ветвь $\operatorname{arg} W(q, 0, \lambda, \nu, 1)$ из точки $(1, 0)$ в точку $(0, 1)$ в плоскости $(\lambda, \nu)$ вдоль отрезка между этими точками. Имеем

$$ \begin{equation*} \begin{aligned} \, W(q,0,\lambda,\nu,1)&= \det \begin{pmatrix} - i E_n & 2q- i \dfrac{\partial \mathcal{E}_\varkappa}{\partial x}(x^0) \\ 2i\nu q^\mathrm{T} & \lambda \end{pmatrix} ={} \\ &= \det[-iE_n] \biggl\{ \lambda - i2\nu q^\mathrm{T}[-iE_n]^{-1} \biggl(2q - i \frac{\partial \mathcal{E}_\varkappa}{\partial x}(x^0)\biggr) \biggr\} ={} \\ &= (-i)^n \biggl\{ \lambda +4 \nu q^2 -i2\nu q^\mathrm{T} \frac{\partial \mathcal{E}_\varkappa}{\partial x}(x^0) \biggr\}, \end{aligned} \end{equation*} \notag $$
где через $E_n$ обозначена единичная матрица размера $n$. Утверждение теоремы следует из того факта, что выражение в фигурных скобках при изменении $(\lambda, \nu)$ остается в правой полуплоскости комплексной плоскости.

При вычислении канонического оператора может быть полезно следующее утверждение.

Предложение 6. При малых положительных значениях $\tau$ справедливо неравенство $J^0(\psi,\tau) < 0$, причем $\arg J^0(\psi,\tau) = -\pi$.

Доказательство. Нам потребуются следующие оценки:

$$ \begin{equation*} \begin{gathered} \, X = x^0 + 2\tau q +O(\tau^2), \qquad P = q+\tau \frac{\partial \mathcal{E}_\varkappa}{\partial x}(x^0) +O(\tau^2), \\ \begin{aligned} \, \frac{\partial X}{\partial \tau} -i\varepsilon \frac{\partial P}{\partial \tau} &= 2q -i\varepsilon \frac{\partial \mathcal{E}_\varkappa}{\partial x}(x^0) + O(\tau), \\ \frac{\partial X}{\partial q} -i\varepsilon \frac{\partial P}{\partial q} &= 2\tau E_n -i\varepsilon E_n + O(\tau^2) = (2\tau - i\varepsilon)\biggl[E_n + O\biggl(\frac{\tau^2}{2\tau -i\varepsilon}\biggr)\biggr]. \end{aligned} \end{gathered} \end{equation*} \notag $$
Имеем
$$ \begin{equation*} \begin{aligned} \, W(q, \tau, 0, 1, \varepsilon) &= \det\Big[ \frac{\partial X}{\partial q}- i\varepsilon \frac{\partial P}{\partial q} \Big] \Big\{ -i2q^\mathrm{T} \Big[ \frac{\partial X}{\partial q}- i\varepsilon \frac{\partial P}{\partial q} \Big]^{-1} \biggl( \frac{\partial X}{\partial \tau}- i\varepsilon \frac{\partial P}{\partial \tau} \biggr) \Big\} ={} \\ &= (2\tau - i\varepsilon)^n \biggl(1 + O\biggl(\frac{\tau^2}{2\tau -i\varepsilon}\biggr)\biggr)\times{} \\ &\qquad\times \biggl\{ \frac{-i2q^\mathrm{T}}{2\tau - i\varepsilon} \biggl[E_n + O\biggl(\frac{\tau^2}{2\tau -i\varepsilon}\biggr)\biggr] \biggl( 2q -i\varepsilon \frac{\partial \mathcal{E}_\varkappa}{\partial x}(x^0) + O(\tau) \biggr) \biggr\} ={} \\ &= (2\tau - i\varepsilon)^n \biggl(1 + O\biggl(\frac{\tau^2}{2\tau -i\varepsilon}\biggr)\biggr)\times{} \\ &\qquad\times \biggl\{ \frac{-i2q^\mathrm{T}}{2\tau - i\varepsilon} \biggl[E_n + O\biggl(\frac{\tau^2}{2\tau -i\varepsilon}\biggr)\biggr] \biggl( 2q -i\varepsilon \frac{\partial \mathcal{E}_\varkappa}{\partial x}(x^0) + O(\tau) \biggr) \biggr\} ={} \\ &= (2\tau - i\varepsilon)^{n-1} \biggl\{ -i4q^2 - 2\varepsilon q^\mathrm{T}\frac{\partial \mathcal{E}_\varkappa}{\partial x}(x^0) + O(\tau) + O\biggl(\frac{\tau^2}{2\tau -i\varepsilon}\biggr) \biggr\}. \end{aligned} \end{equation*} \notag $$
По определению функция $J^0$ вещественна. Из данной формулы видно, что
$$ \begin{equation*} J^0(\psi,\tau) = -i W(Q(\psi), \tau, 0, 1, 0) < 0 \end{equation*} \notag $$
при малых $\tau > 0$. Перенесем ветвь $\operatorname{arg}J^\varepsilon(\psi, \tau) = \operatorname{arg} W(Q(\psi), \tau, 0, 1, \varepsilon) - \pi/2$ из точки $\varepsilon = 1$ в точку $\varepsilon = 0$. При малых $\tau > 0$ аргумент $(2\tau - i\varepsilon)^{n-1}$ получает приращение, близкое к $\pi(n-1)/2$, а приращение аргумента выражения в фигурных скобках, имеющего при всех $\varepsilon \in [0, 1]$ отрицательную мнимую часть, по модулю не превышает $\pi$. В совокупности с предложением 5 это влечет $-2\pi -\pi/2 < \arg J^0(\psi, \tau) < \pi/2$. С учетом отрицательности $J^0$ из этого следует $\arg J^0(\psi,\tau) = -\pi$. Предложение доказано.

Напомним, что действие канонического оператора Маслова определяется сшивкой его локальных представлений в картах некоторого атласа6 $\Lambda_+ = \bigcup_{\beta}V_\beta$ на лагранжевом многообразии. В неособых картах7 результатом действия канонического оператора Маслова является функция типа ВКБ. Выпишем выражение для канонического оператора в произвольной неособой карте $V_{\beta}$. В неособой карте можно однозначно разрешить уравнение $x=\mathcal{X}(\psi,\tau)$. Обозначим его решение через $\psi=\psi_\beta(x)$, $\tau=\tau_\beta(x)$. Пусть $\vartheta_\beta(\psi, \tau)$ – некоторая гладкая функция с носителем в карте $V_{\beta}$. Определим в каждой карте “фазу”, “якобиан” и “амплитуду” как следующие функции от $x$:

$$ \begin{equation*} \begin{gathered} \, S_\beta(x) = S(\psi_\beta(x), \tau_\beta(x)), \qquad J_\beta(x) = J^0(\psi_\beta(x), \tau_\beta(x)), \\ B_\beta(x) = [\vartheta_\beta B](\psi_\beta(x), \tau_\beta(x)). \end{gathered} \end{equation*} \notag $$

Предложение 7. На функциях с носителем в неособой карте $V_{\beta}$ канонический оператор действует в соответствии с выражением

$$ \begin{equation*} K_{\Lambda_+, d\mu_+}^h[\vartheta_\beta B]=e^{-(i/2)\operatorname{arg} J_\beta} \frac{e^{iS_\beta(x)/h} B_\beta(x)}{\sqrt{|J_\beta(x)|}}. \end{equation*} \notag $$
При этом значение ветви $\operatorname{arg} J_\beta$ постоянно в карте. Оно равно приращению аргумента
$$ \begin{equation*} \operatorname{arg} J_\beta = -\pi + \lim_{\varepsilon\to +0}\Delta_\gamma \operatorname{arg} J^\varepsilon(\psi,\tau) \end{equation*} \notag $$
вдоль любого пути $\gamma$ с началом в точке $\Lambda_+$, имеющей достаточно малую координату $\tau > 0$, и с концом в произвольной точке в $V_\beta$.

Выражение для $\mathbf{H}_{\rm sub}(x, p)$ зависит от неявно заданной уравнением (31) функции $\mathcal{E}_\varkappa(x)$. Применение рассуждений, приведенных в замечании 6, позволяет не решать уравнение (31) во всех точках $x$, а ограничиться поиском решения $\mathcal{E}_\varkappa(x^0)$, которое требуется для определения функций $\mathcal{X}(\psi,\tau)$, $\mathcal{P}(\psi,\tau)$, лишь в точке $x=x^0$. Получим не зависящее от $\mathcal{E}_\varkappa(x)$ выражение функции $\mathbf{H}_{\rm sub}(\mathcal{X}(\psi,\tau), \mathcal{P}(\psi,\tau))$, входящей в формулу (33). Введем функцию

$$ \begin{equation*} \mathbf{F}(x,p,\mathcal{E})=i \sum_{j=1}^n p_j \biggl[ (\rho_-^{-1} - \rho_+^{-1}) \frac{\partial D}{\partial x_j} \biggl(\chi^2 \Big|_{ \begin{subarray}{l} z=D(x) \\ A_\pm = a_\pm(x, \mathcal{E}) \end{subarray} } \biggr) + a_-^2 C_- \frac{\partial \rho_-^{-1}}{\partial x_j}+ a_+^2 C_+ \frac{\partial \rho_+^{-1}}{\partial x_j} \biggr], \end{equation*} \notag $$
где $D=D(x)$, $\rho_\pm^{-1} = 1/\rho_\pm(x)$, $a_\pm = a_\pm(x, \mathcal{E})$, $C_\pm = C_\pm(x, \mathcal{E})$, а функция $\chi = \chi(A_+, A_-, x, \mathcal{E}, z)$ задана в (7).

Лемма 3. Справедливо равенство

$$ \begin{equation} \mathbf{H}_{\rm sub}(x, p)=\mathbf{F}(x,p,\mathcal{E}_\varkappa(x)). \end{equation} \tag{34} $$

Доказательство. Разность данных функций равна

$$ \begin{equation} \begin{aligned} \, \notag & \mathbf{H}_{\rm sub}(x, p) -\mathbf{F}(x,p,\mathcal{E}_\varkappa(x)) = -i \sum_{j=1}^n \biggl[ 2 p_j \biggl\langle\chi_\varkappa, \frac{\partial \chi_\varkappa}{\partial x_j} \biggr\rangle +{} \\ &+ (\rho_-^{-1} - \rho_+^{-1}) \frac{\partial D}{\partial x_j} \chi_\varkappa^2(x, D(x)) + \biggl( a_-^2 C_- \frac{\partial \rho_-^{-1}}{\partial x_j}+ a_+^2 C_+ \frac{\partial \rho_+^{-1}}{\partial x_j}\biggr)\Big|_{\mathcal{E} = \mathcal{E}_\varkappa(x)} \biggr]. \end{aligned} \end{equation} \tag{35} $$

Взятие производной по переменной $x_j$ от обеих частей тождества $\langle \chi_\varkappa, \chi_\varkappa\rangle = 1$ дает тождество

$$ \begin{equation} \begin{aligned} \, 2 \biggl\langle \chi_\varkappa, \frac{\partial \chi_\varkappa}{\partial x_j} \biggr\rangle &+ \chi_\varkappa^2(x, D(x)) (\rho_-^{-1} - \rho_+^{-1}) \frac{\partial D}{\partial x_j} +{} \notag \\ &+ \frac{\partial \rho_-^{-1}}{\partial x_j} \int_{z_-}^{D(x)} \chi_\varkappa^2\, dz + \frac{\partial \rho_+^{-1}}{\partial x_j} \int_{D(x)}^{z_+} \chi^2_\varkappa\, dz = 0. \end{aligned} \end{equation} \tag{36} $$
Формулы для интегралов (11) с учетом выражения (15) для функции $\chi_\varkappa$ дают равенства
$$ \begin{equation*} \int_{z_-}^{D(x)} \chi_\varkappa^2\, dz = a_-^2C_-|_{\mathcal{E} = \mathcal{E}_\varkappa(x)}, \qquad \int_{D(x)}^{z_+} \chi^2_\varkappa\, dz = a_+^2C_+|_{\mathcal{E} = \mathcal{E}_\varkappa(x)}. \end{equation*} \notag $$
Подстановка этих равенств в левую часть (36) дает разность (35), которая, таким образом, оказывается равной нулю.

Предложение 8. Справедливо равенство

$$ \begin{equation} \mathbf{H}_{\rm sub}(\mathcal{X}(\psi,\tau), \mathcal{P}(\psi,\tau))= \mathbf{F}(x,p, p^2) \Big|_{ \begin{subarray}{l} x=\mathcal{X}(\psi,\tau) \\ p=\mathcal{P}(\psi,\tau) \end{subarray} }. \end{equation} \tag{37} $$

Доказательство. Данное равенство получается с помощью подстановки в (34) функций $\mathcal{P}(\psi,\tau)$, $\mathcal{X}(\psi,\tau)$ и использования тождества $\mathcal{E}_\varkappa(\mathcal{X}(\psi,\tau)) \equiv \mathcal{P}^2(\psi,\tau)$.

5. Алгоритм построения асимптотики

Приведем алгоритм построения главного члена асимптотического решения уравнения (2).

6. Пример

Пусть $n=2$ и правая часть уравнения имеет вид

$$ \begin{equation} f(\xi, z) = F(\xi) G(z), \qquad F(\xi) = c_1 e^{-(\xi-c_2)^2 / c_3^2}, \qquad G(z) = c_4 e^{-(z-c_5)^2 / c_6^2}. \end{equation} \tag{38} $$
Тогда
$$ \begin{equation*} A(p, z) = G(z) \frac{-i c_1 c_3}{4\pi} e^{-i p c_2} e^{-c_3 p^2/4}. \end{equation*} \notag $$
Для каждой моды $(\mathcal{E}_\varkappa, \chi_\varkappa)$ введем числа
$$ \begin{equation*} A_{\pm,0} = a_\pm(x^0, \mathcal{E}_\varkappa(x^0)), \qquad \mu_{\pm,0} = \mu_\pm(x^0, \mathcal{E}_\varkappa(x^0)) \end{equation*} \notag $$
и функции
$$ \begin{equation*} \begin{aligned} \, \mathbf{P}_\pm(z) ={}& \frac{-i \sqrt{\pi} A_{\pm,0} c_4 c_6}{4\mu_{\pm,0}} \exp\biggl[-i(c_5+z_\pm) \sqrt{\mu_{\pm,0}} - \frac{c_6^2 \mu_{\pm,0}}{4}\biggr] \times {} \\ &\times \biggl\{ \exp[2iz_\pm \sqrt{\mu_{\pm,0}}] \operatorname{erf} \biggl[ \frac{c_5 - z}{c_6} - \frac{i c_6 \sqrt{\mu_{\pm,0}}}{2} \biggr] -{} \\ &\qquad- \exp[2ic_5 \sqrt{\mu_{\pm,0}}] \operatorname{erf} \biggl[ \frac{c_5 - z}{c_6} + \frac{i c_6 \sqrt{\mu_{\pm,0}}}{2} \biggr] \biggr\}, \end{aligned} \end{equation*} \notag $$
в которые входит функция ошибок $\operatorname{erf}(\zeta)$. Справедливо равенство $(d/dz)\mathbf{P}_\pm(z) = G(z) \chi_\pm(x^0, z)$. Следовательно,
$$ \begin{equation*} \begin{aligned} \, A_{\varkappa,0}(p) &= \langle A(p, z), \chi_\varkappa(x^0, z) \rangle = \frac{-i c_1 c_3}{4\pi} e^{-i p c_2} e^{- c_3 p^2/4} \times{} \\ &\hphantom{={}}\times \biggl[\frac{\mathbf{P}_-(D(x^0)) - \mathbf{P}_-(z_-)}{\rho_+(x^0)} + \frac{\mathbf{P}_+(z_+) - \mathbf{P}_+(D(x^0))}{\rho_-(x^0)}\biggr]. \end{aligned} \end{equation*} \notag $$

Рассмотрим следущий пример:

$$ \begin{equation*} \begin{gathered} \, k_+(x)=1, \qquad k_-(x) = \frac{3}{4}, \qquad \rho_+(x) = 1, \qquad \rho_-(x) = 2, \qquad z_+ = 0, \qquad z_- = -5, \\ x^0 = (-3, 0), \qquad D(x) = -1-3 e^{-x_1^2/2}, \end{gathered} \end{equation*} \notag $$
и возьмем правую часть вида (38) с константами
$$ \begin{equation*} c_1 = -i, \qquad c_2 = (0, 0), \qquad c_3 = \sqrt{2}, \qquad c_4 = 1, \qquad c_5 = -1, \qquad c_6 = \frac{1}{\sqrt{8}}. \end{equation*} \notag $$
В точке $x=x^0$ дисперсионное соотношение (14) имеет только одно положительное решение $\mathcal{E}_1(x^0) \approx 0.112918$. На рис. 1 изображен график функции $D(x)$. Для единственной гиперболической моды $\varkappa = 1$ на рис. 2 изображены проекции на пространство $\mathbb{R}^n_x$ координатных линий на лагранжевом многообразии $\Lambda_+$, образованных равномерной прямоугольной сеткой в плоскости координат $(\psi, \tau)$. На рис. 3 изображен график функции $\operatorname{Re} u_\varkappa(x, h)$, $\varkappa=1$, при $h=0.1$. На рис. 4 изображен график вещественной части главного члена асимпотики $\operatorname{Re} u(x, z, h) = \operatorname{Re} [u_\varkappa(x, h) \chi_\varkappa(x, z)]$, $\varkappa=1$, в плоскости $x_2=0$ при $h=0.1$.

Благодарности

Авторы благодарны С. Ю. Доброхотову и В. Е. Назайкинскому за внимание к работе и ценные дискуссии.

Конфликт интересов

Авторы заявляют, что у них нет конфликта интересов.

Список литературы

1. B. Katsnelson, V. Petnikov, J. Lynch, Fundamentals of Shallow Water Acoustics, Springer, New York, 2012  crossref
2. F. B. Jensen, W. A. Kuperman, M. B. Porter, H. Schmidt, Computational Ocean Acoustics, Springer, New York, 2011  crossref  mathscinet
3. G. V. Frisk, Ocean and Seabed Acoustics: A Theory of Wave Propagation, Prentice Hall PTR, Englewood Cliffs, NJ, 1998
4. C. L. Pekeris, “Theory of propagation of explosive sound in shallow water”, Propagation of Sound in the Ocean, Geological Society of America Memoirs, 27, eds. J. L. Worzel; M. Ewing, C. L. Pekeris, Geol. Soc. Amer., New York, 1948, 1–118  crossref
5. P. S. Petrov, M. Yu. Trofimov, A. D. Zakharenko, “Modal perturbation theory in the case of bathymetry variations in shallow-water acoustics”, Russ. J. Math. Phys., 28:2 (2021), 257–262  crossref  mathscinet
6. А. Ю. Аникин, С. Ю. Доброхотов, В. Е. Назайкинский, М. Руло, “Лагранжевы многообразия и конструкция асимптотик для (псевдо)дифференциальных уравнений с локализованными правыми частями”, ТМФ, 214:1 (2023), 3–29  mathnet  crossref  crossref
7. В. М. Бабич, “О коротковолновой асимптотике функции Грина для уравнения Гельмгольца”, Матем. сб., 65(107):4 (1964), 576–630  mathnet  mathscinet  zmath
8. В. М. Бабич, В. С. Булдырев, Асимптотические методы в задачах дифракции коротких волн, Наука, М., 1972  mathscinet
9. В. В. Белов, С. Ю. Доброхотов, В. П. Маслов, Т. Я. Тудоровский, “Обобщенный адиабатический принцип для описания динамики электрона в искривленных наноструктурах”, УФН, 175:9 (2005), 1004–1010  mathnet  crossref  crossref  adsnasa
10. V. V. Belov, S. Yu. Dobrokhotov, T. Ya. Tudorovskiy, “Operator separation of variables for adiabatic problems in quantum and wave mechanics”, J. Eng. Math., 55:1–4 (2006), 183–237  crossref  mathscinet
11. В. П. Маслов, Операторные методы, Наука, М., 1973  mathscinet  mathscinet  zmath  zmath
12. R. P. Feynman, “An operator calculus having applications in quantum electrodynamics”, Phys. Rev., 84:1 (1951), 108–128  crossref  mathscinet
13. А. Ю. Аникин, С. Ю. Доброхотов, В. Е. Назайкинский, М. Руло, “Канонический оператор Маслова на паре лагранжевых многообразий и асимптотика решений стационарных уравнений с локализованными правыми частями”, Докл. РАН, 475:6 (2017), 624–628  crossref  crossref  zmath
14. С. Ю. Доброхотов, Д. С. Миненков, М. Руло, “Принцип Мопертюи–Якоби для гамильтонианов вида $f(x,|p|)$ в некоторых двумерных стационарных квазиклассических задачах”, Матем. заметки, 97:1 (2015), 48–57  mathnet  crossref  crossref  mathscinet  zmath
15. С. Ю. Доброхотов, М. Руло, “Квазиклассический аналог принципа Мопертюи–Якоби и его приложение к линейной теории волн на воде”, Матем. заметки, 87:3 (2010), 458–463  mathnet  crossref  crossref  mathscinet  zmath
16. S. Yu. Dobrokhotov, M. Rouleux, “The semi-classical Maupertuis–Jacobi correspondence for quasi-periodic hamiltonian flows with applications to linear water waves theory”, Asymptotic Anal., 74:1–2 (2011), 33–73  crossref  mathscinet
17. А. Ю. Аникин, С. Ю. Доброхотов, А. И. Клевин, Б. Тироцци, “Скаляризация стационарных квазиклассических задач для систем уравнений и приложение к физике плазмы”, ТМФ, 193:3 (2017), 409–433  mathnet  crossref  crossref  adsnasa
18. Б. Р. Вайнберг, “О коротковолновой асимптотике решений стационарных задач и асимптотике при $t\to\infty$ решений нестационарных задач”, УМН, 30:2(182) (1975), 3–55  mathnet  crossref  mathscinet  zmath
19. П. Н. Петров, С. Ю. Доброхотов, “Асимптотика решения уравнения Гельмгольца в трехмерном слое переменной толщины с локализованной правой частью”, Ж. вычисл. матем. и матем. физ., 59:4 (2019), 566–578  mathnet  crossref  crossref  mathscinet  zmath  adsnasa
20. В. П. Маслов, Теория возмущений и асимптотические методы, Изд-во Моск. ун-та, M., 1965  zmath
21. G. Teschl, Mathematical Methods in Quantum Mechanics, Graduate Studies in Mathematics, 157, AMS, Providence, RI, 2014  crossref  mathscinet
22. A. Zettl, Sturm–Liouville Theory, Mathematical Surveys and Monographs, 121, AMS, Providence, RI, 2005  crossref  mathscinet
23. В. П. Маслов, М. В. Федорюк, Квазиклассическое приближение для уравнений квантовой механики, Наука, М., 1976  mathscinet  zmath
24. С. Ю. Доброхотов, В. Е. Назайкинский, А. И. Шафаревич, “Новые интегральные представления канонического оператора Маслова в особых картах”, Изв. РАН. Сер. матем., 81:2 (2017), 53–96  mathnet  crossref  crossref  mathscinet  adsnasa
25. S. Yu. Dobrokhotov, V. E. Nazaikinskii, “Efficient formulas for the Maslov canonical operator near a simple caustic”, Russ. J. Math. Phys., 25:4 (2018), 545–552  crossref  mathscinet
26. С. Ю. Доброхотов, В. Е. Назайкинский, “Лагранжевы многообразия и эффективные формулы для коротковолновых асимптотик в окрестности точки возврата каустики”, Матем. заметки, 108:3 (2020), 334–359  mathnet  crossref  crossref  mathscinet

Образец цитирования: А. Ю. Аникин, А. И. Клевин, “Асимптотика решений уравнения Гельмгольца в двухслойной среде с локализованной правой частью”, ТМФ, 216:1 (2023), 148–168; Theoret. and Math. Phys., 216:1 (2023), 1036–1054
Цитирование в формате AMSBIB
\RBibitem{AniKle23}
\by А.~Ю.~Аникин, А.~И.~Клевин
\paper Асимптотика решений уравнения Гельмгольца в двухслойной среде с локализованной правой частью
\jour ТМФ
\yr 2023
\vol 216
\issue 1
\pages 148--168
\mathnet{http://mi.mathnet.ru/tmf10421}
\crossref{https://doi.org/10.4213/tmf10421}
\mathscinet{http://mathscinet.ams.org/mathscinet-getitem?mr=4619872}
\adsnasa{https://adsabs.harvard.edu/cgi-bin/bib_query?2023TMP...216.1036A}
\transl
\jour Theoret. and Math. Phys.
\yr 2023
\vol 216
\issue 1
\pages 1036--1054
\crossref{https://doi.org/10.1134/S0040577923070103}
\scopus{https://www.scopus.com/record/display.url?origin=inward&eid=2-s2.0-85165561345}
Образцы ссылок на эту страницу:
  • https://www.mathnet.ru/rus/tmf10421
  • https://doi.org/10.4213/tmf10421
  • https://www.mathnet.ru/rus/tmf/v216/i1/p148
  • Эта публикация цитируется в следующих 1 статьяx:
    Citing articles in Google Scholar: Russian citations, English citations
    Related articles in Google Scholar: Russian articles, English articles
    Теоретическая и математическая физика Theoretical and Mathematical Physics
    Статистика просмотров:
    Страница аннотации:152
    PDF полного текста:13
    HTML русской версии:81
    Список литературы:22
    Первая страница:12
     
      Обратная связь:
     Пользовательское соглашение  Регистрация посетителей портала  Логотипы © Математический институт им. В. А. Стеклова РАН, 2024