0.5
0.4
0.3
0.2
0.1
■ \\\ 1C
■\W\ \ \ 20
; \ у \л\ 4С\ У/ 1
" 1А X. . . . 1 . . , Ï : : : I— ï
0.2
0.4
0.6
0.8
Td
Рис. 7. Карта синхронных режимов системы (8). Кривая /2 разделяет область синхронизации (слева) и область трансверсальной неустойчивости (справа)
лент и обратный каскад бифуркаций удвоения периода, заканчивающийся циклом периода 1. Кривая /2 на рис. 7 ограничивает область полной синхронизации. Вне этой области эти режимы все еще существуют, но уже неустойчивы.
Заключение
Представленные в данной работе результаты свидетельствуют о том, что введение задержки в канале связи дискретной хаотической системы приводит к потере полной синхронизации хаоса [7]. Однако для систем с непрерывным временем существует интервал значений Т л, при которых синхронное хаотическое множество остается устойчивым. Более того, существует типичный сценарий, по которому синхронизация хаоса сменяется синхронизацией периодических ко-
лебаний через обратный каскад бифуркаций удвоения периода (см. рис. 7). При превышении величины Т d определенного порогового значения происходит потеря синхронизации периодических колебаний.
Библиографический список
1. Raddy D.V.R., Sen A., Johnston G.L. Time Delay Induced Death in Coupled Limit Cycle Oscillators // Phys. Rev. Lett. 1998. Vol. 80. P. 5109.
2. Raddy D.V.R., Sen A., Johnston G.L. Experimetnal Evidence of Time-Delay-Induced Death in Coupled Limit-Cycle Oscillators // Phys. Rev. Lett. 2000. Vol. 85. P. 3381.
3. Yeung M.K.S., Strogatz S.H. Time Delay in the Kuramoto Model of Coupled Oscillators // Phys. Rev. Lett. 1999. Vol. 82. P. 648.
4. Chung T.-H., Kim S. Spatio-Temporal Dynamics in Locally Coupled Ginzburg-Landau Oscillator Chain with Time Delay // Stochastic Dynamics and Pattern Formation in Biological and Complex Systems / Ed. by S. Kim, K.J. Lee, W. Sung. Berlin: Springer-Verlag, 2000. P. 67-74.
5. Astakhov V., Shabunin A., Klimshin A., Anishchenko V Inphase and Anti-phase Complete Chaotic Synchronization in Symmetrically Coupled Discrete Maps // Discrete Dynamics in Nature and Society. 2002. Vol. 7. P. 215.
6. Астахов В.В., Шабунин А.В., Анищенко В.С. Механизмы разрушения хаотической синхронизации в системе связанных кубических отображений // Изв. вузов. Прикладная нелинейная динамика. 1999. Т. 7, № 2. С. 3.
7. Астахов В.В., Неходцева Е.И., Астахов С.В., Шабунин А.В. Влияние задержки в канале связи на режимы полной синхронизации хаотических систем с дискретным временем // Изв. вузов: Прикладная нелинейная динамика. 2007. Т. 15, № 5. С. 61.
8. Anishchenko V. Dynamical Chaos: Models and Experiments. Singapore: World Scientific, 1995. 400 p.
УДК 539.194
СТОХАСТИЧЕСКАЯ ДИССОЦИАЦИЯ БРИЗЕРОВ УРАВНЕНИЯ СИНУС-ГОРДОНА ВО ВНЕШНИХ ПОЛЯХ
А.Г. Лавкин, М.Б. Мысенко, И.Н. Антонов
Саратовский государственный университет, кафедра теоретической и математической физики E-mail: TMF@SGU.RU
В работе исследуется методом максимальных ляпуновских показателей распад бризера Синус-Гордона в адиабатическом приближении во внешних полях. Определена зависимость между скоростью диссоциации бризера и значением максимального ля-пуновского показателя. В работе также рассмотрены другие численные методы исследования хаотизации поведения различных динамических систем.
Ключевые слова: бризер уравнения Синус-Гордона; стохастическая диссоциация; адиабатическое приближение; максимальный показатель Ляпунова; хаотичность.
The Sine-Gordon Brizer Stochastic Dissoziation in an External Fields A.G. Lavkin, M.B. Mysenko, I.N. Antonov
It is investigated, the Sine-Gordon brizer stochastic dissoziation in the external fields by the maximal Lyapunov parameters method in adiabatic approach. Dependence between the Sine-Gordon brizer dissoziation speed and the maximal Lyapunov parameters values is established.
© А.Г. Лавкин, М.Б. Мысенко, И.Н. Антонов, 2008
It is discussed also, the efficiency of the various numerical methods of research of the determined chaos of dynamic systems.
Key words: Sine-Gordon brizer; stochastic dissoziation; adiabatic approach; maximal Lyapunov parameter; chaotic character.
Введение
Весьма широкие физические приложения солитонов уравнения Синус-Гордона (СГ)
Ф - Фхх + sin (Ф) = 0
(1)
хорошо известны. Это математическое моделирование самонаведенной прозрачности при распространении ультракоротких лазерных импульсов через двухуровневую систему, динамики флюксонов - квантов потока магнитного поля в джозефсоновской линии, динамики дислокаций в кристаллах и многие другие применения данной полевой модели.
Возмущенное уравнение СГ имеет вид [1]
Ф - Фхх + sin (Ф) = £g.
(2)
Параметры солитона СГ в адиабатическом и следующем приближениях могут быть определены с помощью метода обратной задачи рассеяния (МОЗР).
Хаотичность динамики систем из солитонов СГ во внешнем монохроматическом поле численно моделировалась в работе [2] в максимально упрощенной длинноволновой модели. Методом сечений Пуанкаре и частотным методом показано, что переход системы к хаосу с ростом амплитуды и частоты внешнего поля происходит через удвоения периода внешнего возмущения.
В работе [3] моделировался случай гармонического внешнего поля методами качественной теории дифференциальных уравнений. С помощью функции Мельникова показано, что диссоциации бризера СГ соответствует развитая хаотичность его динамики и трансверсальное пересечение сепаратрис фазового пространства системы.
Хаотичность диссоциации бризера СГ в монохроматическом однородном внешнем поле моделировалась в [4] методом перекрытия резонансов в рамках простой квазичастичной модели. Была численно показана стабилизирующая роль диссипативных процессов (потери энергии бризером на излучение и т.п.) и возможность диссоциации бризера СГ при параметрическом его возбуждении.
В нашей работе рассмотрена динамика бризе-ра СГ под действием неоднородных возмущений, моделирующих, например, влияние на бризер дефектов конденсированных сред. Показано, что бризер СГ может диссоциировать в случайных
полях, а шумовые возмущения солитонных ансамблей приводят к хаотическим резонансам и к потерям энергии на излучение.
Хаотизация динамики бризера СГ внешнем полем
Рассмотрим [5] хаотизацию динамики бризера СГ внешнем полем. Возмущенное уравнение СГ (2) в этом случае имеет вид
Ф« - Фхх + sin (Ф) = eR (Ф), (3)
где &R (Ф) и есть возмущение. Параметры бризер-ного решения уравнения (3)
ф( X, t ) = -4arctg
V
nch(z )
(4)
в адиабатической модели подчиняются следующей системе уравнений:
Y = V = -^I
4cos y 4cos y
sH(I3 - VI4tgу)
(2sin у )
e = 4H cos у - &4H(VI3ctgу +
2 2 +HI4cos у- I5)(4sin у cos у )
(5)
которая получается при решении уравнения (3) методом МОЗР [4]. Здесь точкой обозначена производная по времени, у=аг^ (у/щ) - характерная амплитуда СГ-бризера, V - его скорость (Н=1-¥2), сС6/Ж - частота собственных пульсаций бризера,
= ф + nVZ, Z =
v(X___X°( )) , xo(t) - координа-
v Jn + iv|H ’
та центра бризера. Интегралы I1-I¡ записываются следующим образом (ß=ch2Z+tg2ycos2ф):
+œ ch(z)sin фпґ-, \\ ,
Il = J „ Д(ф( z))dz,
-œ
B
r +œ ch( z)sin фп..,
I2 = J 1 ^(ф( z ))dz,
-œ
B
13 zchMîi!^w z „dz,
B
-œ
14 =« m ( z ))dz,
B
(6)
I5 = 7 w z ))dz.
-œ
B
Рассмотрим динамику бризера СГ в случае, когда возмущающий член в уравнении (3) не зависит от пространственной координаты и имеет вид
єЯ(Ф) = - F (t) - аФ t, (7)
где F (t) - внешнее пространственно-однородное
2
nF sin 0 tgу sin 0
поле, а а характеризует интенсивность диссипативных процессов (потери энергии бризером на излучение и т.п.) [5]. В этом случае интегралы 12 и /3тождественно равны нулю, а V (^=соп^. Иными словами, возмущение (7) не выводит бризер из состояния покоя. Система уравнений (5) при ^0)=0 существенно упрощается (М=1+tg2y соб2#):
Y
- а-
cos y
M
ґл ctgYarsh(tgY cos 0Y (1 + ,—----------),
VM cos 0
(8)
0 = cosy -
nF
4sin y
cos 0
4M
cos2 y
- ctg(y)arsh(tgy cos 0)
cos 0 sin 0 tg (y )arsh(tg y cos 0) . 3^
- а----- ---а^^—¡=^--------- sin3 0.
Особые точки и предельные циклы фазовых траекторий (8) на (у, в) - плоскости численно исследовались в работе [6] в случае ступенчатой зависимости Е ($ внешнего поля от времени при мгновенном и адиабатическом его включении. Случай гармонического внешнего поля Е (^ = =А-Бт (ш$ изучался в [3] с помощью качественной теории дифференциальных уравнений. В данной работе динамика бризера СГ в ступенчатом и гармоническом внешних полях численно моделируется в адиабатической модели на основе уравнений (8) методом максимальных ляпуновских показателей. Кроме того, показывается, что величиной максимального ляпуновского показателя характеризуется скорость развала (диссоциации) бризера СГ во внешнем поле. Отметим еще раз, что диссоциация бризера СГ имеет место и при параметрическом его возбуждении [4].
Определенные заключения о динамике бризера СГ под действием возмущения (7) в адиабатической модели можно сделать на основе упрощенного квазичастичного подхода [7], рассматривая бризер как ангармонический осциллятор с пои (г) = -8 Ы "
тенциалом U (г)--ое -2F{t)г + а|Р0>тн|, где
r - расстояние между солитоном и антисолитоном,
У~
скорость их относительного движения.
Энергия бризера СГ отрицательна (Е= =-2tg2y), и в отсутствии внешнего поля солитон-антисолитон образуют связанную систему,
4Ё
и амплитудои
осциллирующую с частотой
Г 8 1
~1п -—- . Внешнее поле, превышающее крити-
и щ)
ческое значение Екр, определяемое уравнением
E = - F.
ln-
nFK
16
-1
, вызывает диссоциацию
бризера. Например, для Е ~ -0.16 FKp ~ 0.015, что
M cos Y
д/M
согласуется с результатами непосредственного решения уравнения Синус-Гордона (3).
Диссипативные процессы (а Ф 0) увеличивают Fкр и тем самым стабилизируют бризер от развала внешним полем. При достаточно интенсивной диссипации бризер СГ коллапсирует, не успев диссоциировать. Это типичная ситуация - стабилизирующий и дестабилизируюший факторы могут меняться ролями.
Моделирования нерегулярной динамики бризеров
Рассмотрим динамику бризера СГ, когда внешнее поле F (t) зависит от времени ступенчатым образом, используя метод максимальных ляпуновских показателей. При этом, как и в [6], ограничимся начальными условиями:
ctg2y (0) = 0.0785, 0 (0) = ^K +1 jrc , K = 1. (9)
Соответственно бризер имеет начальную энергию Е ~ -0.16 и находится в высоковозбужденном состоянии.
На рис. 1 представлена временная зависимость инкремента неустойчивости lmax (t) бризера СГ при мгновенном и адиабатическом включении в момент времени t=0 внешнего поля при а=0. Видно, что для одной и той же амплитуды F (t) максимальные ляпуновские показатели при мгновенном и адиабатическом включении внешнего поля совпадают. Это свидетельствует о равенстве F' в обоих случаях и согласуется с данными работы [6] при целых К в уравнении (9). Вследствие чего ниже моделируется только случай мгновенного включения внешнего поля.
На рис. 2 показана зависимость максимального ляпуновского показателя от амплитуды ступенчатого внешнего поля в отсутствии и при наличии диссипации (а=0.005). Видно, что
(оо ^ ;)
Рис. 1. Зависимость Х ((): а - для ступенчатого внешнего поля при а=0 (1, 2 - мгновенное, 3 -адиабатическое включение внешнего поля с амплитудой, равной 0.0136; 0.005; 0.005 соответственно, 4 - Е(0 при адиабатическом включении внешнего поля), б - для гармонического внешнего поля Е(0=Л-вт(а0, а=0.005 (1 - А=0.0165, ю=0.15; 2 - А=0.02, ю=0.5; 3 - А=0.03, ю=1)
Рис. 2. Зависимость X (^ю) от амплитуды внешнего поля: • - а=0, о - а=0.005
максимальные ляпуновские показатели положительны и динамика бризера СГ хаотична. Причем переход к развитой хаотичности, при которой Хтса (?—>^)~0.5 и скорость диссоциации бризера высока, происходит в отсутствие диссипации при Гкр ~ 0.015. При наличии диссипации (а=0.005) соответствующее значение Гкр ~ 0.02 заметно больше. Эти выводы находятся в согласии с данными работ [4, 6].
Из рис. 2 также следует, что динамика бризера СГ хаотична и при Г < Гкр. Здесь максимальные ляпуновские показатели определенно положительны, но по величине значительно меньше, чем при Г > Гкр. Последнее означает, что скорость диссоциации бризера чрезвычайно мала.
Учитывая регулярность динамики стабильных систем [7], для которых 1тса (—ю)=0, можно ожидать обратную зависимость времени диссоциации т бризера СГ от величины максимального ляпуновского показателя:
1
т*------------------. (ю)
|Лшах(^ ^ )Г
Скорость диссоциации бризера можно оценить [8, 9] по времени т, когда расстояние между солитоном и антисолитоном в момент развала
1, отн.еа.
бризера достигает определенной величины Ь. На рис. 3, а представлена зависимость т от амплитуды ступенчатого внешнего поля (Ь=2.75) в отсутствии и при наличии диссипации [10]. Аппроксимирующие кривые (10) на этом рисунке проведены с использованием данных, представленных на рис. 2. Видно, что удовлетворительное описание при а=0 достигается, если х=0.28, а при а=0.005, если х=0.64.
Теперь промоделируем с помощью метода максимальных ляпуновских показателей динамику бризера СГ в гармоническом внешнем поле Г (¿)=А-8т (Ш). Этот вопрос изучался в работе [3] методами качественной теории дифференциальных уравнений с использованием функции Мельникова и анализа трансверсальности пересечения сепаратрис фазовых траекторий.
На рис. 3, б плавной кривой показана граница между областью хаотичной динамики бризера СГ (выше кривой), когда он достаточно быстро диссоциирует, и областью его регулярной динамики (ниже кривой) [3]. На рис. 1, б изображена зависимость Хтах ($ при а=0.005 в трех точках граничной кривой, показанной на рис. 3, б: А1 = 0.0165, <^=0.15; А2=0.02, ш2=0.5; А3=0.03, ш3 = 1. Видно, что соответствующие
о
А/а
б
Рис. 3. Зависимость времени диссоциации бризера СГ от амплитуды внешнего поля: а - точки • - а=0, о - а=0.005, кривая 1 - т~А,хта1(г^сю) при х=-0.28, кривая 2 - т~Ххта1(г^да) при х=-0.64; б - граница между областями регулярной и хаотичной динамики бризера СГ в гармоническом внешнем поле А-8т(ш0 при наличии диссипации
значения максимального ляпуновского показателя положительны и малы по величине 0.06; 0.04; 0.13, причем в третьей точке его величина заметно отличается от первых двух. Это говорит о том, что на данной кривой динамика бризе-ра СГ слабохаотична, а выбранный в работе [3] критерий развала (степени хаотичности) бризера СГ в области частот ш ~1 не вполне эквивалентен критерию, использованному там при более низких частотах.
Заключение
В данной работе методом максимальных ляпуновских показателей численно исследована диссоциация бризера СГ во внешних однородных полях в отсутствие и при наличии диссипации. Показана хаотичность процесса диссоциации, установлена связь скорости диссоциации бризера СГ с величиной максимального ляпуновского показателя.
Библиографический список
1. Парментье Р. Солитоны в действии / Под ред. К. Лон-грена, Э. Скотта. М.: Мир, 1991. С. 185.
2. Yeh W., Kao Y. Universal Scaling and Chaotic Behavior of a Josephson-Junction Analog // Physical Review Letters. 1982. Vol. 49. P. 1888.
3. Ariyasu J., Bishop A. Determination of the onset of spatiotemporal complexity induced by soliton competition in a perturbed sine-Gordon chain // Physical Review. 1989. Vol. A39. P. 6409.
4. Абдуллаев Ф.Х. Динамический хаос солитонов. Ташкент: ФАН, 1990.
5. Лавкин А.Г. Исследование диссоциации бризера синус-Гордона во внешнем поле // Ядерная физика. 1991. Т. 54. С. 1290.
УДК 539.194
НЕРЕГУЛЯРНАЯ ДИНАМИКА ПОЛЕЙ ЯНГА-МИЛЛСА В ТЕРМОСТАТЕ
А.Г. Лавкин, М.Б. Мысенко, И.Н. Антонов
Саратовский государственный университет, кафедра теоретической и математической физики E-mail: TMF@sgu.ru
Исследуется хаотичность калибровочных полей в модели Янга-Миллса в формализме ланжевеновских источников в термостате. Методом максимальных ляпуновских показателей показано, что с ростом температуры возрастает и радиус корреляции (конфайн-мента), то есть имеет место фазовый переход адронная материя - кваркглюонная плазма. Причиной этого фазового перехода является стабилизация динамики полей глюонов флуктуациями термостата.
Ключевые слова: калибровочное поле, поле Янга-Миллса, конфайнмент, квантовая хромодинамика, кварк-глюонная плазма, метод источников Ланжевена.
The Irregular Dynamics of the Yang-Mills Fields in Thermostat A.G. Lavkin, M.B. Mysenko, I.N. Antonov
Dynamics of the homogeneous SU (2) Yang-Mills fields in thermostats and at their stochastic quantization are investigated by the Langeven's sources method.
It is shown, that with growth of temperature the phase transition confinement-deconfinement in chromodynamics takes place, and quantization of the Yang-Mills fields results in increase of radius of correlation (confinement) colors.
Key words: gage field, Yang-Mills field, confinement, quantum chromodynamics, quark-gluon plasma, Langeven's sources method.
Актуальность исследования хаотичности калибровочных полей в модели Янга-Миллса (ЯМ) связана с возможностью нетрадиционной интерпретации конфайнмента носителей цветных
6. Карпман В.И. Особые точки и предельные циклы фазовых траекторий бризеров уравнения Синус-Гордона // Журн. эксперим. и теорет. физики. 1983. Т. 84. С. 289.
7. BarashenkovI.V., Smirnov Yu.S., AlexeevaN.V Bifurcation to multisoliton complexes in the ac-driven, damped nonlinear Schrodinger equation // Physical Review. 1998. Vol. Е57. P. 2350.
8. Лавкин А.Г. Исследование диффузионной диссоциации классической гетерополярной молекулы в поле поляризованного ИК излучения // Опт. и спектр. 1991. Т. 70. С. 1021.
9. Матрасулов Д. У. Хаос в сверхкритическом атоме // Ядерная физика. 2001. Т. 64. С. 299.
10. Karpman V. Preprint IZMIRAN. № 52A, 1981.
зарядов (кварков и глюонов) в физике высоких энергий.
Наиболее доступны для исследования на хаотичность пространственно - однородные и плоские калибровочные поля ЯМ. Причем плоские волны ЯМ могут быть сведены к однородному случаю, так как они получаются из однородных с помощью преобразования Лоренца. Это преобразование всегда существует, в отличие, например, от линейной безмассовой электродинамики, именно благодаря нелинейности свободных неабелевых калибровочных полей ЯМ [1].
Несмотря на то, что эти модели полей ЯМ являются крайним упрощением для реальных пространственно-неоднородных полей, их исследования представляют значительный интерес вследствие следующих обстоятельств.
• Условие однородности означает длинноволновое приближение (или сильные поля) X >> >> А'1 (X, А - безразмерные длина волны и амплитуда поля) и можно надеяться, что изучение таких моделей полезно для получения информации об инфракрасном режиме квантовой хромодинамики (КХД) - наиболее неясном ее пункте.
• Однородные модели важны при рассмотрении проблемы квантования калибровочных полей в модели ЯМ. Поскольку хаотичность
© А.Г. Лавкин, М.Б. Мысенко, И.Н. Антонов, 2008