close

Вход

Забыли?

вход по аккаунту

?

Pq)-резонанcные численно точные мобильные бризеры в модели молекулы ДНК Пейрара-Бишопа.

код для вставкиСкачать
ФИЗИКА
408
УДК 539.2
(P/Q)-РЕЗОНАНCНЫЕ ЧИСЛЕННО ТОЧНЫЕ МОБИЛЬНЫЕ БРИЗЕРЫ
В МОДЕЛИ МОЛЕКУЛЫ ДНК ПЕЙРАРА-БИШОПА
© М. И. Фахретдинов, Ф. К. Закирьянов
Башкирский государственный университет
Россия, Республика Башкортостан, 450076 г. Уфа, ул. Заки Валиди, 32.
Тел./факс: +7 (347) 273 93 25.
Email: fmi106tf@gmail.com
Найдены решения в виде устойчивых (p/q)-резонанcных численно точных мобильных дискретных бризеров для модели ДНК Пейрара-Бишопа. Эти решения можно получить только
при превышении величиной параметра межчастичного взаимодействия некоторого порога.
При значениях величины параметра межчастичного взаимодействия, близких к этому порогу,
полученные численно точные мобильные бризеры теряют энергию при движении и останавливаются.
Ключевые слова: (p/q)-резонанcные численно точные мобильные бризеры, модель ДНК
Пейрара-Бишопа.
Введение
Дискретные бризеры (ДБ) – это колебания, локализованные в пространстве и периодические по времени [1]. Они возникают в дискретных нелинейных
системах, в частности в квазилинейных молекулярных
цепочках. Одним из интересных примеров такого рода
объектов является модель ДНК Пейрара-Бишопа (ПБмодель) [2] и ее модификации [3–6]. В этих моделях
ДБ предшествуют возникновению «денатурационных
глазков» (denaturation bubbles) – областей разделения
комплементарных цепочек в процессе плавления ДНК
[7–11].
При определенных условиях такие локализованные осцилляции могут двигаться – обычно их называют мобильными бризерами (МБ). Существует систематический метод поиска решения дискретных
нелинейных систем в виде приближенных МБ [12].
Приближенные мобильные дискретные бризеры имеют достаточно большое время жизни, но при движении они постепенно теряют энергию из-за излучения
фононов. Однако численные исследования показывают существование в нелинейных дискретных решетках так называемых численно точных мобильных дискретных бризеров, передвигающихся без потери энергии и изменения своей формы.
Численно точный мобильный бризер (ЧТМБ) это
локализованное в пространстве периодическое решение дискретной системы, точно повторяющее свой
профиль при сдвиге на p узлов через некоторый период T.

 y n  p t  T   y n t 


 y n  p t  T   y n t 
(1)
Впервые решения такого типа упоминаются в
работе [13]. Обри и Кретеньи получили ЧТМБ для
уравнения Клейна-Гордона с потенциалом Морзе [13].
Гомез-Гарденез с соавторами исследовали численно
точные МБ в решетках, описываемых дискретным
нелинейным уравнением Шредингера [14]. Йошимура
с соавторами получили ЧТМБ для уравнения ФПУ
[15]. Они рассматривали так называемые (p/q)резонанcные ЧТМБ. То есть ЧТМБ имеющие внутреннюю частоту ωb и повторяющие свой профиль при
сдвиге на p узлов через время T = q (2π / ωb), где q –
натуральное число. Имеются также работы где T –
период ЧТМБ не задается, а вычисляется наряду с
профилем ЧТМБ методами минимизации (работа
группы Чечина по ЧТМБ в уравнении ФПУ [16]). В
дальнейшем под ЧТМБ мы будем иметь в виду (p/q)резонанcные ЧТМБ.
Работ по ЧТМБ пока очень мало. Сам алгоритм
получения ЧТМБ нигде формально не описан и есть
только намеки на его реализацию авторами статей
[13–15]. В них говорится, что для получения ЧТМБ
использовался метод Ньютона. Однако подчеркивается, что метод Ньютона, использующийся повсеместно
для получения стационарных ДБ, использовать нельзя, так как матрица Якоби оказывается сингулярной.
Но можно использовать метод Ньютона с изменениями, включающими сдвиг по координате и скорости на
p узлов и вместо решения СЛАУ JΔx = B для нахождения очередного приращения к решению в методе
Ньютона x = x + Δx, используется псевдообратная
(pseudo-inverse) матрица Якоби J–1, то есть Δx = BJ–1.
Псевдообратная матрица получается методом сингулярного разложения (SVD, Singular Value Decomposition) [17]. Также данный метод очень чувствителен к
начальным условиям и сходится только при выборе
«хороших» начальных условий, которыми могут быть
приближенные МБ.
Метод расчета ЧТМБ
Основным уравнением модели ДНК ПейрараБишопа, описывающим растяжение водородных связей оснований ДНК является уравнение:
yn  S  yn 1  2 yn  yn 1  



 2 2 exp  2 yn 1  exp  2 yn
,
(2)
где S – безразмерный параметр, характеризующий
стекинг-взаимодействие молекулы ДНК [2].
Задача нахождения решения (2) в виде ЧТМБ
сводится к нахождению начальных условий
y n 0   n , y n 0   n , при которых выполняются
условия (1). Для этого построим функции F и G следующего вида:
Fn ξ, ν   y n  p T ; ξ, ν    n
Gn ξ, ν   y n  p T ; ξ, ν   v n
(3)
Если начальные условия соответствуют ЧТМБ,
то функции F и G равны нулю. Дадим приращения Δξ
и Δν такие, что:
Fn ξ  ξ, ν  ν   0
Gn ξ  ξ, ν  ν   0
(4)
Разложим функции (4) в ряд Тейлора, отбрасывая квадратичные члены и члены более высокого порядка:
Вестник Башкирского университета. 2015. Т. 20. №2
ISSN 1998-4812
Fn ξ  ξ; ν  ν  
Fn
F
ξ  n ν
ξ
ν
Gn ξ  ξ; ν  ν  
 Fn ξ; ν  
 Gn ξ; ν  
(5)
Gn
Gn
ξ 
ν
ξ
ν
Подставляя (3) в (5) и приравнивая получившееся выражение нулю получим:
 yn  p T ; ξ, ν 


 I ξ 
ξ


yn  p T ; ξ, ν 

ν  ξ n  yn  p
ν
y n  p T ; ξ, ν 
ξ 
ξ
,
(6)
где I – единичная матрица порядка N, где N равно
числу узлов бризера.
Выражение (6) это СЛАУ для приращений Δν,
Δξ. Решить СЛАУ можно методом наименьших квадратов, что эквивалентно поиску решения методом
нахождения псевдообратной матрицы, как описано
выше. Тогда очередное приближение к решению
находится по формуле (7):
v n 1  v n  v n
пользуется как начальное условие для интегрирования
уравнений движения. МБ помещается в узел с номером k (значение k бралось как середина цепочки yn).
Исследуемая система интегрируется до времени
T = q Tb и считается норма разности между решением
и приближенным МБ, сдвинутым на p узлов вправо.
Этот процесс повторяется с постепенным увеличением множителя μ на небольшую величину до тех пор,
пока норма разности не достигнет минимального значения. Приближенный МБ, полученный при этом значении μ и принимается за начальное условие для
нахождения ЧТМБ.
Результаты и их обсуждение
 y n  p T ; ξ, ν 

 
 I ν  vn  y n  p

ν


ξ n 1  ξ n  ξ n
409
На рис. 1 показан профиль ЧТМБ в начальный
момент времени. Скорость ЧТМБ равна υb = 1/10 узлов, константа межчастичного взаимодействия S = 0.6.
Характерной особенностью ЧТМБ является
наличие ненулевых малоамплитудных осцилляций
распространяющихся по всех цепочке, названных
«хвостами» (tails) [13] (см. рис. 2). Хвосты зависят от
скорости ЧТМБ и параметров модели и имеют амплитуду от почти нулевой до величины порядка амплитуды ЧТМБ.
(7)
Собственно нахождение ЧТМБ делится на 2 этапа. Первый – это нахождение хорошего начального
приближения к методу Ньютона, а второй – использование метода Ньютона для нахождения ЧТМБ c
найденными на первом этапе приближенными МБ в
качестве начального приближения.
Нахождение начального приближения. Для получения хорошего начального приближения к методу
получения ЧТМБ вычисляется стационарный ДБ
ynst 0 c частотой ωb для выбранного значения константы межчастичного взаимодействия S. Cкоростная
часть ДБ задается нормализованным локализованным
собственным вектором пиннинг-моды или дискретным градиентом y n 0    y nst1 0  y nst1 0 2 .
Для нахождения приближения к ЧТМБ со скоростью υb = p/(q Tb) скоростная часть ДБ умножается на
множитель μ и получившийся приближенный МБ ис-
Рис. 2. Хвосты ЧТМБ, υb = 1/10, S = 0.6.
Найденные решения в виде ЧТМБ передвигаются по цепочке без изменения своей формы, и для них
строго выполняется условие (1). На рис. 3 изображено
прохождение ЧТМБ со скоростью υb = 1/10 по цепочке с периодическими граничными условиями. Показано смещение всех узлов цепочки. ЧТМБ проходит 500
периодов (500T = 17453.3) без возмущения окружения.
ЧТМБ можно получить только при определенных значениях параметра межчастичного взаимодействия S. При S < 0.44 ДБ становятся сильно локализо-
Рис. 1. Профиль ЧТМБ, в начальный момент времени координатная а) и скоростная часть б) υb = 1/10, S = 0.6.
ФИЗИКА
410
Рис. 4. Хвосты ЧТМБ, υb = 1/15, S = 0.44.
ванными, теряют мобильность и их трудно привести в
движение методом возмущений. Приближенные МБ
полученные при S < 0.44 осциллируют около центра
своего запуска. Метод Ньютона поиска ЧТМБ не сходится.
При 0.44 < S < 0.5 сходимость метода Ньютона
поиска ЧТМБ плохая. При этом часто полученные
ЧТМБ имеют хвосты большой амплитуды (см. рис. 4).
Они теряют энергию при движении, не могут преодолеть барьер Пайерлса-Набарро [1] и останавливаются
(см. рис. 5).
Выводы
Рассмотрен метод Ньютона нахождения (p/q)резонанcных ЧТМБ и поиска начального приближения к нему. Получены решения в виде (p/q)резонантных ЧТМБ в модели ДНК Пейрара-Бишопа.
ЧТМБ можно получить только при определенных значениях параметра межчастичного взаимодействия S.
При малых S ЧТМБ получить не удается. При
0.44 < S < 0.5 сходимость метода Ньютона поиска
ЧТМБ плохая, часто полученные ЧТМБ теряют энергию при движении, не могут преодолеть барьер Пайерлса-Набарро и останавливаются. При больших значениях S получение ЧТМБ для большого интервала
скоростей не представляет трудности.
Рис. 5. Временная эволюция ЧТМБ, υb = 1/15, S = 0.44.
ЛИТЕРАТУРА
1.
2.
3.
4.
5.
6.
7.
8.
Рис. 3. Временная эволюция ЧТМБ, υb = 1/10, S = 0.6.
9.
Flach S., Gorbach A. V. Discrete breathers – advances in
theory and applications // Physics Reports. 2008. Vol. 467.
No. 1. P. 1–116.
Peyrard M., Bishop A. R. Statistical mechanics of a nonlinear
model for DNA denaturation // Physical review letters. 1989.
Vol. 62. No. 23. P. 2755–2758.
Dauxois T. Dynamics of breather modes in a nonlinear “helicoidal” model of DNA // Physics Letters A. 1991. Vol. 159.
No. 8. P. 390-395.
Forinash K., Cretegny T., Peyrard M. Local modes and localization in a multicomponent nonlinear lattice // Physical Review E. 1997. Vol. 55. No. 4. P. 4740–4756.
Zhravković S., Satarić M. V. The impact of viscosity on the
DNA dynamics // Physica Scripta. 2001. Vol. 64. No. 6. P.
612–619.
Фахретдинов М. И., Закирьянов Ф. К. Конформационный
переход в модифицированной модели ДНК ПейярдаБишопа с трением // Вестник Башкирского университета.
2010. Т. 15. №1. С. 7–9.
Cuevas J. et al. Moving breathers in bent DNA with realistic
parameters // Modern Physics Letters B. 2004. Vol. 18. No.
25. P. 1319-1326.
Peyrard M., Cuesta-Lopez S., James G. Nonlinear analysis of
the dynamics of DNA breathing // Journal of biological physics. 2009. Vol. 35. No. 1. P. 73–89.
Maniadis P. et al. Feigenbaum cascade of discrete breathers in
a model of DNA // Physical Review E. 2011. Vol. 83. No. 1.
P. 011904.
ISSN 1998-4812
10.
11.
12.
13.
Вестник Башкирского университета. 2015. Т. 20. №2
Фахретдинов М. И., Закирьянов Ф. К. Дискретные бризеры в модели ДНК Пейрара-Бишопа // Журнал технической
физики. 2013. Т. 83. №7. C. 1–5.
Фахретдинов М. И., Закирьянов Ф. К., Екомасов Е. Г. Дискретные бризеры и мультибризеры в модели ДНК ПейрараБишопа // Нелинейная динамика. 2015. Т. 11. №1. С. 77–87.
Chen D., Aubry S., Tsironis G. P. Breather Mobility in Discrete φ4 Nonlinear Lattices // Physical review letters. 1996.
Vol. 77. No. 23. P. 4776–4779.
Aubry S., Cretegny T. Mobility and reactivity of discrete
breathers // Physica D: Nonlinear Phenomena. 1998. Vol. 119.
No. 1. P. 34–46.
14.
15.
16.
17.
411
Gómez-Gardenes J. et al. Nonintegrable Schrödinger discrete
breathers // Chaos: An Interdisciplinary Journal of Nonlinear
Science. 2004. Vol. 14. No. 4. P. 1130–1147.
Yoshimura K., Doi Y. Moving discrete breathers in nonlinear
lattice: Resonance and stability // Wave Motion. 2007. Vol.
45. No. 1. P. 83–99.
Лобзенко И. П., Чечин Г. М. Численное моделирование
движущихся дискретных бризеров в моноатомных цепочках // Вестник Нижегородского государственного университета. 2013. №4(1). С. 67–69.
Голуб Д., Ван Лоун Ч. Матричные вычисления. М.: Мир,
1999. 548 с.
Поступила в редакцию 28.05.2015 г.
ФИЗИКА
412
(P/Q)-RESONANT NUMERICALLY EXACT MOBILE BREATHERS
IN THE PEYRARD-BISHOP DNA MODEL
© F. M. Fakhretdinov*, F. K. Zakiryanov
Bashkir State University
32 Zali Validi St., 450076 Ufa, Republic of Bashkortostan, Russia.
Phone: +7 (347) 273 93 25.
*Email: fmi106tf@gmail.com
Discrete breathers are spatially localized time periodic oscillations. They appear in discrete nonlinear systems particularly in
the quasi-linear molecular chains. Interesting examples of such objects is the Peyrard-Bishop model of DNA and its modifications.
In the models, discrete breathers precede denaturation bubbles (the regions of separation of complementary chains in the course of
the DNA melting). Under certain conditions such localized oscillations can move and they are known as mobile breathers. There is a
systematic method for finding solutions of nonlinear discrete systems in the form of approximate mobile breathers. Approximate
mobile breathers have quite a long lifetime, but they move slowly losing energy due to the emission of phonons. However, in nonlinear discrete lattices numerical studies show the existence of so-called numerically exact mobile breathers moving without loss of
energy and not changing in shape. The present paper deals with these numerically exact mobile breathers in Peyrard-Bishop model of
DNA. A method of Newton for finding numerically exact mobile breather and finding the initial approximation to it is considered.
After qTb periods (Tb = 2π / ωb where ωb is internal frequency of mobile breather), these solutions repeat the same profile but displaced by p lattice sites. Numerically exact mobile breather can be obtained only at certain values of the interparticle interaction of S.
Keywords: (p/q)-resonant numerically exact mobile breather, Peyrard-Bishop DNA model.
Published in Russian. Do not hesitate to contact us at bulletin_bsu@mail.ru if you need translation of the article.
REFERENCES
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
Flach S., Gorbach A. V. Physics Reports. 2008. Vol. 467. No. 1. Pp. 1–116.
Peyrard M., Bishop A. R. Physical review letters. 1989. Vol. 62. No. 23. Pp. 2755–2758.
Dauxois T. Physics Letters A. 1991. Vol. 159. No. 8. Pp. 390-395.
Forinash K., Cretegny T., Peyrard M. Physical Review E. 1997. Vol. 55. No. 4. Pp. 4740–4756.
Zhravković S., Satarić M. V. Physica Scripta. 2001. Vol. 64. No. 6. Pp. 612–619.
Fakhretdinov M. I., Zakir'yanov F. K. Vestnik Bashkirskogo universiteta. 2010. Vol. 15. No. 1. Pp. 7–9.
Cuevas J. et al. Moving breathers in bent DNA with realistic parameters. Modern Physics Letters B. 2004. Vol. 18. No. 25. Pp. 13191326.
Peyrard M., Cuesta-Lopez S., James G. Journal of biological physics. 2009. Vol. 35. No. 1. Pp. 73–89.
Maniadis P. et al. Feigenbaum cascade of discrete breathers in a model of DNA. Physical Review E. 2011. Vol. 83. No. 1. Pp. 011904.
Fakhretdinov M. I., Zakir'yanov F. K. Zhurnal tekhnicheskoi fiziki. 2013. Vol. 83. No. 7. Pp. 1–5.
Fakhretdinov M. I., Zakir'yanov F. K., Ekomasov E. G. Nelineinaya dinamika. 2015. Vol. 11. No. 1. Pp. 77–87.
Chen D., Aubry S., Tsironis G. P. Physical review letters. 1996. Vol. 77. No. 23. Pp. 4776–4779.
Aubry S., Cretegny T. Physica D: Nonlinear Phenomena. 1998. Vol. 119. No. 1. Pp. 34–46.
Gómez-Gardenes J. et al. Nonintegrable Schrödinger discrete breathers. Chaos: An Interdisciplinary Journal of Nonlinear Science.
2004. Vol. 14. No. 4. Pp. 1130–1147.
Yoshimura K., Doi Y. Wave Motion. 2007. Vol. 45. No. 1. Pp. 83–99.
Lobzenko I. P., Chechin G. M. Vestnik Nizhegorodskogo gosudarstvennogo universiteta. 2013. No. 4(1). Pp. 67–69.
Golub D., Van Loun Ch. Matrichnye vychisleniya [Matrix calculations]. Moscow: Mir, 1999.
Received 28.05.2015.
Документ
Категория
Без категории
Просмотров
4
Размер файла
697 Кб
Теги
бишопа, пейрара, резонанcные, молекули, бризеры, точных, модель, днк, мобильный, численного
1/--страниц
Пожаловаться на содержимое документа