За достаточно большое время резонансы могут оказать сильное влияние на диффузию даже под действием слабого шума. Так, например, если в системе есть бсльшие (неперекрывающиеся) резонансы, то малое внешнее возмущение может перевести траекторию с нерезонансной инвариантной кривой на резонансную. Пусть время между внешними «толчками» велико по сравнению с периодом фазовых колебаний. Тогда следующий толчок может произойти уже
1) Пример подобной диффузии вследствие ошибок округления приведен в работе [73, §5]. Любопытно отмегить, что скорость диффузии оказалась на два порядка меньше, чем для стучайных ошибок той же величины.Прим. ред.

на другой стороне резонанса. Прн этом характерный коэффициент диффузии в окрестности резонанса значительно возрастет $\mathbf{1}$ ):
\[
D-(\Delta I)^{2 / \tau}
\]

где $\Delta I$ – ширина резонанса, а $\tau$ – время между двумя случайными толчками. Однако если значительные области фазового пространства не содержат больших резонансов, то, как будет видно, малая скорость нерезонансной дхффузии существенно подавляет глобальную диффузию.
Скорость диффузии. В присутствли слабого случайного возмущения скорость диффузии для стандартного отображения была получена аналитически Речестером и Уайтом [345] и Речестером и др. [346]. В последней работе был введен метод фурье-траекторий, описанный в предыдущем параграфе. ГІоажем, как следует видоизменить метод п. 5.4 г, чтобы учесть случайное возмущение. Введем в стандартное отображение случайное изменение фазы $\xi$ :
\[
\begin{array}{c}
I_{n+1}=I_{n}+K \sin \theta_{n}, \\
\theta_{n+1}=\theta_{n}+I_{n+1}+\xi_{n},
\end{array}
\]

где $\xi$ имеет гауссово распределение с дисперсией $\sigma$ :
\[
p(\xi)=\frac{1}{\sqrt{2 \pi \sigma}} \exp \left(-\frac{\xi^{2}}{2 \sigma}\right) .
\]

Тогда вероятность перехода (5.4.25) принимает вид
\[
\begin{array}{c}
W\left(I, \theta, n \mid I^{\prime}, \theta^{\prime}, n-1\right)=\int d \xi p(\xi) \delta\left(I-I^{\prime}-K \sin \theta^{\prime}\right) \times \\
\times \delta\left(\theta-\theta^{\prime}-I^{\prime}-K \sin \theta^{\prime}-\xi\right) .
\end{array}
\]

Подставив это выражение в (5.4.33) и проинтегрировав по $\theta$ и $I^{\prime}$, получим (5.4.34) с дополнительным множителем
\[
\int_{-\infty}^{\infty} e^{-i m \xi} p(\xi) d \xi=\exp \left(-\frac{1}{2} m^{2} \sigma\right) \text {. }
\]

В результате вместо (5.4.37б) приходим к следующему рекуррентному соотношению:
\[
a_{n}\left(m_{n}, q_{n}\right)=\sum_{l_{n}} \mathcal{g}_{l_{n}}\left(\left|q_{n-1}\right| K\right) \exp \left(-\frac{1}{2} m_{n}^{2} \sigma\right) a_{n-1}\left(m_{n-1}, q_{n-1}\right) .
\]
1) Нижеследующее выражение в тсчности совпадает с оценкой Будкера (см. [509], с. 50), который первым обратил внимание на такой механизм динамического усиления диффузии. Впоследствии это явление получило название неоклассической диффузии, тесрия которой была развита Галеевым и Сагдеевым [510].- Прим. ред.

Используя фурье-траектории, изображенные на рис. $5.12,6-e$, получаем коэффициент диффузии
\[
D_{n}=\frac{K^{2}}{2}\left[\frac{1}{2}-\mathcal{F}_{2}(K) e^{-\sigma}+\mathcal{F}_{2}^{2}(K) e^{-2 \sigma}+\mathscr{F}_{3}^{2}(K) e^{-3 \sigma}-\mathcal{J}_{1}^{2}(K) e^{-\sigma}\right],
\]

который обобщает соотношение (5.4.45), учитывая влияние внешнего шума при $K \gg 1$.

Рис. 5.16. Зависимость скорости диффузии $D / D_{1}$ от $K$ для различных условий внешнего шума о (по данным рабэты [346]).
Сплошные кривые – мстод фурье-траекторий; пунктирная кривая – оценка (5.5.14); кружки – численный счет $\left(\sigma=10^{-2}\right)$.

Более интересным является случай $K \ll 1$. Речестер и др. [346] показали, что в низшем порядке по $K$ вклад в коэффициент диффузии дают только фурье-траектории с $l=0$ и $l=1$, причем для

Рис. 5.17. Фазовая плоскость стандартного отображения.
Стрелками показана типичная траектория, триводящая к усилению внешней диффузии при $K \leqslant 1$.
$n \gg(2 / K)^{2}$ эти траектории должны начинаться и заканчиваться в точке $m=0, q=0$. В результате они получили
\[
\frac{D_{n}}{D_{1}}=\text { th }\left(\frac{\sigma}{2}\right),
\]

где $D_{1}=K^{2} / 4$ – квазилинейный коэффициент диффузии. Для определения поправок к $D_{n}$ было проведено численное суммирование по многим фурье-траекториям. На рис. 5.16 эти результаты (сплошные кривые) сравниваются с данными численного моделирования (кружки). Для малых $K$ значения $D$ сходятся к (5.5.9), а для больших $K$ величина $D$ не зависит от $\sigma$. Для нас представляет интерес область $K \leqslant 1$, где инвариантные кривые препятствуют развитию внутренней глобальной стохастичности.

Качественно эти результаты можно получить из следующих простых соображений (рис. 5.17). Прежде всего заметим, что изменение фазы $\theta$ на величину $\xi \ll 1$ приводит, как это следует из (5.5.3), к изменению действия
\[
\delta I=\left|I_{m} \cos \theta\right| \xi .
\]

Усреднение по распределению (5.5.4) и по $\theta$ дает
\[
\left\langle(\delta I)^{2}\right\rangle=\frac{1}{2} I_{m}^{2} \sigma .
\]

Принимая $I_{m} \approx K / 2$, получаем коэффициент диффузии
\[
D \approx\left\langle(\delta I)^{2}\right\rangle=\frac{K^{2}}{8} \sigma,
\]

что совпадает при малых $\sigma$ и $K$ с (5.5.9) ${ }^{1}$ ).
Вообще говоря, необходимо еще учесть усиление диффузии на резонансах. Пусть, например, траектория попадает под действием шума в стохастический слой целого резонанса в точке $A$, затем идет вдоль слоя $A \rightarrow B\left(B^{\prime}\right) \rightarrow C^{\prime}$ и покидает резонанс в точке $C^{\prime}$. Если полупериод фазовых колебаний $T$ (3.5.23) мал по сравнению с $1 / \sigma$ (что справедливо для малых $\sigma$ и $K \sim T \sim 1$ ), то область фазового пространства, в которой идет диффузия, сокращается на ширину резонанса. Заметим, что диффузия траектории $A \rightarrow C^{\prime}$ может идти либо наружу от резонанса, как рассмотрено выше, либо внутрь резонанса. В последнем случае средняя скорость диффузии падает. Таким образом, имеются две группы траекторий: быстрые, проходящие резонанс, и медленные, которые захватываются в резонанс. В качестве простой оценки примем, что отношение средних скоростей диффузии для двух групп траекторий обратно пропорционально квадрату интервала диффузии
\[
\frac{D_{f}}{D_{\mathrm{S}}}=\frac{(2 \pi)^{2}}{\left(2 \pi-2 K^{12}\right)^{2}},
\]

где, согласно (4.1.29), ширина целого резонанса на рис. 5.17 равна $2 K^{1 / 2}$, и мы пренебрегаем влиянием других резонансов. Примем далее, что половина траекторий, попадающих на границу резонанса (в точке $A$ ), совершает полпериода фазовых колебаний (от $A$ до $C^{\prime}$ ) и половина из них выходит из резонанса (в точке $C^{\prime}$ ). Таким образом, четверть всех траекторий проходит резонанс и попадает в группу быстрых траекторий. Эффективный коэффициент диффузии равен
\[
D_{\text {эфф }}=\frac{1}{4} D_{f}+\frac{3}{4} D_{s} .
\]

Используя для $D_{s}$ выражение (5.5.12), получаем с учетом (5.5.13) результат, показанный на рис. 5.16 пунктирной кривой. Несмотря на грубость оценки, она дает разумные значения даже для $K \geq 1$,
1) Приведенная оценка справедлиза лишь по порядку величины, и ее совпадение с (5.5.9) случайно [см. ниже (5.5.18) \}. – Прим. ред.

где внутренняя стохастичность становится глобальной. Связано это с тем, что если $K$ незначительно превышает критическое значение $K=1$, то средняя скорості диффузии по-прежнему определяется шумом, поскольку внутренняя диффузия является очень медленной. Другой более эффективный механизм усиления внешней диффузии резонансами будет рассмотрен в п. 6.3б.
Усреднение коэффициента диффузии. Согласно оценке (5.5.14), средняя скорость диффузии близка к медленной между резонансами, если только они не занимают значительную часть фазового пространства. Это связано с непрерывностью диффузионного потока, так что градиент плотности больше там, где локальная скорость диффузии меньше, и обратно. Рассмотрим некоторый стационарный поток $\Gamma$ в области $\Delta I$ с переменным коэффициентом диффузии. В каждом небольшом интервале с постоянной скоростью диффузии имеет место следующее соотношение для функции распределения $P$ :
\[
-\Gamma=D(I) \frac{d P}{d I} \equiv\langle D\rangle\left\langle\frac{d P}{d I}\right\rangle .
\]

Последнее равенство определяет средние по области значения $D$ и $d P / d I$. Таким образом,
\[
\frac{1}{\langle D\rangle} \frac{d P}{d I}=\frac{1}{D(I)}\left\langle\frac{d P}{d I}\right\rangle
\]

и усреднение по интервалу $\Delta I$ дает
\[
\frac{1}{\langle D\rangle} \frac{1}{\Delta I} \int \frac{d P}{d I} d I=\left\langle\frac{d P}{d I}\right\rangle \frac{1}{\Delta I} \int \frac{d I}{D(I)},
\]

или
\[
\frac{1}{\langle D\rangle}=\frac{1}{\Delta I} \int_{\Delta I} \frac{d I}{D(I)} .
\]

Это соотношение показывает, что $\langle D\rangle$ определяется в основном теми областями, в которых величина $1 / D$ максимальна, т. е. $D(I)$ минимальна.

Возвращаясь к оценке скорогти медленной диффузии (5.5.12), мы видим, что ее можно сделать более аккуратно, если найти зависимость $I_{m}(I)$ и усреднить по $I$ :
\[
\frac{1}{\left\langle D_{s}\right\rangle}=\frac{1}{2 \pi} \int_{0}^{2 \pi} \frac{d I}{\frac{1}{2} I_{m}^{2}(I) \sigma} .
\]

В заключение подчеркнем различие между усреднением диффузии по различным траекториям и усреднением по фазовому пространству. В первом случае усредняется сам коэффициент диффузии $D$, а во втором – его обратное значение $1 / D$.