Аннотация:
Пусть f0 и f∞ – формальные степенные ряды в нуле и в бесконечности соответственно, и пусть Pn/Qn, deg(Pn),deg(Qn)≤n, – рациональная функция, которая одновременно интерполирует f0 в нуле с порядком n и f∞ в бесконечности с порядком n+1. Если f0,f∞ – ростки многозначных функций с конечных числом точек ветвления, то (как было показано В. И. Буслаевым) существует единственный компакт F, в дополнении к которому такие рациональные аппроксимации сходятся по емкости к приближаемым функциям. Множество F может разбивать или не разбивать плоскость. Мы изучаем равномерную сходимость аппроксимаций для геометрически простейшего случая множеств F, которые разбивают плоскость.
Библиография: 26 названий.
Исследование выполнено при поддержке Simons Foundation (грант № CGM-354538) и Московского центра фундаментальной и прикладной математики грантом Министерства науки и высшего образования РФ (соглашение № 075-15-2019-1623).
Пусть f(z)=∑∞k=0fkz−k – формальный степенной ряд в бесконечности, и пусть Mn/Nn – такая рациональная функция, что deg(Mn),deg(Nn)⩽n и
(Nnf−Mn)(z)=O(z−n−1)при z→∞.
Известно, что такая рациональная функция единственна и называется классической диагональной аппроксимацией Паде для f в бесконечности. Следующая теорема1[x]1В § 1 приводятся лишь усеченные формулировки теорем, которые достаточны для наших нужд. суммирует один из наиболее фундаментальных вкладов Герберта Шталя в теорию комплексных приближений; см. [18]–[23].
Теорема 1.1 (Г. Шталь). Предположим, что росток f в бесконечности можно аналитически продолжить вдоль любого пути в ¯C, не проходящего через точки фиксированного полярного множества2[x]2Необходимые понятия теории потенциала можно найти в [17]., и что продолжения этого ростка по различным путям хотя бы в одну точку вне упомянутого полярного множества дают хотя бы два различных ростка (т.е. в эту точку росток продолжается не однозначно).
Тогда существует компакт F такой, что:
(i)F не разбивает плоскость и f имеет голоморфное однозначное продолжение в область D:=¯C∖F;
(ii)F состоит из открытых непересекающихся аналитических дуг Ji, их концевых точек и некоторых особых точек функции f, и 3[x]3Такие множества F теперь называются симметричными компактами или S-компактами.
∂gF(z,∞)∂n+=∂gF(z,∞)∂n−
в каждой точке z∈⋃iJi, где gF(⋅,∞) – функция Грина области D с особенностью в бесконечности, а ∂/∂n± обозначают односторонние нормальные производные;
(iii) для любого компакта V⊂D при любом ϵ>0 выполняется соотношение
lim
где \mathrm{cap}(\cdot) – логарифмическая емкость.
Более общо, если выбрать ветвь f, голоморфную в окрестности некоторого замкнутого множества, а также набор 2n+1 не обязательно различных точек этого множества, то можно определить многоточечные аппроксимации Паде, интерполирующие f в этих точках. Аналог теоремы Шталя для многоточечных аппроксимаций Паде был доказан А. А. Гончаром и Е. А. Рахмановым (см. [13]). Однако существование множества F, удовлетворяющего условию (i) и соответствующему обобщению условия (ii), в [13] не показывалось, а предполагалось; из существования такого множества затем выводилось утверждение (iii) (в [16], [2], [25] есть результаты о существовании таких компактов, симметричных во внешнем поле).
Компакты, симметричные во внешнем поле, характеризуются тем, что они минимизируют определенную логарифмическую емкость во внешнем поле. Одно из главных препятствий, мешающих доказать общую теорему об их существовании, состоит в том, что минимизирующий компакт может разбивать плоскость. В [5] В. И. Буслаев интерпретировал эту возможность не как помеху, а как важное свойство. Более точно, пусть f_0 и f_\infty – формальные степенные ряды в нуле и в бесконечности соответственно, т.е.
Рациональная функция P_n/Q_n есть двухточечная аппроксимация Паде типа (n_1,n_2), n_1+n_2=2n+1, пары (f_0,f_\infty), если \deg(P_n),\deg(Q_n)\leqslant n и
Как и в случае классических аппроксимаций Паде, отношение P_n/Q_n всегда единственно. В [5; теорема 1] (см. также [6] для обобщений на m-точечные аппроксимации Паде) В. И. Буслаев доказал следующее.
Теорема 1.2 (В. И. Буслаев). Пусть ростки f_0 и f_\infty в (1.1) можно аналитически продолжить вдоль любого пути в \overline{\mathbb{C}}, не проходящего через некоторое конечное множество точек, причем одна из этих точек является точкой ветвления для f_0, а еще одна – точкой ветвления для f_\infty.
Тогда существует компакт F такой, что:
(i)\overline{\mathbb{C}}\setminus F =D_0\cup D_\infty, где области D_0\ni 0 и D_\infty \ni\infty либо не пересекаются, либо совпадают, и f_e имеет однозначное голоморфное продолжение в D_e, e\in\{0,\infty\} (если D_0=D_\infty=D, то f_0 и f_\infty являются аналитическими продолжениями друг друга в области D);
(ii)F состоит из открытых аналитических дуг и их концевых точек, и в каждой точке этих дуг выполняется равенство
где g_F(z,0) – функция Грина области D_0 с особенностью в нуле, а g_F(z,\infty) – функция Грина области D_\infty с особенностью в бесконечности;
(iii) если последовательность индексов n_1+n_2=2n+1 в (1.2) такова, что \lim_{n\to\infty}n_i/n=1, i\in\{1,2\}, и если \partial D_0\not\subset\partial D_\infty и \partial D_\infty\not\subset\partial D_0, то для любого компакта V \subset {\mathbb{C}}\setminus F и любого \epsilon>0
\begin{equation*}
\lim_{n\to\infty}\mathrm{cap}\biggl\{z\in V \colon \biggl|\biggl(f-\frac{P_n}{Q_n}\biggr)(z)\biggr|^{1/n}\geqslant e^{-g_F(z,0)-g_F(z,\infty)}+\epsilon\biggr\} = 0,
\end{equation*}
\notag
где \mathrm{cap}(\cdot) – логарифмическая емкость и f=f_e в D_e, e\in\{0,\infty\}.
Пусть f – росток в бесконечности, который можно аналитически продолжить вдоль любого пути в \overline{\mathbb{C}}, не проходящего через конечное множество точек, причем одна из них является точкой ветвления для некоторого продолжения ростка. Положим f_\infty = f, и пусть f_0 – один из аналитических элементов f в нуле. Результат Буслаева утверждает, что независимо от выбора f_0 двухточечные аппроксимации Паде всегда сходятся в приближаемой паре ростков. Однако в зависимости от выбора f_0 множество сходимости может состоять из одной или двух связных компонент.
Утверждение (iii) теоремы Буслаева описывает так называемую слабую асимптотику, или асимптотику корня n-й степени аппроксимаций Паде. Цель настоящей работы – получить формулы сильной асимптотики для подходящих подпоследовательностей индексов. В случае, когда F не разбивает плоскость, подобные результаты были получены в [15], [24], [3], [4], [25] при различных предположениях о приближаемой функции. Ниже мы изучаем простейший геометрический случай, для которого области D_0 и D_{\infty} различны (рис. 1).
§ 2. Основные результаты
Пусть \mathbb{D} – открытый единичный круг с центром в нуле. Для заданной точки a\in\mathbb{D}\setminus\{0\} обозначим через Kкомпакт Чеботарёва для тройки точек \{-1,1,\mathfrak j(a)\}, где \mathfrak j(z):=(z+z^{-1})/2 – функция Жуковского. Компакт K состоит из критических траекторий квадратичного дифференциала
\begin{equation*}
-\frac{(\zeta-\mathfrak j(b))\,\mathrm d \zeta^2}{(\zeta^2-1)(\zeta-\mathfrak j(a))},
\end{equation*}
\notag
где точка b\in\mathbb{D} определяется однозначно (точка \mathfrak j(b) называется точкой Чеботарёва). Если \operatorname{Im}(a)\neq 0, то множество K состоит из трех аналитических дуг, выходящих из \mathfrak j(b) и заканчивающихся в каждой из точек -1,1,\mathfrak j(a). Если a\in(-1,0), то b=-1 и K=[-\mathfrak j(a),1], а если a\in (0,1), то b=1 и K=[-1,\mathfrak j(a)]. Компакт Буслаева F, соответствующий точке a, в этом случае задается формулой
\begin{equation*}
F := \{z\colon \mathfrak j(z)\in K\}
\end{equation*}
\notag
(рис. 1–3). Очевидно, F разбивает плоскость на две односвязные компоненты, одна из которых содержит нуль, а вторая – окрестность бесконечности; мы обозначаем их через D_0 и D_\infty соответственно. Введем обозначение
\begin{equation*}
F := F^\circ\cup E, \qquad E:=\{a,a^{-1},b,b^{-1}\},
\end{equation*}
\notag
где F^\circ состоит из четырех (трех, если a вещественно) открытых непересекающихся жордановых дуг; см. рис. 1. Пусть F_{\pm1} (соответственно F_{\pm1}^\circ) – замкнутые (соответственно открытые) жордановы дуги, соединяющие b с b^{-1} и содержащие точки \pm1, а F_{a^{\pm1}} (соответственно F_{a^{\pm1}}^\circ) – замкнутые (соответственно открытые) жордановы дуги, соединяющие a^{\pm1} с b^{\pm1} (при этом F_{-1}=\varnothing при a\in (-1,0) и F_1=\varnothing при a\in(0,1)). Выберем ориентацию этих дуг так, чтобы объединение F_1\cup F_{-1} было замкнутой жордановой кривой, ориентированной против часовой стрелки, а F_{a^{-1}}\cup F_1\cup F_a было жордановой дугой, ориентированной от a^{-1} к a; см. рис. 1. Простой подстановкой можно проверить, что множество F^\circ состоит из критических траекторий квадратичного дифференциала
\begin{equation}
{-}\frac{(z-b)(z-b^{-1})}{(z-a)(z-a^{-1})}\,\frac{\mathrm d z^2}{z^2}.
\end{equation}
\tag{2.1}
Начиная с этого момента, мы фиксируем точку a\in \mathbb{D}\setminus\{0\} и соответствующий компакт Буслаева F. В настоящей работе мы будем изучать приближения пар (1.1) вида
для некоторых классов весов \rho, которые будут указаны позже. Таким образом, в зависимости от ситуации мы будем говорить об аппроксимациях (1.2) либо для пары функций, либо для одной функции, имея в виду в обоих случаях пару (2.2). Как можно увидеть из рис. 1, структура множества F для вещественных и для невещественных значений a качественно различается. Поэтому не удивительно, что описание асимптотики аппроксимаций также будет зависеть от того, является параметр a вещественным или нет.
Напомним, что в этом случае F есть объединение единичной окружности \mathbb{T} и отрезка, соединяющего a и a^{-1} (всегда ориентированного от a^{-1} к a), который мы обозначим через [a^{-1},a]. Начнем с описания класса весов \rho, которые мы рассматриваем. Пусть
обозначает голоморфную ветвь многозначной функции в правой части в дополнении к отрезку [a^{-1},a] так, что w(z) = z + \mathcal O(1) при z\to\infty. Для заданной функции \rho на F нам будет удобно положить4[x]4Для функции g, голоморфной в дополнении к ориентированной жордановой дуге или кривой J, мы обозначаем через g_+ (соответственно g_-) предельные значения функции g на положительной (соответственно отрицательной) стороне J относительно ориентации этой кривой.
где b:=a/|a|. Мы будем говорить, что функция \rho принадлежит классу \mathcal W_1, если h продолжается до голоморфной и не обращающейся в нуль функции в некоторой окрестности множества F\setminus\{a,a^{-1}\} с нулевым приращением аргумента вдоль единичной окружности и если существуют такие вещественные константы \alpha,\beta>-1, что
здесь функция \widetilde h голоморфна и не обращается в нуль в некоторой окрестности отрезка [a^{-1},a], а ветви степенных функций голоморфны в окрестности интервала (a^{-1},a). Нам также будет удобно выделить следующий подкласс класса \mathcal W_1: будем считать, что \rho\in \mathcal W_2 \subset \mathcal W_1, если \alpha=\beta=-1/2, т.е. если h продолжается до голоморфной и не обращающейся в нуль функции в некоторой окрестности компакта F. В частности, если мы определим \rho из соотношения (2.4), положив h(s)\equiv 2, то приближаемая пара функций будет иметь вид (1/w,-1/w).
Для краткости обозначим \sqrt x := \mathrm i \sqrt{|x|} для x<0. Как мы увидим ниже, асимптотическое поведение корня n-й степени из аппроксимаций описывается следующей функцией:
которая голоморфна вне отрезка [a^{-1},a], обращается в нуль в точке z=0 и имеет простой полюс в бесконечности. Кроме того, поскольку w имеет чисто мнимые односторонние предельные значения на отрезке [a^{-1},a] и w^2(e^{\mathrm i t}) = 2e^{\mathrm i t}(\cos(t)-\mathfrak j(a)), то можно легко проверить, что
Для того чтобы лучше описать поведение аппроксимаций, определим следующие функции. Пусть \log h обозначает любую гладкую ветвь логарифма h на единичной окружности (напомним, что функция h непрерывна на \mathbb{T} и ее аргумент имеет нулевое приращение по этой окружности). Положим
\begin{equation}
D(z) := \exp\biggl\{\frac{1}{2\pi\mathrm i}\int_{\mathbb{T}}\frac{\log h(s)\,\mathrm d s}{s-z}\biggr\};
\end{equation}
\tag{2.8}
тогда D(z) – голоморфная и не обращающаяся в нуль функция в дополнении к \mathbb{T}, которая с каждой из сторон гладко продолжается на единичную окружность, причем эти продолжения удовлетворяют соотношению D_+(s)=D_-(s)h(s) (см. [12; гл. I]). Поэтому функция
голоморфна и не обращается в нуль вне отрезка [a^{-1},a], а ее предельные значения с двух сторон этого отрезка удовлетворяют соотношению S_+(s)S_-(s)=\widehat h(s). Используя введенные функции (2.6), (2.8) и (2.9), определим функцию
которая голоморфна, не обращается в нуль в {\mathbb{C}}\setminus F и имеет нуль порядка n в бесконечности (легко видеть, что произведение функций \mathcal R_n и {\mathcal Q}_n кусочно постоянно, и поэтому второе определение может показаться избыточным; однако для невещественных a соотношение между этими двумя функциями не будет столь простым, а мы предпочитаем рассматривать случаи вещественного и невещественного параметра в едином стиле).
Теорема 2.1. Обозначим через P_n/Q_n двухточечную аппроксимацию Паде типа (n,n+1) для (2.2) с \rho\in \mathcal W_1. Положим
\begin{equation}
R_n(z) := z^{-n}(Q_nf_\rho-P_n)(z), \qquad z\in\overline{\mathbb{C}}\setminus F.
\end{equation}
\tag{2.12}
Тогда для всех достаточно больших n полином Q_n имеет степень n и может быть нормирован так, чтобы его старший коэффициент был единицей. При такой нормировке выполняются соотношения
в которых для функций \upsilon=\upsilon_{{\mathcal Q}_n},\upsilon_{\mathcal R_n} локально равномерно в \overline{\mathbb{C}}\setminus F выполняются оценки
с некоторыми константами C>0 и c<1, а постоянные \gamma_n выбираются так, чтобы выполнялось предельное соотношение \lim_{z\to\infty}\gamma_n{\mathcal Q}_n(z)z^{-n} = 1.
Как и в п. 2.1, мы начнем с определения классов весов \rho, которые будем рассматривать.
Мы будем говорить, что функция \rhoпринадлежит классу\mathcal W_1, если выполнены следующие условия:
1) ограничения функции \rho на дуги F_1^\circ,F_{-1}^\circ продолжаются до голоморфных, не обращающихся в нуль функций в окрестности F_1 и F_{-1} соответственно;
2) существуют константы \alpha,\beta>-1 такие, что ограничения \rho на F_a^\circ и F_{a^{-1}}^\circ имеют вид
где \widetilde\rho – голоморфная и не обращающаяся в нуль функция в окрестности F_a или F_{a^{-1}} соответственно, и ветви степенных функций голоморфны в окрестности F_a^\circ и F_{a^{-1}}^\circ;
в некоторой окрестности точки b (верхнее соотношение) и в некоторой окрестности точки b^{-1} (нижнее соотношение), где, допуская некоторую вольность в обозначениях, мы обозначаем через \rho_{|F_e}(z) не только ограничение функции \rho на F_e, но и его аналитическое продолжение.
В отличие от случая вещественного параметра класс \mathcal W_2 теперь не будет пересекаться с классом \mathcal W_1. Пусть
обозначает голоморфную ветвь многозначной функции в правой части равенства, голоморфную в {\mathbb{C}}\setminus(F_{a^{-1}}\cup F_a) и такую, что w(z) = z^2 + \mathcal O(z) при z\to\infty.
Мы будем говорить, что функция \rhoпринадлежит классу\mathcal W_2, если функция
продолжается до голоморфной и не обращающейся в нуль функции в некоторой окрестности компакта F. Легко проверить, что, если h(s)\equiv2, мы опять получаем пару (1/w,-1/w).
Преимущество класса \mathcal W_2 состоит в том, что в аналоге соотношения (2.13) оценки ошибок опять будут экспоненциально убывать. Однако это не очень естественный класс, так как для того, чтобы определить функции из этого класса, нужно знать точку b, а ее нахождение в общем случае является трансцендентной задачей. Класс \mathcal W_1 более естественный; например, он содержит пары (c_0/w_a, c_\infty/w_a) для некоторых констант c_0\neq\pm c_\infty, где w_a(z) := \sqrt{(z-a)(z-a^{-1})} обозначает ветвь радикала, голоморфную вне F_{a^{-1}} \cup F_1 \cup F_a.
Аналоги функций \varphi, D и S из (2.6), (2.8) и (2.9) теперь выглядят сложнее, поскольку в данном случае их нужно определять с помощью различных дифференциалов на римановой поверхности
которая имеет род 1. Более того, этих функций самих по себе не достаточно, чтобы определить аналоги функций {\mathcal Q}_n,\mathcal R_n. Поэтому мы изберем другой подход.
Поверхность \boldsymbol{\mathfrak{R}} можно реализовать как две копии \overline{\mathbb{C}}, разрезанные вдоль F_{a^{-1}}\cup F_a и затем склеенные между собой крест-накрест вдоль соответствующих дуг. Мы будем обозначать эти копии через \boldsymbol{\mathfrak{R}}^{(0)} и \boldsymbol{\mathfrak{R}}^{(1)}. Пусть \pi\colon \boldsymbol{\mathfrak{R}}\to\overline{\mathbb{C}} обозначает естественную проекцию \pi(\boldsymbol z) =z. Для точки z\not\in F_{a^{-1}}\cup F_a мы также будем обозначать через z^{(i)} ее прообраз при проекции, лежащий на \boldsymbol{\mathfrak{R}}^{(i)}. Положим
где дуга \boldsymbol\beta :=\pi^{-1}(F_{a^{-1}}) ориентирована так, что \boldsymbol{\mathfrak{R}}^{(0)}\setminus\boldsymbol\Delta остается справа, когда \boldsymbol\beta проходится в положительном направлении, дуга \boldsymbol\delta:=\pi^{-1}(F_a) ориентирована так, что \boldsymbol{\mathfrak{R}}^{(0)}\setminus\boldsymbol\Delta остается слева, когда \boldsymbol\delta проходится в положительном направлении, а \boldsymbol\alpha:=\pi^{-1}(F_1) и \boldsymbol\gamma:=\pi^{-1}(F_{-1}) ориентированы так, что их положительные направления на \boldsymbol{\mathfrak{R}}^{(0)} совпадают с положительными направлениями дуг F_1 и F_{-1}.
Теорема 2.2. Пусть h(s) – функция на F, для которой существуют такие вещественные константы \alpha(e), e\in E=F\setminus F^\circ, и такие ветви функций (z-e)^{\alpha(e)}, что произведение
продолжается до не обращающейся в нуль гёльдеровой функции на F.
Тогда для любого n\in\mathbb{N} существуют мероморфная на \boldsymbol{\mathfrak{R}}\setminus\boldsymbol\Delta функция \Psi_n и точка \boldsymbol z_n \in \boldsymbol{\mathfrak{R}} такие, что:
(1)\Psi_n имеет непрерывные предельные значения на \pi^{-1}(F^\circ), которые удовлетворяют соотношению
(2)\Psi_n имеет полюс порядка n в точке \infty^{(1)} (порядка n-1, если \boldsymbol z_n\,{=}\,\infty^{(1)}), простой полюс в точке \infty^{(0)} (точку голоморфности, если \boldsymbol z_n=\infty^{(0)}), нуль кратности n в точке 0^{(1)} (кратности n+1, если \boldsymbol z_n=0^{(1)}), простой нуль в точке \boldsymbol z_n\not\in\boldsymbol\Delta\cup\{\infty^{(0)},\infty^{(1)},0^{(1)}\} (если \boldsymbol z_n\in\boldsymbol\Delta, то мы используем явное представление (3.16) для \Psi_n, чтобы рассматривать \boldsymbol z_n как нуль и \Psi_{n+}, и \Psi_{n-}), а в остальных точках \Psi_{n} не обращается в нуль и конечна;
для k\in\{0,1\}, где m_n(e):=1, если \pi(\boldsymbol z_n)=e, и m_n(e):=0 в противном случае5[x]5Запись |g_1(z)|\sim|g_2(z)| при z\to z_0 означает, что существует константа C>1 такая, что C^{-1} \leqslant |(g_1/g_2)(z)| \leqslant C в некоторой окрестности точки z_0..
Обратно, если \Psi – функция с нулем кратности не меньше n в точке 0^{(1)}, с полюсом порядка не выше n в точке \infty^{(1)}, не более чем с простым полюсом в точке \infty^{(0)}, не имеющая других полюсов и удовлетворяющая соотношениям (2.21) и
то \Psi получается из \Psi_n умножением на константу.
Как отмечалось выше, функции \Psi_n можно явно выразить через различные дифференциалы на \boldsymbol{\mathfrak{R}} так же, как и через тета-функции Римана; см. формулу (3.16) ниже.
При нашем асимптотическом анализе мы будем интересоваться только индексами n, для которых точки \boldsymbol z_n отделены от \infty^{(1)}. Как утверждается в следующем предложении, таких индексов бесконечно много.
Предложение 2.1. Для заданного \varepsilon>0 пусть D_\varepsilon(z) – круг радиуса \varepsilon в сферической метрике с центром z\in\overline{\mathbb{C}}, и пусть D_\varepsilon(\boldsymbol z) – связная компонента \pi^{-1}(D_\varepsilon(z)), содержащая \boldsymbol z. Положим
Тогда для всех достаточно малых \varepsilon и для всех натуральных чисел n>1 либо n, либо n-1 принадлежит множеству \mathbb{N}_\varepsilon.
В соответствии с нашими обозначениями для любого множества B положим B^{(i)}:= \pi^{-1}(B)\cap \boldsymbol{\mathfrak{R}}^{(i)}. Теперь мы готовы определить аналоги функций (2.10) и (2.11) для невещественных a. Пусть
– это кусочно голоморфная функция в \overline{\mathbb{C}}\setminus F с нулем кратности n-1 в бесконечности6[x]6Опять же, эти порядки могут измениться в зависимости от положения точки \boldsymbol z_n.. Из-за специфики анализа Римана–Гильберта, который мы будем использовать для изучения поведения аппроксимаций Паде, нам понадобятся еще следующие функции. Пусть \Upsilon_n – рациональная функция на \boldsymbol{\mathfrak{R}}, которая конечна всюду, кроме двух простых полюсов в точках \infty^{(0)} и \boldsymbol z_n, и имеет простой нуль в точке 0^{(1)} (такая функция единственна с точностью до скалярного множителя). Положим
и определим {\mathcal Q}_n^\star и \mathcal R_n^\star соответственно соотношениями (2.24) и (2.25), однако при этом в (2.25) заменим \Psi_n на \Psi_n^\star и z^n на z^{n+1}. Полученные функции голоморфны в {\mathbb{C}}\setminus F, причем {\mathcal Q}_n^\star имеет полюс порядка n в бесконечности, в то время как \mathcal R_n^\star имеет в бесконечности нуль кратности n-1.
Теорема 2.3. Обозначим через P_n/Q_n двухточечную аппроксимацию Паде типа (n,n+1) для пары (2.2) с \rho\in \mathcal W_1 \cup \mathcal W_2. Пусть R_n задается формулой (2.12). Далее, пусть {\mathcal Q}_n и \mathcal R_n задаются формулами (2.24) и (2.25) с функцией \Psi_n, определенной, как в теореме 2.2, но по функции h, задаваемой равенством (2.19).
Тогда для любого \varepsilon> 0 и для всех достаточно больших n\in\mathbb{N}_\varepsilon полином Q_n имеет степень n и может быть нормирован так, что его старший коэффициент равняется единице. При такой нормировке выполняются соотношения
где константа \gamma_n такова, что \lim_{z\to\infty}\gamma_n{\mathcal Q}_n(z)z^{-n} = 1, а функции \upsilon=\upsilon_{nj} локально равномерно в \overline{\mathbb{C}}\setminus F удовлетворяют оценкам (2.14).
Из-за наличия блуждающего нуля \boldsymbol z_n из этой формулы непосредственно не следует локальная равномерная сходимость аппроксимаций к f_\rho. Действительно, если
то {\mathcal Q}_n(z_n)=0, и поэтому аппроксимация имеет полюс в окрестности точки z_n. Следующие результаты способствуют дальнейшему прояснению ситуации.
Теорема 2.4. Функции \Psi_n можно нормировать так, что для любого замкнутого множества B \subset \overline{\mathbb{C}}\setminus F и для любого \delta>0 существуют положительные константы C(B) и C_\delta(B) такие, что
(из принципа минимума для супергармонических функций следует, что 2g(z)-\log|z|\,{>}\,0 в {\mathbb{C}}\setminus F). Более того, константы C(B) и C_\delta(B) можно выбрать так, что
Функции \Upsilon_n можно нормировать так, чтобы неравенства (2.27) и (2.29) выполнялись также для {\mathcal Q}_n^\star и \mathcal R_n^\star. Более того, для любого \delta>0 существует константа C_\delta такая, что
В силу теоремы Буслаева ясно, что 2g(z)-\log|z|=g_F(z,0)+g_F(z,\infty).
Благодарность
Автор выражает благодарность Андрею Мартинес-Финкельштейну за неоднократные ценные обсуждения.
§ 3. Доказательства теорем 2.2, 2.4 и предложения 2.1
Пусть \boldsymbol{\mathfrak{R}} – риманова поверхность, определенная в (2.20). Будем считать, что каждое из множеств \boldsymbol{\mathfrak{R}}^{(i)} – замкнутое подмножество в \boldsymbol{\mathfrak{R}}, т.е. оно содержит циклы \boldsymbol\beta и \boldsymbol\delta. Определим конформную инволюцию поверхности \boldsymbol{\mathfrak{R}} формулой \boldsymbol z=(z,w) \mapsto \boldsymbol z^* =(z,-w). Легко видеть, что пара циклов (\boldsymbol\alpha,\boldsymbol\beta) образует базис гомологий поверхности \boldsymbol{\mathfrak{R}}. В частности, область \boldsymbol{\mathfrak{R}}_{\boldsymbol\alpha,\boldsymbol\beta}:= \boldsymbol{\mathfrak{R}}\setminus(\boldsymbol\alpha\cup\boldsymbol\beta) односвязна.
3.1. Дифференциал Наттолла
Пусть w(\boldsymbol z) := (-1)^iw(z), \boldsymbol z\in \boldsymbol{\mathfrak{R}}^{(i)}, где w(z) – ветвь, определенная в (2.18). Для удобства положим
и v(z):= v(\boldsymbol z) для \boldsymbol z\in\boldsymbol{\mathfrak{R}}^{(0)}\setminus\{\boldsymbol\beta\cup\boldsymbol\gamma\cup\boldsymbol\delta\}. Заметим, что v(0)=v(\infty)=1. Дифференциал
\begin{equation}
\mathcal N(\boldsymbol s) = \frac{1-v(\boldsymbol s)}{2s}\,\mathrm d s
\end{equation}
\tag{3.1}
голоморфен всюду, кроме двух простых полюсов в точках 0^{(1)} и \infty^{(1)} с вычетами соответственно 1 и -1. Более того, можно легко проверить, используя (2.1), что все периоды \mathcal N(\boldsymbol s) чисто мнимые, и поэтому мы можем определить вещественные константы:
\begin{equation}
\omega := -\frac1{2\pi\mathrm i}\oint_{\boldsymbol\beta}\mathcal N, \qquad \tau:=\frac1{2\pi\mathrm i}\oint_{\boldsymbol\alpha}\mathcal N.
\end{equation}
\tag{3.2}
Лемма 3.1. Пусть \sqrt a – главное значение квадратного корня. Положим
где путь интегрирования полностью лежит в \boldsymbol{\mathfrak{R}}_{\boldsymbol\alpha,\boldsymbol\beta}.
Тогда функция \Phi(\boldsymbol z) голоморфна и не обращается в нуль в \boldsymbol{\mathfrak{R}}_{\boldsymbol\alpha,\boldsymbol\beta}\setminus\{\infty^{(1)},0^{(1)}\}, имеет простой полюс в точке \infty^{(1)} и простой нуль в точке 0^{(1)}\in \boldsymbol{\mathfrak{R}}^{(1)}. Выполняется равенство \Phi(\boldsymbol z)\Phi(\boldsymbol z^*)=z, а предельные значения \Phi(\boldsymbol z) удовлетворяют соотношениям
Более того, |\Phi(\boldsymbol z)|^2<|z| при z\in D_{\mathcal R} и |\Phi(\boldsymbol z)|^2>|z| при z\in D_{\mathcal Q}.
Доказательство. Свойства голоморфности \Phi непосредственно следуют из свойств голоморфности \mathcal N. Поскольку v(\boldsymbol z) = -v(z) при \boldsymbol z\in\boldsymbol{\mathfrak{R}}^{(1)}\setminus\{\boldsymbol\beta\cup\boldsymbol\gamma\}, то выполняется соотношение
\begin{equation}
\Phi(\boldsymbol z)\Phi(\boldsymbol z^*) = a\exp\biggl\{\int_a^z\frac{(1-v(t))\,\mathrm d t}{2t}+\int_a^z\frac{(1+v(t))\,\mathrm d t}{2t}\biggr\} = a\exp\biggl\{\int_a^z\frac{\mathrm d t}{t}\biggr\} = z.
\end{equation}
\tag{3.5}
Далее, на \boldsymbol \alpha и \boldsymbol\beta мы получаем соответственно
что дает (3.4); см. (3.2). Наконец, напомним, что F состоит из критических траекторий квадратичного дифференциала -(v(z)\,\mathrm dz/z)^2; см. (2.1). Значит, интеграл от v(t)\,\mathrm d t/t по любой дуге, содержащейся в F, чисто мнимый, и поэтому
где знак минус нужно брать при \boldsymbol s\in\boldsymbol{\mathfrak{R}}^{(0)}, а знак плюс – при \boldsymbol s\in\boldsymbol{\mathfrak{R}}^{(1)}. Теперь последнее утверждение леммы следует из принципа максимума модуля. Лемма доказана.
3.2. Голоморфные дифференциалы
Можно легко проверить, что
\begin{equation*}
\mathcal H(\boldsymbol s) := \frac{C\,\mathrm d s}{w(\boldsymbol s)}, \quad\text{где }\ C:=\biggl(\oint_{\boldsymbol\alpha}\frac{\mathrm d s}{w(\boldsymbol s)}\biggr)^{-1},
\end{equation*}
\notag
является голоморфным дифференциалом на \boldsymbol{\mathfrak{R}} (единственным с точностью до мультипликативной константы). Обозначим
\begin{equation*}
\mathbf B:= \oint_{\boldsymbol \beta}\mathcal H.
\end{equation*}
\notag
Известно, что \operatorname{Im}(\mathbf B) >0. Доказательство следующей леммы полностью аналогично доказательству леммы 3.1.
Лемма 3.2. Для данного \sigma\in{\mathbb{C}} положим
где путь интегрирования полностью лежит в \boldsymbol{\mathfrak{R}}_{\boldsymbol\alpha,\boldsymbol\beta}.
Тогда функция A_\sigma(\boldsymbol z) голоморфна и не обращается в нуль в \boldsymbol{\mathfrak{R}}_{\boldsymbol\alpha,\boldsymbol\beta}. Выполняется равенство A_\sigma(\boldsymbol z) A_\sigma(\boldsymbol z^*)\equiv1, а предельные значения A_\sigma(\boldsymbol z) на \boldsymbol\alpha,\boldsymbol\beta удовлетворяют соотношениям
Обозначим через \mathcal C_{\boldsymbol z} единственный мероморфный дифференциал, который имеет два простых полюса в точках \boldsymbol z и \boldsymbol z^* с вычетами 1 и -1 соответственно и нулевой \boldsymbol\alpha-период. Если \pi(\boldsymbol z)\in{\mathbb{C}}, то можно легко проверить, что
Тогда функция S_h(\boldsymbol z) голоморфна и не обращается в нуль в \boldsymbol{\mathfrak{R}}\setminus\boldsymbol\Delta. Выполняется равенство S_h(\boldsymbol z)S_h(\boldsymbol z^*)\equiv1, а предельные значения функции S_h(\boldsymbol z) удовлетворяют соотношению
Доказательство. Пусть \boldsymbol\nu – симметричный относительно инволюции цикл на \boldsymbol{\mathfrak{R}}, проходящий через точки ветвления \boldsymbol p_1, \boldsymbol p_2, и пусть \lambda_{\boldsymbol\nu} – такая инвариантная относительно инволюции функция на \boldsymbol\nu, что для некоторых вещественных констант \alpha_1, \alpha_2 функция
Поскольку \mathcal C_{\boldsymbol z^*} = - \mathcal C_{\boldsymbol z}, имеем \Lambda_{\boldsymbol\nu}(\boldsymbol z) + \Lambda_{\boldsymbol\nu}(\boldsymbol z^*) \equiv 0. Более того, известно (см. [26; формулы (2.7)–(2.9)]), что \Lambda_{\boldsymbol\nu}(\boldsymbol z) – голоморфная функция в \boldsymbol{\mathfrak{R}}\setminus(\boldsymbol\nu\cup\boldsymbol\alpha) с непрерывными предельными значениями на (\boldsymbol\nu\setminus\boldsymbol\alpha)\setminus\{\boldsymbol p_1,\boldsymbol p_2\} и (\boldsymbol\alpha\setminus\boldsymbol\nu)\setminus\{\boldsymbol p_1,\boldsymbol p_2\}, которые удовлетворяют соотношению
где мы использовали, что \lambda_{\boldsymbol\nu}(\boldsymbol s) = \lambda_{\boldsymbol\nu}(\boldsymbol s^*); на дугах, входящих в \boldsymbol\nu \cap \boldsymbol\alpha, скачки нужно складывать. При отсутствии логарифмических особенностей, т.е. когда все \alpha(e) равны нулю, утверждения леммы можно получить, суммируя \Lambda_{\boldsymbol\nu} по всем \boldsymbol\nu\in\{\boldsymbol\alpha,\boldsymbol\beta,\boldsymbol\gamma,\boldsymbol\delta\} и взяв \lambda_{\boldsymbol\nu} := \lambda_{h|\boldsymbol\nu}.
Пусть \Lambda :=\sum_{\boldsymbol\nu=\boldsymbol\alpha,\boldsymbol\beta,\boldsymbol\gamma,\boldsymbol\delta}\Lambda_{\boldsymbol\nu}. Если точка ветвления \boldsymbol p такова, что \pi(\boldsymbol p)\in\{a,a^{-1}\}, то ровно один цикл из цепи \boldsymbol\Delta проходит через \boldsymbol p (либо \boldsymbol\beta, либо \boldsymbol\delta). Более того, поскольку \boldsymbol\beta и \boldsymbol\delta разбивают \boldsymbol{\mathfrak{R}} на листы \boldsymbol{\mathfrak{R}}^{(0)}\setminus\{\boldsymbol\beta \cup \boldsymbol\delta\} и \boldsymbol{\mathfrak{R}}^{(1)}\setminus\{\boldsymbol\beta \cup \boldsymbol\delta\}, то можно использовать анализ, изложенный в [1; п. 5.2], и получить, что для k\in\{0,1\}
где ветви функции \log(p\,{-}\,{\cdot}) голоморфны в некоторой окрестности точки p, разрезанной вдоль дуги F_p. Ситуация, когда p\in\{b,b^{-1}\}, опять очень похожа на ту, которая обсуждается в [1; п. 5.2]. Ясно, что структура особенности функции в точке p определяется первым слагаемым в (3.8). Как объясняется в [12; п. I.8.5], для того, чтобы изучить эту структуру, достаточно найти функцию, которая имеет логарифмическую особенность в p и такие же скачки на циклах, составляющих \boldsymbol\Delta. Таким образом, можно проверить, что
где знак плюс выбирается, если z\in D_0, а знак минус, если z\in D_\infty, и у ветви \log({\cdot}\,{-}\,p) есть скачок на дуге F_a при p=b и на F_{a^{-1}} при p=b^{-1}.
3.4. Задача обращения Якоби
Мы определяем отображение Абеля на \boldsymbol{\mathfrak{R}} формулой
где путь интегрирования целиком лежит в \boldsymbol{\mathfrak{R}}_{\boldsymbol\alpha,\boldsymbol\beta}, и полагаем \mathfrak a(\boldsymbol z) := \mathfrak a_+(\boldsymbol z) при \boldsymbol z\in\boldsymbol\alpha \cup \boldsymbol\beta. Поскольку поверхность \boldsymbol{\mathfrak{R}} имеет род 1, любая задача обращения Якоби на \boldsymbol{\mathfrak{R}} однозначно разрешима. В частности, для заданных функции h и целого числа n\in\mathbb{N} существуют единственная точка \boldsymbol z_n=\boldsymbol z_n(h)\in\boldsymbol{\mathfrak{R}} и единственные числа j_n,m_n\in\mathbb{Z} такие, что
где \mathcal M_{\boldsymbol z_1,\boldsymbol z_2} – мероморфный дифференциал, имеющий два простых полюса в точках \boldsymbol z_1 и \boldsymbol z_2 с вычетами 1 и -1 соответственно и нулевой период по циклу \boldsymbol\alpha. Заметим, что
Из непрерывности \mathfrak a(\boldsymbol z) и единственности решения задачи обращения Якоби теперь следует, что если \boldsymbol z_n\to\infty^{(1)} по некоторой подпоследовательности индексов \mathbb{N}', то \boldsymbol z_{n-1}\to 0^{(1)} при \mathbb{N}'\ni n\to\infty, а из этого следует как бесконечность множества \mathbb{N}_\varepsilon, так и то, что при всех достаточно малых \varepsilon либо n, либо n-1 попадает в \mathbb{N}_\varepsilon. Те же рассуждения доказывают последнее утверждение леммы, поскольку сходимость \boldsymbol z_{n-1}^\star \to \infty^{(0)} по некоторой подпоследовательности индексов влечет сходимость \boldsymbol z_{n-1} \to 0^{(1)} и соответственно сходимость \boldsymbol z_n\to\infty^{(1)} по той же подпоследовательности. Лемма доказана.
3.5. Тета-функция Римана
Напомним, что тета-функция, ассоциированная с числом \mathbf B, – это целая трансцендентная функция, определяемая формулой
Функция \Theta_n мероморфна в \boldsymbol{\mathfrak{R}}_{\boldsymbol\alpha,\boldsymbol\beta}, имеет простой нуль в точке \boldsymbol z_n, простой полюс в точке \infty^{(0)}, а в остальных точках конечна и не обращается в нуль. Заметим, что она голоморфно продолжается через \boldsymbol\beta и
используя (3.3), (3.6), (3.9) и (3.13). Свойства мероморфности этой функции напрямую следуют из лемм 3.1–3.3 и 3.5, а свойство (2.21) получается комбинацией (3.4), (3.7), (3.10) и (3.14). Асимптотика (2.22) вблизи точек ветвления \boldsymbol{\mathfrak{R}} есть прямое следствие (3.11). Если теперь \Psi – функция со свойствами, описанными в формулировке теоремы, то отношение \Psi/\Psi_n – рациональная функция на \boldsymbol{\mathfrak{R}}, которая может иметь полюс только в точке \boldsymbol z_n. Поскольку риманова поверхность \boldsymbol{\mathfrak{R}} имеет род 1, на ней не существует рациональных функций с единственным полюсом. Таким образом, отношение \Psi/\Psi_n с необходимостью постоянно. Теорема и предложение доказаны.
Заметим, что образ отображения Абеля \mathfrak a(\boldsymbol z) ограничен. Напомним также, что \operatorname{Im}(\mathbf B)>0. Поэтому из (3.12) следует, что последовательность \{ n\tau+m_n \} ограничена. Следовательно, для любого \delta>0 найдется такая константа C_\delta\,{>}\,1, что
для точек \boldsymbol z вне круговых окрестностей “радиуса” \delta каждой из точек ветвления. Далее, из компактности \boldsymbol{\mathfrak{R}} и непрерывности отображения Абеля следует, что семейство функций \{ \Theta(\,\cdot\,;\boldsymbol p ) \}, где
также компактно и поэтому обязательно имеет равномерно ограниченные сверху модули в точках \boldsymbol z вдали от \infty^{(0)}. Аналогично можно увидеть, что семейство функций \{\Theta(\,\cdot\,;\boldsymbol p)/\Phi\} также имеет равномерно ограниченные сверху модули вдали от точки 0^{(1)}. Двух последних наблюдений достаточно, чтобы завершить доказательство верхних оценок в (2.27) и (2.29). Ясно, что для получения нижних оценок нужно оценить снизу модули функций \Theta(\boldsymbol z;\boldsymbol p) и \Theta(\boldsymbol z;\boldsymbol p)/\Phi вне D_\delta(\boldsymbol p). Наличие таких оценок для каждой из функций при каждом значении \boldsymbol p очевидно, а положительность их точной нижней грани (по всем \boldsymbol p) опять следует из компактности.
Далее заметим, что второй нуль функции \Upsilon_n – это точка \boldsymbol z_n^\star по лемме 3.4. Поэтому из (3.12) получаем, что
Из этого представления мы можем получить оценки (2.27) и (2.29) в точности так же, как выше. Наконец, заметим, что отношение \Psi_{n-1}^\star/\Psi_n равно \mathcal Q_{n-1}^\star/{\mathcal Q}_n в D_{\mathcal Q} и равно \mathcal R_{n-1}^\star/\mathcal R_n в D_{\mathcal R}. Значит, нам просто нужно оценить модуль |\Psi_{n-1}^\star/\Psi_n| на \boldsymbol{\mathfrak{R}}. Из (3.16) и (3.18), очевидно, следует, что
Аналогично (3.17) мы можем утверждать, что первый множитель в этом произведении равномерно ограничен сверху по всем n и на всей поверхности \boldsymbol{\mathfrak{R}}. Второй множитель представляет собой функцию, не зависящую от n, с единственным полюсом в точке \infty^{(0)}. Наконец, последнее отношение имеет единственный полюс в точке \boldsymbol z_n и поэтому равномерно ограничено сверху в \boldsymbol{\mathfrak{R}}\setminus D_\delta(\boldsymbol z_n) для всякого \delta> 0 по приведенному выше соображению компактности. Теорема 2.4 доказана.
§ 4. Доказательство теоремы 2.1 в случае \rho\in\mathcal W_2
Чтобы изучить асимптотическое поведение полиномов Q_n и линеаризованных функций остатка R_n, мы будем использовать матричный подход Римана–Гильберта, впервые предложенный А. С. Фокасом, А. Р. Итсом и А. В. Китаевым (см. [10], [11]), и нелинейный метод наискорейшего спуска, развитый П. Дейфтом и С. Чжоу (см. [9]). Для дальнейшего нам будет удобно положить
Пусть \Omega – ограниченная кольцеобразная область, содержащая F и не содержащая 0, граница которой состоит из двух гладких замкнутых жордановых кривых. Предполагая, что \partial \Omega положительно ориентирована, мы получим
\begin{equation}
0 = \int_{\partial\Omega} R_n(z)z^k\,\mathrm d z = \int_{\partial \Omega} Q_n(z)f_\rho(z)z^{k-n}\,\mathrm d z = -\int_FQ_n(s)s^{k-n}\rho(s)\,\mathrm d s
\end{equation}
\tag{4.3}
для всех k\in\{0,\ldots,n-1\}, где первое равенство следует из (4.2) и теоремы Коши, примененной к внешности \Omega, второе получается, если применить теорему Коши внутри \Omega, а последнее следует из (2.2), теоремы Фубини–Тонелли и интегральной формулы Коши. Аналогично (4.3) мы получаем
\begin{equation}
0 = \int_FQ_{n-1}^\star(s)s^{k-n}\rho(s)\,\mathrm d s, \qquad k\in\{0,\ldots,n-2\}.
\end{equation}
\tag{4.4}
Таким образом, a_n бесконечно тогда и только тогда, когда Q_{n-1}^\star также удовлетворяет (4.4) с k=n-1. Однако если это так, то Q_{n-1}^\star удовлетворяет (4.3). Обратно, если существует полином степени не выше n-1, удовлетворяющий (4.3), то он автоматически удовлетворяет (4.4), и поэтому коэффициент при z^{-n} в разложении R_{n-1}^\star в бесконечности должен быть нулевым. В целом a_n конечен, если и только если \deg(Q_n)=n, где Q_n – полином наименьшей степени, удовлетворяющий (4.3).
Тогда эта матричнозначная функция является решением следующей задачи Римана–Гильберта (\mathrm{RHP}\text{-}\boldsymbol Y). Найти (2\times2)-матричнозначную функцию \boldsymbol Y такую, что:
Действительно, можно прямой проверкой установить, что \boldsymbol Y удовлетворяет \mathrm{RHP}\text{-}\boldsymbol Y(a), если \deg(Q_n)=n, из чего также следует, что a_n конечно. Пусть Q обозначает либо Q_n, либо Q_{n-1}^\star, а R обозначает либо R_n, либо R_{n-1}^\star. Тогда из (2.12), (4.1) и формул Сохоцкого–Племеля (см. [12; п. I.4.2]) мы получаем, что
где предельные значения \rho_{|F_{-b}}(b\pm) вычисляются в соответствии с ориентацией F_{-b} (также следует помнить, что отрезок [a^{-1},a ] всегда ориентирован от a^{-1} к a). Таким образом, свойство \mathrm{RHP}\text{-}\boldsymbol Y(c) вытекает из известного поведения интегралов Коши вблизи точек разрыва веса (см. [12; пп. I.8.1–4]), а тот факт, что второй столбец не имеет логарифмической особенности в b, напрямую следует из (4.6). Утверждение, что решение задачи \mathrm{RHP}\text{-}\boldsymbol Y, если оно существует, должно иметь вид (4.5), к настоящему времени стандартно; см., например, [14; лемма 2.3] или [1; лемма 1]. Таким образом, мы доказали следующую лемму.
Лемма 4.1. Если решение задачи \mathrm{RHP}\text{-}\boldsymbol Y существует, то оно единственно и задается равенством (4.5), причем \deg(Q_n)=n. Обратно, если \deg(Q_n)=n, то (4.5) решает задачу \mathrm{RHP}\text{-}\boldsymbol Y.
4.3. Раскрытие линз
Пусть \Gamma_0 и \Gamma_\infty – две положительно ориентированные жордановы кривые, которые лежат в D_0 и D_\infty соответственно. Предположим далее, что эти кривые достаточно близки к F, так что h(z) голоморфна и не обращается в нуль в кольцеобразной области, ограниченной этими двумя кривыми. Обозначим через \Omega_0 и \Omega_\infty пересечения этой кольцеобразной области с D_0 и с D_\infty соответственно. Пусть
Лемма 4.2. Задача \mathrm{RHP}\text{-}\boldsymbol X разрешима тогда и только тогда, когда разрешима задача \mathrm{RHP}\text{-}\boldsymbol Y. Если решения задач \mathrm{RHP}\text{-}\boldsymbol X и \mathrm{RHP}\text{-}\boldsymbol Y существуют, то они единственны и связаны между собой соотношением (4.7).
Чтобы решить задачу \mathrm{RHP}\text{-}\boldsymbol N, вспомним определения (2.10) и (2.11) функций {\mathcal Q}_n и \mathcal R_n. Поскольку D_+=D_-h на \mathbb{T}, S_+S_-=D^2/h на F_a и S_+S_-=D^2h на F_{a^{-1}}, можно легко проверить, что
где соглашение о корнях из отрицательных чисел такое же, как в (2.6). Можно показать аналогично тому, как это было сделано для \varphi(z), что функция \phi(z) голоморфна вне [a^{-1},a], имеет простой нуль в точке z=0 и удовлетворяет условию \phi_-(s)\phi_+(s)\,{=}\,s при s\in[a^{-1},a]. Далее, используя (2.10) и (2.11), положим
– это голоморфная и не обращающаяся в нуль функция в {\mathbb{C}}\setminus F с нулем кратности n-1 в бесконечности. Очевидно, что функции \mathcal Q_{n-1}^\star и \mathcal R_{n-1}^\star также удовлетворяют условию (4.8). Тогда можно легко проверить, что матричнозначная функция
решает задачу \mathrm{RHP}\text{-}\boldsymbol N, где \gamma_{n-1}^\star такова, что \lim_{z\to\infty}\gamma_{n-1}^\star z^{n-1}\mathcal R_{n-1}^\star(z)\,{=}\,1. Заметим, что функция \boldsymbol M удовлетворяет условию \mathrm{RHP}\text{-}\boldsymbol Y(c) и \det(\boldsymbol N)\equiv1, потому что \det(\boldsymbol N)(z) – голоморфная функция вне точек \{a,a^{-1}\}, в которых она имеет особенности не более чем порядка квадратного корня, причем эта функция имеет значение 1 в бесконечности. Также верно, что
\begin{equation}
\det(\boldsymbol M) = (\gamma_n\gamma_{n-1}^\star)^{-1} = -\frac4{2b+a+a^{-1}}.
\end{equation}
\tag{4.13}
где c(\Gamma_0,\Gamma_\infty)\in (0,1); последнее равенство следует из (2.7) и принципа максимума модуля для голоморфных функций. Утверждение леммы теперь следует из тех же соображений, что и в [7; следствие 7.108]. Лемма доказана.
4.6. Асимптотика
Пусть \boldsymbol Z – решение задачи \mathrm{RHP}\text{-}\boldsymbol Z, даваемое леммой 4.3, и пусть \boldsymbol C,\boldsymbol M определены формулами (4.12). Тогда можно легко проверить, что функция \boldsymbol X := \boldsymbol{CZM} решает задачу \mathrm{RHP}\text{-}\boldsymbol X, и поэтому решение задачи \mathrm{RHP}\text{-}\boldsymbol Y получается из (4.7).
Для заданного замкнутого множества B\subset\overline{\mathbb{C}}\setminus F выберем \Omega_0 и \Omega_\infty так, что B\subset \overline{\mathbb{C}} \setminus (\overline\Omega_0 \cup \overline\Omega_\infty). Тогда \boldsymbol Y = \boldsymbol{CZM} на B. Значит, если ввести обозначения \upsilon_{n1}, \upsilon_{n2}, полагая, что первая строка матрицы \boldsymbol Z имеет вид \begin{pmatrix} 1+\upsilon_{n1} & \upsilon_{n2} \end{pmatrix}, мы можем записать
вследствие (4.5) и (4.12). Теперь равенства (2.13) вытекают из (2.10) и (4.10), а также из (2.11) и (4.11), поскольку мы знаем из леммы 4.3, что |\upsilon_{nk}|\leqslant c^n равномерно в \overline{\mathbb{C}} (\upsilon_{nk}(\infty)=0, так как \boldsymbol{Z}(\infty)=\boldsymbol{I}).
§ 5. Доказательство теоремы 2.3 при \rho\in\mathcal W_2
Можно проверить напрямую, что все, написанное в пп. 4.1–4.3, остается справедливым и в этом случае, за исключением условия \mathrm{RHP}\text{-}\boldsymbol Y(c), которое теперь имеет вид
Далее, формулировка задачи \mathrm{RHP}\text{-}\boldsymbol N также остается той же самой. Чтобы решить эту задачу, заметим, что условие (4.8) по-прежнему выполняется (с заменой F_{-b}^\circ на F_{-1}^\circ\cup F_1^\circ), если функции {\mathcal Q}_n и \mathcal R_n теперь определить равенствами (2.24) и (2.25). Действительно, для s\in F_a^\circ имеем
как и утверждается. Доказательство равенства (4.8) на оставшихся дугах совершенно аналогично (нужно только следить за выбранной ориентацией циклов \boldsymbol\alpha,\boldsymbol\beta,\boldsymbol\gamma,\boldsymbol\delta). Более того, функции \mathcal Q_{n-1}^\star и \mathcal R_{n-1}^\star, определенные перед формулировкой теоремы 2.3, также удовлетворяют условию (4.8), что можно показать аналогичным способом, поскольку \Psi_{n-1}^\star, очевидно, удовлетворяет (2.21). Таким образом, легко проверить, используя (4.8), что задача \mathrm{RHP}\text{-}\boldsymbol N решается для каждого такого n, что \boldsymbol z_n\neq\infty^{(1)}, с помощью формул (4.12), в которых константы \gamma_n и \gamma_{n-1}^\star опять определяются условиями
Эти константы определены корректно в силу определения множества \mathbb{N}_\varepsilon и леммы 3.4. Опять заметим, что \boldsymbol M имеет такое же поведение, как и \boldsymbol Y, вблизи E. Более того, \det(\boldsymbol N)\equiv1 по тем же причинам, что и выше, и поэтому \det(\boldsymbol M) = (\gamma_n\gamma_{n-1}^\star)^{-1}. По решению задачи \mathrm{RHP}\text{-}\boldsymbol N мы опять можем построить решение задачи \mathrm{RHP}\text{-}\boldsymbol Z. Очевидно, скачок \boldsymbol Z равен левой части (4.15). Поскольку
то из определения \mathbb{N}_\varepsilon, леммы 3.4 и соображений компактности из доказательства теоремы 2.4 следует, что последовательность \{|\gamma_n\gamma_{n-1}^\star|\}_{n\in\mathbb{N}_\varepsilon} ограничена сверху (константой, не зависящей от \varepsilon). Поэтому утверждение леммы 4.3 опять выполняется, но только для достаточно больших n\in\mathbb{N}_\varepsilon и с константой c=c_\varepsilon; для доказательства равенства в (4.15) нам нужно использовать соотношения (2.29) и (2.28), а также принцип максимума для гармонических функций. Наконец, доказательство равенств (2.26) абсолютно такое же, как и в случае теоремы 2.1. Теорема 2.3 доказана.
§ 6. Доказательство теоремы 2.1 при \rho\in\mathcal W_1
6.1. Начальная задача Римана–Гильберта
Все изложенное в п. 4.1 остается справедливым и в этом случае. Единственное изменение, которое нужно внести в рассуждения из п. 4.2, – это заменить условие \mathrm{RHP}\text{-}\boldsymbol Y(c) на условие
Здесь мы выбираем \Gamma_0 и \Gamma_\infty, как в п. 4.3, за исключением того, что теперь мы требуем, что \Gamma_0 имеет с F общую точку a, а \Gamma_\infty имеет с F общую точку a^{-1}. Мы опять определяем \Omega_0 и \Omega_\infty, как в п. 4.3, но, однако, теперь они больше не являются кольцеобразными областями. Далее, мы по-прежнему определяем \boldsymbol X формулой (4.7) с функцией \rho(s), продолженной на \Omega_0\cup\Omega_\infty, как отношение h(z)/w(z) (мы предполагаем, что множества разреза ветвей (z-a)^{\alpha+1/2} и (z-a^{-1})^{\beta+1/2} в (2.5) проходят вне пересечения некоторых окрестностей a и a^{-1} с \overline\Omega_0\,{\cup}\,\overline\Omega_\infty). Задача Римана–Гильберта \mathrm{RHP}\text{-}\boldsymbol X остается той же самой, за исключением условия \mathrm{RHP}\text{-}\boldsymbol X(c), которое заменяется в \Omega_0\cup\Omega_\infty на следующее:
при \Omega_0\cup\Omega_\infty\,{\ni}\, z\to a; и аналогичное изменение нужно сделать вблизи a^{-1}. С учетом описанных выше изменений лемма 4.2 остается справедливой.
6.3. Модельная и локальная задачи Римана–Гильберта
Модельная задача Римана–Гильберта \mathrm{RHP}\text{-}\boldsymbol N формулируется и решается в точности так же, как в случае \rho\in\mathcal W_2. Более того, опять верно, что \det(\boldsymbol N)\equiv 1 (особенности элементов \boldsymbol N в точках a,a^{-1} сокращаются при вычислении определителя). Пусть теперь U_a и U_{a^{-1}} – открытые множества, содержащие соответственно точки a и a^{-1}. Положим
где первое соотношение выполняется для s\in (\Gamma_0\cup \Gamma_\infty)\setminus ( \overline U_a \cup \overline U_{a^{-1}}), а второе для s\in \partial U_e\setminus (\Gamma_0\cup F\cup \Gamma_\infty), e\in\{a,a^{-1}\}.
Имеет место следующая лемма.
Лемма 6.1. Для всех достаточно больших n решение задачи \mathrm{RHP}\text{-}\boldsymbol Z существует и обладает следующим свойством: \boldsymbol{Z}\,{=}\,\boldsymbol{I}\,{+}\,\mathcal O(1/n) равномерно в \overline{\mathbb{C}}.
Доказательство. Доказательство того, что скачок \boldsymbol Z экспоненциально мал на (\Gamma_0\cup \Gamma_\infty)\setminus ( \overline U_a \cup \overline U_{a^{-1}}), такое же, как в случае \rho\in\mathcal W_2. Далее, имеем
на \partial U_a, и аналогичная формула верна на \partial U_{a^{-1}}. В любом случае это фиксированная матрица, не зависящая от n. Таким образом, скачок \boldsymbol Z имеет порядок \boldsymbol I + \mathcal O(1/n) на \partial U_a \cup \partial U_{a^{-1}}. Утверждение леммы теперь можно получить так же, как в случае \rho\in\mathcal W_2.
6.5. Асимптотика
Формулы (2.13) теперь получаются в точности так же, как в случае \rho\in\mathcal W_2.
6.6. Решение задачи \mathrm{RHP}\text{-}{\boldsymbol P}_e, e\in\{a,a^{-1}\}
Мы построим только матрицу \boldsymbol P_a, поскольку построение \boldsymbol P_{a^{-1}} полностью аналогично.
6.6.1. Модельная задача
Ниже мы всегда будем предполагать, что вещественная прямая и ее интервалы ориентированы слева направо. Положим
и будем считать, что лучи I_\pm ориентированы в направлении к началу координат. Для данного \alpha>-1 пусть \boldsymbol\Psi_\alpha – матричнозначная функция такая, что:
где \displaystyle \boldsymbol S(\zeta) := \frac{\zeta^{-\sigma_3/4}}{\sqrt2}\begin{pmatrix} 1 & \mathrm i \\ \mathrm i & 1 \end{pmatrix}, и мы выбираем главную ветвь функции \zeta^{1/4}.
Явное построение такой матричнозначной функции можно найти в [14] (оно использует модифицированные функции Бесселя и Ганкеля). Заметим, что
где функция \varphi задана в (2.6). Из (2.7) следует, что функция \zeta_a голоморфно продолжается через F_a. Также из явного представления для \varphi следует, что \zeta_a обращается в нуль в точке a. Более того, поскольку
этот нуль в точке a обязательно простой. Заметим также, что |\varphi_+|=|\varphi_-| на [a^{-1},a], и поэтому |s/\varphi_\pm^2(s)| \equiv 1 на этом отрезке в соответствии с (2.7). Таким образом, функция \zeta_a отображает F_a во множество отрицательных вещественных чисел. Также несложно проверить, что оставшиеся вещественные числа в U_a отображаются функцией \zeta_a в положительные вещественные числа. Положим
Из предыдущего обсуждения понятно, что \zeta_a(U_a^\pm) \subset \{\pm\operatorname{Im}(z)>0\}. Пусть \Gamma_0^\pm : = \Gamma_0\cap U_a^\pm. Заметим, что в соответствии с выбранной ориентацией кривой \Gamma_0 дуга \Gamma_0^+ ориентирована по направлению к a, а дуга \Gamma_0^- ориентирована по направлению от a. Поскольку мы имели некоторую свободу в выборе кривой \Gamma_0, выберем ее так, что \zeta_a(\Gamma_0^\pm) \subset I_\pm.
Наконец, в дальнейшем мы будем понимать под \zeta_a^{1/2} ветвь, задаваемую выражением в скобках в (6.6) со скачком вдоль F_a. В частности, тогда
где матрица \boldsymbol D была определена в (6.3). Аналогично будем обозначать через \zeta_a^{1/4} ветвь, которая отображает U_a в сектор |\arg(z)|<\pi/4. Например, на F_a выполняется равенство \zeta_{a+}^{1/4} = \mathrm i \zeta_{a-}^{1/4}.
6.6.3. Матрица \boldsymbol P_a
При условиях, наложенных на функции класса \mathcal W_1, верно, что
где функция h_* голоморфна и не обращается в нуль в U_a, \alpha>-1, а для степеней, содержащих \alpha, выбираются главные значения. Напомним также, что \rho на F_a определяется как предельное значение на F_a функции \rho_{|U_a^+}. Приведенную выше формулу для \rho можно эквивалентным образом переписать в виде
где опять выбираются главные ветви степенных функций. Тогда r_a – голоморфная и не обращающаяся в нуль функция в U_a\setminus F_a, для которой выполняются соотношения
удовлетворяет \mathrm{RHP}\text{-}{\boldsymbol P}_e(a, b, c), если \boldsymbol E_a – произвольная голоморфная матричнозначная функция (отметим, что ориентация дуги \zeta_a(\Gamma_0^-) противоположна ориентации луча I_-). Далее, из \mathrm{RHP}\text{-}\boldsymbol N(b), (2.7) и (6.5) вытекает, что функция
функция \boldsymbol E_a действительно голоморфна в U_a. Наконец, условие \mathrm{RHP}\text{-}{\boldsymbol P}_e(d) следует из (6.7) и \mathrm{RHP}\text{-}{\boldsymbol \Psi_\alpha}(d).
§ 7. Доказательство теоремы 2.3 при \rho\in\mathcal W_1
Как обычно, результаты пп. 4.1, 4.2 буквально переносятся на данный случай, если заменить условие \mathrm{RHP}\text{-}\boldsymbol Y(c) на (6.1).
7.1. Раскрытие линз
Мы выбираем кривые \Gamma_0 и \Gamma_\infty так же, как в п. 4.3, однако требуем, чтобы \Gamma_0 имела с F общую точку a, а \Gamma_\infty имела с F общую точку a^{-1} (рис. 4).
Более того, мы также выбираем открытые ориентированные дуги \Gamma_{v,1}, \Gamma_{v,-1}, \Gamma_v, соединяющие v с \Gamma_0\cup\Gamma_\infty, v\in\{b,b^{-1}\}, как показано на рис. 4. Кроме внутренней области \Gamma_0 и внешней области \Gamma_\infty, объединение выбранных дуг, которое мы назовем \Gamma, вместе с F ограничивает восемь областей, которые мы помечаем, как на рис. 4. Заметим, что \rho имеет голоморфные и не обращающиеся в нуль продолжения в каждую из этих восьми областей (если нужно, мы можем для этого выбрать дуги \Gamma_0 и \Gamma_\infty ближе к F). Мы предполагаем, что все введенные дуги гладкие. Положим
где \Omega_+ := \Omega_{\infty,1}\cup\Omega_{\infty,-1}\cup\Omega_{b,-1}\cup\Omega_{b^{-1},1} и \Omega_- := \Omega_{0,1}\cup\Omega_{0,-1}\cup\Omega_{b,1}\cup\Omega_{b^{-1},-1}.
Тогда задачу Римана–Гильберта для \boldsymbol X можно поставить следующим образом:
при s\in(\Gamma_0\cup\Gamma_\infty)\,{\setminus}\,\{a,a^{-1}\}, где знак минус следует выбирать для части \Gamma_0, ограничивающей \Omega_{b,-1}, и для части \Gamma_\infty, ограничивающей \Omega_{b^{-1},-1}, и
При описанных выше изменениях лемма 4.2 выполняется с заменой (4.7) на (7.1).
7.2. Модельная и локальная задачи Римана–Гильберта
Модельная задача Римана–Гильберта \mathrm{RHP}\text{-}\boldsymbol N формулируется и решается в точности так же, как в случае \rho\in\mathcal W_2 для n\in \mathbb{N}_\varepsilon.
Пусть теперь U_e – открытое множество, содержащее e\in E. В случае U_{b^{\pm1}} мы будем предполагать, что эти множества не пересекают множество \Gamma_0\,{\cup}\,\Gamma_\infty и полностью содержат \Gamma_{b^{\pm1}},\Gamma_{b^{\pm1},1} и \Gamma_{b^{\pm1},-1} (рис. 5). Положим \boldsymbol D(z) := \Phi(\boldsymbol z)^{n\sigma_3}, \boldsymbol z\in D_{\mathcal Q}, где \Phi задается формулой (3.3), а открытое множество D_{\mathcal Q} было определено перед формулировкой теоремы 2.4. Как и в п. 6.3, нам понадобится решить задачу \mathrm{RHP}\text{-}{\boldsymbol P}_e для всех e\in E. Как и в п. 6.3, мы отложим построение матричнозначных функций, решающих эту задачу, до конца этого параграфа.
Лемма 7.1. Для достаточно больших n\in\mathbb{N}_\varepsilon решение задачи \mathrm{RHP}\text{-}\boldsymbol Z существует и удовлетворяет условию \boldsymbol{Z}=\boldsymbol{I}\,{+}\, \mathcal O( 1/n), где \mathcal O(\cdot) равномерно на \overline{\mathbb{C}} и зависит от \varepsilon.
Доказательство. Доказательство того, что скачок \boldsymbol Z экспоненциально мал на (\Gamma_0\cup \Gamma_\infty)\setminus\bigcup_{e\in E} \overline U_e, такое же, как в случае \rho\in\mathcal W_2. Далее,
на \partial U_e. Из (2.24), (2.25), (3.16) и равенства \Phi(\boldsymbol z)\Phi(\boldsymbol z^*)=z (см. лемму 3.1) следует, что первая строка матрицы \boldsymbol{MD}^{-1} равна
В процессе доказательства теоремы 2.4 было показано, что эти функции имеют модули, равномерно ограниченные сверху на компактных подмножествах {\mathbb{C}}. Подобным же образом можно показать, что то же самое верно и для второй строки матрицы \boldsymbol{MD}^{-1}. Поскольку
\begin{equation*}
\det(\boldsymbol{MD}^{-1}) = \det(\boldsymbol M) = (\gamma_n\gamma_{n-1}^\star)^{-1}
\end{equation*}
\notag
и константы |\gamma_n\gamma_{n-1}^\star| равномерно ограничены сверху для n\in\mathbb{N}_\varepsilon, мы получаем, что скачок \boldsymbol Z имеет вид \boldsymbol I + \mathcal O(1/n) на \partial U_e для каждого e\in E. Теперь утверждение леммы выводится так же, как в случае \rho\in\mathcal W_2. Лемма доказана.
7.4. Асимптотика
Формулы (2.26) теперь выводятся так же, как в случае \rho\in\mathcal W_2.
7.5. Решение задачи \mathrm{RHP}\text{-}{\boldsymbol P}_e для e\in\{a,a^{-1}\}
Как и в п. 6.6, мы построим только матрицу \boldsymbol P_a. Построение опять основано на матричнозначной функции \boldsymbol\Psi_\alpha, решающей задачу \mathrm{RHP}\text{-}{\boldsymbol \Psi_\alpha}. Мы опять начинаем с того, что определяем специальное конформное отображение в окрестности точки a.
7.5.1. Конформное отображение
Используя те же обозначения, что и в п. 3.1, положим
Поскольку v_+=-v_- на F_a^\circ, функция \zeta_a голоморфна в U_a. Более того, поскольку v имеет особенность типа квадратного корня в точке a, функция \zeta_a имеет в точке a простой нуль. Таким образом, можно выбрать U_a столь малым, чтобы \zeta_a была конформна в \overline U_a. Напомним, что дифференциал v(s)\,\mathrm d s/s чисто мнимый на F_a^\circ (см. последнюю часть доказательства леммы 3.1). Поэтому \zeta_a отображает точки F_a в отрицательные вещественные числа. Поскольку мы имели некоторую свободу в выборе кривой \Gamma_0, теперь мы выберем эту кривую в U_a так, чтобы часть \Gamma_0, ограничивающая \Omega_{b,1} (скажем, \Gamma_0^+), отображалась в I_+, а часть, ограничивающая \Omega_{b,-1} (скажем, \Gamma_0^-), отображалась в I_-. Отметим, что ориентация дуги \zeta_a(\Gamma_0^-) противоположна ориентации луча I_-.
В дальнейшем под \zeta_a^{1/2} мы будем понимать ветвь, задаваемую выражением в скобках в (7.2). Уравнения (3.1) и (3.3) показывают, что
Так как \Phi(\boldsymbol z)\Phi(\boldsymbol z^*)\equiv z, соотношение (6.7) остается справедливым для \zeta_a, как и выше, и \boldsymbol D(z)=\Phi(\boldsymbol z)^{n\sigma_3}, \boldsymbol z\in D_\mathcal{Q}. Наконец, как и в случае вещественного a, на дуге F_a выполняется равенство \zeta_{a+}^{1/4} = \mathrm i \zeta_{a-}^{1/4}.
7.5.2. Матрица \boldsymbol P_a
Пусть J_a – такая дуга в U_a, выходящая из a, что \zeta_a(J_a)\,{\subset} [0,\infty). В силу условий, наложенных на функции класса \mathcal W_1, выполняется равенство
где \rho_* – функция, голоморфная и не обращающаяся в нуль в U_a, а (z-a)^\alpha обозначает ветвь, голоморфную в U_a\setminus J_a. Пусть U_a^\pm – связные компоненты множества U_a\setminus(F_a\cup J_a), содержащие \Gamma_0^\pm. Определим функцию
равномерно в {\mathbb{C}}\setminus(I_+\cup I_-\cup(-\infty,\infty)), где функция \boldsymbol S(\zeta) была определена в \mathrm{RHP}\text{-}{\boldsymbol \Psi_\alpha}(d).
Такая матричнозначная функция была построена в [8] с помощью функций Эйри.
7.6.2. Конформное отображение
Используя те же обозначения, что и в п. 3.1, положим
Поскольку v_+=-v_- на F_{a^{-1}}^\circ, функция \zeta_{b^{-1}}^3 голоморфна в U_{b^{-1}}. Более того, поскольку v стремится к нулю при z\to b^{-1} как квадратный корень, функция \zeta_{b^{-1}}^3 имеет нуль третьего порядка в b^{-1}, и поэтому \zeta_{b^{-1}} голоморфна в U_{b^{-1}}. Размер окрестности U_{b^{-1}} можно подобрать так, чтобы функция \zeta_{b^{-1}} была конформна в U_{b^{-1}}. Напомним, что дифференциал v(s)\,\mathrm d s/s на F чисто мнимый. Таким образом, мы можем выбрать такую ветвь выражения (7.3), что
Далее мы понимаем под \zeta_{b^{-1}}^{3/2} ветвь, задаваемую выражением в скобках в (7.3), и выбираем ветвь \zeta_{b^{-1}}^{1/4} со скачком вдоль F_{a^{-1}}, для которой \zeta_{b^{-1}+}^{1/4}=\mathrm i \zeta_{b^{-1}-}^{1/4}.
Пусть точка z\in U_{b^{-1}}\setminus (F_{a^{-1}}\cup F_1) принадлежит компоненте, содержащей F_{-1}, скажем, U_{b^{-1}}^1 (см. рис. 4 и 5). Пусть \gamma – путь из a в b^{-1}, а \gamma_z – путь из b^{-1} в z, которые полностью лежат в U_{b^{-1}}^1. Как обычно, обозначим через B^{(i)} поднятие множества B на \boldsymbol{\mathfrak{R}}^{(i)}. Тогда пути
идут соответственно из \boldsymbol a в z^{(0)} и в z^{(1)} и лежат в \boldsymbol{\mathfrak{R}}_{\boldsymbol\alpha,\boldsymbol\beta} (формально циклы \boldsymbol\alpha, \boldsymbol\beta в (7.4) нужно деформировать в гомологичные им циклы в \boldsymbol{\mathfrak{R}}_{\boldsymbol\alpha,\boldsymbol\beta}). Тогда, используя пути (7.4) в формуле (3.3) и вспоминая определение (3.2), мы получим
Пусть теперь точка z \in U_{b^{-1}}\setminus (F_{a^{-1}}\cup F_1) принадлежит компоненте, не содержащей F_{-1}, скажем U_{b^{-1}}^2. Выберем путь \gamma_z, лежащий в этой компоненте. Тогда пути
соединяют \boldsymbol a с z^{(0)} и z^{(1)} соответственно и лежат в \boldsymbol{\mathfrak{R}}_{\boldsymbol\alpha,\boldsymbol\beta} (с той же оговоркой, что и ранее). Таким образом, из (7.5), (3.3) и (3.2) следует, что
удовлетворяет условиям \mathrm{RHP}\text{-}{\boldsymbol P}_e(a, b, c) для e=b^{-1}, если \boldsymbol E_{b^{-1}} – произвольная голоморфная матричнозначная функция в U_{b^{-1}}. Таким образом, остается только выбрать матрицу \boldsymbol E_{b^{-1}} так, чтобы выполнялось условие \mathrm{RHP}\text{-}{\boldsymbol P}_e(d). Положим
Напомним, что \boldsymbol D(z)=\Phi^{n\sigma_3}(\boldsymbol z) при \boldsymbol z\in D_{\mathcal Q}. Обозначим элемент матрицы \boldsymbol D(z) с индексом (1,1) через d(z). Тогда
в силу свойств \Phi(\boldsymbol z)\Phi(\boldsymbol z^*) = z и (3.4); здесь предельные значения d(z) берутся на дугах в комплексной плоскости, а предельные значения \Phi(\boldsymbol z) берутся на циклах на поверхности \boldsymbol{\mathfrak{R}}. Используя \mathrm{RHP}\text{-}\boldsymbol N(b), (6.5), (7.7) и явные определения для r_{b^{-1}}, \boldsymbol J и \boldsymbol K, можно прямой проверкой выяснить, что функция \boldsymbol E_{b^{-1}} голоморфна в U_{b^{-1}}\setminus \{b^{-1}\}. Кроме того, из (2.22) и поведения функции \boldsymbol S в нуле следует, что
Это показывает, что фактически эта функция голоморфна на всем открытом множестве U_{b^{-1}}. Соотношение \mathrm{RHP}\text{-}{\boldsymbol P}_e(d) следует из \mathrm{RHP}\text{-}{\boldsymbol \Psi}(d) и (7.6).
Список литературы
1.
A. I. Aptekarev, M. L. Yattselev, “Padé approximants for functions with branch points – strong asymptotics of Nuttall–Stahl polynomials”, Acta Math., 215:2 (2015), 217–280
2.
L. Baratchart, H. Stahl, M. Yattselev, “Weighted extremal domains and best rational approximation”, Adv. Math., 229:1 (2012), 357–407
3.
L. Baratchart, M. Yattselev, “Convergent interpolation to Cauchy integrals over analytic arcs”, Found. Comput. Math., 9:6 (2009), 675–715
4.
L. Baratchart, M. Yattselev, “Convergent interpolation to Cauchy integrals over analytic arcs with Jacobi-type weights”, Int. Math. Res. Not. IMRN, 2010:22 (2010), 4211–4275
5.
В. И. Буслаев, “О сходимости многоточечных аппроксимаций Паде кусочно аналитических функций”, Матем. сб., 204:2 (2013), 39–72; англ. пер.: V. I. Buslaev, “Convergence of multipoint Padé approximants of piecewise analytic functions”, Sb. Math., 204:2 (2013), 190–222
6.
В. И. Буслаев, “О сходимости m-точечных аппроксимаций Паде набора многозначных аналитических функций”, Матем. сб., 206:2 (2015), 5–30; англ. пер.: V. I. Buslaev, “Convergence of m-point Padé approximants of a tuple of multivalued analytic functions”, Sb. Math., 206:2 (2015), 175–200
7.
P. Deift, Orthogonal polynomials and random matrices: a Riemann–Hilbert approach, Courant Lect. Notes Math., 3, Courant Inst. Math. Sci., New York; Amer. Math. Soc., Providence, RI, 2000, viii+273 pp.
8.
P. Deift, T. Kriecherbauer, K. T.-R. McLaughlin, S. Venakides, X. Zhou, “Strong asymptotics of orthogonal polynomials with respect to exponential weights”, Comm. Pure Appl. Math., 52:12 (1999), 1491–1552
9.
P. Deift, X. Zhou, “A steepest descent method for oscillatory Riemann–Hilbert problems. Asymptotics for the MKdV equation”, Ann. of Math. (2), 137:2 (1993), 295–368
10.
A. S. Fokas, A. R. Its, A. V. Kitaev, “Discrete Panlevé equations and their appearance in quantum gravity”, Comm. Math. Phys., 142:2 (1991), 313–344
11.
A. S. Fokas, A. R. Its, A. V. Kitaev, “The isomonodromy approach to matrix models in 2D quantum gravitaty”, Comm. Math. Phys., 147:2 (1992), 395–430
12.
Ф. Д. Гахов, Краевые задачи, 3-е изд., Наука, М., 1977, 640 с. ; англ. пер. 2-го изд.: F. D. Gakhov, Boundary value problems, Pergamon Press, Oxford–New York–Paris; Addison-Wesley Publishing Co., Inc., Reading, MA–London, 1966, xix+561 с.
13.
А. А. Гончар, Е. А. Рахманов, “Равновесные распределения и скорость рациональной аппроксимации аналитических функций”, Матем. сб., 134(176):3(11) (1987), 306–352; англ. пер.: A. A. Gonchar, E. A. Rakhmanov, “Equilibrium distributions and degree of rational approximation of analytic functions”, Math. USSR-Sb., 62:2 (1989), 305–348
14.
A. B. Kuijlaars, K. T.-R. McLaughlin, W. Van Assche, M. Vanlessen, “The Riemann–Hilbert approach to strong asymptotics for orthogonal polynomials on [-1,1]”, Adv. Math., 188:2 (2004), 337–398
15.
G. López, “Szegő's theorem for polynomials orthogonal with respect to varying measures”, Orthogonal polynomials and their applications (Segovia, 1986), Lecture Notes in Math., 1329, Spinger, Berlin, 1988, 255–260
16.
E. A. Rakhmanov, “Orthogonal polynomials and S-curves”, Recent advances in orthogonal polynomials, special functions, and their applications, Contemp. Math., 578, Amer. Math. Soc., Providence, RI, 2012, 195–239
17.
T. Ransford, Potential theory in the complex plane, London Math. Soc. Stud. Texts, 28, Cambridge Univ. Press, Cambridge, 1995, x+232 pp.
18.
H. Stahl, “Extremal domains associated with an analytic function. I”, Complex Variables Theory Appl., 4:4 (1985), 311–324
19.
H. Stahl, “Extremal domains associated with an analytic function. II”, Complex Variables Theory Appl., 4:4 (1985), 325–338
20.
H. Stahl, “The structure of extremal domains associated with an analytic function”, Complex Variables Theory Appl., 4:4 (1985), 339–354
21.
H. Stahl, “Orthogonal polynomials with complex-valued weight function. I”, Constr. Approx., 2:3 (1986), 225–240
22.
H. Stahl, “Orthogonal polynomials with complex-valued weight function. II”, Constr. Approx., 2:3 (1986), 241–251
23.
H. Stahl, “The convergence of Padé approximants to functions with branch points”, J. Approx. Theory, 91:2 (1997), 139–204
24.
H. Stahl, “Strong asymptotics for orthonormal polynomials with varying weights”, Acta Sci. Math. (Szeged), 66:1-2 (2000), 147–192
25.
M. L. Yattselev, “Symmetric contours and convergent interpolation”, J. Approx. Theory, 225 (2018), 76–105
26.
Э. И. Зверович, “Краевые задачи теории аналитических функций в гёльдеровских классах на римановых поверхностях”, УМН, 26:1(157) (1971), 113–179; англ. пер.: È. I. Zverovich, “Boundary value problems in the theory of analytic functions in Hölder classes on Riemann surfaces”, Russian Math. Surveys, 26:1 (1971), 117–192
Образец цитирования:
М. Л. Ятцелев, “Сходимость двухточечных аппроксимаций Паде к кусочно голоморфным функциям”, Матем. сб., 212:11 (2021), 128–164; M. L. Yattselev, “Convergence of two-point Padé approximants to piecewise holomorphic functions”, Sb. Math., 212:11 (2021), 1626–1659
С. П. Суетин, “Прямое доказательство теоремы Шталя для некоторого класса алгебраических функций”, Матем. сб., 213:11 (2022), 102–117; S. P. Suetin, “A direct proof of Stahl's theorem for a generic class of algebraic functions”, Sb. Math., 213:11 (2022), 1582–1596