Идеальная цепь

редактировать

Идеальная цепь (или свободно соединенная цепь ) - это простейшая модель для описывают полимеры, такие как нуклеиновые кислоты и белки. Он предполагает только полимер как случайное блуждание и не учитывает любые взаимодействия между мономерами. Хотя он прост, его общность дает представление о физике полимеров.

В этой модели мономеры представляют собой жесткие стержни фиксированной длины l, и их ориентация полностью не зависит от ориентации и положения соседних мономеров в той степени, в которой два мономера могут сосуществовать в одном месте. В некоторых случаях мономер имеет физическую интерпретацию, например, аминокислота в полипептиде. В других случаях мономер - это просто сегмент полимера, который можно моделировать как дискретную, свободно соединенную единицу. Если это так, l - это длина Куна. Например, хроматин моделируется как полимер, в котором каждый мономер представляет собой сегмент длиной примерно 14-46 т.п.н.

Содержание
  • 1 Модель
  • 2 Общность модели
  • 3 Энтропийная эластичность идеальной цепи
    • 3.1 Идеальная цепь при ограничении длины
    • 3.2 Идеальная длина обмена полимера с резервуаром
    • 3.3 Идеальная цепь при ограничении постоянной силы - расчет
  • 4 См. Также
  • 5 Внешние ссылки
  • 6 Ссылки
Модель

N мер образуют полимер, общая развернутая длина которого составляет:

L = N l {\ displaystyle L = N \, l}L = N \, l , где N - количество меров.

В этом очень простом подходе, где взаимодействия между мерами не рассматриваются, энергия полимера считается независимой от его формы, что означает, что при термодинамических равновесие, все его конфигурации формы с одинаковой вероятностью возникнут, поскольку полимер колеблется во времени, в соответствии с распределением Максвелла – Больцмана.

. Позвольте нам называть R → {\ displayst yle {\ vec {R}}}{\ vec {R}} полный сквозной вектор идеальной цепи и r → 1,…, r → N {\ displaystyle {\ vec {r}} _ { 1}, \ ldots, {\ vec {r}} _ {N}}{\ vec r} _ {1}, \ ldots, {\ vec r} _ {N} векторы, соответствующие отдельным мерам. Эти случайные векторы имеют компоненты в трех направлениях пространства. Большинство выражений, приведенных в этой статье, предполагают, что количество мер N велико, поэтому применяется центральная предельная теорема . На рисунке ниже показан эскиз (короткой) идеальной цепи.

Случайное блуждание по идеальной цепочке. Svg

Два конца цепочки не совпадают, но они колеблются вокруг друг друга, так что, конечно же:

⟨R →⟩ = Σ i = 1 N ⟨r → i⟩ = 0 → {\ displaystyle \ langle {\ vec {R}} \ rangle = \ Sigma _ {i = 1} ^ {N} \ langle {\ vec {r}} _ {i} \ rangle = {\ vec {0}} ~}\ langle {\ vec {R}} \ rangle = \ Sigma _ {{i = 1}} ^ { N} \ langle {\ vec r} _ {i} \ rangle = {\ vec 0} ~

На протяжении всей статьи скобки ⟨⟩ {\ displaystyle \ langle \ rangle}\ langle \ rangle будут использоваться для обозначения среднего (значений, взятых с течением времени) случайной величины или случайный вектор, как указано выше.

Поскольку r → 1,…, r → N {\ displaystyle {\ vec {r}} _ {1}, \ ldots, {\ vec {r}} _ {N}}{\ vec r} _ {1}, \ ldots, {\ vec r} _ {N} независимы, из Центральной предельной теоремы следует, что R → {\ displaystyle {\ vec {R}}}{\ vec {R}} распределяется согласно нормальному распределению (или гауссовскому распределению): точно, в 3D, R x, R y, {\ displaystyle R_ {x}, R_ {y},}R_ {x}, R_ {y}, и R z {\ displaystyle R_ {z}}R_ {z} распределяются согласно нормальному распределению из среднего 0 и дисперсия :

σ 2 знак равно ⟨р Икс 2⟩ - ⟨р Икс⟩ 2 = ⟨р Икс 2⟩ - 0 {\ displaystyle \ sigma ^ {2} = \ langle R_ {x} ^ {2} \ rangle - \ langle R_ {x} \ rangle ^ {2} = \ langle R_ {x} ^ {2} \ rangle -0}\ sigma ^ {2} = \ langle R_ {x} ^ {2} \ rangle - \ langle R_ {x} \ rangle ^ {2} = \ langle R_ {x} ^ {2} \ rangle -0
⟨R x 2⟩ = ⟨R y 2⟩ = ⟨R z 2⟩ = N l 2 3 {\ displaystyle \ langle R_ {x} ^ {2} \ rangle = \ langle R_ {y} ^ {2} \ rangle = \ langle R_ {z} ^ {2} \ rangle = N \, {\ frac {l ^ {2}} {3}}}\ langle R_ {x} ^ {2} \ rangle = \ langle R_ {y} ^ {2} \ rangle = \ langle R_ {z} ^ {2} \ rangle = N \, {\ frac {l ^ {2}} {3}}

Итак, ⟨R 2 →⟩ = N l 2 = L l {\ displaystyle \ langle {\ vec {R ^ {2}}} \ rangle = N \, l ^ {2} = L \, l ~}\ langle {\ vec {R ^ {2}}} \ rangle = N \, l ^ {2} = L \, l ~ . Сквозной вектор цепочки распределяется согласно следующей функции плотности вероятности :

P (R →) = (3 2 π N l 2) 3/2 e - 3 R → 2 2 N l 2 {\ Displaystyle P ({\ vec {R}}) = \ left ({\ frac {3} {2 \ pi Nl ^ {2}}} \ right) ^ {3/2} e ^ {- {\ frac {3 {\ vec {R}} ^ {2}} {2Nl ^ {2}}}}}P ({\ vec R}) = \ left ({\ frac {3} {2 \ pi Nl ^ {2}}} \ right) ^ {{3/2}} e ^ {{- {\ frac {3 {\ vec R} ^ {2}} {2Nl ^ {2}}}}}

Среднее расстояние от конца до конца полимера:

⟨R 2 →⟩ = N l = L l {\ displaystyle {\ sqrt {\ langle {\ vec {R ^ {2}}} \ rangle}} = {\ sqrt {N}} \, l = {\ sqrt {L \, l}} ~}{\ sqrt {\ langle {\ vec {R ^ {2}}} \ rangle}} = {\ sqrt N} \, l = {\ sqrt {L \, l}} ~

В физике полимеров часто используется величина радиус инерции :

RG = N l 6 {\ displaystyle {\ mathit {R}} _ {G} = {\ frac {{\ sqrt {N}} \, l} {{\ sqrt {6}} \}}}{\ mathit {R}} _ {G} = {\ frac { {\ sqrt N} \, l} {{\ sqrt 6} \}}

Следует отметить, что указанное выше среднее расстояние от конца до конца, которое в случае этой простой модели также является типичным амплитуда колебаний системы становится незначительной по сравнению с общей развернутой длиной полимера N l {\ displaystyle N \, l}N \, l при термодинамическом пределе. Этот результат является общим свойством статистических систем.

Математическое замечание: строгая демонстрация выражения плотности вероятности P (R →) {\ displaystyle P ({\ vec {R}})}P ({\ vec R}) не является так же прямо, как показано выше: из применения обычной (1D) центральной предельной теоремы можно вывести, что R x {\ displaystyle R_ {x}}R_ {x} , R y {\ displaystyle R_ {y}}R_ {y} и R z {\ displaystyle R_ {z}}R_ {z} распределяются согласно центрированному нормальному распределению дисперсии N l 2/3 {\ Displaystyle N \, l ^ {2} / 3}N \, l ^ {2} / 3 . Тогда приведенное выше выражение для P (R →) {\ displaystyle P ({\ vec {R}})}P ({\ vec R}) не единственное, которое совместимо с таким распределением для R x {\ displaystyle R_ {x}}R_ {x} , R y {\ displaystyle R_ {y}}R_ {y} и R z {\ displaystyle R_ {z}}R_ {z} . Однако, поскольку компоненты векторов r → 1,…, r → N {\ displaystyle {\ vec {r}} _ {1}, \ ldots, {\ vec {r}} _ {N}}{\ vec r} _ {1}, \ ldots, {\ vec r} _ {N} некоррелированы для рассматриваемого нами случайного блуждания, из этого следует, что R x {\ displaystyle R_ {x}}R_ {x} , R y {\ displaystyle R_ {y} }R_ {y} и R z {\ displaystyle R_ {z}}R_ {z} также некоррелированы. Это дополнительное условие может быть выполнено, только если R → {\ displaystyle {\ vec {R}}}{\ vec {R}} распределено согласно P (R →) {\ displaystyle P ({\ vec {R}})}P ({\ vec R}) . В качестве альтернативы, этот результат также может быть продемонстрирован путем применения многомерного обобщения центральной предельной теоремы или с помощью аргументов симметрии.

Общность модели

Хотя описанная выше элементарная модель полностью не адаптирована к описанию реальных полимеров в микроскопическом масштабе, она действительно демонстрирует некоторую актуальность в макроскопическом масштабе в случае полимера в растворе, мономеры которого образуют идеальную смесь с растворителем (в этом случае взаимодействия между мономером и мономером, молекулой растворителя и молекулой растворителя, а также между мономером и растворителем идентичны, и энергия системы может считаться постоянной, что подтверждает гипотезы модели).

Однако актуальность модели ограничена даже в макроскопическом масштабе тем фактом, что она не учитывает какой-либо исключенный объем для мономеров (или, говоря химическими терминами, пренебрегает) стерические эффекты ).

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

Энтропическая эластичность идеальной цепи

Если два свободных конца идеальной цепи прикреплены к какому-либо микроманипуляционному устройству, то это устройство испытывает силу, оказываемую полимером. Энергия идеальной цепи постоянна, и, следовательно, ее среднее по времени, внутренняя энергия, также постоянна, что означает, что эта сила обязательно проистекает из чисто энтропийного эффекта.

Эта энтропийная сила очень похожа на давление, испытываемое стенками коробки, содержащей идеальный газ. внутренняя энергия идеального газа зависит только от его температуры, а не от объема контейнера, в котором он находится, поэтому это не эффект энергии, который имеет тенденцию чтобы увеличить объем коробки, как это делает газ давление. Это означает, что давление идеального газа имеет чисто энтропийное происхождение.

Каково микроскопическое происхождение такой энтропийной силы или давления? Самый общий ответ состоит в том, что эффект тепловых флуктуаций имеет тенденцию приближать термодинамическую систему к макроскопическому состоянию, которое соответствует максимуму в количестве микроскопических состояний (или микросостояний), совместимых с этим макроскопическим состоянием. Другими словами, тепловые флуктуации приводят систему к макроскопическому состоянию максимальной энтропии.

. Что это означает в случае идеальной цепи? Во-первых, для нашей идеальной цепи микроскопическое состояние характеризуется суперпозицией состояний r → i {\ displaystyle {\ vec {r}} _ {i}}\ vec r_i каждого индивидуального мономера ( где i изменяется от 1 до N). В своем растворителе идеальная цепь постоянно подвергается ударам от движущихся молекул растворителя, и каждый из этих ударов переводит систему из ее текущего микроскопического состояния в другое, очень похожее микроскопическое состояние. Для идеального полимера, как будет показано ниже, существует больше микроскопических состояний, совместимых с коротким сквозным расстоянием, чем есть микроскопических состояний, совместимых с большим сквозным расстоянием. Таким образом, для идеальной цепи максимизация ее энтропии означает уменьшение расстояния между двумя ее свободными концами. Следовательно, сила, которая стремится разрушить цепь, действует идеальная цепь между двумя ее свободными концами.

В этом разделе будет выведено среднее значение этой силы. Затем будет обсуждаться общность выражения, полученного при термодинамическом пределе.

Идеальная цепь при ограничении длины

В этом подразделе будет рассмотрен случай идеальной цепи, два конца которой прикреплены к фиксированным точкам. Вектор R → {\ displaystyle {\ vec {R}}}{\ vec {R}} , соединяющий эти две точки, характеризует макроскопическое состояние (или макросостояние) идеальной цепи. Каждому макросостоянию соответствует определенное количество микросостояний, которые мы будем называть Ω (R →) {\ displaystyle \ Omega ({\ vec {R}})}\ Omega ({\ vec R}) (микросостояния определены во введении к этому разделу). Поскольку энергия идеальной цепи постоянна, каждое из этих микросостояний имеет равную вероятность возникновения. энтропия, связанная с макросостоянием, таким образом, равна:

S (R →) = k B log ⁡ (Ω (R →)) {\ displaystyle S ({\ vec {R} }) = k_ {B} \ log (\ Omega ({\ vec {R}}))}S ({\ vec R}) = k_ {B} \ log (\ Omega ({\ vec R})) , где k B {\ displaystyle k_ {B}}k_ {B} равно константа Больцмана

Вышеприведенное выражение дает абсолютную (квантовую) энтропию системы. Для точного определения Ω (R →) {\ displaystyle \ Omega ({\ vec {R}})}\ Omega ({\ vec R}) потребуется квантовая модель для идеальной цепи, что выходит за рамки данной статья. Однако мы уже вычислили плотность вероятности P (R →) {\ displaystyle P ({\ vec {R}})}P ({\ vec R}) , связанную со сквозным вектором безусловного идеала. цепь, выше. Поскольку вероятность возникновения всех микросостояний идеальной цепочки одинакова, P (R →) {\ displaystyle P ({\ vec {R}})}P ({\ vec R}) пропорционально Ω (R →) {\ Displaystyle \ Omega ({\ vec {R}})}\ Omega ({\ vec R}) . Это приводит к следующему выражению для классической (относительной) энтропии идеальной цепи:

S (R →) = k B log ⁡ (P (R →)) + C st {\ displaystyle S ({\ vec {R}}) = k_ {B} \ log (P ({\ vec {R}})) + C_ {st}}S ({\ vec R}) = k_ {B} \ журнал (P ({\ vec R})) + C _ {{st}} ,

где C st {\ displaystyle C_ {st }}C _ {{st}} - фиксированная константа. Назовем F → {\ displaystyle {\ vec {F}}}{\ vec {F}} силу, прилагаемую цепью к точке, к которой прикреплен ее конец. Из приведенного выше выражения энтропии мы можем вывести выражение этой силы. Предположим, что положение двух концов идеальной цепочки теперь не фиксируется, а контролируется оператором. Оператор управляет эволюцией сквозного вектора R → {\ displaystyle {\ vec {R}}}{\ vec {R}} . Если оператор меняет R → {\ displaystyle {\ vec {R}}}{\ vec {R}} на крошечную величину d R → {\ displaystyle {\ vec {dR}}}{\ vec {dR}} , тогда изменение внутренней энергии цепочки равно нулю, поскольку энергия цепочки постоянна. Это условие можно записать как:

0 = d U = δ W + δ Q {\ displaystyle 0 = dU = \ delta W + \ delta Q ~}0 = dU = \ delta W + \ delta Q ~

δ W {\ displaystyle \ delta W}\ delta W определяется как элементарное количество механической работы, переданное оператором идеальной цепочке, а δ Q {\ displaystyle \ delta Q}\ delta Q определяется как элементарное количество тепла, передаваемое растворителем идеальной цепи. Теперь, если мы предположим, что преобразование, наложенное оператором на систему, является квазистатическим (т. Е. Бесконечно медленным), то преобразование системы будет обратимым во времени, и мы можем предположить, что во время ее перехода из макросостояния R → {\ displaystyle {\ vec {R}}}{\ vec {R}} в макросостояние R → + d R → {\ displaystyle {\ vec {R}} + {\ vec {dR}}}{\ vec R} + {\ vec {dR}} , система проходит через серию термодинамического равновесия макросостояний. Это имеет два следствия:

  • во-первых, количество тепла, полученного системой во время преобразования, может быть связано с изменением ее энтропии :
δ Q = T d S {\ displaystyle \ delta Q = TdS ~}\ delta Q = TdS ~ , где T - температура цепи.
  • секунд, чтобы преобразование оставалось бесконечно медленным, среднее усилие, прилагаемое оператор на конечных точках цепи должен уравновесить среднее усилие, оказываемое цепью на ее конечные точки. Вызов f → op {\ displaystyle {\ vec {f}} _ {op}}{\ vec f} _ {{op}} силы со стороны оператора и f → {\ displaystyle {\ vec {f}} }\ vec f сила, действующая со стороны цепи, мы имеем:
δ W = ⟨f → op⟩ ⋅ d R → = - ⟨f →⟩ ⋅ d R → {\ displaystyle \ delta W = \ langle {\ vec {f}} _ {op} \ rangle \ cdot {\ vec {dR}} = - \ langle {\ vec {f}} \ rangle \ cdot {\ vec {dR}} ~}\ delta W = \ langle {\ vec f} _ {{op}} \ rangle \ cdot {\ vec {dR}} = - \ langle {\ vec f} \ rangle \ cdot {\ vec {dR}} ~

Таким образом, мы приходим к:

⟨е →⟩ = T d S d R → = k BTP (R →) d P (R →) d R → {\ displaystyle \ langle {\ vec {f}} \ rangle = T {\ frac {dS} {\ vec {dR}}} = {\ frac {k_ {B} T} {P ({\ vec {R}})}} {\ frac {dP ({\ vec { R}})} {\ vec {dR}}} ~}\ langle {\ vec f} \ rangle = T {\ frac {dS} {{\ vec {dR}}}} = {\ frac {k_ {B} T} {P ({\ vec R})}} {\ frac {dP ({\ vec R})} {{\ vec {dR}}}} ~
⟨е →⟩ = - k BT 3 R → N l 2 {\ displaystyle \ langle {\ vec {f}} \ rangle = -k_ { B} T {\ frac {3 {\ vec {R}}} {Nl ^ {2}}} ~}\ langle {\ vec f} \ rangle = -k_ {B} T {\ frac {3 {\ vec R}} {Nl ^ {2}}} ~

Вышеприведенное уравнение является уравнением состояния идеальной цепи. Поскольку выражение зависит от центральной предельной теоремы , оно является точным только в пределах полимеров, содержащих большое количество мономеров (то есть термодинамический предел ). Это также справедливо только для небольших расстояний от конца до конца по отношению к общей длине контура полимера, где поведение похоже на пружину. Поведение в более широких диапазонах сил можно смоделировать, используя канонический ансамбль, идентичный намагничиванию парамагнитных спинов. Для произвольных сил зависимость силы растяжения будет выражаться функцией Ланжевена L {\ displaystyle {\ mathcal {L}}}{\ displaystyle {\ mathcal {L}}} :

RN l = coth ⁡ (flk BT) - k BT fl = L (flk BT), {\ displaystyle {\ frac {R} {Nl}} = \ coth \ left ({\ frac {fl} {k_ {B} T}} \ right) - {\ frac { k_ {B} T} {fl}} = {\ mathcal {L}} \ left ({\ frac {fl} {k_ {B} T}} \ right),}{\ displaystyle {\ frac {R} {Nl}} = \ coth \ left ({\ frac {fl} {k_ {B} T} } \ right) - {\ frac {k_ {B} T} {fl}} = {\ mathcal {L}} \ left ({\ frac {fl} {k_ {B} T}} \ right),}

где расширение равно R = | R → | {\ displaystyle R = | {\ vec {R}} |}{\ displaystyle R = | {\ vec {R}} |} .

Для произвольных расширений зависимость сила-расширение может быть аппроксимирована следующим образом:

flk BT = L - 1 (RN l) ≈ 3 RN l + 1 5 (RN l) 2 грех ⁡ (7 R 2 N l) + (RN l) 3 1 - RN l {\ displaystyle {\ frac {fl} {k_ {B} T}} = {\ mathcal {L} } ^ {- 1} \ left ({\ frac {R} {Nl}} \ right) \ приблизительно 3 {\ frac {R} {Nl}} + {\ frac {1} {5}} \ left ({ \ frac {R} {Nl}} \ right) ^ {2} \ sin \ left ({\ frac {7R} {2Nl}} \ right) + {\ frac {\ left ({\ frac {R} {Nl }} \ right) ^ {3}} {1 - {\ frac {R} {Nl}}}}}{\ displaystyle {\ frac {fl} {k_ {B} T} } = {\ mathcal {L}} ^ {- 1} \ left ({\ frac {R} {Nl}} \ right) \ приблизительно 3 {\ frac {R} {Nl}} + {\ frac {1} {5}} \ left ({\ frac {R} {Nl}} \ right) ^ {2} \ sin \ left ({\ frac {7R} {2Nl}} \ right) + {\ frac {\ left ( {\ frac {R} {Nl}} \ right) ^ {3}} {1 - {\ frac {R} {Nl}}}}} ,

где L - 1 {\ displaystyle {\ mathcal {L}} ^ {- 1 }}{\ displaystyle {\ mathcal {L}} ^ { -1}} - это обратная функция Ланжевена, N - количество связей в молекуле (поэтому, если в молекуле есть N связей, она имеет N + 1 мономеров, составляющих молекулу.).

Наконец, модель может быть расширена до еще более широких диапазонов сил путем включения модуля упругости по длине контура полимера. То есть, позволяя длине каждого звена цепи упруго реагировать на приложенную силу.

Исправлен фиксированный резервуар идеальной цепи 2.JPG

Идеальная длина замены полимера с резервуаром

В этом подразделе, как и в предыдущем, два конца полимера прикреплены к устройству для микроманипуляции. Однако на этот раз устройство не поддерживает два конца идеальной цепи в фиксированном положении, а поддерживает постоянную тянущую силу f → op {\ displaystyle {\ vec {f}} _ {op} }{\ vec f} _ {{op}} на идеальной цепочке. В этом случае два конца полимера колеблются около среднего положения ⟨R →⟩ {\ displaystyle \ langle {\ vec {R}} \ rangle}\ langle {\ vec R} \ rangle . Идеальная цепь реагирует с постоянной противоположной силой f → = - f → op {\ displaystyle {\ vec {f}} = - {\ vec {f}} _ {op}}{\ vec f} = - {\ vec f} _ {{ op}}

Для идеальной цепи меняя длину с резервуаром, макросостояние системы характеризуется вектором f → {\ displaystyle {\ vec {f}}}\ vec f .

Переход между идеальной цепью фиксированной длины и идеальной цепью в контакте с резервуаром длины очень похоже на изменение между микроканоническим ансамблем и каноническим ансамблем (см. статью Статистическая механика об этом). Переход от состояния, в котором фиксированное значение накладывается на определенный параметр, к состоянию, в котором система может обмениваться этим параметром с внешним миром. Рассматриваемый параметр - это энергия для микроканонического и канонического описаний, тогда как в случае идеальной цепи параметром является длина идеальной цепи.

Как и в микроканоническом и каноническом ансамблях, два описания идеальной цепи различаются только тем, как они трактуют колебания системы. Таким образом, они эквивалентны при термодинамическом пределе. Уравнение состояния идеальной цепи остается прежним, за исключением того, что R → {\ displaystyle {\ vec {R}}}{\ vec {R}} теперь подвержено колебаниям:

е → = - к BT 3 ⟨R →⟩ N l 2 {\ displaystyle {\ vec {f}} = - k_ {B} T {\ frac {3 \ langle {\ vec {R}} \ rangle} { Nl ^ {2}}} ~}{\ vec f} = - k_ {B} T {\ frac {3 \ langle {\ vec R} \ rangle} {Nl ^ { 2}}} ~ .

Идеальная цепь при ограничении постоянной силы - расчет

Схема идеальной цепи, ограниченной постоянной силой.

Рассмотрим свободно соединенную цепь из N связей длины l {\ displaystyle l}l , подверженный постоянной силе удлинения f, приложенной к его концам вдоль оси z, и температуре окружающей среды T {\ displaystyle T}T . Примером может служить цепь с двумя противоположными зарядами + q и -q на концах в постоянном электрическом поле E → {\ displaystyle {\ vec {E}}}{\ vec {E}} применяется вдоль оси z {\ displaystyle z}z , как показано на рисунке справа. Если прямое кулоновское взаимодействие между зарядами игнорируется, то на двух концах существует постоянная сила f → {\ displaystyle {\ vec {f}}}{\ vec {f}} .

Различные конформации цепи не равновероятны, потому что они соответствуют разной энергии цепи во внешнем электрическом поле.

U = - q E → ⋅ R → = - f → ⋅ R → = - f R z {\ displaystyle U = -q {\ vec {E}} \ cdot {\ vec {R}} = - { \ vec {f}} \ cdot {\ vec {R}} = - fR_ {z}}{\ displaystyle U = -q {\ vec {E}} \ cdot {\ vec {R}} = - {\ vec {f}} \ cdot {\ vec {R}} = - fR_ {z}}

Таким образом, разные конформации цепи имеют разные статистические факторы Больцмана exp (- U / k BT) {\ displaystyle exp (-U / k_ {B} T)}{\ displaystyle exp (-U / k_ {B} T)} .

функция распределения :

Z = ∑ statesexp (- U / k BT) = ∑ statesexp (f R zk BT) {\ Displaystyle Z = \ сумма _ {состояния} ехр (-U / k_ {B} T) = \ sum _ {состояния} ехр ({fR_ {z} \ над k_ {B} T})}{\ displaystyle Z = \ sum _ {состояния} exp (-U / k_ {B} T) = \ sum _ {состояния} exp ({fR_ {z} \ over k_ {B} T})}

Каждое соединение мономера в цепочке характеризуется вектором ri → {\ displaystyle {\ vec {r_ {i}}}}{\ displaystyle {\ vec {r_ {i}}}} длины l {\ displaystyle l}l и углы θ i, φ i {\ displaystyle \ theta _ {i}, \ varphi _ {i}}{\ displaystyle \ theta _ {i}, \ varphi _ {i}} в сферическая система координат. Сквозной вектор может быть представлен как: R z = ∑ i = 1 N l cos ⁡ θ i {\ displaystyle R_ {z} = \ sum _ {i = 1} ^ {N} l { \ text {}} \ cos \ theta _ {i}}{\ displaystyle R_ {z } = \ sum _ {i = 1} ^ {N} l {\ text {}} \ cos \ theta _ {i}} . Следовательно:

Z = ∫ exp ⁡ (flk BT ∑ i = 1 N cos ⁡ θ i) ∏ i = 1 N sin ⁡ θ id θ id φ i = [∫ 0 π 2 π sin ⁡ θ i exp ⁡ ( flk BT cos ⁡ θ i) d θ i] N = [2 π fl / (k BT) [exp ⁡ (flk BT) - exp ⁡ (- flk BT)]] N = [4 π sinh ⁡ (fl / ( KBT)) fl / (к BT)] N {\ displaystyle {\ begin {align} Z = \ int \ exp ({fl \ over k_ {B} T} \ sum _ {i = 1} ^ {N} \ cos \ theta _ {i}) \ prod _ {i = 1} ^ {N} \ sin \ theta _ {i} d \ theta _ {i} d \ varphi _ {i} = \\ \ left [ \ int _ {0} ^ {\ pi} 2 \ pi {\ text {}} \ sin \ theta _ {i} {\ text {}} \ exp ({fl \ over k_ {B} T} \ cos \ theta _ {i}) d \ theta _ {i} \ right] ^ {N} = \\ \ left [{2 \ pi \ over fl / (k_ {B} T)} [\ exp ({fl \ над k_ {B} T}) - \ exp (- {fl \ over k_ {B} T})] \ right] ^ {N} = \\ \ left [{4 \ pi {\ text {}} \ sinh (fl / (K_ {B} T)) \ over fl / (k_ {B} T)} \ right] ^ {N} \ end {align}}}{\ displaystyle {\ begin {align} Z = \ int \ exp ({fl \ over k_ { B} T} \ sum _ {i = 1} ^ {N} \ cos \ theta _ {i}) \ prod _ {i = 1} ^ {N} \ sin \ theta _ {i} d \ theta _ {i} d \ varphi _ {i} = \\ \ left [\ int _ {0} ^ {\ pi } 2 \ pi {\ text {}} \ sin \ theta _ {i} {\ text {}} \ exp ({fl \ over k_ {B} T} \ cos \ theta _ {i}) d \ theta _ {i} \ right] ^ {N} = \\ \ left [{2 \ pi \ over fl / (k_ {B} T)} [\ exp ({fl \ over k_ {B} T}) - \ exp (- {fl \ over k_ {B} T})] \ right] ^ {N} = \\ \ left [{4 \ pi {\ text {}} \ sinh (fl / (K_ {B} T)) \ over fl / (k_ {B} T)} \ right] ^ {N} \ end {align}}}

Свободная энергия Гиббса G можно напрямую вычислить из статистической суммы:

G (T, f, N) = - k BT ln Z (T, f, N) = - N k BT [ln (4 π sinh (flk BT)) - пер (flk BT)] {\ displaystyle G (T, f, N) = - k_ {B} T { \ text {}} ln {\ text {}} Z (T, f, N) = - Nk_ {B} T [ln (4 \ pi {\ text {}} sinh ({fl \ over k_ {B} T })) - ln ({fl \ over k_ {B} T})]}{\ displaystyle G (T, f, N) = - k_ {B} T {\ text {}} ln {\ text {}} Z (T, f, N) = - Nk_ {B} T [ln (4 \ pi {\ text {}} sinh ({fl \ over k_ {B } T})) - ln ({fl \ over k_ {B} T})]}

Здесь используется свободная энергия Гиббса, потому что ансамбль цепочек соответствует постоянной температуре T {\ displaystyle T}T и постоянная сила f {\ displaystyle f}f (аналог изотермино-изобарического ансамбля, который имеет постоянные температура и давление).

Среднее расстояние от конца до конца, соответствующее данной силе, можно получить как производную от свободной энергии:

< R>= - ∂ G ∂ f = N l [coth (flk BT) - 1 fl / (k BT)] {\ displaystyle = - {\ partial G \ over \ partial f} = Nl [coth ({fl \ over k_ {B} T}) - {1 \ over fl / (k_ {B} T)}]}{\displaystyle <R>= - {\ partial G \ over \ partial f} = Nl [coth ({fl \ over k_ {B} T}) - {1 \ over fl / ( k_ {B} T)}]}

Это выражение является функцией Ланжевена L {\ displaystyle {\ mathcal {L}}}{\ mathcal {L}} , также упомянутой в предыдущих абзацах:

среднее расстояние ⟨R →⟩ {\ displaystyle \ left \ langle {\ vec {R}} \ right \ rangle}{\ displaystyle \ left \ langle {\ vec {R}} \ right \ rangle} цепи как функция от α {\ displaystyle \ alpha }\ alpha .

L (α) = coth (α) - 1 α {\ displaystyle {\ mathcal {L}} (\ alpha) = coth (\ alpha) - {1 \ over \ alpha}}{\ displaystyle {\ mathcal {L}} (\ alpha) = coth (\ alpha) - {1 \ over \ alpha}} где, α = flk BT {\ display style \ alpha = {fl \ over k_ {B} T}}{\ displaystyle \ alpha = {fl \ over k_ {B} T}} .

Для малых относительных удлинений (⟨R⟩ ≪ R max = l N {\ displaystyle \ left \ langle R \ right \ rangle \ ll R_ {max} = lN}{\ displaystyle \ left \ langle R \ right \ rangle \ ll R_ {max} = lN} ) зависимость приблизительно линейна,

L (α) ≅ α 3 {\ displaystyle {\ mathcal {L}} (\ alpha) \ cong {\ alpha \ over 3}}{\ displaystyle {\ mathcal {L }} (\ альфа) \ cong {\ alpha \ over 3}} для α ≪ 1 {\ displaystyle \ alpha \ ll 1}{\ displaystyle \ alpha \ ll 1}

и следует закону Гука, как показано в предыдущих абзацах:

f → = KBT 3 ⟨R →⟩ N l 2 {\ displaystyle {\ vec {f}} = K_ {B} T {3 \ left \ langle {\ vec {R}} \ right \ rangle \ over Nl ^ {2}} }{\ displaystyle {\ vec {f}} = K_ {B} T {3 \ left \ langle {\ vec {R}} \ right \ rangle \ over Nl ^ {2}}} .

См. Также
Внешние ссылки
Ссылки
Последняя правка сделана 2021-05-23 10:28:31
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте