Научная статья на тему 'Непрерывность вероятности перколяции бернуллиевских случайных полей на однородных древесных графах'

Непрерывность вероятности перколяции бернуллиевских случайных полей на однородных древесных графах Текст научной статьи по специальности «Теория меры, представления булевых алгебр, динамические системы»

CC BY
472
50
Поделиться
Ключевые слова
вероятность перколяции / однородное дерево кэли / бернуллиевское случайное поле
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Текст научной работы на тему «Непрерывность вероятности перколяции бернуллиевских случайных полей на однородных древесных графах»

УДК 517.987

НЕПРЕРЫВНОСТЬ ВЕРОЯТНОСТИ ПЕРКОЛЯЦИИ БЕРНУЛЛИЕВСКИХ СЛУЧАЙНЫХ ПОЛЕЙ НА ОДНОРОДНЫХ ДРЕВЕСНЫХ ГРАФАХ Е.С. Антонова, Ю.П. Вирченко

Белгородский государственный университет,

ул. Победы, 85, Белгород, 308015, Россия, e-mail: antonova_e_s@mail.ru

Аннотация. В работе рассматривается «задача узлов» дискретной теории перколяции на графах типа однородного дерева Кэли. Изучается вероятность перколяции P(c) однородного бернуллиевского случайного поля из корневой вершины. Доказано, что при любом индексе s € {2, 3, 4,...} ветвления дерева функция P(с) непрерывна.

Ключевые слова: вероятность перколяции, однородное дерево Кэли, бернуллиевское случайное поле.

1. Введение. Теория перколяции изучает задачи, связанные с понятием связности случайных множеств (см., например, [1]-[3]). Она подразделяется на два самостоятельных направления - дискретная и непрерывная теории перколяции. Несмотря на важность теории для приложений, в особенности в физике твёрдого тела, ввиду чрезвычайной сложности возникающих в ней задач, имеется не так уж много точно установленных математических фактов. Это касается как результатов качественного характера, так и количественных оценок основных перколяционных характеристик с гарантированной точностью. Основные математические достижения в теории перколяции относятся к дискретной её части, которая имеет дело со случайными множествами на бесконечных графах, которые при заданной структуре смежности исходного графа порождают на нём случайные подграфы.

В настоящей работе изучается самая простая модель дискретной теории перколяции, в которой случайное множество индуцируется однородным бернуллиевским случайным полем и бесконечный граф является однородным деревом с постоянным порядком ветвления. Эта модель является исторически первой перколяционной моделью, возникшей ещё до появления термина перколяция и формулировки основных понятий теории, а именно, своим появлением она обязана теории ветвящихся марковских случайных процессов с дискретным временем (см., например, [4]), а именно, теории процесса Гальтона-Ватсона. Несмотря на это, до сих пор имеется ряд связанных с рассматриваемой моделью невыясненных принципиальных вопросов. Статья носит, скорее, методический характер и служит введением в проблематику, относящуюся к указанной перколяцион-ной модели.

2. Перколяция на однородном древесном графе. Будем рассматривать бесконечные неориентированные графы (V, Ф), не содержащие петель, где V - множество вершин графа и Ф С V х V - симметричное подмножество, не содержащее точек диагонали {(х, х); X £ V}, которое определяет отношение смежности вершин. В дальнейшем

вершины графа мы будем обозначать жирными строчными буквами латинского алфавита. Если пара {х, у) Є Ф, то они называются смежными и отношение их смежности мы будем обозначать х^у.

Всякую последовательность {х, Хі, х2,...) конечную или бесконечную, для которой имеет место ХкірХк+і, к = 0,1, 2,..., называется путём на графе {V, Ф). При этом первая компонента последовательности называется началом пути. Пути на графе мы будем обозначать греческой буквой 7. Число компонент этой последовательности 7 = {х, х1, х2,...) называется длиной пути и мы будем обозначать эту величину посредством І71. В частности, если последовательность бесконечна, то |71 = то. Путь называется несамопересекающимся, если для любого к Є М+ имеет место хк = хк+і, І Є N хо = х.

Определение 1. Граф {V, Ф) называется однородным деревом с порядком ветвления в Є N\ {1}, если множество его вершин V состоит из вершины 0, которая называется корневой, и множества вершин г, которое координатизируется множеством конечных последовательностей {іі, ...,]п), где і к =1 ^ в, к =1 ^ п для всех п Є N. Таким образом, V = {0}и{£ : г = {іі,...,іп)}, п Є N. Множество Ф смежности этого графа определяется следующими отношениями смежности: 0^{і) и

{Іі,--,Зп)^{Зі,---,Зп,І) , І = 1 ^ в

для любой вершины, определяемой последовательностью {іі, ...,іп) с ік = 1^в, к = 1^п, п Є N.

Определение 2. Случайное поле {с(х);хЄ V) на графе {V, Ф) называется однородным бернуллиевским, если все случайные величины с(х), хЄ V статистически независимы в совокупности, принимают только два значения {0,1} и одинаково распределены.

Бернулливское однородное поле полностью определяется значением вероятности c = Pr{c(x)= 1}, которую мы будем называть концентрацией. Бернулливское случайное поле на графе индуцирует случайное множество с пространством погружения V, случайные реализации которого определяются формулой C = { х: с(х) = 1}. В дальнейшем эти реализации мы называем конфигурациями. Каждая случайная конфигурация определяет на графе (V, Ф) случайный подграф, множеством вершин которого служит C, а отношение смежности индуцируется отношением смежности Ф исходного графа. При этом пути на исходном графе устанавливают отношение связанности на любом случайном подграфе. Легко понять, что это отношение связности является отношением эквивалентности. Поэтому каждая случайная конфигурация C разбивается дизъюнктивным образом на связные компоненты, согласно введеному понятию эквивалентности.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Обозначим посредством Со связную компоненту конфигурации, которая содержит корневую вершину 0. Если 0 € С, то положим С0 = 0. Вероятность Pr{|C0| = то} = P(c) называется вероятностью перколяции из корневой вершины (здесь и далее символ | • | обозначает число элементов множества). Если она положительна, то говорят, что из неё имеется перколяция. Очевидно, что данное определение вероятности перколяции

эквивалентно следующему

Р(с) = Рг{3(7 С Со : |71 = то)} .

Далее, посредством 7 мы обозначаем пути на С, начинающиеся в корневой вершине. Если путь начинается в другой вершине, то этот факт мы отмечаем соответствующими индексами. В том случае, когда длина пути 7 не меньше, чем п, то будем говорить, что имеется перколяция на расстояние п от корневой вершины. Обозначим вероятность этого случайного события посредством Рп(с). Таким образом,

Р„(с) = Рг{3(7 С Со : |71 ^ п)} .

При этом положим

Ро(с) = Рг{Со = 0} = с.

Замечание 1. Имеется связь между перколяцией на однородных деревьях и марковскими ветвящимися случайными процессами с дискретным временем [4]. Наличие перколяции с точки зрения теории таких процессов означает, что для них реализуется т.н. надкритический режим.

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

Обозначим посредством Рп(с) вероятность перколяции из корневой вершины на расстояние п, то есть вероятность того, что на случайной конфигурации С = {ж: с(х) = 1} существует путь 7 из корневой вершины, имеющий длину не менее, чем п. Она является функцией от концентрации Рг{с(х) = 1} = с. Таким образом, вероятность Рп(с) определяется формулой

Р„(с) = Рг{3(7 С Со : |71 ^ п)} .

Непосредственно из этого определения, замечаем, что верна

Теорема 1. Последовательность функций (Рп(с); п £ М+) такова, что при каждом с £ [0,1] числовая последовательность (Рп(с); п £ М) является невозрастающей и существует предел Иш Рп(с) .

□ Очевидно, что имеет место включение случайных событий

{3(7 С Со : ы ^ п + 1)} С {3(7 С Со : |71 ^ п)} .

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Поэтому для вероятностей этих событий выполняется неравенство, утверждаемое в формулировке теоремы,

Рп+1 (с) = Рг{3(7 С Со : |71 ^ п + 1)} ^ Рг{3(7 С Со : |71 ^ п)} = Рга(с).

Так как Рп(с) ^ 0 для любого с £ [0,1], то для всех значений с невозрастающие числовые последовательности (Рп(с); п £ М) ограничены снизу, что приводит к утверждению теоремы. Ж

Следствием этой теоремы является связь между вероятностью перколяции Р(с) и вероятностями Рп(с), п € М+.

Теорема 2. Имеет место

Р(с) = Пт Рга(с). (1)

n

n—

□ Формула (1) доказывается использованием свойства непрерывности вероятности относительно теоретико-множественного предельного перехода,

lim Pr {3(y С C0 : |y| ^ n)} = Pr f lim {3(7 С C0 : |y| ^ n)}) =

n—^ L J \n—<^ /

= Pr{3(y С Co : |yI = <^)} = P(c) ■ В

Определение З. Число

c* = inf{c : P(c) > О}

называется порогом перколяции однородного бернуллиевского поля на графе (V, Ф).

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Введём вероятности Qn(c) = І — Pn(c) того, что на конфигурации C = {x: c(x) = І} все пути y из корневой вершины графа имеют длину менее n,

Qn(c)=Pr{V(Y С Co : |y| < n)} ■ (2)

Теорема 3. В каждой точке с Є [0,1] последовательность функций (Рга(с); п Є М+) не убывает и имеет предел, равный Р(с) = 1 — Р(с).

□ Первая часть утверждения следует из связи Рга(с) = 1 — Рп(с) значений функций

и Рп и из монотонного невозрастания последовательности (Рп(с); п Є М+) при каждом с Є [0,1]. Монотонно неубывающая числовая последовательность (Рга(с); п Є М+), с Є [0,1] всегда имеет предел. На основании (1) это предельное значение Р(с) равно

Р(с) = Ііт р„(с) = Ііт (1 — Рп(с)) = 1 — ІІШ Рп(с) = 1 — Р(с). ■ (3)

Свойство невозрастания каждой из функций Рга(с), п Є М+ приводит к тому, что справедлива

Теорема 4. Функция Р(с) не возрастающая (Р(с) не убывающая) на [0,1].

□ Сформулированное утверждение следует из того, что все функции Рга(с), п Є М+ невозрастающие по с, и поточечный предел последовательности (Рга(с); п Є М+) является невозрастающей функцией. ■

Особенностью перколяции на древесных графах является то, что совокупность вероятностей Рга(с), п Є М+ или совокупность вероятностей Рп(с), п Є М+ удовлетворяют соответствующему разностному уравнению, на основе которого они могут быть вычислены, последовательно переходя от значения п к (п +1).

Теорема 5. Если однородное дерево имеет порядок ветвления равный в € {2, 3, 4,...}, вероятности Рга(с), п € М+ связаны уравнением

Рга+1(с) = 1 - с + С РП(с) . (4)

□ Воспользуемся дизъюнктивным разложением

{У(7 С Со : |71 < п + 1)} = {Со = 0} и {Со = 0, У(т С Со : |71 < п + 1)} .

Второе событие в этом разложении запишем в следующей форме {Со = 0,У(7 С Со : |71 < п + 1)} =

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

= {Со = 0,У(71 С Со : Ы < п),У(7з С Со : Ы < п)} =

= П{Со = 0^(Т' С Со : |7^-1 < п)} ,

'=1

где 7^ - пути, начинающиеся соответственно в вершинах (3), 3 = 1 ^ в и несодержащие корневую вершину. События

{Со = 0,У(т,- С Со : ^1 < п)} , з = 1 ^ в

условно независимы относительно условия {Со = 0}. Это следует из того, что компонента Со случайной конфигурации представляется в виде дизъюнктивного объединения путей Ту, з = 1 ^ в и одноточечного множества {0}. Тогда

Рг{Со = 0,У(т С Со : |т| < п + 1)} =

ПРг{Со = 0.У(7, С Со : Ь| < п)|Со = 0} ) Рг{Со = 0} .

4.7=1 /

Так как Рг{Со = 0} = с и все условные вероятности Рг{Со = 0, У(т? С Со : |т? | <

п)|Со = 0}, 3 = 1 ^ в совпадают, ввиду однородности дерева, и равны безусловной

вероятности Рг{У(т С Со : |т| < п)}, то мы получаем

Рг{Со = 0,У(т С Со : |т| < п + 1)} = РП(с) .

Учитывая, кроме того, что Рг{Со = 0} = 1 — с, находим

Рг{Со = 0, У(т С Со : Ы < п +1)} = 1 — с + сРП(с). ■

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Следствие 1. Вероятности Рп(с), п € М+ подчинены уравнению

Рп+1(с) = с (1 — (1 — Рга(с))5) •

На основе уравнения (4) вероятности Фп(с) вычисляются однозначным образом на при использования "начального" условия Фо(с) = 1 — с.

Теорема 6. Все функции Фп(с), п € М+ - убывающие.

□ Из уравнения (4) имеем

фП+1(с) = фП(с) — 1 + вс 1(сШс) (5)

и, кроме того, Фо(с) = —1 < 0. Так как ФП(с) ^ 1, то индукцией по п € М+, использующей соотношение (5) в качестве индукционного шага получаем утверждение теоремы. ■

Из рекуррентного соотношения (4) сразу следует Теорема 7. Вероятность Ф(с) является решением уравнения

Ф(с) = 1 — с + с Ф5(с) , (6)

которое подчинено неравенствам 0 ^ Ф(с) ^ 1. Точно также вероятность Р(с) является решением уравнения

Р(с) = с(1 — (1 — Р(с))5) , 0 ^ Р(с) ^ 1.

□ Доказательство получается переходом к пределу п то в (4) с использованием

(3). Уравнение для вероятности Р(с) получается из (6) заменой ф(с) = 1 — Р(с). ■

4. Основная теорема. На основании следующего утверждения будет дан ответ на вопрос о поведении вероятности Ф(с), как функции от концентрации с.

Лемма. Пусть значения каждой из совокупности функций {Л-Д-, с) : [0,1] ^ К; с € [0,1], в € {2, 3,...}} определяются формулой

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Л5(г, с) = сг5 — г + 1 — с, г € [0,1]. (7)

Тогда уравнение

Л8(г,с) = 0 (8)

имеет:

1) единственный корень г =1 при условии св ^ 1;

2) два корня {1, г*(с)}, г*(с) < 1 при условии св > 1.

□ Функция Л5(г,с) выпукла по г, и поэтому уравнение (2) имеет не более двух корней. При любых значениях с € [0,1] это уравнение имеет корень г =1, Л5(1, с) = 0.

Так как Л5(0, с) = 1 — с > 0, то для того чтобы существовал на интервале (0,1) единственный корень, необходимо и достаточно чтобы минимум функции ЛД-, с) достигался на (0, 1) и был отрицателен.

Точкой минимума функции Л5(^с) является = (вс)-1/(5-1). Из этого выражения видно, что € (0,1) тогда и только тогда, когда св > 1. Таким образом, при св ^ 1

функция Л8(-,с) не имеет минимума на (0,1), и поэтому Л5(г,с) > 0, г € (0,1), что доказывает утверждение 1) леммы.

Пусть св > 1. Тогда

Поэтому из выражения (9) следует, что Л5(гт, с) < 0. Это означает, что в рассматриваемом случае существует единственный корень г*(с) уравнения (8) на интервале (0,1).

Так как Л5(0,с) > 0, то из утверждения 2) леммы следует, что необходимым и достаточным условием для выполнимости неравенства Л5(г,с) > 0 является условие г < г*(с). ■

Следующая теорема выявляет строение функции Ф(с), которая является решением уравнения (6). При её доказательстве, казалось бы, естественно было использовать теорему о неявной функции. Однако, оказывается, гораздо проще воспользоваться ограничением, которое возникает вследствие утверждения доказанной леммы, благодаря предельному соотношению (3) и разностному уравнением (4).

Теорема 8. Функция Ф(с) (Р(с)) определяется следующими формулами:

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

1) ф(с) = 1 (Р(с) = 0) при с ^ в-1;

2) ф(с) = г*(с) (Р(с) = 1 — г*(с)) при с > в-1.

□ Утверждение 1) тривиально, так как решение уравнения (6) единственно и равно 1 при с ^ в-1.

Ввиду того, что при каждом с € [0,1] выполняется Фп+1(с) ^ Фп(с), используя (4), получаем неравенство

то есть Л5(фп(с),с) > 0, что эквивалентно неравенству (см. емма 3) Фп(с) ^ г*(с) < 1. Переходя к пределу п ^ то в этом неравенстве, на основании (2), имеем Ф(с) ^ г* (с) <

только либо z*(c), либо 1, то Q(c) = z*(c). ■

Заметим, что при доказательстве теоремы не было использовано a priori свойство непрерывности функции Q(c), которое само нуждается в отдельном доказательстве (см. Следствие 3).

Следствие 2. Точка с = s-1 является порогом перколяции случайного однородного бернуллиевского поля на однородном древесном графе с порядком ветвления s.

Следствие 3. Функция Q(c) вещественно-аналитична на (1/s, 1) и непрерывна в точке с = 1/s.

Определим значения функции д : К ^ К формулой д(у) = у(1 — у5-1). Эта функция достигает максимума на [0,1] в точке в-1/(5-1). Следовательно, полагая у = с1/(5-1), имеем д (с1/(5-1)) < д (в-1/(,5-1)) при с = 1/в или, что эквивалентно,

с Qn — Qn(c) + 1 — с ^ 0 ,

1. С другой стороны, так как при каждом с € [0,1] функция Ф(с) может быть равна

□ Функция Ф(с) является вещественно-аналитической на (1/в, 1), так как она в этом случае определяется как ветвь аналитической функции, являющейся решением уравнения Л5(г,с) = 0, коэффициенты которого аналитическим образом зависят от с.

Функция Ф(с) непрерывна при с € [0,1], так как она непрерывна на [0,1/в) и (1/в, 1]. Кроме того, г*(с) ^ 1 при с ^ 1/в + 0, так как в этом случае ^ г*(с) и Л5(гт, с) ^

0. ■

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Непрерывность функции в Ф(с) в точке с = 1/в, где у неё имеется особенность, с физической точки зрения означает, что в этой точке происходит фазовый переход второго рода.

Следствие 4. При с > 1/в функция Ф(с) подчинена неравенству

\в-1

вс ф5 (с) < 1. (10)

□ Так как функция ф(с) - вещественно-аналитическая при с > 1/в, то |Ф;(с)| < то при этих значениях с, и поэтому из уравнения (6) следует

ф'(с) (1 — всфв-1(с)) = фв(с) — 1 . (11)

Это указывает на то, что 1 = всф5-1(с). С другой стороны, при ф(1) = 0, то есть

неравенство (10) выполняется в точке с =1. Поэтому, по непрерывности, ввиду необра-

щаемости в нуль функции (1 — всф5-1(с)) на (1/в, 1], получаем, что (10), действительно, имеет место на этом полуинтервале. ■

Замечание 2. Неравенство (10) можно доказать непосредственно. Так как функция ЛД-, с) выпукла, то есть дЛ5(г, с)/дг монотонно возрастает по г € К, имеем неравенство

^ {дк^,с)\ ^

дг V дг / *—*.(<;)

при г > г*(с). Положим г = гт. Так как в точке г = реализуется отрицатель-

ный минимум, а в точке г* (с) происходит пересечение нулевого уровня слева от точки минимальности, то > г*(с), и поэтому неравенство (12) выполняется. Заметив, что

дЬ8{г,с)\ __а_,

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

-1-0,

получим из (16) неравенство всг^-1(с) — 1 < 0, эквивалентное (10).

Теорема 9. При с > 1/в производная ф;(с) определяется формулой

°'<с) = 1 ) <13> 1 — свф5 1(с)

и строго отрицательна, то есть ф(с) монотонно убывает.

□ Отрицательность следует из того ф(с) < 1 при с > 1/в и неравенства (10). ■

Замечание 3. Монотонное убывание функции ф(с) следует также из монотонного убывания каждой из функций фп(с), п € М+ и предельного соотношения (3), так как предел последовательности монотонно убывающих функций является монотонно невозрастающей функцией.

Следствие 5. Имеет место формула

Ит 0'(с) =--------. (14)

С— 1/5+0 в 1

□ Так как

Пт ^ ~1 = /,

с—— 1/5+0 ф(с) 1

то, используя значение этого предела при I = в, в — 1,

. Ф5(с) — 1 ф(с) — 1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

А = Ит <3 (с) = Ит ----------------- _ = ^ Ит —

с—1/5+0 с—1/5+0 1 — вс ф5 1(с) с—1/5+0 1 — вс ф5 1(с)

1 — св Ит-------------------5 + 1

с—1/5+0 ф(с) — 1

Предел, стоящий в знаменателе последнего выражения, на основе правила Лопиталя, сводится к

1 — вс в

1т1

с—1/5+0 ф(с) — 1 в/А + (в — 1) Отсюда следует уравнение (А = 0)

в

1 = -

в + (в — 1)А

решение которого даёт формулу (14). ■

Конечность правой производной функции ф(с) в точке 1/в, с физической точки зрения, означает, что возникновение новой фазы в точке фазового перехода происходит без появления "критического режима", что является нетипичным в статистической физике.

Теорема 10. При с > 1/в производная ф/;(с) равна

0/7 ч зСГ-2(с)(1 - <3(с))

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

<3<с)= (1-8еСГЧе»> Н‘(С)- <15)

где значения функции Н5 : [1/в, 1] ^ К определяются формулой

Н5(с) = 2вс + (в — 1)ф(с) — (в + 1) •

Она строго положительна, то есть ф(с) выпукла.

в

□ Дифференцируя согласно формуле (13), получаем

О" (с) =_____вСГ2(с)______х

(1 — зс<3Л’-1(с))2

х [ф'(с) (с(в — 1)(Ф5(с) — 1) + Ф(с)(1 — всФ5-1(с))) + ф(с)(ф5(с) — 1)] _ вф5-2(с)(1 — Ф(с))

(1 — всф5-1(с))3

вф5-2(с)(1 — Ф(с))

[с(в — 1)(1 — ф5(с)) — 2ф(с)(1 — вс ф5 1(с))] [2вс + (в — 1)ф(с) — (в + 1)] ,

(1 — всф5-1(с))3

где мы воспользовались уравнением (6) для преобразования выражения в квадратных скобках. Таким образом, получаем формулу (15). Тогда

#5(с) =(в — 1)ф'(с) +2в, #5'(с) = (в — 1)ф''(с) • (16)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Так как (5(1) = 0, ф'(1) = —1, то Н5(1) = #5(1) = (в + 1) и, проинтегрировав (16), получаем интегральное представление для функции Н5,

1 1

#5(с) = (в + 1)(2 — с) + (в — 1) ^ д"(п)^п. (17)

С ?

Пусть с' и с" точки, ближайшие к 1 слева от неё, в которых соответственно Н5(с) и ф''(с) обращаются в нуль. Тогда из (17) следует, что с' < с", и поэтому на интервале (с/,с/') выполняется ф''(с) < 0. С другой стороны, из (15) следует, что это невозможно. Таким образом, Н5(с) > 0 и ф''(с) > 0 одновременно на полуинтервале (1/в, 1]. ■

Пример. Рассмотрим отдельно случай в = 2. Последовательность (фп(с); п € М+) определяется разностным уравнением (см. (4))

фп+1 (с) = 1 — с + с (с)

и начальным условием ф0(с) = 1 — с. На основе этого рекуррентного соотношения, индукцией по п € М+, устанавливается неравенство фп(с) < (2с)-1. А именно, так как ф0(с) = 1 — с < (2с)-1, ввиду с(1 — с) ^ 1/4 при с € [0,1], используя индукционный шаг, находим

Фп+1(с) ^ 1 — с + (4с) 1 ^ (2с) 1 •

Так как при каждом фиксированном с € (0,1) функция фп(с) стремится к вероятности ф(с), то ф(с) ^ (2с)-1.

Таким образом, из двух решений

<3±(с) = 1(1 ± Л/1/4-с(1-с)) = 1(1 ± |1 - 2с\),

{Q+(c) = 1, Q-(c) = (1 — c)/c} уравнения (6) с s = 2,

Q(c) = 1 — c + c Q2(c) ,

при c > 1/2 нужно выбрать то, которое подчинено неравенству Q(c) < (2c)-1 при с > 1/2. Этим свойством обладает только решение Q-(c) = c-1(1 — c). То же самое следует из условия неубывания последовательности (Qn(c); n £ N+). Неравенство Qn+i(c) ^ Qn(c) даёт

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

c Qnn(c) — Qn(c) + 1 — c ^ 0.

Это неравенство может быть удовлетворено только при Qn(c) ^ Q-(c), либо Qn(c) ^ Q+(c) = 1. Второе невозможно, так как Qn(1) = 0.

Непосредственное вычисление

Q'(c) = —1/c2, Q'(1/2 + 0) = —4 , Q"(c) = 2/c3

указывает на то, что Q(c) монотонно убывает и выпукла на [1/2,1].

5. Заключение. В статье математически корректно доказаны основные свойства вероятности перколяции P(c), как функции от концентрации c, для однородного бернул-лиевского случайного поля (или, в другой терминологии, в задаче узлов) на однородных бесконечных древесных графах. Единственной нерешённой проблемой для графов указанного типа остаётся доказательство того факта, что функции Pn(c) для любого n £ N

имеют ровно одну точку перегиба на (0,1), то есть уравнения P^(c) = 0, n £ N имеют

при каждом n £ N на (0,1) ровно одно решение. Этот факт усматривается в компьютерных экспериментах при значениях n, меньших 10, однако, общее доказательство его правильности отсутствует.

Все доказательства в этой работе построены на рекуррентном соотношении (4) для вероятностей Qn(c). Его наличие тесно связано с тем, что граф имеет древесную структуру. Для более сложных графов, в которых имеются циклы, соотношения, подобные

(4), отсутствуют. Поэтому доказательства общих качественных свойств вероятности перколяции таких, как непрерывность, наличие одной точки c*, в которой P(c) не является аналитической, вогнутость P(c) при c > c* в настоящее время отсутствуют.

Литература

1. Kesten H. Percolation Theory for Mathematicians/ H. Kesten. - Boston: Birkhauser, 1982.

2. Grimmett G. Percolation. 2nd Edition/ G. Grimmet. - New York: Springer-Verlag, 1999.

3. Тарасевич Ю.Ю. Перколяция: теория, приложения, алгоритмы/ Ю.Ю. Тарасевич. -Москва: Эдиториал УРСС, 2002.

4. Севастьянов Б.А. Ветвящиеся процессы / Б.А. Севастьянов. - М.: Наука, 1971. - 436 с.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

CONTINUITY OF BERNOULLI RANDOM FIELD PERCOLATION PROBABILITY OF UNIFORM TREE GRAPHS E.S. Antonova Yu.P. Virchenko

Belgorod State University,

Pobedy St., 85, Belgorod, 308015, Russia, e-mail: antonova_e_s@mail.ru

Belgorod State University,

Pobedy St., 85, Belgorod, 308015, Russia, e-mail: virch@bsu.edu.ru

Abstract. The cite problem for arbitrary infinite graphs A(V, $) in discrete percolation theory is considered. Graphs under consideration are so-called uniform Caley trees. The percolation probability P(c) of the uniform Bernoulli random field from the root vertex is studied. It is proved that the function P(c) is continuous at any branching index s € {2, 3, 4,...} of the tree.

Key words: percolation probability, uniform Caley tree, Bernoulli random field.