Теория поля полимера

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

A теория поля полимера - это статистическая теория поля, описывающая статистическое поведение нейтрального или заряженного полимерная система. Его можно получить, преобразовав статистическую сумму из его стандартного многомерного интегрального представления по степеням свободы частиц в представление функционального интеграла по вспомогательному полю функция, используя либо преобразование Хаббарда – Стратоновича, либо дельта-функциональное преобразование. Компьютерное моделирование, основанное на теориях полимерного поля, показало полезные результаты, например, для расчета структур и свойств растворов полимеров (Baeurle 2007, Schmid 1998), расплавов полимеров (Schmid 1998, Matsen 2002, Fredrickson 2002) и термопласты (Baeurle 2006).

Содержание

  • 1 Канонический ансамбль
    • 1.1 Частичное представление канонической статистической суммы
    • 1.2 Теоретико-полевое преобразование
    • 1.3 Базовое теоретико-полевое представление канонической статистической суммы
  • 2 Большой канонический ансамбль
    • 2.1 Базовое теоретико-полевое представление большой канонической статистической суммы
  • 3 Аппроксимация среднего поля
  • 4 Поправки высшего порядка
  • 5 Методы перенормировки
    • 5.1 Теория ренормгруппы
    • 5.2 Перенормировка Хартри
    • 5.3 Перенормировка головастика
  • 6 Численное моделирование
    • 6.1 Представление среднего поля
    • 6.2 Гауссовское эквивалентное представление
  • 7 Ссылки
  • 8 Внешние ссылки

Канонический ансамбль

Частичное представление канонического Статистическая сумма

Стандартная континуальная модель гибких полимеров, введенная Эдвардсом (Edwards, 1965), рассматривает раствор, состоящий из n {\ displaystyle n}n линейных монодисперсных гомополимеров как систему из крупнозернистого полимера s, в котором статистическая механика цепей описывается моделью непрерывной гауссовой нити (Baeurle 2007), а растворитель учитывается неявно. Модель гауссовой нити можно рассматривать как непрерывный предел модели дискретной гауссовой цепи, в которой полимеры описываются как непрерывные линейно упругие нити. Каноническая статистическая сумма такой системы, сохраненная при обратной температуре β = 1 / k BT {\ displaystyle \ beta = 1 / k_ {B} T}{\ displaystyle \ beta = 1 / k_ {B} T} и ограниченная объемом V {\ displaystyle V}V, может быть выражено как

Z (n, V, β) = 1 n! (λ T 3) N N ∏ J знак равно 1 N ∫ D rj ехр ⁡ (- β Φ 0 [r] - β Φ ¯ [r]), (1) {\ displaystyle Z (n, V, \ beta) = {\ frac {1} {n! (\ lambda _ {T} ^ {3}) ^ {nN}}} \ prod _ {j = 1} ^ {n} \ int D \ mathbf {r} _ {j } \ exp \ left (- \ beta \ Phi _ {0} \ left [\ mathbf {r} \ right] - \ beta {\ bar {\ Phi}} \ left [\ mathbf {r} \ right] \ right), \ qquad (1)}{\ displaystyle Z (n, V, \ beta) = {\ frac {1} {n! (\ Lambda _ {T} ^ {3}) ^ {nN}}} \ prod _ {j = 1} ^ {n} \ int D \ mathbf {r} _ {j} \ exp \ left (- \ beta \ Phi _ {0} \ left [\ mathbf {r} \ right] - \ beta {\ bar {\ Phi}} \ left [\ ma thbf {r} \ right] \ right), \ qquad (1)}

где Φ ¯ [r] {\ displaystyle {\ bar {\ Phi}} \ left [\ mathbf {r} \ right]}{\ displaystyle {\ bar {\ Phi}} \ left [\ mathbf {r} \ right]} - потенциал средней силы, определяемый как

Φ ¯ [r] = N 2 2 ∑ j = 1 n ∑ k = 1 n ∫ 0 1 ds ∫ 0 1 ds ′ Φ ¯ (| rj (s) - rk (s ') |) - 1 2 N N Φ ¯ (0), (2) {\ displaystyle {\ bar {\ Phi}} \ left [\ mathbf {r} \ right] = {\ гидроразрыв {N ^ {2}} {2}} \ sum _ {j = 1} ^ {n} \ sum _ {k = 1} ^ {n} \ int _ {0} ^ {1} ds \ int _ {0} ^ {1} ds '{\ bar {\ Phi}} \ left (\ left | \ mathbf {r} _ {j} (s) - \ mathbf {r} _ {k} (s') \ right | \ right) - {\ frac {1} {2}} nN {\ bar {\ Phi}} (0), \ qquad (2)}{\displaystyle {\bar {\Phi }}\left[\mathbf {r} \right]={\frac {N^{2}}{2}}\sum _{j=1}^{n}\sum _{k=1}^{n}\int _{0}^{1}ds\int _{0}^{1}ds'{\bar {\Phi }}\left(\left|\mathbf {r} _{j}(s)-\mathbf {r} _{k}(s')\right|\right)-{\frac {1}{2}}nN{\bar {\Phi }}(0),\qquad (2)}

представляет опосредованные растворителем несвязанные взаимодействия между сегментами, а Φ 0 [r] {\ displaystyle \ Phi _ {0} [\ mathbf {r}]}{\ displaystyle \ Phi _ {0} [\ mathbf {r}]} представляет h армоническая энергия связи цепей. Последний вклад в энергию можно сформулировать как

Φ 0 [r] = 3 k B T 2 N b 2 ∑ l = 1 n ∫ 0 1 d s | d r l (s) d s | 2, {\ displaystyle \ Phi _ {0} [\ mathbf {r}] = {\ frac {3k_ {B} T} {2Nb ^ {2}}} \ sum _ {l = 1} ^ {n} \ int _ {0} ^ {1} ds \ left | {\ frac {d \ mathbf {r} _ {l} (s)} {ds}} \ right | ^ {2},}{\ displaystyle \ Phi _ {0} [\ mathbf {r}] = {\ frac {3k_ {B} T} {2Nb ^ {2}}} \ sum _ {l = 1} ^ {n} \ int _ {0} ^ {1} ds \ left | {\ frac {d \ mathbf {r} _ {l} (s)} {ds}} \ right | ^ {2},}

где b {\ displaystyle b}b - длина статистического сегмента, а N {\ displaystyle N}N - индекс полимеризации.

Теоретико-полевое преобразование

Чтобы вывести базовое теоретико-полевое представление канонической статистической суммы, ниже вводится оператор сегментной плотности полимерной системы

ρ ^ (r) = N ∑ j = 1 n ∫ 0 1 ds δ (r - rj (s)). {\ Displaystyle {\ шляпа {\ rho}} (\ mathbf {r}) = N \ sum _ {j = 1} ^ {n} \ int _ {0} ^ {1} ds \ delta \ left (\ mathbf {r} - \ mathbf {r} _ {j} (s) \ right).}{\ displaystyle {\ hat {\ rho}} (\ mathbf {r}) = N \ sum _ {j = 1} ^ {n} \ int _ {0} ^ {1} ds \ delta \ left (\ mathbf {r} - \ mathbf {r} _ {j} (s) \ right).}

Используя это определение, можно переписать уравнение. (2) как

Φ ¯ [r] = 1 2 ∫ dr ∫ dr ′ ρ ^ (r) Φ ¯ (| r - r ′ |) ρ ^ (r ′) - 1 2 n N Φ ¯ (0). (3) {\ displaystyle {\ bar {\ Phi}} \ left [\ mathbf {r} \ right] = {\ frac {1} {2}} \ int d \ mathbf {r} \ int d \ mathbf { r} '{\ hat {\ rho}} (\ mathbf {r}) {\ bar {\ Phi}} (\ left | \ mathbf {r} - \ mathbf {r}' \ right |) {\ hat { \ rho}} (\ mathbf {r} ') - {\ frac {1} {2}} nN {\ bar {\ Phi}} (0). \ qquad (3)}{\displaystyle {\bar {\Phi }}\left[\mathbf {r} \right]={\frac {1}{2}}\int d\mathbf {r} \int d\mathbf {r} '{\hat {\rho }}(\mathbf {r}){\bar {\Phi }}(\left|\mathbf {r} -\mathbf {r} '\right|){\hat {\rho }}(\mathbf {r} ')-{\frac {1}{2}}nN{\bar {\Phi }}(0).\qquad (3)}

Затем выполняется преобразование в теорию поля с помощью преобразования Хаббарда-Стратоновича или дельта-функционального преобразования

∫ D ρ δ [ρ - ρ ^] F [ρ] = F [ρ ^], ( 4) {\ displaystyle \ int D \ rho \; \ delta \ left [\ rho - {\ hat {\ rho}} \ right] F \ left [\ rho \ right] = F \ left [{\ hat {\ rho}} \ right], \ qquad (4)}{\ displaystyle \ int D \ rho \; \ delta \ left [\ rho - {\ hat {\ rho}} \ right] F \ left [\ rho \ right] = F \ left [{\ hat {\ rho}} \ right], \ qquad (4)}

где F [ρ ^] {\ displaystyle F \ left [{\ hat {\ rho}} \ right]}{\ displaystyle F \ left [{\ hat {\ rho}} \ right]} - это функционал, а δ [ρ - ρ ^] {\ displaystyle \ delta \ left [\ rho - {\ hat {\ rho}} \ right]}{\ displaystyle \ delta \ left [\ rho - {\ hat {\ rho}} \ right]} - дельта-функционал, задаваемый

δ [ρ - ρ ^] знак равно ∫ D wei ∫ drw (r) [ρ (r) - ρ ^ (r)], (5) {\ displaystyle \ delta \ left [\ rho - {\ hat {\ rho}} \ right] = \ int Dwe ^ {i \ int d \ mathbf {r} w (\ mathbf {r}) \ left [\ rho (\ mathbf {r}) - {\ hat {\ rho}} (\ mathbf {r}) \ right]}, \ qquad (5)}{\ displaystyle \ delta \ left [\ rho - {\ hat {\ rho}} \ right] = \ int Dwe ^ {i \ int d \ mathbf {r} w (\ mathbf {r}) \ left [ \ rho (\ mathbf {r}) - {\ hat {\ rho}} (\ mathbf {r}) \ right]}, \ qquad (5)}

с w (r) = ∑ г вес (G) ехр ⁡ [я г р] {\ displaystyle w (\ mathbf {r}) = \ sum \ nolimits _ {\ mathbf {G}} w (\ mathbf {G}) \ exp \ left [ i \ mathbf {G} \ mathbf {r} \ right]}{\ displaystyle w (\ mathbf {r}) = \ sum \ nolimits _ {\ mathbf {G}} вес (\ mathbf {G}) \ exp \ left [я \ mathbf {G} \ mathbf {r} \ right]} , представляющий функцию вспомогательного поля. Здесь мы отмечаем, что расширение функции поля в ряд Фурье означает, что периодические граничные условия применяются во всех направлениях и что G {\ displaystyle \ mathbf {G}}\ mathbf {G} -векторы обозначают векторы обратной решетки сверхъячейки.

Базовое теоретико-полевое представление канонической статистической суммы

Используя уравнения. Используя уравнения (3), (4) и (5), мы можем преобразовать каноническую статистическую сумму в уравнение (1) в теоретико-полевом представлении, которое приводит к

Z (n, V, β) = Z 0 ∫ D w exp ⁡ [- 1 2 β V 2 ∫ drdr ′ w (r) Φ ¯ - 1 ( р - r ') вес (r')] Q N [iw], (6) {\ displaystyle Z (n, V, \ beta) = Z_ {0} \ int Dw \ exp \ left [- {\ frac { 1} {2 \ beta V ^ {2}}} \ int d \ mathbf {r} d \ mathbf {r} 'w (\ mathbf {r}) {\ bar {\ Phi}} ^ {- 1} ( \ mathbf {r} - \ mathbf {r} ') w (\ mathbf {r}') \ right] Q ^ {n} [iw], \ qquad (6)}{\displaystyle Z(n,V,\beta)=Z_{0}\int Dw\exp \left[-{\frac {1}{2\beta V^{2}}}\int d\mathbf {r} d\mathbf {r} 'w(\mathbf {r}){\bar {\Phi }}^{-1}(\mathbf {r} -\mathbf {r} ')w(\mathbf {r} ')\right]Q^{n}[iw],\qquad (6)}

где

Z 0 = 1 п! (ехр ⁡ (β / 2 N Φ ¯ (0)) Z ′ λ 3 N (T)) n {\ displaystyle Z_ {0} = {\ frac {1} {n!}} \ left ({\ frac { \ exp \ left (\ beta / 2N {\ bar {\ Phi}} (0) \ right) Z '} {\ lambda ^ {3N} (T)}} \ right) ^ {n}}{\displaystyle Z_{0}={\frac {1}{n!}}\left({\frac {\exp \left(\beta /2N{\bar {\Phi }}(0)\right)Z'}{\lambda ^{3N}(T)}}\right)^{n}}

может интерпретироваться как статистическая сумма для идеального газа невзаимодействующих полимеров и

Z ′ = ∫ DR exp ⁡ [- β U 0 (R)] (7) {\ displaystyle Z '= \ int D \ mathbf { R} \ exp \ left [- \ beta U_ {0} (\ mathbf {R}) \ right] \ qquad (7)}{\displaystyle Z'=\int D\mathbf {R} \exp \left[-\beta U_{0}(\mathbf {R})\right]\qquad (7)}

- интеграл по путям свободного полимера в нулевом поле с упругой энергией

U 0 [R] = k BT 4 R g 0 2 ∫ 0 1 ds | d R (s) d s | 2. {\ displaystyle U_ {0} [\ mathbf {R}] = {\ frac {k_ {B} T} {4R_ {g0} ^ {2}}} \ int _ {0} ^ {1} ds \ left | {\ frac {d \ mathbf {R} (s)} {ds}} \ right | ^ {2}.}{\ displaystyle U_ { 0} [\ mathbf {R}] = {\ frac {k_ {B} T} {4R_ {g0} ^ {2}}} \ int _ {0} ^ {1} ds \ left | {\ frac {d \ mathbf {R} (s)} {ds}} \ right | ^ {2}.}

В последнем уравнении невозмущенный радиус вращения цепи R g 0 = N б 2 / (6) {\ displaystyle R_ {g0} = {\ sqrt {Nb ^ {2} / (6)}}}{\ displaystyle R_ {g0} = {\ sqrt {Nb ^ {2} / (6)}}} . Более того, в формуле. (6) статистическая сумма одного полимера, подчиненного полю w (R) {\ displaystyle w (\ mathbf {R})}{\ displaystyle w (\ mathbf {R})} , задается как

Q [iw ] = DR exp ⁡ [- β U 0 [R] - i N ∫ 0 1 dsw (R (s))] ∫ DR exp exp [- β U 0 [R]]. (8) {\ displaystyle Q [iw] = {\ frac {\ int D \ mathbf {R} \ exp \ left [- \ beta U_ {0} [\ mathbf {R}] -iN \ int _ {0} ^ {1} ds \; w (\ mathbf {R} (s)) \ right]} {\ int D \ mathbf {R} \ exp \ left [- \ beta U_ {0} [\ mathbf {R}] \ right]}}. \ qquad (8)}{\ displaystyle Q [iw] = {\ frac {\ int D \ mathbf {R} \ exp \ left [- \ beta U_ {0} [\ mathbf {R}] -iN \ int _ {0} ^ {1} ds \; w (\ mathbf {R} (s)) \ right]} {\ int D \ mathbf {R} \ exp \ left [- \ beta U_ {0} [\ mathbf {R}] \ right]}}. \ Qquad (8)}

Большой канонический ансамбль

Базовое теоретико-полевое представление большой канонической статистической суммы

Для вывода большой канонической статистической суммы мы используем его стандартное термодинамическое отношение к канонической статистической сумме, заданное формулой

Ξ (μ, V, β) = ∑ n = 0 ∞ e β μ n Z (n, V, β), {\ displaystyle \ Xi (\ mu, V, \ beta) = \ sum _ {n = 0} ^ {\ infty} e ^ {\ beta \ mu n} Z (n, V, \ beta),}{\ displaystyle \ Xi (\ mu, V, \ beta) = \ sum _ {n = 0} ^ {\ infty} e ^ {\ beta \ mu n} Z (n, V, \ beta),}

где μ {\ displaystyle \ mu}\ mu - химический потенциал, а Z (n, V, β) {\ displaystyle Z (n, V, \ beta)}{\ displaystyle Z (n, V, \ beta)} определяется формулой. (6). Суммирование дает теоретико-полевое представление большой канонической статистической суммы,

Ξ (ξ, V, β) = γ Φ ¯ ∫ D w exp ⁡ [- S [w]], {\ displaystyle \ Xi (\ xi, V, \ beta) = \ gamma _ {\ bar {\ Phi}} \ int Dw \ exp \ left [-S [w] \ right],}{\ displaystyle \ Xi (\ xi, V, \ beta) = \ gamma _ {\ bar {\ Phi}} \ int Dw \ exp \ left [-S [w] \ right],}

где

S [w ] Знак равно 1 2 β В 2 ∫ drdr ′ вес (г) Φ ¯ - 1 (г - г ') вес (г') - ξ Q [iw] {\ Displaystyle S [ш] = {\ гидроразрыва {1} { 2 \ beta V ^ {2}}} \ int d \ mathbf {r} d \ mathbf {r} 'w (\ mathbf {r}) {\ bar {\ Phi}} ^ {- 1} (\ mathbf { r} - \ mathbf {r} ') w (\ mathbf {r}') - \ xi Q [iw]}{\displaystyle S[w]={\frac {1}{2\beta V^{2}}}\int d\mathbf {r} d\mathbf {r} 'w(\mathbf {r}){\bar {\Phi }}^{-1}(\mathbf {r} -\mathbf {r} ')w(\mathbf {r} ')-\xi Q[iw]}

- это большое каноническое действие с Q [iw] {\ displaystyle Q [iw] }{\ displaystyle Q [iw]} определяется формулой. (8) и константа

γ Φ ¯ = 1 2 ∏ G (1 π β Φ ¯ (G)) 1/2. {\ displaystyle \ gamma _ {\ bar {\ Phi}} = {\ frac {1} {\ sqrt {2}}} \ prod _ {\ mathbf {G}} \ left ({\ frac {1} {\ pi \ beta {\ bar {\ Phi}} (\ mathbf {G})}} \ right) ^ {1/2}.}{\ displaystyle \ gamma _ {\ bar {\ Phi}} = {\ frac {1} {\ sqrt {2}}} \ prod _ {\ mathbf {G } } \ left ({\ frac {1} {\ pi \ beta {\ bar {\ Phi}} (\ mathbf {G})}} \ right) ^ {1/2}.}

Кроме того, параметр, связанный с химическим потенциалом, задается как

ξ знак равно ехр ⁡ (β μ + β / 2 N Φ ¯ (0)) Z ′ λ 3 N (T), {\ displaystyle \ xi = {\ frac {\ exp \ left (\ beta \ mu + \ beta / 2N {\ bar {\ Phi}} (0) \ right) Z '} {\ lambda ^ {3N} (T)}},}{\displaystyle \xi ={\frac {\exp \left(\beta \mu +\beta /2N{\bar {\Phi }}(0)\right)Z'}{\lambda ^{3N}(T)}},}

где Z ′ {\ displaystyle Z'}Z'обеспечивается формулой. (7).

Приближение среднего поля

Стандартная стратегия приближения для теорий поля полимеров - это приближение среднего поля (MF), которое заключается в замене члена взаимодействия многих тел в действие по члену, в котором все тела системы взаимодействуют со средним эффективным полем. Этот подход сводит любую проблему с несколькими телами к эффективной задаче с одним телом, предполагая, что интеграл статистической суммы модели определяется конфигурацией одного поля. Основное преимущество решения проблем с помощью приближения МП или его численной реализации, обычно называемой теорией самосогласованного поля (ССПП), заключается в том, что оно часто дает полезные сведения о свойствах и поведении сложных систем многих тел при относительно низкая стоимость вычислений. Успешные применения этой стратегии приближения могут быть найдены для различных систем полимеров и сложных жидкостей, таких как, например, сильно сегрегированные блок-сополимеры с высокой молекулярной массой, высококонцентрированные растворы нейтральных полимеров или высококонцентрированные растворы блок-полиэлектролита (PE) (Schmid 1998, Matsen 2002, Fredrickson 2002). Однако существует множество случаев, когда SCFT дает неточные или даже качественно неверные результаты (Baeurle 2006a). К ним относятся растворы нейтрального полимера или полиэлектролита в разбавленных и полуразбавленных режимах концентрации, блок-сополимеры вблизи их перехода порядок-беспорядок, смеси полимеров вблизи фазовых переходов и т. Д. В таких ситуациях интеграл статистической суммы, определяющий теоретико-полевую модель, не полностью определяется Конфигурация одиночного МП и далекие от него конфигурации поля могут внести важный вклад, который требует использования более сложных методов расчета, выходящих за пределы уровня приближения МП.

Поправки высшего порядка

Одна из возможностей решения проблемы - вычислить поправки высшего порядка к приближению МП. Цончев и др. разработал такую ​​стратегию, включающую поправки на флуктуации ведущего (однопетлевого) порядка, что позволило по-новому взглянуть на физику замкнутых растворов ПЭ (Цончев, 1999). Однако в ситуациях, когда приближение МП является плохим, для получения желаемой точности необходимо много требующих вычислений поправок к интегралу более высокого порядка.

Методы перенормировки

Альтернативный теоретический инструмент для решения проблем сильных флуктуаций, возникающих в теориях поля, был предоставлен в конце 1940-х годов концепцией перенормировки, которая изначально использовалась был разработан для вычисления функциональных интегралов, возникающих в квантовых теориях поля (КТП). В QFT стандартная стратегия аппроксимации заключается в расширении функциональных интегралов в степенной ряд по константе связи с использованием теории возмущений. К сожалению, как правило, большинство членов расширения оказываются бесконечными, что делает такие вычисления невозможными (Ширков 2001). Один из способов избавиться от бесконечностей из QFT - это использовать концепцию перенормировки (Baeurle 2007). В основном он заключается в замене голых значений параметров связи, например, электрические заряды или массы, перенормированными параметрами связи и требуя, чтобы физические величины не изменялись при этом преобразовании, тем самым приводя к конечным членам в разложении возмущений. Простую физическую картину процедуры перенормировки можно получить на примере классического электрического заряда, Q {\ displaystyle Q}Q , помещенного в поляризуемую среду, например, в раствор электролита. На расстоянии r {\ displaystyle r}r от заряда из-за поляризации среды его кулоновское поле будет эффективно зависеть от функции Q (r) {\ displaystyle Q (r)}{\ displaystyle Q (r) } , т.е. эффективный (перенормированный) заряд, вместо чистого электрического заряда, Q {\ displaystyle Q}Q . В начале 1970-х годов К.Г. Вильсон стал пионером в силе концепций перенормировки, разработав формализм теории ренормгруппы (RG) для исследования критических явлений статистических систем (Wilson, 1971).

Теория ренормгруппы

Теория RG использует серию преобразований RG, каждое из которых состоит из шага грубой зернистости, за которым следует изменение масштаба (Wilson 1974). В случае статистико-механических задач эти шаги реализуются путем последовательного исключения и изменения масштаба степеней свободы в сумме разбиений или интеграле, который определяет рассматриваемую модель. Де Женн использовал эту стратегию, чтобы установить аналогию между поведением классической векторной модели с нулевыми компонентами ферромагнетизма вблизи фазового перехода и самопроизвольного случайного блуждания полимерной цепи бесконечной длины на решетке, для расчета показателей исключенного объема полимера (de Gennes, 1972). Адаптация этой концепции к теоретико-полевым функциональным интегралам подразумевает систематическое изучение того, как модель теории поля изменяется при удалении и изменении масштаба определенного числа степеней свободы из интеграла статистической суммы (Wilson 1974).

перенормировка Хартри

Альтернативный подход известен как приближение Хартри или самосогласованное однопетлевое приближение (Amit 1984). Он использует поправки на гауссовские флуктуации для вклада МФ 0-го {\ displaystyle 0 ^ {th}}0 ^ {th} -порядка, чтобы перенормировать параметры модели и самосогласованным образом извлечь доминирующую длину масштаб колебаний концентрации в режимах критической концентрации.

Перенормировка головастика

В более поздней работе Ефимов и Ноговицин показали, что альтернативный метод перенормировки, происходящий из QFT, основанный на концепции перенормировки головастика, может быть очень эффективным подходом для вычисления функциональных интегралов. возникающие в статистической механике классических систем многих частиц (Ефимов 1996). Они продемонстрировали, что основной вклад в классические интегралы статистической суммы вносят диаграммы Фейнмана типа головастиков низкого порядка, которые учитывают расходящиеся вклады из-за самодействия частицы . Процедура перенормировки, выполняемая в этом подходе, влияет на вклад самодействия заряда (например, электрона или иона) в результате статической поляризации, индуцированной в вакууме из-за наличия этого заряда (Baeurle 2007). Как свидетельствуют Ефимов и Ганболд в более ранней работе (Ефимов, 1991), процедура перенормировки головастика может быть очень эффективно использована для устранения расхождений в действии основного теоретико-полевого представления статистической суммы и приводит к альтернативному функциональному интегралу представление, называемое гауссовским эквивалентным представлением (GER). Они показали, что процедура обеспечивает функциональные интегралы со значительно улучшенными свойствами сходимости для аналитических расчетов возмущений. В последующих работах Baeurle et al. разработали эффективные недорогие методы аппроксимации, основанные на процедуре перенормировки головастика, которые показали полезные результаты для прототипов полимеров и растворов ПЭ (Baeurle 2006a, Baeurle 2006b, Baeurle 2007a).

Численное моделирование

Другой возможностью является использование алгоритмов Монте-Карло (MC) и выборка полного интеграла статистической суммы в теоретико-полевой формулировке. Полученная процедура затем называется полимером теоретико-полевым моделированием. В недавней работе, однако, Баерле продемонстрировал, что выборка MC в сочетании с базовым теоретико-полевым представлением неосуществима из-за так называемой проблемы числового знака (Baeurle 2002). Сложность связана со сложным и колеблющимся характером получаемой функции распределения, что приводит к плохой статистической сходимости средних значений по ансамблю желаемых термодинамических и структурных величин. В таких случаях необходимы специальные аналитические и численные методы для ускорения статистической сходимости (Baeurle 2003, Baeurle 2003a, Baeurle 2004).

Представление среднего поля

Чтобы сделать методологию пригодной для вычислений, Баёрл предложил сместить контур интегрирования интеграла статистической суммы через однородное решение МП, используя интегральную теорему Коши, обеспечивая так называемое представление среднего поля. Эта стратегия ранее успешно применялась Baer et al. в теоретико-полевых расчетах электронной структуры (Baer 1998). Баерле смог продемонстрировать, что этот метод обеспечивает значительное ускорение статистической сходимости средних значений ансамбля в процедуре выборки MC (Baeurle 2002, Baeurle 2002a).

Гауссово эквивалентное представление

В последующих работах Baeurle et al. (Baeurle 2002, Baeurle 2002a, Baeurle 2003, Baeurle 2003a, Baeurle 2004) применили концепцию перенормировки головастика, приводящую к гауссовскому эквивалентному представлению интеграла статистической суммы, в сочетании с передовыми методами МК в большом каноническом ансамбле. Они могли убедительно продемонстрировать, что эта стратегия обеспечивает дальнейшее повышение статистической сходимости желаемых средних ансамблевых значений (Baeurle 2002).

Ссылки

  • Baeurle, S.A.; Ноговицин, Е.А. (2007). «Сложные законы масштабирования гибких полиэлектролитных растворов с эффективными концепциями перенормировки». Полимер. 48 (16): 4883. doi : 10.1016 / j.polymer.2007.05.080.
  • Fredrickson, G.H.; Ganesan, V.; Дроле, Ф. (2002). "Теоретико-полевые методы компьютерного моделирования полимеров и сложных жидкостей". Макромолекулы. 35 : 16. Bibcode : 2002MaMol..35... 16F. doi : 10.1021 / ma011515t.
  • Baeurle, S.A.; Усами, Т.; Гусев, А.А. (2006). «Новый подход многомасштабного моделирования для прогнозирования механических свойств полимерных наноматериалов». Полимер. 47 (26): 8604. doi : 10.1016 / j.polymer.2006.10.017.
  • de Gennes, P.G. (1972). «Экспоненты для проблемы исключенного объема, полученные методом Вильсона». Phys. Lett. 38 A : 339.
  • Ефимов Г.В.; Ноговицин, Е.А. (1996). «Статистические суммы классических систем в гауссовском эквивалентном представлении функциональных интегралов». Physica A. 234 : 506. Bibcode : 1996PhyA..234..506V. doi : 10.1016 / S0378-4371 (96) 00279-8.
  • Baeurle, S.A.; Ефимов, Г.В.; Ноговицин, Е.А. (2006b). «О новой теории самосогласованного поля для канонического ансамбля». J. Chem. Phys. 124 (22): 224110. Bibcode : 2006JChPh.124v4110B. doi : 10.1063 / 1.2204913. PMID 16784266.
  • Baeurle, S.A.; Шарло, М.; Ноговицин Е.А. (2007a). «Грандиозные канонические исследования прототипных моделей полиэлектролитов за пределами приближения среднего поля». Phys. Ред. E. 75 : 011804. Bibcode : 2007PhRvE..75a1804B. doi : 10.1103 / PhysRevE.75.011804.
  • Baeurle, S.A. (2002). "Метод гауссовского эквивалентного представления: новый метод уменьшения знаковой проблемы функциональных интегральных методов". Phys. Rev. Lett. 89 (8): 080602. Bibcode : 2002PhRvL..89h0602B. doi : 10.1103 / PhysRevLett.89.080602. PMID 12190451.
  • Baeurle, S.A. (2003a). «Метод Монте-Карло вспомогательного поля стационарной фазы: новая стратегия уменьшения проблемы знака в методологиях вспомогательного поля». Comput. Phys. Commun. 154 (2): 111. Bibcode : 2003CoPhC.154..111B. doi : 10.1016 / S0010-4655 (03) 00284-4.
  • Baeurle, S.A. (2004). «Большое каноническое вспомогательное поле Монте-Карло: новый метод моделирования открытых систем с высокой плотностью». Comput. Phys. Commun. 157 (3): 201. Bibcode : 2004CoPhC.157..201B. doi : 10.1016 / j.comphy.2003.11.001.
  • Baer, ​​R.; Хед-Гордон, М.; Нойхаузер, Д. (1998). "Вспомогательное поле со смещенным контуром Монте-Карло для ab initio электронной структуры: преодоление знаковой проблемы". J. Chem. Phys. 109 (15): 6219. Bibcode : 1998JChPh.109.6219B. doi : 10,1063 / 1,477300.
  • Baeurle, S.A.; Martonak, R.; Парринелло, М. (2002a). «Теоретико-полевой подход к моделированию в классическом каноническом и большом каноническом ансамбле». J. Chem. Phys. 117 (7): 3027. Bibcode : 2002JChPh.117.3027B. doi : 10.1063 / 1.1488587.

Внешние ссылки

Последняя правка сделана 2021-06-02 10:33:22
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте