Динамика файла

редактировать

Термин динамика файла - это движение множества частиц в узком канале.

В науке: в химии, физике, математике и родственных областях, динамике файлов (иногда называют, динамика одного файла ) - это диффузия N (N → ∞) идентичных броуновских твердых сфер в квазиодномерном канале длиной L (L → ∞), так что сферы не прыгают одна на другую, а средняя плотность частиц приблизительно фиксирована. Наиболее известные статистические свойства этого процесса заключаются в следующем: среднеквадратичное смещение (MSD) частицы в файле: MSD ≈ t 1 2 {\ displaystyle \ mathrm {MSD} \ приблизительно t ^ {\ frac {1} {2}}}{\ displaystyle \ mathrm {MSD} \ приблизительно t ^ {\ frac {1} {2}}} , и его функция плотности вероятности (PDF) равна гауссову в позиции с дисперсией MSD.

Результаты в файлах, которые обобщают базовый файл, включают:

  • В файлах с законом плотности, который не фиксирован, но затухает по степенному закону с показателем a с расстоянием от начала координат, частица в начале координат имеет МСД, который масштабируется как МСД ≈ t 1 + a 2 {\ displaystyle MSD \ приблизительно t ^ {\ frac {1 + a} {2}}}{\ displaystyle MSD \ приблизительно t ^ {\ frac {1 + a} {2}}} с гауссовским PDF.
  • Когда, кроме того, коэффициенты диффузии частиц распределены по степенному закону с показателем γ (вокруг начала координат), MSD следует, MSD ≈ t 1 - γ 2 / (1 + a) - γ {\ displaystyle MSD \ приблизительно t ^ {\ frac {1- \ gamma} {2 / (1 + a) - \ gamma}}}{\ displaystyle MSD \ приблизительно t ^ { \ frac {1- \ gamma} {2 / (1 + a) - \ gamma}}} , с гауссовским PDF.
  • В аномало us файлы, которые обновляются, а именно, когда все частицы пытаются совершить прыжок вместе, однако, со временем перехода, взятым из распределения, которое затухает по степенному закону с показателем −1 - α, MSD масштабируется как MSD соответствующей нормальной файл в степени α.
  • В аномальных файлах независимых частиц MSD работает очень медленно и масштабируется как, MSD ≈ log 2 (t) {\ displaystyle MSD \ приблизительно log ^ {2 } (t)}{\ displaystyle MSD \ приблизительно log ^ {2} (t)} . Что еще более интересно, частицы образуют кластеры в таких файлах, определяя динамический фазовый переход. Это зависит от степени аномалии α: процент частиц в кластерах ξ следует, ξ ≈ 1 - α 3 {\ displaystyle \ xi \ приблизительно {\ sqrt {1- \ alpha ^ {3}}}}{\ displaystyle \ xi \ приблизительно {\ sqrt {1- \ alpha ^ {3}}}} .
  • Другие обобщения: когда частицы могут обходить друг друга с постоянной вероятностью при встрече, наблюдается усиленная диффузия. Когда частицы взаимодействуют с каналом, наблюдается более медленная диффузия. Файлы, встроенные в двухмерные, демонстрируют аналогичные характеристики файлов в одном измерении.

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

Содержание
  • 1 Математическая формулировка
    • 1.1 Простые файлы
    • 1.2 Гетерогенные файлы
    • 1.3 Обновленные, аномальные, гетерогенные файлы
    • 1.4 Аномальные файлы с независимыми частицами
  • 2 Математический анализ
    • 2.1 Простые файлы
    • 2.2 Гетерогенные файлы
    • 2.3 Обновление аномальных разнородных файлов
    • 2.4 Аномальные файлы с независимыми частицами
      • 2.4.1 Законы масштабирования для аномальных файлов независимых частиц
      • 2.4.2 Численные исследования аномальных файлов независимых частиц
  • 3 См. Также
  • 4 Ссылки
Математическая формулировка

Простые файлы

В простых броуновских файлах P (x, t ∣ x 0) {\ displaystyle P (\ mathbf {x}, t \ mid \ mathbf {x_ {0}})}{\ displaystyle P (\ mathbf {x}, t \ mid \ mathbf {x_ {0}})} , совместная вероятность den Функция sity (PDF) для всех частиц в файле подчиняется уравнению нормальной диффузии:

∂ t P (x, t ∣ x 0) = D Σ j = - MM ∂ xj 2 P (x, t ∣ х 0). {\ displaystyle \ partial _ {t} P (\ mathbf {x}, t \ mid \ mathbf {x_ {0}}) = D \ Sigma _ {j = -M} ^ {M} \ partial _ {x_ { j}} ^ {2} P (\ mathbf {x}, t \ mid \ mathbf {x_ {0}}).}{\ displaystyle \ partial _ {t} P (\ mathbf {x}, t \ mid \ mathbf {x_ {0}}) = D \ Sigma _ {j = -M} ^ {M} \ partial _ {x_ {j}} ^ {2} P (\ mathbf {x}, t \ mid \ mathbf {x_ {0}}).}

(1)

In P (x, t ∣ x 0) {\ displaystyle P (\ mathbf {x}, t \ mid \ mathbf {x_ {0}})}{\ displaystyle P (\ mathbf {x}, t \ mid \ mathbf {x_ {0}})} , x = {x - M, x - M + 1,…, x M} {\ displaystyle \ mathbf {x} = \ {x _ {- M}, x _ {- M + 1}, \ ldots, x_ {M} \}}{\ displaystyle \ mathbf {x} = \ {x _ {- M}, x _ {- M + 1}, \ ldots, x_ {M} \}} - это набор положений частиц в момент t {\ displaystyle t}t и x 0 {\ displaystyle \ mathbf {x_ {0}}}\ mathbf {x_ {0}} - набор начальных положений частиц в начальное время t 0 {\ displaystyle t_ {0}}t_ {0} (установить на ноль). Уравнение (1) решается с соответствующими граничными условиями, которые отражают твердосферный характер файла:

(D ∂ xj P (x, t ∣ x 0)) xj = xj + 1 = (D ∂ xj + 1 P (x, t ∣ x 0)) xj + 1 = xj; j = - M,…, M - 1, {\ displaystyle {\ big (} D \ partial _ {x_ {j}} P (\ mathbf {x}, t \ mid \ mathbf {x_ {0}}) { \ big)} _ {x_ {j} = x_ {j + 1}} = {\ big (} D \ partial _ {x_ {j + 1}} P (\ mathbf {x}, t \ mid \ mathbf { x_ {0}}) {\ big)} _ {x_ {j + 1} = x_ {j}}; \ qquad j = -M, \ ldots, M-1,}{\ displaystyle {\ big (} D \ partial _ {x_ {j}} P (\ mathbf { x}, t \ mid \ mathbf {x_ {0}}) {\ big)} _ {x_ {j} = x_ {j + 1}} = {\ big (} D \ partial _ {x_ {j + 1) }} P (\ mathbf {x}, t \ mid \ mathbf {x_ {0}}) {\ big)} _ {x_ {j + 1} = x_ {j}}; \ qquad j = -M, \ ldots, M-1,}

(2)

и с соответствующим начальным условием:

P (x, t → ∞ ∣ x 0) = Π j = - MM δ (xj - x 0, j). {\ displaystyle P (\ mathbf {x}, t \ rightarrow \ infty \ mid x_ {0}) = \ Pi _ {j = -M} ^ {M} \ delta (x_ {j} -x_ {0, j }).}{\ displaystyle P (\ mathbf {x}, т \ rightarrow \ infty \ mid x_ {0}) = \ Pi _ {j = -M} ^ {M} \ delta (x_ {j} -x_ {0, j}).}

(3)

В простом файле начальная плотность фиксирована, а именно: x 0, j = j Δ {\ displaystyle x_ {0, j} = j \ Delta}{\ displaystyle x_ {0, j} = j \ Delta} , где Δ {\ displaystyle \ Delta}\ Delta - параметр, представляющий микроскопическую длину. Координаты PDF-файлов должны соответствовать порядку: x - M ≤ x - M + 1 ≤ ⋯ ≤ x M {\ displaystyle x _ {- M} \ leq x _ {- M + 1} \ leq \ cdots \ leq x_ {M}}{\ displaystyle x _ {- M} \ leq x _ {- M + 1} \ leq \ cdots \ leq x_ {M}} .

Гетерогенные файлы

В таких файлах уравнение движения следует,

∂ t P (x, t ∣ x 0) = Σ j = - MMD j ∂ xj 2 P (х, t ∣ х 0). {\ Displaystyle \ partial _ {t} P (\ mathbf {x}, t \ mid \ mathbf {x_ {0}}) = \ Sigma _ {j = -M} ^ {M} D_ {j} \ partial _ {x_ {j}} ^ {2} P (\ mathbf {x}, t \ mid \ mathbf {x_ {0}}).}{\ displaystyle \ partial _ {t} P (\ mathbf {x}, t \ mid \ mathbf {x_ {0}}) = \ Sigma _ {j = -M} ^ {M} D_ {j} \ partial _ {x_ {j}} ^ {2} P (\ mathbf {x}, т \ mid \ mathbf {x_ {0}}).}

(4)

с граничными условиями:

( D j ∂ xj P (x, t ∣ x 0)) xj = xj + 1 = (D j + 1 ∂ xj + 1 P (x, t ∣ x 0)) xj + 1 = xj; j = - M,…, M - 1, {\ displaystyle {\ big (} D_ {j} \ partial _ {x_ {j}} P (\ mathbf {x}, t \ mid \ mathbf {x_ {0}) }) {\ big)} _ {x_ {j} = x_ {j + 1}} = {\ big (} D_ {j + 1} \ partial _ {x_ {j + 1}} P (\ mathbf {x }, t \ mid \ mathbf {x_ {0}}) {\ big)} _ {x_ {j + 1} = x_ {j}}; \ qquad j = -M, \ ldots, M-1,}{\ displaystyle {\ big (} D_ {j} \ partial _ {x_ {j}} P ( \ mathbf {x}, t \ mid \ mathbf {x_ {0}}) {\ big)} _ {x_ {j} = x_ {j + 1}} = {\ big (} D_ {j + 1} \ частичное _ {x_ {j + 1}} P (\ mathbf {x}, t \ mid \ mathbf {x_ {0}}) {\ big)} _ {x_ {j + 1} = x_ {j}}; \ qquad j = -M, \ ldots, M-1,}

(5)

и с начальным условием, уравнение. (3), где начальные положения частиц подчиняются:

x 0, j = s i g n (j) Δ ∣ j ∣ 1 / (1 - a); 0 ≤ a ≤ 1. {\ displaystyle x0, j = sign (j) \ Delta \ mid j \ mid ^ {1 / (1-a)}; \ qquad 0 \ leq {\ mathit {a}} \ leq 1.}{\ Displaystyle x0, j = знак (j) \ Delta \ mid j \ mid ^ {1 / (1-a)}; \ qquad 0 \ leq {\ mathit {a}} \ leq 1.}

(6)

Коэффициенты распространения файла взяты независимо от PDF,

W (D) = 1 - γ Λ (D / Λ) - γ, 0 ≤ γ ≤ 1, {\ displaystyle W (D) = {\ frac {1- \ gamma} {\ Lambda}} (D / \ Lambda) ^ {- \ gamma}, \ qquad 0 \ leq \ gamma \ leq 1,}{\ displaystyle W (D) = {\ frac {1- \ gamma} {\ Lambda}} (D / \ Lambda) ^ {- \ gamma}, \ qquad 0 \ leq \ gamma \ leq 1,}

(7)

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

Обновление, аномальные, гетерогенные файлы

В файлах с аномальным возобновлением случайный период берется независимо от функции плотности вероятности времени ожидания (WT-PDF; см. Марков с непрерывным временем процесс для получения дополнительной информации) вида: ψ α (t) ∼ k (kt) - 1 - α, 0 < α ≤ 1 {\displaystyle \psi _{\alpha }(t)\sim k(kt)^{-1-\alpha },0<\alpha \leq 1}{\ displaystyle \ psi _ {\ альфа} (т) \ сим к (кт) ^ {- 1- \ альфа}, 0 <\ альфа \ leq 1} , где k - параметр. Затем все частицы в файле остаются неподвижными на этот случайный период, после чего все частицы пытаются прыгнуть в соответствии с правилами файла. Эта процедура повторяется снова и снова. Уравнение движения для PDF частиц в файле с аномальным обновлением получается при свертывании уравнения движения для броуновского файла с ядром k α (t) {\ displaystyle k _ {\ alpha} (t)}{\ displaystyle k _ {\ alpha} (t)} :

∂ t P (x, t ∣ x 0) = Σ j = - MMD j ∂ xj 2 ∫ 0 tk α (t - u) P (x, u ∣ x 0) du. {\ Displaystyle \ partial _ {t} P (\ mathbf {x}, t \ mid \ mathbf {x_ {0}}) = \ Sigma _ {j = -M} ^ {M} D_ {j} \ partial _ {x_ {j}} ^ {2} \ int _ {0} ^ {t} k _ {\ alpha} (tu) P (\ mathbf {x}, u \ mid \ mathbf {x_ {0}}) \, du.}{\ displaystyle \ partial _ {t} P (\ mathbf {x}, t \ mid \ mathbf {x_ {0}}) = \ Sigma _ {j = -M} ^ {M} D_ {j} \ partial _ {x_ {j}} ^ {2} \ int _ {0} ^ {t} k _ {\ alpha} (tu) P (\ mathbf {x}, u \ mid \ mathbf {x_ {0}}) \, du. }

(8)

Здесь ядро ​​k α (t) {\ displaystyle k _ {\ alpha} (t)}{\ displaystyle k _ {\ alpha} (t)} и WT-PDF ψ α (t) {\ displaystyle \ psi _ {\ alpha} (t)}{\ displaystyle \ psi _ {\ alpha} (t)} связаны в пространстве Лапласа, k ¯ α (s) = s ψ ¯ α (s) 1 - ψ ¯ α (s) {\ Displaystyle {\ bar {k}} _ {\ alpha} (s) = {\ frac {s {\ bar {\ psi}} _ {\ alpha} (s)} {1- {\ bar {\ psi}} _ {\ alpha} (s)}}}{\ displaystyle {\ bar {k}} _ {\ alpha} (s) = {\ frac {s {\ bar {\ psi}) } _ {\ alpha} (s)} {1 - {\ bar {\ psi}} _ {\ alpha} (s)}}} . (Преобразование Лапласа функции f (t) {\ displaystyle f (t)}f (t) читает, f ¯ (s) = ∫ 0 ∞ f (t) e - stdt { \ displaystyle {\ bar {f}} (s) = \ int _ {0} ^ {\ infty} f (t) e ^ {- st} \, dt}{\ displaystyle {\ bar {f}} (s) = \ int _ {0} ^ {\ infty} f (t) e ^ {- st } \, dt} .) Отражающие граничные условия сопровождал уравнение. (8) получаются при свертывании граничных условий броуновского файла с ядром k α (t) {\ displaystyle k _ {\ alpha} (t)}{\ displaystyle k _ {\ alpha} (t)} , где здесь и в броуновский файл начальные условия идентичны.

Аномальные файлы с независимыми частицами

Когда каждой частице в аномальном файле назначается свое собственное время перехода, нарисованная форма ψ α (t) {\ displaystyle \ psi _ {\ alpha } (t)}{\ displaystyle \ psi _ {\ alpha} (t)} (ψ α (t) {\ displaystyle \ psi _ {\ alpha} (t)}{\ displaystyle \ psi _ {\ alpha} (t)} одинаково для всех частиц), аномальный файл не является файлом обновления. Базовый динамический цикл в таком файле состоит из следующих шагов: частица с самым быстрым временем перехода в файле, скажем, ti {\ displaystyle t_ {i}}t_ {i} для частицы i, попытки прыжок. Затем время ожидания для всех остальных частиц корректируется: мы вычитаем t i {\ displaystyle t_ {i}}t_ {i} из каждой из них. Наконец, для частицы i отображается новое время ожидания. Наиболее существенное различие между аномальными файлами обновления и аномальными файлами, которые не являются обновлением, заключается в том, что, когда каждая частица имеет свои собственные часы, частицы фактически связаны также во временной области, и результатом является дальнейшее замедление работы системы (доказано в основной текст). Уравнение движения для PDF в аномальных файлах независимых частиц имеет следующий вид:

∂ ti P (x, t ∣ x 0) = D i ∂ xi 2 ∫ 0 tik α (ti - ui) P (x, t ′ (i), ui ∣ x 0) dui; - M ≤ i ≤ M. {\ displaystyle \ partial _ {t_ {i}} P (\ mathbf {x}, \ mathbf {t} \ mid \ mathbf {x_ {0}}) = D_ {i} \ partial _ {x_ {i}} ^ {2} \ int _ {0} ^ {t_ {i}} k _ {\ alpha} (t_ {i} -u_ {i}) P (\ mathbf {x}, \ mathbf {t} ^ {'( i)}, u_ {i} \ mid \ mathbf {x_ {0}}) \, du_ {i}; \ qquad -M \ leq i \ leq M.}{\displaystyle \partial _{t_{i}}P(\mathbf {x},\mathbf {t} \mid \mathbf {x_{0}})=D_{i}\partial _{x_{i}}^{2}\int _{0}^{t_{i}}k_{\alpha }(t_{i}-u_{i})P(\mathbf {x},\mathbf {t} ^{'(i)},u_{i}\mid \mathbf {x_{0}})\,du_{i};\qquad -M\leq i\leq M.}

(9)

Обратите внимание, что аргумент времени в PDF P (x, t ∣ x 0) {\ displaystyle P (\ mathbf {x}, \ mathbf {t} \ mid \ mathbf {x_ {0}})}{\ displaystyle P (\ mathbf {x}, \ mathbf {t} \ mid \ mathbf {x_ {0}})} - вектор времен: t = {ti} i = - MM {\ displaystyle \ mathbf {t} = \ {t_ {i} \} _ {i = -M} ^ {M}}{\ displaystyle \ mathbf {t} = \ {t_ {i} \} _ {i = -M} ^ {M }} и t ′ (i) = {tc} c = - M, c ≠ i M {\ displaystyle \ mathbf {t} ^ {'(i)} = \ {t_ {c} \ } _ {c = -M, c \ neq i} ^ {M}}{\displaystyle \mathbf {t} ^{'(i)}=\{t_{c}\}_{c=-M,c\neq i}^{M}}. Сложение всех координат и выполнение интегрирования сначала в порядке меньшего времени (порядок определяется случайным образом из равномерного распределения в пространстве конфигураций) дает полное уравнение движения в аномальных файлах независимых частиц (усреднение уравнения по всем поэтому требуется дополнительная конфигурация). Действительно, даже уравнение. (9) очень сложно, и усреднение еще больше усложняет ситуацию.

Математический анализ

Простые файлы

Решение уравнений. (1) - (2) - полный набор перестановок всех начальных координат, входящих в гауссианы,

P (x, ∣ x 0) = 1 c N Σ {p} e - 1/4 D t Σ j = - MM (xj - x 0, j (p)) 2. {\ Displaystyle P (\ mathbf {x}, \ mid \ mathbf {x_ {0}}) = {\ frac {1} {c_ {N}}} \ Sigma _ {\ {p \}} e ^ {- 1 / 4Dt \ Sigma _ {j = -M} ^ {M} (x_ {j} -x_ {0, j} (p)) ^ {2}}.}{\ displaystyle P (\ mathbf {x}, \ mid \ mathbf {x_ {0}}) = {\ frac {1} {c_ {N}}} \ Sigma _ {\ {p \}} e ^ {- 1 / 4Dt \ Sigma _ {j = -M} ^ {M} (x_ {j} -x_ {0, j} (p)) ^ {2}}.}

(10)

Здесь, индекс p {\ displaystyle p}p проходит по всем перестановкам начальных координат и содержит N! {\ displaystyle N!}N! перестановки. Из уравнения. (10), PDF-файл помеченной частицы в файле, P (r, t ∣ r 0) {\ displaystyle P (r, t \ mid r_ {0})}{\ displaystyle P (r, t \ mid r_ {0})} , вычисляется

P (r, t ∣ r 0) ∼ 1 c N e - R d 2/2 τ, {\ displaystyle P (r, t \ mid r_ {0}) \ sim {\ frac {1 } {c_ {N}}} e ^ {- R_ {d} ^ {2} / {\ sqrt {2 \ tau}}},}{\ displaystyle P (r, t \ mid r_ {0}) \ sim {\ frac {1} {c_ {N}}} e ^ {- R_ {d} ^ {2} / {\ sqrt {2 \ tau}}},}

(11)

В уравнении. (11), R d = rd Δ {\ displaystyle R_ {d} = r_ {d} \ Delta}{\ displaystyle R_ { d} = r_ {d} \ Delta} , rd = r - r 0 {\ displaystyle r_ {d} = r-r_ { 0}}{\ displaystyle r_ {d} = r -r_ {0}} (r 0 {\ displaystyle r_ {0}}r_ {0} - начальное состояние помеченной частицы) и τ = Δ - 2 D t {\ displaystyle \ tau = \ Дельта ^ {- 2} Dt}{\ displaystyle \ тау = \ Delta ^ {- 2} Dt} . MSD для меченой частицы получается непосредственно из уравнения. (11):

R d 2⟩ ∼ τ. {\ displaystyle \ langle R_ {d} ^ {2} \ rangle \ sim {\ sqrt {\ tau}}.}{\ displaystyle \ langle R_ {d} ^ {2} \ rangle \ sim {\ sqrt {\ tau}}.}

(12)

Гетерогенные файлы

Решение уравнений (4) - (7) аппроксимируется выражением

P (x, t ∣ x 0) ≈ 1 c N Σ {p} e - Σ j = - MM (xj - x 0, j (p)) 2/4 t D j. {\ Displaystyle P (х, т \ середина х_ {0}) \ приблизительно {\ гидроразрыва {1} {c_ {N}}} \ Sigma _ {\ {p \}} e ^ {- \ Sigma _ {j = -M} ^ {M} (x_ {j} -x_ {0, j} (p)) ^ {2} / 4tD_ {j}}.}{\ displaystyle P (x, t \ mid x_ {0}) \ приблизительно {\ frac {1} {c_ {N}}} \ Sigma _ {\ {p \}} e ^ { - \ Sigma _ {j = -M} ^ {M} (x_ {j} -x_ {0, j} (p)) ^ {2} / 4tD_ {j}}.}

(13)

Начиная с уравнения. (13) следует PDF-файл помеченной частицы в гетерогенном файле:

P (r, t ∣ r 0) ∼ 1 c N e - R d 2/4 τ τ (1 - a)) / (2 - γ (1 + a)); τ = Δ - 2 Λ t. {\ displaystyle P (r, t \ mid r_ {0}) \ sim {\ frac {1} {c_ {N}}} e ^ {- R_ {d} ^ {2} / 4 \ tau \ tau ^ { (1-a)) / (2- \ gamma (1 + a))}}; \ qquad \ tau = \ Delta ^ {- 2} \ Lambda t.}{\ displaystyle P (r, t \ mid r_ {0}) \ sim {\ frac {1} {c_ {N}}} e ^ {- R_ {d } ^ {2} / 4 \ tau \ tau ^ {(1-a)) / (2- \ gamma (1 + a))}}; \ qquad \ tau = \ Delta ^ {- 2} \ Lambda t. }

(14)

СКО Помеченная частица в гетерогенном файле взята из уравнения. (14):

⟨R d 2⟩ = 2 τ (1 - γ) / (2 c - γ), c = 1 / ((1 + a)). {\ Displaystyle \ langle R_ {d} ^ {2} \ rangle = 2 \ tau ^ {(1- \ gamma) / (2c- \ gamma)}, \ qquad c = 1 / ((1 + a)). }{\ displaystyle \ langle R_ {d} ^ {2} \ rangle = 2 \ tau ^ {(1- \ gamma) / (2c- \ gamma)}, \ qquad c = 1 / ((1 + a)).}

(15)

Обновление аномальных разнородных файлов

Результаты обновления аномальных файлов просто выводятся из результатов броуновских файлов. Во-первых, PDF в уравнении. (8) записывается в терминах PDF, который решает несверточное уравнение, то есть уравнение броуновского файла; это соотношение осуществляется в пространстве Лапласа:

P ¯ (x, s ∣ x 0) = 1 k ¯ α (s) P ¯ nrml (x, s / k ¯ α (s) ∣ x 0). {\ displaystyle {\ bar {P}} (x, s \ mid x_ {0}) = {\ frac {1} {{\ bar {k}} _ {\ alpha} (s)}} {\ bar { P}} _ {\ text {nrml}} (x, s / {\ bar {k}} _ {\ alpha} (s) \ mid x_ {0}).}{\ displaystyle {\ bar {P}} (x, s \ mid x_ {0}) = {\ frac {1} {{\ bar {k}} _ { \ alpha} (s)}} {\ bar {P}} _ {\ text {nrml}} (x, s / {\ bar {k}} _ {\ alpha} (s) \ mid x_ {0}).}

(16)

( Нижний индекс nrml обозначает нормальную динамику.) Из уравнения. (16), легко связать MSD броуновских разнородных файлов и неоднородных файлов с аномальным обновлением,

⟨r ¯ 2 (s)⟩ = 1 k ¯ α (s) ⟨r ¯ 2 (s / k ¯ α (s))⟩ nrml. {\ displaystyle \ langle {\ bar {r}} ^ {2} (s) \ rangle = {\ frac {1} {{\ bar {k}} _ {\ alpha} (s)}} \ langle {\ bar {r}} ^ {2} (s / {\ bar {k}} _ {\ alpha} (s)) \ rangle _ {\ text {nrml}}.}{\ displaystyle \ langle {\ bar {r}} ^ {2} (s) \ rangle = {\ frac {1} {{\ bar {k}} _ {\ alpha} (s)}} \ langle {\ bar {r}} ^ {2} (s / {\ bar {k}} _ {\ alpha} (s)) \ rangle _ {\ text {nrml}}.}

(17)

От Уравнение (18), обнаруживается, что MSD файла с нормальной динамикой в ​​степени α {\ displaystyle \ alpha}\ alpha является MSD соответствующего файла с аномальным обновлением,

⟨R 2 (t)⟩ ∼ r 2 (t)⟩ nrml α. {\ displaystyle \ langle r ^ {2} (t) \ rangle \ sim \ langle r ^ {2} (t) \ rangle _ {\ text {nrml}} ^ {\ alpha}.}{\ displaystyle \ langle r ^ {2} (t) \ rangle \ sim \ langle r ^ {2 } (t) \ rangle _ {\ text {nrml}} ^ {\ alpha}.}

(19)

Аномальные файлы с независимыми частицами

Уравнение движения для аномальных файлов с независимыми частицами (9) очень сложно. Решения для таких файлов достигаются при выводе законов масштабирования и с помощью численного моделирования.

Законы масштабирования для аномальных файлов независимых частиц

Во-первых, мы записываем закон масштабирования для среднего абсолютного смещения (MAD) в файл обновления с постоянной плотностью:

⟨∣ r ∣⟩ ∼ ⟨∣ r ∣⟩ бесплатно / n. {\ displaystyle \ langle \ mid r \ mid \ rangle \ sim \ langle \ mid r \ mid \ rangle _ {\ text {free}} / n.}{\ displaystyle \ langle \ mid r \ mid \ rangle \ sim \ langle \ mid r \ mid \ rangle _ {\ text {free}} / n.}

(20)

Здесь n {\ displaystyle n}n - количество частиц в покрытой длине ⟨∣ r ∣⟩ {\ displaystyle \ langle \ mid r \ mid \ rangle}{\ displaystyle \ langle \ mid r \ mid \ rangle} , и ⟨∣ r ∣⟩ free {\ displaystyle \ langle \ mid r \ mid \ rangle _ {\ text {free}}}{\ displaystyle \ langle \ mid r \ mid \ rangle _ {\ text {free}}} - MAD свободной аномальной частицы, ⟨∣ р ∣⟩ бесплатно ∼ t α / 2) {\ displaystyle \ langle \ mid r \ mid \ rangle _ {\ text {free}} \ sim t ^ {\ alpha / 2})}{\ displaystyle \ langle \ mid r \ mid \ rangle _ {\ text {free}} \ sim t ^ {\ alpha / 2})} . В формуле. (20), n {\ displaystyle n}n входит в вычисления, поскольку все частицы на расстоянии ⟨∣ r ∣⟩ {\ displaystyle \ langle \ mid r \ mid \ rangle}{\ displaystyle \ langle \ mid r \ mid \ rangle} от помеченной частицы должна двигаться в том же направлении, чтобы помеченная частица достигла расстояния ⟨∣ r ∣⟩ {\ displaystyle \ langle \ mid r \ mid \ rangle}{\ displaystyle \ langle \ mid r \ mid \ rangle} из исходного положения. На основании уравнения. (20), мы пишем обобщенный закон масштабирования для аномальных файлов независимых частиц:

⟨| г | ⟩ ∼ ⟨∣ r ∣⟩ бесплатно n f (n); 0 < f ( n) < 1. {\displaystyle \langle |r|\rangle \sim {\frac {\langle \mid r\mid \rangle _{\text{free}}}{n}}f(n);\qquad 0{\ displaystyle \ langle | r | \ rangle \ sim {\ frac {\ langle \ mid r \ mid \ rangle _ {\ text {free}}} {n}} f (n); \ qquad 0 <f (n) <1.}

(21)

Первый член в правой части уравнения. (21) также появляется в файлах продления; тем не менее, член f (n) единственен. f (n) - это вероятность, которая объясняет тот факт, что для перемещения n аномальных независимых частиц в одном направлении, когда эти частицы действительно пытаются прыгнуть в одном направлении (выражается с помощью члена (⟨∣ r ∣⟩ free / n) {\ displaystyle \ langle \ mid r \ mid \ rangle _ {\ text {free}} / n)}{\ displaystyle \ langle \ mid r \ mid \ rangle _ {\ text {free}} / n)} ), частицы на периферии должны двигаться первыми, так что частицы в середине файла будет иметь свободное пространство для перемещения, что потребует более быстрого перехода для тех, кто находится на периферии. f (n) появляется, потому что нет типичной шкалы времени для скачка в аномальных файлах, а частицы независимы, и поэтому конкретная частица может стоять на месте в течение очень долгого времени, существенно ограничивая возможности прогресса для частиц вокруг него., в течение этого времени. Ясно, что 0 < f ( n) < 1 {\displaystyle 0{\ displaystyle 0 <f (n) <1} , где f (n) = 1 для файлов обновления, поскольку частицы прыгают вместе, но также и в файлах независимых частиц с α>1 {\ displaystyle \ alpha>1}\alpha>1 , так как в таких файлах есть типичная шкала времени для скачка, рассматриваемая как время для синхронизированного скачка. Мы вычисляем f (n) из числа конфигураций, в которых порядок времени скачка частиц допускает движение; то есть порядок, в котором более быстрые частицы всегда находятся к периферии. Для n частиц существует n! различных конфигураций, из которых одна конфигурация является оптимальной; поэтому (1 / n!) ≤ f (n) {\ displaystyle (1 / n!) \ leq f (n)}{ \ Displaystyle (1 / п!) \ Leq е (п)} . Тем не менее, хотя и не оптимально, распространение также возможно во многих других конфигурациях; когда m - количество движущихся частиц, тогда

f (n) ∼ (нм) (п - т)! 1 / п!, {\ disp Laystyle f (n) \ sim {\ dbinom {n} {m}} (n-m)! 1 / n !,}{\ displaystyle f (n) \ sim {\ dbinom {n} {m}} (nm)! 1 / n !,}

(22)

где (n m) (n - m)! {\ displaystyle {\ dbinom {n} {m}} (n-m)!}{\ displaystyle {\ dbinom {n} {m}} (нм)!} подсчитывает количество конфигураций, в которых m частиц вокруг помеченной частицы имеют оптимальный порядок перехода. Теперь, даже когда m ~ n / 2, f (n) ∼ e - n / 2 {\ displaystyle f (n) \ sim e ^ {- n / 2}}{\ displaystyle f (n) \ sim e ^ {- n / 2}} . Использование в формуле. (21), е (n) ∼ e - n / n 0 {\ displaystyle f (n) \ sim e ^ {- n / n_ {0}}}{\ displaystyle f (n) \ sim e ^ {-n / n_ {0}}} (n 0 {\ displaystyle n_ {0}}n_ {0} небольшое число больше 1), мы видим,

MSD ∼ (α n 0) 2 ln 2 ⁡ (t). {\ displaystyle \ mathrm {MSD} \ sim \ left ({\ frac {\ alpha} {n_ {0}}} \ right) ^ {2} \ ln ^ {2} (t).}{\ displaystyle \ mathrm {MSD} \ sim \ left ({\ frac {\ alpha} {n_ {0 }}} \ right) ^ {2} \ ln ^ {2} (t).}

(23)

(В уравнении (23) мы используем MSD ∼∣ MAD ∣ 2 {\ displaystyle MSD \ sim \ mid MAD \ mid ^ {2}}{\ displaystyle MSD \ sim \ mid MAD \ mid ^ {2}} .) Уравнение (23) показывает, что асимптотически частицы очень медленные в аномальных файлах независимых частиц.

Численные исследования аномальных файлов независимых частиц

Рис. 1 Траектории моделирования 501 аномальной независимой частицы с α = 0.9, 0.1 {\ displaystyle \ alpha = 0.9, 0.1}{\ displaystyle \ альфа = 0,9,0.1} (рекомендуется: открыть файл в новом окне)

При численных исследованиях видно, что аномальные файлы независимых частиц образуют кластеры. Это явление определяет динамический фазовый переход. В установившемся режиме процент частиц в кластере, ξ (α) {\ displaystyle \ xi (\ alpha)}{\ displaystyle \ xi (\ alpha)} , следует,

ξ (α) ≈ 1 - α 3 { \ displaystyle \ xi (\ alpha) \ приблизительно {\ sqrt {1- \ alpha ^ {3}}}}{\ displaystyle \ xi (\ alpha) \ приблизительно {\ sqrt {1- \ alpha ^ {3}}}}

(24)

На рисунке 1 мы показываем траектории от 9 частиц в файле из 501 частицы. (Рекомендуется открывать файл в новом окне). На верхних панелях показаны траектории для α = 0,9 {\ displaystyle \ alpha = 0.9}{\ displaystyle \ alpha = 0.9} , а на нижних панелях показаны траектории для α = 0,1 {\ displaystyle \ alpha = 0,1}\ alpha = 0,1 . Для каждого значения α {\ displaystyle \ alpha}\ alpha показаны траектории на ранних этапах моделирования (слева) и на всех этапах моделирования (справа). Панели демонстрируют явление кластеризации, когда траектории притягиваются друг к другу, а затем перемещаются в значительной степени вместе.

См. Также
Ссылки
Последняя правка сделана 2021-05-20 03:41:50
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте