Кристаллоструктурный метод расчета ионности связей и эффективных зарядов атомов в минералах

1.2.3. Кристаллоструктурный метод расчета ионности связей и эффективных зарядов атомов в минералах

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

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

Конкретная процедура оценки ионности (эффективных зарядов атомов) минералов производится в ряд этапов (Зуев, 1990):

1. 
   Находится общее число валентных электронов минерала, равное суммарному заряду остовов с обратным знаком. Проблема оценки числа участвующих в химической связи электронов у атомов в соединениях подробно рассмотрена в предыдущем разделе 1.2.2.
   
2. 
   Суммарное число валентных электронов минерала распределяется по соответствующим двухцентровым (межостовным) связям, т. е. оцениваются их электронные заряды.
   
3. 
   Электронные заряды связей делятся на доли электронного облака, локализованные у атомных остовов пропорционально их силовым характеристикам.
   
4. 
   Суммированием зарядов остовов и относящихся к ним долям электронных облаков вычисляются эффективные заряды атомов (ионов).
   

В качестве примеров приведем некоторые ОЭФ минералов и молекул, в которых по сравнению с монографиями (Зуев, 1990; Зуев, Денисов, Мочалов и др., 2000) внесены (в отдельных случаях) некоторые изменения (уточнения) в отношении фактической валентности неметаллов:

Насколько правомерно или, другими словами, отвечает реальности представление кристаллов в виде ассоциации атомных остовов и связующих (а также неподеленных) электронов-электридов, что отражают предлагаемые остовно-электронные формулы? Подробное обоснование этого подхода имеется в уже достаточно обширной литературе (Bent, 1970; Гиллеспи, 1975; Cohen, 1984; Зуев, 1990; Гиллеспи, Харгиттаи, 1992; Семенов, 2001). Поэтому ограничимся лишь некоторыми дополнительными аргументами.

Обычно ковалентную двухэлектронную связь трактуют в виде спаривания электронов у нейтральных атомов: например, Н-Н или Н^¯Н^­ в простейшем случае молекулы водорода Н₂. Однако эти схемы явно упрощают реальную картину ковалентной связи, природу которой раскрывают квантово-химические методы (молекулярных орбит, валентных связей), в рамках которых необходимо учитывать обмен электронами у вступающих в связь атомов. Так, по Полингу (Полинг, 1974, с.140) молекулу водорода можно записать в виде резонанса двух структур:

H_A^↓H_B^↑­H и _A^­↑ H_B^↓

Это значит, что электрон с положительным спином (↓) и электрон с отрицательным спином (­↑) меняются местами, т. е. происходит обмен электронами у атомов водорода. Понятно, что обмен электронами у первоначально нейтральных атомов водорода Н_А и Н_В предполагает отрыв у них электронов (с затратами энергий ионизации Ry = 13,6 эВ), образование остовов [H_A⁺] и [H_B⁺] и локализацию в межостовном пространстве ставшей общей (связующей) пары неразличимых электронов, которые попадают под действие силовых полей обоих остовов. Описанную ситуацию как раз и отражает остовно-электронная модель молекулы водорода - |[H⁺](2e^-)[H⁺]|.[1] Аналогичным образом можно обосновать и другие более сложные остовно-электронные модели молекул и кристаллов.

Имеются специальные формулы для оценки q в кристаллах (Зуев, 1990). Так, в случае простых (бинарного состава МхХу) минералов:

а в случае сложных минералов электронный заряд связи металл-неметалл находится по формуле:

где Z_M, Z_X - числа валентных электронов атома, соответствующих зарядам остовов; ΣЭОl_M, ΣЭОl_X - суммы электроотрицательностей лигандов вокруг атомов М и Х в соответствующих атомных полиэдрах кристаллической структуры. Необходимые для таких расчетов координационные ЭО атомов имеются в (Зуев, 1990).

Как показано ранее в разделе 1.2.1, исходя из сферической формы остовов и связывающих электронных облаков (электридов), обладающих соответствующими размерами и зарядами (рис. 1.3), можно строить остовно-электронные модели структур минералов. На рис. 1.4-1.8 даны примеры таких моделей для рутила, алмаза, кварца, сфалерита и периклаза. Важно подчеркнуть, что эти модели фактически демонстрируют новый остовно-электронный способ изображения кристаллических структур, в какой-то мере адекватный картам разностной электронной плотности (рис. 1.11-1.12), а также расчетам зарядовой электронной плотности в кристаллах по методу псевдопотенциала (Cohen, 1984). В последнем случае на соответствующих схемах, например, для кристалла Si, непосредственно фигурируют как связывающие электронные облака, так и остовы кремния.

Остовно-электронные модели структур можно представлять в двух используемых в кристаллохимии аспектах (Зуев, 1990):

• 
  как постройки из сферических остовов и электридов, обладающих соответствующими зарядами и радиусами (рис. 1.4-1.8);
  
• 
  как ассоциации остовных координационных полиэдров (рис. 1.13-1.16).
  

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

Найденные по формулам (1.25) и (1.26) электронные заряды связей q(M-X) делятся между остовами пропорционально силовым параметрам F_M и F_X, в качестве которых можно использовать величины:

F = Z/(КЧxr_c или 2)F = Z/(Lr_c2), (1.27)

где L- общая лигандность остова, т. е. его координация с учетом неподеленных электронов.

В монографии (Зуев, 1990) приведены многочисленные расчеты Q атомов в минералах с использованием параметров Z и соответственно F. Однако понятно, что одинаковые параметры Z слабо отражают индивидуальные химические особенности таких разных катионных остовов, как [Mg2+] и [Fe2+]. Поэтому более корректно, по-видимому, в формулах (1.27) вместо номинального заряда Z использовать суммарный эффективный заряд ядра, действующий на валентные (сверх остова) электроны атома - ΣZ_i*. Тогда эти формулы пре-образуются в:

F* = (ΣZ_i*)/(КЧxr_c или 2)F* = (ΣZ_i*)/(Lr_c2), (1.28)

где Z_i* - действующий на валентный электрон эффективный заряд ядра, вычисляемый по правилам С. С. Бацанова (Бацанов и Звягина, 1966; Бацанов, 1986), либо рассчитываемый из потенциалов ионизации атома по формуле (1.12). Именно второй способ оценки эффективных зарядов ядер, как более простой, принят в данной работе.

Использованием эффективных зарядов ядер в формулах (1.28) удается в явной форме, хотя и в суммарно усредненном виде, учесть силовые характеристики атома для всех его сверхостовных электронов. Так, в случае муассанита SiC, остовно-электронная формула которого |[Si⁴⁺] (2e^-)₄ [C⁴⁺]|, совершенно очевидна необходимость учета эффективных зарядов ядер для всех четырех валентных электронов атомов Si и C при оценке их силовых параметров.

Вычисленные по формулам (1.27) или (1.28) силовые параметры используются для расчета ионного (полярного) характера связей М-Х кристалла:

λ(М-Х) = (F_X - F_M)/(F_X + F_M) (1.29)

λ(М-Х) = (F_X* - F_M*)/(F_X* + F_M*) (1.30)

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

Схема полярной связи М-Х приведена на рис. 1.17, где чисто ионная связь (λ = 1) является лишь условной абстракцией. На самом деле прообразом такой связи является, например, неподеленная пара электронов при остове кислорода в кварце (рис. 1.6). В этом случае ввиду отсутствия второго (катионного) остова можно принять
F_M = F_M* = 0, тогда по формулам (1.29) и (1.30) λ = 1, т. е. неподеленную пару при остове можно отождествить с химической двухэлектронной связью 100%-полярности. Поскольку в реальных связях М-Х с двумя остовами всегда F_M > 0 и F_M* > 0, то чисто ионные связи в гетероатомных химических соединениях (минералах) невозможны, соответственно нереальны предельно-ионные структуры типа Na⁺Cl^-, Mg²⁺O^2-, Al₂³⁺O₃^2- и др. Впрочем, это не означает, что по распределению зарядов реальные кристаллические соединения не могут быть весьма близки к предельно-ионным моделям типа указанных, о чем свидетельствуют современные данные расчетов электронной плотности.

Заключительная процедура вычисления эффективных зарядов атомов в минералах производится по следующим формулам (Зуев, 1990):

где N - порядок связи М-Х, равный q/2.

Возможен расчет Q атомов в минералах на основе величин параметров q, минуя параметры λ. Обозначим через с - электронную долю связи, локализованную на катионном компоненте, а через q-c соответственно - на анионном компоненте связи М-Х. Тогда из пропорции с/(q-c) = F_M*/F_X* имеем:

 

с = (qF_M*)/(F_M*+F_X *),q-c = (qF_X*)/(F_M*+F_X*). (1.32)

Зная КЧ атомов М и Х в структуре минерала, приходим к следующим формулам:

которые дают те же самые результаты, что и формулы (1.31).

Для сопоставления вычисленных по разработанной методике Q атомов в кристаллических соединениях с соответствующими данными их диэлектрической ионности по Филлипсу-Ван Фехтену рекомендуется выражение:

f_i = Q/Q_и, (1.34)

представляющее собой степень ионизации атомов, равную отношению эффективных зарядов к предельно-ионным зарядам атомов. В монографии (Зуев, 1990) даны примеры хорошего соответствия оценок эффективных зарядов атомов в кристаллах в рамках двух указанных подходов.

Разработанная методика оценки ионности (эффективных зарядов атомов) универсальна в том смысле, что применима к соединениям различного типа, будь то молекулы, кристаллы, комплексы и т. д. Соответственно она была успешно апробирована на большом числе минералов и искусственных соединений. Для расчетов были использованы эффективные заряды ядер для валентных электронов атомов согласно (Zhang, 1982), электронные заряды связей в соединениях рассчитывались по приведенным формулам, либо непосредственно заимствовались из книги (Зуев, 1990). Ковалентные радиусы неметаллических компонентов r_c(X) были приняты с учетом данных Слэтера, Урусова, Лебедева (Лебедев, 1969; Урусов, 1975; Зуев, 1990), а металлических - находились по разности r_c(M) = d(M-X) - r_c(X), где d(M-X) - среднее межатомное расстояние катион-анион в структуре.

В качестве примеров рассчитаем эффективные заряды атомов в молекуле NaCl, галите NaCl, cфалерите ZnS, периклазе MgO, кварце SiO₂, форстерите Mg₂SiO₄ и бишофите MgCl₂(H₂O)₆.

Остовно-электронная формула молекулы NaCl - |[Na⁺]2e^-[(Cl⁷⁺)6e^-]|, она обладает ординарной полярной связью Na-Cl. Размерные характеристики молекулы: d(Na-Cl) = 2,36 Å, r_c(Cl) = 0,7 Å, r_c(Na) = 1,66 Å. Эффективные заряды ядер для остовов ΣZ_i*(Na⁺) = 1,78, ΣZ_i*(Cl⁺) = 2,82, по формулам (1.28) F*(Na⁺) = 1,78/1,662 = 0,646, F*(Cl+) = 2,82/0,72 = 5,755, по формуле (1.30) l(Na-Cl) = (5,755 - 0,646)/(5,755 + 0,646) = 0,8, т. е. получаем Na+0.8Cl-0.8, что соответствует ионности, вычисленной из дипольного момента молекулы NaCl (Бацанов, 1986).

Остовно-электронная формула галита была приведена выше, q(Na-Cl) = 1e-, d(Na-Cl) = 2,82 Å, r_c(Cl) = 1,0 Å, r_c(Na) = 1,82 Å, в этом случае F*(Na⁺) = 1,78/(6x1,822) = 0,09, F*(Cl⁵⁺) = 23,76/(6x12) = 3,96. По формуле (1.30) λ(Na-Cl) = 0,96 и далее по формулам (1.31): Q(Na) = 1 - 6x0,5(1 - 0,96) = +0,88, Q(Cl) = 5 - 6x0,5(1 + 0,96) = -0,88, т. е. для галита имеем Na^+0.88Cl^-0.88, что соответствует диэлектрической ионности кристалла NaCl (Бацанов, 1982). Обратим внимание на заметный рост эффективных зарядов атомов при переходе от молекулы к кристаллу.

В сфалерите ZnS КЧ атомов равны 4, q(Zn-S) = 2e^-, d(Zn-S) = 2,35 Å, F*(Zn²⁺) = 6,83/(4x1,35²) = 0,94, F*(S⁶⁺) = 30,35/(4x1²) = 7,59, по формуле (1.30) λ(Zn-S) = 0,78 и далее по формулам (1.31), в которых N = 1: Q(Zn) = 2 - 4(1 - 0,78) = +1,12 и Q(S) = 6 - 4(1 + 0,78) = -1,12, т. е. для сфалерита получаем Zn^+1.12S^-1.12.

В периклазе MgO (рис. 1.8) КЧ обоих атомов равны 6, q(Mg-O) = 1e^-, d(Mg-O) = 2,10 Å, r_c(O) = 0,55 Å, r_c(Mg) = 1,55 Å, SZ_i*(Mg²⁺) =

В случае кварца SiO₂ (рис. 1.6) КЧ(Si) = 4(O), КЧ(О) = 2(Si), q(Si-O) = 3e-, d(Si-O) = 1,61 Å, r_c(O) = 0,55 Å, r_c(Si) = 1,06 Å, SZ_i*(Si⁴⁺) = 15,205, SZ_i*(O⁴⁺) = 13,94, F*(Si4+) = 15,205/(4x1,06²) = 3,383, F*(O4+) = 13,94/(2x0,55²) = 23,04. По формуле (1.30) λ(Si-O) = 0,744 и далее по формулам (1.31) находим: Q(Si) = 4 - 4x1,5(1 - - 0,744) = +2,46, Q(O) = 4 - 2x1,5(1 + 0,744) = -1,23. Таким образом, для кварца имеем Si^+2.46O₂^-1.23, что соответствует его диэлектрической ионности f_i = 0,6 (Бацанов, 1982).

В форстерите Mg₂SiO₄, структура которого характеризуется КЧ(Mg) = 6(O), КЧ(Si) = 4(O) и КЧ(О) = 3Mg+Si, усредненные электронные заряды связей, межатомные расстояния и ковалентные радиусы атомов: q(Mg-O) = 1,23e^-, q(Si-O) = 2,31e^-, d(Mg-O) = 2,12 Å, d(Si-O) = 1,63 Å, r_c(Mg) = 1,57 Å, r_c(Si) = 1,08 Å, r_c(O) = 0,55 Å. Суммарные эффективные заряды ядер для валентных электронов атомов - ΣZ_i*(Mg²⁺) = 5,21, ΣZ_i*(Si⁴⁺) = 15,205, ΣZ_i*(O⁴⁺) = 13,94, силовые параметры остовов - F*(Mg²⁺) = 5,21/(6x1,57²) = 0,352, F*(Si⁴⁺) = 15,205/(4x1,08²) = 3,259, F*(O⁴⁺) = 13,94/(4x0,55²) = 11,521. По формуле (1.30) λ(Mg-O) = 0,941, λ(Si-O) = 0,559 и далее по формулам (1.31): Q(Mg) = 2 - 6x0,615(1 - 0,941) = +1,782, Q(Si) = 4 - 4x1,155(1 - 0,559) = +1,963, Q(O) = 4 - 3x0,615(1 + 0,941)-1x1,155(1 + 0,559) = -1,382, т. е. приходим к Mg₂^+1.78Si^+1.96O₄^-1.38. Найденные в форстерите эффективные заряды оказались удивительно близкими к полученным японскими учеными из прецизионных расчетов электронной плотности - Mg₂^+1.75Si^+2.1O₄^<-1.4> (Fujino, Sasaki, Takeuchi, Sadanada, 1981).

В качестве заключительных примеров расчета Q атомов рассмотрим комплексные кристаллические соединения.

Бишофит, как показали структурные исследования, имеет молекулярную структуру, состоящую из октаэдрических комплексов [Mg(H₂O)₆]²⁺, нейтрализованных двумя ионами Cl^-. Соответственно формула минерала - [Mg(H₂O)₆]²⁺Cl₂^-, а структурная формула отдельного иона [Mg(H₂O)₆]²⁺:

Здесь атом кислорода, имея плоско-треугольную координацию - КЧ(О) = Mg+2H, находится, по-видимому, в состоянии 3s2p²-гибри­дизации с неподеленной парой 2s²-электронов, и соответственно с остовом [O⁴⁺]. Межатомные расстояния: d(Mg-O) = 2,08 Å, d(H-O) = 1,0 Å, ковалентные радиусы: r_c(Mg) = 1,53 Å, r_c(H) = 0,43 Å, r_c(O) = 0,55 Å. Исходное распределение заряда в комплексе - [Mg2+(H₂O)₆], т. е. число валентных электронов у атома магния равно нулю. Тогда расчет по формуле (1.26) дает q(Mg-O) = -[(0x3,5)/(6x3,5) + (4x1,3)/(1,3 + 2x2,1)] = 0,945e^-, q(H-O) = -[(1x3,5)/(1x3,5) + (4x2,1)/(1,3 + 2x2,1)] = 2,527e^-. Силовые параметры остовов: F*(Mg²⁺) = 5,21/(6x1,53²) = 0,37, F*(H⁺) = 0,85/(1x0,43²) = 4,6, F*(O⁴⁺) = 13,94/(3x0,55²) = 15,36. Воспользовавшись в данном случае формулами (1.32), для связи Mg-O находим с = 0,0224e^-, q-c = 0,9231e^-, а для связи H-O c = 0,5824e^-, q-c = 1,9446e^-. Далее по формулам (1.33) находим Q(Mg) = 2 - 6x0,0224 = +1,866, Q(H) = 1 - 1x0,5824 = +0,417, Q(O) = 4 - 1x0,9231 - 2x1,9446 = -0,812. Результирующее распределение зарядов атомов в бишофите оказалось следующим - [Mg^+1.866x(H₂^+0.417O^-0.812)₆]²⁺Cl₂^-. Баланс зарядов в молекуле воды приводит к (H₂O)^+0.022, что позволяет понять природу примеси ковалентной составляющей связи в данном и других аква-комплексах.[2] Как это очевидно, дело заключается в незначительной перекачке электронной плотности от молекул воды на свободные орбитали центрального катиона. Результатом этого является некоторое уменьшение формального заряда центрального катиона и соответственно - появление на молекулах воды указанного избыточного положительного заряда. Другими словами, от каждой из шести первоначально электронейтральных молекул воды на свободные орбитали центрального катиона (Mg²⁺) переходит 0,022e-, и он приобретает эффективный заряд +2 - 0,022x6 = +1,868.

В случае бианкита ZnSO₄x6H₂O конституцию этого минерала можно рассматривать как соединение из октаэдрических и тетраэдрических комплексных ионов - [Zn(H₂O)₆]²⁺ и [SO₄]^2-. Расчет приводит к следующим эффективным зарядам атомов в октаэдрическом комплексе [Zn^+1.83(H₂^+0.418O^-0.807)₆]²⁺. Эффективные заряды атомов в комплексе [SO₄]^2- вычислены для структуры с двойными связями и следующими параметрами: d(S-O) = 1,50 Å, r_c(O) = 0,5 Å, r_c(S) = 1,0 Å, q(S-O) = 4,0e^-, F*(S⁶⁺) = 30,35/(4x1²) = 7,59, F*(O²⁺) = 5,19/(1x0,52) = 20,76. Использование этих параметров приводит к следующим эффективным зарядам в сульфатном комплексе– [S^+1.72O₄^-0.93]^2-, а общее распределение зарядов в бианките будет таковым: [Zn^+1.83(H₂^+0.418xO^-0.807)₆]²⁺[S^+1.72O₄^-0.93]^2-. Для водного сульфата FeSO₄×6H₂O аналогичные расчеты приводят к следующему распределению эффективных зарядов - [Fe^+1.86(H₂^+0.411O^-0.80)₆]²⁺[S^+1.72O₄^-0.93]^2-.

Для октаэдрического аква-комплекса Fe(III) расчеты дали такие эффективные заряды атомов - [Fe^+2.69(H₂^+0.424O^-0.796)₆]³⁺.

Для безводных сульфатов M(II)SO₄, структура которых характеризуется следующей взаимной координацией атомов: КЧ(М) = 2О_А + 4О_В, КЧ(S) = 2O_A + 2O_B, КЧ(О_А) = M + S, КЧ(О_В) = 2М + S, получились в порядке уменьшения Q(M) такие данные:

Mn^+1.843[S^+1.418O_2(A)^-0.662O_2(B)^-1.969]

Обратим внимание на антибатный ход изменения эффективных зарядов нерадикальных катионов M(II) и радикальных катионов S(VI) в рассмотренном ряду сульфатов и некоторое уменьшение Q(S) при переходе от свободного радикала [SO₄]^-2 к связанному в кристаллической решетке безводному сульфату.

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

Таблица 1.32

Эффективные заряды атомов в некоторых комплексах и молекулах, рассчитанные в предположении двойных связей «катион-кислород»

Комплекс (молекула)	d(M-On), Å	Эффективные заряды атомов	Комплекс	d(M-O4), Å	Эффективные заряды атомов
[NO3]-	1,24	N+0,58O3-0.527	[ClO4]-	1,51	Cl+1.81O4-0.703
[CO3 ]2-	1,28	C+1O3-1	[MnO4]-	1,59	Mn+2.22O4-0.805
[SO4 ]2-	1,50	S+1.72O4-0.93	[CrO4 ]2-	1,65	Cr+2.53O4-1.133
ArO4	1,50	Ar+1.94O4-0.485	[WO4 ]2-	1,80	W+2.8O4-1.2
FeO4	1,50	Fe+2O4-0.5	[MoO4 ]2-	1,83	Mo+3.06O4-1.26

В рамках продемонстрированного подхода были выполнены расчеты Qатомов в тех же минералах, что и в монографии (Зуев, 1990, стр. 55, таблица 3.9), в результате выяснилась близость расчетных эффективных зарядов атомов в обоих случаях и подтвердились три общих закономерности распределения этих зарядов в сложных минералах: 1) явление (или эффект) взаимного влияния катионов (ВВК); 2) различие в эффективных зарядах кристаллохимически нетождественных (по числу или сорту координирующих катионов) анионов, например, кислорода, обозначаемых О_А, О_В, О_Си т. д.; 3) выполнения известного второго правила Л. Полинга о локальной компенсации ионных зарядов в структуре, если в качестве таковых использовать рассчитанные величины Q атомов.

Эффект ВВК проиллюстрируем следующими уравнениями с указанием эффективных зарядов атомов:

Mg^+1.73O^-1.73+Si^+2.46O₂^-1.23 = Mg^+1.79Si^+2.20O₃^<-1.33>;

2Mg^+1.73O^-1.73+Si+2.46O₂^-1.23 = Mg₂^+1.78Si^+1.96O₄^-1.38;

Al₂^+2.38O₃^-1.587+Si^+2.46O₂^-1.23 = Al₂^+2.44Si^+2.14O₅^<-1.404>;

Zr^+3.21O₂^<-1.61>+Si^+2.46O₂^-1.23 = Zr^+3.36Si^+2.14O₄^-1.375,

т. е. в сложном минерале (энстатите, форстерите, кианите, цирконе) по сравнению с простыми составляющими оксидами (периклазом, кварцем, корундом, бадделеитом,) происходит увеличение Q (уменьшение ковалентности связи) более слабого катиона (Mg, Al, Zr) за счет уменьшения Q (роста ковалентности связи) более сильного катиона (Si).

Подробно эффекты ВВК в сложных минералах рассмотрены в монографии (Зуев, 1990), однако по поводу эффектов ВВК необходимо сделать весьма существенное дополнение, упущенное в указанной монографии. В самом деле, как меняется ионность (ковалентность) связей при переходе от исходных простых составляющих к сложному минералу? Для выяснения данного вопроса сравним средние параметры ионности f_i и ковалентности (f_c = 1-f_i) для простых составляющих кристаллов, с одной стороны, и для соответствующего сложного соединения, с другой стороны. Такие сравнительные оценки удобно произвести по анионным компонентам соединений с использованием формулы 1.34, примеры соответствующих оценок приведены в таблице 1.33.

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

Расчеты показали, что аналогичная картина имеет место и в случае содержащих более двух сортов катионов сложных кристаллических соединений, например: Ca₃Al₂Si₃O₁₂ (Δf_i/Df_c = -0,07/+0,07), Mg₃Al₂Si₃O₁₂ (Δf_i/Df_c = -0,04/+0,04), CaTiSiO₅ (Δf_i/Δf_c = -0,05/+0,05), Be₃Al₂Si₆O₁₈ (Δf_i/Δf_c = -0,02/+0,02) и т. д.

В соединениях с разносортными катионами и анионами, например, в топазе, установленная закономерность сохраняется:

2/3Al₂^+2.38O₃^-1.587 + 2/3Al^+2.592F₃^-0.864 + Si^+2.46O₂^-1.23 = Al₂^+2.50Si^+2.193O₄^-1.44F₂^-0.72 (Δf_i/Δf_c = -0,04/+0,04).

В свете установленной закономерности полученные в рамках компьютерного моделирования структур эффективные заряды (Урусов, Оганов, Еремин, 1998) в кристаллах: Al₂^+1.86O₃^-1.24 + Si^+1.92O₂^-0.96 = Al₂^+1.83Si^+2.08O₅^-1.15 с противоположным (по сравнению со всеми предыдущими соответствующими данными таблицы 1.33) изменением величин в параметрах Δf_i/Δf_c = +0,01/-0,01 не могут считаться корректными, поскольку находятся в противоречии с эффектом ВВК. Сказанное подтверждают данные изучения этих кристаллов методом ЭСХА (РФЭС) (таблица 1.34).

Вопрос о том, как меняется ионность и ковалентность при переходе от простых составляющих к сложному кристаллическому соединению, является принципиальным и очень важным, объясняя причину распространенности среди минералов сложных кристаллических соединений. Так, установленный при этом рост ковалентности (снижение ионности) ведет, по нашему мнению, к энергетической стабилизации сложного соединения, термодинамическим выражением которой может, очевидно, служить превышение его энтальпии над суммой энтальпий простых составляющих соединений (Зуев, 1990). Другими словами, отрицательную величину энтальпии образования сложного кристаллического соединения (например, оксида) из простых составляющих оксидов вполне логично и обоснованно можно связать с соответствующим ростом ковалентности первого по сравнению со средней ковалентностью вторых. В частности, получают объяснение: почти двукратное увеличение энтальпии образования из оксидов форстерита Mg₂SiO₄ (Δf_c = +0,05, ΔH^o_ox = -58 кДж/моль) по сравнению с энстатитом MgSiO₃ (Δf_c = +0,03, ΔH^o_ox = -36 кДж/моль), низкая величина энтальпии образования из оксидов кианита Al₂SiO₅ (Δf_c = +0,02, ΔH^o_ox = -8 кДж/моль) и, наоборот, весьма высокая величина энтальпии образования из оксидов шеелита CaWO₄ (Δf_c = +0,18, ΔH^o_ox = -162 кДж/моль) и т. д. Попытка построения соответствующей корреляции (рис. 1.18) привела к следующей эмпирической формуле:

-ΔH^o_ox(кДж/моль) = 40,425Ln(Δf_c) + 212. (1.35)

Следует признать, что эта формула дает очень грубые оценки энтальпии образования сложного оксида из простых составляющих оксидов, гораздо более корректные результаты дают изложенные в других наших работах (Zuev, 1987; Зуев, 1990) соответствующие методики.

В заключение раздела на примере форстерита Mg₂SiO₄ покажем, как влияет принятие различных параметров Z, Z* и r_c на результаты расчетов эффективных зарядов атомов в рамках кристаллоструктурного метода. С этой целью ниже приводятся следующие две серии данных:

Параметры: Z(Mg) = +2 Z(Si) = +4 Z(O) = +4	Параметры: ΣZi*(Mg2+) = +5,206 ΣZi*(Si4+) = +15,205 ΣZi*(O4+) = +13,94
Mg2+1.64Si+1.66O4-1.235 Mg2+1.71Si+2.10O4-1.355 Mg2+1.77Si+2.487O4-1.507	Mg2+1.72Si+1.50O4-1.235 Mg2+1.78Si+1.96O4-1.38 Mg2+1.83Si+2.37O4-1.507

В обоих случаях первый вариант эффективных зарядов (верхние строки) отвечают расчетам при r_c(O) = 0,6 Å, второй вариант (средние строки) - r_c(O) = 0,55 Å, третий вариант (нижние строки) - r_c(O) = 0,5 Å. Соответствующие ковалентные радиусы катионных компонентов (Mg, Si) находились по разности экспериментальных межатомных расстояний и фиксированных величин r_c(O).

Из приведенных данных следует вывод о заметном влиянии на результаты расчетов эффективных зарядов атомов лишь параметров ковалентных радиусов, что было установлено ранее (Зуев, Денисов, Мочалов и др., 2000). Оптимальные заряды получаются, как нам представляется, по второму варианту расчета.

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

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

[1] В использованном Полингом способе изображения - H_A(↓↑­)H_B или H_A(­↑↓)H_B.

[2] Расчет эффективных зарядов в приближении остовно-электронного подхода для свободной молекулы воды дает H₂^+0.34O^-0.68 в хорошем согласии с расчетами из дипольного момента и электроотрицательностей.