САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ
На правах рукописи
УДК 541.5
Специальность 02.00.01 - неорганическая химия
А В Т О Р Е Ф Е Р А Т
диссертации на соискание ученой степени
кандидата химических наук
Санкт-Петербург
1995
Актуальность темы. По мере развития вычислительных методов квантовой химии молекул и кристаллов становится все труднее интерпретировать результаты расчета в терминах локальных свойств электронной структуры, таких как атомные заряды, электронные конфигурации, кратности связей, валентности. В большинстве теоретических работ, посвященных исследованию электронной структуры молекул и, особенно, кристаллов, основное внимание уделяется описанию электронных и энергетических свойств вещества в целом, тогда как природа химической связи и локальные свойства, связанные с состоянием атомов и связей в веществе не изучаются.
Классическое определение валентности атома ориентировано на описание достаточно простых молекул, находящихся в нормальных условиях. Развитие подходов к определению валентности атомов в твердых телах или в возбужденных молекулах неочевидно. Квантовохимическая интерпретация валентности, как величины, характеризующей общее число электронов атома, участвующих в образовании химических связей, позволяет не только реабилитировать это понятие классической химии, но и расширить область его практического применения. В частности, в последнее время делаются попытки объяснить активность молекул в химических реакциях, необычные физико-химические свойства идеальных и дефект-содержащих кристаллов наличием в них атомов, находящихся в необычном валентном состоянии.
Интересными объектами исследования, для которых характерен широкий спектр изменения химической связи (от ионной к ионно-ковалентной и металлической), являются смешанные кристаллические оксиды sp- и sd- элементов. Многие из этих соединений проявляют необычные свойства (высокотемпературная сверхпроводимость, фотопроводимость, активность в каталитических и хемосорбционных процессах на поверхности), что обуславливает постоянный интерес к ним как к технологически перспективным объектам.
Целью настоящей работы является развитие аппарата, применяемого для анализа локальных свойств электронной структуры, построение квантовохимического определения валентности атома в молекулах и кристаллах и изучение природы химической связи в сложных кристаллических оксидах sp и sd элементов: свинца, висмута, лантана, никеля и меди.
Научная новизна. В настоящей работе предлагается новая схема расчета валентности атома, основанная на квантовохимическом анализе распределения электронной плотности в молекуле и кристалле и изучении составных частей химической связи - ионной, ковалентной и металлической. Впервые выполнены расчеты электронной структуры ряда кристаллических оксидов свинца, висмута, никеля и меди (всего более сорока соединений). Исследована природа химической связи в изученных соединениях и выявлены закономерности в изменениях природы химической связи в оксидах sp и sd элементов при изменении формальной степени окисления атомов металла.
Практическая ценность. Сделанное в работе обобщение разработанных в квантовой химии молекул подходов к анализу локальных свойств электронной структуры на случай кристаллических твердых тел, открывает новые возможности для изучения устойчивости поверхности и точечных дефектов в твердых телах. Выполненные расчеты электронной структуры кристаллических оксидов имеют большое значение для интерпретации экспериментальных данных. Особенный интерес представляет исследование электронной структуры и анализ природы химической связи в кристаллах высокотемпературных сверхпроводников и Мотт-Хаббардовских диэлектриков. Предложенная формула определения полной валентности атома реализована в ряде квантовохимических программ расчетов молекул, кластеров и кристаллов.
Публикации. По материалам диссертации опубликованы 14 статей, один обзор и 9 тезисов докладов.
Апробация работы. Материалы диссертации докладывались на II и III Всесоюзном совещании по ВТСП (Харьков 1989, 1991), XIV всесоюзном (Пекаровском) совещании по теории полупроводников (Донецк, 1989), Всесоюзном совещании по квантовой химии (Новосибирск, 1990), I международной конференции "Химия твердого тела" (Одесса 1990), конференции по квантовой химии твердого тела (Рига 1990), Международной конференции по дефектам в диэлектриках ICDIM-92 (Нордкирхен, ФРГ, 1992), II международной конференции по применению компьютеров в материаловедении CAMSE'92 (Йокогама, Япония, 1992), 5-международной конференции по электронной спектроскопии ICES-5 (Киев, Украина, 1993), I Международной конференции по теоретической химической физике (Жирона, Испания, 1993), IV международной конференции механизмы и материалы ВТСП M2S-HTSC-IV (Гренобль, Франция, 1994),
Положения, выносимые на защиту:
Объем и структура диссертации. Работа изложена на 190 страницах машинописного текста, включает 14 рисунков, 18 таблиц и приложение. Основной текст состоит из введения, обзора литературы, оригинальной части, выводов и списка литературы, включающего 129 наименований.
Обзор литературы включает в себя обсуждение следующих вопросов: Как с позиций классической химии определяется валентность атомов и описывается химическая связь в соединениях? Какие подходы применяются в квантовой химии для вычисления распределения электронной плотности в молекулах и кристаллах? Какие характеристики электронной структуры вещества можно сопоставить с атомными зарядами, порядками связей и ковалентностями атомов?
Хотя классическое определение валентности атома, основанное на анализе стехиометрической формулы соединения, далеко не всегда приводит к однозначным результатам, можно выделить достаточно большое количество молекул, для которых применение такого подхода оправдано. Вводя то или иное определение валентности атома, следует анализировать его корректность, сравнивая результаты, вычисленные в рамках нового определения со стехиометрическими значениями для ряда указанных молекул.
При использовании расчетных схем, основанных на приближении Хартри-Фока-Рутаана (метод ЛКАО) удается вычислить матричные элементы матрицы плотности молекулы в ортогональном атомном базисе. Для кристалла, описываемого в рамках зонной модели, находят матрицу плотности, построенную в блоховском базисе. Однако она может быть преобразована к матрице плотности в локализованном атомном базисе.
Квантовохимический подход к анализу природы химической связи в соединении состоит в вычислении атомных зарядов QA (описывающих ионную компоненту химической связи) и индексов Вайберга WAB (описывающих ковалентную составляющую связи) (Lцwdin P.O. Adv. Phys. 5 (1956) 1; Wiberg K.B. Tetrahedron. 24 (1968) 1083).
Эти величины определяются матричными элементами одночастичной матрицы плотности Pab построенной в ортогональном атомном базисе.
(1)
(2)
орбитали a,b центрированы на атомах A и B, соответственно.
Сумма индексов Вайберга по всем атомам, кроме данного, определяет общее число электронов, участвующих в образовании ковалентных связей (Борисова Н.П., Семенов С.Г. Вестн. ЛГУ. N16, (1973) 119; Armstrong D.R., Perkins P.G., Stewart J.J.P. J. Chem. Soc. Dalton Trans. 7 (1973) 838)
(3)
Используя свойство идемпотентности матрицы плотности, выражение (3) можно переписать так, чтобы величина CA определялась только элементами матрицы плотности, относящимися к диагональному атомному блоку A-A:
(4)
В квантовой химии молекул для определения валентности атома используют выражение (3), тогда как в физике твердого тела часто отождествляют QA и валентность.
В настоящей работе показано, что оба этих подхода корректны только в крайних случаях соответственно, ковалентной неполярной и чисто ионной связи. На рисунке 1 представлена зависимость отношения CA к стехиометрической валентности от степени ионности (величины относительного атомного заряда) для ряда молекул. Приведенные результаты взяты из литературных источников или получены в результате расчета молекул (в экспериментальной геометрии) методами нулевого дифференциального перекрывания (CNDO, AM1).
Рис. 1
Только для соединений с малой степенью ионности численные значения CA близки к стехиометрическим валентностям, тогда как для ионных соединений величина CA стремится к нулю. Фактически, CA описывает только ковалентную составляющую химической связи, и ее следует называть ковалентностью атома.
Очевидно, что использование только атомных зарядов или только ковалентностей для описания валентного состояния атомов ограничено крайними случаями соответственно чисто ионной связи и ковалентной неполярной связи. Естественно предположить, что валентность атома VA, как характеристика общего числа электронов данного атома вступивших в образование ионных и ковалентных связей, должна зависеть от величин QA и CA (VA=f(QA,CA)). Для определения этой функциональной зависимости, потребуем чтобы 1) в указанных крайних случаях величина VA совпадала соответственно с QA или CA. 2) величина VA имела единичное значение для атома, валентная оболочка которого содержит один электрон на единственной s орбитали (таким атомом, например, является водород или атом щелочного металла).
Этим условиям удовлетворяет следующее выражение:
(5)
Выражение (5) совпадает с ковалентностью (3) для чисто ковалентных соединений, совпадает с абсолютной величиной атомного заряда (1) для чисто ионных систем и тождественно равно единице для водородо-подобного атома, для которого как следует из (1) и (4) CA=1-QA2. Величина VA введенная в выражении (5) описывает как ионный, так и ковалентный вклады в химическую связь, поэтому ее можно назвать полной валентностью атома.
Для соединений с существенно ионным характером
химической связи
выражение (5)
переходит в
, которое имеет
наглядную интерпретацию: полная валентность атома
складывается из абсолютной величины атомного заряда,
перенесенного в ионной связи, и половины всех электронных
пар, участвующих в образовании ковалентных связей. В общем
случае способ "деления" ковалентных электронных пар носит
более сложный характер, который определяется балансом между
"связывающими" и неподеленными электронами атома.
Предположив, что зависимость VA от QA и CA имеет не более чем квадратичный характер и описывается квадратичной формой с девятью произвольными коэффициентами
, можно строго показать, что выражение
(5) является единственным, удовлетворяющим следующим
условиям:
Корректность введенного определения полной валентности (5) продемонстрирована на ряде модельных систем, для которых изучены границы допустимых значений VA.
Подставляя значения QA и CA , представленные на рисунке 1, в выражение (5) полной валентности атома, находим, что отличия значений полной валентности атома от стехиометрических - незначительны (см. рис. 2).
Рис. 2
Вычисленные по (5) значения квантовохимических валентностей атомов с постоянной (Li(I), F(I), Mg(II), O(II), B(III)) и переменной (C(II)-(IV), Si(II)-(IV), N(III)-(V), P(III)-(V), S(II)-(IV)-(VI)) валентностью в молекулах с различным типом химической связи находятся в хорошем согласии с классическими представлениями.
В работе изучены возможные причины отклонений найденных валентностей атомов от стехиометрических значений: 1) использованная в расчете структура молекулы или кристалла не соответствует равновесной геометрии; 2) приближения, использованные в квантовохимическом расчете, не позволяют удовлетворительно описать электронную структуру системы; 3) природа химической связи в данном веществе носит особый характер, что может проявляться в активности данного соединения в химических реакциях.
Для обобщения на кристаллы разработанных для молекул подходов к анализу локальных свойств электронной структуры прежде всего следует определить элементы матрицы плотности в локализованном атомном базисе. Особая проблема связана с нахождением ковалентности атома, для вычисления которой необходимо проводить суммирование в (3) по всему бесконечному кристаллу.
При изучении металлических соединений возникает проблема описания делокализованной металлической связи в терминах локальных характеристик (заряды, порядки связей). В работе изучено поведение индексов Вайберга при увеличении длины связи в диэлектрических и металлических соединениях. Показано, что обычная ковалентная связь образуется парами катион-анион, ее порядок быстро спадает при увеличении длины связи. В металлах наблюдается немонотонный характер убывания индексов Вайберга при увеличении расстояния.
Для выделения металлической компоненты в соединениях со
смешанным характером связи (ионно- ковалентно-
металлической), предлагается разделить суммарную
ковалентность на гомоатомную, определяемую связями между
атомами одного элемента, и гетероатомную.
(6)
Заключительная глава диссертации посвящена изучению электронного строения и химической связи в кристаллических оксидах. В качестве объектов исследования представлены кристаллические оксиды s-(щелочные, щелочноземельные металлы), sp-(свинец, висмут), и sd-(лантан, иттрий, медь, никель) металлов.
В рамках зонной модели, в приближении CNDO (полное пренебрежение дифференциальным перекрыванием) изучены следующие кристаллические оксиды: I) PbO-a, PbO-b, PbO2-a, PbO2-b, Pb2O3, Pb3O4, BaPbO3, II) Bi2O3-a, Bi2O3-b, Bi2O3-d, KBiO3, BaBiO3, III) La2O3, IV) NiO, La2NiO4-tetr, La2NiO4-orth, Li2NiO2, Na2NiO2, LiNiO2, NaNiO2, LaNiO3, LaNiO2, LaSrNiO3, BaNiO3, V) Cu2O, CuO, Cu4O3, KCuO, YCuO2, LaCuO2, La2CuO4-tetr, La2CuO4-orth, Li2CuO2, KCuO2, NaCuO2, LaCuO3, YBa2Cu3O6, YBa2Cu3O7. Самосогласованная матрица плотности вычислялась интегрированием по 64 точкам зоны Бриллюэна. Для описания соединений d- элементов использовался расширенный по спину метод расчета. Атомные параметры полуэмпирической схемы CNDO калибровались по экспериментальным данным для простых оксидов (PbO, Bi2O3, NiO, Cu2O, La2O3).
Получен обширный материал (полные и парциальные плотности состояний, дисперсионные кривые, энергии оптических переходов), позволяющий установить закономерности в электронном строении оксидов свинца, висмута, никеля и меди. Показано, что использование схемы CNDO позволяет передать основные электронные характеристики соединений sp и sd элементов, находящихся в различных степенях окисления. Изучена природа диэлектрической щели в оксидах меди и никеля, и показано, что эти соединения относятся к диэлектрикам с переносом заряда.
В таблице приведены значения атомных зарядов QA, заселенностей орбиталей Ns (Np, Nd), наибольшего индекса Вайберга связи Me-O, гомоатомной и гетероатомной составляющей ковалентности и полной валентности атомов в некоторых из изученных соединений.
| Атом | Кристалл | Q | Ns | Np/Nd | WMeO | Chetero | Chomo | V |
| Pb(II) | PbO-a | 0.76 | 1.84 | 1.40 | 0.43 | 1.87 | 0.15 | 2.14 |
| Pb(IV) | PbO2-b | 2.06 | 0.75 | 1.19 | 0.48 | 2.78 | 0.07 | 3.87 |
| Bi(III) | Bi2O3-d | 2.22 | 1.99 | 0.79 | 0.21 | 1.30 | 0.08 | 2.97 |
| Bi(V) | KBiO3 | 4.19 | 0.54 | 0.28 | 0.20 | 1.30 | 0.01 | 4.89 |
| Y(III) | YCuO2 | 2.47 | 0.08 | 0.45 | 0.16 | 0.98 | 0.03 | 3.01 |
| La(III) | La2O3 | 2.73 | 0.03 | 0.25 | 0.10 | 0.52 | 0.01 | 3.00 |
| Ni(II) | NiO | 1.77 | 0.07 | 8.16 | 0.07 | 0.42 | 0.02 | 1.99 |
| La2NiO4-t | 1.69 | 0.07 | 8.24 | 0.12 | 0.54 | 0.02 | 1.99 | |
| Ni(I) | LaNiO2 | 0.95 | 0.03 | 9.02 | 0.03 | 0.20 | 0.82 | 1.05 |
| Ni(III) | LaNiO3 | 2.39 | 0.05 | 7.56 | 0.15 | 1.00 | 0.96 | 2.94 |
| Атом | Кристалл | Q | Ns | Nd | WMeO | Chetero | Chomo | V |
| Cu(I) | Cu2O | 0.96 | 0.04 | 9.99 | 0.02 | 0.07 | 0.02 | 1.00 |
| Cu(II) | CuO | 1.66 | 0.07 | 9.26 | 0.13 | 0.50 | 0.04 | 1.93 |
| La2CuO4-t | 1.68 | 0.06 | 9.26 | 0.11 | 0.50 | 0.03 | 1.94 | |
| Cu(III) | KCuO2 | 1.69 | 0.06 | 9.25 | 0.26 | 1.05 | 0.00 | 2.30 |
| LaCuO3 | 1.98 | 0.05 | 8.98 | 0.13 | 1.09 | 0.48 | 2.59 |
При увеличении формальной степени окисления атомов sp и sd элементов степень ионности понижается, а ковалентность растет. Наблюдается рост ионности в ряду оксидов sp - sd - s элементов. Найдено, что простые оксиды обладают наибольшей степенью ионности по сравнению с тетрарными оксидами.
При вычислении полной валентности атома в металлических соединениях необходимо учитывать только гетероатомную часть ковалентности. В таблице приведены значения суммарной ковалентности, ее гомоатомной части (связанной с порядками связей металл-металл, и описывающей металлическую компоненту химической связи) и значения полной валентности, при вычислении которой использовалась либо вся ковалентность CA (Vtotal), либо лишь ее гетероатомная составляющая (Vhetero). (В скобках для сравнения указаны диэлектрические соединения)
| Кристалл | Ctotal | M=Chomo | Vtotal | Vhetero | |
| Ni(III) | LaNiO3 | 1.96 | 0.96 | 3.59 | 2.94 |
| LiNiO2 | 1.92 | 1.15 | 3.70 | 2.98 | |
| Ni(I) | LaNiO2 | 1.02 | 0.82 | 1.62 | 1.05 |
| Ni(I) | (LaSrNiO3) | 0.16 | 0.02 | 1.05 | 1.05 |
| Cu(II) | YBa2Cu3O7^ | 0.91 | 0.33 | 2.10 | 1.88 |
| Cu(II) | (CuO) | 0.54 | 0.04 | 1.96 | 1.93 |
| Cu(III) | LaCuO3 | 1.56 | 0.48 | 2.94 | 2.59 |
| Cu(III) | (KCuO2) | 1.06 | 0.00 | 2.30 | 2.30 |
^ - атом меди в [CuO2] плоскости.
Как следует из приведенной таблицы, для диэлектриков значения Vtotal и Vhetero практически не отличаются. Тогда как для металлов значения Vhetero соответствуют величинам, которые характерны для полной валентности атомов в родственных соединениях.
Особый интерес вызывают кристаллы, в структуре которых атомы металла занимают неэквивалентные позиции Вайкоффа и следовательно, способны проявлять различные валентности.
| Соединение группа симметрии | Позиция Вайкоффа | Q | Ns | Np/Nd | Chetero | Chomo | V |
| Pb2O3 | (4e) | 2.10 | 0.74 | 1.16 | 2.67 | 0.13 | 3.82 |
| P21/a | (4e) | 0.93 | 1.92 | 1.15 | 1.72 | 0.17 | 2.13 |
| Pb3O4 | (4d) | 2.60 | 0.54 | 0.86 | 2.17 | 0.09 | 3.90 |
| P42/mbc | (8h) | 0.80 | 1.85 | 1.36 | 1.85 | 0.11 | 2.14 |
| BaBiO3 | (1a) | 1.44 | 1.95 | 1.61 | 2.31 | 0.13 | 3.01 |
| P2/m | (1c) | 2.80 | 1.23 | 0.97 | 2.46 | 0.11 | 4.29 |
| Cu4O3 | (8c) | 0.97 | 0.05 | 9.98 | 0.08 | 0.05 | 1.02 |
| I41/amd | (8d) | 1.65 | 0.07 | 9.28 | 0.50 | 0.05 | 1.91 |
| YBa2Cu3O6 | (1a) | 0.91 | 0.10 | 9.99 | 0.19 | 0.00 | 1.00 |
| P4/mmm | (2g) ^ | 1.55 | 0.12 | 9.32 | 0.64 | 0.04 | 1.91 |
| YBa2Cu3O7 | (1a) | 1.65 | 0.04 | 9.30 | 0.96 | 0.03 | 2.20 |
| Pmmm | (2q) ^ | 1.56 | 0.04 | 9.40 | 0.58 | 0.33 | 1.88 |
Как следует из результатов, приведенных в таблице, смешанные оксиды Pb2O3, Pb3O4 содержат атомы двух и четырехвалентного свинца. Хотя электронное строение неэквивалентных атомов висмута в структуре BaBiO3 различно, отметим, что величины зарядов на атомах существенно отличаются от таковых, найденных для других соединений Bi(III) и Bi(V). Найденные значения атомных зарядов, ковалентностей и валентностей атомов меди, позволяют сделать однозначный вывод о наличии одно и двухвалентной меди в кристаллах Cu4O3 и YBa2Cu3O6. Найденные значения валентности для обоих типов атомов меди в YBa2Cu3O7 близки к двум. В то же время, следует учесть то обстоятельство, что и в других трехвалентных соединениях меди, вычисленная полная валентность принимает значения, меньшие трех. На основании результатов выполненного расчета электронной структуры YBa2Cu3O7 можно сделать вывод о том, что атом меди (находящийся в цепочках CuO), имеет электронное строение близкое к таковому в других формально трехвалентных соединениях меди (KCuO2, LaCuO3).
Найдено, что для высокоокисленных соединений меди и никеля индекс Вайберга связи O-O увеличивается до ~ 0.3, однако его величина существенно меньше порядков связей в пероксидных соединениях.
Изучены изменения, происходящие в электронном строении атомов sp и sd металла при увеличении его формальной степени окисления. На рисунках представлены зависимости полной валентности, а также ионной и ковалентной составляющей химической связи от формальной степени окисления свинца, висмута, никеля и меди. Ковалентности атомов обозначены звездочками и обведены прямоугольником, атомные заряды обозначены знаком "o" и обведены овалом. Полные валентности атомов обозначены заштpихованными кpугами. Прямая линия f(x)=x приведена для наглядного представления о величине отклонения полной валентности атомов от соответствующих стехиометрических значений.
Рис. 3. Зависимость QA , CA и VA от формальной степени окисления атомов Pb, Bi, Ni, Cu в кристаллических оксидах.
Как следует из приведенных рисунков величины QA и CA изменяются достаточно сильно, тогда как суммарное значение полной валентности атома изменяется в значительно более узких пределах. Приведенные графики позволяют проследить, как для соединений различных элементов реализуется определенное валентное состояние. Для никеля и висмута при увеличении формальной степени окисления ковалентность атома изменяется незначительно и увеличение валентности полностью определяется изменением атомного заряда. Для атомов свинца и меди происходит практически одновременный рост ионности и ковалентности. Введенное в настоящей работе определение валентности атома позволяет установить взаимосвязь между ионной и ковалентной составляющими химической связи, состоящее в перераспределении электронной плотности, локализованной на атомах или связях.
В работе рассмотрены возможные применения развитого аппарата для изучения химической связи в кристаллах, содержащих дефекты и их энергетической стабильности. Как показано в работе на примере большого числа идеальных кристаллов, даже при существенно различном кристаллическом окружении атомы стремятся проявлять определенные характерные для них значения валентностей. Следует ожидать, что после образования дефектов, электронная подсистема кристалла будет стремиться к "невозмущенному" положению, при котором локальные свойства стали бы близки к характерным для данного элемента значениям. Таким образом, анализ локальных свойств электронной структуры поверхности или кристаллов с точечными дефектами может использоваться для оценок устойчивости дефектов и даже поиска возможной равновесной геометрии дефект-содержащего кристалла. В качестве примеров рассмотрены точечные дефекты на поверхности графита (замещение атомов углерода на бор и азот) и изучены различные механизмы релаксации поверхности кристалла ZnO.
ВЫВОДЫ
Основные публикации по теме диссертации:
[1]Верязов В. А., Эварестов Р. А. Квантовохимическое определение валентности атома в молекуле или кристалле. // Вестн. ЛГУ,.-1990.-18,.- p. 38-43.
[2]Evarestov R. A., Veryazov V. A. The electronic structure of copper oxide crystalline compounds. I. LUC-CNDO approach to the electronic structure of Cu2O and CuO crystals. // Phys. Status Solidi(b).-1990.-157.- p. 281-291.
[3]Evarestov R. A., Veryazov V. A. The electronic structure of copper oxide crystalline compounds. II. Chemical bonding in copper-oxygen crystals. // Phys. Status Solidi(b).-1990.-158.- p. 201-212.
[4]Верязов В. А., Эварестов Р. А. Электронная структура дефектов на поверхности графита // Теоретическая и экспериментальная химия.-1990.-26.- p. 63-65.
[5]Верязов В. А., Китаев Ю. Э., Смирнов В. П., Эварестов Р. А. Электронные и фононные состояния в кристаллических высокотемпературных сверхпроводниках, в кн. Высокотемпературная сверхпроводимость / ред. А. А. Киселев, Л., "Машиностроение", 1990, стр. 446-512.
[6]Evarestov R. A., Veryazov V. A. Quantum-chemical definition of the atomic valence in molecules and crystals. // Theor. Chim. Acta.-1991.-81.- p. 95-103.
[7]Evarestov R. A., Veryazov V. A. Use of the cyclic model in the description of local properties of composite crystalline oxides. // Rev. of Solid State Science.-1991.-5.- p. 415-433. (in Quantum mechanical cluster calculations in solid state studies/ ed. by R.W. Grimes, C.R.A. Сatlow, A.L.Shluger, Singapore, World Scientific Publ., 1992, p. 355-373.)
[8]Evarestov R. A., Veryazov V. A. The electronic structure of crystalline lead oxides. I. Crystal structure and LUC-CNDO calculations. // Phys. Status Solidi (b).-1991.-165.- p. 401-410.
[9]Evarestov R. A., Veryazov V. A. The electronic structure of crystalline lead oxides. II. Chemical bonding in the crystalline lead oxides. // Phys. Status Solidi (b).-1991.-165.- p. 411-418.
[10]Верязов В. А., Эварестов Р. А. Химическая связь и локальные свойства кристаллических ВТСП // Физика низких температур,.-1991.-17.- p. 1254-1256.
[11]Evarestov R. A., Leko A. V., Murin I. V., Petrov A. V., Veryazov V. A. Electronic structure of La2O3 and LaF3 crystals. // Phys. Status Solidi (b).-1992.-170.- p. 145-153.
[12]Evarestov R. A., Veryazov V. A., Tupitsyn I. I., Afanasiev V. V. Electronic structure of crystalline nickel oxides, // J. Electronic Spectroscopy and Related Phenomena.-1994.-68.- p. 555-563.
[13]Evarestov R. A., Shapovalov V. O., Veryazov V. A. Electronic structure and chemical bonding in Bi2O3. // Phys. Status Solidi (b).-1994.-183.- p. K15-K17.