О КОНФОРМАЦИОННОЙ ИЗОМЕРИИ ЦИКЛОТРИФОСФАТА НАТРИЯ
Ниязбекова Актоты Болатовна
к.х.н., доцент кафедры химии и химической технологии
Шакиров Тимур Аскарович
старший преподаватель кафедры химии и химической технологии
Западно-Казахстанский аграрно-технический университет им. Жангир хана
г.Уральск, Республика Казахстан
АННОТАЦИЯ
Проведены квантово-химические расчеты термодинамических характеристик конформационных изомеров циклотрифосфата натрия. Рассчитаны константы равновесия, количественное соотношение форм в равновесных системах и энергия активации при конформационных переходах. По результатам расчетов построена энергетическая диаграмма. Расчетные данные находятся в полном согласии с известными закономерностями.
ABSTRACT
Quantum-chemical calculations of the thermodynamic characteristics of conformational isomers of sodium cyclotriphosphate are carried out. Equilibrium constants, the proportion of forms in equilibrium systems, and the activation energy during conformational transitions are calculated. According to the results of calculations, an energy diagram was constructed. Calculated data are in complete agreement with known laws.
Ключевые слова: циклофосфаты, квантово-химические расчёты, конформеры, устойчивость, структура.
Keywords: cyclophosphates, quantum chemical calculations, conformers, stability, structure.
Экспериментальное изучение циклических фосфатов имеет большое практическое значение так как области их применения с каждым днем расширяются. Расширение знаний о пространственной структуре подобных соединений позволит объяснить поведение модифицированных полифосфатных систем и причины их эффективности в качестве удобрений, ингибиторов коррозии и ионообменников. В настоящем исследовании предпринята попытка с применением методов компьютерной химии проанализировать устойчивость различных конформационных изомеров циклотрифосфата натрия, имеющего на сегодняшний день наиболее широкое практическое значение [1,2,3].
Конформационная изомерия циклотрифосфата натрия
Каждый атом в молекуле всегда стремится, чтобы его валентные углы совпадали с углами между его связывающими орбиталями. Любое отклонение от «нормального» значения валентного угла сопровождается возникновением углового напряжения. Несвязанные атомы или группы стремятся занять такие положения, чтобы были обеспечены наиболее благоприятные диполь-дипольные взаимодействия: положения, в которых минимальны диполь-дипольные отталкивания и максимальны диполь-дипольные притяжения. Для того чтобы вывести наиболее устойчивую конформацию молекулы, необходимо, рассмотреть все возможные комбинации валентных и торсионных углов, и даже длин связей и выяснить, какая из комбинаций приведет к меньшему содержанию энергии.
Поскольку молекулы циклических фосфатов имеют кольцевое строение, они представляют собой в свободном состоянии довольно стабильные структуры, сохраняющиеся в том виде, в каком образовались в процессе синтеза. Известно, например, что циклотрифосфат в кристаллическом состоянии находится в виде «кресла» [4], а в растворе может параллельно существовать и в форме «ванны» [5]. Однако отсутствуют какие-либо достоверные данные о количественном соотношении этих форм.
Так как основа структуры молекулы циклотрифосфата представляет собой шестичленное кольцо, образованное чередующимися атомами фосфора и кислорода, то по аналогии с подобными циклическими соединениями, предполагается его существование в виде трех конформационных изомеров: «кресло», «ванна» и «твист» («искаженная ванна»).
Однако существует другой тип пространственных структур одного и того же соединения, которые не могут непосредственно переходить друг в друга по пространственным причинам, и сохраняются в ходе химических превращений в первозданном виде (как, например, оптические D- и L-изомеры органических соединений).
В случае циклических фосфатов речь идет о взаимном расположении атомов кислорода при двойной связи, относительно плоскости циклофосфатного кольца. При этом, структуру в которой все атомы кислорода при двойной связи, расположены по одну сторону кольца следует обозначить как α-, а структуру, в которой один атом кислорода при двойной связи располагается по другую сторону кольца – β (рис.1):
Рисунок 1. α- и β-формы циклотрифосфата натрия
Образование той или иной формы возможно только в процессе термического синтеза. Следовательно, целесообразно рассмотреть существование циклотрифосфата в виде шести различных пространственных структур.
Расчеты энергетических характеристик
Для оценки устойчивости рассматриваемых структур проведены квантово-химические расчеты полуэмпирическим методом РМ3. Первичная оптимизация проводилась методом молекулярной механики в силовом поле ММ+.
Расчет выполнен с применение программного пакета GAMESS [6]. Для удобства ввода структур молекул, формирования расчётного задания и представления полученной информации использовались программы Avogadro 1.2.0 и Facio 19.1.7.
При расчетах были оценены следующие характеристики: энтальпия образования конформеров и переходных форм, энергия активации при конформационных переходах, свободная энергия Гиббса, константа конформационного равновесия и количественные соотношения конформеров в равновесной системе.
Наиболее информативным показателем в такого рода исследованиях является свободная энергия Гиббса, поскольку только процессы, для которых ∆G<0 могут протекать самопроизвольно. К тому же, величина ∆G связана с константой равновесия K соотношением:
ΔG = – RT ln K
При Т = 298,15K различию в свободной энергии находящихся в равновесии изомеров ∆G = 0,5 ккал/моль соответствует процентное соотношение 30:70 в равновесной смеси. При ∆G = 1,0 ккал/моль — это соотношение составит уже 15:85, а при ∆G = 2,0 ккал/моль содержание менее устойчивого изомера составит лишь 3%. Это дает представление о том, какая точность определения энергии необходима для надежного предсказания такого рода соотношений. Тем не менее, для изомерных соединений при стандартных условиях величина T∆S невелика, и для оценки констант равновесия можно принять ∆G ≈ ∆H. При этом погрешность не превышает допустимую и расчет дает хорошее согласование с экспериментальными данными, но только при условии сохранения числа однотипных связей при конформационных переходах [7,8].
Расчетные значения энтальпии образования конформеров циклотрифосфата натрия представлены в таблице 1.
Таблица 1.
Энтальпия образования Н 0298,15 конформеров циклотрифосфата натрия
| Конформация | Н 0298,15, ккал/моль | |
| α-форма | β-форма | |
| «кресло» | -379.27 | -379.00 |
| «ванна» | -378.61 | -377.70 |
| «твист» | -366.78 | -365.13 |
По рассчитанным данным видно, что устойчивость конформеров (не зависимо от формы) уменьшается в ряду:
«кресло» ˃ «ванна» ˃ «твист».
Разница в потенциальной энергии показывает, что при 298,15К преобладающая часть молекул (60,56% α-формы и 83,81% β-формы) находятся в креслообразной форме, характеризующейся более низкой энергией. Эта конформация свободна не только от углового, но также и от торсионного напряжения. Форма лодки появляется в большей концентрации только при значительно более высоких температурах, когда приданная молекуле энергия позволяет ей перейти в состояние с высокой потенциальной энергией. Содержание же твист-конформеров в равновесной системе для α-формы не превышает 6∙10-6 %, а для β-формы 2∙10-6 %.
При расчетах принято, что конформационные переходы осуществляются через переходное состояние, энергетика которого также рассчитана [9]. Данные представлены в таблице 2.
Таблица 2.
Энергия переходных состояний при конформационных переходах циклотрифосфата
| Равновесная система | Еперех.сост., ккал/моль | |
| α-форма | β-форма | |
| «кресло» ↔ «ванна» | -374.68 | -376.30 |
| «ванна» ↔ «твист» | -364.71 | -363.07 |
| «кресло» ↔ «твист» | -359.76 | -360.07 |
Энергия активации (энергетический барьер) оценена как разница между энергией переходного и начального состояния системы [10]. Данные представлены в таблице 3.
Таблица 3.
Энергия активации Еакт конформационных переходов циклотрифосфата
| Равновесная система | Еакт., ккал/моль | |
| α-форма | β-форма | |
| «кресло» ↔ «ванна» | 4.58 | 2.70 |
| «ванна» ↔ «твист» | 7.02 | 5.06 |
| «кресло» ↔ «твист» | 14.63 | 13.90 |
Рассчитанные значения свободной энергии Гиббса и константы равновесия представлены в таблице 4
Таблица 4.
Свободная энергия Гиббса G0298,15 и константы равновесия при конформационных переходах циклотрифосфата (принято, что ∆G ≈ ∆H).
| Равновесная система | α-форма | β-форма | ||||
| — G 0298,15, ккал/моль | lnK | Kравн. | — G 0298,15, ккал/моль | lnK | Kравн. | |
| «кресло» ↔ «ванна» | 0.66 | -1.115 | 3.94∙10-1 | 1.29 | -2.182 | 1.62∙10-1 |
| «ванна» ↔ «твист» | 11.82 | -19.960 | 5.87∙10-8 | 12.57 | -21.226 | 2.04∙10-8 |
| «кресло» ↔ «твист» | 12.49 | -21.075 | 2.32∙10-8 | 13.87 | -23.408 | 3.31∙10-9 |
Расчеты показывают, что конформационные переходы легче протекают у β-формы, и в обоих случаях происходят по схеме:
«кресло» ↔ «ванна» ↔ «твист».
Переход непосредственно «кресло» ↔ «твист», через соответствующее переходное состояние, в обоих случаях энергетически не выгоден, так как для этого необходимо в 3 раза больше энергии, чем для перехода «кресло» ↔ «ванна». При этом рассчитанные значения полностью согласуются с законом Гесса:
| Н | 298,15 | + | Н | 298,15 | = | Н | 298,15 | |||||||
| кресло ванна | ванна твист | кресло твист | ||||||||||||
Расчетные значения изменения энтальпии при конформационных переходах представлены в таблице 5.
Таблица 5.
Изменение энтальпии Н0298,15 при конформационных переходах циклотрифосфата
| Конформационный переход | Н0298,15, ккал/моль | |
| α-форма | β-форма | |
| «кресло» ↔ «ванна» | 0.66 | 1.29 |
| «ванна» ↔ «твист» | 11.82 | 12.57 |
| «кресло» ↔ «твист» | 12.49 | 13.87 |
Как видно, все конформационные переходы являются эндотермическими процессами и для смещения равновесия в сторону прямой реакции требуется затрата энергии.
При этом относительное содержание молекул конформеров в равновесном состоянии равно:
«α-кресло» : «α-ванна» : «α-твист» = 2,62∙107 : 1,70∙107 : 1;
«β-кресло» : «β-ванна» : «β-твист» = 2,55∙1010 : 4,92∙109 : 1.
Наиболее наглядно полученные закономерности выглядят на соответствующей энергетической диаграмме (рис. 2). При этом внешний вид диаграммы для обеих форм циклотрифосфата идентичен. Отличия состоят только в абсолютных значениях рассчитанных энергетических параметров.
Рисунок 2. Энергетическая диаграмма переходов между конформациями циклотрифосфата
Можно также предположить, что конформация «ванна» представляет собой не реально существующий конформер, а лишь переходное состояние в равновесной системе «кресло» ↔ «твист». Однако, на энергетической кривой эта конформация не отвечает максимуму потенциальной энергии и, следовательно, данное предположение является ошибочным. Так же не возможен переход по схеме «кресло» ↔ «твист» ↔ «ванна».
Для сравнения, энергетический барьер внутреннего вращения в молекуле этана (между заторможенной D3d и заслоненной D3h конформацией) составляет 2,9 ккал/моль, а для перехода «кресло» «ванна» в молекуле циклогексана 5,3 ккал/моль [7].
Выводы:
1. существует две формы молекулы циклотрифосфата, α- и β-, между которыми взаимный переход невозможен;
2. обе формы могут существовать в виде трех конформационных изомеров: «кресло», «ванна» и «твист»;
3. устойчивость конформационных изомеров уменьшается в ряду «кресло» ˃ «ванна» ˃ «твист»;
4. взаимопереход конформеров осуществляется по схеме: «кресло» ↔ «ванна» ↔ «твист»;
5. конформация «ванна» не является переходным состоянием в равновесной системе«кресло» ↔ «твист»;
6. относительное содержание молекул конформеров в равновесном состоянии равно:
«α-кресло» : «α-ванна» : «α-твист» = 2,62∙107 : 1,70∙107 : 1;
«β-кресло» : «β-ванна» : «β-твист» = 2,55∙1010 : 4,92∙109 : 1;
7. преобладающая часть молекул (60,56% α-формы и 83,81% β-формы) находятся в форме «кресло»;
8. содержание твист-конформеров в равновесной системе для α-формы не превышает 6∙10-6 %, а для β-формы 2∙10-6 %.
Список литературы
1. Corbridge D.E.C. Phosphorus. Chemistry, biochemistry and technology. 6th edition, CRC press, 2013
2. Д. ван Везер Фосфор и его соединения Пер. с англ. под ред. к.т.н. А.И. Шерешевского. — М.: Изд. ин. лит., 1962. — 688 с.
3. Куанышева Г.С., Макашева Г.Р., Джамансариева К.У. Химия и основы модифицирования конденсированных фосфатов. Методическое пособие. — Алматы.: Қазақ университеті, 2002. – 32с.
4. W. P. Griffith Vibrational spectra of metaphosphates, meta-arsenates, and metavanadates J. Chem. Soc. A, 1967, 905-908
5. S. P. Mishra, S.C. Ojha, L.P.Pathak Isomerism, conformation and molecular symmetry in phosphorus compounds Asian J. of Chem. Vol. 20, No. 6 (2008), 4931-4933
6. M.W.Schmidt, K.K.Baldridge, J.A.Boatz, S.T.Elbert, M.S.Gordon, J.H.Jensen, S.Koseki, N.Matsunaga, K.A.Nguyen, S.J.Su, T.L.Windus, M.Dupuis, J.A.Montgomery J. Comput. Chem. 14, 1347-1363 (1993)
7. Кобычев В. Б. Квантовая химия на ПК: Компьютерное моделирование молекулярных систем: учеб.-метод. пособие / В. Б. Кобычев. – Иркутск: Иркут. гос. ун-т, 2006. – 87 с.
8. Блатов В.А., Шевченко А.П., Пересыпкина Е.В. Полуэмпирические расчетные методы квантовой химии: Учебное пособие. Изд. 2-е. Самара: Изд-во «Универс-групп», 2005. — 32 с.
9. Плехович С.Д., Зеленцов С.В. Расчет переходных состояний методами квантовой химии. Учебно-методическое пособие. – Нижний Новгород: НГУ, 2015. – 21 с.
10. Цышевский Р.В., Гарифзянова Г.Г., Храпковский Г.М. Квантово-химические расчеты механизмов химических реакций: учебно-методическое пособие: — Изд-во КНИТУ, Казань, 2012. – 87 с.