CC BY
47
УДК 539.192
КВАНТОВО-ХИМИЧЕСКИЙ РАСЧЕТ СТРУКТУРНЫХ И ТЕРМОДИНАМИЧЕСКИХ ПАРАМЕТРОВ КЛАСТЕРИЗАЦИИ НАСЫЩЕННЫХ КАРБОНОВЫХ КИСЛОТ НА МЕЖФАЗНОЙ ПОВЕРХНОСТИ ВОДА/ВОЗДУХ
Е.А.Беляева
QUANTUM-CHEMICAL CALCULATION OF SOME STRUCTURAL AND THERMODYNAMIC PARAMETERS OF CLUSTERIZATION OF SATURATED CARBOXYLIC ACIDS ON THE AIR/WATER INTERFACE
Е.А.Belyaeva
Донецкий национальный технический университет, Украина, eabelyaeva@mail.ru
В представленной работе в рамках квантово-химического полуэмпирического программного комплекса Mopac 2012 (РМ3 метод) рассчитаны структурные и термодинамические параметры кластеризации насыщенных карбоновых кислот на межфазной поверхности вода/воздух при образовании монослоя гексагональной структуры. Сопоставление с имеющимися экспериментальными параметрами показало хорошее соответствие рассчитанных и экспериментальных данных.
Ключевые слова: гексагональная кристаллическая решетка, насыщенные карбоновые кислоты, структурные и термодинамические параметры монослоя, поверхность раздела фаз вода/воздух
In this paper structural and thermodynamic parameters of clusterization of fatty carboxylic acids clusters of hexagonal structure at air/water interface were calculated within quantum-chemical semiempirical program complex Mopac 2012 (PM3 method). Comparison of the calculated parameters with the corresponding experimental parameters gave good agreement between them.
Keywords: hexagonal unit cell, fatty carboxylic acids, structural and thermodynamic parameters of the monolayer, air/water interface
Введение
Экспериментально установлено, что карбоновые кислоты образуют монослои гексагональной структуры на межфазной поверхности вода/воздух [15], причем длины сторон элементарной ячейки (а и Ь) гексагональной кристаллической решетки определяются следующим соотношением а = Ь'^3 [6-12].
В данной работе рассчитаны параметры кластеризации насыщенных карбоновых кислот, являющихся элементами структуры монослоя гексагональной структуры. Расчеты проводились в рамках полуэмпирического метода РМ3 (Морас 2012 [13]) для молекул с общей формулой СпН2п+1СООН (где п варьировали от 7 до 16). Как ранее [14], наличие
межфазной поверхности вода/воздух учитывали неявно, через ее растягивающее и ориентирующее действие.
Обсуждение результатов Структурные параметры
В результате проведения конформационного анализа для молекул мономеров насыщенных карбоновых кислот, который заключался в варьировании величины торсионного углов C2QO1H и С3С2С1О2 (см. рис.1) от 0° до 360°, показано, что энергетический минимум для угла С2С101Н соответствует 180°, тогда как углу С3С2С102 соответствуют два минимума: первый — от -75°до -112°, второй — от -35° до -40°.
Рис.1 Фрагмент молекулы насыщенной карбоновой кислоты, R — углеводородный радикал
Для того чтобы определить угол наклона молекул монослоя к нормали относительно поверхности раздела фаз одну из молекул, являющейся стороной элементарной ячейки монослоя (см. рис.2), передвигали относительно второй такой молекулы [19]. При варьировании угла наклона от -15° до 25° на кривой энергии Гиббса образования наблюдаются два минимума: 15°-16° и 1° относительно нормали к межфазной поверхности. В данной работе рассматривали параметры кластеризации для структур с наклоном 15-16° относительно нормали к межфазной поверхности.
На рис.2 представлена структура фрагмента монослоя карбоновых кислот, состоящего из девяти мономеров. Стрелками на рис.2 условно обозначено направление от атома С к атому Н в метильном фрагменте молекулы.
Параметры элементарной ячейки: ^ = 16°, а = 8—8,2 А; Ь = 4,2—4,5 А, соответствующие экспериментальные данные: ^ = 20°, а = 8,4—8,7 А и Ь = 4,9—5,0 А. Параметры оптимизированной структуры воспроизводят экспериментальные.
Термодинамические параметры
Расчет термодинамических параметров кластеризации мономеров и малых кластеров, являющихся элементами представленного монослоя, также проводили в рамках квантово-химического полуэмпирического метода РМ3 (Морас 2012).
Мономеры. Ранее в [14] уже были рассчитаны энтальпия и энергия Гиббса образования мономеров карбоновых кислот данной структуры (см. рис.1), а также их абсолютная энтропия (мономеры 1). В работе [14] расчет был проведен в рамках квантово-химического программного комплекса Морас 2000, и так как в данной работе все расчеты проводились в рамках Морас 2012, возникла необходимость сравнить вновь рассчитанные параметры с соответствующими значениями, полученными ранее в рамках Мо-рас 2000 [16], а также соответствующими экспериментальными данными. Также были рассчитаны термодинамические параметры образования мономер с величиной торсионного угла С3С2С1О2 - 35°, такие мономеры обозначены как мономеры 2. Результаты расчета приведены в табл.1.
Из табл. 1 видно, что параметры мономеров 1 и 2 близки, поэтому в последующих расчетах будут использоваться только характеристики мономеров 1. Результаты расчета в рамках программного комплекса Морас 2012 хорошо воспроизводят как параметры, рассчитанные ранее в Морас 2000, так и экспериментальные данные [14,17-19, 20].
При расчете абсолютной энтропии необходимо учитывать, что в рамках программного комплекса Морас 2012 не учитывается вклад от свободного вращения метиленовых фрагментов. Поэтому, как и ранее (см., например [21-24]), величину данной поправки определяли методом наименьших квадратов. Она составила 7,84 Дж/мольК, что согласуется с аналогичными поправками для спиртов, тиоспиртов и аминов [21-24]. В табл.1 приведены значения абсолютной энтропии и энергии Гиббса образования мономеров карбоновых кислот с учетом данной поправки.
Малые кластеры. На рис.3 представлен фрагмент монослоя карбоновых кислот (вид сверху), на котором выделено три типа межмолекулярных СН' ■ НС взаимодействий, образующихся между углеводородными радикалами. Данные типы взаимодействий условно обозначали как Ка, Кь и Кс (см. рис.3), по аналогии с [21,26,27], где они описаны более подробно. Также необходимо выделить 6 типов взаимодействия между карбоксильными группами, эти взаимодействия представлены на рис.4.
ь
Рис.2. Фрагмент монослоя насыщенных карбоновых кислот на межфазной поверхности вода/воздух: a — вид сверху; Ь фронтальный вид
Таблица 1
Термодинамические параметры ряда мономеров насыщенных карбоновых кислот
Молекула AH°298,mon, кДж/молЬ
Mopac 2000 Mopac 2012 Эксперимент [17,18]
Мономеры 1 Мономеры 1 Мономеры 2
C7H^COOH -535,52 -558,66 -558,64 -556,00
C8HnCOOH -558,20 -581,38 -581,35 -577,30
C9H19COOH -580,86 -604,10 -604,14 -594,30
C10H21COOH -603,54 -626,83 -626,83 -614,60
C11H23COOH -626,22 -649,56 -649,60 -640,00
C12H25COOH -648,90 -672,30 -672,31 -660,20
C13H27COOH -671,58 -695,03 -695,07 -683,00
C14H29COOH -694,26 -717,77 -717,77 -699,00
C15H31COOH -716,94 -740,49 -740,49 -723,00
C16H33COOH -763,25 -763,25 -743,00
S°298, mon, Дж/мольК
C7H15COOH 535,30 528,29 527,81 520,00
C8H17COOH 573,36 566,43 565,60 559,00
C9H19COOH 612,20 604,52 604,40 599,00
C10H21COOH 651,25 641,63 641,61 638,00
C11H23COOH 688,88 679,88 680,59 677,40
C12H25COOH 727,75 717,21 716,62 717,00
C13H27COOH 765,40 754,79 756,88 754,00
C14H29COOH 802,55 792,33 791,41 796,00
C15H31COOH 841,22 829,43 829,21 833,00
C16H33COOH 865,52 865,37 874,00
AG°298,mon, кДж/молЬ
C7H15COOH -331,58 -331,85 -331,69 -325,00
C8H17COOH -324,96 -325,33 -325,06 -317,00
C9H19COOH -318,58 -318,80 -318,81 -305,00
C10H21COOH -312,27 -311,99 -311,98 -296,63
C11H23COOH -305,54 -305,51 -305,75 -293,10
C12H25COOH -299,17 -298,77 -298,60 -284,50
C13H27COOH -292,44 -292,09 -292,75 -278,00
C14H29COOH -285,57 -285,41 -285,14 -266,00
C15H31COOH -279,14 -278,58 -278,52 -260,00
C16H33COOH -271,49 -271,45 -252,00
Димеры 1
Димеры 2
Рис.5. Оптимизированные структуры димеров насыщенных карбоновых кислот
Структуры димеров, термодинамические параметры которых были рассчитаны, приведены на рис.5.
Стрелками на рис.3 обозначено условное направление от атома углерода к кетонному атому ки-
слорода в карбоксильной группе. Вклады от этих взаимодействий различаются вследствие различий в их взаимной ориентации. Для описания монослоя нужно определить вклад в энтальпию и энтропию от каждого СН-НС взаимодействия и каждого взаимодействия между функциональными группами. Для этого в структуре кластера (см. рис.3) были выделены структуры димеров, тримеров и тетрамеров, содержащие в своем составе все рассматриваемые взаимодействия.
Стрелками на рис.5 показаны межмолекулярные СН-НС взаимодействия (взаимодействия «а»-типа для димеров 1 и «Ь»-типа — для димеров 2).
Рассмотрены структуры пяти тримеров (рис.6), двух тетрамеров от (рис.7). Для них рассчитаны термодинамические параметры образования и кластеризации (см. табл.2). Термодинамические параметры кластеризации рассчитывали по следующим формулам: АНС1298 =
= ДН°т,98 - тДН°-
298,тош С1
дбс
= Б
т,298
- тБ'
298,топ,
т,298
ДНС1298 = ДНС1298 - Т^ДБС1298, где ДН°т,298 — энтальпия образования кластера, кДж/моль; Б°т,298 — абсолютная энтропия кластера, Дж/мольК; Т — абсолютная температура, К; т — число мономеров в кластере.
Рис.6. Оптимизированные структуры тримеров насыщенных карбоновых кислот
Таблица 2
Термодинамические параметры димеров, тримеров, тетрамеров насыщенных карбоновых кислот
ДH°m,298, S°m,298, ДНС1298, ДБС1298, ДGC1298,
Молекула кДж/моль Дж/моль-К кДж/моль Дж/моль-К кДж/моль
Димеры 1
С7Н!зСООН -1149,19 792,45 -32,89 -183,43 21,78
С8НпСООН -1202,01 851,46 -40,38 -188,56 15,81
СДдТООН -1249,92 901,92 -42,94 -200,36 16,76
С10Н21СООН -1303,00 952,03 -50,68 -215,86 13,65
С11Н23СООН -1350,91 1002,77 -53,46 -224,26 13,37
С12Н25СООН -1404,04 1052,24 -61,04 -242,08 11,10
С13Н27СООН -1451,90 1103,82 -63,53 -249,93 10,95
С14Н29СООН -1505,09 1153,16 -71,43 -267,15 8,18
С15Н31СООН -1552,98 1206,51 -73,99 -273,99 7,66
С16Н33СООН
Продолжение табл. 2
Димеры 1
С7Н!зСООН -1160,45 708,39 -43,13 -238,43 27,92
С8НПС00Н -1212,05 749,55 -49,29 -257,87 27,55
С9Н19СООН -1263,68 792,77 -55,48 -275,15 26,52
С10Н21СООН -1315,30 832,62 -61,64 -293,84 25,93
С11Н23СООН -1366,93 874,34 -67,81 -312,94 25,44
С12Н25СООН -1418,56 915,94 -73,96 -330,32 24,47
С13Н27СООН -1470,20 957,81 -80,14 -347,93 23,54
С14Н29СООН -1521,84 999,40 -86,30 -365,74 22,69
С15Н31СООН -1573,49 1040,19 -92,51 -383,47 21,77
С16Н33СООН -1625,13 1080,32 -98,63 -403,70 21,67
Тримеры 1
С7Н!5СООН -1769,75 995,44 -93,77 -424,79 32,82
С8Н!7СООН -1849,29 1062,60 -105,15 -448,53 28,51
С9Н19СООН -1935,41 1126,30 -123,11 -475,58 18,62
С10Н21СООН -2015,07 1186,26 -134,58 -503,43 15,44
С11Н23СООН -2101,17 1251,03 -152,49 -529,89 5,42
С12Н25СООН -2180,94 1311,76 -164,04 -557,63 2,13
С13Н27СООН -2267,00 1378,18 -181,91 -580,43 -8,94
С14Н29СООН -2346,82 1437,75 -193,51 -609,96 -11,74
С15Н31СООН -2432,87 1504,80 -211,40 -630,69 -23,46
С16Н33СООН -2512,71 1562,84 -222,96 -663,19 -25,32
Тримеры 2
С7Н!5СООН -1764,71 1008,17 -88,73 -412,06 34,06
С8Н!7СООН -1844,72 1072,27 -100,58 -438,86 30,20
С9Н19СООН -1930,54 1135,81 -118,24 -466,07 20,64
С10Н21СООН -2010,38 1199,00 -129,89 -490,69 16,34
С11Н23СООН -2096,42 1260,71 -147,74 -520,21 7,28
С12Н25СООН -2176,17 1322,45 -159,27 -546,94 3,71
С13Н27СООН -2262,32 1386,25 -177,23 -572,36 -6,66
С14Н29СООН -2342,01 1446,00 -188,70 -601,71 -9,39
С15Н31СООН -2428,23 1510,72 -206,76 -624,77 -20,58
С16Н33СООН -2507,87 1568,83 -218,12 -657,20 -22,27
Тримеры 3
С7Н!5СООН -1757,96 941,86 -81,98 -478,37 60,58
С8Н!7СООН -1838,42 996,65 -94,28 -514,48 59,03
С9Н19СООН -1918,72 1054,05 -106,42 -547,83 56,83
С10Н21СООН -1999,21 1107,07 -118,72 -582,62 54,90
С11Н23СООН -2079,59 1162,29 -130,91 -618,63 53,44
С12Н25СООН -2160,06 1218,11 -143,16 -651,28 50,92
С13Н27СООН -2240,49 1272,87 -155,40 -685,74 48,95
С14Н29СООН -2320,96 1327,83 -167,65 -719,88 46,88
С15Н31СООН -2401,42 1383,71 -179,95 -751,78 44,09
С16Н33СООН -2481,87 1438,03 -192,12 -788,00 42,70
Тримеры 4
С7Н15СООН -1769,04 934,62 -93,06 -485,61 51,65
С8Н17СООН -1846,88 998,77 -102,74 -512,36 49,95
С9Н19СООН -1941,59 1026,13 -129,29 -575,75 42,28
С10Н21СООН -2010,65 1113,18 -130,16 -576,51 41,64
С11Н23СООН -2096,44 1165,37 -147,76 -615,55 35,68
С12Н25СООН -2174,46 1226,73 -157,56 -642,66 33,96
С13Н27СООН -2260,22 1278,69 -175,13 -679,92 27,48
С14Н29СООН -2338,29 1339,83 -184,98 -707,88 25,97
С15Н31СООН -2424,03 1389,99 -202,56 -745,50 19,59
С16Н33СООН -2502,12 1452,06 -212,37 -773,97 18,27
Продолжение табл. 2
Тримеры 5
С7Н15СООН -1764,52 990,38 -88,54 -429,85 39,56
С8Н17СООН -1841,94 1066,87 -97,80 -444,26 34,59
С9Н19СООН -1928,14 1115,08 -115,84 -486,80 29,22
С10Н21СООН -2005,57 1190,92 -125,08 -498,77 23,56
С11Н23СООН -2091,87 1239,13 -143,19 -541,79 18,26
С12Н25СООН -2169,27 1315,69 -152,37 -553,70 12,63
С13Н27СООН -2255,64 1360,51 -170,55 -598,10 7,68
С14Н29СООН -2333,02 1441,24 -179,71 -606,47 1,02
С15Н31СООН -2419,44 1484,60 -197,97 -650,89 -4,00
С16Н33СООН -2496,80 1561,95 -207,05 -664,08 -9,16
Тетрамеры 1
С7Н15СООН -2385,00 1179,07 -150,36 -714,57 62,58
С8Н17СООН -2494,69 1263,82 -169,17 -751,02 54,63
С9Н19СООН -2613,39 1326,88 -196,99 -808,96 44,08
С10Н21СООН -2722,99 1410,01 -215,67 -842,91 35,51
С11Н23СООН -2841,84 1472,74 -243,60 -901,82 25,14
С12Н25СООН -2951,48 1557,47 -262,28 -935,05 16,37
С13Н27СООН -3070,34 1619,30 -290,22 -992,18 5,45
С14Н29СООН -3180,00 1705,67 -308,92 -1024,61 -3,59
С15Н31СООН -3298,86 1766,30 -336,90 -1081,02 -14,75
С16Н33СООН -3408,55 1850,84 -355,55 -1117,20 -22,62
Тетрамеры 2
С7Н15СООН -2379,60 1203,79 -144,96 -689,85 60,61
С8Н17СООН -2489,55 1288,46 -164,03 -726,38 52,43
С9Н19СООН -2607,99 1348,15 -191,59 -787,69 43,15
С10Н21СООН -2717,82 1434,91 -210,50 -818,01 33,27
С11Н23СООН -2836,45 1495,40 -238,21 -879,16 23,77
С12Н25СООН -2946,23 1579,55 -257,03 -912,97 15,03
С13Н27СООН -3064,95 1641,70 -284,83 -969,78 4,16
С14Н29СООН -3174,69 1724,98 -303,61 -1005,30 -4,03
С15Н31СООН -3293,47 1786,17 -331,51 -1061,15 -15,29
С16Н33СООН -3403,19 1870,18 -350,19 -1097,86 -23,02
Тетрамеры 1
Тетрамеры 2
Рис.7. Оптимизированные структуры тетрамеров насыщенных карбоновых кислот
На основе рассчитанных значений термодинамических параметров кластеризации (см. табл.1 и табл.2) были построены корреляционные зависимости энтальпии и энтропии кластеризации от числа межмолекулярных СН-НС взаимодействий между углеводородными радикалами и числа межмолекулярных взаимодействий между карбоксильными группами.
ДНС1 = -(9,66 ± 0,06) Ка - (6,00 ± 0,03>КЬ -- (7,22 ± 0,39)-Кс - (3,55 ± 0,73>п - (0,17 ± 0,40)х х(п2 +п3) + (2,43 ± 0,48)^п4 - (5,34 ± 0,52>п5 -
- (2,68 ± 0,52)^п6, кДж/моль (5)
(R = 0,99995; Б = 1,81 кДж/моль; N = 90)
ДБС1 = -(17,48 ± 0,56)^Ка - (16,97 ± 0,68>КЬ -
- (14,54 ± 0,39)^Кс - (130,29 ± 6,15>п -
- (122,98 ± 7,95)^п2 - (75,87 ± 5,93)-п3 -
- (62,50 ± 5,94)-п4 - (84,81 ± 4,44>п5 -
- (51,48 ± 4,44)^п6, Дж/моль-К (6)
(R = 0,9997; Б = 15,40 Дж/молыК; N = 90) где Ка — число парных межмолекулярных СН'НС взаимодействий «а»-типа; Кь — число парных межмолекулярных СН'НС взаимодействий «Ь»-типа; Кс — число парных межмолекулярных СН'НС взаимодействий «с»-типа; п1-п6 — число межмолекулярных взаимодействий) между карбоксильными группами (при различной ориентации этих групп (см. рис.4).
Значения угловых коэффициентов корреляционных зависимостей (5) и (6) совпадают с рассчитанными ранее аналогичными параметрами для жирных карбоновых кислот других классов замещенных ал-канов, рассчитанных в рамках программного комплекса Морас 2000.
Воспользовавшись известной
формулой
ДGC1298 = ДНС1298 - Т^С1298, получена корреляционная зависимость энергии Гиббса кластеризации насыщенных карбоновых кислот от числа межмолекулярных взаимодействий в кластере:
ДGC1298 = -(4,45 ± 0,23)-Ка - (0,94 ± 0,23>КЬ -- (2,89 ± 0Д6Ж + (35,28 ± 2,56>п + + (36,47 ± 2,64)п + (22,43 ± 2,17)п + + (21,95 ± 1,98>щ + (19,93 ± 1,84)п + + (12,66 ± 1,84)^п6, кДж/моль (7)
Видим, что на одно парное «а»-взаимодействие приходится 4,45 кДж/моль, на «Ь»-взаимодействие — 0,94 кДж/моль, а на «с»-взаимодействие — 2,89 кДж/моль, то СН-НС взаимодействия «а»-типа наиболее энергетически выгодны. Это согласуется с результатами наших предыдущих работ [14, 21-26].
Выводы
Структурные параметры монослоя насыщенных карбоновых кислот «гексагональной» структуры, полученные с помощью квантово-химического полуэмпирического метода РМ3 (Морас 2012), хорошо воспроизводят соответствующие экспериментальные параметры. Так, длины сторон элементарной ячейки составили а = 8—8,2 А; Ь = 4,2—4,5 А (а = 8,4—8,7 А и Ь = 4,9—5,0 А согласно экспериментальным данным). Угол наклона молекул монослоя относительно нормали к поверхности раздела фаз равен 16°, тогда как согласно имеющимся экспериментальным данным эта величина составляет 20°.
Рассчитанные термодинамические параметры образования мономеров карбоновых кислот с высокой степенью точности воспроизводят соответствующие экспериментальные данные. Так, отклонение рассчитанных энтальпии образования и абсолютной энтропии от экспериментальных параметров составляют соответственно 4,74 кДж/моль и 5,51 Дж/мольК.
Рассчитанные термодинамические параметры кластеризации малых кластеров в дальнейшем будут положены в основу общей аддитивной схемы, позволяющей рассчитать соответствующие термодинамические параметры кластеризации кластеров насыщенных карбоновых кислот любой размерности вплоть до монослоев.
1. Kenn R.M., Biihm C., Bib A.M., Peterson I.R., Mohwald H. Mesophases and Crystalline Phases In Fatty Acid Monolayers // Phys. Chem. 1991. V.95. Issue 5. P.2092.
2. Friedenberg M.C., Fuller G.G., Frank C.W., Robertson C.R. Direct Visualization of Flow-Induced Anisotropy in a Fatty Acid Monolayer // Langmuir. 1996. V.12. P.1594.
3. Maruyama T., Lauger J., Fuller G. G., Frank C. W., Robertson C. R. Orientation in a Fatty Acid Monolayer: Effect of Flow Type // Langmuir. 1998. V.14. P.1836.
4. Riviere S., Henon S., Meunier J. Textures and phase transitions in Langmuir monolayers of fatty acids. A comparative Brew-ster angle microscope and polarized fluorescence microscope study // J. Chem. Phys. 1994. V.101. Issue 11. P.10045.
5. Overbeck G.A., Mobius D.A. New Phase in the Generalized Phase Diagram of Monolayer Films of Long-chain Fatty Acids // J. Phys. Chem. 1993. V.97. P.7999.
6. Knobler C.M. Recent developments in the study interface of monolayers at the air-water // Advance in Chemical Physics. 1990. V.LXXVI. P.397.
7. Honig D., Overbeck G.A., Mobius D. Morphology of Penta-decanoic Acid Monolayers at the Air/Water Interface Studied by BAM // Adv Mater. 1992. V.4. P.419.
8. Schwartz D.K., Tsao M.-W., Knobler C.M. Domain morphology in a two-dimensional anisotropic mesophase: Cusps and boojum textures in a Langmuir monolayer // J. Chern. Phys. 1994. V.101. Issue 9. P.8258.
9. Johann R, Brezesinski G, Vollhardt D, Mohwald H. The effect of headgroup interactions on structure and morphology of arachidic acid monolayers // J. Phys. Chem. B. 2001. V.105. P.2958.
10. Weidemann G., Brezesinski G., Vollhardt D., Bringezu F., Mei-jere K. De, Mohwald H. Comparing Molecular Packing and Textures of Langmuir Monolayers of Fatty Acids and Their Methyl and Ethyl Esters // J. Phys. Chem. B. 1998. V.102. P.148.
11. Peterson I.R., Brezesinski G., Struth B., Scalas E. Grazing-Incidene X-ray Diffraction Study of Octadecanoic Acid Monolayers // J. Phys. Chem. B 1998. V.102. P.9437.
12. Durbin M.K., Malik A., Ghaskadvi R., Shih M.C., Zschack P., Dutta P. X-ray Diffraction Study of a Recently Identified Phase Transition in Fatty Acid Langmuir Monolayers // J. Phys. Chem. 1994. V.98. P.1753.
13. James J.P. Stewart, Stewart Computational Chemistry, Colorado Springs, CO, USA, HTTP://OpenMOPAC.net
14. Vysotsky Yu.B., Muratov D.V., Boldyreva F.L., Fainerman V.B., Vollhardt D., Miller R. Quantum Chemical Analysis of the Thermodynamics of 2D Cluster Formation of n-Carboxylic Acids at the Air/Water Interface // J. Phys. Chem. B. 2006. V.110. P.4717.
15. Vysotsky Yu.B., Fomina E.S., Belyaeva E.A., Vollhardt D., Fainerman V.B., Miller R. Quantum Chemical Analysis of the Thermodynamics of 2D Cluster Formation of Aliphatic Amides at the Air/Water Interface // J. Phys. Chem. C. 2012, V.116. Issue 50. P.26358.
16. Stewart J. J. P. MOPAC 2000.00 Manual; Fujitsu Limited: Tokyo, Japan, 1999.
17. Daubert T.E., Dannerm R.P., Sibul H.M., Stebbins C.C. Physical and thermodynamic properties of pure chemicals: Data Compilation, Part 1 - Part 5; Taylor & Francis: Pennsylvania, 1998. 9860 p.
18. Zwolinski B.J, Wilhoit R. Heats of formation and heats of combustion in «American Institute of physics handbook». 3-rd edition by D.E.Gray. Mc. Graw-Hill. N.-Y., 1972. P.4-342.
19. Vysotsky Yu.B., Belyaeva E.A., Fainerman V.B., Aksenenko E.V., Vollhardt D., Miller R. Quantum Chemical Analysis of the Thermodynamics of 2-Dimensional Cluster Formation of Alkylamines at the Air/Water Interface // J. Phys. Chem. C. 2007. V.111. P.15342-15349.
20. Vysotsky Yu.B., Belyaeva E.A., Fainerman V.B., Vollhardt D., Aksenenko E.V., Miller R. Thermodynamics of the Clusterization Process of Cis Isomers of Unsaturated Fatty Acids at the Air/Water Interface // J. Phys. Chem. B. 2009. V.113. P.4347-4359.
21. Vysotsky Yu.B., Bryantsev V.S., Fainerman V.B., Vollhardt D., Miller R. Quantum chemical semi-empirical approach to the thermodynamic characteristics of oligomers and large aggregates of alcohols at the water/air interface // Colloid and Surfaces A. 2002. V.209. P.1-14.
22. Vysotsky Yu. B., Belyaeva E. A., Fainerman V. B., Vollhardt
D., Miller R. Quantum Chemical Analysis of Thermodynamics of 2D Cluster Formation of n-Thioalcohols at the Air/Water Interface J. Phys. Chem., C 2007, 111, 5374-5381.
23. Vysotsky Yu.B., Belyaeva E.A., Fainerman V.B., Aksenenko
E.V., Vollhardt D., Miller R. Quantum Chemical Analysis of the Thermo-dynamics of 2-Dimensional Cluster Formation of Alkylamines at the Air/Water Interface // J. Phys. Chem. C.
2007. V.111. P.15342-15349.
24. Vysotsky Yu.B., Belyaeva E.A., Fainerman V.B., Vollhardt
D., Aksenenko E.V., Miller R. Thermodynamics of the Clusterization Process of Cis Isomers of Unsaturated Fatty Acids at the Air/Water Interface // J. Phys. Chem. B. 2009. V.113. P.4347-4359.
25. Vysotsky Yu.B., Bryantsev V.S., Fainerman V.B. Thermodynamics of 2D Cluster Formation of Odd n-Alcohols at the Air/Water Interface // J. Phys. Chem. B. 2002. V.106. P.11285-11294.
26. Vysotsky Yu.B., Belyaeva E.A., Vollhardt D., Aksenenko
E.V., Miller R. Simplified method of the quantum chemical analysis for determination of thermodynamic parameters of 2D cluster formation of amphiphilic compounds at the air/water interface // Journal of Colloid and Interface Science.
2008. V.326. P.339-346.
