Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Нелинейная физика ДНК..pdf
Скачиваний:
9
Добавлен:
15.11.2022
Размер:
16.39 Mб
Скачать

первого приближения модель идеального газа и переписать корреляци­ онные функции (...), которые являются средними по ансамблю солитонов, в следующем виде

( - - - ) = N 3

(7.29)

где N s — среднее число частиц, а скобки (.. .)х

обозначают усреднение

по состояниям отдельного солитона

 

Li

J dZ J d P ( ...)exp ( - E 3/k BT)

(7.30)

dZ J dP exp ( - E s/kBT)

7.3.Задача рассеяния и нелинейные математические модели

Давайте обсудим теперь, каким образом подход, описанный в пре­ дыдущем разделе, может быть использован для решения задачи рас­ сеяния. В качестве примера рассмотрим, как и в разделе 4.3, задачу о рассеянии нейтронов. Хотя необходимо отметить, что результирующая формула для динамического фактора имеет общий храрактер и может быть применена для любой задачи рассеяния. В качестве примера мы покажем позднее, как эта формула может быть применена для задачи рассеяния света.

Начнем с формулы (4.106) и выражений для динамических факто­ ров когерентного и некогерентного рассеяния (4.107), которые являют­ ся достаточно общими и справедливыми как для линейного, так и для нелинейного случаев. Чтобы рассчитать упругое когерентное рассеяние, воспользуемся однофононным приближением, как мы делали это в раз­ деле 4.3. Тогда результирующую формулу для динамического фактора

можно записать как

 

S coh (x ,w ') = (2тrhN)

1е х р (- 2Щ * ) ^ ^ | exp [- х (Д° - Д°,)]

9оо

П П '

х J dtexp (iw't) (хщ, (t),

(0))}.

(7.31)

Формула (7.31) во многом похожа на формулу (4.121). Однако, мож­ но ожидать, что новые результаты вычислений корреляционных функций (xun(t), жтхп/(0)), которые мы предполагаем получить для нелинейного случая, будут отличаться от старых результатов, полученных в главе 4 для линейного случая. Естественно ожидать также, что эти новые ре­ зультаты будут различны для различных динамических моделей ДНК.

В этой главе мы представим результаты вычислений динамических факторов нейтронного рассеяния, полученные для трех нелинейных мо­ делей: для простой модели синус-Гордона, предложенной Инглэндером, для спиральной версии модели синус-Гордона и для У-модели.

7.3.1. Динамический фактор для простой модели синус-Гордона

Модель синус-Гордона, используемая для описания внутренней ди­ намики молекулы ДНК, была подробно изложена в разделе 5.1, а соот­ ветствующие динамические уравнения, записанные в безразмерном ви­ де, обсуждались еще раз в предыдущем разделе (см. формулу (7.14)).

Чтобы облегчить интерпретацию результатов, давайте вернемся к первоначальным переменным z^t. Тогда вместо уравнения (7.14) мы получим уравнение типа (6.1) (или (6.2)), а соответствующее солитонное решение в виде кинка будет иметь следующий вид

Ч>(z , t) = 4 tg-1 {exp [(7/d) (z - v t - г0)]},

(7.32)

где 7 = (1 v l / c y 1/2; d = al (K / v o f 2

Давайте проиллюстрируем, как вычисляется неупругая компонен­ та когерентного рассеяния, представляющая наибольший интерес в ис­ следованиях нелинейной динамики ДНК. Начнем с общего выражения для динамического фактора неупругого когерентного рассеяния, который определяется формулой (4.121). Эта формула может быть преобразована к более простому виду (7.31), так как в случае модели синус-Гордона индекс j можно опустить.

Чтобы вычислить динамический фактор, нам необходимо знать вы­ ражение для вектора и^, описывающего смещения оснований ДНК из положений равновесия. Это выражение легко получить, если учесть сходство между ДНК и механическим аналогом системы синус-Гордон. В этом случае смещения оснований эквивалентны смещениям маятников и имеют вид

где <рп (t ) — угловое смещение n-го маятника; I — его длина. А формула для динамического фактора неупругого когерентного рассеяния прини­ мает вид

 

5 coh (x,w') = (ж,г*/) + 5ц (ж, гг/),

(7.34)

где продольная и поперечная компоненты определяются формулами

5ц (ж, u/) =

(12X 1/2'X KKN) exp (—2WX) ЕЕ exp [—xza (n —n')] x

 

n

nf

 

 

 

4-00

 

 

 

x

J dt exp (iw't) ((1 - cos <pn (t )), (1 - cos tpn>(0)));

—OO

 

 

 

S_L (ж, гг/) =

(l2Xy/27rhN) exp (—2 W X) EE exp [—xza (n — n')] x

 

n

n'

 

 

 

4-oo

 

 

 

x

J dt exp (iw't) (sin <pn (t ), sin y?n/ (0))

 

—OO

 

 

(7.35)

 

 

 

 

В континуальном пределе формулы (7.35) преобразуются к виду

 

4-оо

4-оо

4-оо

 

5ц (*,«/) = (l2Xx/2nhN) exp(—2Wx) J dz j

dz' J

dtx

 

—oo

—OO

—OO

 

x{exp [—ixz (z —z)\ exp (iw't) ((1 —cos ip(z,t) ) , (1 —cosy? (г', 0)))};

 

4-oo

4-oo 4-oo

 

S±(x,u/) = (Pa*/2KhN)exp(-2Wx) J dz J

dz' J

dtx

 

—oo

—oo

—oo

 

x{exp[—ixz (z —г')] exp (iw't) (sirup (z^t) ,siny? (z\ 0))}.

(7.36) Теперь давайте воспользуемся приближением идеального газа,

в рамках которого предполагается, что

(7.37)

Здесь Ns — среднее число солитонов. Для величины Ns возьмем значе­ ние, полученное Куррие и соавторами [286]

Ж = {2Na/d) (Е0/2жкв Т ) ехр ( - Е0/кв Т )

(7.38)

Для

возьмем формулу

 

 

+оо

 

 

<•••>! = {2М0С0К 1(Ео/кв Т ) } - 1 J dpz (...)exp(Es/ k BT),

(7.39)

 

—оо

 

где

— функция МакДональда.

 

 

Подставляя (7.32) и (7.38) в (7.36) и учитывая формулы

(7.39)

и (7.34), найдем окончательное выражение для динамического факто­ ра неупругого когерентного рассеяния

Scoh (*,«>') = {2l2adl0 (E0/2kBT)1/2/hC0xzK,(E0/kBT) } e x p (- 2Wx) x

xexp (- E 0/kBT ) {x2 [(nxzd/-y0) / sh (■Kxzd/2'y0)}2 +

+x2y [(пхЛ/'уо) / Ch (■ Kxzd/2'yo)]2exp (-Eolo/keT)},

 

(1 - w,2/x 2zC$)

(7.40)

где 7o =

1/2

Результат (7.40) был получен для общего («релятивистского») слу­

чая. В

случае низких

температур (Т < Ео/кв) и малых скоростей

((wf/x z) С Со), формулу (7.40) можно привести к более простому виду

Scoh (x,wf) = {4l2adEo/hCoxzkBT } exp (—2WX) x

xexp ( - E 0/kBT ){ x l [(7rxzd/'yo) / s h {'Kxzd/2^o)f +

+ x2 [('KXzd/w) /ch (7rxzd/2jo)]2 } exp ( - M OW'Q/2х\кв Т)

(7.41) Заметим, что температура, соответствующая обычным физиологи­ ческим условиям, так же как и комнатная температура, принадлежат к так называемым «низким» температурам, то есть для них справедли­

во

условие Т

(Ео/кв ).

Действительно, если взять значение Ео =

=

6 ккал/моль,

даваемое

авторами работы [15], и значение кв =

=

1,38х10-23 джК-1 , то

(Ео/кв ) = ЗхЮ3 К. Последняя величи­

на

на порядок

превышает

комнатную температуру. Второе условие,

(w'/xz) «С Со, означает, что формула (7.41) справедлива для солитонных волн, распространяющихся в ДНК со скоростями, много меньшими скорости акустической волны Со (Со = 2х 103 м/сек [231]).

В заключение давайте покажем, что результаты (7.40) и (7.41) мо­ гут быть использованы и для других задач рассеяния. В качестве приме­ ра рассмотрим задачу рассеяния молекулой ДНК инфракрасного света.