Цепочка Тоды

Шаблон:Значения Цепо́чка То́ды (Шаблон:Lang-en) — система дискретных нелинейных уравнений, описывающих динамику взаимосвязанных нелинейных осцилляторов. Имеет важное значение в теории колебаний кристаллических решёток.

Система в общем случае имеет вид:

$m{\ddot{r}}_n=2f(r_n)-f(r_{n+1})-f(r_{n-1})$

где $r_n(t)$ имеет смысл величины отклонения n-го осциллятора от положения равновесия, а $f(r_i)$ — нелинейная функция, имеющая смысл возвращающей силы, действующей на i-ый осциллятор. Точки означают взятие операции дифференцирования.

Впервые предложена и проанализирована для случая $f(t)=-\alpha (1-\exp (-\beta t))$ Морикадзу Тодой в 1967 году.

Эквивалентная форма

Уравнение цепочки Тоды удобно анализировать в эквивалентной форме следующего вида

${\dot{s}}_n=f(r_n)$

$m{\dot{r}}_n=2s_n-s_{n+1}-s_{n-1}$

Решения

Можно показать, что уравнения, описывающие динамику цепочки Тоды, имеют решения в виде стационарных бегущих волн, имеющих вид

$s_n=S(\theta )=S(\omega t-pn)$

где функция $S(\theta )$ в случае, если $f(r_i)=-\alpha (1-\exp (-\beta t))$, удовлетворяет уравнению

$\frac{m{\omega }^2}{\beta }\frac{S^{″}}{\alpha +\omega S^{\prime }}=S(\theta +p)+S(\theta -p)-S(\theta )$

Решение этого уравнения выражается через эллиптические функции Якоби:

$S(\theta )=bZ(2K\theta )$

где

$Z(\theta )=E(\theta )-\theta \frac{E(K)}{K}$ — дзета-функция Якоби, имеющая период 2K

$E(\zeta )=\underset{0}{\overset{\zeta }{\int }}{\mathrm{d}\mathrm{n}}^{2}zdz$

Здесь K — полный эллиптический интеграл первого рода. Связь коэффициентов b и $\omega$ с параметрами $\alpha$, $\beta$ и m достаточно сложна, однако упрощается в предельных случаях.

Функция $r_n$ находится из соотношения

$-\alpha (1-e^{-\beta r_n})={\dot{s}}_n=2K\omega b({\mathrm{d}\mathrm{n}}^{2}2K\theta -\frac{E}{K})$

Особым решением является уединённое локализованное решение солитонного типа. Оно может быть получено в пределе $K\to \mathrm{\infty }$, при одновременном выполнении условий:

$2K\omega \to \mathrm{\Omega }$

$2Kp\to P$

В этом случае эллиптические функции переходят в гиперболические, и решение принимает вид

$-\alpha (1-\beta e^{-r_n})=\frac{m{\mathrm{\Omega }}^2}{\beta {\cosh }^2[\mathrm{\Omega }t-Pn]}$

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

Литература

• Шаблон:Статья
• Шаблон:Книга