2025-11-14T14:22:18.492353

Functional limit theorems for elephant random walks on general periodic structures

Shibata
This paper investigates functional limit theorems for the Elephant Random Walk (ERW) on general periodic structures, extending the Bertenghi's results on $\mathbb{Z}^d$. Our results reveal new structure-dependent quantities that do not appear in the classical setting $\mathbb{Z}^d$, highlighting how the underlying structure affects the asymptotic behavior of the walk.
academic

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

Основная информация

  • ID статьи: 2511.10347
  • Название: Functional limit theorems for elephant random walks on general periodic structures
  • Автор: Shuhei Shibata (Kyushu University)
  • Классификация: math.PR (Теория вероятностей)
  • Дата публикации: 13 ноября 2025 г. (препринт arXiv)
  • Ссылка на статью: https://arxiv.org/abs/2511.10347

Аннотация

В данной работе исследуются функциональные предельные теоремы для слоновьих случайных блужданий (Elephant Random Walk, ERW) на общих периодических структурах, что расширяет результаты Bertenghi на стандартной целочисленной решётке Zd\mathbb{Z}^d. Исследование выявляет новые структурно-зависимые величины, которые не появляются в классической постановке Zd\mathbb{Z}^d, подчёркивая, как базовая структура влияет на асимптотическое поведение блуждания.

Исследовательский контекст и мотивация

Постановка проблемы

  1. Основной вопрос: Исследование асимптотического поведения случайных блужданий с долгосрочной памятью. Слоновьи случайные блуждания (ERW) были введены Schütz и Trimper в 2004 году для изучения эффектов памяти в одномерном дискретном случайном блуждании, демонстрируя фазовый переход от диффузивного к супердиффузивному поведению.
  2. Значимость проблемы:
    • ERW является важной моделью для исследования эффектов памяти с полной историей
    • Модель демонстрирует существенно различное поведение в разных параметрических областях (диффузивное, критическое, супердиффузивное)
    • Понимание влияния структуры на поведение случайного блуждания имеет теоретическую и прикладную ценность
  3. Ограничения существующих методов:
    • Большинство существующих работ сосредоточены на стандартной решётке Zd\mathbb{Z}^d
    • Bercu и Laulin, Bertenghi и др. установили асимптотическую теорию на Zd\mathbb{Z}^d
    • Отсутствует систематическое исследование более общих периодических структур (треугольная решётка, шестиугольная решётка, кирпичная решётка и т.д.)
  4. Исследовательская мотивация:
    • Обобщение теории ERW на более общие периодические структуры
    • Открытие новых структурно-зависимых величин, раскрывающих влияние базовой геометрической структуры на асимптотическое поведение
    • Установление единого аналитического каркаса через метод урн Пóльи

Основные вклады

  1. Теоретическое расширение: Обобщение функциональных предельных теорем Bertenghi с Zd\mathbb{Z}^d на общие периодические структуры, включая треугольную решётку, шестиугольную решётку и кирпичную решётку.
  2. Открытие новых структурных величин: Выявление и анализ структурно-зависимых величин, не появляющихся в классической постановке Zd\mathbb{Z}^d, таких как матрицы ковариации Σ(U)\Sigma(U) и Σ(W)\Sigma(W), кодирующие информацию о базовой геометрической структуре.
  3. Единый аналитический каркас: Установление единого аналитического каркаса через метод урн Пóльи, применимого к:
    • ERW типа I (одноцветное множество вершин, U=WU=W)
    • ERW типа II (двухцветное множество вершин, UWU \neq W)
  4. Полный анализ фазовой диаграммы: Установление сильного закона больших чисел и функциональных предельных теорем во всех параметрических областях (диффузивной, критической, супердиффузивной) с точным критическим значением pcm=m+12mp_c^m = \frac{m+1}{2m}.
  5. Вычисления для конкретных примеров: Предоставление явных вычислений ключевых величин для нескольких типичных структур (стандартная решётка, треугольная решётка, шестиугольная решётка, кирпичная решётка и т.д.).

Подробное описание методов

Определение задачи

Объект исследования — слоновьи случайные блуждания {Sn}n=0\{S_n\}_{n=0}^{\infty}, определённые на общей периодической структуре Γ\Gamma:

Входные данные:

  • Множества векторов шагов U={u1,,um}RdU = \{u_1, \ldots, u_m\} \subset \mathbb{R}^d и W={w1,,wm}RdW = \{w_1, \ldots, w_{m'}\} \subset \mathbb{R}^d
  • Параметр памяти p(0,1)p \in (0,1)
  • Начальные векторы шагов ui0u_{i_0} и wj0w_{j_0}

Выходные данные:

  • Сильный закон больших чисел: Snn12(uˉ+wˉ)\frac{S_n}{n} \to \frac{1}{2}(\bar{u} + \bar{w}) п.н.
  • Функциональные предельные теоремы: предельные распределения при надлежащей нормализации в различных параметрических областях

Ограничения:

  • Γ\Gamma должна быть решёткой в Rd\mathbb{R}^d
  • Для ERW типа II требуется выполнение условия чередования (двудольная структура)

Архитектура модели

1. Конструкция пространства состояний

Случай U = W (ERW типа I): Γ={i=1mkiui:kiN{0}}\Gamma = \left\{\sum_{i=1}^m k_i u_i : k_i \in \mathbb{N} \cup \{0\}\right\} Одноцветное множество вершин, все вершины структурно эквивалентны.

Случай U ≠ W (ERW типа II): Сначала определяем: Γ0={i=1mkiui+j=1mljwj:i=1mki=j=1mlj,ki,ljN{0}}\Gamma_0 = \left\{\sum_{i=1}^m k_i u_i + \sum_{j=1}^{m'} l_j w_j : \sum_{i=1}^m k_i = \sum_{j=1}^{m'} l_j, k_i, l_j \in \mathbb{N} \cup \{0\}\right\}

Затем: Γ=Γ0(Γ0+U)\Gamma = \Gamma_0 \sqcup (\Gamma_0 + U)

Определяем классы вершин: ZU:=Γ0Z_U := \Gamma_0 и ZW:=Γ0+UZ_W := \Gamma_0 + U, образующие двудольную структуру.

2. Определение ERW

ERW типа II (UWU \neq W):

  • Определение позиции: S2n=i=1n(σi+τi),S2n1=S2(n1)+σnS_{2n} = \sum_{i=1}^n (\sigma_i + \tau_i), \quad S_{2n-1} = S_{2(n-1)} + \sigma_n
  • Механизм выбора шага: В момент времени n+1n+1, при заданной истории σ1,,σn\sigma_1, \ldots, \sigma_n и τ1,,τn\tau_1, \ldots, \tau_n: P(σn+1=σUn)=p,P(σn+1=σ)=1pm1 для σU{σUn}P(\sigma_{n+1} = \sigma_{U_n}) = p, \quad P(\sigma_{n+1} = \sigma) = \frac{1-p}{m-1} \text{ для } \sigma \in U \setminus \{\sigma_{U_n}\} где UnU_n равномерно распределён на {1,2,,n}\{1,2,\ldots,n\}.

ERW типа I (U=WU = W): Sn=i=1nσiS_n = \sum_{i=1}^n \sigma_i Механизм выбора шага аналогичен, но выборка производится из всей истории.

Технические инновации

1. Связь с моделью урн Пóльи

Установление распределительной эквивалентности между ERW и процессом урны:

ERW типа II: {S2n}n=0=d{i=1mXniui+j=1mYnjwj}n=0\{S_{2n}\}_{n=0}^{\infty} \stackrel{d}{=} \left\{\sum_{i=1}^m X_n^i u_i + \sum_{j=1}^{m'} Y_n^j w_j\right\}_{n=0}^{\infty}

ERW типа I: {S2n}n=0=d{i=1mX2niui}n=0\{S_{2n}\}_{n=0}^{\infty} \stackrel{d}{=} \left\{\sum_{i=1}^m X_{2n}^i u_i\right\}_{n=0}^{\infty}

где XnX_n и YnY_n — независимые процессы урн Пóльи.

2. Анализ матрицы замены

Ключевая матрица замены размером m×mm \times m: A=1pm1Jm+mp1m1ImA = \frac{1-p}{m-1}J_m + \frac{mp-1}{m-1}I_m

Собственные значения:

  • λ1=1\lambda_1 = 1 (кратность 1)
  • λ2=mp1m1\lambda_2 = \frac{mp-1}{m-1} (кратность m1m-1)

Критическое значение определяется из λ2/λ1=1/2\lambda_2/\lambda_1 = 1/2: pcm=m+12mp_c^m = \frac{m+1}{2m}

3. Матрицы ковариации, зависящие от структуры

Определение ключевых матриц ковариации размером d×dd \times d: Σ(U)=1mi=1m(uiuˉ)(uiuˉ)\Sigma(U) = \frac{1}{m}\sum_{i=1}^m (u_i - \bar{u})(u_i - \bar{u})^{\top}Σ(W)=1mj=1m(wjwˉ)(wjwˉ)\Sigma(W) = \frac{1}{m'}\sum_{j=1}^{m'} (w_j - \bar{w})(w_j - \bar{w})^{\top}

Эти матрицы кодируют геометрическую конфигурацию векторов шагов. В случае Zd\mathbb{Z}^d они упрощаются до Id/dI_d/d, но в общих структурах могут быть недиагональными, отражая корреляцию диффузии по различным координатным осям.

4. Связь с общей теорией Янсона

Статья искусно применяет общую теорию Янсона 2004 о многотипных ветвящихся процессах и обобщённых урнах Пóльи (особенно Теоремы 3.21, 3.24, 3.31), преобразуя асимптотическое поведение процесса урны в асимптотическое поведение ERW через теорему о непрерывном отображении.

Экспериментальная установка

Теоретическая верификация

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

Вычисления для типичных примеров

В разделе 6 статья предоставляет вычисления ключевых величин для 6 конкретных примеров:

  1. Стандартная решётка Zd\mathbb{Z}^d: m=2dm=2d, pcm=2d+14dp_c^m = \frac{2d+1}{4d}, Σ(U)=Id/d\Sigma(U) = I_d/d
  2. Треугольная решётка: m=6m=6, pcm=7/12p_c^m = 7/12, Σ(U)=I2/2\Sigma(U) = I_2/2
  3. Шестиугольная решётка: m=m=3m=m'=3, pcm=2/3p_c^m = 2/3, Σ(U)=I2/2\Sigma(U) = I_2/2
  4. Кирпичная решётка: m=m=3m=m'=3, pcm=2/3p_c^m = 2/3, Σ(U)=29(3001)\Sigma(U) = \frac{2}{9}\begin{pmatrix}3 & 0\\0 & 1\end{pmatrix}
  5. Смешанная структура 1: U={±u1,±u2}U = \{\pm u_1, \pm u_2\}, W={±e1,±e2}W = \{\pm e_1, \pm e_2\}, m=m=4m=m'=4
  6. Смешанная структура 2: U={±e1,±e2,(1,2)}U = \{\pm e_1, \pm e_2, (1,2)^{\top}\}, W={±e1,±e2}W = \{\pm e_1, \pm e_2\}, m=5,m=4m=5, m'=4

Основные результаты

1. Сильный закон больших чисел (Теорема 5.1)

Для всех p(0,1)p \in (0,1): Snn12(uˉ+wˉ)п.н. при n\frac{S_n}{n} \to \frac{1}{2}(\bar{u} + \bar{w}) \quad \text{п.н. при } n \to \infty

Ключевые особенности:

  • Предел не зависит от числа рёбер m,mm, m'
  • Зависит только от среднего значения векторов шагов
  • Справедлив для ERW типа I и типа II

2. Диффузивная область (Теорема 5.3)

Условие: 0<p<pcmpcm0 < p < p_c^m \leq p_c^{m'}

Результат: {S2ntnt(uˉ+wˉ)n}t0{Wt}t0\left\{\frac{S_{\lfloor 2nt \rfloor} - nt(\bar{u} + \bar{w})}{\sqrt{n}}\right\}_{t \geq 0} \Rightarrow \{W_t\}_{t \geq 0}

где {Wt}t0\{W_t\}_{t \geq 0} — центрированный Rd\mathbb{R}^d-значный непрерывный гауссовский процесс со структурой ковариации: E[WsWt]=Cas(ts)aΣ(U)+Cas(ts)aΣ(W)\mathbb{E}[W_s W_t^{\top}] = C_a s\left(\frac{t}{s}\right)^a \Sigma(U) + C_{a'} s\left(\frac{t}{s}\right)^{a'} \Sigma(W)

где a=mp1m1a = \frac{mp-1}{m-1}, Ca=112aC_a = \frac{1}{1-2a}.

Инновационные аспекты:

  • Матрицы ковариации Σ(U),Σ(W)\Sigma(U), \Sigma(W) воплощают структурную зависимость
  • В случае Zd\mathbb{Z}^d вырождается в результат Bertenghi
  • Когда Σ(U),Σ(W)\Sigma(U), \Sigma(W) недиагональны, диффузия по различным координатным осям коррелирована

3. Критическая область (Теорема 5.5)

Условие: 0<pcm=p=pcm0 < p_c^m = p = p_c^{m'}

Результат: {S2ntnt(uˉ+wˉ)nt/2logn}t0{Wt}t0\left\{\frac{S_{\lfloor 2nt \rfloor} - nt(\bar{u} + \bar{w})}{n^{t/2}\sqrt{\log n}}\right\}_{t \geq 0} \Rightarrow \{W_t\}_{t \geq 0}

Структура ковариации: E[WsWt]=sΣ(U,W)\mathbb{E}[W_s W_t^{\top}] = s\Sigma(U, W)

где Σ(U,W)=Σ(U)+Σ(W)\Sigma(U, W) = \Sigma(U) + \Sigma(W).

Особенности:

  • Требуется дополнительный множитель нормализации logn\sqrt{\log n}
  • Предельный процесс может быть представлен как Wt=Σ(U,W)1/2BtW_t = \Sigma(U,W)^{1/2}B_t, где BtB_t — стандартное броуновское движение
  • В случае Zd\mathbb{Z}^d процесс {Wt/d}\{W_t/\sqrt{d}\} является стандартным броуновским движением

4. Супердиффузивная область (Теорема 5.8)

Условие: 0<pcm=pcm<p0 < p_c^m = p_c^{m'} < p

Результат для ERW типа II: {S2ntnt(uˉ+wˉ)na}t0{taL}t0\left\{\frac{S_{\lfloor 2nt \rfloor} - nt(\bar{u} + \bar{w})}{n^a}\right\}_{t \geq 0} \Rightarrow \{t^a L\}_{t \geq 0}

Результат для ERW типа I: {Sntntuˉna}t0{taL~}t0\left\{\frac{S_{\lfloor nt \rfloor} - nt\bar{u}}{n^a}\right\}_{t \geq 0} \Rightarrow \{t^a \tilde{L}\}_{t \geq 0}

где L,L~L, \tilde{L} — ненулевые Rd\mathbb{R}^d-значные случайные векторы.

Важные наблюдения:

  • ERW типа I и типа II имеют различные предельные распределения в супердиффузивной области
  • Предельное распределение зависит от выбора начального шага
  • При предположениях uˉ=wˉ=0\bar{u} = \bar{w} = 0 и равномерном выборе начального шага: E[LL]=1(2a1)Γ(2a)Σ(U,W)\mathbb{E}[LL^{\top}] = \frac{1}{(2a-1)\Gamma(2a)}\Sigma(U,W)

5. Смешанные параметрические области (Замечания 5.7, 5.9)

Статья также обсуждает случаи pcmpcmp_c^m \neq p_c^{m'}:

  • Когда 0<pcm=p<pcm0 < p_c^m = p < p_c^{m'}, процесс XX доминирует, предел зависит только от Σ(U)\Sigma(U)
  • Когда 0<pcm<pcm<p0 < p_c^m < p_c^{m'} < p, доминирующий член определяется в зависимости от соотношения aa и aa'

Связанные работы

Исследования одномерного ERW

  • Schütz и Trimper 2004: Введение модели ERW, открытие фазового перехода при p=3/4p=3/4
  • Baur и Bertoin 2016: Установление связи ERW с урнами Пóльи
  • Bercu 2017: Анализ ERW методом мартингалов
  • Coletti и др. 2017: Центральная предельная теорема
  • Kubota и Takei 2019: Гауссовские флуктуации в супердиффузивной области

Исследования многомерного ERW

  • Bercu и Laulin 2019: Исследование асимптотического поведения MERW методом мартингалов
  • Bertenghi 2022: Установление функциональных предельных теорем на Zd\mathbb{Z}^d (работа, непосредственно обобщаемая в данной статье)
  • González-Navarrete 2020: Многомерные блуждания со случайным смещением
  • Chen и Laulin 2023: Многомерное ERW с плавным забыванием
  • Curien и Laulin 2024: Возвратность плоского ERW
  • Qin 2025: Возвратность и транзиентность многомерного ERW

Проблемы столкновения

  • Roy, Takei и Tanemura 2024: Проблема столкновения двух ERW на Z\mathbb{Z}
  • Shibata и Shirai 2025: Проблемы столкновения и асимптотика расстояния при различных параметрах памяти

Теория урн Пóльи

  • Janson 2004: Функциональные предельные теоремы для многотипных ветвящихся процессов и обобщённых урн Пóльи (основной теоретический инструмент статьи)
  • Athreya и Karlin 1968: Вложение модели урны в непрерывный ветвящийся процесс
  • Chauvin и др. 2011: Предельные распределения для больших урн Пóльи

Преимущества данной работы

  1. Общность: Не ограничивается Zd\mathbb{Z}^d, охватывает широкий класс периодических структур
  2. Структурная зависимость: Раскрывает новые структурно-зависимые величины
  3. Единый каркас: Единообразная обработка ERW типа I и типа II
  4. Полнота: Охватывает полную теорию для всех параметрических областей

Заключение и обсуждение

Основные выводы

  1. Успешное теоретическое расширение: Успешное обобщение результатов Bertenghi с Zd\mathbb{Z}^d на общие периодические структуры, демонстрирующее универсальность метода урн Пóльи.
  2. Количественное описание влияния структуры: Точное количественное описание влияния базовой геометрической структуры на асимптотическое поведение ERW посредством матриц ковариации Σ(U)\Sigma(U) и Σ(W)\Sigma(W).
  3. Сохранение фазовых переходов: Критическое значение pcm=m+12mp_c^m = \frac{m+1}{2m} определяется собственными значениями матрицы замены, явление фазового перехода сохраняется в общих структурах.
  4. Различие между типом I и типом II: За исключением супердиффузивной области, оба типа ERW имеют одинаковые предельные распределения; в супердиффузивной области влияние начальных условий приводит к различным пределам.

Ограничения

  1. Ограничения структуры:
    • Рассматриваются только периодические структуры, представимые в виде решёток
    • Для многодольных графов (l3l \geq 3), таких как решётка кагоме, требуются различные аналитические методы
    • Исключены тривиальные случаи при m=1m=1
  2. Начальные условия:
    • Предполагается начало из начала координат
    • Первые два шага детерминированы
    • Результаты в супердиффузивной области зависят от распределения начальных шагов
  3. Диапазон параметров:
    • Рассматривается только p(0,1)p \in (0,1)
    • p=1p=1 соответствует тривиальному случаю (детерминированное блуждание)
    • Случай p=0p=0 не обсуждается
  4. Детали распределения:
    • Распределение предельной случайной величины LL в супердиффузивной области не полностью описано
    • Даны только выражения для первых двух моментов

Направления будущих исследований

  1. Расширение на более общие структуры:
    • Непериодические структуры
    • Многодольные графы (l3l \geq 3), такие как решётка кагоме
    • ERW в случайной среде
  2. Тонкий анализ предельных распределений:
    • Полное распределение LL в супердиффузивной области
    • Применение метода неподвижной точки Guérin и др. 2023, 2025
  3. Проблемы столкновения:
    • Проблема столкновения двух ERW на общих структурах
    • Случаи с различными параметрами памяти
  4. Другие свойства:
    • Полная характеризация возвратности и транзиентности
    • Принципы больших отклонений
    • Исследование локального времени

Глубокая оценка

Достоинства

  1. Математическая строгость:
    • Доказательства полные и строгие, адекватно используется общая теория Янсона
    • Классификация различных параметрических областей ясна и полна
    • Технические детали обработаны надлежащим образом (например, вычисления ковариации)
  2. Теоретический вклад:
    • Первое систематическое исследование ERW на общих периодических структурах
    • Открытие новых структурно-зависимых величин Σ(U),Σ(W)\Sigma(U), \Sigma(W)
    • Установление единого каркаса для ERW типа I и типа II
  3. Ясность изложения:
    • Рациональная организация статьи, переход от простого к сложному
    • Наглядные диаграммы (треугольная решётка, шестиугольная решётка, кирпичная решётка)
    • Вычисления конкретных примеров в разделе 6 очень полезны
  4. Полнота:
    • Охватывает все параметрические области (диффузивную, критическую, супердиффузивную)
    • Обсуждаются различные смешанные параметрические случаи
    • Связь с существующей литературой ясна

Недостатки

  1. Ограниченность примеров:
    • Хотя теория общая, приведено только 6 конкретных примеров
    • Отсутствует обсуждение некоторых интересных структур (например, решётка кагоме)
    • Можно было бы добавить больше примеров нестандартных структур
  2. Физическая интуиция:
    • Недостаточно физической или геометрической интуиции для структурной зависимости
    • Смысл недиагональных элементов Σ(U),Σ(W)\Sigma(U), \Sigma(W) не полностью объяснён
    • Механизм, посредством которого различные структуры приводят к различному поведению, недостаточно обсуждён
  3. Супердиффузивная область:
    • Характеризация предельного распределения LL неполна
    • Даны только вторые моменты, высшие моменты или полное распределение неизвестны
    • Глубокие причины различия между типом I и типом II недостаточно выяснены
  4. Обсуждение приложений:
    • Отсутствует обсуждение практических сценариев применения
    • Не указано, какие выводы эти результаты дают для понимания реальных систем

Влияние

  1. Вклад в область:
    • Открывает новое направление в исследовании ERW (общие структуры)
    • Предоставляет методологию, применимую к другим случайным блужданиям с памятью
    • Обогащает приложения теории урн Пóльи
  2. Практическая ценность:
    • Теория сильна, но может служить теоретической основой для изучения диффузии в материалах, случайных процессов на сетях и т.д.
    • Формулы в разделе 6 можно непосредственно использовать для конкретных вычислений
  3. Воспроизводимость:
    • Доказательства полные и могут быть проверены
    • Формулы для вычисления ключевых величин явны
    • Теоретические результаты применимы к новым структурам

Сценарии применения

  1. Математические исследования:
    • Предельные теоремы в теории вероятностей
    • Теория случайных процессов
    • Комбинаторные случайные структуры
  2. Физические системы:
    • Процессы диффузии на кристаллических решётках
    • Явления переноса с памятью
    • Исследование фазовых переходов
  3. Сетевая наука:
    • Случайные блуждания на графах
    • Модели распространения информации
    • Алгоритмы исследования сетей
  4. Статистическая физика:
    • Немарковские процессы
    • Системы с долгосрочными корреляциями
    • Аномальная диффузия

Избранные ссылки

  • 2 Baur & Bertoin (2016): Elephant random walks and their connection to Pólya-type urns. Physical Review E.
  • 5 Bercu & Laulin (2019): On the multi-dimensional elephant random walk. J. Stat. Phys.
  • 7 Bertenghi (2022): Functional limit theorems for the multi-dimensional elephant random walk. Stoch. Models.
  • 17 Janson (2004): Functional limit theorems for multitype branching processes and generalized Pólya urns. Stochastic Process. Appl. (основной теоретический инструмент)
  • 23 Schütz & Trimper (2004): Elephants can always remember: Exact long-range memory effects in a non-Markovian random walk. Physical Review E. (происхождение ERW)

Общая оценка: Это высококачественная теоретическая статья по теории вероятностей, успешно обобщающая теорию ERW на общие периодические структуры и раскрывающая глубокое влияние структуры на асимптотическое поведение. Математически строга и полна, с существенным теоретическим вкладом, закладывающим основу для дальнейших исследований в этой области. Основные недостатки заключаются в недостаточном обсуждении физической интуиции и сценариев применения результатов, а также в неполной характеризации предельного распределения в супердиффузивной области.