Аннотация:
Рассматривается задача Коши для гамильтоновой системы, состоящей из поля Клейна–Гордона и бесконечного гармонического кристалла. Динамика системы трансляционно-инвариантна относительно дискретной подгруппы $\mathbb{Z}^d$ группы $\mathbb{R}^d$. Начальные данные задаются как случайная функция, которая близка к двум пространственно-однородным (относительно подгруппы $\mathbb{Z}^d$) процессам при $\pm x_1>a$ с некоторым $a>0$. Изучается распределение $\mu_t$ решения в момент времени $t\in\mathbb{R}$ и доказывается слабая сходимость $\mu_t$ к гауссовой мере $\mu_\infty$ при $t\to\infty$. Кроме того, доказано существование предела корреляционных функций и выведены явные формулы для ковариации предельной меры $\mu_\infty$. Дается приложение к случаю гиббсовских мер.
Ключевые слова:
поле Клейна–Гордона, взаимодействующее с гармоническим кристаллом, преобразование Зака, случайные начальные данные, гауссовы и гиббсовские меры, слабая сходимость мер.
Поступило в редакцию: 30.05.2023 После доработки: 30.05.2023
Данная статья является продолжением работ [1]–[6], посвященных изучению долговременно́го поведения распределений решений задачи Коши для бесконечномерных гамильтоновых систем. В качестве модели мы здесь рассматриваем линейную гамильтонову систему, состоящую из вещественного скалярного поля Клейна–Гордона $\psi(x)$ и его импульса $\pi(x)$, $x\in\mathbb{R}^d$, которое взаимодействует с гармоническим кристаллом, описываемым отклонениями $u(k)\in\mathbb{R}^n$ частиц (атомов, молекул, ионов и т. п.) от их положения равновесия и их скоростями $v(k)\in\mathbb{R}^n$, $k\in\mathbb{Z}^d$.
Функция взаимодействия $R(x)$ имеет значения в $\mathbb{R}^n$, $m_0>0$, точка обозначает скалярное произведение в $\mathbb{R}^n$ (ниже мы используем это обозначение также для скалярного произведения в $\mathbb{R}^d$), $V$ – вещественная матрица взаимодействия, $V(k)\in\mathbb{R}^n\times\mathbb{R}^n$, размерности $d,n\geqslant 1$. Вычисляя вариационные производные от $\mathrm H(\psi,u,\pi,v)$, приходим к следующей системе динамических уравнений:
Условия на матрицу $V$ и функцию $R$ приведены в п. 2.2 ниже. Система (1.2) может служить моделью движения блоховских электронов в периодической среде, порожденной ионными ядрами. Такое движение широко изучалось в литературе (см., например, труды [7]–[10] и библиографию в них).
Для системы (1.2) рассматривается задача Коши с начальными данными
Мы предполагаем, что начальные данные $Y_0$ задаются как случайная функция с распределением $\mu_0$. Обозначим через $\mathbb{E}$ математическое ожидание по мере $\mu_0$. На меру $\mu_0$ налагаются условия S1–S3 (см. п. 2.3 ниже). В частности, мера $\mu_0$ имеет нулевое среднее значение, а ее корреляционная матрица $Q_0(p,p'):=\mathbb{E}(Y_0(p)\otimes Y_0(p'))$, $p,p'\in\mathbb{P}^d$, имеет вид
Здесь $p= (p_1,\dots,p_d)$, $ p'= (p'_1,\dots,p'_d)\in\mathbb{P}^d$, $Q_{\pm}(p,p')$ – корреляционные матрицы некоторых мер $\mu_{\pm}$, трансляционно-инвариантных относительно группы $\mathbb{Z}^d$, т. е. для любого $k\in\mathbb{Z}^d$
Заметим, что $\mu_0$ не является инвариантной относительно группы $\mathbb{Z}^d$, если $Q_{-}\neq Q_{+}$.
Обозначим через $\mu_t$, $t\in\mathbb{R}$, вероятностную борелевскую меру, которая является распределением решения $Y(t)$ задачи (1.5) со случайными начальными данными $Y_0$. Мы изучаем асимптотику меры $\mu_t$ при $t\to\pm\infty$. Наша цель – доказать сходимость корреляционных функций меры $\mu_t$,
Мы показываем, что предельная мера $\mu_{\infty}$ является гауссовой и инвариантной относительно сдвигов в $\mathbb{Z}^d$, и выводим явные формулы для ковариации $Q_\infty$ меры $\mu_\infty$ (см. формулы (3.3)–(3.7) ниже).
В качестве примеров мы рассматриваем две начальные меры $\mu_0$. В первом случае мы предполагаем, что $\mu_{\pm}\equiv g_{\beta_{\pm}}$ (где $\beta_\pm=T_\pm^{-1}$) – это гиббсовские меры с различными температурами $T_\pm>0$ (определение см. в п. 4.1 ниже). Во втором случае предполагается, что начальные данные для поля и кристалла независимы при $t=0$, и, кроме того, в начальный момент времени “левые” и “правые” части поля и кристалла находятся в тепловом равновесии с различными температурами. В обоих случаях мы доказываем существование гиббсовских мер и получаем явные формулы для предельной корреляционной матрицы $Q_\infty$.
Для одномерных цепочек гармонических осцилляторов сходимость (1.9) была доказана в работах [11], [12]. Для одномерных цепочек негармонических осцилляторов, связанных с тепловыми резервуарами, соответствующие результаты были получены в [13], [14]. Для бесконечных многомерных гармонических кристаллов сходимость (1.9) была доказана в [2], [5], а для волновых уравнений и уравнений Клейна–Гордона – в [1], [4], [6]. В случае, когда начальные меры инвариантны относительно сдвигов в $\mathbb{Z}^d$, сходимость (1.9) была доказана в [3]. В настоящей статье мы обобщаем результаты работы [3] на более общий класс начальных мер $\mu_0$ и матриц взаимодействия $V$ в кристалле.
2. Постановка задачи
Введем фазовое пространство $\mathcal E^{s,\alpha}$ начальных данных $Y_0$.
Определение 2.1. 1. Пусть $H^{s,\alpha}=H^{s,\alpha}(\mathbb{R}^d)$, $s\in\mathbb{R}$, $\alpha\in\mathbb{R}$, – гильбертово пространство распределений $\psi\in S'(\mathbb{R}^d)$ с нормой $\|\psi\|_{s,\alpha}\equiv\|\langle x\rangle^\alpha\Lambda^s\psi\|_{L^2(\mathbb{R}^d)}<\infty$. Здесь
и $\hat\psi:=F \psi$ – преобразование Фурье функции $\psi$ умеренного роста. Для $\psi\in C_0^\infty(\mathbb{R}^d)$ полагаем $\hat\psi(\xi)=\int e^{\mathrm i\xi\cdot x}\psi(x)\,dx$. При $s=0,1,2,\dots{}$ пространство $H^{s,\alpha}$ – это гильбертово пространство вещественных функций $\psi(x)$ с нормой
2. Пусть $L^\alpha\equiv [\ell^2_\alpha(\mathbb{Z}^d)]^n$, $\alpha\in\mathbb{R}$, – гильбертово пространство векторнозначных функций $u(k)\in\mathbb{R}^n$, $k\in\mathbb{Z}^d$, с нормой
3. Пусть $\mathrm{E}^{s,\alpha}:=H^{s,\alpha}\oplus L^\alpha$, $s,\alpha\in\mathbb{R}$, – гильбертово пространство пар $Y^0\equiv(\psi,u)$ с нормой $\|Y^0\|^2_{s,\alpha}=\|\psi\|^2_{s,\alpha}+\|u\|^2_{\alpha}<\infty$.
4. Пусть $\mathcal E^{s,\alpha}:=\mathrm{E}^{s+1,\alpha}\oplus\mathrm{E}^{s,\alpha}\equiv H^{1+s,\alpha}\oplus L^\alpha\oplus H^{s,\alpha}\oplus L^\alpha$ – гильбертово пространство функций $Y=(Y^0,Y^1)$, где $Y^0=(\psi,u)$, $Y^1=(\pi,v)$, с нормой
Заметим, что $\mathcal E^{s_1,\alpha_1}\subset\mathcal E^{s_2,\alpha_2}$ для любых $s_2<s_1$ и $\alpha_2<\alpha_1$, и это вложение компактно по теореме вложения Соболева. Зафиксируем $\alpha<-d/2$. Ниже рассматриваются две задачи в $\mathcal E^{s,\alpha}$: с $s=0$ и $s<-d/2$. В случае $s=0$ доказывается сходимость (1.9) при условиях S1–S4 (см. п. 2.3 ниже) на начальные (вообще говоря, негауссовы) меры $\mu_0$. В случае $s<-d/2$ мы предполагаем, что $\mu_0$ – гауссова мера, удовлетворяющая условиям S1 и S3, где $\mu_\pm$ (см. формулу (1.6)) являются гиббсовскими мерами с температурами $T_\pm>0$.
Предположим сначала, что начальные данные $Y_0$ принадлежат фазовому пространству $\mathcal E:=\mathcal E^{0,\alpha}$.
2.1. Преобразование Блоха–Флоке–Зака
Динамика системы (1.5) инвариантна относительно пространственных сдвигов в $\mathbb{Z}^d$. Поэтому, применяя преобразование Зака, мы сведем систему к задаче Блоха (см., например, [7]) на торе. Для этого сначала представим $x\in\mathbb{R}^d$ как $x=l+y$, $l\in\mathbb{Z}^d$, $y\in\mathrm K_1^d:=[0,1]^d$, и применим дискретное преобразование Фурье к решениям $\psi(l+y,t)$, $u(k,t)$, $\pi(l+y,t)$, $v(k,t)$:
преобразования $\tilde\pi$ и $\tilde v$ определяются аналогичным образом. Функции $\tilde\psi$, $\tilde\pi$ периодические по $\theta$ и квазипериодические по $y$, т. е. $\tilde\psi(\theta,y+m,t)=e^{-\mathrm im\cdot\theta}\psi(\theta,y,t)$, $\tilde\pi(\theta,y+m,t)=e^{-\mathrm im\cdot\theta}\pi(\theta,y,t)$, $m\in\mathbb{Z}^d$.
Применим к решениям $Y(\,{\cdot}\,,t)$ вариант преобразования Блоха–Флоке–Зака (или, для краткости, преобразования Зака, см. работу [10]), которое определяется следующим образом:
Заметим, что $\tilde\psi_e$, $\tilde\pi_e$ являются периодическими функциями по $y$ и квазипериодическими по $\theta$, т. е. $\tilde\psi_e(\theta+2\pi e_j,y,t)=e^{\mathrm i2\pi y_j}\tilde\psi_e(\theta,y,t)$, где $e_j\in\mathbb{Z}^d$ – вектор с координатами $e_j^i=\delta_{ij}$ и $\delta_{ij}$ – символ Кронекера. Обозначим через $\mathbb{T}_1^d:=\mathbb{R}^d/\mathbb{Z}^d$ $d$-мерный тор, пусть $\mathcal T^d:=\mathbb{T}_1^d\cup\{0\}$. Далее, $\widetilde Y_e(\theta,t)\equiv \widetilde Y_e(\theta,r,t)$, $r\in\mathcal T^d$. Тогда задача (1.5) сводится к задаче Блоха на торе $\mathbb{T}_1^d$ с параметром $\theta\in\mathrm K^d\equiv [-\pi,\pi]^d$:
Условие V1. Существуют положительные константы $C$ и $\gamma$, такие что $\|V(k)\|\leqslant Ce^{-\gamma|k|}$ при $k\in\mathbb{Z}^d$, где $\|V\|$ обозначает матричную норму.
Условие V2. Имеет место равенство $V^{\mathrm T}(-k)=V(k)$, $k\in\mathbb{Z}^d$, где $V^{\mathrm T}$ обозначает транспонирование матрицы $V$.
Из условий V1 и V2 вытекает, что $\widetilde V(\theta)$ – вещественно-аналитическая эрмитова матричнозначная функция от $\theta\in\mathbb{T}^d$, где $\mathbb{T}^d$ – $d$-мерный тор $\mathbb{R}^d/(2\pi\mathbb{Z})^d$.
Условие V3. Матрица $\widetilde V(\theta)$ положительно определена для каждого $\theta\in\mathbb{T}^d$, т. е. существует константа $\nu_0>0$, такая что $\bar u\cdot\widetilde V(\theta)u\geqslant \nu^2_0|u|^2$, $u\in\mathbb{C}^n$.
Замечание 2.1. Условия V1–V3 выполнены, например, для кристалла со взаимодействием в соседних точках, для которого
и выполнено условие V3 с $\nu_0^2:=\min_l\nu_l^2>0$.
На вещественную векторнозначную функцию взаимодействия $R(x)$ налагаются следующие условия.
Условие R1. Функция $R\in C^{\infty}(\mathbb{R}^d)$ и $|R(x)|\leqslant Ae^{-\varepsilon |x|}$ с некоторыми константами $\varepsilon>0$ и $A<\infty$.
Введем обозначение $\mathrm H_1^s:=\mathrm H^s(\mathbb{T}^d_1)\oplus\mathbb{C}^n$, $s\in\mathbb{R}$, где $\mathrm H^s(\mathbb{T}^d_1)$ – пространство Соболева. Пусть $(\,{\cdot}\,,\,{\cdot}\,)$ обозначает скалярное произведение в $\mathrm H_1^0$,
Здесь $\omega_{\mathrm{min}}^2(\theta)$ – минимальное собственное значение матрицы $\widetilde V(\theta)$, $\omega_{\mathrm{min}}^2(\theta)\,{\geqslant}\,\nu_0^2\,{>0}\,$, где $\nu_0^2$ взято из условия V3. В свою очередь, оценка (2.5) справедлива, если выполнено следующее условие (см. замечание 9.3 в [3]).
2. Условие R2' выполнено для функций $R$, удовлетворяющих R1 с $A\varepsilon^{-d}\ll 1$.
Утверждение 2.1. Пусть выполнены условия V1–V3 и R1, R2. Тогда для любого $Y_0\in\mathcal E$ существует, и притом единственное, решение $Y(t)\in C(\mathbb{R},\mathcal E)$ задачи Коши (1.5). Оператор $W(t)\colon Y_0\mapsto Y(t)$ непрерывен в $\mathcal E$ для любого $t\in\mathbb{R}$,
Это утверждение может быть доказано аналогично предложению 2.7 в [3].
Как было показано в работе [3], из условий R1 и R2 вытекает, что для фиксированного $\theta\in\mathrm K^d$ оператор $\widetilde{\mathcal H}(\theta)$ положительно определен, самосопряжен в $\mathrm H_1^0$ и имеет дискретный спектр. Определим эрмитов положительно определенный оператор
Оператор $\Omega(\theta)$ имеет собственные значения $0<\omega_1(\theta)\leqslant\omega_2(\theta)\leqslant\cdots\leqslant\omega_\sigma(\theta)\leqslant\cdots{}$ и соответствующие спектральные проекторы $\Pi_\sigma(\theta)$ с кратностью $r_\sigma=\operatorname{tr}\Pi_\sigma(\theta)$.
В физике твердого тела собственные значения $\omega_\sigma(\theta)$ называются (энергетическими) зонами Блоха. Функции $\omega_\sigma(\,{\cdot}\,)$ и $\Pi_\sigma(\,{\cdot}\,)$ являются вещественно-аналитическими вне множества точек, где зоны пересекаются. Однако эти функции, вообще говоря, не являются гладкими в точках пересечения. Следующая лемма уточняет свойства функций $\omega_\sigma$ и $\Pi_\sigma$ (см. также раздел XIII.16 в [15]).
Лемма 2.1 (см. лемму 2.8 в [3]). Пусть выполнены условия R1 и R2. Тогда существует замкнутое подмножество $\mathcal C_*\subset \mathrm K^d$ нулевой лебеговой меры, такое что справедливы следующие утверждения.
1. Для любой точки $\Theta\in\mathrm K^d\setminus\mathcal C_*$ и каждого числа $N\in\mathbb{N}$ существует окрестность $\mathcal O(\Theta)\subset \mathrm K^d\setminus\mathcal C_*$, такая что каждая из функций $\omega_\sigma(\theta)$ и $\Pi_\sigma(\theta)$, $\sigma=1,2,\ldots,N$, может быть выбрана вещественно-аналитической в области $\mathcal O(\Theta)$.
2. Собственные значения $\omega_\sigma(\theta)$ имеют постоянную кратность $r_\sigma$ в $\mathcal O(\Theta)$.
3. Для любого $\theta\in\mathcal O(\Theta)$ имеет место спектральное разложение
где $\Pi_\sigma(\theta)$ – ортогональный проектор в $\mathrm H^0_1$.
Предполагается, что система (1.2) удовлетворяет следующим условиям.
Условие R3. Для любого $\Theta\in\mathrm K^d\setminus\mathcal C_*$ и $\sigma=1,2,\ldots{}$ функции $\partial_{\theta_1}\omega_\sigma(\theta)$ не равны тождественно нулю в $\mathcal O(\Theta)$, где $\mathcal O(\Theta)$ введена в лемме 2.1.
Для доказательства сходимости (1.8) достаточно наложить условие R3. Однако для доказательства (1.9) нам потребуется более сильное условие.
Замечание 2.3. 1. В частности, условие R4 выполняется, если $p_{\pm,\sigma\sigma'}^{ij}=0$ при $\sigma\neq\sigma'$, $i=0,1$, а это, в свою очередь, справедливо для гиббсовских мер, рассмотренных в п. 4.1 (см. формулы (4.3) ниже).
2. Вместо условия R4 можно наложить более сильное, но простое условие.
Условие R4'. Если $\sigma\ne \sigma'$, то функции $\omega_\sigma(\theta)\pm\omega_{\sigma'}(\theta)$ в области $\mathcal O(\Theta)$ не равны тождественно ненулевым константам.
2.3. Условия на начальную меру $\mu_0$
Напомним, что $Y_0$ в (1.5) является измеримой случайной функцией, $\mu_0$ обозначает борелевскую вероятностную меру в $\mathcal E$, которая является распределением случайной функции $Y_0$, и $\mathbb{E}$ обозначает математическое ожидание по мере $\mu_0$.
Введем корреляционные функции меры $\mu_0$ следующим образом (см. (1.4)):
где $p=(p_1,\ldots,p_d)$, $p'=(p'_1,\ldots,p'_d)$, $Q_\pm(p,p')$ – корреляционные матрицы некоторых мер $\mu_\pm$, которые имеют нулевое среднее значение и инвариантны относительно сдвигов в $\mathbb{Z}^d$, т. е. выполнено условие (1.7) для почти всех (п.в.) $p,p'\in\mathbb{P}^d$. Следовательно, $Q_\pm(p,p'):=\mathbb{E}_\pm(Y(p)\otimes Y(p'))$ и
Здесь $\mathbb{E}_\pm$ обозначает интеграл по мере $\mu_\pm$, $k=[x]$, $r=x-[x]$, если $p=x\in\mathbb{R}^d$, и $r=0$, если $p=k\in\mathbb{Z}^d$. Аналогично для $k'$ и $r'$.
Пусть $L^2(\mathbb{P}^d):= L^2(\mathbb{R}^d)\oplus [l^2(\mathbb{Z}^d)]^n$ и $\langle\,{\cdot}\,,\,{\cdot}\,\rangle$ обозначает скалярное произведение в $L^2(\mathbb{P}^d)\otimes\mathbb{R}^N$ для различных $N=1,2,\ldots{}$:
Мера $\mu$ называется гауссовой (с нулевым средним значением), если ее характеристический функционал имеет вид $\hat\mu(Z)=e^{-\mathcal Q(Z,Z)/2}$, $Z \in\mathcal D$, где $\mathcal Q$ – вещественная неотрицательная квадратичная форма в $\mathcal D$.
Для доказательства сходимости (1.8) налагаются условия S1–S3. Однако для доказательства утверждения (1.9) (в случае негауссовых начальных мер) нам нужно более сильное условие, чем S2. Чтобы сформулировать его, введем условие перемешивания.
Определение 2.3. Пусть $\sigma (\mathcal A)$ для $\mathcal A\subset\mathbb{P}^d$ обозначает $\sigma$-алгебру в $\mathcal E$, порожденную линейными функционалами $Y\mapsto\langle Y,Z\rangle$, где $Z\in\mathcal D$ с $\operatorname{supp}Z\subset\mathcal A$. Определим коэффициент перемешивания Розенблатта [16] вероятностной меры $\mu_0$ на $\mathcal E$ следующим образом:
Здесь верхняя грань берется по всем множествам $A\in\sigma(\mathcal A)$, $B\in\sigma(\mathcal B)$ и всем открытым выпуклым подмножествам $\mathcal A,\mathcal B\subset\mathbb{P}^d$ с расстоянием $\rho(\mathcal A,\mathcal B)\geqslant r$. Говорят, что мера $\mu_0$ удовлетворяет условию сильного перемешивания Розенблатта, если $\alpha(r)\to 0$ при $r\to\infty$.
для $k\in\mathbb{Z}^d$ и п.в. $x\in\mathbb{R}^d$. Более того, $\mu_0$ удовлетворяет условию сильного перемешивания Розенблатта и $r^{d-1}\alpha^{p}(r)\in L^1(0,+\infty)$ с $p:=\min(\delta/(2+\delta),1/2)$.
Заметим, что из условий S1 и S4 вытекает S2 с $h(s):=10e_\delta^{2/(2+\delta)}\alpha^{\delta/(2+\delta)}(s)$, где величина $e_\delta$ введена в (2.14). Это следует из леммы 17.2.3 в [17].
2.4. Пример начальной меры
Построим пример начальной меры $\mu_0$, которая удовлетворяет условиям S1–S4. Сначала запишем корреляционные матрицы $Q_\pm(p,p')=(Q_\pm^{ij}(p,p'))_{i,j=0,1}$ мер $\mu_\pm$ в виде
для $i,j=0,1$, $p,p'\in\mathbb{P}^d$. Пусть $\mu_\pm$ – гауссовы меры в $\mathcal E$ с корреляционными функциями $Q_\pm^{ij}(p,p')$, $i,j=0,1$. Пусть, кроме того, $Q_\pm^{ij}(p,p')=0$ при $i\neq j$, а при $i=j$
Тогда в силу теоремы Минлоса [18] на $\mathcal E\equiv\mathcal E^{0,\alpha}$, $\alpha<-d/2$, существуют борелевские гауссовы меры $\mu_{\pm}$ с корреляционными матрицами $Q_\pm(p,p')$, так как
Далее введем пару $(Y_{-},Y_{+})$ как единичную случайную функцию на вероятностном пространстве $(\mathcal E\times\mathcal E, \mu_{-}\times\mu_{+})$, где $Y_\pm=(Y_\pm^0,Y_\pm^1)$, $Y_\pm^i=(\psi^i_\pm(x),u^i_\pm(k))$. Тогда $Y_{\pm}$ – гауссовы независимые векторы в $\mathcal E$. Наконец, определим борелевскую вероятностную меру $\mu_0$ как распределение случайной функции
где $\zeta_\pm$ введены в условии S3. Тогда корреляционные функции меры $\mu_0$ имеют вид (2.12). Следовательно, мера $\mu_0$ удовлетворяет условиям S1 и S3.
а $f$ – четная неотрицательная последовательность, такая что $f\in\ell^1$ и $\Delta_{\mathrm L} f(k)\equiv f(k+1)-2f(k)+f(k-1)\geqslant 0$ для любого $k\geqslant 1$. Тогда $\tilde f(\theta)\geqslant 0$ в силу теорем 4.1 и 2.7 в [19] и справедливы условия (2.16). Если дополнительно $|f(k)|\leqslant C(1+|k|)^{-N}$ с $N>d$, то справедливо условие (2.18).
Тогда $\tilde f(\theta)=\frac{1-\cos N_0\theta}{1-\cos\theta}$, $\theta\in\mathbb{T}\equiv \mathbb{R}/(2\pi\mathbb{Z})$, и условие (2.16) выполнено. Более того, в этом случае из (2.19) и (2.20) вытекает, что
Тогда выполнено условие S4 с $\alpha(r)=0$ при $r\geqslant r_0$ и справедливы неравенства (2.18) с $h(r)=0$ при $r\geqslant r_0$.
3. Основные результаты
Определение 3.1. 1. Пусть $\mu_t$, $t\in\mathbb{R}$, – борелевская вероятностная мера в $\mathcal E$, являющаяся распределением решения $Y(t)$ задачи Коши (1.5), $\mu_t(B)=\mu_0(W(-t)B)$ для всех $B\in\mathcal B(\mathcal E)$, $t\in\mathbb{R}$. Здесь и далее $\mathcal B(X)$ обозначает $\sigma$-алгебру борелевских множеств топологического пространства $X$.
2. Корреляционные функции меры $\mu_t$, $t\in\mathbb{R}$, определяются следующим образом:
для всех $Z_1,Z_2\in D_{\mathrm F}\oplus D_{\mathrm L}$. Обозначим через $\mathcal Q_t$, $t\in\mathbb{R}$, квадратичную форму с интегральным ядром $Q_t(p,p')$.
Цель настоящей статьи – доказать слабую сходимость (1.9) мер $\mu_t$ к предельной мере $\mu_\infty$ при $t\to\infty$ в пространствах $\mathcal E^{s,\beta}$ для любых $s<0$ и $\beta<\alpha<-d/2$. Это означает сходимость интегралов
Здесь $(\,{\cdot}\,,\,{\cdot}\,)$ – скалярное произведение в пространстве $[\mathrm H_1^0]^2$ (ср. с формулой (2.4)), $\tilde q_\infty(\theta)=\operatorname{Op}(\tilde q_\infty(\theta,r,r'))$ обозначает интегральный оператор с ядром $\tilde q_\infty(\theta,r,r')$, где
Здесь $\partial_{\theta_1}\equiv\partial/\partial\theta_1$, $\tilde{\mathbf q}^\pm(\theta):=(\tilde q_{+}(\theta)\pm\tilde q_{-}(\theta))/2$, $\Pi_\sigma(\theta)$ – спектральный проектор, введенный в п. 3 леммы 2.1, и $L$ определяется следующим образом:
Здесь и далее $C^*$ обозначает матрицу, эрмитово-сопряженную к матрице $C$.
Замечание 3.1. В силу формул (3.6), (3.7) матрица $\tilde q_\infty(\theta)=\bigl(\tilde q^{ij}_\infty(\theta)\bigr)_{i,j=0}^1$ удовлетворяет “условию равновесия”, т. е.
2. Пусть выполнены условия S1, S2, S4, R1, R2, R3', R4. Тогда для любых $s<0$, $\beta<-d/2$ в пространстве $\mathcal E^{s,\beta}$ имеет место сходимость (1.9). При этом предельная мера $\mu_\infty$ является гауссовой в $\mathcal E$. Кроме того, мера $\mu_\infty$ инвариантна относительно сдвигов в $\mathbb{Z}^d$ и ее характеристический функционал имеет вид
3. Мера $\mu_\infty$ стационарна по времени, т. е. $[W(t)]^*\mu_\infty=\mu_\infty$, $t\in\mathbb{R}$. Более того, группа $W(t)$ удовлетворяет условию перемешивания относительно $\mu_\infty$,
где $Z_\beta$ – нормирующий множитель, функция Гамильтона $\mathrm H(Y)$ определена в (1.1), $T$ – соответствующая абсолютная температура. Введем гиббсовские меры $g_\beta$ как гауссовы борелевские вероятностные меры $g_\beta(dY)=g^0_\beta(dY^0)\times g^1_\beta(dY^1)$ в пространствах $\mathcal E^{s,\alpha}\equiv\mathrm{E}^{s+1,\alpha}\oplus\mathrm{E}^{s,\alpha}$, $s,\alpha<-d/2$, с характеристическими функционалами
Следовательно, $g^1_\beta(dY^1)\equiv g^1_\beta(d\pi\,dv)$ является гауссовой мерой в $\mathrm{E}^{s,\alpha}$, $s,\alpha<-d/2$, с корреляционной матрицей (см. (2.15))
где $\mathrm I$ – единичная матрица. Ниже через $\mathrm I$ мы обозначаем также тождественный оператор в $\mathrm H_1^0$. Так как $Q^{11}_\beta(k+r,k'+r')=:q^{11}_\beta(k-k',r,r')$, в преобразовании Зака
В свою очередь, мера $g^0_\beta(dY^0)\equiv g^0_\beta(d\psi\,du)$ является гауссовой в $\mathrm{E}^{s+1,\alpha}$, $s,\alpha<-d/2$, с корреляционной матрицей $Q_\beta^{00}(p,p')$, где $Q_\beta^{00}(k+r,k'+r')=:q_\beta^{00}(k-k',r,r')$. Применяя преобразование Зака, получаем
Таким образом, мера $g_\beta$ является гауссовой в $\mathcal E^{s,\alpha}$ с корреляционной матрицей $Q_\beta(p,p')=\bigl(Q^{ij}_\beta(p,p')\bigr)_{i,j=0,1}$, где $Q^{ij}_\beta=0$ при $i\neq j$ и
Заметим, что меры $g_\beta$ инвариантны относительно динамики $W(t)$: мы имеем $W(t)^*g_\beta=g_\beta$, $t\in\mathbb{R}$.
Пример 4.1. Пусть $\mu_0=g_0$ есть борелевская вероятностная мера в пространстве $\mathcal E^{s,\alpha}$, $s,\alpha<-d/2$, задающая распределение случайной функции $Y_0$ вида (2.17), где $Y_\pm$ – гауссовы независимые векторы в $\mathcal E^{s,\alpha}$ с распределениями $\mu_{\pm}\equiv g_{\beta_\pm}$, $\beta_\pm=T_\pm^{-1}>0$. Тогда меры $\mu_{\pm}$ имеют корреляционные матрицы $Q_\pm^{ij}(k+r,k'+r'):=q_{\beta_\pm}^{ij}(k-k',r,r')$, определяемые формулами (4.1) и (4.2).
В случае гиббсовских мер $\mu_\pm\equiv g_{\beta_\pm}$ имеем
Следовательно, условие R4 выполнено. Кроме того, в данном случае не требуется условие перемешивания (см. условие S4), поскольку $g_0$ – гауссова мера. Заметим, что $g_0$ удовлетворяет условиям S1 и S3, но не S2.
Обозначим через $g_t$ распределение случайной функции $W(t)Y_0$, $t\in\mathbb{R}$. Существование решений задачи (1.5) в случае начальных данных $Y_0\in\mathcal E^{s,\alpha}$, $s\in\mathbb{R}$, можно доказать аналогично утверждению 2.1, если вместо условия R2 предположить, что (см. оценку (2.5))
Теорема 4.1. Пусть $s,\alpha<-{d}/{2}$. Тогда существует борелевская вероятностная мера $g_\infty$, такая что $g_t\rightharpoondown g_\infty$ при $t\to\infty$ в пространстве $\mathcal E^{s,\alpha}$. Более того, $g_\infty$ – гауссова мера в $\mathcal E^{s,\alpha}$ с нулевым средним значением и с корреляционной матрицей $Q_{\infty}(p,p')$ вида
Доказательство. Воспользуемся рассуждениями, которые использовались при доказательстве п. 1 теоремы 3.1. Однако в случае, когда начальная мера $\mu_0$ удовлетворяет условию S2, справедливы оценки из леммы А.1. Эти оценки позволяют применить теорему Лебега–Римана к осциллирующим интегралам в (А.16) (детали см. в приложении A). Если $\mu_0=g_0$, то условие S2 не выполнено. В этом случае мы применим явные выражения для $\mathcal Q_t$ через корреляционные матрицы гиббсовских мер $g_{\beta_\pm}$.
А именно, сначала запишем квадратичную форму $\mathcal Q_t(Z,Z)$ в виде (А.16). Далее упростим подынтегральное выражение в правой части (А.16), используя (А.17) и (А.18). Для этого заметим, что $p_{\pm,\sigma\sigma'}^{ij}=0$ при $\sigma\neq \sigma'$ в силу (4.3). Следовательно, $M_{\sigma\sigma'}=0$ при $\sigma\neq \sigma'$, и тогда $A^{\pm}_{\sigma\sigma'}=B^{\pm}_{\sigma\sigma'}=0$. При $\sigma=\sigma'$ из (4.3) вытекает, что
Более того, в силу (4.3) и (4.4) $A^{+}_{\sigma\sigma}(\theta)=\Pi_\sigma(\theta) \tilde q_\infty(\theta)\Pi_\sigma(\theta)$ с $\tilde q_\infty(\theta)=\bigl(\tilde q^{ij}_\infty(\theta)\bigr)_{i,j=0}^{1}$. Таким образом, из (А.19) вытекает, что $\mathcal Q_t(Z,Z)=\mathcal Q_\infty(Z,Z)+o(1)$, $t\to\infty$.
Пример 4.2. Построим еще одну гауссову меру $\mu_0$ на $\mathcal E^{s,\alpha}$, $s,\alpha<-d/2$, следующим образом. Предположим, что случайные начальные данные для поля и кристалла $\Psi_0(x)=(\psi^0_0,\psi^1_0)\equiv (\psi_0,\pi_0)$ и $U_0(k)=(u^0_0,u_0^1)\equiv(u_0,v_0)$ независимы. Кроме того, при $t=0$ точки поля с координатами $x_1<-a$ и $x_1>a$ находятся в тепловом равновесии с температурами $T_{-}^{\mathrm F}$ и $T_{+}^{\mathrm F}$. Аналогично в начальный момент времени частицы кристалла с координатами $k_1<-a$ и $k_1>a$ имеют гиббсовские распределения с температурами $T_{-}^{\mathrm L}$ и $T_{+}^{\mathrm L}$. Точнее, введем гиббсовские меры для поля как гауссовы борелевские меры
в пространстве $\mathcal H^{s,\alpha}_{\mathrm F}:=H^{s+1,\alpha}(\mathbb{R}^d)\oplus H^{s,\alpha}(\mathbb{R}^d)$, $s,\alpha<-d/2$, с характеристическими функционалами вида
в пространстве $\mathcal L^\alpha_{\mathrm L}:=L^\alpha(\mathbb{Z}^d)\oplus L^\alpha(\mathbb{Z}^d)$, $\alpha<-d/2$, с характеристическими функционалами вида
где $\mathrm I$ – единичная матрица, $\mathcal E_{\mathrm L}(k)$ – фундаментальное решение оператора $\mathcal V$, т. е. в дискретном преобразовании Фурье
Введем независимые случайные векторы $\Psi_\pm=(\psi^0_\pm,\psi^1_\pm)$ и $U_\pm=(u^0_\pm,u^1_\pm)$ в вероятностных пространствах $(\mathcal H_{\mathrm F}^{s,\alpha},g^{\mathrm F}_{\beta_\pm^{\mathrm F}})$ и $(\mathcal L_{\mathrm L}^\alpha,g^{\mathrm L}_{\beta_\pm^{\mathrm L}})$ соответственно; здесь $\beta_\pm^{\mathrm F}:=(T_\pm^{\mathrm F})^{-1}$ и $\beta_\pm^{\mathrm L}:=(T_\pm^{\mathrm L})^{-1}$. Тогда $Y_\pm:=(\Psi_\pm(x),U_\pm(k))$ – случайные функции в вероятностном пространстве $(\mathcal E^{s,\alpha},g^{\mathrm F}_{\beta_\pm^{\mathrm F}}\times g^{\mathrm L}_{\beta_\pm^{\mathrm L}})$. Наконец, введем борелевскую вероятностную меру $\mu_0$ как распределение случайной функции $Y_0$ вида (2.17). Тогда $\mu_0$ – гауссова мера в пространстве $\mathcal E^{s,\alpha}$, $s,\alpha<-d/2$, ее корреляционная матрица $Q_0$ имеет вид (2.12), где $Q^{ij}_\pm=0$ при $i\neq j$ и
Обозначим через $g_t$ распределение случайной функции $W(t)Y_0$, $t\in\mathbb{R}$. Тогда справедливо следующее утверждение.
Лемма 4.3. Меры $g_t$ слабо сходятся к предельной мере $g_\infty$ при $t\to \infty$ в пространстве $\mathcal E^{s,\alpha}$, $s,\alpha<-d/2$. Более того, $g_\infty$ – гауссова мера в $\mathcal E^{s,\alpha}$ с нулевым средним значением и с корреляционной матрицей $Q_{\infty}(p,p')$ вида (3.3). В преобразовании Зака корреляционные операторы $\tilde q^{ij}_{\infty}(\theta)=\operatorname{Op}(\tilde q^{ij}_{\infty}(\theta,r,r'))$ равны
Для доказательства введем вспомогательные обозначения и докажем необходимые оценки для начальных корреляционных функций.
А.1. Оценки для начальной ковариации
Введем разложение $p=k+r$, где $k\in\mathbb{Z}^d$ и $r\in\mathrm K_1^d\cup \{0\}$. Другими словами, $r=x-[x]\in\mathrm K_1^d$, если $p=x\in\mathbb{R}^d$, и $r=0$, если $p=k\in\mathbb{Z}^d$. Из (1.7) и (2.15) вытекает, что
Квадратичные формы с матричными ядрами $Q_\pm(p,p')$ также непрерывны в $\mathbf L^2$.
2. Пусть $i,j=0,1$, $\sigma,\sigma'\in\mathbb{N}$, функции $p^{ij}_{\pm,\sigma\sigma'}$ определены в (2.9) и $k,p\in\{-1;0\}$ для любых $i,j$ или $k=1$, если $i=0$, и $p=1$, если $j=0$. Тогда операторы $\Omega^{k}(\theta)\widetilde{p}^{ij}_{\pm,\sigma\sigma'}(\theta) \Omega^{p}(\theta)$ удовлетворяют следующей оценке:
3. Квадратичная форма $\mathcal Q_\infty(Z,Z)$ непрерывна в $\mathbf L^2$.
Первое утверждение следствия вытекает из оценок (А.4) с использованием либо теста Шура (см., например, [20]), либо неравенства Юнга. Утверждения 2 и 3 доказываются аналогично (см. работу [3], следствие 3.3).
А.2. Разложение ковариации $Q_t$
Лемма А.2. Функции $\zeta_\pm$, введенные в условии S3, допускают следующие представления (в преобразовании Зака):
где $H(l)$ – функция Хевисайда интервала $[0,+\infty)$. Применяя преобразование Зака, получаем (А.5). Аналогично, если мы положим $\alpha^{\mathrm d}_\pm(k):=\zeta^{\mathrm d}_\pm(k)-\zeta^{\mathrm d}_\pm(k\mp1)$, $k\in\mathbb{Z}$, то получим (А.5) для $\zeta^{\mathrm d}_\pm$.
Применим преобразование Зака к начальной корреляции $Q_0(p,p')$:
Через $\widetilde Q_0(\theta,\theta')$ обозначим интегральный оператор с ядром $\widetilde Q_0(\theta,\theta',r,r')$. Применяя равенство $\widetilde{fg}(\theta)=(2\pi)^{-2d}\tilde f*\tilde g(\theta)$ для $\theta\in\mathbb{R}^{2d}$, представление (2.12) и равенство (А.2), получаем
Здесь $\mathrm I$ – единичная матрица, а матрица $C(\theta)$ определена в (3.7). Применяя преобразование Зака к матрице $Q_t(p,p')$, определенной в (3.1), получаем
С учетом $\widetilde{{\mathcal G}}^{\kern1.5pt\mathrm T}_t(-\theta')=\widetilde{{\mathcal G}}^*_t(\theta')$, используя преобразование (А.7), находим, что при $Z\in\mathcal D$
Заметим, что $\operatorname{supp}\widetilde Z_e\cap (\mathcal C\cup\partial \mathrm K^d)=\varnothing$. Следовательно, для любой точки $\Theta\in\operatorname{supp}\widetilde Z_e$ существует окрестность $\mathcal O(\Theta)\subset \mathrm K^d\setminus(\mathcal C \cup \partial \mathrm K^d) $ со свойствами из леммы 2.1. Таким образом, $\operatorname{supp} \widetilde Z_e\subset\bigcup_{m=1}^M \mathcal O(\Theta_m)$, где $\Theta_m\in\operatorname{supp} \widetilde Z_e$. Следовательно, существует конечное разбиение единицы
Здесь $g_m$ – неотрицательные функции из $C_0^\infty(\mathrm K^d)$ с $\operatorname{supp}g_m\subset\mathcal O(\Theta_m)$, а собственные значения $\omega_\sigma(\theta)$ и проекторы $\Pi_\sigma(\theta)$ являются вещественно-аналитическими функциями от $\theta\in\operatorname{supp}g_m$ для каждого $m$ (мы не нумеруем эти функции индексом $m$, чтобы не усложнять обозначения). Пусть $Z\in\mathcal D_N$ с некоторым $N\in\mathbb{N}$. Используя разложение (2.7) и формулы (А.8)–(А.10), получаем
Следующая лемма вытекает из предложения А.4(i) в [11].
Лемма А.4. Пусть $\chi\in C^1(\mathrm K^d)$ и $\partial_{\theta_1} \omega_\sigma(\theta)\neq 0$ при $\theta\in\operatorname{supp}\chi\subset (\mathrm K^d\setminus\mathcal{C})$. Тогда для $\theta\in\operatorname{supp}\chi$ при $t\to+\infty$ мы имеем
Осциллирующие интегралы в (А.16) с $\omega^\pm_{\sigma\sigma'}(\theta)\not\equiv\text{const}$ стремятся к нулю при $t\to\infty$ по теореме Лебега–Римана, поскольку все подынтегральные выражения в (А.16) суммируемы в силу п. 2 следствия А.1. Более того, тождества $\omega^\pm_{\sigma\sigma'}(\theta)\equiv\text{const}_\pm$ с $\text{const}_\pm\neq 0$ невозможны в силу условия R4' (см. п. 2 замечания 2.3). Если выполнено условие R4, то при $\omega^\pm_{\sigma\sigma'}(\theta)\equiv\text{const}_\pm$ (с $\text{const}_\pm\ne 0$) матричные коэффициенты $A^\mp_{\sigma\sigma'}$ и $B^\pm_{\sigma\sigma'}$ в (А.17) равны нулю. Таким образом, только интегралы с $\omega^{-}_{\sigma\sigma'}(\theta)\equiv 0$ вносят вклад в предел, поскольку из $\omega^{+}_{\sigma\sigma'}(\theta)\equiv 0$ вытекает, что $\omega_\sigma(\theta)\equiv\omega_{\sigma'}(\theta)\equiv 0$, что невозможно в силу условия R2 (а также R3). Поэтому при $t\to\infty$
Наконец, применяя (2.9) и формулу для $C_\sigma$ из (А.13), перепишем матрицу $A^{+}_{\sigma\sigma}$ в виде $A^{+}_{\sigma\sigma}(\theta) =\Pi_\sigma(\theta)\tilde q_\infty(\theta)\Pi_\sigma(\theta)$, где $\tilde q_\infty(\theta)$ определено в (3.6). Утверждение 1 теоремы 3.1 доказано.
Приложение Б. Существование гиббсовских мер
Обозначим через $\mathbb{E}_\beta$ математическое ожидание по мере $g_\beta$. По теореме Минлоса для доказательства леммы 4.1 достаточно проверить, что для всех $s,\alpha<-d/2$
Пусть $Y(p)=(\psi(x),u(k),\pi(x),v(k))$ – случайная функция с распределением $g_\beta$. Обозначим через $C^s_\beta$ корреляционный оператор случайной функции
Тогда $e^s_\beta$ равно следу оператора $C^s_\beta$ в $\mathcal H_0$. Заметим, что $C^s_\beta=\operatorname{Op}(q^s_\beta(0,r,r'))$ является интегральным оператором с ядром $q^s_\beta(0,r,r')$, где
Обозначим через $\tilde q^s_\beta(\theta):=\operatorname{Op}(\tilde q^s_\beta(\theta,r,r'))$ интегральный оператор с ядром $\tilde q^s_\beta(\theta,r,r')$, где $\tilde q^s_\beta(\theta)=\bigl(\tilde q^{s,ij}_\beta(\theta)\bigr)_{i,j=0,1}$. Тогда $\tilde q^{s,ij}_\beta(\theta)=0$ при $i\neq j$ и
где $\tilde q^{ii}_\beta(\theta)$ введены в (4.1) и (4.2). Оператор $\widetilde{\mathbf{\Lambda}}^s(\theta)\colon\mathrm H_1^s\to \mathrm H^0_1$ определяется соотношениями
Оператор $e^{-\mathrm ir\cdot\theta}\colon(\psi(y),u)\to (e^{-\mathrm iy\cdot\theta}\psi(y),u)$ ограничен в $\mathrm H_1^0$. В свою очередь, операторы $\tilde q^{s,ii}_\beta(\theta)$, $i=0,1$, неотрицательны и самосопряжены в $\mathrm H_1^0$, так как
ограничен в $\mathrm H_1^0$ (равномерно по $\theta\in\mathrm K^d$), так как $\|\widetilde{\mathcal H}^{-1}(\theta)Y\|_{\mathrm H_1^1}\leqslant C\|Y\|_{\mathrm H_1^{-1}}$ для всех $Y\in\mathrm H_1^{-1}$ (см. лемму 9.2 в [3]). При этом
Тогда $e^s_\beta\leqslant C\beta^{-1}<\infty$ в силу теоремы 1.6 в [21]. Отсюда вытекает оценка (Б.1).
Конфликт интересов
Автор заявляет, что у нее нет конфликта интересов.
Список литературы
1.
T. V. Dudnikova, A. I. Komech, H. Spohn, “On a two-temperature problem for wave equation”, Markov Process. Related Fields, 8:1 (2002), 43–80, arXiv: math-ph/0508044
2.
T. V. Dudnikova, A. I. Komech, N. J. Mauser, “On two-temperature problem for harmonic crystals”, J. Stat. Phys., 114:3–4 (2004), 1035–1083, arXiv: math-ph/0211017
3.
T. V. Dudnikova, A. I. Komech, “On the convergence to a statistical equilibrium in the crystal coupled to a scalar field”, Russ. J. Math. Phys., 12:3 (2005), 301–325, arXiv: math-ph/0508053
4.
Т. В. Дудникова, А. И. Комеч, “О двухтемпературной проблеме для уравнения Клейна–Гордона”, Теория вероятн. и ее примен., 50:4 (2005), 675–710
5.
T. V. Dudnikova, “Convergence to stationary states and energy current for infinite harmonic crystals”, Russ. J. Math. Phys., 26:4 (2019), 428–453
6.
Т. В. Дудникова, “Сходимость к стационарным неравновесным состояниям для уравнений Клейна–Гордона”, Изв. РАН. Сер. матем., 85:5 (2021), 110–131
7.
Н. Ашкрофт, Н. Мермин, Физика твердого тела, Мир, М., 1979
8.
G. Panati, H. Spohn, S. Teufel, “Effective dynamics for Bloch electrons: Peierls substitution and beyond”, Commun. Math. Phys., 242:3 (2003), 547–578, arXiv: math-ph/0212041
9.
G. Panati, H. Spohn, S. Teufel, “Motion of electrons in adiabatically perturbed periodic structures”, Analysis, Modeling and Simulation of Multiscale Problems, ed. A. Mielke, Springer, Berlin, Heidelberg, 2006, 595–618, arXiv: 0712.4365
10.
J. Zak, “Dynamics of electrons in solids in external fields”, Phys. Rev., 168:3 (1968), 686–695
11.
C. Boldrighini, A. Pellegrinotti, L. Triolo, “Convergence to stationary states for infinite harmonic systems”, J. Stat. Phys., 30:1 (1983), 123–155
12.
H. Spohn, J. L. Lebowitz, “Stationary non-equilibrium states of infinite harmonic systems”, Commun. Math. Phys., 54:2 (1977), 97–120
13.
J.-P. Eckmann, C.-A. Pillet, L. Rey-Bellet, “Non-equilibrium statistical mechanics of anharmonic chains coupled to two heat baths at different temperatures”, Commun. Math. Phys., 201:3 (1999), 657–697
14.
L. Rey-Bellet, L. E. Thomas, “Exponential convergence to non-equilibrium stationary states in classical statistical mechanics”, Commun. Math. Phys., 225:2 (2002), 305–329
15.
М. Рид, Б. Саймон, Методы современной математической физики, т. 4, Анализ операторов, Мир, М., 1982
16.
M. A. Rosenblatt, “A central limit theorem and a strong mixing condition”, Proc. Nat. Acad. Sci. USA, 42:1 (1956), 43–47
17.
И. А. Ибрагимов, Ю. В. Линник, Независимые и стационарно связанные величины, Наука, М., 1965
18.
И. И. Гихман, А. В. Скороход, Теория случайных процессов, т. 1, Наука, М., 1971
19.
Y. Katznelson, An Introduction in Harmonic Analysis, Cambridge Univ. Press, Cambridge, 2004
20.
Н. К. Никольский, Лекции об операторе сдвига, Наука, М., 1980
21.
B. Simon, Trace Ideals and Their Applications, London Mathematical Society Lecture Note Series, 35, Cambridge Univ. Press, Cambridge, 1979
Образец цитирования:
Т. В. Дудникова, “Стабилизация статистических решений при больших временах для гармонической решетки, взаимодействующей с полем Клейна–Гордона”, ТМФ, 218:2 (2024), 280–305; Theoret. and Math. Phys., 218:2 (2024), 241–263