Одним из семейств сложных и необычных углеродных кластеров являются тороидальные формы углерода. Такого вида молекулы были экспериментально зафиксированы в 1992 году [1]. Они представляют собой замкнутые поверхности, состоящие из углеродных многоугольников (пяти-, шести- и семиугольники). Наличие в их составе пятиугольников обуславливает положительное искривление поверхности, в то время как семиугольники, наоборот, обеспечивают отрицательную кривизну. Подобные молекулы имеют довольно сложную форму (рисунок 8.1, 8.2). Они могут содержать различное количество атомов (от 80 до нескольких тысяч).
Рисунок 8.1. Единичные структуры исследуемых тороидальных молекул, изображенных на Рис. 10.2.
a/) С120 (elong.) (С120elong_MNDO.hin),
b/)C120 (short) (С120short_MNDO.hin),
c/) C240 (С240_MNDO.hin). Неэквивалентные атомы обозначены буквами,а связи между ними пронумерованы.
Рисунок 8.2. Тороидальные молекулы и их комплексы с литием, единичные структуры которых изображены на Рис. 10.1
a/) С120 (elong.) (С120elong_MNDO.hin),
b/) C120 (short) (С120short_MNDO.hin),
c/) C240 (С240_MNDO.hin),
d/) Li2@C120(elong.) (Li2@С120elong_MNDO.hin),
e/) Траектория движения атомов лития в процессе молекулярного моделирования комплекса Li2@C120(elong.) (Li2@С120ih_300K_Lang.hin),
f/) Li4@C120(elong.) (Li4@С120elong_MNDO.hin).
В период с 1992 года по настоящее время в печати опубликована целая серия работ, посвященная теоретическим исследованиям тороидальных углеродных молекул. Так С. Ито [2-4], а также К. Джонсон [5] в своих работах изучали зависимости энергии кластеров и их электронной структуры от количества атомов и геометрических параметров объекта исследования (симметрия, размеры внутреннего и внешнего диаметров и т.д.). С. Ито провел подобные исследования для изомеров С240 [6]. В ходе своих исследований С. Ито для оптимизации геометрии кластеров использовал молекулярную динамику, а для изучения электронной структуры метод Хюккеля. В.Данлоп [7] и В.Мейнер [8] осуществляли изучение атомной и электронной структуры отдельных тороидальных молекул. Так, Данлопом рассматривались углеродные торы содержащие 240, 540 и 576 атомов. По результатам расчетов были получены значения ширины запрещенной зоны 1, 0.04 и 0.02 эВ соответственно для каждой из представленных структур. Мейнером был исследована молекула С1960, для которой он получил ширину запрещенной щели 0.05 эВ. При этом оба автора использовали в своих работах для оптимизации геометрии вычислительные методы, основанные на эмпирических потенциалах, а для получения электронной структуры модель сильной связи. В общем случае отличия тороидальных молекул в электронной структуре обуславливаются их различиями в их атомных структурах.
Однако, гораздо больший интерес с точки зрения особенностей электронной структуры представляют производные тороидальных молекул, в частности, их металлокомплексы. Подобная задача до сих пор не рассматривалась, несмотря на то, что тороидальные молекулы, имеющие во внутренней полости атомы металлов, их электронная структура, а так же поведение внедренных атомов чрезвычайно интересны.
В данной Главе проводилось теоретическое изучение методами квантовой химии и молекулярной динамики тороидальных форм углерода (содержащих 120 и 240 атомов углерода (рисунок 10.1, 10.2)) и их электрондопированных производных (эндоэдральных комплексов с литием), с целью выявления закономерностей в пространственной и электронный структуре данных молекул. В случае соединения С120 рассматривались два различных изомера, отличные друг от друга не только расположением в них пяти- и семиугольных циклов, но и своими размерами (внутренним и внешним диаметрами). Структура соединений С240 (C240_MNDO.hin) и С120 (elongate) (C120elong_MNDO.hin) были предложены в работах [2-3]. В отличие от вышеупомянутого изомера С120, другой был назван нами (short) (C120short_MNDO.hin). Координаты неэквивалентных атомов данного изомера представлены в Таблице 1.
Таблица 1
Координаты неэквивалентных атомов изомера С120 (short), рассчитанные методом MNDO, симметрия молекулы – D5d (Å) (C120short_MNDO.hin).
X
Y
Z
0.0000000000
4.6992476708
1.2477305295
0.0000000000
6.0552438387
1.6057198707
0.7582005602
7.1114336400
0.2071982406
0.6784827003
5.1688426852
1.4971372533
1.5186476126
6.4116043234
-1.2346751306
1.2622942952
6.7170404340
1.1064264943
1.2650400899
4.1822375962
0.7031065629
Помимо самих тороидальных структур изучались их металлокомплексы с литием (рисунок 8.2 D, F). Рассматривались соединения с двумя и четырьмя атомами металла внутри углеродных остовов (Li2@С120elong_MNDO.hin, Li2@С120short_MNDO.hin, Li4@С120elong_MNDO.hin, Li4@С120short_MNDO.hin, Li2@С240_MNDO.hin, Li4@С240_MNDO.hin). В данном случае внимание уделялось положению атомов металла и их поведению при различных температурах (в молекулярно-динамических расчетах моделировалась температура 300К), а так же влиянию внедренных атомов на электронную структуру углеродных торов.
Расчеты осуществлялись полуэмпирическим (MNDO) (С120elong_MNDO.hin, С120short_MNDO.hin, С240_MNDO.hin, Li2@С120elong_MNDO.hin, Li2@С120short_MNDO.hin, Li4@С120elong_MNDO.hin, Li4@С120short_MNDO.hin, Li2@С240_MNDO.hin, Li4@С240_MNDO.hin) и неэмпирическим (базис 3-21G) (С120elong.hin и С120shot.hin) методами. Причем последний применялся только в случае расчетов изомеров С120, вследствие большого числа базисных функций для крупных молекул. Расчет тороидальных соединеий осуществлялся с использованием симметрийных ограничений С120(elong.), С240 – С5, С120(short) – D5d, а в случае их металлокомплексов, подобные ограничения не накладывались, что позволило получить координацию атомов металла, соответствующую структуре с минимальной энергией. При этом все расчеты проводились как для синглетных, так и для триплетных состояний. В качестве проверки достижения минимума энергии рассчитывалась матрица гессиана (набор вторых производных по координатам ядер) и колебательный спектр соединений (минимум энергии характеризуется тем, что все частоты колебаний действительные, если был найден не минимум, а седловая точка, одна из частот становится мнимой). Все вышеупомянутые расчеты осуществлялись с помощью программы GAMESS [9].
Поведение атомов лития в молекулах рассматривалось методом молекулярной динамики. В настоящее время молекулярная динамика, не требующая введения эмпирических межмолекулярных и межатомных потенциалов для расчетов, широко используется для исследований динамических свойств молекулярных систем (см. Главу II). В данной работе для моделирования поведения атомов лития в комплексе использовалась демонстрационная версия программы HyperChem 5.02, позволяющая проводить расчеты с использованием полуэмпирических (MNDO, PM3) потенциалов.
Прежде всего, хотелось бы отметить отличия в электронной и пространственной структуре исследуемых молекул при использовании полуэмпирических и неэмпирических методов. Исходя из того, что в случае бакминстерфуллерена, наиболее точные значения длин связей были получены с помощью полуэмпирических методов (PM3, MNDO), можно сделать заключение, о том, что использование в неэмпирических расчетах базиса 3-21G приводит к заниженным значениям длинны связи (как впрочем, и в С60). При этом наблюдалась определенная градация изменений: наиболее подвержены им, оказались связи c порядком близким ко второму (1.7) (изменение длинны около 0.04 Å), и наоборот, в случае если порядок связи приблизительно равен единице, то ее длина практически не менялась (таблица 2, рисунок 10.1).
Таблица 2
Межатомные расстояния между неэквивалентными атомами в различных тороидальных молекулах (Å)
Номер связи
С120elong
MNDO
C120elong
3-21G
C120short
MNDO
C120short
3-21G
C240
1
1.477
1.474
1.471
1.454
1.433
2
1.378
1.333
1.417
1.364
1.413
3
1.485
1.472
1.513
1.507
1.469
4
1.386
1.341
1.357
1.312
1.371
5
1.466
1.447
1.523
1.531
1.448
6
1.420
1.385
1.457
1.425
1.489
7
1.479
1.464
1.516
1.539
1.423
8
1.477
1.461
1.461
1.415
1.427
9
1.477
1.458
1.510
1.520
1.450
10
1.411
1.374
1.434
1.405
1.388
11
1.495
1.490
1.402
1.346
1.478
12
1.466
1.447
1.483
13
1.485
1.472
1.498
14
1.431
15
1.456
16
1.483
17
1.470
18
1.407
19
1.494
20
1.468
21
1.482
22
1.395
Изменения в электронной структуре оказались не менее существенными. В полуэмпирических расчетах электронные уровни рассматриваемых систем оказались существенно ближе друг к другу, нежели в неэмпирических (Рисунок 8.3). Однако, несмотря на подобный сдвиг, их структура в области потолка валентной зоны осталась практически неизменной. Тем не менее, это привело к различной оценке ширины запрещенной щели (таблица 3) полуэмпирическим и неэмпирическим методами (отличие для обоих изомеров С120 составило 1.35 эВ).
Рисунок 8.3. Плотности состояний тороидальных структур.
a/) С120 (elong) C120 (short) – MNDO,
b/)С120 (elong) C120 (short) – 3-21G,
c/) С120 (elong) – MNDO, 3-21G (оптимизация происходила методом MNDO), 3-21G,
d/)d) C120 (short) – MNDO, 3-21G (оптимизация происходила методом MNDO), 3-21G.
Не менее интересным оказался сравнительный анализ различных изомеров С120 между собой (таблица 3). Структура С120 short оказалась энергетически менее выгодной (на 66.8 эВ по оценке MNDO и 58 эВ для 3-21G), вследствие меньшего диаметра внутренней полости, и как результат, более напряженной структуры. При этом, в молекуле С120 short значения длин связей в соответствующих многоугольниках оказались существенно выше по сравнению с С120 elongate.
Таблица 3
Энергия на атом и ширина запрещенной щели для различных тороидальных молекул и их комплексов с атомами лития
Объект исследования
C120(elong.)MNDO
C120(elong.)3-21G
C120(elon.)3-21G_optMNDO
Li2@C120(elong.)
Li4@C120(elong.)
C120(short)MNDO
C120(short)3-21G
C120(short)3-21G_optMNDO
Li2@C120(short)
Li4@C120(short)
C240
Li2@C240
Li4@C240
E/атом. эВ
127.18
1024.34
1024.31
125.19
123.27
126.62
1023.85
1023.81
124.57
122.58
127.27
126.29
125.31
DEho-lu эВ
6.454
7.809
7.138
1.611
2.114
5.828
7.181
6.201
2.686
4.634
4.422
4.304
3.603
Однако, несмотря на это, пропорциональные отношения между связями, принадлежащими различным циклам, в обоих изомерах оказались очень близки. В случае С240 семиугольники и пятиугольники непосредственно не соединяются как в С120, между ними расположены шестичленные циклы, вследствие этого длина связи в некоторых многоугольниках существенно изменилась (таблица 2, рисунок 8.1). Некоторые особенности, характерные для межатомных расстояний в С120, сохранились и в С240. Так наиболее коротким оказались связь между семи- и шестиугольниками (связь 4 у всех молекул), а так же направленная от вершины пятиугольника (связь 10 в С120 и 22 в С240).
В соединениях с двумя литиями, после процедуры оптимизации геометрии, атомы внедренных металлов расположились на максимально удаленном расстоянии друг от друга. В изомере С120 (elong.) это расстояние составило 7.79 Å, при этом они координировались напротив связи 5 (рисунок 10.1, рисунок 10.2 D) углеродного остова. В С120 (short) – 10.9 Å, атомы лития разместились в полости около связи 11. В комплексе Li2@С240 расстояние между металлами оказалось равным 10.5 Å с расположением их возле связи 22. Необходимо отметить, что во всех случаях литии находились практически на равном удалении от стенок углеродного части комплекса.
В случае четырех внедренных атомов лития, их координация в двух изомерах С120 была различной. Так в комплексе Li4@C120(elong.) щелочные металлы оказались сгруппированы попарно, расстояние между атомами в паре составляет 4.48 Å, а между парами – 6.44 Å (рисунок 2 F). Для изомера Li4@C120(short) наиболее характерной стала конфигурация с максимально удаленными друг от друга литиями (8 Å между ближними атомами и 11.6 Å между противоположными). В Li4@C240 атомы металла так же удалились на максимальное расстояние (7.5 Å между атомами – соседями и 10.39 Å между противолежащими литиями).
Расчеты показали, что введение внутрь углеродной структуры атомов – доноров электронов, какими являются щелочные металлы, вызывает искажение самого остова. Поскольку валентный электрон лития практически полностью переходит на углеродную часть комплекса, речь, в данном случае, идет о влиянии на связи углерод – углерод двух факторов. Появление дополнительных электронов, внесенных в систему щелочными металлами и воздействие положительно заряженных ионов лития (заряд ~0.6). Увеличение длин углеродных связей в металлокомплексах по сравнению с чистыми тороидальными молекулами указывает на антисвязывающую природу вакантного состояния, на котором располагаются электроны металлов. Максимально увеличивается длинна связей с большим порядком. При этом изменение длин связей имеет локальный характер, вследствие большой протяженности тороидального остова. По мере удаления от мест координации щелочных металлов изменения в углеродных связях сходят на нет.
Процесс интеркалляции щелочными металлами приводит и к изменениям электронной структуре тороидальных форм углерода. В отличие от чистых углеродных торов, для которых более низкими по энергии являются синглеты, для комплексных соединений С120 с двумя или четырьмя интеркалированными атомами лития, за исключением Li4@C120(short), триплетные состояния оказались, согласно проведенным расчетам, более низкими по энергии. При этом разница между триплетом и синглетом для комплекса Li2@C120(elong.) составила 1.255 эВ, для Li4@C120(elong.) – 0.054 эВ, Li2@C120(short) – 0.061 эВ и Li4@C120(short) – 0.343эВ.
Рисунок 8.4. Структура электронных уровней валентной зоны молекул a)С120 (elong.), b)C120 (short), c)C240.
Для того чтобы понять причину подобного явления, необходимо рассмотреть вакантные уровни изомеров С120. Нижним вакантным состоянием является А, а следующим за ним по энергии – двукратно вырожденное Е1 (рисунок 10.4 А, В). В случаях комплексов с двумя внедренными металлами происходит инверсия данных уровней и два электрона, внесенных атомами лития, заполняют электронные состояния углеродного каркаса комплекса, располагаясь на двух орбиталях, образовывавших в чистом углеродном торе нижнее вакантное состояние Е1 (формально двукратно вырожденного уровней не существует поскольку расчет производится в симметрии С1), что в соответствии с правилом Хунда приводит к появлению триплета. Процесс инверсии уровней наблюдается и для Li4@C120(short). Как и в предыдущем случае, роль внедренного состояния играет пара уровней, на которых размещаются четыре электрона, принадлежащих атомам металла, образуя синглет. В комплексе Li4@C120(elong.) подобного изменения в расположении уровней не происходит, и два из четырех электронов заполняют А, а два находятся на паре уровней (бывших Е1) делая основным триплетное состояние. Состояния с более высокой мультиплетностью не могут существовать из-за особенностей электронной структуры.
В металлокомплексах С240 расчеты показали, что синглетные состояния более выгодны по энергии в сравнении с триплетными (на 0.377 эВ для Li2@С240 и 0.101 эВ для Li4@С240). Таким образом, в данном случае предпочтительнее оказывается последовательное заполнение однократно вырожденных уровней А (рисунок 10.4 С).
Во всех изучаемых комплексах, уровни оказавшиеся заполненными, в результате появления электронов лития, располагаются в запрещенной щели. Таким образом, переход на углеродную часть комплекса электронов от металла приводит к появлению в запрещенной щели внедренного состояния (рисунок 8.5). Подобный процесс характерен и для металлокомплексов С60 [10]. Сдвиг электронных уровней под влиянием положительных зарядов атомов лития, в тороидальных молекулах менее заметен чем в С60, вследствие большого размера углеродной составляюшей металлокомплекса. Однако можно предположить, что инверсия уровней, приводящая к триплетному состоянию появляется как раз в результате воздействия данного фактора. Это подтверждается расчетами изоэлектронного иона , в котором инверсии не наблюдается.
Рисунок 8.5. Плотности состояний металлокомплексов Lin@Cm.
a/) С120 (elong.), Li2@C120(elong.), Li4@C120(elong.),
b/) C120 (short), Li2@C120(short), Li4@C120(elong.),
c/) Электронная структура комплекса Li2@C120(elong.) соответствующая максимальному и минимальному расстоянию между атомами лития в процессе молекулярного моделирования.
d/) C240, Li2@C240, Li4@C240.
Как видно из рисунка 8.5 (А, В, D) положение внедренного состояния зависит от количества атомов металла. При этом в комплексах С240 (рисунок 5 D) по мере увеличения количества интеркаллированых атомов внедренное состояние сдвигалось по направлению к вакантным уровням. Подобное движение внедренного состояния приводит к уменьшению запрещенной щели (таблица 3). Для С120(short) наблюдается обратная картина (рисунок 8.5В, таблица 3): внедренное состояние располагается ближе к вакантным уровням в соединении Li2@C120(short). В комплексах С120 (elong.) как с двумя атомами лития, так и с четырьмя (рисунок 5А, таблица 3), наблюдается сдвиг внедренного состояния, приводящий к практически полному замыванию запрещенной щели (1.8 эВ). Таким образом, появление в электронной структуре внедренного состояния является общим результатом для всех литийинтеркаллированых тороидальных соединений, изучаемых в работе.
Поведение атомов щелочных металлов внутри комплексов изучалось методом молекулярной динамики. Расчет производился для комплексов с двумя и четырьмя литиями, причем как для триплетных, так и для синглетных состояний.
В случае изомера С120 (short) (Li2_c120nash_300K.hin) при любом количестве внедренных атомов, а так же каждом из значений мультиплетности, атомы – гости совершали незначительное движение в потенциальной яме возле минимума энергии, найденного при оптимизации геометрии. При этом они не покидали пространства ограниченного шестиугольником d, e, f (рисунок 8.1). Данный факт легко объясняется малым внутренним диаметром углеродного остова данных комплексов. Кроме того, тороидальная молекула подобного вида состоит из своего рода полостей, с диаметром ~ 2.8 Å, которые разделены перемычками, диаметр которых составляет всего 2.2 Å Данные перемычки и играют роль барьера, удерживающего атомы лития в строго ограниченной области внутри углеродной молекулы.
В изомере С120 (elong.) (синглетное основное состояние) (Li2@C120ih_300K_Lang.hin, расчет проведен в приближении ланжевеновской механики) картина, описывающая поведение атомов – гостей существенно меняется. Поскольку данный изомер имеет существенно больший внутренний диаметр (~4 Å) и на его внутренней поверхности отсутствуют узлы как в С120 (short), внедренные щелочные металлы могут передвигаться по полости углеродного остова без существенных помех. Тем не менее, внутренняя поверхность молекулы имеет два вида потенциальных минимумов. Один из них соответствует положению атомов лития, полученному в процессе оптимизации геометрии. Второй минимум, более высокий по энергии, находится в центре шестиугольника d, e, f (рисунок 8.1). При моделировании поведения атомов лития в комплексах Li2@C120 в триплетном состоянии (Li2@c120ih_M3_MD.hin) и Li4@C120 (Li4@c120ih_300K.hin) атомы – гости переходили из одного вида минимумов в другой. Подобные переходы имели место только между соседними минимумами, и существенных изменений в положении внедренных атомов не наблюдалось.
Однако в случае синглетного состояния Li2@C120 (Li2@c120ih_300K_lang.hin) наблюдалась принципиально иная картина. Внедренные щелочные металлы, переходя из одного минимума в другой, передвинулись друг к другу и образовали пару с расстоянием между атомами ~ 3Å. В дальнейшем они стали совершать совместные переходы между соседними минимумами, мигрируя вдоль внутренней полости. При этом в паре происходили колебания расстояния Li – Li, величина которого при максимальном удалении атомов друг от друга не превышала 6.3Å. Подобное поведение и как следствие возникновение ненулевого орбитального момента ионов в молекуле при повышении температуры объясняется изменением полного электронного орбитального момента. Действительно, полный момент системы, который является совокупностью электронного орбитального момента, ядерного момента и спина, должен сохраняться. Однако при повышении температуры изменяются эффективные межъядерные расстояния, что приводит к изменению и электронной волновой функции системы: , так как при изменении температуры изменяется набор координат ядер (), которые входят в полноэлектронную волновую функцию в виде параметров. Следовательно, вся система должна компенсировать изменение электронного орбитального момента изменением ионного орбитального момента системы (поскольку ее мультиплетность не меняется), что и проявляется в движении атомов лития в синглетном состоянии. Отличие спиновой компоненты от нуля должно приводить к изменению пропорции между электронным и ионным орбитальными моментами и должно приводить к изменению температуры начала движения атомов внутри полости. Это объясняет отсутствие движения литиев в комплексе Li2@C120 (триплет) (Li2@c120ih_M3_MD.hin). Подобный эффект наблюдается и при изменении количества внедренных атомов.
Таким образом, в случае Li2@C120 (синглет) при различных температурах наблюдаются различные атомные конфигурации, которые будут иметь и отличную друг от друга электронную структуру. Так на рисунке 10.5С представлены электронные структуры соответствующие максимальному и минимальному расстояниям между атомами лития. В случае температурных перестроек в электронных спектрах должна наблюдаться картина суперпозиции электронных структур, отвечающих различным атомным конфигурациям, и как следствие некоторые спектральные особенности должны иметь размытую форму. В частности, это происходит с внедренным состоянием. При удалении атомов лития друг от друга внедренное состояние сдвигается к вакантным орбиталям, и наоборот, при сближении атомов металла внедренное состояние перемещается от вакантных уровней. Разница в положении внедренного состояния составляет ~ 2 эВ.
Однако поскольку без учета корреляционных эффектов мы не можем достоверно сказать о мультиплетности системы, то, интерпретируя результаты температурной зависимости экспериментальных фотоэлектронных спектров можно судить о спиновом состоянии комплекса.
В комплексах углеродного тора С240 как в синглетном, так и в триплетном состоянии и с различным количеством атомов лития (два и четыре), существенного перемещения внедренных атомов не наблюдалось (Li2@c240_300K.hin и Li2_non_diag_@c240_M1_300K.hin). Это возможно объяснить существованием на внутренней поверхности углеродного каркаса комплекса глубоких потенциальных минимумов, отсутствующих в С120, вследствие меньшего диаметра и как следствие более напряженной структуры.
Таким образом, мы увидели, что тороидальные углеродные молекулы между собой существенно отличаются как атомной, так и электронной структурой. Такие факторы как количество атомов в кластере, его симметрия, размеры (внутренний и внешний диаметры), порядок соединения циклов между собой, обуславливают строение электронных уровней и ширину запрещенной щели.
Введение во внутреннюю полость углеродной молекулы атомов щелочного металла приводит к родственным для всех изучаемых в работе кластеров изменениям в атомной и электронной структуре. В присутствии атомов лития увеличиваются межатомные расстояния в углеродном остове комплекса. В электронной структуре подобных металлокомплексов наблюдается внедренное состояние в запрещенной щели. Его положение (приближение или удаление от вакантных орбиталей) зависит от большого количества параметров (количества атомов, строения исходной углеродной молекулы, мультиплетности и т.д.). Однако, в целом, влияние атомов металла на электронную структуру комплексов можно свести к уменьшению ширины запрещенной щели за счет индуцирования внедренного состояния. При этом, как показывает результат молекулярной динамики, при соблюдении некоторых условий (достаточный размер внутренней полости, отсутствие на ее поверхности высоких барьеров или глубоких минимумов, а так же находжение системы в синглетном состоянии), в экспериментальных спектрах будет наблюдаться размывание пика, соответствующего внедренным уровням, за счет движения ионов металлов внутри углеродных торов.
Другие возможные торы можно посмотреть в директории Databank/Toruses файлы 300.hin, 320.hin, 360_prim.hin, 400.hin. Подвижность атомов гелия в одном из торов С360 иллюстрируется молекулярно-динамическим расчетом Tor_He8_MM+_300K.hin. Молекулярно-динамическое моделирование проводилось в модельном потенциале молекулярной механики.
СПИСОК ЛИТЕРАТУРЫ К ГЛАВЕ III.8.
1. B.M. Terrones, Phil. Trans. R. Soc. Lond. A, 354, 2025 (1992).
2. S. Itoh, S. Ihara, Phys.Rev., 47,1703 (1993).
3. S. Itoh, S. Ihara, Phys.Rev., 47, 12908 (1993).
4. S. Itoh, S. Ihara, Phys.Rev., 48, 8323 (1993).
5. J.K. Johnson, Phys.Rev., 50, 17575 (1994).
6. S. Itoh, S. Ihara, Phys.Rev., 49, 13970 (1994).
7. B.I. Dunlap, Phys.Rev., 46, 1933 (1992).
8. V. Menner, Phys.Rev., 57, 14886 (1998).
9. M.W. Schmidt, K.K. Baldridge, J.A. Boatz. J. Comp. Chem. 14, 6, 1347 (1993)
10. С.А. Варганов, П.В. Аврамов, С.Г. Овчинников. Физика твердого тела, т. 55, 378 (2000).
Subscribe to:
Post Comments (Atom)
No comments:
Post a Comment