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

Кинетическая и потенциальные энергии будут иметь тогда следую­ щий вид

т =

+ М < 2/2 } ;

(4.77)

 

п

 

V = К К ,1

- «п-1д|2 /2 +

К К а - «2-i.ll2 /2 +

(4.78)

п

 

N

 

+ ] Г > К

1 - г Ч 2 |2/2,

 

 

П

 

 

 

где М — масса маятника, К — жесткость продольных пружинок и /с — жесткость поперечных пружинок.

Внося (4.76) в (4.77)-(4.78), мы получим модельный гамильтони­ ан Ht

Ht = T + V =

= Y . M i 2 {Ф1Л +</^)2} / 2 + Х / Х } ^ 2

/ 2+

п

п j

 

+ £ w

2 (y?n,1 - ¥42)2 / 2,

 

(4.79)

где переменные (pnj — независимые.

Мы можем записать теперь уравнения движения, соответствующие гамильтониану (4.79),

1фп,1 К 1 2 (<рп + 1,1

2<Рп,1 “Ь ф п —l,l)

М (<Рп,1

\

1 ф п ,2 = К р (^„+1,2 -

2у?п,2 +

2) -

k l 2 (y>„)2 -

(4.80)

<рп > 1 ) .

Здесь I = Ml2

4.2.2.1.Общее решение модельных уравнений

Чтобы найти решения модельных уравнений (4.80), сделаем преоб­ разование от переменных (fnj (t) к переменным Qqi9 (£), которые обычно называют нормальными координатами. Для этой цели давайте предпо­ ложим, что угловые смещения <pnj (t) имеют временную зависимость

Тогда уравнения (4.80) приобретут следующий вид

^ п , 1 = - к ?

(ёп+ 1л - 2ёпл+ ёп - 1,1)+ ы 2

;

I ^ V n ,2 =

(^п+1,2 - 2^ п ,2 + ^ п —1,2) + Ы 2 (<Pn,2 “

^ n . l )

Система 2iV линейных уравнений (4.82) имеет нетривиальное реше­ ние, если

det \1п>26П}П,6м, - Ап'П'ш* | = 0.

(4.83)

Здесь ненулевые коэффициенты

равны

 

АП}П-1}1 = An,n^2 = 2Kl2 + kl2-

 

Ап,п\1,2 = AnjTl]2,i =

—kl2]

(4.84)

An,n+l\jJ An^n—i-jj = —Kl ,

 

а корни уравнения (4.83) представляют собой так называемые частоты нормальных мод.

Учитывая трансляционную симметрию рассматриваемых уравнений, полезно сделать следующую замену

<п,з = ( ^ j / jl/2) ехР (i(lR n) = (Щ/11/2) ехР {щпа) ,

(4.85)

где q = {ях\Яу\Яг} = {0;0; gr} — волновой вектор, а его значения лежат в первой зоне Бриллюэна.

Теперь уравнения на собственные значения для частот нормальных

колебаний можно записать как

 

Iw2ipj =

exp [iq (га' —га) а] <^-.

(4.86)

п '

j '

 

Принимая во внимание соотношения (4.84), мы можем переписать уравнения (4.86) следующим образом

ь?Трх = {[2Kl2(1 - cos да) + kl2] / / } ^ i —{kl2/I) Тр2]

w2^ = {kl2/I) + { [2K12 (1 —cos да) + kl2] /1 } щ .

Для каждого значения g мы находим два решения для w2