Теория полимера поля является статистической теорией поля , описывающая статистического поведение нейтральной или заряженной полимерной системы. Он может быть получен путем преобразования статистической суммы из ее стандартного многомерного интегрального представления по степеням свободы частиц в функциональное интегральное представление по вспомогательной полевой функции, используя либо преобразование Хаббарда – Стратоновича, либо дельта-функциональное преобразование. Компьютерное моделирование на основе теорий поля полимеров показали полезные результаты, например, для расчета структур и свойств растворов полимеров (Baeurle 2007, Schmid 1998), расплавов полимеров (Schmid 1998, Matsen 2002, Fredrickson 2002) и термопластов (Baeurle 2006). .
Канонический ансамбль
Частичное представление канонической статистической суммы
Стандартная модель континуума гибких полимеров, введенная Эдвардсом (Edwards 1965), рассматривает раствор, состоящий из линейные монодисперсные гомополимеры как система крупнозернистых полимеров, в которой статистическая механика цепей описывается моделью непрерывной гауссовой нити (Baeurle 2007), а растворитель учитывается неявно. Модель гауссовой нити можно рассматривать как непрерывный предел модели дискретной гауссовой цепи, в которой полимеры описываются как непрерывные линейно упругие нити. Каноническая статистическая сумма такой системы при обратной температуре и заключен в объем , можно выразить как
где является потенциалом средней силы , заданной,
представляющие опосредованные растворителем несвязанные взаимодействия между сегментами, в то время как представляет собой гармоническую энергию связи цепей. Последний энергетический вклад можно сформулировать как
где - длина статистического сегмента и индекс полимеризации.
Теоретико-полевые преобразования
Чтобы вывести основное теоретико-полевое представление канонической статистической суммы, в дальнейшем вводится оператор сегментной плотности полимерной системы
Используя это определение, можно переписать уравнение. (2) как
Затем модель преобразуется в теорию поля, используя преобразование Хаббарда-Стратоновича или дельта-функциональное преобразование.
где это функциональный и - дельта-функционал, задаваемый формулой
с участием представляющий функцию вспомогательного поля. Здесь отметим, что расширение функции поля в ряд Фурье означает, что периодические граничные условия применяются во всех направлениях и что-векторами обозначены векторы обратной решетки сверхъячейки.
Основное теоретико-полевое представление канонической статистической суммы
Используя уравнения. Используя уравнения (3), (4) и (5), мы можем преобразовать каноническую статистическую сумму в уравнение (2). (1) в теоретико-полевом представлении, что приводит к
где
можно интерпретировать как статистическую сумму для идеального газа невзаимодействующих полимеров и
- интеграл по траекториям свободного полимера в нулевом поле с упругой энергией
В последнем уравнении невозмущенный радиус вращения цепи . Более того, в формуле. (6) статистическая сумма одного полимера, подверженного действию поля, дан кем-то
Большой канонический ансамбль
Основное теоретико-полевое представление большой канонической статистической суммы
Чтобы вывести большую каноническую статистическую сумму, мы используем ее стандартную термодинамическую связь с канонической статистической суммой, заданной формулой
где химический потенциал и дается формулой. (6). Суммирование дает теоретико-полевое представление большой канонической статистической суммы,
где
это великое каноническое действие с определяется формулой. (8) и постоянная
Кроме того, параметр, связанный с химическим потенциалом, имеет вид
где обеспечивается формулой. (7).
Приближение среднего поля
Стандартная стратегия приближения для теорий поля полимеров - приближение среднего поля (МП), которое заключается в замене члена взаимодействия многих тел в действии на член, в котором все тела системы взаимодействуют со средним эффективным полем. Этот подход сводит любую проблему с несколькими телами к эффективной задаче с одним телом, предполагая, что интеграл статистической суммы модели определяется конфигурацией одного поля. Основное преимущество решения проблем с помощью приближения МП или его численной реализации, обычно называемой теорией самосогласованного поля (ССПП), заключается в том, что оно часто дает полезные сведения о свойствах и поведении сложных систем многих тел при относительно низкие вычислительные затраты. Успешные применения этой стратегии приближения можно найти для различных систем полимеров и сложных жидкостей, таких как, например, сильно сегрегированные блок-сополимеры с высокой молекулярной массой, высококонцентрированные растворы нейтральных полимеров или высококонцентрированные растворы блок- полиэлектролитов (PE) (Schmid 1998, Matsen 2002, Фредриксон 2002). Однако существует множество случаев, когда SCFT дает неточные или даже качественно неверные результаты (Baeurle 2006a). К ним относятся растворы нейтрального полимера или полиэлектролита в разбавленных и полуразбавленных режимах концентрации, блок-сополимеры вблизи их перехода порядок-беспорядок, смеси полимеров вблизи их фазовых переходов и т. Д. В таких ситуациях интеграл статистической суммы, определяющий теоретико-полевую модель, не полностью определяется Конфигурация одиночного МП и далекие от него конфигурации поля могут внести важный вклад, который требует использования более сложных методов расчета, выходящих за пределы уровня приближения МП.
Поправки высшего порядка
Одна из возможностей решить эту проблему - вычислить поправки более высокого порядка к приближению МП. Цончев и др. разработал такую стратегию, включающую поправки на флуктуации ведущего (однопетлевого) порядка, что позволило по-новому взглянуть на физику замкнутых растворов ПЭ (Цончев, 1999). Однако в ситуациях, когда приближение МП является плохим, для получения желаемой точности необходимо много требующих вычислений поправок к интегралу более высокого порядка.
Техники перенормировки
Альтернативный теоретический инструмент для решения проблем сильных флуктуаций, возникающих в теориях поля, был предоставлен в конце 1940-х годов концепцией перенормировки , которая первоначально была разработана для вычисления функциональных интегралов, возникающих в квантовых теориях поля (КТП). В КТП стандартная стратегия приближения заключается в расширении функциональных интегралов в степенной ряд по константе связи с использованием теории возмущений . К сожалению, обычно большинство членов разложения оказываются бесконечными, что делает такие вычисления невозможными ( Ширков 2001). Способ убрать бесконечности из КТП - использовать концепцию перенормировки (Baeurle 2007). Он в основном состоит в замене голых значений параметров связи, таких как, например, электрических зарядов или масс, перенормированными параметрами связи и требования, чтобы физические величины не изменялись при этом преобразовании, что приводит к конечным членам в разложении возмущений. Простую физическую картину процедуры перенормировки можно составить на примере классического электрического заряда:, вставленный в поляризуемую среду, например, в раствор электролита. На расстоянии от заряда, обусловленного поляризацией среды, ее кулоновское поле будет эффективно зависеть от функции , т. е. эффективный (перенормированный) заряд вместо затравочного электрического заряда, . В начале 1970-х годов К.Г. Уилсон стал пионером в силе концепций перенормировки, разработав формализм теории ренормгруппы (РГ) для исследования критических явлений статистических систем (Wilson, 1971).
Теория ренормгруппы
Теория RG использует серию преобразований RG, каждое из которых состоит из шага грубой зернистости, за которым следует изменение масштаба (Wilson 1974). В случае статистико-механических задач эти шаги реализуются путем последовательного исключения и изменения масштаба степеней свободы в сумме разбиений или интеграле, который определяет рассматриваемую модель. Де Женн использовал эту стратегию, чтобы установить аналогию между поведением классической векторной модели ферромагнетизма с нулевыми компонентами вблизи фазового перехода и случайным блужданием полимерной цепи бесконечной длины по решетке с самоизбеганием, чтобы вычислить исключенный объем полимера. экспоненты (de Gennes, 1972). Адаптация этой концепции к теоретико-полевым функциональным интегралам подразумевает систематическое изучение того, как модель теории поля изменяется при удалении и изменении масштаба определенного числа степеней свободы из интеграла статистической суммы (Wilson 1974).
Перенормировка Хартри
Альтернативный подход известен как приближение Хартри или самосогласованное однопетлевое приближение (Amit 1984). Он использует поправки на гауссовские флуктуации длявклад МП-порядка для перенормировки параметров модели и самосогласованного извлечения доминирующего масштаба длин флуктуаций концентрации в режимах критической концентрации.
Перенормировка головастика
В более поздней работе Ефимов и Ноговицин показали, что альтернативная техника перенормировки, происходящая из КТП, основанная на концепции перенормировки головастика , может быть очень эффективным подходом для вычисления функциональных интегралов, возникающих в статистической механике классических систем многих частиц (Ефимов 1996) . Они продемонстрировали, что основной вклад в классические интегралы статистической суммы вносят диаграммы Фейнмана типа головастиков низкого порядка , которые учитывают расходящиеся вклады из-за самодействия частиц . Процедура перенормировки, выполняемая в этом подходе, влияет на вклад самодействия заряда (например, электрона или иона) в результате статической поляризации, индуцированной в вакууме из-за присутствия этого заряда (Baeurle 2007). Как свидетельствуют Ефимов и Ганболд в более ранней работе (Ефимов, 1991), процедура перенормировки головастика может быть очень эффективно использована для устранения расхождений в действии основного теоретико-полевого представления статистической суммы и приводит к альтернативному функциональному интегралу представление, называемое гауссовским эквивалентным представлением (GER). Они показали, что эта процедура обеспечивает функциональные интегралы со значительно улучшенными свойствами сходимости для аналитических вычислений возмущений. В последующих работах Baeurle et al. разработали эффективные недорогие методы аппроксимации, основанные на процедуре перенормировки головастика, которые показали полезные результаты для прототипов полимеров и растворов полиэтилена (Baeurle 2006a, Baeurle 2006b, Baeurle 2007a).
Численное моделирование
Другой возможностью является использование алгоритмов Монте-Карло (МК) и выборка полного интеграла статистической суммы в теоретико-полевой формулировке. Полученная процедура называется теоретико-полевым моделированием полимеров . Однако в недавней работе Баерле продемонстрировал, что выборка MC в сочетании с базовым теоретико-полевым представлением неосуществима из-за так называемой проблемы числового знака (Baeurle 2002). Трудность связана со сложным и колеблющимся характером получаемой функции распределения, что приводит к плохой статистической сходимости средних по ансамблю требуемых термодинамических и структурных величин. В таких случаях необходимы специальные аналитические и численные методы для ускорения статистической сходимости (Baeurle 2003, Baeurle 2003a, Baeurle 2004).
Среднее представление поля
Чтобы сделать методологию пригодной для вычислений, Бэрле предложил сместить контур интегрирования интеграла статистической суммы через однородное решение МП, используя интегральную теорему Коши , обеспечивая его так называемое представление среднего поля . Эта стратегия ранее успешно применялась Baer et al. в теоретико-полевых расчетах электронной структуры (Baer 1998). Баерле смог продемонстрировать, что этот метод обеспечивает значительное ускорение статистической сходимости средних значений ансамбля в процедуре выборки МС (Baeurle 2002, Baeurle 2002a).
Гауссово эквивалентное представление
В последующих работах Baeurle et al. (Baeurle 2002, Baeurle 2002a, Baeurle 2003, Baeurle 2003a, Baeurle 2004) применили концепцию перенормировки головастика, приводящую к гауссовскому эквивалентному представлению интеграла статистической суммы, в сочетании с передовыми методами MC в большом каноническом ансамбле. Они могли убедительно продемонстрировать, что эта стратегия обеспечивает дальнейшее повышение статистической конвергенции желаемых средних значений по ансамблю (Baeurle 2002).
Рекомендации
- Baeurle, SA; Ноговицин Е.А. (2007). «Сложные законы масштабирования гибких полиэлектролитных растворов с эффективными концепциями перенормировки». Полимер . 48 (16): 4883. DOI : 10.1016 / j.polymer.2007.05.080 .
- Шмид, Ф. (1998). «Теории самосогласованного поля для сложных жидкостей». J. Phys .: Condens. Материя . 10 (37): 8105. arXiv : cond-mat / 9806277 . Bibcode : 1998JPCM ... 10.8105S . DOI : 10.1088 / 0953-8984 / 10/37/002 .
- Мацен, MW (2002). «Стандартная модель Гаусса для расплавов блок-сополимеров». J. Phys .: Condens. Материя . 14 (2): R21. Bibcode : 2002JPCM ... 14R..21M . DOI : 10.1088 / 0953-8984 / 14/2/201 .
- Fredrickson, GH; Ganesan, V .; Дроле, Ф. (2002). "Теоретико-полевые методы компьютерного моделирования полимеров и сложных жидкостей". Макромолекулы . 35 : 16. Bibcode : 2002MaMol..35 ... 16F . DOI : 10.1021 / ma011515t .
- Baeurle, SA; Усами, Т .; Гусев, А.А. (2006). «Новый подход многомасштабного моделирования для прогнозирования механических свойств полимерных наноматериалов». Полимер . 47 (26): 8604. DOI : 10.1016 / j.polymer.2006.10.017 .
- Эдвардс, С.Ф. (1965). «Статистическая механика полимеров с исключенным объемом». Proc. Phys. Soc . 85 (4): 613. Bibcode : 1965PPS .... 85..613E . DOI : 10.1088 / 0370-1328 / 85/4/301 .
- Baeurle, SA; Ефимов, Г.В. Ноговицин, Е.А. (2006а). «Расчет теории поля за пределами уровня среднего поля». Europhys. Lett . 75 (3): 378. Bibcode : 2006EL ..... 75..378B . DOI : 10,1209 / EPL / i2006-10133-6 .
- Цончев, С .; Коулсон, РД; Дункан, А. (1999). "Статистическая механика заряженных полимеров в растворах электролитов: подход теории поля решетки". Phys. Rev. E . 60 (4): 4257. arXiv : cond-mat / 9902325 . Bibcode : 1999PhRvE..60.4257T . DOI : 10.1103 / PhysRevE.60.4257 .
- Ширков, Д.В. (2001). «Пятьдесят лет ренормгруппе» . ЦЕРН Курьер . 41 : 14.
- Уилсон, KG (1971). "Ренормализационная группа и критические явления. II. Фазово-пространственный анализ критического поведения" . Phys. Rev. B . 4 (9): 3184. Bibcode : 1971PhRvB ... 4.3184W . DOI : 10.1103 / PhysRevB.4.3184 .
- Wilson, KG; Когут Дж. (1974). «Ренормализационная группа и ε-разложение». Phys. Rep . 12 (2): 75. Bibcode : 1974PhR .... 12 ... 75W . DOI : 10.1016 / 0370-1573 (74) 90023-4 .
- де Женн, PG (1972). «Экспоненты для проблемы исключенного объема, полученные методом Вильсона». Phys. Lett . 38 А : 339.
- Амит, диджей (1984). «Теория поля, ренормгруппа и критические явления» . Сингапур, World Scientific . ISBN 9812561196.
- Ефимов, Г.В. Ноговицин Е.А. (1996). «Статистические суммы классических систем в гауссовском эквивалентном представлении функциональных интегралов». Physica . 234 : 506. Bibcode : 1996PhyA..234..506V . DOI : 10.1016 / S0378-4371 (96) 00279-8 .
- Ефимов, Г.В. Ганболд, Г. (1991). «Функциональные интегралы в режиме сильной связи и собственная энергия полярона». Physica Status Solidi . 168 : 165. Bibcode : 1991PSSBR.168..165E . DOI : 10.1002 / pssb.2221680116 . ЛВП : 10068/325205 .
- Baeurle, SA; Ефимов, Г.В. Ноговицин, Е.А. (2006б). «О новой теории самосогласованного поля для канонического ансамбля». J. Chem. Phys . 124 (22): 224110. Bibcode : 2006JChPh.124v4110B . DOI : 10.1063 / 1.2204913 . PMID 16784266 .
- Baeurle, SA; Шарло, М .; Ноговицин Е.А. (2007а). «Грандиозные канонические исследования прототипных моделей полиэлектролитов за пределами приближения среднего поля». Phys. Rev. E . 75 : 011804. Bibcode : 2007PhRvE..75a1804B . DOI : 10.1103 / PhysRevE.75.011804 .
- Baeurle, SA (2002). "Метод гауссовского эквивалентного представления: новый метод уменьшения знаковой проблемы функциональных интегральных методов". Phys. Rev. Lett . 89 (8): 080602. Bibcode : 2002PhRvL..89h0602B . DOI : 10.1103 / PhysRevLett.89.080602 . PMID 12190451 .
- Baeurle, SA (2003). «Расчет в рамках подхода вспомогательного поля». J. Comput. Phys . 184 (2): 540. Bibcode : 2003JCoPh.184..540B . DOI : 10.1016 / S0021-9991 (02) 00036-0 .
- Baeurle, SA (2003a). «Метод Монте-Карло вспомогательного поля со стационарной фазой: новая стратегия уменьшения проблемы знака в методологиях вспомогательного поля». Comput. Phys. Commun . 154 (2): 111. Bibcode : 2003CoPhC.154..111B . DOI : 10.1016 / S0010-4655 (03) 00284-4 .
- Baeurle, SA (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, SA; Martonak, R .; Парринелло, М. (2002a). «Теоретико-полевой подход к моделированию в классическом каноническом и большом каноническом ансамбле». J. Chem. Phys . 117 (7): 3027. Bibcode : 2002JChPh.117.3027B . DOI : 10.1063 / 1.1488587 .
Внешние ссылки
- Исследовательская группа Регенсбургского университета по теории и расчетам перспективных материалов