Известно, что создание Li-ионных батарей было успешным в большой степени благодаря высокой ионной проводимости в структуре LiCoO2 [1]. Предсказание свойств материалов с помощью компьютерных вычислений позволяет значительно сократить время разработки новых материалов. В данном исследовании с использованием кристаллохимического анализа (программный пакет TOPOS) были проанализированы 2367 неорганических Na-содержащих соединений. Выявлено, что 192 соединения способны обладать ионной проводимостью при комнатной температуре. Из них 43–1D ионные проводники, 55–2D, 94–3D. Результаты могут быть использованы в дальнейших расчетах количественных характеристик ионной проводимости из первых принципов и создании новых ионных проводников.
Ключевые слова: кристаллохимический анализ, твердые электролиты, суперионные проводники, электрохимические батареи, суперконденсаторы, твердотельные газовые сенсоры.
Ионные проводники представляют интерес при создании суперконденсаторов [2–4], топливных элементов, электрохимических батарей [5–8]. Выбор материалов катодов электрохимических батарей для каждого конкретного применения нетривиален — он определяется балансом удельной запасаемой энергии, удельной мощности, срока службы, безопасности, стоимости [9]. К примеру, несмотря на то, что такой известный материал как LiFePO4 обладает приемлемыми электрохимическими характеристиками, безопасностью, термической стабильностью, у него имеются такие недостатки как меньшая в сравнении с другими аналогами удельная энергия и снижение характеристик при высокой и низкой температурах. LiNiCoAlO2, напротив, несмотря на то, что менее широко используется на потребительском рынке, обладает высокой удельной запасаемой энергией, удельной мощностью и длительным сроком службы, чем привлекает внимание производителей автомобильной продукции. Именно LiNiCoAlO2 используется в электромобилях Tesla. Требования, предъявляемые к катодным материалам, разнообразны, не удовлетворены в полной степени, а некоторые важные свойства известных материалов не понятны, потому поисковые исследования в данном направлении продолжаются. Наноматериалы позволяют раскрыть высокую производительность материалов, как это было в случае с LiFePO4 [10–15]. Также имеются компьютерные ресурсы, помогающие инженерам и ученым в поиске новых материалов, например AFLOWLIB.org консорциум [16] и Materials Project [17]. В качестве дополнительного инструмента, позволяющего определять возможность существования ионной проводимости, может рассматриваться метод кристаллохимического анализа с использованием программного пакета TOPOS [18]. Данный подход позволяет доступными компьютерными средствами оценить возможность существования ионной проводимости в исследуемых соединениях и структурах. Определение количественных характеристик, влияющих на величину ионной проводимости возможно с использованием более точных подходов метода теории функционала плотности. Первичный анализ данным методом, тем не менее, позволяет исключить структуры, в которых не может наблюдаться ионная проводимость, сэкономить имеющиеся у исследователя компьютерные ресурсы. Наличие в ресурсе Materials Project информации о величине запрещенной зоны, положению уровня Ферми для интересующего материала позволяет также, не прибегая к ресурсоемким расчетам в первом приближении оценить порядок электронной проводимости. Высокая электронная проводимость необходима для свободного перемещения электронов в область протекания электрохимической реакции, а также быстрого установления равновесия между данными областями, отсутствие которого может отразиться на величине саморазряда [19]. Данные рассуждения верны для материалов катодов электрохимических батарей и электродов суперконденсаторов с псевдоемкостным накоплением заряда [20]. Следует отметить, что твердых электролитов твердотельных электрохимических батарей величина электронной проводимости, напротив, должна быть минимальной.
Тем не менее, имеются исключения из правил: LiFePO4 несмотря на невысокие электронную и ионную проводимости (для монокристаллического материала) демонстрирует высокие скорости заряда-разряда и длительный срок службы (в виде нанопорошков, при допировании материала и покрытии их углеродом). Предполагается, что данный факт можно объяснить существованием фазы метастабильного твердого раствора на границе LiFePO4/FePO4 [21].
Данный метод анализа не является достаточным — некоторые из катионов могут быть неподвижными, формировать вместе с анионами многогранники и не участвовать в ионной проводимости. Такие многогранники способны соединяться с другими и образовывать устойчивый каркас как, например, это имеет место в случае с Li10GeP2S12 [22]. Возможно улучшение данной методики — определение неподвижных катионов Li и исключение позиций, занимаемых ими из определения путей миграции. Данный способ мог бы быть реализован, например, построением доменов Дирихле вокруг каждого катиона Li и определение неподвижных на основе измерения телесных углов данных доменов. Мы предполагаем продемонстрировать перспективность данного подхода на примере соединения Li10GeP2S12 в дальнейших публикациях.
Анализируя результаты расчетов (табл. 1) видно, что с увеличением размерности путей миграции катионов число соединений, в которых может наблюдаться ионная проводимость при комнатной температуре, увеличивается (43 для 1D, 55 для 2D, 94 для 3D). Для пространственной группы P1, напротив, уменьшается, (14 для 1D, 16 для 2D, 7 для 3D): вероятно с уменьшением симметрии кристаллической решетки наблюдение 3D ионной проводимости все менее вероятно. Однако следует отметить, что пример наличия высокой ионной проводимости у 2D ионного проводника LiCoO2 доказывает, что размерность каналов не является определяющей для высокой ионной проводимости. Только использование более точных (и ресурсоемких) компьютерных вычислений методами теории функционала электронной плотности позволяет определять энергии активации ионной проводимости и в данном случае метод кристаллохимического анализа также не является достаточным.
Для соединений, в которых возможна ионная проводимость, помимо замещения среди переходных и редкоземельных металлов, распространены случаи замещения P5+ на V5+ (3 случая — см. табл. 1). Данные соединения при замещении приобретают различную электронную структуру согласно информации ресурса Materials Project: ширина запрещенной зоны, как и для других случаев, различается существенно — 3,279 и 2,395 эВ для Na3P(SO)2 и Na3V(SO)2 соответственно. На примере данных соединений допирование позволило бы понижать ширину запрещенной зоны и повышать электронную проводимость, необходимую для материалов катодов электрохимических батарей. Использование допирования для повышения электронной проводимости используется для такого известного материала как LiFePO4. Как и в случае с Li-ионными проводниками, метод кристаллохимического анализа с использованием программного пакета TOPOS позволяет производить экспериментальный поиск новых ионных проводников путем гетеровалентного замещения катионов каркаса [23]. Термодинамически нестабильные ионные проводники при этом можно получить, например, методом ионного обмена термодинамически стабильных.
Таблица 1
Размерность путей миграции ионов Na в неорганических Na-содержащих материалах
Состав |
Число соединений |
Пространственная группа |
1D путь миграции |
||
Na5WO5 |
1 |
C2/c |
Na2Co2(B4O7)3 |
1 |
C2/c |
Na2Tl2B10O17 |
1 |
C2/c |
Na2FeO3 |
1 |
Ccmm |
Na4TiO4 |
1 |
Cmcm |
Na4CoAuO5 |
1 |
Cmcm |
NaCr4O8 |
1 |
I4/m |
Na2DyMoPO8 M = Y, Tb, Dy, Yb |
8 |
Ibca |
Na2MO2 M = Ni, Pt |
2 |
Immm |
Na25Zr11Sc5Si12(PO8)12 |
1 |
P1 |
Na3Zr2Si2PO12 |
2 |
P1 |
Na10Zr6Si7(PO18)2 |
2 |
P1 |
Na20Zr11Si10(PO9)8 |
2 |
P1 |
Na7Bi3P12(Pb5O24)2 |
1 |
P1 |
Na17Zr12Si11P7O72 |
2 |
P1 |
Na4Al2Si6ClO24 |
2 |
P1 |
Na2Zn3(SiO4)2 |
1 |
P1 |
Na10Zr6Si7(PO18)2 |
1 |
P1 |
Na6Si8O19 |
1 |
P21/c |
Na2Pr4Cl9O2 |
1 |
P21/m |
Na2V3O7 |
1 |
P31c |
Na3TmSi2O7 |
1 |
P63 |
NaP3(PbO3)4 |
1 |
P63 |
NaSb3(PO5)2 |
1 |
Pbc21 |
NaASbO5 A = Si, Ge |
2 |
Pc21n |
Na5AsO5 |
1 |
Pcmb |
Na5Lu(SiO3)4 |
1 |
R32 |
Na2V3O7 |
1 |
|
Na4Ti5O12 |
1 |
|
2D путь миграции |
||
Na5Bi2As(CO4)4 |
1 |
C2/c |
Ba3NaBiO6 |
1 |
C2/c |
Na3Sc2(PO4)3 |
1 |
Cc |
NaCo2H3(SeO5)2 |
1 |
Cm |
Na3MnO3 |
1 |
Cmce |
Na5MoNO4 |
1 |
Cmcm |
Na5Bi2P(CO4)4 |
1 |
F222 |
Na4Bi2C4SO16 |
1 |
F222 |
Na4Bi2As(CO4)4 |
1 |
Fddd |
Na2Ti2Sb2O |
1 |
I4/mmm |
Na3Zr2Si2PO12 |
10 |
P1 |
Na11Zr8Si7P5O48 |
2 |
P1 |
Na4Ga3Si3NO14 |
2 |
P1 |
Na13Zr8Si9(PO16)3 |
2 |
P1 |
Na4Al3Ge3NO15 |
2 |
P2 |
NaBiPCO7 M = Zr, Ce, Pr, Nd, Bi |
5 |
P21 |
NaBiACO7 A = S, As |
2 |
P21 |
Na8M2O7 M = Fe, Ga |
2 |
P21/c |
Na3FeO3 |
1 |
P21/c |
Na2FeO2 |
1 |
P21/c |
Na2Ti2O5 |
1 |
P21/c |
Na3B5SO8 |
1 |
P21/c |
Na3NbAs2O9 |
1 |
P21/c |
Na2CoGeO4 |
1 |
P21nm |
Na3Fe2P2O8F3 |
2 |
P42/mbc |
Na2ZnP2O7 |
1 |
P42/mnm |
NaHg2IO2 |
1 |
P6222 |
NaAsO2 |
1 |
Pcab |
Na3V2P2O8F3 |
2 |
Pnnm |
NaAl11O17 |
1 |
R3m |
Na2SO4 |
1 |
|
Na2CrO4 |
1 |
|
NaTl3(SO3)2 |
1 |
|
NaNiAsO4 |
1 |
|
3D путь миграции |
||
Na2CdO2 |
1 |
C2 |
Na7(CoO3)2 |
1 |
C2/c |
Na2Sb4O7 |
1 |
C2/c |
NaW2(Br3O)2 |
1 |
C2/m |
Na4FeO3 |
1 |
Cc |
Na4BeAlSi4ClO12 |
1 |
Cc |
Na8FeO6 |
1 |
Cc |
Na3AS3O A = P, V |
2 |
Ccm21 |
Na22In5O15 |
1 |
Cm |
Na5WNO4 |
1 |
Cmc21 |
Na2Bi2C4SO16 |
1 |
Fddd |
Na2Sn2O3 |
1 |
I213 |
Na5M(MoO4)4 M = Y, La, Gd, Tb, Er, Lu |
6 |
I41/a |
Na5M(WO4)4 M = Y, La, Nd, Tb, Lu |
5 |
I41/a |
Na4Th(MoO4)4 M = Th, U, Pu |
3 |
I41/a |
Na2Sc4Zn2(MoO4)9 |
2 |
P1 |
Na2Li5Al7Si29O72 |
1 |
P1 |
Na13Zr7Si5P7O48 |
2 |
P1 |
Na4CoO4 |
1 |
P1 |
Na8Be3Si9(ClO12)2 |
1 |
P1 |
Na2Pr4Br9NO |
2 |
P21 |
Na14Al4O13 |
1 |
P21/c |
Na2Al2S3(Cl4O3)2 |
2 |
P21/c |
Na3VO3 |
1 |
P21/c |
Na2TiS2O |
1 |
P21/c |
Na3Si2PO8 |
1 |
P212121 |
Na2MoSeO6 |
1 |
P213 |
Na3MoNO3 |
1 |
P21nm |
Na8Al6Si6SO28 |
1 |
P23 |
Na5MSO2 M = Co, Ni |
2 |
P4/mmm |
Na3Nd2(NO3)9 |
1 |
P4132 |
Na3NiO2 |
1 |
P42/mnm |
Na6MO4 M = Co, Ni |
2 |
P42/nmc |
NaAl3Tl2(SiO4)3 |
2 |
P4mm |
Na3LuSi2O7 |
1 |
P63/m |
Na3IO |
1 |
P63/mmc |
Na8AO6 A = Ti, Co, Sn |
3 |
P63cm |
Na6CoO4 |
1 |
P63mc |
Na5FeO4 |
1 |
Pbca |
Na9Fe2O7 |
1 |
Pca21 |
Na5FeO4 A = Fe, Co, Ni, Ga, Tl |
5 |
Pcab |
Na3A(SO)2 A = P, V |
2 |
Pcab |
Na3TlO2 |
1 |
Pcmn |
Na8TiO6 |
1 |
Pna21 |
Na3VS3O |
1 |
Pnma |
Na7SbO6 |
1 |
R3 |
Na3ASO3 A = P, V |
2 |
R3c |
NaSn4(PO4)3 |
1 |
R3c |
Na8BeAl4Si7(AO12)2 A = Cl, Br |
2 |
|
Na3VO4 |
1 |
|
Na3BiO3 |
1 |
|
NaSbO3 |
1 |
|
Na14Ni2O9 |
1 |
|
Na2CrO4 |
1 |
|
Na5CrO4 |
1 |
|
Na4M(NO)2 M = Mo, W |
2 |
|
Na3SO4F |
1 |
|
NaLa6OsI12 |
1 |
|
Na3Sc2(BO3)3 |
1 |
|
Na4Al3A3IO12 A = Si, Ge |
4 |
|
Na4Ga3Si3IO12 |
2 |
С использованием кристаллохимического анализа (программный пакет TOPOS) проанализированы 2367 неорганических Na-содержащих соединений на способность обладать ионной проводимостью, определена размерность ион-проводящих каналов. Выявлено, что 192 соединения способны обладать ионной проводимостью при комнатной температуре. Метод позволяет провести топологический анализ структур и ион-проводящих каналов с использованием минимальных требований к компьютерным ресурсам и таким образом произвести подготовку к дальнейшим расчетам из первых принципов.
Литература:
1. Guyomard D., in New Trends in Electrochemical Technology: Energy Storage Systems in Electronics, ed. Osaka T. and Datta M., Gordon and Breach Publishers, Philadelphia, 2000, ch. 9.
2. Коваленко А. С., Шилова О. А., Морозова Л. В., Калинина М. В., Дроздова И. А., Арсентьев М. Ю. Особенности синтеза и исследование нанокристаллической кобальто-никелевой шпинели // Физика и химия стекла, 2014, т. 40, № 1, С. 135–145.
3. Арсентьев М. Ю., Тихонов П. А., Калинина М. В., Цветкова И. Н., Шилова О. А. Cинтез и физико-химические свойства электродных и электролитных нанокомпозитов для суперконденсаторов // Физика и химия стекла, 2012, т. 38, № 5, С. 653–664.
4. Шилова О. А., Антипов В. Н., Тихонов П. А., Кручинина И. Ю., Арсентьев М. Ю., Панова Т. И., Морозова Л. В., Московская В. В., Калинина М. В., Цветкова И. Н. Керамические нанокомпозиты на основе оксидов переходных металлов для ионисторов // Физика и химия стекла, 2013, т. 39, № 5, С. 803–815.
5. Тихонов П. А., Попов В. П., Арсентьев М. Ю., Подзорова Л. И., Ильичева А. С., Чернышева И. В., Андреева Н. С. Керамика и тонкие слои на основе трех-компонентной системы ZrO2-CeO2-Al2O3 как перспективные твердые электролиты // Огнеупоры и техническая керамика, 2009, № 6, С. 7–11.
6. Тихонов П. А., Калинина М. В., Арсентьев М. Ю., Пугачев К. Э. Протонопроводящая керамика и тонкие пленки на основе цирконатов La и Sm // Физика и химия стекла, 2012, т. 38, № 4, С. 553–564.
7. Тихонов П. А., Арсентьев М. Ю., Калинина М. В. Наноразмерные пленки на основе диоксидов циркония и церия // Физика и химия стекла, 2010, т. 36, № 2, С. 289–296.
8. Арсентьев М. Ю., Калинина М. В., Тихонов П. А., Морозова Л. В., Коваленко А. С., Ковалько Н. Ю., Хламов И. И., Шилова О. А. Синтез и свойства сенсорных оксидных наноразмерных пленок в системе ZrO 2-CeO 2 // Физика и химия стекла, 2014, т. 40, № 3, С. 478–484.
9. http://batteryuniversity.com/learn/article/types_of_lithium_ion
10. Грачева И. Е., Мошников В. А., Абрашова Е. В. Обобщение результатов анализа величины фрактальной размерности золь-гель пористых иерархических структур // Материаловедение, 2013, № 6, с. 13–22.
11. Abrashova E. V., Gracheva I. E., Moshnikov V. A. Functional nanomaterials based on metal oxides with hierarchical structure // Journal of Physics: Conference Series, 2013, Vol. 461, conference 1, P. 012019.
12. Abrashova E. V., Gracheva I.E, Moshnikov V. A. Metal oxide SnO2–ZnO–SiO2 films prepared by sol-gel // Smart Nanocomposites, 2014, Vol. 4, Iss. 2, pp. 1–7.
13. Абрашова Е. В., Барановский М. В. Получение и анализ спектральных характеристик нанокомпозитов на основе широкозонных проводящих металлоксидов системы ZnO-SnO2-SiO2 // Известия СПбТЭТУ «ЛЭТИ», 2013, т. 5, с. 16–21.
14. Шмигель А. В., Тихонов П. А., Арсентьев М. Ю., Дроздова И. А. Получение и исследование наночастиц металлического серебра методом анодного окисления // Физико-химические аспекты изучения кластеров, наноструктур и наноматериалов, 2013, № 5, с. 390–396.
15. Шмигель А. В., Тихонов П. А., Арсентьев М. Ю., Лапшин А. Е. Электрохимическое получение наноструктурированных материалов на основе никеля // Физико-химические аспекты изучения кластеров, наноструктур и наноматериалов, 2013, № 5, с. 385–389.
16. Curtarolo S. et al. AFLOWLIB.ORG: A distributed materials properties repository from high-throughput ab initio calculations // Comp. Mater. Sci., 2012, Vol. 58, pp. 227–235.
17. Jain A. et al. A high-throughput infrastructure for density functional theory calculations // Comp. Mater. Sci., 2011, Vol. 50, Iss. 8, pp. 2295–2310.
18. Blatov V. A. Multipurpose crystallochemical analysis with the program package TOPOS // IUCr CompComm Newsletter, 2006, pp. 4–38., Anurova N. A., Blatov V. A., Ilyushin G. D., Blatova O. A., Ivanov-Schitz A. K., Dem’yanets L. N. Migration maps of Li+ cations in oxygen-containing compounds // Solid State Ionics, 2008, Vol. 179, Iss. 39, pp. 2248–2254.
19. Liu Tongchang, Pell W. G., Conway B. E. Self-discharge and potential recovery phenomena at thermally and electrochemically prepared RuO2 supercapacitor electrodes // Electrochimica Acta, 1997, Vol. 42, Iss. 23–24, pp. 3541–3552.
20. Chung S.-Y., Bloking J. T., and Chiang Y.-M. Electronically Conductive Phospho-Olivines as Lithium Storage Electrodes // Nature Materials, 2002, Vol. 1, Iss. 2, pp. 123–128.
21. http://www.greencarcongress.com/2014/06/20140609-mit.html
22. Kayama N., Homma K., Yamakawa Y., Kanno R., Yonemura M., Kamiyama T., Kato Y., Hama S., Kawamoto K., Mitsui A. A lithium superionic conductor // Nature Materials, 2011, Vol. 31, Iss. 10, pp. 682–686.
23. Anurova N. A., Blatov V. A., Ilyushin G. D., Blatova O. A., Ivanov-Schitz A. K., Dem'yanets L. N. Migration maps of Li+ cations in oxygen-containing compounds // Solid State Ionics, 2008, Vol. 179, Iss. 39, pp. 2248–2254.