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

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

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



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






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


Теоретическая и математическая физика, 2024, том 219, номер 3, страницы 440–461
DOI: https://doi.org/10.4213/tmf10649
(Mi tmf10649)
 

Определитель Ганкеля для полуклассического унитарного ансамбля Лагерра, уравнения Пенлеве IV и Гойна

Дань Ван

School of Computer Science and Artificial Intelligence, Changzhou University, China
Список литературы:
Аннотация: Проведен анализ асимптотического поведения определителя Ганкеля, порожденного полуклассическим весом Лагерра. Для этого применяются лестничные операторы и исследуется эволюция параметров. Установлено, что вспомогательная величина, связанная с полуклассическим весом Лагерра, удовлетворяет уравнению Пенлеве IV, в котором выполнено соответствующее преобразование переменных. С помощью метода кулоновской жидкости получено разложение при больших $n$ логарифма определителя Ганкеля. Это позволяет разобраться в скейлинге и флуктуациях определителя, а также способствует более глубокому пониманию его поведения в полуклассическом ансамбле Лагерра. Исследуется асимптотическая эволюция ортогональных многочленов по полуклассическому весу Лагерра, рассматривается частный случай. Изучаются свойства и характеристики этих полиномов, связанные с ансамблем. Объясняется связь между дифференциальными уравнениями второго порядка, которым удовлетворяют ортогональные многочлены по полуклассическому весу Лагерра, и триконфлюэнтными уравнениями Гойна или биконфлюэнтными уравнениями Гойна.
Ключевые слова: Ганкеля, асимптотика, уравнение Пенлеве IV, уравнение Гойна.
Финансовая поддержка Номер гранта
Changzhou University ZMF22020116
Д. Ван выражает благодарность университету Чанчжоу за грант № ZMF22020116.
Поступило в редакцию: 30.11.2023
После доработки: 15.01.2024
Дата публикации: 30.05.2024
Англоязычная версия:
Theoretical and Mathematical Physics, 2024, Volume 219, Issue 3, Pages 913–932
DOI: https://doi.org/10.1134/S0040577924060035
Реферативные базы данных:
Тип публикации: Статья
MSC: 15B52, 42C05

1. Введение

В заданном пространстве эрмитовых $(n \times n)$-матриц $\mathcal{H}_n$ [1] введем плотность вероятности для $M=(M_{ij}) \in \mathcal{H}_n$:

$$ \begin{equation*} \begin{gathered} \, \textsf{P}(M)\,dM \propto e^{-\operatorname{Tr}(\mathrm{v}(M))}\operatorname{vol}(dM), \\ \operatorname{vol}(dM)=\prod_{i=1}^{n}dM_{ii}\prod_{1\leqslant j<k \leqslant n}d(\operatorname{Re} M_{jk})\,d(\operatorname{Im} M_{jk}), \end{gathered} \end{equation*} \notag $$
где $\operatorname{vol}(dM)$ – элемент объема [2], а $\mathrm{v}(M)$ – матричная функция [3], заданная жордановой канонической формой. При разложении по собственным значениям и собственным векторам унитарный ансамбль индуцирует совместную функцию плотности вероятности на наборе собственных значений $\{z_k\}_{k=1}^{n} \in \mathbb{R}^n$ по формуле
$$ \begin{equation*} \frac{1}{D_n[w]}\frac{1}{n!}\prod_{1\leqslant j<i\leqslant n}(z_i-z_j)^2\prod_{k=1}^{n}w(z_k)\,dz_k, \end{equation*} \notag $$
где $w(z)$ – весовая функция с носителем на $[a,b]$ и конечными моментами
$$ \begin{equation*} \mu_j:=\int_a^b z^{j}w(z)\,dz,\qquad j \in \{0,1,2,\ldots\}. \end{equation*} \notag $$
Нормировочную константу $D_n[w]$ в совместной функции плотности вероятности можно вычислить точно – как определитель матрицы Ганкеля на основе тождества Андреева или Гейне [4]:
$$ \begin{equation*} D_n[w]:=\frac{1}{n!}\int_{[a,b]^n }\prod_{1\leqslant j<i\leqslant n}(z_i-z_j)^2\prod_{k=1}^{n}w(z_k)\,dz_k =\operatorname{det}(\mu_{i+j})|_{i,j=0}^{n-1}. \end{equation*} \notag $$

Для обсуждаемой в данной работе задачи формула для вероятности $\mathbb{P}_n(s,t)$ того, что собственные числа унитарного ансамбля больше $s$, основывается на разложении по собственным числам и векторам,

$$ \begin{equation*} \begin{aligned} \, \mathbb{P}_n(s,t)={}&\frac{1}{n!}\int_{(s,\infty)^n}\prod_{1\leqslant l<k\leqslant n}(z_k-z_l)^2\prod_{j=1}^nw(z_j,t)\,dz_j \times{} \\ &\qquad\qquad\qquad\qquad\times \biggl\{\frac{1}{n!} \int_{(0,\infty)^n}\prod_{1\leqslant l<k\leqslant n}(z_k-z_l)^2\prod_{j=1}^nw(z_j,t)\,dz_j\biggr\}^{-1}={} \\ ={}&\operatorname{det}\biggl(\int_s^{\infty}z^{l+k}w(z,t)\,dz\biggr)\bigg|_{l,k=0}^{n-1}\bigg/\operatorname{det}\biggl(\int_0^{\infty}z^{l+k}w(z,t)\,dz\biggr)\bigg|_{l,k=0}^{n-1}, \end{aligned} \end{equation*} \notag $$
где $w(z,t)$ – положительная весовая функция с носителем $[0,\infty)$, которая показывает, что определитель Ганкеля связан не только с распределением собственных чисел унитарных ансамблей, но и с вероятностью ошибки беспроводных коммуникационных систем с многоканальной принимающей и передающей антенной. Среди прочих применений отметим применения, указанные в работах [2], [5]–[8]. Фундаментальный подход к изучению определителя Ганкеля включает представление его в виде произведения квадратов $L^2$-норм нестандартных ортогональных полиномов и использование условий совместности с лестничными операторами для ортогональных полиномов по данному весу. С помощью этого метода можно вывести ряд разностных и дифференциальных уравнений, характеризующих определитель Ганкеля. Отметим результаты, полученные в работах [3], [9]–[11] для смещенного или классического веса Якоби, в работах [12]–[14] для классического или сингулярно возмущенного веса Лагерра, в работах [15], [16] для деформации классического гауссова веса.

В настоящей работе мы также используем указанный выше подход при изучении определителя Ганкеля

$$ \begin{equation*} D_n(s,t)=\frac{1}{n!}\int_{[s,\infty)^n}\prod_{1\leqslant j<k\leqslant n}(x_j-x_k)^2\prod_{l=1}^nw(x_l,t)\,dx_l=\operatorname{det}\biggl(\int_s^{\infty}z^{i+j}w(z,t)\,dz\biggr)\bigg|_{i,j=0}^{n-1}, \end{equation*} \notag $$
где момент
$$ \begin{equation*} \mu_k(s,t)=\int_s^{\infty}z^kw(z,t)\,dz \end{equation*} \notag $$
связан с полуклассическим весом Лагерра с двумя параметрами $s$ и $t$, которые называются соответственно “щель” и “время”,
$$ \begin{equation} w(z,t)=e^{-z^2+tz},\qquad z\in[s,\infty),\,\, t>0,\,\, s>0. \end{equation} \tag{1.1} $$

Замечание 1. Чтобы определить вероятность того, что все собственные значения унитарного ансамбля превосходят некоторое пороговое значение $s$, рассмотрим связанный с ансамблем вес $w(z,t)$. Эту вероятность можно представить как отношение $D_n(s,t)$ к его значению при $s=0$, где $D_n(s,t)$ является определителем $n$-мерных матриц Ганкеля, порожденных моментами веса $w(z,t)$ для $z$ в интервале $[s,\infty)$. Исследование, проведенное в работе [14], было сосредоточено исключительно на определителе Ганкеля $D_n(s,t)$, связанном с сингулярно возмущенным унитарным ансамблем Лагерра, который характеризуется весом

$$ \begin{equation} w(z,t)=z^{\alpha}e^{-z-t/z},\qquad z\in[s,\infty)\subset(0,\infty),\,\, t\geqslant 0,\,\, \alpha>-1. \end{equation} \tag{1.2} $$
Это частное исследование способствовало более глубокому пониманию сути проблемы и является одной из мотиваций настоящей работы.

Другая мотивация, побуждающая рассмотреть задачу полуклассического унитарного ансамбля Лагерра (1.1), связана с работой [17], в которой объяснена связь между рекуррентными коэффициентами ортогональных полиномов по полуклассическому весу Лагерра

$$ \begin{equation} w(z,t)=z^{\alpha} e^{-z^2+tz},\qquad z\in[0,\infty),\,\, t\in\mathbb{R},\,\, \alpha>-1, \end{equation} \tag{1.3} $$
и классическими решениями уравнения Пенлеве IV. Эта связь между ортогональными полиномами и уравнениями Пенлеве является интересной областью математической физики, ее изучение внесло большой вклад в исследования и ортогональных полиномов, и уравнений Пенлеве. На основе указанных мотиваций мы сосредотачиваемся на определителе Ганкеля для полуклассического веса Лагерра (1.1) с параметрами “щель” и “время”.

Ниже исследуется определитель Ганкеля в скейлинговых пределах с привлечением двух параметров и порядка матрицы Ганкеля. Эти исследования вносят значительный вклад в понимание поведения и свойств определителя Ганкеля.

Пусть $\{P_n(z;s,t)\}$ – многочлены степени $n$ по весу (1.1) на интервале $[s,\infty)$,

$$ \begin{equation} \int_s^{\infty}P_j(z;s,t)P_k(z;s,t)w(z,t)\,dz=h_j(s,t)\delta_{j,k},\qquad j,k=0,1,2,\dots, \end{equation} \tag{1.4} $$
при этом
$$ \begin{equation} P_n(z;s,t)=z^n+p(n,s,t)z^{n-1}+\cdots+P_n(0;s,t), \end{equation} \tag{1.5} $$
где $h_j(s,t)$ – квадрат $L^2$-нормы многочлена $P_n(z;s,t)$, а $\delta_{j,k}$ – символ Кронекера. Следовательно, $D_n(s,t)$ можно записать в другом представлении:
$$ \begin{equation} D_n(s,t)=\prod_{j=0}^{n-1}h_j(s,t). \end{equation} \tag{1.6} $$

Отсюда следует трехчленное рекуррентное соотношение

$$ \begin{equation} zP_n(z;s,t)=P_{n+1}(z;s,t)+\alpha_n(s,t)P_n(z;s,t)+\beta_n(s,t)P_{n-1}(z;s,t), \end{equation} \tag{1.7} $$
$n=0,1,2,\dots$, с начальными условиями
$$ \begin{equation*} \beta_0(s,t)P_{-1}(z;s,t)=0,\qquad P_0(z;s,t)=1. \end{equation*} \notag $$
Простым следствием формул (1.4), (1.5) и (1.7) являются соотношения
$$ \begin{equation} \alpha_n(s,t) =p(n,s,t)-p(n+1,s,t), \end{equation} \tag{1.8} $$
$$ \begin{equation} \beta_n(s,t) =\frac{h_n(s,t)}{h_{n-1}(s,t)}, \end{equation} \tag{1.9} $$
причем телескопическая сумма выражений (1.8) имеет вид
$$ \begin{equation} \sum_{j=0}^{n-1}\alpha_j(s,t)=-p(n,s,t), \end{equation} \tag{1.10} $$
где $p(0,s,t)=0$.

Статья имеет следующую структуру. В разделе 2 применяется метод лестничных операторов вместе с соответствующими дополнительными условиями для введения двух вспомогательных величин, а именно $R_n(s,t)$ и $r_n(s,t)$, которые близко связаны с весом (1.1). Установлена разностная связь между $R_n(s,t)$, $r_n(s,t)$ и рекуррентными коэффициентами в (1.7). В разделе 3 уравнение (1.4) продифференцировано по $s$ и $t$ для получения системы уравнений в частных производных, в которые входят $\alpha_n(s,t)$, $\beta_n(s,t)$, $R_n(s,t)$ и $r_n(s,t)$. Показано, что уравнения в частных производных второго порядка, управляющие вспомогательной величиной $R_n(s,t)$, при подходящей замене переменных связаны с уравнением Пенлеве IV. Кроме того, рекуррентные коэффициенты выражаются через самих себя с помощью формул, включающих частные производные. Раздел 4 посвящен применению метода жидкости Дайсона–Кулона к анализу асимптотического поведения логарифмической формы определителя Ганкеля при $n \to \infty$. Для оценки асимптотического разложения при больших $n$ применяется система нелинейных дифференциальных уравнений, в которую входят рекуррентные коэффициенты. В разделе 5 численно представлены две асимптотические характеристики, полученные на основе результатов раздела 4. Мы аппроксимируем выражение для $P_n(0;s,t)$ и показываем, что дифференциальные уравнения второго порядка, которым удовлетворяют ортогональные полиномы по весу (1.1), приблизительно эквивалентны триконфлюэнтному уравнению Гойна в скейлинговом пределе при $n \to \infty$. В разделе 6 представлены выводы.

2. Предварительные сведения и разностные уравнения

Соотношения для лестничных операторов (см., например, [3], [6], [7], [11], [18], [19]), которым удовлетворяют ортогональные многочлены $\{P_n(z;s,t)\}$ по весу (1.1) на интервале $[s,\infty)$, имеют вид

$$ \begin{equation} \biggl(\frac{\partial}{\partial z}+B_n(z;s,t)\biggr)P_n(z;s,t) =\beta_n(s,t)A_n(z;s,t)P_{n-1}(z;s,t), \end{equation} \tag{2.1} $$
$$ \begin{equation} \biggl(\frac{\partial}{\partial z}-B_n(z;s,t)-\frac{\partial v(z,t)}{\partial z}\biggr)P_{n-1}(z;s,t) =-A_{n-1}(z;s,t)P_n(z;s,t), \end{equation} \tag{2.2} $$
где
$$ \begin{equation*} \begin{aligned} \, A_n(z;s,t)={}&\frac{P_n^2(y;s,t)w(y,t)}{(y-z)h_n(s,t)}\bigg|_{y=s}^{y=\infty}+{} \\ &+\frac{1}{h_n(s,t)}\int_s^{\infty}\frac{\partial v(z,t)/\partial z-\partial v(y,t)/\partial y}{z-y}P_n^2(y;s,t)w(y,t)\,dy, \\ B_n(z;s,t)={}&\frac{P_n(y;s,t)P_{n-1}(y;s,t)w(y,t)}{(y-z)h_{n-1}(s,t)}\bigg|_{y=s}^{y=\infty}+{} \\ &+\frac{1}{h_{n-1}(s,t)}\int_s^{\infty}\frac{\partial v(z,t)/\partial z-\partial v(y,t)/\partial y}{z-y}P_n(y;s,t)P_{n-1}(y;s,t)w(y,t)\,dy, \end{aligned} \end{equation*} \notag $$
$v(z,t)=-\ln w(z,t)$, а соответствующие условия совместности имеют вид
$$ \begin{equation} B_n(z;s,t)+B_{n+1}(z;s,t)=(z-\alpha_n(s,t))A_n(z;s,t)-\frac{\partial v(z,t)}{\partial z}, \end{equation} \tag{$\mathrm{S}_1$} $$
$$ \begin{equation} 1+(z-\alpha_n(s,t))(B_{n+1}(z;s,t)-B_n(z;s,t))=\beta_{n+1}(s,t)A_{n+1}(z;s,t)-{} \nonumber \end{equation} \tag{$\mathrm{S}_2$} $$
$$ \begin{equation} \qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad-\beta_n(s,t)A_{n-1}(z;s,t). \end{equation} \notag $$
Комбинация соотношений ($\mathrm{S}_1$) и ($\mathrm{S}_2$) приводит к выражению, в которое входят $\sum_{j=0}^{n-1}A_j(z;s,t)$, $B_n(z;s,t)$ и $\partial v(z,t)/\partial z$:
$$ \begin{equation} B_n^2(z;s,t)+\sum_{j=0}^{n-1}A_j(z;s,t)+B_n(z;s,t)\frac{\partial v(z,t)}{\partial z}=\beta_n(s,t)A_n(z;s,t)A_{n-1}(z;s,t). \end{equation} \tag{$\mathrm{S}'_2$} $$

Замечание 2. Для упрощения обозначений мы будем использовать $\partial_xf(x,y)$ и $\partial_{xx}f(x,y)$ вместо $\partial f(x,y)/\partial x$ и $\partial^2 f(x,y)/\partial x^2$ соответственно. Более полное объяснение и вывод соотношений (2.1), (2.2), ($\mathrm{S}_1$), ($\mathrm{S}_2$) и ($\mathrm{S}'_2$) можно найти в приложении.

Для веса (1.1) имеем

$$ \begin{equation*} \frac{\partial_z v(z,t)-\partial_y v(y,t)}{z-y}=2, \end{equation*} \notag $$
а из определения $A_n(z;s,t)$ и $B_n(z;s,t)$ найдем
$$ \begin{equation} A_n(z;s,t) =\frac{R_n(s,t)}{z-s}+2, \end{equation} \tag{2.3} $$
$$ \begin{equation} B_n(z;s,t) =\frac{r_n(s,t)}{z-s}, \end{equation} \tag{2.4} $$
где вспомогательные величины $R_n(s,t)$ и $r_n(s,t)$ вводятся по формулам
$$ \begin{equation*} R_n(s,t)=\frac{P_n^2(s;s,t)w(s,t)}{h_n(s,t)},\qquad r_n(s,t)=\frac{P_n(s;s,t)P_{n-1}(s;s,t)w(s,t)}{h_{n-1}(s,t)}. \end{equation*} \notag $$

Замечание 3. Пара лестничных операторов и условия совместности интенсивно обсуждаются в литературе (см., например, [18], [20]–[24]). Тождества ($\mathrm{S}_1$), ($\mathrm{S}_2$) и ($\mathrm{S}'_2$) справедливы для $z\in \mathbb{C}\cup\{\infty\}$.

Отметим, что зависимость от $t$ и $s$ через вес $e^{-z^2+tz}$, $z\in [s,\infty)$, $t>0$, $s>0$, индуцирует зависимость от $t$ и $s$ в полиномах $P_n(z;s,t)$, $h_j(s,t)$ и в связанных с ними величинах. Для краткости ниже мы не всегда явно указываем зависимость от $t$ и $s$, если в этом нет необходимости. Мы считаем, что это не приведет к недоразумениям.

Имеют место следующие соотношения:

$$ \begin{equation} r_0(s,t)=0,\qquad R_0(s,t)=\frac{P_0^2(s;s,t)w(s,t)}{h_0(s,t)}=\frac{2\pi^{-1/2}e^{-s^2+ts-t^2/4}}{1+\operatorname{erf}[(t-2s)/2]}, \end{equation} \tag{2.5} $$
где $\operatorname{erf}(z)$ – функция ошибок.

Подставляя (2.3) и (2.4) в ($\mathrm{S}_1$) и ($\mathrm{S}_2$), получим

$$ \begin{equation} R_n =2\alpha_n-t, \end{equation} \tag{2.6} $$
$$ \begin{equation} r_n+r_{n+1} =(s-\alpha_n)R_n, \end{equation} \tag{2.7} $$
$$ \begin{equation} r_{n+1}-r_n =2\beta_{n+1}-2\beta_n-1, \end{equation} \tag{2.8} $$
$$ \begin{equation} (\alpha_n-s)(r_{n+1}-r_n) =\beta_nR_{n-1}-\beta_{n+1}R_{n+1}. \end{equation} \tag{2.9} $$
Аналогично, подставляя (2.3) и (2.4) в ($\mathrm{S}'_2$), получим
$$ \begin{equation} r_n =2\beta_n-n, \end{equation} \tag{2.10} $$
$$ \begin{equation} \sum_{j=0}^{n-1}R_j =2\beta_nR_{n-1}+2\beta_nR_n+tr_n-2sr_n, \end{equation} \tag{2.11} $$
$$ \begin{equation} r_n^2 =\beta_nR_nR_{n-1}. \end{equation} \tag{2.12} $$
Из телескопической суммы (2.6) с учетом (1.10) получаем
$$ \begin{equation} \sum_{j=0}^{n-1}R_j=-2p(n)-nt. \end{equation} \tag{2.13} $$
Используя (2.6) и (2.7), чтобы исключить $\alpha_n$, получим
$$ \begin{equation} r_n+r_{n+1}=\biggl(s-\frac{t+R_n}{2}\biggr)R_n. \end{equation} \tag{2.14} $$

Предложение 1. Рекуррентные коэффициенты $\alpha_n$ и $\beta_n$ удовлетворяют системе нелинейных дифференциальных уравнений

$$ \begin{equation} 2\beta_n-2n =(s-\alpha_n)(2\alpha_n-t)-2\beta_{n+1}+1, \end{equation} \tag{2.15а} $$
$$ \begin{equation} (\alpha_n-s)(2\beta_{n+1}-2\beta_n-1) =(2\alpha_{n-1}-t)\beta_n-(2\alpha_{n+1}-t)\beta_{n+1} \end{equation} \tag{2.15б} $$
с начальными условиями
$$ \begin{equation} \alpha_0={} \frac{\pi^{-1/2}e^{-s^2+ts-t^2/4}}{1+\operatorname{erf}[(t-2s)/2]}+\frac{t}{2},\qquad \beta_0=0, \end{equation} \tag{2.16} $$
$$ \begin{equation} \alpha_1={} \frac{t}{2}+\biggl(s-\frac{t}{2}-\frac{\pi^{-1/2}e^{-s^2+ts-t^2/4}}{1+\operatorname{erf}[(t-2s)/2]}\biggr)\biggl\{1-\biggl[\biggl(\frac{\pi^{-1/2}e^{-s^2+ts-t^2/4}}{1+\operatorname{erf}[(t-2s)/2]}\biggr)\times {} \nonumber \end{equation} \notag $$
$$ \begin{equation} \qquad\qquad\qquad\qquad\qquad\times \biggl(s-\frac{t}{2}-\frac{\pi^{-1/2}e^{-s^2+ts-t^2/4}}{1+\operatorname{erf}[(t-2s)/2]}\biggr)+1\biggr]^{-1}\biggr\}, \end{equation} \tag{2.17} $$
$$ \begin{equation} \beta_1={} \biggl(s-\frac{t}{2}-\frac{\pi^{-1/2}e^{-s^2+ts-t^2/4}}{1+\operatorname{erf}[(t-2s)/2]}\biggr)\biggl(\frac{\pi^{-1/2}e^{-s^2+ts-t^2/4}}{1+\operatorname{erf}[(t-2s)/2]}\biggr)+\frac{1}{2}. \end{equation} \tag{2.18} $$

Доказательство. Вычитая (2.7) из (2.8), получаем

$$ \begin{equation} 2r_n=(s-\alpha_n)R_n-2\beta_{n+1}+2\beta_n+1. \end{equation} \tag{2.19} $$
Тогда, заменяя $r_n$ и $R_n$ соответственно выражениями (2.10) и (2.6), из (2.19) получим (2.15а).

Исключая $r_{n+1}-r_n$ из (2.8) и (2.9),

$$ \begin{equation} (\alpha_n-s)(2\beta_{n+1}-2\beta_n-1)=\beta_nR_{n-1}-\beta_{n+1}R_{n+1}, \end{equation} \tag{2.20} $$
после подстановки $2\alpha_{n-1}-t$ и $2\alpha_{n+1}-t$ вместо $R_{n-1}$ и $R_{n+1}$ в (2.20) в силу (2.6) приходим к (2.15б).

Из (2.6) и (2.10) при $n=0$ с помощью (2.5) получим (2.16). В частном случае $n=0$ уравнения (2.15а) и (2.15б) принимают вид

$$ \begin{equation*} \begin{aligned} \, &(s-\alpha_0)(2\alpha_0-t)-2\beta_1+1=0,\\ &(\alpha_0-s)(2\beta_1-1)+(2\alpha_1-t)\beta_1=0, \end{aligned} \end{equation*} \notag $$
откуда, комбинируя полученную систему уравнений с (2.16), получим (2.17) и (2.18).

3. Дифференциальные уравнения и уравнение Пенлеве IV

Начнем с вычисления частных производных $h_n(s,t)$ по $t$ и $s$ исходя из (1.4) при $j=k=n$,

$$ \begin{equation*} h_n(s,t)=\int_s^{\infty}P_n^2(z;s,t)w(z,t)\,dz, \end{equation*} \notag $$
что приводит соответственно к уравнениям
$$ \begin{equation*} \partial_t h_n(s,t)=\alpha_n(s,t)h_n(s,t),\qquad \partial_s h_n(s,t)=-R_n(s,t)h_n(s,t). \end{equation*} \notag $$
Следовательно, имеем
$$ \begin{equation} \partial_t\ln h_n(s,t) =\alpha_n(s,t), \end{equation} \tag{3.1} $$
$$ \begin{equation} \partial_s\ln h_n(s,t) =-R_n(s,t). \end{equation} \tag{3.2} $$
Отсюда следуют соотношения
$$ \begin{equation} \partial_s\ln \beta_n(s,t) =-R_n(s,t)+R_{n-1}(s,t), \end{equation} \tag{3.3} $$
$$ \begin{equation} \partial_t\ln \beta_n(s,t) =\alpha_n(s,t)-\alpha_{n-1}(s,t), \end{equation} \tag{3.4} $$
где использовано выражение (1.9).

Снова взяв частные производные по $t$ и $s$ от обеих частей равенства

$$ \begin{equation*} \int_s^{\infty}P_n(z;s,t)P_{n-1}(z;s,t)w(z,t)\,dz=0, \end{equation*} \notag $$
получим
$$ \begin{equation} \partial_t p(n,s,t) =-\beta_n(s,t), \end{equation} \tag{3.5} $$
$$ \begin{equation} \partial_s p(n,s,t) =r_n(s,t), \end{equation} \tag{3.6} $$
откуда в силу (1.8) находим
$$ \begin{equation} \partial_s\alpha_n(s,t) =r_n(s,t)-r_{n+1}(s,t), \end{equation} \tag{3.7} $$
$$ \begin{equation} \partial_t\alpha_n(s,t) =-\beta_n(s,t)+\beta_{n+1}(s,t). \end{equation} \tag{3.8} $$

Теперь мы готовы доказать следующую лемму.

Лемма 1. Вспомогательные величины $R_n(s,t)$ и $r_n(s,t)$ удовлетворяют четырем уравнениям в частных производных первого порядка:

$$ \begin{equation} \partial_sR_n(s,t) =4r_n(s,t)-(2s-t-R_n(s,t))R_n(s,t), \end{equation} \tag{3.9а} $$
$$ \begin{equation} \partial_tR_n(s,t) =\biggl(s-\frac{t+R_n(s,t)}{2}\biggr)R_n(s,t)-2r_n(s,t), \end{equation} \tag{3.9б} $$
$$ \begin{equation} \partial_sr_n(s,t) =\frac{2r_n^2(s,t)}{R_n(s,t)}-(r_n(s,t)+n)R_n(s,t), \end{equation} \tag{3.9в} $$
$$ \begin{equation} \partial_tr_n(s,t) =\frac{1}{2}(r_n(s,t)+n)R_n(s,t)-\frac{r_n^2(s,t)}{R_n(s,t)}. \end{equation} \tag{3.9г} $$

Доказательство. Продифференцируем (2.6) по $s$ и используем (3.7):

$$ \begin{equation} \partial_sR_n(s,t)=2r_n(s,t)-2r_{n+1}(s,t). \end{equation} \tag{3.10} $$
Комбинируя (3.10) и (2.14), чтобы исключить $r_{n+1}(s,t)$, приходим к (3.9а).

Аналогично, дифференцируя (2.6) по $t$ и подставляя (3.8), получим

$$ \begin{equation} \partial_tR_n(s,t)=2\partial_t\alpha_n(s,t)-1=2\beta_{n+1}(s,t)-2\beta_n(s,t)-1. \end{equation} \tag{3.11} $$
Используя далее (2.10), получим
$$ \begin{equation} \partial_tR_n(s,t)=r_{n+1}(s,t)-r_n(s,t). \end{equation} \tag{3.12} $$
Избавляясь от $r_{n+1}(s,t)$ с помощью (3.12) и (2.14), приходим к (3.9б).

Взяв производные (2.10) по $s$ и $t$, найдем

$$ \begin{equation} \partial_sr_n(s,t) =2\partial_s\beta_n(s,t), \end{equation} \tag{3.13} $$
$$ \begin{equation} \partial_tr_n(s,t) =2\partial_t\beta_n(s,t). \end{equation} \tag{3.14} $$
Подставляя (3.13) в (3.3), имеем
$$ \begin{equation} \partial_sr_n(s,t)=2\beta_n(s,t)R_{n-1}(s,t)-2\beta_n(s,t)R_n(s,t). \end{equation} \tag{3.15} $$
Согласно (3.4) и (2.6) из (3.14) следует
$$ \begin{equation} \partial_tr_n(s,t)=2(\alpha_n(s,t)-\alpha_{n-1}(s,t))\beta_n(s,t)=\beta_n(s,t)R_n(s,t)-\beta_n(s,t)R_{n-1}(s,t). \end{equation} \tag{3.16} $$
Исключая $\beta_n(s,t)R_{n-1}(s,t)$ и $\beta_n(s,t)$ из (3.15) и (3.16) с помощью (2.12) и (2.10), получим соответственно (3.9в) и (3.9г).

Предложение 2. Очевидно, из леммы 1 следует, что $R_n$ и $r_n$ удовлетворяют следующему дифференциальному уравнению:

$$ \begin{equation} \partial_sf(s,t)+2\partial_t f(s,t)=0. \end{equation} \tag{3.17} $$

Опираясь на лемму 1, покажем, что вспомогательная величина $R_n(s,t)$ удовлетворяет трем дифференциальным уравнениям в частных производных второго порядка.

Теорема 1. При $n\in\{0,1,2,\dots\}$ вспомогательная величина $R_n(s,t)$ удовлетворяет уравнениям

$$ \begin{equation} \partial_{st}R_n+\frac{(\partial_sR_n)^2}{4R_n}+\frac{3}{4}R_n^3+(t-2s)R_n^2+\biggl(\frac{1}{4}t^2-st+s^2-2n-1\biggr)R_n=0, \end{equation} \tag{3.18а} $$
$$ \begin{equation} \partial_{ss}R_n-\frac{(\partial_sR_n)^2}{2R_n}-\frac{3}{2}R_n^3+(4s-2t)R_n^2-\biggl(\frac{1}{2}t^2-2st+2s^2-4n-2\biggr)R_n=0, \end{equation} \tag{3.18б} $$
$$ \begin{equation} \partial_{tt}R_n-\frac{(\partial_tR_n)^2}{2R_n}-\frac{3}{8}R_n^3+\biggl(s-\frac{1}{2}t\biggr)R_n^2-\biggl(\frac{1}{8}t^2-\frac{1}{2}st+\frac{1}{2}s^2-n-\frac{1}{2}\biggr)R_n=0. \end{equation} \tag{3.18в} $$

Доказательство. Выразим $r_n$ из (3.9а):

$$ \begin{equation} r_n=\frac{1}{4}[\partial_sR_n+(2s-t-R_n)R_n]. \end{equation} \tag{3.19} $$
Вывод уравнений (3.18а)(3.18в) заключается в подстановке (3.19) в (3.9г) и (3.9в) и использовании уравнения (3.17) для упрощения итоговых выражений. Уравнение (3.18а) получим подстановкой (3.19) в (3.9г) при использовании (3.17) для упрощения выражения. Уравнение (3.18б) получим подстановкой (3.19) в (3.9в). Уравнение (3.18в) получим исключением $r_n$ с помощью комбинации (3.9б) и (3.9г).

Замечание 4. Подстановка выражения для $r_n$, полученного из (3.9б), в (3.9в) также приводит к (3.18а). Уравнения (3.9а)(3.9г) также можно преобразовать с помощью (3.17).

Обсуждение дифференциальных уравнений второго порядка, которым удовлетворяет $R_n(s,t)$, можно продолжить после введения новых переменных. Как указано в замечании 3, значение величины $R_n(s,t)$ зависит как от $s$, так и от $t$. Следующая теорема получена путем контроля этих переменных, при этом вводятся величины $R_n(s)$ и $R_n(t)$. Здесь $R_n(s)$ представляет случай, когда в $R_n(s,t)$ переменной считается $s$, а значение $t$ фиксировано. Для удобства используем в этом случае сокращенное обозначение $R_n(s)$ для $R_n(s,t)$. Аналогично вводится обозначение $R_n(t)$.

Теорема 2. Введем $S:=t/2-s$. Для удобства обозначим $\mathbb{R}_n(S)$ через $\mathbb{R}_n$. Эта величина удовлетворяет следующему дифференциальному уравнению второго порядка:

$$ \begin{equation} \partial_{SS}\mathbb{R}_n-\frac{(\partial_S\mathbb{R}_n)^2}{2\mathbb{R}_n}-\frac{3}{2}\mathbb{R}_n^3-4S\mathbb{R}_n^2-2(S^2-2n-1)\mathbb{R}_n=0, \end{equation} \tag{3.20} $$
которое известно как уравнение Пенлеве IV, точнее, $P_\mathrm{IV}(2n+1,0)$. В этом уравнении $\mathbb{R}_n(S)$ эквивалентно $R_n(s)$. Аналогично, вводя $T:=t/2-s$, получим величину $\mathbf{R}_n(T)$, удовлетворяющую уравнению Пенлеве IV$(2n+1,0)$:
$$ \begin{equation} \begin{aligned} \, \partial_{TT}\mathbf{R}_n-\frac{(\partial_T\mathbf{R}_n)^2}{2\mathbf{R}_n}-\frac{3}{2}\mathbf{R}_n^3-4T\mathbf{R}_n^2-2(T^2-2n-1)\mathbf{R}_n=0, \end{aligned} \end{equation} \tag{3.21} $$
где $\mathbf{R}_n(T)$ соответствует $R_n(t)$.

Доказательство. Для вывода уравнений (3.20) и (3.21) подставим $s=t/2-S$ и $t=2T+2s$ в (3.18б) и (3.18в) соответственно. Тогда, проделав вычисления и упрощения, получим искомые результаты.

Замечание 5. Два параметра $s$ и $t$ в (1.1) оказывают влияние как на величину $\mathbb{R}_n(S)$, так и на $\mathbf{R}_n(T)$, которые удовлетворяют уравнению Пенлеве IV$(2n+1,0)$ благодаря замене переменных. Говоря конкретнее, заменив $s$ на $t/2-S$, получим $\mathbb{R}_n(S)$, а заменив $t$ на $2T+2s$, получим $\mathbf{R}_n(T)$. Эти замены позволяют описать связь между $s$ и $t$ как “зеркальные изображения” друг друга, поскольку зависимость $\mathbb{R}_n(S)$ от $s$ и $t$ отражается в зависимость $\mathbf{R}_n(T)$ от $t$ и $s$ соответственно. Эта связь играет решающую роль при анализе асимптотического поведения определителя Ганкеля, что обсуждается в разделе 4.

В конце настоящего раздела мы используем некоторые из полученных уравнений, включающие в себя рекуррентные коэффициенты, для получения следующих результатов.

Предложение 3. Уравнения в частных производных, которым удовлетворяют $\alpha_n$ и $\beta_n$, можно выразить через $\alpha_n$ и $\beta_n$, т. е.

$$ \begin{equation} (s\partial_s+t\partial_t)\alpha_n=(2s-t)\biggl[2\beta_n-\frac{1}{2}(s-\alpha_n)(2\alpha_n-t)-\frac{1}{2}-n\biggr]+s, \end{equation} \tag{3.22а} $$
$$ \begin{equation} (s\partial_s+t\partial_t)\ln \beta_n(s,t)=\biggl(s-\frac{t}{2}\biggr)\biggl[\frac{(2\beta_n-n)^2}{(2\alpha_n-t)\beta_n}+t-2\alpha_n\biggr], \end{equation} \tag{3.22б} $$
$$ \begin{equation} (s^2\partial_{ss}+2st\partial_{st}+t^2\partial_{tt})\ln \beta_n\!=\biggl\{\frac{1}{2}\biggl[s-\frac{t}{2}-\frac{(2\beta_n-n)^2}{2(2\alpha_n-t)\beta_n}\biggr]\frac{(2\beta_n-n)^2}{(2\alpha_n-t)\beta_n}+{} \nonumber \end{equation} \notag $$
$$ \begin{equation} \qquad\qquad\qquad\qquad\qquad\qquad \qquad +2n-4\beta_n+\frac{1}{2}(s-\alpha_n)(2\alpha_n-t)\biggr\}(2s-t)^2. \end{equation} \tag{3.22в} $$

Доказательство. Применим $s\partial_s+t\partial_t$ к обеим частям (1.8) и заменим $p(n,s,t)$ на $p(n+1,s,t)$ с помощью (3.5) и (3.6), получим

$$ \begin{equation} \begin{aligned} \, (s\partial_s+t\partial_t)\alpha_n(s,t)&=s(r_n(s,t)-r_{n+1}(s,t))+t(\beta_{n+1}(s,t)-\beta_n(s,t))={} \notag\\ &=(2s-t)(\beta_n(s,t)-\beta_{n+1}(s,t))+s, \end{aligned} \end{equation} \tag{3.23} $$
где второе равенство следует из (2.10). Исключая $\beta_{n+1}$ из (3.23) и (2.15а), получим (3.22а).

Повторяя те же шаги, что и при выводе (3.22а), используя (3.3) и (3.4), получим

$$ \begin{equation*} \begin{aligned} \, (s\partial_s+t\partial_t)\ln\beta_n(s,t)&=s\partial_s\ln\beta_n(s,t)+t\partial_t\ln\beta_n(s,t)={}\\ &=s(R_{n-1}(s,t)-R_n(s,t))+t(\alpha_n(s,t)-\alpha_{n-1}(s,t)). \end{aligned} \end{equation*} \notag $$
Затем, исключая $R_{n-1}(s,t)$ и $R_n(s,t)$ с помощью (2.6), приходим к уравнению
$$ \begin{equation} (s\partial_s+t\partial_t)\ln \beta_n(s,t)=(2s-t)(\alpha_{n-1}-\alpha_{n}). \end{equation} \tag{3.24} $$
Из (2.6), подставляя $n-1$ вместо $n$, получим
$$ \begin{equation*} \alpha_{n-1}=\frac{1}{2}(R_{n-1}+t). \end{equation*} \notag $$
Заменяя теперь $R_{n-1}$ на $r_n^2/\beta_nR_n$ (см. (2.12)), получим
$$ \begin{equation} \alpha_{n-1}=\frac{r_n^2}{2\beta_nR_n}+\frac{t}{2}=\frac{(2\beta_n-n)^2}{2(2\alpha_n-t)\beta_n}+\frac{t}{2}, \end{equation} \tag{3.25} $$
где мы использовали (2.10) и (2.6). После подстановки (3.25) в (3.24) приходим к (3.22б).

Снова применяя $s\partial_s+t\partial_t$ к (3.24), после упрощения получим

$$ \begin{equation} \begin{aligned} \, (s^2\partial_{ss}+2st\partial_{st}+t^2\partial_{tt})\ln\beta_n(s,t)={}&(2s-t)(s\partial_s\alpha_{n-1}(s,t)+t\partial_t\alpha_{n-1}(s,t))-{} \notag\\ &-(2s-t)(s\partial_s\alpha_n(s,t)+t\partial_t\alpha_n(s,t)). \end{aligned} \end{equation} \tag{3.26} $$
Непосредственно подставляя (3.23) в (3.26), имеем
$$ \begin{equation} (s^2\partial_{ss}+2st\partial_{st}+t^2\partial_{tt})\ln \beta_n(s,t)=(2s-t)^2(\beta_{n-1}-2\beta_{n}+\beta_{n+1}). \end{equation} \tag{3.27} $$
Из (2.15а) с индексом $n-1$ следует
$$ \begin{equation} \beta_{n-1}=\frac{1}{2}(s-\alpha_{n-1})(2\alpha_{n-1}-t)+n-\beta_n-\frac{1}{2}, \end{equation} \tag{3.28} $$
а после подстановки (3.25) в (3.28), чтобы исключить $\alpha_{n-1}$, получим
$$ \begin{equation} \beta_{n-1}=\frac{1}{2}\biggl[s-\frac{t}{2}-\frac{(2\beta_n-n)^2}{2(2\alpha_n-t)\beta_n}\biggr]\frac{(2\beta_n-n)^2}{(2\alpha_n-t)\beta_n}+n-\frac{1}{2}-\beta_n. \end{equation} \tag{3.29} $$
Подставляя (3.29) и (2.15а) в (3.27) и исключая тем самым $\beta_{n-1}$ и $\beta_{n+1}$, приходим к (3.22в).

Замечание 6. Как и при выводе уравнений (3.18а)(3.18в), подставляя (3.9в) и (3.9г) в (3.9а) и (3.9б), получим систему уравнений в частных производных второго порядка для величины $r_n(s,t)$. Однако мы не приводим здесь эти уравнения, поскольку они не связаны напрямую с основными результатами и обсуждениями настоящей работы.

4. Асимптотики и определитель Ганкеля

В данном разделе применяется метод кулоновской жидкости для формулировки асимптотического поведения определителя Ганкеля, порождаемого весом (1.1), в логарифмическом виде. Метод кулоновской жидкости является мощным инструментом для изучения асимптотик определителей в теории случайных матриц и в связанных с ней областях.

Более подробную информацию о кулоновской жидкости, о связях между кулоновской жидкостью и ортогональными полиномами и о равновесной плотности можно найти в работах [6], [18], [19], [25]. Согласно работе [19] равновесная плотность $\sigma(x)$ с носителем на $J$ определяется минимизацией функционала свободной энергии

$$ \begin{equation*} F[\sigma]:=\int_{J}\sigma(x)v(x)\,dx-\int_{J}\int_{J}\sigma(x)\ln|x-y|\sigma(y)\,dx\,dy \end{equation*} \notag $$
при условии
$$ \begin{equation} \int_{J}\sigma(x)\,dx=n, \end{equation} \tag{4.1} $$
где $v(x)=-\ln w(x)$ – потенциал.

Отсюда следует, что плотность $\sigma(x)$ удовлетворяет интегральному уравнению

$$ \begin{equation} A=v(x)-2\int_J\ln|x-y|\sigma(y)\,dy,\qquad x\in J, \end{equation} \tag{4.2} $$
где $A$ – множитель Лагранжа для условия (4.1). Напомним, что $A$ является постоянной, не зависящей от $x$, но она зависит от $n$ и параметров потенциала $v(x)$. Взяв производную по $x$ от уравнения (4.2), получим сингулярное интегральное уравнение
$$ \begin{equation} v^{\prime}(x)-2\,\mathrm{P.V.}\int_J\frac{\sigma(y)}{x-y}\,dy=0,\qquad x\in J, \end{equation} \tag{4.3} $$
где $\mathrm{P.V.}$ означает главное значение в смысле Коши. В этом случае полученная для веса (1.1) на основе теории сингулярных интегральных уравнений [26] плотность $\sigma(x)$ с носителем на $[s,b]$ имеет вид
$$ \begin{equation*} \sigma(x)=\frac{1}{2\pi^2}\sqrt{\frac{b-x}{x-s}} \,\mathrm{P.V.} \int_s^b\frac{v^{\prime}(y)}{y-x}\,\sqrt{\frac{y-s}{b-y}}\,dy,\qquad x\in(s,b), \end{equation*} \notag $$
и является решением уравнения (4.3), при этом условие нормировки (4.1) принимает вид
$$ \begin{equation} \int_s^bv^{\prime}(x)\sqrt{\frac{x-s}{b-x}}\,dx=2\pi n. \end{equation} \tag{4.4} $$

Заметив для нашей задачи, что

$$ \begin{equation} v(x)=-\ln w(x)=x^2-tx, \end{equation} \tag{4.5} $$
получим
$$ \begin{equation*} \sigma(x)=\frac{1}{2\pi}(b-s+2x-t)\sqrt{\frac{b-x}{x-s}},\qquad x\in(s,b). \end{equation*} \notag $$
Подстановка (4.5) в (4.4) приводит к уравнению
$$ \begin{equation*} 3b^2-2(s+t)b-s^2+2ts-8n=0, \end{equation*} \notag $$
которое имеет интересующее нас положительное решение при $n\to\infty$:
$$ \begin{equation} \begin{aligned} \, b={}&2\,\sqrt{\frac{2n}{3}}+\frac{1}{3}(t+s)+\frac{1}{12\,\sqrt{6}}(t^2-4st+4s^2)n^{-1/2}-\frac{1}{1152\,\sqrt{6}}\times{} \notag\\ &\times (t^2-4st+4s^2)^2n^{-3/2}+\frac{1}{55296\,\sqrt{6}}(t^2-4st+4s^2)^3n^{-5/2}+\mathcal{O}(n^{-7/2}). \end{aligned} \end{equation} \tag{4.6} $$
Прямыми вычислениями получим
$$ \begin{equation} \frac{b+s}{2}={} \sqrt{\frac{2n}{3}}+\frac{2s}{t}+\frac{t}{6}+\frac{1}{24\,\sqrt{6}}(t-2s)^2n^{-1/2}-\frac{1}{2304\,\sqrt{6}}(t-2s)^4n^{-3/2}+{} \nonumber \end{equation} \tag{4.7} $$
$$ \begin{equation} +\frac{1}{110592\,\sqrt{6}}(t-2s)^6n^{-5/2}+\mathcal{O}(n^{-7/2}), \end{equation} \notag $$
$$ \begin{equation} \biggl(\frac{b-s}{2}\biggr)^2={} \frac{n}{6}+\frac{1}{6\,\sqrt{6}}(t-2s)n^{1/2}+\frac{1}{72}(t-2s)^2+\frac{1}{288\,\sqrt{6}}(t-2s)^3n^{-1/2}-{} \nonumber \end{equation} \tag{4.8} $$
$$ \begin{equation} -\frac{1}{27648\,\sqrt{6}}(t-2s)^5n^{-3/2}+\frac{1}{1327104\,\sqrt{6}}(t-2s)^7n^{-5/2}+\mathcal{O}(n^{-3}). \end{equation} \notag $$

Лемма 2. Множитель Лагранжа (4.2), связанный с весом (1.1), имеет вид

$$ \begin{equation} A=\frac{1}{8}(3b^2+2bs+3s^2)-\frac{1}{2}(b+s)t-2n\ln\biggl(\frac{b-s}{4}\biggr). \end{equation} \tag{4.9} $$
В пределе $n\to\infty$ получим
$$ \begin{equation} \begin{aligned} \, A={}&(1+\ln6-\ln n)n+\frac{\sqrt{6}}{3}(2s-t)n^{1/2}+\frac{1}{12}(8s^2-8st-t^2)+\frac{1}{72\,\sqrt{6}}(8s^3+{} \notag\\ &+6st^2-12s^2t-t^3)n^{-1/2}+\frac{1}{11520\,\sqrt{6}}(t^5-10st^4+40s^2t^3-80s^3t^2+{} \notag\\ &+80s^4t-32s^5)n^{-3/2}+\mathcal{O}(n^{-5/2}). \end{aligned} \end{equation} \tag{4.10} $$

Доказательство. Исходя из (4.2) имеем

$$ \begin{equation*} v(x)-2\int_s^b\sigma(y)\ln|x-y|\,dy=A,\qquad x\in(s,b). \end{equation*} \notag $$
Умножая обе части этого уравнения на $1/\sqrt{(b-x)(x-s)}$ и интегрируя по $x$ от $s$ до $b$, получим
$$ \begin{equation} \int_s^b\frac{(x^2-tx)}{\sqrt{(b-x)(x-s)}}\,dx-2\int_s^b dy\, \sigma(y)\int_s^b\frac{\ln|x-y|}{\sqrt{(b-x)(x-s)}}\,dx=A\pi, \end{equation} \tag{4.11} $$
где
$$ \begin{equation*} \int_s^b\frac{dx}{\sqrt{(b-x)(x-s)}}=\pi. \end{equation*} \notag $$

Заметим, что $A$ – константа, которая не зависит от $x$ при $x\in[s,b]$, но зависит от $n$, поэтому можно заменить $y$ на $b$,

$$ \begin{equation} \begin{aligned} \, \int_s^b\frac{\ln|x-y|}{\sqrt{(b-x)(x-s)}}\,dx&=\int_s^b\frac{\ln(b-x)}{\sqrt{(b-x)(x-s)}}\,dx={} \notag\\ &=\int_0^1\frac{\ln[(b-s)(1-X)]}{\sqrt{(1-X)X}}\,dX=\pi\ln\biggl(\frac{b-s}{4}\biggr), \end{aligned} \end{equation} \tag{4.12} $$
где при вычислении использовалось преобразование переменной $x:=(b-s)X+s$.

С другой стороны, проделав некоторые вычисления, преобразуем первый интеграл в (4.11) в следующее выражение:

$$ \begin{equation} \int_s^b\frac{x^2-tx}{\sqrt{(b-x)(x-s)}}\,dx=\frac{1}{8}\pi(3b^2+2bs+3s^2)-\frac{1}{2}\pi(b+s)t. \end{equation} \tag{4.13} $$
Подставляя (4.12), (4.13) и (4.1) в (4.11), после упрощения получим (4.9).

Подставляя (4.6) в (4.9), получим разложение для $A$ при $n\to\infty$, т. е. формулу (4.10).

Теорема 3. Асимптотические представления рекуррентных коэффициентов $\alpha_n(s,t)$ и $\beta_n(s,t)$ имеют следующий вид:

$$ \begin{equation} \alpha_n(s,t) =\sqrt{\frac{2n}{3}}+\frac{1}{6}(4s+t)+\frac{1}{24\,\sqrt{6}}[12+(t-2s)^2]n^{-1/2}+\mathcal{O}(n^{-3/2}), \end{equation} \tag{4.14а} $$
$$ \begin{equation} \beta_n(s,t) =\frac{n}{6}+\frac{\sqrt{6}}{36}(t-2s)n^{1/2}+\frac{1}{72}(t-2s)^2+\frac{1}{288\,\sqrt{6}}(t-2s)^3n^{-1/2}+\mathcal{O}(n^{-1}). \end{equation} \tag{4.14б} $$

Доказательство. Используя (4.7), (4.8) и (4.10) и опираясь на метод, изложенный в работе [19], заметим, что $\alpha_n(s,t)$ и $\beta_n(s,t)$ допускают следующие разложения при $n\to\infty$:

$$ \begin{equation} \alpha_n(s,t) =\sqrt{\frac{2n}{3}}+\sum_{j=0}^{\infty}\frac{a_j(s,t)}{n^{j/2}}, \end{equation} \tag{4.15а} $$
$$ \begin{equation} \beta_n(s,t) =\frac{n}{6}+\sum_{j=-1}^{\infty}\frac{b_j(s,t)}{n^{j/2}}. \end{equation} \tag{4.15б} $$
Подставляя (4.15а) и (4.15б) в систему разностных уравнений, которой удовлетворяют $\alpha_n(s,t)$ и $\beta_n(s,t)$, а именно в (2.15а) и (2.15б), и взяв предел при $n\to\infty$, можно рекуррентным образом определить коэффициенты разложения $a_j(s,t)$ и $b_j(s,t)$, сопоставляя коэффициенты перед степенями $n$:
$$ \begin{equation*} \begin{gathered} \, a_0(s,t)=\frac{1}{6}(4s+t),\qquad a_1(s,t)=\frac{1}{24\,\sqrt{6}}[12+(t-2s)^2],\qquad a_2(s,t)=0, \\ b_{-1}(s,t)=\frac{\sqrt{6}}{36}(t-2s),\qquad b_0(s,t)=\frac{1}{72}(t-2s)^2,\qquad b_1(s,t)=\frac{1}{288\,\sqrt{6}}(t-2s)^3, \end{gathered} \end{equation*} \notag $$
и т. д. Это завершает доказательство.

Теорема 4. Логарифм определителя Ганкеля, порожденного весом (1.1), имеет следующее разложение при больших $n$:

$$ \begin{equation} \begin{aligned} \, \ln D_n(s,t)={}&\frac{2\,\sqrt{6}}{9}(t-2s)n^{3/2}+\biggl(\frac{t^2}{12}+\frac{2}{3}st-\frac{2s^2}{3}\biggr)n+{}\notag\\ &+\sqrt{6}\bigg(\frac{t^3}{216}-\frac{st^2}{36}+\frac{s^2t}{18} -\frac{s^3}{27}\bigg)n^{1/2}-\frac{(t-2s)^4}{864}-{}\notag\\ &-\sqrt{6}\bigg(\frac{s^5}{1080}+\frac{s}{72}-\frac{t}{144}-\frac{s^4t}{432}+\frac{s^3t^2}{432} -\frac{s^2t^3}{864}+\frac{st^4}{3456}-\frac{t^5}{34560}\bigg)n^{-1/2}+{}\notag\\ &+C(s,t)+\mathcal{O}(n^{-1}), \end{aligned} \end{equation} \tag{4.16} $$
где $C(s,t)$ – произвольная функция от $s$, $t$ и $n$.

Доказательство. Применяя оператор $s\partial_s+t\partial_t$ к прологарифмированному уравнению (1.6)

$$ \begin{equation*} \ln D_n(s,t)=\sum_{j=0}^{n-1}\ln h_j(s,t), \end{equation*} \notag $$
получим
$$ \begin{equation} \begin{aligned} \, (s\partial_s+t\partial_t)\ln D_n(s,t)&=s\sum_{j=0}^{n-1}\partial_s\ln h_j(s,t)+t\sum_{j=0}^{n-1}\partial_t\ln h_j(s,t)={} \notag \\ &=-s\sum_{j=0}^{n-1}R_j(s,t)+t\sum_{j=0}^{n-1}\alpha_j(s,t) \end{aligned} \end{equation} \tag{4.17} $$
в силу (3.1) и (3.2).

Определим

$$ \begin{equation} H_n(s,t):=(s\partial_s+t\partial_t)\ln D_n(s,t), \end{equation} \tag{4.18} $$
тогда (4.17) можно преобразовать, используя (2.6) и (2.11), чтобы исключить $\alpha_j(s,t)$ и затем $\sum_{j=0}^{n-1}R_j(s,t)$, к следующему виду:
$$ \begin{equation} \begin{aligned} \, H_n(s,t)={}&\frac{1}{2}nt^2+\biggl(\frac{t}{2}-s\biggr)\sum_{j=0}^{n-1}R_j={} \notag\\ ={}&\frac{1}{2}nt^2+\biggl(\frac{t}{2}-s\biggr)[2\beta_nR_{n-1}+2\beta_nR_n+(t-2s)r_n]. \end{aligned} \end{equation} \tag{4.19} $$
Тогда, учитывая (2.12), (2.6) и (2.10), получим выражение для $H_n(s,t)$ через $\alpha_n$ и $\beta_n$:
$$ \begin{equation} \begin{aligned} \, H_n(s,t)=\frac{1}{2}nt^2+\biggl(\frac{t}{2}-s\biggr)\biggl[\frac{2(2\beta_n-n)^2}{2\alpha_n-t}+2(2\alpha_n-t)\beta_n+(t-2s)(2\beta_n-n)\biggr]. \end{aligned} \end{equation} \tag{4.20} $$
Подставляя (4.14а) и (4.14б) в (4.20), при $n\to\infty$ получим
$$ \begin{equation} \begin{aligned} \, H_n(s,t)={}&\frac{2}{3}\,\sqrt{\frac{2}{3}}(t-2s)n^{3/2}+\frac{1}{6}(t^2+8st-8s^2)n+\frac{\sqrt{6}}{72}(t-2s)^3\sqrt{n}+{} \notag\\ &+\frac{1}{216}(t-2s)(t^3-6st^2+12s^2t-8s^3)+\frac{\sqrt{6}}{6912}(t-2s)(16s^4-{} \notag\\ &-32s^3t+24s^2t^2-8st^3+t^4+48)n^{-1/2}+\mathcal{O}(n^{-1}). \end{aligned} \end{equation} \tag{4.21} $$
Принимая во внимание (4.18), после интегрирования (4.21) получим искомый результат (4.21).

Замечание 7. Заметим, что для обсуждаемой нами задачи вес (1.1) при $s=0$ соответствует полуклассическому весу Лагерра $w(x;t)=x^{\lambda}e^{-x^2+tx}$, $x\in\mathbb{R^{+}}$, при $\lambda= 0$ (см. [27]). Применяя теорему 1.8 из работы [27], выразим произвольную функцию, фигурирующую в (4.16), через асимптотическое представление для $t=s=0$:

$$ \begin{equation*} C(0,0)\sim\biggl(\frac{1}{2}n^2-\frac{1}{6}\biggr)\ln n-\frac{1}{4}(3+2\ln 6)n^2+n\ln(2\pi)+2\xi^{\prime}(-1)-\frac{1}{6}\ln 2+\frac{1}{8}\ln 3, \end{equation*} \notag $$
где $\xi(x)$ – дзета-функция Римана.

Предложение 4. Используя (2.12) и (2.10), чтобы исключить $\beta_nR_{n-1}$ и $\beta_n$ из (4.19), выразим $H_n(s,t)$ через $R_n(s,t)$ и $r_n(s,t)$ следующим образом:

$$ \begin{equation} H_n(s,t)=\frac{nt^2}{2}+\biggl(\frac{t}{2}-s\biggr)\biggl[\frac{2r_n^2}{R_n}+(r_n+n)R_n+(t-2s)r_n\biggr]. \end{equation} \tag{4.22} $$
Для функции $H_n(s,t)$ есть еще одно представление в терминах $R_n(s,t)$ и $\partial_sR_n(s,t)$:
$$ \begin{equation} \begin{aligned} \, H_n(s,t)={}&\frac{t(\partial_sR_n)^2}{16R_n}-\frac{s(\partial_sR_n)^2}{8R_n}+\biggl(\frac{s}{8}-\frac{t}{16}\biggr)R_n^3+\biggl(\frac{st}{2}-\frac{t^2}{8}-\frac{s^2}{2}\biggr)R_n^2+{} \notag\\ &+\biggl(\frac{3st^2}{8}-\frac{t^3}{16}-\frac{3s^2t}{4}+\frac{s^3}{2}+\frac{nt}{2}-ns\biggr)R_n+\frac{nt^2}{2}, \end{aligned} \end{equation} \tag{4.23} $$
если подставить (3.9а) в (4.22) и исключить таким образом $r_n$.

5. Другие численные характеристики

Опираясь на асимптотические характеристики, которые обсуждались в разделе 4, можно вычислить приближенное выражение для $P_n(0;s,t)$ при $n\to\infty$ и упростить дифференциальное уравнение второго порядка, которому удовлетворяют ортогональные полиномы, порожденные весом (1.1).

Ортогональные многочлены по весу (1.1) можно аппроксимировать формулой

$$ \begin{equation} P_n(z;s,t)\sim\exp[-F_1(z;s,t)-F_2(z;s,t)],\qquad z\notin [s,b],\quad n\to\infty, \end{equation} \tag{5.1} $$
где
$$ \begin{equation} \exp[-F_1(z;s,t)]=\frac{1}{2}\biggl[\biggl(\frac{z-b}{z-s}\biggr)^{1/4}+\biggl(\frac{z-s}{z-b}\biggr)^{1/4}\biggr], \end{equation} \tag{5.2а} $$
$$ \begin{equation} F_2(z;s,t)=-n\ln\biggl(\frac{\sqrt{z-s}+\sqrt{z-b}}{2}\,\biggr)^{\!2}+{} \nonumber \end{equation} \notag $$
$$ \begin{equation} \qquad\qquad\qquad\qquad+\frac{1}{2\pi}\int_s^b\frac{v(x)}{\sqrt{(b-x)(x-s)}}\biggl(\frac{\sqrt{(z-s)(z-b)}}{x-z}+1\biggr)\,dx, \end{equation} \tag{5.2б} $$
которая получена с помощью теории логарифмического потенциала с внешним полем, а также с помощью линейной статистической формулы (см. [28]–[31]). Следовательно, можно использовать эти уравнения для оценки $P_n(0;s,t)$ при больших $n$.

Теорема 5. Выражение для оценки полиномов $P_n(z;s,t)$ по весу (1.1) при $z=0$ имеет следующий вид:

$$ \begin{equation} \begin{aligned} \, (-1)^n & P_n(0;s,t)\sim{} 2^{-2n-5/8}3^{-1/8}s^{-1/4}(2^{3/4}3^{-1/4}n^{1/4}+s^{1/2})^{2n}\times{} \notag\\ &\times \exp\biggl[-\frac{n}{2}-2^{1/4}3^{-3/4}s^{1/2}n^{3/4}+\frac{(t-2s)n^{1/2}}{2\,\sqrt{6}}+2^{-9/4}3^{3/4}(t-s)s^{1/2}n^{1/4}+{} \notag\\ &\qquad\qquad\qquad\qquad\qquad+\frac{1}{24}(t^2+8st-8s^2)+\frac{1}{8}\ln n\biggr],\qquad n\to\infty. \end{aligned} \end{equation} \tag{5.3} $$

Доказательство. Очевидно,

$$ \begin{equation} \exp[-F_1(0;s,t)]=\frac{1}{2}\biggl[\biggl(\frac{b}{s}\biggr)^{\!1/4}+\biggl(\frac{s}{b}\biggr)^{\!1/4}\,\biggr]. \end{equation} \tag{5.4} $$
Вычислим интеграл:
$$ \begin{equation*} \begin{aligned} \, \int_s^b&\frac{x^2-tx}{\sqrt{(b-x)(x-s)}}\biggl(\frac{\sqrt{bs}}{x}+1\biggr)\,dx={} \\ &\qquad\qquad\qquad\qquad=\frac{\pi}{2}\,\sqrt{bs}(b+s)-t\pi\,\sqrt{bs}-\frac{t\pi}{2}(b+s)+\frac{\pi}{8}(3b^2+2bs+3s^2). \end{aligned} \end{equation*} \notag $$
Тогда из (5.2б) получим
$$ \begin{equation} \begin{aligned} \, -F_2(0;s,t)={}&\ln\biggl(\frac{\sqrt{-s}+\sqrt{-b}}{2}\biggr)^{2n}-\frac{1}{4}\,\sqrt{bs}(b+s)+\frac{t}{2}\,\sqrt{bs}+{} \notag \\ &+\frac{t}{4}(b+s)-\frac{1}{16}(3b^2+2bs+3s^2). \end{aligned} \end{equation} \tag{5.5} $$
Взяв ветвь $-s=s e^{i\pi}$ и $-b=b e^{i\pi}$ в (5.5) и подставляя (5.5) и (5.4) в (5.1), при $z=0$ получим
$$ \begin{equation} \begin{aligned} \, (-1)^n & P_n(0;s,t)\sim{} 2^{-2n-1}(\sqrt{s}+\sqrt{b})^{2n}\biggl[\biggl(\frac{b}{s}\biggr)^{\!1/4}+\biggl(\frac{s}{b}\biggr)^{\!1/4} \,\biggr]\times{} \notag\\ &\times\exp\biggl[\frac{t}{4}(b+s)+\frac{t}{2}\,\sqrt{bs}-\frac{1}{16}(3b^2+2bs+3s^2)-\frac{\sqrt{bs}}{4}(b+s)\biggr]. \end{aligned} \end{equation} \tag{5.6} $$

Подставляя (4.6) в (5.6) и устремляя $n$ к бесконечности, сразу получим (5.3).

Возвращаясь к паре лестничных операторов, исключим $P_{n-1}(z;s,t)$ из (2.1) и (2.2) с помощью ($\mathrm{S}'_2$). Это приводит к дифференциальному уравнению второго порядка, которому удовлетворяют ортогональные полиномы по весу (1.1) и которое имеет вид

$$ \begin{equation} \partial_{zz}P_n(z;s,t)+G_n(z;s,t)\partial_zP_n(z;s,t)+g_n(z;s,t)P_n(z;s,t)=0, \end{equation} \tag{5.7} $$
где
$$ \begin{equation} G_n(z;s,t)=-\partial_zv(z,t)-\frac{\partial_zA_n(z;s,t)}{A_n(z;s,t)}, \end{equation} \tag{5.8} $$
$$ \begin{equation} q_n(z,t_1,t_2)=\partial_zB_n(z;s,t)-\frac{B_n(z;s,t)\partial_zA_n(z;s,t)}{A_n(z;s,t)}+\sum_{j=0}^{n-1}A_j(z;s,t). \end{equation} \tag{5.9} $$
Подставляя (1.1), (2.3) и (2.4) в (5.8) и (5.9), используя затем (2.11), (2.6) и (2.10), чтобы исключить $\sum_{j=1}^{n-1}R_j$, $R_n$ и $r_n$, получим выражения
$$ \begin{equation} G_n(z;s,t)={} t-2z+\frac{2\alpha_n-t}{(s-z)(2s+t-2z-2\alpha_n)}, \end{equation} \tag{5.10} $$
$$ \begin{equation} g_n(z;s,t)={} 2n+\frac{2(2\alpha_{n-1}-t)\beta_n+2(2\alpha_n-t)\beta_n+(t-2s)(2\beta_n-n)}{z-s}-{} \nonumber \end{equation} \tag{5.11} $$
$$ \begin{equation} -\frac{2(2\beta_n-n)}{(s-z)(2s+t-2z-2\alpha_n)}, \end{equation} \notag $$
каждое из которых представлено через $\alpha_{n}$, $\alpha_{n-1}$ и $\beta_n$.

Теорема 6. В скейлинговом пределе ортогональные полиномы $P_n(z;s,t)$ по весу (1.1) удовлетворяют следующим дифференциальным уравнениям второго порядка: уравнению

$$ \begin{equation} \partial_{zz}\widetilde{P}_n(z;s,t)+\biggl(\frac{1}{z-s}+t-2z\biggr)\partial_z\widetilde{P}_n(z;s,t)+\frac{4\,\sqrt{6}n^{3/2}}{9(z-s)}\widetilde{P}_n(z;s,t)=0 \end{equation} \tag{5.12} $$
при $n\to\infty$, которое при $s=0$ эквивалентно биконфлюэнтному уравнению Гойна [32]
$$ \begin{equation*} \frac{d^2u(z)}{dz^2}+\biggl(\frac{\gamma}{z}+\delta+\varepsilon z\biggr)\frac{du(z)}{dz}+\frac{\alpha z-q}{z}u(z)=0, \end{equation*} \notag $$
где $\widetilde{P}_n(z;s,t)=u(z)$, $\gamma=1$, $\delta=t$, $\varepsilon=-2$, $\alpha=0$, $q=-4\,\sqrt{6}n^{3/2}/9$, и уравнению
$$ \begin{equation} \begin{aligned} \, \partial_{zz}\widetilde{P}_n(z;s,t)&+\biggl[\biggl(t-\frac{1}{s}\biggr)-\biggl(2+\frac{1}{s^2}\biggr)z-\frac{z^2}{s^3}\biggr]\partial_z\widetilde{P}_n(z;s,t)+{}\notag\\ &+\biggl[-\frac{4nz}{3s} -\frac{4\,\sqrt{6}n^{3/2}}{9s}+\frac{2}{3s}(s+t)n\biggr]\widetilde{P}_n(z;s,t)=0 \end{aligned} \end{equation} \tag{5.13} $$
при $n\to\infty$, $s\to\infty$, $t\to 0$, которое эквивалентно триконфлюэнтному уравнению Гойна [32]
$$ \begin{equation*} \frac{d^2u(z)}{dz^2}+(\gamma+\delta z+\varepsilon z^2)\frac{du(z)}{dz}+(\alpha z-q)u(z)=0, \end{equation*} \notag $$
где
$$ \begin{equation*} \begin{aligned} \, \widetilde{P}_n(z;s,t)&=u(z), \qquad \gamma=t-\frac{1}{s}, \qquad \delta=-2-\frac{1}{s^2}, \qquad \varepsilon=-\frac{1}{s^3}, \\ \alpha&=-\frac{4n}{3s}, \qquad q=\frac{4\,\sqrt{6}n^{3/2}}{9}-\frac{2}{3s}(s+t)n. \end{aligned} \end{equation*} \notag $$

Доказательство. Подставляя (4.14а) и (4.14б) в (5.10) и (5.11) и взяв предел при $n\to\infty$, получим (5.12), где опущены остаточные члены. Аналогично, рассматривая предел при $n\to\infty$, $s\to\infty$ и $t\to 0$ и проделав те же шаги, получим уравнение (5.13).

6. Выводы

Целью нашего исследования является изучение определителя Ганкеля, порожденного полуклассическим весом Лагерра с “щелевым” параметром $s$ и “временны́м” параметром $t$, как описано в выражении (1.1). Эту проблему можно считать расширением или продолжением исследования [14]. Используя лестничные операторы и связанные с ними условия совместности, мы ввели две вспомогательные величины, а именно $R_n(s,t)$ и $r_n(s,t)$. Эти величины введены в связи с полуклассическим ансамблем Лагерра с “щелевым” параметром $s$ и “временны́м” параметром $t$. Проведя анализ лестничных операторов и условий совместности, мы получили ряд нелинейных дифференциальных уравнений, которым удовлетворяют $R_n(s,t)$, $r_n(s,t)$ и рекуррентные коэффициенты в уравнении (1.7).

Вычисление частных производных уравнения (1.4) по $t$ и $s$ приводит к соответствующим уравнениям в частных производных первого порядка, которым удовлетворяют вспомогательные величины. Эти уравнения в частных производных в комбинации с полученными ранее нелинейными разностными уравнениями приводят к более глубокому пониманию поведения вспомогательной величины $R_n(s,t)$. После преобразования переменных дифференциальные уравнения второго порядка, котором удовлетворяет $R_n(s,t)$, оказались связанными с уравнением Пенлеве IV. Конкретные дифференциальные уравнения второго порядка, которым удовлетворяет $R_n(s,t)$, можно представить в виде (3.20) и (3.21). Интересно отметить, что уравнения (3.20) и (3.21) проявляют связь, описанную как “зеркальное изображение” друг друга. Кроме того, уравнения в частных производных, которым удовлетворяют $\alpha_n$ и $\beta_n$, можно выразить через них же.

Применяя метод кулоновской жидкости, мы получили асимптотики рекуррентных коэффициентов $\alpha_n(s,t)$ и $\beta_n(s,t)$ в уравнении (1.7) при больших $n$. При выводе асимптотических выражений рассматривался множитель Лагранжа, связанный с (1.1). Система нелинейных разностных уравнений (2.15а) и (2.15б) играет решающую роль в получении асимптотик рекуррентных коэффициентов при больших $n$. Эти разностные уравнения являются основой для изучения поведения и соотношений между коэффициентами, их использование в анализе увеличивает точность и надежность асимптотических результатов. Кроме того, найдено асимптотическое выражение для логарифмического представления определителя Ганкеля, порожденного весом (1.1), при больших $n$.

Путем комбинации результатов работы [31] с полученными ранее асимптотическими поведениями вычислены ортогональные многочлены $P_n(z;s,t)$, связанные с весом (1.1), при $z=0$. Показано, что линейное дифференциальное уравнение второго порядка для ортогональных многочленов в специальном скейлинговом пределе сводится к аналогу биконфлюэнтного уравнения Гойна и триконфлюэнтного уравнения Гойна. Сведение линейного дифференциального уравнения второго порядка для ортогональных многочленов к биконфлюэнтному уравнению Гойна и триконфлюэнтному уравнению Гойна позволяет провести дальнейший анализ и исследование свойств и поведения этих полиномов. Решения этих специальных дифференциальных уравнений углубляет понимание структуры и характеристик ортогональных многочленов, связанных с весом (1.1). Связь проведенного анализа ортогональных многочленов с результатами работы [31] и с асимптотическими поведениями, полученными ранее, дает более глубокое понимание полуклассического ансамбля Лагерра и связанных с ним ортогональных полиномов.

Приложение

Приведем вывод понижающего оператора (2.1) (см., например, [3], [6], [7], [11], [18], [19]).

Начнем с соотношения

$$ \begin{equation*} \partial_z P_n(z;s,t)=\sum_{k=0}^{n-1}c_{n,k}P_k(z;s,t), \end{equation*} \notag $$
где $c_{n,k}$ определяется из соотношений ортогональности
$$ \begin{equation*} c_{n,k}=\frac{1}{h_k(s,t)}\int_s^{\infty}P_k(y;s,t)[\partial_y P_n(y;s,t)]w(y,t)\,dy. \end{equation*} \notag $$
Непосредственно получим
$$ \begin{equation} \partial_z P_n(z;s,t)=\sum_{k=0}^{n-1}\frac{P_k(z;s,t)}{h_k(s,t)}\int_s^{\infty}P_k(y;s,t)[\partial_y P_n(y;s,t)]w(y,t)\,dy. \end{equation} \tag{П.1} $$
Интегрируя по частям, имеем
$$ \begin{equation} \begin{aligned} \, \int_s^{\infty}&P_k(y;s,t)[\partial_y P_n(y;s,t)]w(y,t)\,dy={} \notag \\ ={}&P_n(y;s,t)P_k(y;s,t)w(y,t)\big|_{y=s}^{y=\infty}-\int_s^{\infty}P_n(y;s,t)P_k(y;s,t)[\partial_y w(y,t)]\,dy={} \notag\\ ={}&P_n(y;s,t)P_k(y;s,t)w(y,t)\big|_{y=s}^{y=\infty}+\int_s^{\infty}P_n(y;s,t)P_k(y;s,t)w(y,t)[\partial_y v(y,t)]\,dy={} \notag\\ ={}&P_n(y;s,t)P_k(y;s,t)w(y,t)\big|_{y=s}^{y=\infty}-\int_s^{\infty}P_n(y;s,t)P_k(y;s,t)w(y,t)[\partial_z v(z,t)-{} \notag\\ &-\partial_y v(y,t)]\,dy, \end{aligned} \end{equation} \tag{П.2} $$
где используется равенство $v(z,t)=-\ln w(z,t)$. Подставляя (П.2) в (П.1), получим
$$ \begin{equation} \begin{aligned} \, \partial_z{}& P_n(z,s,t)=P_n(z;s,t)\sum_{k=0}^{n-1}\frac{P_k(z;s,t)P_k(y;s,t)w(y,t)}{h_k(s,t)}\bigg|_{y=s}^{y=\infty}-{} \notag\\ &-\int_s^{\infty}P_n(y;s,t)\sum_{k=0}^{n-1}\frac{P_k(z;s,t)P_k(y;s,t)}{h_k(s,t)}[\partial_zv(z,t)-\partial_yv(y,t)]w(y,t)\,dy={} \notag\\ ={}&\frac{P_n(z;s,t)P_n(y;s,t)P_{n-1}(y;s,t)w(y,t)}{(z-y)h_{n-1}(s,t)}\bigg|_{y=s}^{y=\infty}-\frac{P_n^2(y;s,t)P_{n-1}(z;s,t)w(y,t)}{(z-y)h_{n-1}(s,t)}\bigg|_{y=s}^{y=\infty}-{} \notag\\ &-\frac{P_n(z;s,t)}{h_{n-1}(s,t)}\int_s^{\infty}\frac{P_n(y;s,t)P_{n-1}(y;s,t)w(y,t)[\partial_z v(z,t)-\partial_y v(y,t)]}{z-y}\,dy+{} \notag\\ &+\frac{P_{n-1}(z;s,t)}{h_{n-1}(s,t)}\int_s^{\infty}\frac{P_n^2(y;s,t)w(y,t)[\partial_z v(z,t)-\partial_y v(y,t)]}{z-y}\,dy, \end{aligned} \end{equation} \tag{П.3} $$
где применена формула Кристоффеля–Дарбу [33]
$$ \begin{equation*} \sum_{k=0}^{n-1}\frac{P_k(z;s,t)P_k(y;s,t)}{h_k(s,t)}=\frac{P_n(z;s,t)P_{n-1}(y;s,t)-P_n(y;s,t)P_{n-1}(z;s,t)}{(z-y)h_{n-1}(s,t)}. \end{equation*} \notag $$
Далее, используя равенство $\beta_n(s,t)=h_n(s,t)/h_{n-1}(s,t)$, из (П.3) получим уравнение (2.1). Выполнив аналогичные действия, можно вывести уравнение (2.2).

Из определения $B_n(z;s,t)$ и рекуррентного соотношения получим

$$ \begin{equation} \begin{aligned} \, B_n(z;s,t)+{}&B_{n+1}(z;s,t) =\frac{[y-\alpha_n(s,t)]P_n^2(y;s,t)w(y,t)}{(y-z)h_n(s,t)}\bigg|_{y=s}^{y=\infty}+{} \notag\\ &+\int_s^{\infty}\frac{[y-\alpha_n(s,t)]P_n^2(y;s,t)w(y,t)[\partial_zv(z,t)-\partial_yv(y,t)]}{(z-y)h_n(s,t)}\,dy={} \notag\\ ={}&[z-\alpha_n(s,t)]A_n(z;s,t)-\partial_zv(z,t)+\frac{P_n^2(y;s,t)w(y,t)}{h_n(s,t)}\bigg|_{y=s}^{y=\infty}+{} \notag\\ &+\frac{1}{h_n(s,t)}\int_s^{\infty}P_n^2(y;s,t)w(y,t)[\partial_yv(y,t)]\,dy, \end{aligned} \end{equation} \tag{П.4} $$
где
$$ \begin{equation*} y-\alpha_n(s,t)=(z-\alpha_n(s,t))-(z-y). \end{equation*} \notag $$

Выпишем (П.2) при $k=n$:

$$ \begin{equation*} \int_s^{\infty}P_n(y;s,t)[\partial_y P_n(y;s,t)]w(y,t)\,dy=0. \end{equation*} \notag $$
Согласно предпоследнему равенству из (П.2) имеем
$$ \begin{equation} \int_s^{\infty}P_n^2(y;s,t)w(y,t)[\partial_y v(y,t)]\,dy+P_n^2(y;s,t)w(y,t)\big|_{y=s}^{y=\infty}=0. \end{equation} \tag{П.5} $$
Подстановка (П.5) в (П.4) завершает доказательство условия ($\mathrm{S}_1$). Аналогично доказываются условия ($\mathrm{S}_2$) и ($\mathrm{S}'_2$).

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

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

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

1. M. L. Mehta, Random Matrices, Elsevier, Amsterdam, 2004  mathscinet
2. E. Basor, Y. Chen, L. Zhang, “PDEs satisfied by extreme eigenvalues distributions of GUE and LUE”, Random Matrices Theory Appl., 1:1 (2012), 1150003, 21 pp.  crossref  mathscinet
3. E. Basor, Y. Chen, T. Ehrhardt, “Painlevé V and time-dependent Jacobi polynomials”, J. Phys. A: Math. Theor., 43:1 (2010), 015204, 25 pp.  crossref  mathscinet
4. Г. Сегё, Ортогональные многочлены, Физматгиз, М., 1962  mathscinet  zmath  zmath
5. Y. Chen, N. S. Haq, M. R. McKay, “Random matrix models, double-time Painlevé equations, and wireless relaying”, J. Math. Phys., 54:6 (2013), 063506, 55 pp.  crossref  mathscinet
6. Y. Chen, M. R. McKay, “Coulomb fluid, Painlevé transcendents and the information theory of MIMO systems”, IEEE Trans. Inform. Theory, 58:7 (2012), 4594–4634  crossref  mathscinet
7. S. Lyu, Y. Chen, E. Fan, “Asymptotic gap probability distributions of the Gaussian unitary ensembles and Jacobi unitary ensembles”, Nucl. Phys. B, 926 (2018), 639–670  crossref  mathscinet
8. C. Min, Y. Chen, “Gap probability distribution of the Jacobi unitary ensemble: an elementary treatment, from finite $n$ to double scaling”, Stud. Appl. Math., 140:2 (2018), 202–220  crossref  mathscinet
9. Y. Chen, D. Dai, “Painelvé V and a Pollaczek–Jacobi type orthogonal polynomials”, J. Approx. Theory, 162:12 (2010), 2149–2167  crossref  mathscinet
10. Y. Chen, L. Zhang, “Painlevé VI and the unitary Jacobi ensembles”, Stud. Appl. Math., 125:1 (2010), 91–112  crossref  mathscinet
11. D. Dai, L. Zhang, “Painlevé VI and Hankel determinants for the generized Jacobi weight”, J. Phys. A: Math. Theor., 43:5 (2010), 055207, 14 pp.  crossref  mathscinet
12. E. Basor, Y. Chen, “Painlevé V and the distribution function of a discontinuous linear statistics in the Laguerre unitary ensembles”, J. Phys. A: Math. Theor., 42:3 (2009), 035203, 18 pp.  crossref  mathscinet
13. E. Basor, Y. Chen, “Perturbed Laguerre unitary ensembles, Hankel determinants, and information theory”, Math. Methods Appl. Sci., 38:18 (2015), 4840–4851  crossref  mathscinet
14. S. Lyu, J. Griffin, Y. Chen, “The Hankel determinant associated with a singularly perturbed Laguerre unitary ensemble”, J. Nonlinear Math. Phys., 26:1 (2019), 24–53  crossref  mathscinet
15. Y. Chen, M. V. Feigin, “Painlevé IV and degenerate Gaussian unitary ensembles”, J. Phys. A: Math. Gen., 39:40 (2006), 12381–12393  crossref  mathscinet
16. Y. Chen, G. Pruessner, “Orthogonal polynomials with discontinuous weights”, J. Phys. A: Math. Gen., 38:12 (2005), L191–L198  mathscinet
17. P. A. Clarkson, K. Jordaan, “The relationship between semiclassical Laguerre polynomials and the fourth Painlevé equation”, Constr. Approx., 39:1 (2014), 223–254  crossref  mathscinet
18. Y. Chen, M. E. H. Ismail, “Ladder operators and differential equations for orthogonal polynomials”, J. Phys. A: Math. Gen., 30:22 (1997), 7817–7829  crossref  mathscinet
19. Y. Chen, M. E. H. Ismail, “Thermodynamic relations of the Hermitian matrix ensembles”, J. Phys. A: Math. Gen., 30:19 (1997), 6633–6654  crossref  mathscinet
20. W. C. Bauldry, “Estimate of the asymmetric Freud polynpomials on the real line”, J. Approx. Theory, 63:2 (1990), 225–237  crossref  mathscinet
21. S. Belmehdi, A. Ronveaux, “Laguerre–Freud's equations for the recurrence coefficients of semi-classical orthogonal polynomials”, J. Approx. Theory, 76:3 (1994), 351–368  crossref  mathscinet
22. S. Bonan, D. S. Clark, “Estimates of the orthogonal polynomials with weight $\exp(-x^m)$, $m$ an even positive integer”, J. Approx. Theory, 46:4 (1986), 408–410  crossref  mathscinet
23. S. Bonan, D. S. Clark, “Estimates of the Hermite and the Freud polynomials”, J. Approx. Theory, 63:2 (1990), 210–224  crossref  mathscinet
24. S. Bonan, P. Nevai, “Orthogonal polynomials and their derivatives. I”, J. Approx. Theory, 40:2 (1984), 134–147  crossref  mathscinet
25. F. J. Dyson, “Statistical theory of the energy levels of complex systems. I”, J. Math. Phys., 3:1 (1962), 140–156  crossref; “II”, 157–165  crossref; “III”, 166–175  crossref  mathscinet
26. E. B. Saff, V. Totik, Logarithmic Potentials with External Fields, Grundlehren der mathematischen Wissenschaften, 316, Springer, Berlin, Heidelberg, 1997  crossref  mathscinet
27. C. Min, Y. Chen, “Painlevé IV, Chazy II, and asymptotics for recurrence coefficients of semi-classical Laguerre polynomials and their Hankel determinants”, Math. Methods Appl. Sci., 46:14 (2023), 15270–15284  crossref
28. Y. Chen, D. S. Lubinsky, “Smallest eigenvalues of Hankel matrices for exponential weights”, J. Math. Anal. Appl., 293:2 (2004), 476–495  crossref  mathscinet
29. Y. Chen, N. Lawrence, “Small eigenvalues of large Hankel matrices”, J. Phys. A: Math. Gen., 32:42 (1999), 7305–7315  crossref  mathscinet
30. Y. Chen, N. Lawrence, “On the linear statistics of Hermitian random matrices”, J. Phys. A: Math. Gen., 31:4 (1998), 1141–1152  crossref  mathscinet
31. M. Chen, Y. Chen, “Singular linear statistics of the Laguerre unitary ensemble and Painlevé. III. Double scaling analysis”, J. Math. Phys., 56:6 (2015), 063506, 14 pp.  crossref  mathscinet
32. A. Ronveaux (ed.), Heun's Differential Equations, Oxford Univ. Press, Oxford, 1995  mathscinet  zmath
33. M. E. H. Ismail, Classical and Quantum Orthogonal Polynomials in One Variable, Encyclopedia of Mathematics and its Applications, 98, Cambridge Univ. Press, Cambridge, 2005  crossref  mathscinet

Образец цитирования: Дань Ван, “Определитель Ганкеля для полуклассического унитарного ансамбля Лагерра, уравнения Пенлеве IV и Гойна”, ТМФ, 219:3 (2024), 440–461; Theoret. and Math. Phys., 219:3 (2024), 913–932
Цитирование в формате AMSBIB
\RBibitem{Wan24}
\by Дань~Ван
\paper Определитель Ганкеля для полуклассического унитарного ансамбля Лагерра, уравнения Пенлеве IV и Гойна
\jour ТМФ
\yr 2024
\vol 219
\issue 3
\pages 440--461
\mathnet{http://mi.mathnet.ru/tmf10649}
\crossref{https://doi.org/10.4213/tmf10649}
\mathscinet{https://mathscinet.ams.org/mathscinet-getitem?mr=4767965}
\adsnasa{https://adsabs.harvard.edu/cgi-bin/bib_query?2024TMP...219..913W}
\transl
\jour Theoret. and Math. Phys.
\yr 2024
\vol 219
\issue 3
\pages 913--932
\crossref{https://doi.org/10.1134/S0040577924060035}
\scopus{https://www.scopus.com/record/display.url?origin=inward&eid=2-s2.0-85196786001}
Образцы ссылок на эту страницу:
  • https://www.mathnet.ru/rus/tmf10649
  • https://doi.org/10.4213/tmf10649
  • https://www.mathnet.ru/rus/tmf/v219/i3/p440
  • Citing articles in Google Scholar: Russian citations, English citations
    Related articles in Google Scholar: Russian articles, English articles
    Теоретическая и математическая физика Theoretical and Mathematical Physics
     
      Обратная связь:
     Пользовательское соглашение  Регистрация посетителей портала  Логотипы © Математический институт им. В. А. Стеклова РАН, 2025