CC BY
36
®
УДК 538.975 ББК 22.3
НАНОТУБУЛЯРНЫЕ ФОРМЫ БОРА: ОСОБЕННОСТИ
ЭЛЕКТРОННО-ЭНЕРГЕТИЧЕСКОГО СТРОЕНИЯ И ПРОВОДЯЩИХ СВОЙСТВ
Е.В. Борознина, И.В. Запороцкова, С.В. Борознин, Д.И. Поликарпов
В данной работе представлены результаты исследования влияния структурных модификаций борных нанотруб на электронное строение и проводящие характеристики, а также результаты расчетов поверхностных, в том числе сорбционных свойств однослойных борных нанотубуле-нов. Обсуждаются особенности механизмов внешней адсорбции различных атомов на поверхности борной нанотрубки малого диаметра и вызванные ими изменения проводящего состояния поверхностно-модифицированных тубуленов. Расчеты выполнены в рамках моделей ионно-встро-енного ковалентно-циклического и молекулярного кластеров с помощью полуэмпирической схемы MNDO.
Ключевые слова: борные нанотрубки, полуэмпирические расчеты, адсорбция, протонная проводимость, молекулярный кластер.
о
(N «
«
о
я «
§
с pq
О я
СП
о &
РЧ
pq
К «
§
я о а о
00 pq
W л я
о &
W
©
Введение
В 2004 г. ученым удалось синтезировать чистые однослойные борные нанотрубки [4], которые могут быть рассмотрены как новый класс топологической структуры бора, обладающий уникальными физико-химическими и проводящими свойствами, открывающий широчайшие перспективы его использования в нано- и микроэлектронике. Ожидается, что борные нанотубулярные системы станут функциональными блоками и элементами (транзисторами, диодами и т. п.) для следующего поколения наноэлектронных устройств.
Однако до настоящего времени не существует однозначного мнения о точной структуре нанотубулярного бора, практически ничего не известно о его электронно-энергетическом строении, физических и химических свойствах. Поэтому исследования новой конфигурации бора, а именно его нанотубулярной модификации, в настоящее время являются
чрезвычайно актуальными ввиду ожидаемых перспектив его применения.
Исследование влияния структурных модификаций борных нанотруб на электронное строение и проводящие характеристики
Основным расчетным методом в работе является метод ИВ-КЦК в рамках схемы MNDO [2]. Ниже представлены матричные элементы оператора Хартри-Фока-Рутана для двухцентрового взаимодействия и одноцент-рового взаимодействия:
(A)
Faa =8 U +у P , ,
ЦУ ЦУ ЦЦ / J Ц v
(цу|ц'У) - 2 <ЦЦ'1 уУ)
(I) B( *A)
(B)
; zb •(h sBsI
X, a
+8„v[vA») - V«],
Цу A
1 (A) (B)
fA=ß^ - 2 ;; p- • ич
+
,,
+
где ц, V, ц', V' - атомные орбитали 5-,р- и d- типа, локализованные на атоме А; s, 1 - атомные орбитали, локализованные на атоме В;
игг - кинетическая и потенциальная энергия
электрона в атоме; 5 - символ Кронекера; Р является матрицей плотностей (матрицы порядков связей); вц1 - двухцентровой одноэлектронный остов-
ный резонансный интеграл; <цц'|пп>, <цп|1а> - одноцентровые и двух-центровое интегралы отталкивания электронов, соответственно; яв - атомные орбитали s-типа атома В; VJI} - маделунговский потенциал в точке А, создаваемый только теми атомами, которые находятся в области I; VJr} - маделунговский потенциал, создаваемый остальной областью нанотрубки. Были исследованы электронно-энергетические характеристики чистых однослойных нанот-руб на основе гексагонального бора. Рассмотрены фрагменты однослойных борных нанотрубок типа (п, п) и (п, 0) (п = 4, 5, 6, 9, 11, 12). В качестве геометрических моделей изучаемых тубуленов выбраны кластеры, содержащие п шестиатомных борных циклов по периметру трубки и два и более элементарных слоя вдоль ее оси. На рисунке 1 в качестве примера представлены расширенные элементарные ячейки (далее - РЭЯ) нанотрубок (6, 6) и (6, 0). На выбранные РЭЯ накладывались циклические граничные условия вдоль оси трубок. Длины связи между соседними атомами бора, полученные в результате оптимизации геометрии систем, получились равными 1,44 Е для всех рассмотренных типов нанотрубок.
б
Рис. 1. РЭЯ борных нанотруб: а - (6, 6)-тип; б - (6, 0)-тип
Результаты расчетов основных характеристик борных нанотруб приведены в таблицах 1 и 2. Анализ ширины запрещенной зоны DEg, позволил сделать вывод, что все они являются узкозонными полупроводниками, независимо от диаметра. Были вычислены энергии деформации Едеф как разность энергий ква-зипланарной РЭЯ из атомов бора и РЭЯ соответствующей нанотрубки, полученной в результате скручивания гексагонального бора.
Анализ результатов показал, что энергия деформации уменьшается с увеличением диаметра для борных трубок типа (п, п), в отличие от борных нанотруб (п, 0) типа (рис. 2). Это позволило сделать вывод, что образование нанот-руб типа (п, 0) путем скручивания квазипланар-ного гексагонального бора энергетически невыгодно и, соответственно, менее вероятно. Анализ зарядового состояния борных тубуленов установил их устойчивую симметрию - заряды на атомах бора практически нулевые.
Можно рассчитать концентрацию носителей заряда в собственных полупроводниках, к которым относятся рассмотренные нами БНТ, опираясь на известный
Таблица 1
Основные электронно-энергетические характеристики борных нанотрубок (п, п)
а
n Количество трансляций Удельная энергия, эВ Потенциал ионизации, эВ AEg,, эВ
4 4 68,01 6,43 0,07
5 4 67,91 6,07 0,04
6 4 68,14 6,32 0,90
9 3 68,97 7,55 0,02
11 2 69,34 7,73 0,35
12 2 69,33 7,84 0,22
Таблица 2
Основные электронно-энергетические характеристики борных нанотрубок (n, 0)-типа
n Количество трансляций Удельная энергия, эВ Потенциал ионизации, эВ ДЕ&,эВ
4 4 67,85 6,50 0,27
5 4 67,80 7,18 0,00
6 4 67,88 6,24 0,20
8 3 67,92 6,69 0,01
12 2 69,03 7,41 0,02
а б
Рис. 2. Зависимость энергии деформации Е от диаметра dтубулена: а - (n, п)-тиш; б - (n, 0)-типа
формализм теории твердого тела [1] и используя формулу:
n=^ й1 )з/2( ef - )з
3 h
где m * - эффективная масса электрона; h = 6,63 *110-34 ДжХс - постоянная Планка; Ef - энергия уровня Ферми; EC - энергия дна зоны проводимости.
Как известно, эффективная масса определяется кривизной зависимости Е^). В нашем расчетном методе энергетический спектр строится для точки k = 0, когда эффективная масса электрона фактически совпадает с массой свободного электрона. В результате для борных нанотруб (далее -БНТ) типа «arm-chair» нами получено значение концентрации основных носителей: n ~ 1020 см-3. Получаем, что БНТ малого диаметра в смысле проводимости являют-
ся системами, близкими к металлам, для которых концентрация носителей ~ 6х1022см-3.
Поверхностная модификация борных нанотруб как способ создания элементов аноэлектроники с регулируемой проводимостью
Известно, что дефекты существенно изменяют проводящие свойства углеродных на-нотруб. Интересно проследить за изменением свойств борных тубуленов с дефектами.
Изучено электронное строение и энергетические характеристики борной трубки типа (6, 6), включающей дефекты замещения. Дефект находился в центре кластера (см. рис. 3). В качестве дефекта выбирался либо нейтральный атом углерода (С), либо положительно (С+) и отрицательно (С ) заряженные ионы углерода.
Рис. 3. Молекулярный кластер борной нанотрубки типа (6, 6), содержащей одиночный дефект замещения атома бора поверхности атомом углерода
Анализ электронно-энергетической структуры показал, что дефект вносит дополнительные уровни в валентную зону, что приводит к ее расширению по сравнению с валентной зоной (далее - ВЗ) бездефектной борной нанотрубки. Величина запрещенной щели не изменяется, хотя меняется положение нижней вакантной и верхней заполненной орбиталей.
Итак, можно утверждать, что характер проводимости В-нанотрубки не изменяется при введении одиночного дефекта замещения.
Далее была исследована борная нанот-рубка (6, 6), содержащая вакансию. Процесс образования вакансии моделировался пошаговым удалением (с шагом 0,1 А) атома В от поверхности тубулена до момента его отрыва. В результате выполненных расчетов был построен профиль потенциальной энергии процесса образования вакансии, представленный на рисунке 4.
Рис. 4. Профиль потенциальной энергии процесса образования вакансии в борной нанотрубке типа (6, 6), нормированной к энергии отрыва атома (значение энергии системы ЕиЬВ = 0 на расстоянии RtuЬ_в = N
Из графика видно, что при отрыве атом В преодолевает потенциальный барь-
ер, отождествляемый с энергией активации дефекта Еакт=0,68 эВ. Найденное значение Еакт позволило оценить относительную долю вакансий при температуре Т= 1000 К по известной формуле:
n/N ~ exp[-Exp[-kT]
Она оказалась равной 0,53. Таким образом, при температуре 1000 К на 100 узлов борной нанотубки может приходиться 53 вакансии, при этом система должна сохранить свою устойчивость [1].
Анализ одноэлектронных спектров борных нанотруб с вакансией не обнаружил изменения ширины запрещенной зоны последних, что свидетельствует о низкой чувствительности проводящих свойств борных тубуленов к одиночным точечным дефектам.
Далее была исследована адсорбция атома водорода как способ создания носителя заряда на поверхности борного тубулена. Рассмотрены три варианта ориентации адсорбирующегося атома Н (адатома) относительно поверхности нанотрубки: I) над атомом бора, II) над центром связи В-В, III) над центром гексагона (рис. 5).
Рис. 5. Три варианта ориентации адатома (Н) относительно поверхности борного тубулена: I - над атомом бора; II - над центром связи В-В;
III - над центром гексагона
Процесс адсорбции моделировался пошаговым приближением адсорбирующихся атомов к поверхности борного тубулена. Выполненные расчеты позволили построить профили поверхности потенциальной энергии процессов адсорбции (см. рис. 6).
Оказалось, что во всех рассмотренных положениях атом Н активно адсорбируется на поверхности борной нанотрубки. Анализ зарядового состояния системы установил, что при адсорбции атома водорода имеет место перенос электронной плотности с атома Н на поверхность борной нанотрубки, то есть фактически образуется протон Н+.
Рис. 6. Профили поверхности потенциальной энергии взаимодействия атома Н с поверхностью борной нанотрубки типа (6, 6) для трех вариантов расположения его относительно борного тубулена:
I - над атомом бора; II - над центром связи В-В;
Ш - над центром борного гексагона
Таким образом, было обнаружено, что присоединение атома водорода к поверхности борного тубулена приводит к созданию внешнего носителя положительного заряда -протона, что позволяет предположить возможность использования нанотрубок на основе бора в качестве элементов наноэлектроники с поверхностной протонной проводимостью.
Для исследования мханизма протонной проводимости на поверхности борных нанот-руб были выбраны тубулены типа (6, 6) и (8, 8). Рассмотрены два варианта миграции одиночного протона Н+ вдоль поверхности нанотрубки (рис. 7):
1) так называемый «прыжковый» механизм, когда протон Н+ движется от одного атома бора поверхности до другого над двумя следующими друг за другом гексагонами (путь I);
2) «эстафетный» механизм, когда протон Н+ перемещается от одного атома бора до другого вдоль соединяющей их связи (путь II).
Рис. 7. Пути миграции протона по поверхности нанотруб типа «arm-chair» варианты I и II
Для обоих вариантов построены кривые энергий процессов миграции протона
(см. рис. 8). Установлено, что во всех случаях на кривых имеются максимумы, которые отождествляются с энергией активации (Е ).
Используя полученные значения величин потенциальных барьеров, преодолеваемых протоном при миграции вдоль поверхности, можно оценить подвижность носителей поверхностного положительного заряда. Для этого была выполнена следующая расчетная процедура.
Для начала мы произвели оценку доли частиц, обладающих достаточной энергией для преодоления потенциального барьера классическим способом по формуле:
а = exp
E
акт
' kT
Далее была оценена скорость поверхностной миграции протонов с использованием формулы:
где п - концентрация протонов.
В нашем случае скорость поверхностной миграции протона, отнесенная к значению концентрации протонов п, может считаться дрейфовой и использоваться для оценки подвижности протона, являющегося внешним носителем положительного заряда по формуле
¡л =
_др_ E
где Е - напряженность электрического поля.
Значения подвижности протонов для рассматриваемых борных нанотруб представлены в таблице 3.
Для сравнения подвижность электронов в кремнии по различным данным составляет (0,14...0,19) м2/(ВЧс), а в арсениде галлия -(0,93 ...1,1) м2/(ВЧс). Подвижность дырок в кремнии и арсениде галлия составляет (0,04...0,05) м2/(ВЧс) и 0,045 м2/(ВЧс), соответственно. Таким образом, подвижность протона на поверхности борных нанотрубок сравнима с подвижностью основных носителей типичных полупроводников.
v =
s
б
а
в г
Рис. 8. Потенциальная энергия миграции протона по поверхности нанотрубки: а - (6, 6) - путь II; б - (6, 6) - путь I; в - (8, 8) - путь II; г - (8, 8) - путь I
Таблица 3
Основные характеристики процесса продольной миграции протона по поверхности
борных нанотрубок типов (6, 6) и (8, 8) *
БНТ (6, 6) (8, 8)
Пути миграции протона I II I II
а 1,3х10-4 7,8 х10-2 1,7х10-2 2,0х10-2
Vs, с-1 м-2 1,5 х10-2 8,9 1,9 2,3
ц, м2/Вхс 1,5 х10-3 0,89 0,19 0,23
* Доля частиц (протонов), преодолевающих потенциальный барьер классическим способом а; скорость поверхностной миграции протона вдоль БНТ V, подвижность протона ц.
Далее была исследована регулярная гидрогенизация и оксидирование борной нанотрубки. Выполнены расчеты двух вариантов присоединения атомов водорода и кислорода к поверхности нанотрубки, как показано на рисунке (см. рис. 9). Оказалось, что с энергетической точки зрения группе атомов водорода и кислорода наиболее выгодно расположиться над противолежащими вершинами гек-сагонов, а не в виде цепочки.
Анализ полученных одноэлектронных спектров гидридов и оксидов борной нанотруб-ки и чистой БНТ, представленных на рисунке 10, позволил сделать вывод, что наличие атомов водорода и кислорода на поверхности изменяет проводящие характеристики В-тубуле-на, а именно, гидрогенизированная и оксидированная борная нанотрубка обладает металлической проводимостью: энергетическая щель для гидридов обоих вариантов равна нулю.
Рис. 9. Развернутые в плоскость РЭЯ борных тубуленов (6, 6) с указанием положений атомов водорода на поверхности: а) 1 вариант; б) 2 вариант
о -|-
Е. эв а) б) в)
0.0 — 2.0 0.0 — 2.0
-6--—-¡5-——
-8--=-
0.0
■10--=-20
-14 -1-
Рис. 10. Одноэлектронные энергетические спектры гидридов борных тубуленов (6, 6) для двух вариантов присоединения водорода:
а) вариант 1 несмещенного расположения атомов водорода относительно поверхности; б) вариант 2 смещенных друг относительно друга
колец атомов водорода; в) чистая борная нанотрубка. Показаны дважды заполненные и вакантные уровни
Итак, можно утверждать, что БНТ могут выступать в роли аккумуляторов водорода, наряду с углеродными нанотрубками, а гидриды и оксиды борных нанотруб могут использоваться в качестве элементов наноэлек-троники с переходом «узкозонный полупро-
водник - металл» (элементы с барьером Шотки), если полупроводящую борную на-нотрубку с одного конца насытить атомами водорода или кислорода, а другой конец оставить чистым.
Существенная модификация свойств тубуленов может быть достигнута при их покрытии разнообразными металлами и соединениями. Поэтому мы исследовали возможность присоединения атомов щелочных металлов К) к внешней поверхности борной нанотрубки типа (6, 6) с целью выяснения возможности создания металлофазных композитов на основе борной нанотрубки. Положения адсорбирующихся атомов изображены на рисунке 11.
I II
Рис. 11. Варианты ориентации адатомов относительно поверхности борного тубулена: I - над атомом бора, II - над центром гексагона
Анализ энергетических кривых показал, что в первом случае (адсорбция над поверхностным атомом В) все выбранные атомы металлов активно адсорбируются на поверхности В-тубулена (рис. 12).
г^ А
Рис. 12. Профили поверхности потенциальной энергии взаимодействия атомов Li, К, № с поверхностью борной нанотрубки типа (6, 6); вариант I - над атомом бора
Во втором случае (см. рис. 13) адсорбция наблюдается лишь для атомов калия и натрия.
Рис. 13. Профили поверхности потенциальной
энергии взаимодействия атомов Li, К, № с поверхностью борной нанотрубки типа (6, 6); вариант II - над центром гексагона
Значения энергий адсорбции Еа атомов Li, К, № и соответствующих расстояний адсорбции Rа представлены в таблице 4.
Таблица 4
Основные характеристики процессов адсорбции атомов лития Li, калия К и натрия ^ на поверхности борной нанотрубки (6, 6) *
Вариант ориентации адатомов (I) Вариант ориентации адатомов (II)
Е^эВ Яд, А Q ЕаД,эВ Яд, А Q
Li 2,79 2,3 0,73 - - -
K 8,30 2,0 0,74 1,56 2,1 0,63
Na 4,17 2,2 0,60 2,30 2,9 0,37
* Здесь Eaд - энергия адсорбции, эВ, Яад -расстояние адсорбции; Е, Q - заряды на атомах щелочных металлов; вариант ориентации I - над атомом бора поверхности, вариант ориентации II - над центром борного гексагона.
Проанализировано зарядовое перераспределение в полученных нанокомпозитах. Установлен факт переноса электронной плотности с атомов металлов на борную трубку, что увеличивает число основных носителей заряда в тубулене, а следовательно, усиливает металлические свойства борных трубок, которые мы ранее отнесли к узкощелевым полупроводникам.
Далее было исследовано регулярное расположение атомов щелочных металлов (далее - ЩМ) в позициях, показанных на рисунке 14.
Данные системы, выполненные атомами щелочных металлов, можно условно считать
металлической сверхрешеткой прямоугольного (рис. 14 (I)) и ромбического (рис. 14 (II)) вида. При выбранных положениях атомов металлов постоянные сверхрешеток оказываются сравнимыми с экспериментальными значениями постоянных кристаллических решеток ЩМ [3].
Рис. 14. Расширенная элементарная ячейка борной нанотрубки типа (6, 6) с указанными положениями атомов металлов (Ь^ К, №), образующих над поверхностью трубки сверхрешетки I - прямоугольного типа, II - ромбического типа
В результате проведенного исследования мы получили значения ширины запрещенной зоны ДЕg , представленые в таблице 5. Для варианта модифицирования I ширина запрещенной зоны борных нанотруб, модифицированных атомами Ь^ по сравнению с чистым борным тубуленом в пределах заданной точности не изменяется. Для тубуленов, модифицированных атомами К и происходит незначительное увеличение запрещенной зоны, но тем не менее, данные структуры по-прежнему можно классифицировать как узкозонные полупроводники. Для варианта модифицирования II ширина запрещенной зоны для БНТ, модифицированных атомами Li и увеличивается незначительно, а в борном тубулене с атомами К энергетическая щель не изменяется по сравнению с энергетической щелью чистой БНТ.
Анализ зарядового состояния для обоих вариантов расположения атомов металлов над поверхностью БНТ установил, что так же, как и в случае одиночной адсорбции, происходит перенос электронной плотности с атомов щелочных металлов на поверхностные атомы бора (см. табл. 5).
Таблица 5
Энергетические характеристики модифицированных борных тубуленов типа (6, 6) *
Вариант модифицирования б Е эВ А Е, эВ А, эВ
Чистая БНТ - - -207,00 0,15 -
БНТ+12Ы I 0,96 -212,68 0,16 1,37
II 0,95 -211,31 0,39
БНТ+12К I 0,96 -207,88 0,35 1,44
II 0,94 -206,44 0,14
БНТ+12№ I 0,92 -212,88 0,28 1,49
II 0,97 -211,39 0,20
* I, II - варианты модифицирования; Q - заряды на атомах щелочных металлов; Е - энергия системы, приходящейся на один гексагон, эВ; АЕ, - ширина запрещенной зоны, эВ; А = - разность полных энергий систем, эВ
Таким образом, расчеты показывают, что поверхностное модифицирование борных нанотруб атомами металлов не меняет проводящий характер борных нанотруб, но приводит к возникновению поверхностной проводимости за счет появления свободных носителей заряда на поверхности трубки, тем самым обеспечивая появление перехода «узкозонный полупроводник - металл» в полученных таким образом металлофазных композитах борных нанотруб.
Заключение
1. Определены основные характеристики квазипланарного гексагонального бора и соответствующих нанотубулярных форм. Доказана возможность преимущественного образования (п, п) борных нанотруб, причем БНТ малого диаметра (до 2 нм) являются узкозонными полупроводниками. Установлено, что БНТ малого диаметра являются эффективными в смысле проводимости системами, близкими к металлам по числу носителей - электронов.
2. Одиночные дефекты не влияют на характер проводимости борных нанотруб. Это позволит при конструировании приборов наноэлектроники широко использовать борные нанотрубки произвольно выбранного диаметра и «чистоты» поверхности без
потери искомых проводящих свойств последних.
3. Исследован механизм адсорбции атома водорода на внешней поверхности борных нанотруб. Выяснено, что при адсорбции происходит перенос электронной плотности с атома Н на поверхность тубулена, что фактически свидетельствует об образовании свободного носителя положительного заряда - протона Н+.
4. Исследованы «эстафетный» и «прыжковый» механизмы миграции протона Н+ по внешней поверхности борной на-нотрубки. Проведена оценка подвижности протонов и установлено, что она сравнима с подвижностью основных носителей типичных полупроводников. Доказано, что БНТ могут быть отнесены к классу протонпро-водящих материалов, что обеспечит новые интересные перспективы их использования в электронной технике.
5. Исследованы механизмы адсорбции атомов газовой фазы на поверхности борных нанотруб и определены наиболее предпочтительные места их присоединения. Тип проводимости получившихся фторида и хлорида борного тубулена не изменяется по сравнению с чисто борной нанотрубкой и сохраняет полупроводящий характер.
6. Изучен процесс регулярной адсорбции атомов водорода и кислорода на поверхности В-нанотруб. Обнаружено изменение
проводимости получившихся композитов в сторону металлизации по сравнению с полупроводящим состоянием чистой БНТ.
7. Исследован механизм адсорбции атомов щелочных металлов на поверхности борной нанотрубки и установлено наиболее предпочтительное место их присоединения. Доказано, что модификация поверхности БНТ атомами металлов приводит к возникновению металлической сверхрешетки и появлению перехода «узкозонный полупроводник - металл» в полученных таким образом металлофазных композитах борных нанотруб.
СПИСОК ЛИТЕРАТУРЫ
1. Китель, Ч. Введение в физику твердого тела / Ч. Китель. - М. : Наука. - 1978. - 79 с.
2. Литинский, А. О. Модель ионно-встроен-ного ковалентно-циклического кластера в MNDO-расчетах межмолекулярных взаимодействий в гетерогенных системах / А. О. Литинский, Н. Г. Лебедев, И. В. Запороцкова // Журнал физической химии. - 1995. - Т. 69. - №№1 - С. 189.
3. Пенкаля, Т. Очерки кристаллохимии / Т. Пен-каля. - Л. : Химия. - 1972. - 496 с.
4. Ciuparu, D. Synthesis of Pure Boron SingleWall Nanotubes / D. Ciuparu [et al.] // J. Phys. Chem. B. - 2004. - V 108. - P. 3967-3969.
NANOTUBULYARNYE OF THE PINE FOREST FORM: FEATURES OF THE ELECTRONIC AND POWER STRUCTURE AND CARRYING-OUT PROPERTIES
E. V. Boroznina, I. V. Zaporotskova, S. V. Boroznin, D.I. Polikarpov
In this work results of research of influence of structural modifications of boric nanopipes on an electronic structure and carrying-out characteristics, and also results of calculations superficial, including sorption properties single-layer boric HaHOTy6yjieHOB are presented. Features of mechanisms of external adsorption of various atoms on surfaces of a boric nanotube of small diameter and the changes of a carrying-out condition caused by them superficial modified Ty6yjieHOB are discussed. Calculations are executed within models of ion-built-in covalent and cyclic and molecular clusters by means of the semi-empirical scheme of MNDO.
Key words: boric nanotubes, semi-empirical calculations, adsorption, proton conductivity, molecular cluster.
