работа студента 417 группы Башевого К. В (стр. 1 из 2)
Московский Государственный Университет им. М.В. Ломоносова
Физический факультет
Кафедра Физики Полимеров и Кристаллов
Курсовая работа студента 417 группы Башевого К.В.
Компьютерное моделирование процесса ультраметрической диффузии
научные руководители:
доцент Иванов В.А.
д.ф.-м.н. Аветисов В.А.
Заведующий кафедрой:
профессор Хохлов А.Р.
Москва, 2006
1. Введение
1.1. Актуальность проблемы
Ультраметрическая диффузия есть модельное описание динамики системы, энергетические состояния которой представлены многомерной сильно пересеченной поверхностью с большим числом локальных минимумов и седловых точек. В современной литературе такие системы часто называют "сложными системами", а соответствующие энергетические поверхности – многомерными сильно пересеченными энергетическими ландшафтами.
Развитие теоретических и компьютерных методов исследования динамики на сильно пересеченных ландшафтах стало особенно актуальным после того, как язык энергетических ландшафтов стал использоваться в изучении динамических и релаксационных процессов в неупорядоченных конденсированных средах, макромолекулярных структурах и биополимерах.
Сложность и трудоемкость подобных исследований обусловлена многомерностью и чрезвычайной изрезанностью энергетического ландшафта, что затрудняет применение вычислительных средств и, в конечном итоге, приводит к не прозрачным результатам. Использование теоретических моделей, в частности, модели ульраметрической диффузии, дает только качественную оценку поведения сложной системы, однако такая оценка служит естественным ориентиром в детальных исследованиях с применением методов компьютерного моделирования.
К настоящему времени показано, что модели ультраметрической диффузии могут быть использованы при описании конформационной динамики белковой макромолекулы и могут оказаться полезными для объяснения некоторых экспериментально наблюдаемых явлений, в частности, кинетики связывания СО миоглобином и явлении аномальной спектральной диффузии в белках в широкой области температур 0.1-300 К.
1.2. Цель работы
Целью курсовой работы являлась проверка теоретических моделей ультраметрических случайных процессов двух типов с помощью компьютерного эксперимента. В моделях первого типа, "частица" случайно блуждает в ультраметрическом пространстве, а время является обычной вещественной переменной. Такой случайный процесс называется ультраметрической диффузией. В моделях второго типа, "частица" случайно блуждает в пространстве с обычной (евклидовой) метрикой, а ультраметричность отражена во временной последовательности, на которой происходят случайные изменения положения "частицы". Такая временная последовательность вводится с помощью временных интервалов между моментами изменения случайной величины, длительность и статистика появления которых непосредственно связана с ультраметрической диффузией.
1.3. Основы теоретической модели ультраметрической диффузии
Теоретическая модель, которая исследуется в настоящей работе, была разработана научной группой В.А.Аветисова и состоит из двух взаимосвязанных частей. Первая часть – это случайное блуждание на прямой с дискретно заданными положениями блуждающей частицы. От обычной одномерной диффузии наша модель отличается тем, что смещение частицы происходит в моменты времени (тоже дискретные) разделенные различными временными промежутками. Вторая часть данной модели представляет собой генератор тех самых моментов времени, в которые осуществляется случайное смещение блуждающей частицы. В эти моменты времени с равной вероятностью осуществляется случайное смещение блуждающей частицы на единичное расстояние в положительном или отрицательном направлении оси x.
Проводилось моделирование двух вариантов одномерного случайного блуждания с ультраметричной временной последовательностью: "модель с задержками" и "модель активных бассейнов". В модели с задержками, временные интервалы между двумя последовательными смещениями блуждающей частицы определялись временами ожидания и их вероятностным распределением взятыми для "прыжков" другой, виртуальной "частицы", блуждающей в ультраметрическом пространстве.
В модели активных бассейнов изменения положения блуждающей частицы происходили только в те моменты времени, когда виртуальная "частица", блуждающая на ультраметрической решетке, попадала в определенные, наперед "помеченные" узлы.
2. Обзор литературы
Понятие ультраметрического пространства было введено в начале XX-го века вместе с понятием метрического пространства. Дело в том, что наряду с вещественными числами, в математике были уже введены и, так называемые, р-адические числа. Это было сделано Куртом Хенселем (Kurt Hensel) в1897, и именно р-адические числа требовали введения специального ультраметрического пространства.
В физической литературе термин «ультраметричность» появился относительно недавно, в 1979 году (Parisi, 1980), в связи с описанием фазовых переходов в неупорядоченных конденсированных системах в терминах спонтанного нарушения симметрии. Важным толчком к развитию таких представлений стало исследование спиновых стекол, а открытие ультраметрической структуры состояний спинового стекла было неожиданностью. Спиновые стекла оказались наиболее богатой моделью для исследований в этой области. Эти системы обладают каскадом фазовых переходов (и, следовательно, иерархическим типом организации при низких температурах), который никогда ранее в физике не наблюдался.
Одна из самых современных моделей ультраметрической диффузии разработана и описана в статье В.А. Аветисова [1]. Она позволила наиболее полно исследовать класс точно решаемых задач ультраметрической диффузии и использовать такие задачи для описания реальной кинетики связывания CO миоглобином. Более ранние исследования задачи об ультраметрической диффузии были ограничены лишь несколькими частными моделями иерархического типа. Они проводились в середине восьмидесятых годов и опубликованы в работах Огиельского и Штейна [7], Паладина [8], Гроссмана [4]. Следует также отметить работу Киркпатрика [5], в которой был разработан, так называемый, "метод вынужденного спуска", позволявший достаточно быстро "отыскивать нужное состояние" на сложном энергетическом ландшафте. Все эти работы дали толчок также и в решении ряда проблем биологии и химии.
3. Основная часть
3.1. Постановка задачи
Наша основная задача состоит в следующем. Пусть задана иерархия активационных барьеров с заданной зависимостью высоты барьеров от уровня иерархии. Также, задаем масштабный параметр для высоты барьеров α и выбираем модель (с задержками или случайных бассейнов). Требуется промоделировать процесс случайного блуждания с временными задержками и определить следующие характеристики процесса: зависимости среднего квадрата смещения частицы от времени при статистике временных задержек заданной генератором ультраметрической диффузии выбранного типа, распределение вероятности преодоления барьера i-того уровня в ультраметрическом генераторе, распределение времен задержек между моментами смещения блуждающей частицы. Также, проверить правильность поведения модели по сравнению с предсказаниями теории и выявить основные закономерности для такого случайного процесса.
3.2. Теория моделирования ультраметрических процессов
Для реализации вышеупомянутой второй части нашей модели, прежде всего, был построен генератор ультраметрического случайного блуждания. Он задается ультраметрической решеткой с иерархически вложенными друг в друга бассейнами состояний и матрицей вероятностей переходов между соответствующими бассейнами. На рис.1 показан отрезок прямой, на котором точками показаны дискретные состояния, в которых может находиться система. Эти состояния представляют собой конечные узлы некоторого древообразного графа с параметром ветвления p (на рисунке 1 этот параметр выбран равным p=2). Чтобы попасть из одной состояния в другое, система должна «подняться» сначала до той точки ветвления, откуда возможен «спуск» в новое состояние, а потом спуститься именно в одно из состояний, достижимых после преодоления точки ветвления. Можно переформулировать описание такого процесса на языке энергетических барьеров. Под состояниями системы понимаются квазиравновесные конфигурации, отвечающие минимумам энергетического ландшафта сложной системы. Каждое состояние отделено от любого другого состояния некоторым энергетическим барьером. Барьеры представляют собой иерархическую систему. Обозначим уровень иерархии барьера g, его высоту Hγ, а максимальный уровень иерархии r=γmax (см. рис.1). Сопоставим каждому состоянию индекс i, i=1,…,pγ. Множество состояний, каждому из которых можно приписать такой индекс i, назовем бассейном Bγ. Поделим Bγ на p взаимно непересекающихся подмножеств (бассейнов) Bγ-1(a1), каждый из которых включает pγ-1 состояний,
. Вероятность переходов, связанных с преодолением барьера Hγ, обозначим через qγ. Поделим каждый бассейн Bγ-1(a1) на p меньших бассейнов Bγ-2 (a1a2), a2=1, …, p, включающих в себя pγ-2 состояний, 
. Величину активационных барьеров между различными бассейнами уровня γ – 1 примем равной H γ-1. При этом вероятность перехода между любыми двумя состояниями 
и 
из различных бассейнов Bγ-2 есть q γ -1 (зависит от H γ -1) для a1=a|1, и q γ (зависит от H γ) для a1≠a|1. Данная процедура деления бассейнов продолжается до уровня, на котором бассейны B1(a1a2…aγ) включают по одному состоянию.