Ultrametric diffusion in a strong centrally symmetric field

Abstract

A random process at the boundary of a finite regularly branching tree encapsulated in the central-symmetric external field is considered with respect to introduced ultrametricity. We demonstrate an explicit procedure of reduction of dimensionality of the problem. In addition, we consider the strongfield-limit and show that in this case the problem can be solved exactly. The exact solution of the strong-field-limit problem related to the case of linearly growing hierarchy of barriers is exemplified and supplemented by estimations of the transition kinetics into the ground state.

Full Text

Введение. В последние два десятилетия наблюдается рост интереса к случайным процессам над неархимедовыми полями. Если до начала 80-х годов эта тема звучала почти исключительно в математической литературе (см. [2, 3] и цитированную там литературу), то к настоящему моменту ультраметрические конструкции занимают видное место в физике конденсированных сред с беспорядком [4-6], динамике белковых молекул [7, 8], математическом моделировании и дизайне молекулярных машин [9-12], описании компактной укладки ДНК в ядре клетки [13, 14], и наконец, теории эволюции [15]. Складывается впечатление, что всякий раз, когда теория сталкивается с многомасштабной архитектурой, иерархической упорядоченностью и самоподобием, возникает желание использовать ультраметрическое пространство для описания состояний системы и её динамики. В таких случаях часто говорят о «сложных системах», управляющая функция которых (потенциальная или свободная энергия физических систем, или «ландшафт приспосабливаемости» эволюционных систем) представляет собой сильно пересечённую многомерную гиперповерхность с огромным числом локальных © 2015 Самарский государственный технический университет. Образец для цитирования С и з о в а О. М. Ультраметрическая диффузия в сильном центрально-симметричном поле // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2015. T. 19, № 1. С. 87-104. doi: 10.14498/vsgtu1389. ∗ Настоящая статья представляет собой расширенный вариант доклада [1], сделанного автором на Четвёртой международной конференции «Математическая физика и её приложения» (Россия, Самара, 25 августа - 1 сентября 2014). 87 С и з о в а О. М. экстремумов [16, 17]. Интересно, что обсуждение стохастической динамики на сильно пересечённых многомерных ландшафтах и попытки связать её с блужданиями на древообразных графах можно найти в работах генетиковэволюционистов 30-х годов прошедшего столетия (см., например, [18]), пытавшихся проложить путь от комбинаторной природы генетики к древообразному (таксономическому) описанию видового разнообразия. Новую волну интереса к этим идеям можно увидеть в работах середины 80-х годов прошедшего столетия, мотивированных желанием связать ультраметрическую картину фазовых состояний конденсированных систем с беспорядком с характерной для них аномальной (многомодовой) релаксационной кинетикой и использовать этот подход для описания динамики белка и кинетики ферментативной реакции [19-24]. К этому же периоду относится небольшая заметка [25], в которой была поставлена и решена задача об ультраметрическом случайном блуждании на дискретной решётке, оснащённой бинарным деревом ультраметрических расстояний. Эта заметка представляется удобной отправной точкой для постановки задачи настоящей статьи. В работе [25] было рассмотрено случайное блуждание на границе регулярного бинарного дерева из r уровней и 2r листьев, занумерованных от 0 до 2r - 1 (см. рис. 1), с заданной на них ультраметрикой. А именно, расстояние между любыми двумя листьями бинарного дерева (состояниями системы) определяется только номером γ уровня, на котором находится вершина минимального поддерева, содержащего эти листья. При этом считается, что переход из одного состояния в другое происходит в результате преодоления барьера Hγ , высота которого есть функция ультраметрического расстояния между состояниями. Обозначим за fi (t) переходную вероятность случайного блуждания на границе дерева, т. е. вероятность обнаружить блуждающую «частицу» в состоянии i в момент времени t. Будем считать, что константа перехода через барьер Hγ равна e-Hγ , что физически согласуется с законом Аррениуса. Тогда, учитывая энтропийный фактор, определим константу перехода между двумя состояниями, разделёнными барьером Hγ , как qγ = e-Hγ /µ(γ), где µ(γ) = 2γ-1 - мера поддерева, в которое осуществляется переход через данный барьер. Кинетическое уравнение, описывающее этот случайный процесс, имеет вид ∂f (t) = Af (t), ∂t Рис. 1. Бинарное дерево при r = 3 [Figure 1. Binary tree, when r = 3] 88 Ультраметрическая диффузия . . . где f (t) = (f0 (t), . . . , f2r -1 (t)) и A - матрица переходов вида d q1 q2 q2 A = q3 q3 q3 q3 . . . q1 d q2 q2 q3 q3 q3 q3 . . . q2 q2 d q1 q3 q3 q3 q3 . . . q2 q2 q1 d q3 q3 q3 q3 . . . q3 q3 q3 q3 q3 q3 q3 q3 d q1 q1 d q2 q2 q2 q2 . . . . . . q3 q3 q3 q3 q2 q2 d q1 . . . q3 q3 q3 q3 q2 q2 q1 d . . . ... ... ... ... ... , ... ... ... .. . (1) где d = - q1 + 2q2 + 22 q3 + · · · + 2r-1 qr - сумма элементов по столбцу со знаком минус. Заметим, что матрица (1) симметрична, вырождена и имеет блочно-иерархическую структуру. Как было показано в [25], собственные вектора такой матрицы не зависят от элементов qγ . Более того они разбиваются на группы, образуя собственные подпространства Aγ , γ = 0, . . . , r с собственными значениями λγ следующего вида: λ0 = 0 : A0 =   1 1 . ; . . 1 λ1 = -2 q1 + q2 + 2q3 + . . . + 2r-2 qr :       1 0 0 -1  0   .       .  . 0 1       0 A1 =  0  , -1 , . . . ,  0  ; 0 0       0  .   .     .   .  1 . . 0 -1 0 λ2 = -22 q2 + q3 + 2q4 + . . . + 2r-3 qr :       0 0 1 0 1 0 -1  0  0       -1  0  0       0 1 0       0 1  .  A2 =   ,   , . . . ,  .  ;  .   0  -1 0        0  -1 1       0 0 1  .   .  -1  .   .  . . -1 0 0 ... 89 С и з о в а О. М.  λr = -2r qr : Ar =  1 1    .   .  .   1 -1 .   -1    .   .  . -1 Pазлагая начальное условие по собственным подпространствам Aγ , получаем решение задачи Коши: f (0) = a0 + a1 + · · · + ar , λ0 t f (t) = a0 e + a1 e λ1 t a γ ∈ Aγ , + · · · + ar eλr t . Отсюда видно, что стационарное решение представляет собой собственный вектор a0 , соответствующий λ0 = 0. Отметим, что эта задача была также описана в терминах p-адического анализа в [26]. Задачу, рассмотренную в работе [25], нетрудно обобщить на случай дерева с произвольным индексом ветвления m. В этом случае система имеет mr состояний, а матрица переходов A принимает следующий вид: d q1 . . . q1 q2 q2 . . . q2 . . . q2 q2 . . . q2 . . . q1 d . . . q1 q2 q2 . . . q2 . . . q2 q2 . . . q2 . . . ··· ··· .. . ··· ··· ··· . . . ··· . . . ··· ··· . . . ··· . . . q1 q1 . . . d q2 q2 . . . q2 . . . q2 q2 . . . q2 . . . q2 q2 . . . q2 d q1 . . . q1 . . . q2 q2 . . . q2 . . . q2 q2 . . . q2 q1 d . . . q1 . . . q2 q2 . . . q2 . . . · · · q2 · · · q2 . ··· . . · · · q2 · · · q1 · · · q1 . .. . . . ··· d . . . . . . · · · q2 · · · q2 . . . . . . · · · q2 . . . . . . · · · q2 · · · q2 . ··· . . · · · q2 · · · q2 · · · q2 . ··· . . · · · q2 . .. . . . ··· d · · · q1 . . ··· . · · · q1 . . . . . . q2 q2 . . . q2 q2 q2 . . . q2 . . . q1 d . . . q1 . . . · · · q2 · · · q2 . ··· . . · · · q2 · · · q2 · · · q2 . ··· . . · · · q2 . . . . . . · · · q1 · · · q1 . .. . . . ··· d . . . . . . ··· ··· ··· ··· ··· ··· ··· ··· . (2) ··· ··· ··· ··· ··· .. . Здесь d = -(m - 1) q1 + mq2 + . . . + mr-1 qr . Эта матрица также имеет r + 1 собственное подпространство, которые строятся аналогично. А именно, собственный вектор, соответствующий нулевому собственному значению, это вектор (1, 1, . . . , 1). Далее, чтобы построить подпространство Aγ , следует разбить вектор на блоки длины mγ , выбрать один из них, а остальные заполнить нулями. В выбранном блоке необходимо выделить m подблоков длины mγ-1 , каждый i-тый из них заполнить числом xi , так чтобы m xi = 0. i=1 90 Ультраметрическая диффузия . . . Таким способом, можно построить (m-1)mr-γ линейно независимых векторов, порождающих собственное подпространство Aγ с собственным значением λγ = -mγ qγ - (m - 1)mγ (qγ+1 + mqγ+2 + . . . + mr-γ-1 qr ). Аналогично, получаем решение задачи Коши. (Описание этой задачи в терминах m-адического анализа см. в [27].) Следует отметить, что ультраметрическое случайное блуждание с симметричной матрицей переходов соотносится только с весьма узким классом физических систем, у которых все состояния энергетически вырождены. В ряде случаев, как, например, при спектральной диффузии в низкотемпературных белках [8], ультраметрическая диффузия на сильно вырожденном энергетическом ландшафте хорошо описывает наблюдаемые свойства. Однако это скорее счастливое исключение, чем правило. Гораздо чаще энергетические состояния макро- и супрамолекулярных структур не вырождены и имеют хорошо определенное основное состояние (глобальный минимум) c большим бассейном притяжения (см., например, [28, 29]). Динамику таких систем тоже можно пытаться описывать ультраметрическим случайным блужданием при том понимании, что в невырожденном случае вся древообразная структура активационных барьеров помещена в «охватывающий» потенциал, выделяющий бассейн притяжения. Хотя уравнения для таких систем (в p-адической форме) обсуждались, например, в [30-32], существенного продвижения в части методов их решения и получения самих решений добиться пока не удалось. В данной статье мы рассматриваем уравнение для задачи с определённым видом потенциала, а именно, потенциальной ямой с одним основным состоянием, и показываем, что в предельном случае сильного поля удается найти точное аналитическое решение соответствующей задачи Коши. Ультраметрическое блуждание в центрально-симметричном поле. Видоизменим описанное выше симметричное дерево переходов, поместив его в потенциальную яму с минимумом энергии в нулевом состоянии, которое далее будем называть основным состоянием системы (см. рис. 2). Назовём уровнем γ(i) состояния i ультраметрическое расстояние до основного состояния. Обозначим высоту состояния уровня γ относительно нулевого уровня за hγ , а высоту переходной вершины между состояниями i и j за Hij . Тогда барьер перехода из i в j есть Hij - hγ(i) . Положим Hγ(i) = H0i , тогда константа Рис. 2. Нерегулярное дерево при m = 3, r = 3 [Figure 2. Irregular tree, when m = 3, r = 3] 91 С и з о в а О. М. перехода aγ из основного состояния в состояние уровня γ aγ = e-Hγ . (m - 1)mγ-1 Заметим, что барьеры переходов между состояниями одного уровня не меняются по отношению к задаче с симметричным деревом, поэтому для переходов, не увеличивающих уровень состояния, соответствующие константы остаются прежними и имеют вид qγ = e-(Hγ -hγ ) . (m - 1) mγ-1 Вычислим константы переходов для переходов i → j, увеличивающих уровень состояния: qi→j = e-Hγ(j) e-(Hij -hγ(i) ) = ehγ(i) = aj si , (m - 1)mγ(i)-1 (m - 1)mγ(i)-1 где si = ehγ(i) . Запишем матрицу переходов: d0 a1 . . . a1 a2 . . . a2 . . . a2 . . . a2 a3 . . . q1 d1 . . . q1 a2 s 1 . . . a2 s 1 . . . a2 s 1 . . . a2 s 1 a3 s 1 . . . ··· ··· .. . ··· ··· . . . ··· . . . ··· . . . ··· ··· . . . q1 q2 q1 q2 . . . . . . d1 q2 a2 s1 d2 . . . . . . a2 s1 q1 . . . . . . a2 s1 q2 . . . . . . a2 s1 q2 a3 s 1 a3 s 2 . . . . . . ··· ··· ··· ··· ··· .. . ··· . . . ··· . . . ··· ··· . . . q2 q2 . . . q2 q1 . . . d2 . . . q2 . . . q2 a3 s 2 . . . ··· ··· ··· ··· ··· ··· ··· .. . ··· ··· ··· ··· . . . q2 q2 . . . q2 q2 . . . q2 . . . d2 . . . q1 a3 s2 . . . ··· ··· ··· ··· ··· ··· ··· . . . ··· .. . ··· ··· . . . q2 q3 q2 q3 . . . . . . q2 q3 q2 q3 . . . . . . q2 q3 . . . . . . q1 q3 . . . . . . d2 q3 a3 s2 d3 . . . . . . ··· ··· ··· ··· ··· ··· ··· . ··· ··· (3) ··· ··· ··· .. . Здесь также d0 , d1 , . . ., dr - суммы недиагональных элементов по соответствующим столбцам, взятые со знаком минус. Построим собственные векторы такой матрицы. Заметим, что матрица (3) имеет блочно-иерархическую структуру, аналогичную матрице (2), поэтому каждый собственный вектор матрицы (3) можно также разбить на блоки, кратные индексу ветвления m. Более того, некоторые собственные векторы матрицы (3) совпадают с собственными векторами матрицы (2). А именно, это векторы, содержащие нули в первом блоке. Легко найти и некоторые другие векторы следующим образом. Возьмём первый блок размера mk , разобьём его на m подблоков размера k-1 и первый из них заполним нулями. Остальные m - 1 - числами x , m i 92 Ультраметрическая диффузия . . . i = 2, . . . , m, соответственно, так что m xi = 0. Как нетрудно убедиться, i=2 эти векторы также будут собственными для матрицы (3). Остаётся ещё r +1 собственный вектор. В силу введённого в задачу потенциала эти собственные векторы отвечают переходам из основного состояния системы. Поскольку переходы 0 → i зависят только от γ(i), объединим состояния, находящиеся на одном уровне, и построим матрицу переходов для вектора заселённостей объединённых состояний, который также будем обозначать f . Для этого сначала сложим строки матрицы (3), имеющие одинаковые элементы в первом столбце, заменяя соответствующие блоки строк строкой-суммой, а затем удалим повторяющиеся столбцы. В результате по˜ лучим матрицу A размерности (r + 1) × (r + 1):       d0 (m - 1)a1 (m - 1)ma2 . . . q1 d1 (m - 1)ma2 s1 . . . q2 (m - 1)q2 d2 . . .  ... qr . . . (m - 1)qr  . . . (m - 1)mqr  .   . ..  . . . (m - 1)mr-1 ar (m - 1)mr-1 ar s1 (m - 1)mr-1 ar s2 . . . dr (4) В этой матрице также сумма элементов по столбцам равна 0. Таким образом, задача сводится к нахождению собственных векторов и собственных значений матрицы (4). Ниже мы рассмотрим один физически интересный предельный случай, для которого такая задача может быть решена аналитически. Случай сильного поля. Рассмотрим случай, когда запрещены переходы, увеличивающие уровень состояния. Это означает, что a1 = a2 = . . . = ar = 0. ˜ В этом случае матрица A примет вид 0 q1 q2 q3 0 -q1 (m - 1)q2 (m - 1)q3  -mq2 (m - 1)mq3 0 0 0 0 0 -m2 q3  . . . . . . . . . . . . 0 0 0 0   ... qr ... (m - 1)qr   . . . (m - 1)mqr  . . . (m - 1)m2 qr  .   . .. .  . . ... (5) -mr-1 qr В такой задаче разумно взять начальное условие f0 = (0, . . . , 0, 1) и исследовать кинетику заселённости основного состояния (т. е. первую координату вектора f ). Матрица (5) треугольная, поэтому её диагональные элементы есть собственные значения. Собственный вектор, соответсвующий λ0 = 0, находится сразу: v0 = (1, 0, . . . , 0). Рассматривая структуру матрицы (5), нетрудно понять, что достаточно вычислить собственный вектор vr , соответствующий λr = -mr-1 qr , т. е. найти ненулевое решение матричного уравнения ∂vr ˜ = A - λr E vr . ∂t (6) 93 С и з о в а О. М. Остальные векторы vk находятся аналогичным образом заменой r = k и заполнением координат xi , i > k, собственного вектора нулями. Решение уравнения (6) находится последовательно методом Гаусса: xr-2 xr = 1, (m - 1)qr , xr-1 = - mqr - qr-1 m(m - 1)qr (qr - qr-1 ) =- , (mqr - qr-1 ) (m2 qr - qr-2 ) ... . Докажем, что при k = 0 xr-k = - mk-1 (m - 1)qr (qr - qr-1 ) (mqr - qr-2 ) . . . mk-2 qr - qr-k+1 . (mqr - qr-1 ) (m2 qr - qr-2 ) . . . (mk qr - qr-k ) (7) Из соответствующей строки матрицы следует xr-k = - mr-k-1 (m - 1) [qr xr + qr-1 xr-1 + . . . + qr-k+1 xr-k+1 ] . mr-1 qr - mr-k-1 qr-k Вычислим последовательно сумму в квадратных скобках. mqr (qr - qr-1 ) , mqr - qr-1 m2 qr (qr - qr-1 ) (mqr - qr-2 ) , = (mqr - qr-1 ) (m2 qr - qr-2 ) ... . qr xr + qr-1 xr-1 = qr xr + qr-1 xr-1 + qr-2 xr-2 По индукции получим: qr xr + . . . + qr-k+1 xr-k+1 = mk-1 qr (qr - qr-1 ) . . . mk-2 qr - qr-k+1 . (mqr - qr-1 ) . . . (mk-1 qr - qr-k+1 ) Подставляя сумму в выражение (7), получим xr-k . Аналогично можно вывести формулу для x0 : x0 = - (qr - qr-1 ) (mqr - qr-2 ) . . . mr-2 qr - q1 . (mqr - qr-1 ) (m2 qr - qr-2 ) . . . (mr-1 qr - q1 ) Для решения задачи Коши остаётся разложить начальное условие по построенному базису размерности r + 1, решив матричное уравнение Bc = f0 (8) относительно вектора коэффициентов разложения c, где B - матрица, состоящая из столбцов vk . Это уравнение также решается методом Гаусса. 94 Ультраметрическая диффузия . . . При произвольном выборе величин qi для матрицы (5) затруднительно получить компактные выражения для вектора c, но их можно получить в случае линейного роста энергий Hγ - hγ от уровня γ. Случай сильного линейного поля. Зададим барьеры переходов Hγ - hγ = ˜ = γ. Тогда матрица A и собственный вектор vr примут вид   1 -3 1 1 1 -2 e . . . mr-2 e-(r-1) mr-1 e-r 0 e-1 me m2   0 -e-1 m-1 e-2 m-1 e-3 . . . m-1 e-(r-1) m-1 e-r  m m2 mr-2 mr-1    m-1 -(r-1) m-1 -3 m-1 -r  -2 0 0 -e . . . mr-3 e e  m e mr-2    m-1 -r  m-1 -(r-1) -3 0 e , 0 0 -e . . . mr-4 e mr-3    . . . . . . ..  . . . . . . . . . . . .  .    m-1 -r  e  0 0 0 ... -e-(r-1) 0 m 0 0 0 0  1 - mr-1 ... -e-r 0 (mer-1 -1)...(me-1)  (er-1 -1)...(e-1)    m-1 (mer-2 -1)...(me-1)   r-1 r-1   m (e -1)(er-2 -1)...(e-1)    .   .   .   mek-1 -1)...(me-1)   m-1 (  mk ek -1 ek-1 -1 ...(e-1)  ( )( )    . vr =  .   . .      m-1 (me2 -1)(me-1)    m3 (e3 -1)(e2 -1)(e-1)     me-1 m-1   m2 (e2 -1)(e-1)     m-1 1   m e-1 1 Для получения окончательного решения задачи Коши разложим начальное условие по собственному базису. Запишем систему (8) в матричном виде и решим её:  2 r-1 -1 ...(me-1)  ) 1 (me -1)(me-1) 1 (me 1 . . . - mr-1 (er-1 -1)...(e-1) 1 -1 - m me-1 - m2 (e2 -1)(e-1) e-1    r-2 -1 ...(me-1)  ) m-1 me-1 m-1 (me m-1 1 0 1 . . . mr-1 (er-1 -1)...(e-1)    m e-1 m2 (e2 -1)(e-1)    mer-3 -1)...(me-1)  m-1 1 m-1 ( 0 0 1 . . . mr-2 (er-2 -1)...(e-1)  m e-1   r-4 -1 ...(me-1)   ) × m-1 (me 0 0 0 1 . . . mr-3 (er-3 -1)...(e-1)    .  . . . . ..  . . . . . .  . . . . .     m-1 1 0 0 ...  0 0 m e-1 0 0 0 0 ... 1 95 С и з о в а О. М.    0 a0  a1  0      a2  0     ×  a3  = 0 . (9)  .  .  .  .  .  . ar-1  0 1 ar  Будем находить последовательно значения ai : ar = 1, ar-1 = - m-1 1 , m e-1 ar-2 = - m-1 e-m , m2 (e2 - 1) (e - 1) ... . Докажем, что для k < r ar-k = - ek-1 - m . . . (e - m) m-1 . mk (ek - 1) (ek-1 - 1) . . . (e - 1) Предположим, что это верно для ar-i при всех i < k. Покажем, что равенство верно и для i = k. Из матрицы (9) находим: ek-2 - m . . . (e - m) m-1 1 m-1 · k-1 k-1 + m e-1 m (e - 1) (ek-2 - 1) . . . (e - 1) ek-3 - m . . . (e - m) m-1 me - 1 m-1 + · k-2 k-2 + m2 (e2 - 1) (e - 1) m (e - 1) (ek-3 - 1) . . . (e - 1) mek-2 - 1 . . . (me - 1) m-1 m-1 1 + . . . + k-1 k-1 · - k-2 - 1) . . . (e - 1) m e-1 m (e - 1) (e mek-1 - 1 . . . (me - 1) m-1 - . mk (ek - 1) (ek-1 - 1) . . . (e - 1) ar-k = Дальнейшее доказательство является прямым следствие следующей леммы. Лемма. Для любых x = 1, m ∈ R и k ∈ N верно равенство mxk-1 - 1 . . . (mx - 1) xk-1 - m . . . (x - m) = k + (xk - 1) (xk-1 - 1) . . . (x - 1) (x - 1) (xk-1 - 1) . . . (x - 1) xk-2 - m . . . (x - m) 1 + (m - 1) · k-1 + x - 1 (x - 1) (xk-2 - 1) . . . (x - 1) xk-3 - m . . . (x - m) mx - 1 + (m - 1) 2 · k-2 + (x - 1) (x - 1) (x - 1) (xk-3 - 1) . . . (x - 1) mxk-2 - 1 . . . (mx - 1) 1 + . . . + (m - 1) k-1 · . (x - 1) (xk-2 - 1) . . . (x - 1) x - 1 Д о к а з а т е л ь с т в о. Докажем лемму по индукции. В правой части равенства - - k слагаемых, которые обозначим bi , 1 i k. Сложим первые два слагаемых: 96 Ультраметрическая диффузия . . . xk-2 - m . . . (x - m) xk-1 - m m - 1 + = x-1 (xk-1 - 1) (xk-2 - 1) . . . (x - 1) xk - 1 xk-2 - m . . . (x - m) xk-1 - 1 (mx - 1) = k-1 · = (x - 1) (xk-2 - 1) . . . (x - 1) (xk - 1) (x - 1) 1 mx - 1 xk-2 - m . . . (x - m) = k · · . x -1 x-1 (xk-2 - 1) . . . (x - 1) b1 + b2 = Пусть для суммы первых n слагаемых верна формула b1 + . . . + bn = mxn-1 - 1 . . . (mx - 1) xk-n - m . . . (x - m) 1 · · , (xn-1 - 1) . . . (x - 1) xk - 1 (xk-n - 1) . . . (x - 1) докажем е¨ для n + 1 слагаемых. Имеем е bn+1 = (m-1) mxn-1 - 1 . . . (mx - 1) xk-n-1 - m . . . (x - m) ; · k-n (xn - 1) (xn-1 - 1) . . . (x - 1) (x - 1) (xk-n-1 - 1) . . . (x - 1) b1 + . . . + bn+1 = mxn-1 - 1 . . . (mx - 1) xk-n-1 - m . . . (x - m) · k-n × (xn-1 - 1) . . . (x - 1) (x - 1) (xk-n-1 - 1) . . . (x - 1) xk-n - m m-1 × + n = x -1 xk - 1 mxn-1 - 1 . . . (mx - 1) xk-n-1 - m . . . (x - m) = · k-n × (xn-1 - 1) . . . (x - 1) (x - 1) (xk-n-1 - 1) . . . (x - 1) xk-n - 1 (mxn - 1) × = (xk - 1) (xn - 1) 1 (mxn - 1) . . . (mx - 1) xk-n-1 - m . . . (x - m) = k · · . (xn - 1) . . . (x - 1) x -1 (xk-n-1 - 1) . . . (x - 1) = Теперь, полагая n = k, получим b1 + . . . + bk = mxk-1 - 1 . . . (mx - 1) . (xk - 1) (xk-1 - 1) . . . (x - 1) Что и требовалось доказать. Итак, остаётся найти a0 . Складывая все строки матрицы (9), получим строку (1, 0, . . . , 0) и единицу в столбце свободных членов. Отсюда a0 = 1. Теперь можно выписать решение задачи Коши в терминах заселённостей уровней:     1 -1 0 1     0 m - 1  0  -e-1 t er-2 - m . . . (e - m) 0 -  e f (t) =   - r-1 (er-1 - 1) (er-2 - 1) . . . (e - 1)  0  m .  .  .  .  . . 0 0 97 С и з о в а О. М.  1 - m me-1 e-1     er-3 - m . . . (e - m) m-1  - r-2 r-2 r-3 - 1) . . . (e - 1)  m (e - 1) (e    m-1 1 m e-1      -e-2 t e -     1 0 . . . 0   2 1 (me -1)(me-1) - m2 (e2 -1)(e-1)   m-1 me-1   m2 (e2 -1)(e-1)      -e-3 t er-4 - m . . . (e - m) m-1 m-1 1  e - r-3 r-3 - ... m e-1  r-4 - 1) . . . (e - 1)  m (e - 1) (e   1     . .   . 0   mer-2 -1)(mer-3 -1)...(me-1) 1 ( - mr-2 (er-2 -1)(er-3 -1)...(e-1)     (mer-3 -1)...(me-1) m-1   r-2 (er-2 -1)(er-3 -1)...(e-1)   m   m-1 1   -e-(r-1) t (mer-4 -1)...(me-1) m-1 + ... -  e mr-3 (er-3 -1)(er-4 -1)...(e-1)  m e-1   .   . .     1 0   r-1 -1 mer-2 -1 ...(me-1) )( ) 1 (me -  mr-1 (er-1 -1)(er-2 -1)...(e-1)    (mer-2 -1)...(me-1) m-1     mr-1 (er-1 -1)(er-2 -1)...(e-1)   r-3 -1 ...(me-1)   (me ) m-1  -e-r t  mr-2 (er-2 -1)(er-3 -1)...(e-1) . + e   .   . .       m-1 1   m e-1 1 Кинетика перехода в основное состояние. Рассмотрим первую координату вектора f (t) - зависимость заселённости основного состояния от времени, обозначив её S(t): S(t) = 1 + er-2 - m . . . (e - m) m-1 -1 · e-e t + r-1 r-1 - 1) (er-2 - 1) . . . (e - 1) m (e + 98 er-3 - m . . . (e - m) me - 1 -e-2 t ·e + ... r-2 - 1) (er-3 - 1) . . . (e - 1) e - 1 (e Ультраметрическая диффузия . . . ... + mer-2 - 1 mer-3 - 1 . . . (me - 1) -e-(r-1) t 1 ·e - e-1 (er-2 - 1) (er-3 - 1) . . . (e - 1) - 1 mr-1 mer-1 - 1 mer-2 - 1 . . . (me - 1) -e-r t ·e . (er-1 - 1) (er-2 - 1) . . . (e - 1) При помощи численных расчётов в программе Wolfram Mathematica была построена асимптотика функции S (t) на временном промежутке e1 , er-1 . Как выяснилось, S (t) ∼ Cta , где a = ln m + 1. Это означает, что в случае сильного поля, т.е. достаточно глубокого основного состояния, кинетика перехода имеет степенное поведение на промежутке характерных времён. График поведения функции S (t) при m = 3, r = 15 показан на рис. 3, а на на рис. 4 изображён график S (t) в логарифмической шкале, наложенный на график соответствующей степенной функции. Рис. 3. График поведения релаксационной функции S(t) при m = 3, r = 15 [Figure 3. The behaviour of the relaxation function S(t), when m = 3, r = 15] Рис. 4. График функции S(t) (сплошная линия) в сравнении с функцией 3 · 10-14 tln 3+1 (штриховая линия) в логарифмической шкале [Figure 4. A log-log plot of the function S(t) (solid line) against the function 3 · 10-14 tln 3+1 (dashed line)] Заключение. Данный результат может представлять интерес для понимания особенностей кинетики коллапса полимерной цепи из клубкового состояния в глобулярное, в частности, для фолдинга белка. Действительно, следуя 99 С и з о в а О. М. общепринятым представлениям (см., например, [33, c. 334], время перехода клубок-глобула можно характеризовать максимальным барьером, лежащим на пути перехода. Однако по тем же представлениям сам переход может пониматься как движение системы на многомерном ландшафте с сильно пересеченным энергетическим рельефом. Как показано выше, такие представления не противоречат друг другу, только если основное состояние достаточно глубокое. Благодарности. Автор выражает благодарность Владику Аванесовичу Аветисову за обсуждение задачи и полезные замечания.