INTERRELATION OF COMPOSITION, STRUCTURE AND CHEMICAL PROPERTIES OF REFRACTORY INTERSTITUAL PHASES. PART I. REGULARITIES OF REACTIONS OF REFRACTORY INTERSTITIAL PHASES WITH SOLID REAGENTS

Abstract


The interrelation of composition, structure and chemical properties cubic (type NaCl) refractory interstitial phases (RIP) – carbides, nitrides and monoxides transitive metals IV, V groups is investigated. Regularities of reactions of RIP with various solid reagents (transitional metals of IV, V groups, their oxides, carbon) are discussed. It is established that chemical activity of RIP components is function their thermodynamic activity. It is shown that chemical activity of RIP components is function their thermodynamic activity. It is noted that regularities of change of the RIP chemical activity in the field of their homogeneity and in the rank MeC–MeN–MeO are in good agreement with the solid solution nature of these phases.

Full Text

В предыдущей статье [1] нами проанализирована взаимосвязь состава, структуры и физических свойств кубических (типа NaCl) тугоплавких фаз внедрения (ТФВ) – карбидов, нитридов и монооксидов переходных металлов IV, V групп. На многочисленных примерах показано, что непротиворечивая интерпретация различных физических свойств ТФВ (структурных, термодинамических, физико-механических, теплофизических) возможна лишь в том случае, если исходить из твердорастворной природы этого класса веществ, т.е. из представления о решающей роли в них состояния Me–Me-связей, формирующих металлический каркас решетки. Предложено соответствующее модельное представление, качественно верно описывающие твердорастворную природу ТФВ, характер существующих в них Me–X- и Ме–Ме-связей, их относительную прочность и закономерности изменения с составом. С другой стороны, если представление о твердорастворной природе ТФВ (и следовательно, относительной автономности в них подрешеток металла и неметалла) действительно справедливо, мы вправе ожидать сохранения каждым компонентом, входящим в состав ТФВ, своей химической индивидуальности, или, другими словами, проявления ими селективной активности в различных химических процессах. С целью подтверждения этого предположения нами был разработан метод так называемого химического зондирования подрешеток металла и неметалла в ТФВ различными твердыми реагентами [2]. В качестве «зонда» в этом методе используются специально подобранные вещества, способные как избирательно, так и одновременно взаимодействовать с каждым компонентом, входящим в состав ТФВ. Откликом исследуемой ТФВ на конкретное химическое возмущение является специфичный для каждого состава этой фазы продукт реакции (набор продуктов реакции). Необходимая для решения вопросов информация извлекается из анализа функциональных зависимостей между составом исследуемой ТФВ, с одной стороны, температурой начала реакции и составом первичного продукта взаимодействия, с другой. В первой части настоящей статьи излагаются основные результаты изучения начальных стадий взаимодействия порошкообразных ТФВ с различными твердыми реагентами (переходные металлы IV, V и VIII групп, оксиды переходных металлов, углерод), полученные с использованием изложенной методологии. Основное внимание при этом уделялось выяснению следующих дискуссионных вопросов химического материаловедения ТФВ: 1. Какой из возможных процессов – реакционный или диффузионный – реализуется в начальный период взаимодействия ТФВ с реагентом? 2. Чем определяется состав первичного продукта химической реакции в смесях ТФВ – реагент при наличии альтернативных возможностей? 3. Существует ли корреляция между химической и термодинамической активностью компонентов, входящих в состав ТФВ? Результаты исследования здесь и далее представлены в виде схем протекания соответствующих химических процессов. Над формульным обозначением ТФВ указан параметр ее кубической решетки. Реакционное взаимодействие (собственно химическая реакция) обозначается стрелкой ®, а диффузионное – стрелкой . В качестве основных аналитических методов в этом исследовании были использованы химический анализ, рентгенофазовый анализ (дифрактометры ДРОН-2,0 и Stadi P), термический анализ (ДТА-ТГ анализатор Thermoflex), сканирующая электронная микроскопия и рентгеноспектральный микроанализ (установка JCXA-733). Системы ТФВ – переходный металл. Начальные стадии взаимодействия карбидов с переходными металлами IV, V групп, которые характеризуются высоким сродством к элементам внедрения, изучены методами ДТА-ТГ, РФА и химического анализа. Методические особенности экспериментов заключались в следующем. Брикетированные смеси исходных компонентов (dср = 1…3 мкм), взятых в эквимолярном соотношении, нагревались со скоростью 20 град/мин в токе высокочистого аргона до температуры начала реакции, фиксируемой по отклонению кривой ДТА от нулевой линии. В качестве эталона сравнения использовалась сама исследуемая ТФВ, что позволяло определить начало реакции по изменению температуропроводности смеси. После 5-минутной выдержки при температуре начала реакции образцы охлаждались со скоростью 200 град/мин. Продукты реакции анализировались методом РФА. Если количество образующейся фазы и/или степень ее гомогенности были недостаточными для ее надежной идентификации, то после тщательного перетирания и брикетирования смеси процедура повторялась. Характерные примеры изученных реакций и рентгенограмм, на основании которых были составлены схемы их протекания, приведены ниже (реакции (1) и (2), рис. 1 и 2). Контрольные реакции свободного углерода с теми же переходными металлами (помечены звездочкой) показаны для сравнения. Подпись: 4,68 Å 4,469 Å 4,68 Å NbC0,97 + Zr → NbC0,97–x + ZrC~0,6, (1) * C + Zr → ZrC~0,6. ZrC0,97 + Nb → ZrC0,97–x + Nb2C, (2) * C + Nb → Nb2C. Аналогичные химические реакции протекают в начальный период взаимодействия и в смесях карбидов титана и ванадия с одноименными и разноименными металлами (реакции (3)–(6)), а также в контрольных смесях этих металлов со свободным углеродом: VC0,87 + Ti → VC0,87–x + TiC~0,6, (3) 4,165 Å 4,31 Å TiC0,96 + Ti → TiC0,96–x + TiC~0,6, (4) * C + Ti → TiC~0,6. TiC0,96 + V → TiC0,96–x + V2C, (5) VC0,87 + V → VC0,87–x + V2C, ( 4,326 Å 4,31 Å 6) * C + V → V2C. а б Рис. 1. Рентгенограммы смеси NbC0,97 + Zr: а – до реакции; б – после реакции (750 °С, 3´5 мин) а б Рис. 2. Рентгенограммы смеси ZrC0,97 + Nb: а – до реакции; б – после реакции (920 °С, 2´5 мин) Кратко прокомментируем приведенные результаты. Отметим, во-первых, что начальной стадией взаимодействия карбидов с переходными металлами IV, V групп является химическая реакция углерода карбидной фазы с металлом (непосредственное диффузионное взаимодействие реагентов на начальном этапе процесса не реализуется). Этот результат представляется вполне закономерным, так как энергия активации твердофазной химической реакции должна быть сопоставима с энергией активации поверхностной диффузии атомов реагирующих компонентов (T ≈ 0,3Тпл), которая в полтора-два раза меньше объемной (T ≥ 0,5Тпл) [3, 4]. Бытующее в литературе мнение, что энергия активации твердофазной химической реакции обычно больше, чем объемной диффузии (см., например, [5, 6]), является, по нашему мнению, следствием неоправданного распространения закономерностей, установленных для газофазных реакций (где это положение справедливо), на твердофазные процессы. Как было показано еще В.И. Архаровым [7], аррениусовский формализм, разработанный для оценки энергии активации реакций в газовой фазе, использовать для изучения твердофазных химических реакций в общем случае некорректно. Из приведенных результатов также следует, что первичным продуктом взаимодействия в рассматриваемых системах является карбидная фаза, термодинамически наименее стабильная в данных условиях (фаза с минимальным содержанием углерода). Причем для металлов IV группы это кубический карбид, близкий по составу к нижней границе области гомогенности, а в случае металлов V группы – гексагональный карбид Me2C. Эти результаты тоже не являются неожиданными. Они хорошо согласуются с известным правилом последовательности превращений Оствальда, уточненным позднее Н.В. Танцовым [8, 9] (принцип минимального изменения энтропии). Особый интерес представляет тот факт, что начальные стадии взаимодействия свободного углерода и углерода карбидной фазы MeC~1,0 с одним и тем же переходным металлом протекают химически однотипно, т.е. сопровождаются образованием карбидов практически одного и того же состава. Это обстоятельство свидетельствует о слабой связи углерода с решеткой близких к стехиометрии карбидов и хорошо согласуется с его высокой термодинамической активностью [10–13]. Наконец, важным результатом исследования является впервые обнаруженная корреляция между температурой начала взаимодействия ТФВ с переходными металлами IV,V групп и степенью ионности Me–X-связей, оцениваемой по энтальпии образования ТФВ (рис. 3–6). Рис. 3. Зависимость температуры начала реакции TiC, TiN и TiO с титаном от энтальпии образования ТФВ Рис. 4. Зависимость температуры начала реакции VC, VN и VO с ванадием от энтальпии образования ТФВ Рис. 5. Зависимость температуры начала реакции TiC, TiN и TiO с цирконием от энтальпии образования ТФВ Рис. 6. Зависимость температуры начала реакции VC, VN и VO с цирконием от энтальпии образования ТФВ Наличие подобных корреляций свидетельствует о том, что Me–X-связи в ТФВ являются преимущественно ионными. Согласуются с таким заключением и относительно высокие стандартные энтальпии образования ТФВ. Например, для TiC, TiN и TiO они равны соответственно –46, –81 и –124 ккал/моль [14], что сопоставимо со стандартной энтальпией образования изоструктурного им NaCl (–98 ккал/моль [15]). Ввиду того, что этот вывод находится в противоречии с общепринятыми сегодня теоретическими представлениями о характере Me–X-связей в ТФВ, рассмотрим суть этого противоречия подробнее. Согласно результатам современных квантово-химических расчетов и основанной на них интерпретации рентгеноспектральных данных (см., например, источники [16–20]) высокие прочностные свойства ТФВ обусловлены наличием в них сильных, преимущественно ковалентных Me–X-связей. На этом основании делается вывод о том, что причиной снижения температур плавления, твердости и прочности ТФВ в ряду MeC–MeN–MeO, где Me = Ti или V, является уменьшение в этом направлении степени ковалентности Me–X-связей. Приведенные на рис. 3–6 экспериментальные результаты не согласуются с таким представлением. Действительно, в рядах TiC–TiN–TiO и VC–VN–VO температуры начала взаимодействия ТФВ с одноименными (Ti, V) и с разноименными (Zr) переходными металлами не снижаются, как следовало бы ожидать, если исходить из ковалентной природы Me–X-связей, а напротив, повышаются. При этом наблюдается четкая корреляция с энтальпией образования ТФВ: чем она больше по абсолютной величине (чем выше степень ионности Me–X-связей [15]), тем выше температура начала взаимодействия ТФВ с металлом. Это означает, что Me–X-связи в близких к стехиометрии ТФВ имеют преимущественно ионный, а не ковалентный характер. Противоречие снимается, и наблюдаемые корреляции получают разумное объяснение только в том случае, если исходить из представления о твердорастворной природе ТФВ и конкурентном характере существующих в них Me–Me- и Me–X-связей. Тогда естественно, что усиление Me–X-связей в ряду MeC–MeN–MeO, т.е. их сокращение, обусловливает вынужденное сближение соседних атомов металлов, что сопровождается ростом кинетической энергии электронов Me–Me-связей и их делокализацией. Следствием этого эффекта является рост плотности Me3d-состояний вблизи уровня Ферми и снижение устойчивости фазы (в соответствии с критерием J. Yamashita [21]) вплоть до появления вакансий в подрешетке металла ТФВ, снимающих избыточную энергию системы. Итоговый результат – деградация прочностных свойств ТФВ. Твердорастворная природа ТФВ проявляется и при их взаимодействии с переходными металлами VIII группы, которые характеризуются низким сродством к элементам внедрения. Покажем это на примере систем MeCx–Ni. Особенности проведения экспериментов заключались в следующем. Брикетированные смеси исходных компонентов крупностью 50–63 мкм, содержащие 20, 50 и 80 % никеля, отжигались в среде особо чистого аргона при 1150 °С в течение 5–100 ч и анализировались методами РФА и РСМА. Для уточнения промежуточных составов тугоплавкой и металлической фаз проводилась процедура фазового разделения реагентов (металлическая фаза образца отделялась от тугоплавкой фазы путем ее растворения в кипящей соляной кислоте) с последующим химическим анализом ТФВ и раствора. На рис. 7 и 8 в качестве примера приведены зависимости «состав ТФВ – продолжительность отжига», полученные для смесей MeCx+80%Ni. Видно, что процесс растворения карбидов в никеле сопровождается изменением их химического состава. Установлено, что системы MeCx (x→1)–Ni не являются квазибинарными. Об этом свидетельствует появление в образцах при достижении ими равновесного состояния от 0,4 до 0,9 % свободного углерода, что в несколько раз превышает его содержание в исходных смесях. Неквазибинарность систем MeCx (x→1)–Ni и инконгруэнтный характер растворения ТФВ в никеле являются, очевидно, следствием твердорастворной природы ТФВ. Таким образом, что механизм взаимодействия карбидов с переходными металлами IV, V и VIII групп тесно связан с термодинамической активностью (химическим потенциалом) компонентов, входящих в их состав. Действительно, если активность углерода в карбиде приближается к единице (MeCx при x → 1), его связь с решеткой ТФВ ослабляется. В результате углерод карбидной фазы либо вступает с переходным металлом в химическую реакцию с образованием соответствующего карбида (переходные металлы IV, V групп), либо преимущественно растворяется в нем (переходные металлы VIII группы). Рис. 7. Кинетика изменения состава ТФВ в смесях TiCx + Ni и ZrCx + Ni Рис. 8. Кинетика изменения состава ТФВ в смесях VCx + Ni и NbCx + Ni С уменьшением x в MeCx термодинамическая активность атомов углерода снижается, а металла – соответственно растет [10–13]. В карбидах с минимальным содержанием углерода активность металлического компонента ТФВ уже на несколько порядков (от трех в VCx до девяти в ZrCx) превышает соответствующую величину для углерода. Следствием повышенной активности металлического компонента ТФВ является его преимущественное растворение в никеле, что сопровождается изменением состава растворяющейся карбидной фазы (увеличением содержания в ней углерода). Системы ТФВ – оксид переходного металла. Селективный характер проявления химической активности компонентами, входящими в состав ТФВ, взаимосвязь химической и термодинамической активности компонентов четко проявляются и при зондировании подрешеток металла и неметалла ТФВ кислородом различных оксидных фаз. Два примера таких реакций приведены ниже. Рассмотрим сначала особенности начальных стадий взаимодействия карбида титана с кислородом оксида никеля (реакции (7)–(9)): TiC0,96 + NiO → TiC0,96–x + CO↑ + Ni, (7) TiC0,80 + NiO → TiO2 + CO↑ + Ni, (8) TiC0,52 + NiO → Ti2O3 + C + Ni. (9) Из приведенных результатов видно, что если состав карбида титана близок к стехиометрии, то в реакцию с кислородом оксида никеля вступает преимущественно углерод (несмотря на то, что химическое сродство титана к кислороду намного выше, чем углерода), а если далек – то титан (с выделением углерода в свободном виде). В случае промежуточного состава карбидной фазы, реакция (8), оба ее компонента одновременно участвуют в акте восстановления оксида никеля. Соответственно изменяется состав продуктов взаимодействия. Особый интерес, как с теоретической, так и практической точки зрения (карботермическое восстановление диоксида титана) представляет реакция карбида титана с одноименным монооксидом. Дело в том, что изоструктурные карбид и монооксид титана характеризуются полной взаимной растворимостью [22]. Традиционно считается, что твердые растворы между ними формируются путем взаимной диффузии компонентов без участия каких-либо промежуточных химических реакций. Между тем такое участие не только нельзя исключать, но следует предполагать, поскольку, как мы убедились, подрешетки металла и неметалла в карбидах относительно автономны, а углерод и кислород характеризуются высоким химическим сродством. С целью проверки этого предположения мы провели подробное исследование механизма образования непрерывного твердого раствора в эквимолярной смеси TiC0,96 + TiO1,03. Оказалось, что, действительно, диффузионной стадии взаимодействия TiC0,96 и TiO1,03 предшествует химическая реакция углерода карбидной фазы с кислородом оксидной фазы, сопровождающаяся диспропорционированием последней на Ti2O3 и металлический титан. При этом фаза Ti2O3 рентгенографически фиксируется в продуктах взаимодействия уже после первичной термообработки смеси (1100 °C, 5 мин), а фаза на основе металлического титана – только после дополнительной двукратной обработки той же смеси (рис. 9). В результате протекания этого процесса существенно (от 0,77 до 1,08 м2/г) увеличивается удельная поверхность смеси. Дальнейшее изучение механизма образования твердого раствора в рассматриваемой системе проводилось на образцах, полученных термообработкой в вакуумной печи в изотермических условиях (1200–1700 °C). Для этого брикетированные смеси исходных компонентов отжигались по ступенчатому (через 100 °C) режиму в течение 1–5 ч. Составы промежуточных фаз определялись методом РФА и уточнялись путем моделирования каждой стадии контрольными реакциями (10)–(13). а б Рис. 9. Рентгенограммы смеси TiC0,96 + TiO1,03: а – до реакции; б – после реакции (1100 °С, 3´5 мин) Общая схема механизма образования твердого раствора в эквимолярной смеси TiC0,96 и TiO1,03 приведена на рис. 10. Рис. 10. Механизм образования твердого раствора в эквимолярной смеси TiC0,96 и TiO1,03 Контрольные реакции: C + TiO1,03 → CO↑ + Ti2O3 + α-Ti, (10) α-Подпись: 4,326 Å 4,26 Å Подпись: 4,326 Å 4,31 Å Ti + TiC0,96 → TiC0,6 + TiC0,96–x , (11) Ti2O3 + TiC0,96 → TiC0,3O0,7 + CO↑, (12) Подпись: 4,300 Å 4,21 Å Ti2O3 + TiC0,52 → TiC0,1O0,9 + CO↑. (13) Дополнительными экспериментами установлено, что непосредственное диффузионное взаимодействие TiCx и TiOy наблюдается лишь в том случае, если составы обеих ТФВ близки к нижней границе области гомогенности, т.е. когда активность элементов внедрения в них недостаточна для протекания реакции [C] + [O] → CO↑. Взаимное растворение TiCx и TiOy протекает в этом случае без заметного изменения массы смеси и сопровождается монотонным изменением параметров решетки обеих фаз (реакция (14), рис. 11). 4,314 Å 4,192 Å 4,256 Å TiC0,65 + TiO0,78 TiC0,65–xOy + TiO0,78–yCx TiC0,32O0,39. (14) Рис. 11. Кинетика образования твердого раствора в эквимолярной смеси TiC0,65 + TiO0,78 Системы ТФВ – углерод. Корреляция между химической и термодинамической активностью компонентов ТФВ обнаруживается и при зондировании подрешетки металла ТФВ углеродом. Показательна в этом плане реакция порошкообразного (dср = 1…3 мкм) нитрида титана с углеродом (Sуд = = 19,0 м2/г). Методом ДТА-ТГ установлено, что с уменьшением z в TiNz (т.е. с ростом активности титана) температура начала реакции [Ti] + C монотонно снижается (рис. 12). Как и следовало ожидать, диффузия углерода в решетку нитрида титана на начальной стадии процесса не реализуется. По данным РФА первичным продуктом взаимодействия в смесях TiNz+C, как и в контрольной смеси Ti+C, является дефектный карбид титана TiC~0,6 (а = 4,31 Å). Азот практически не принимает участия в этой реакции и выделяется в газовую фазу. Подпись: Рис. 12. Влияние состава нитрида титана на температуру начала реакции с углеродом Сделаем следующие выводы. Исследована взаимосвязь состава, структуры и химических свойств кубических (типа NaCl) тугоплавких фаз внедрения (ТФВ) – карбидов, нитридов и монооксидов переходных металлов IV, V групп. На примере твердофазных реакций ТФВ с переходными металлами IV, V, VIII групп, их оксидами и углеродом продемонстрирован селективный характер проявления химической активности компонентами, входящими в их состав. Установлено, что химическая активность компонентов ТФВ является функцией их термодинамической активности. Показано, что закономерности изменения химической активности ТФВ в области их гомогенности и в ряду MeC–MeN–MeO хорошо согласуются с твердорастворной природой этих фаз.

About the authors

Victor Aleksandrovich Zhilyaev

Institute of Solid State Chemistry Ural Branch of RAS, Yekaterinburg

Email: zhilyaev@ihim.uran.ru
620990, Yekaterinburg, Pervomayskaya st., 91 Doctor of Technical Sciences, leading research associate, Institute of Solid State Chemistry, Ural Branch of RAS

References

  1. Жиляев В.А. Взаимосвязь состава, структуры и физических свойств тугоплавких фаз внедрения // Вестник ПГТУ. Машиностроение, материаловедение. – Пермь, 2011. – Т. 13, № 3. – С. 106–116.
  2. Жиляев В.А., Могилевская И.C., Аскарова Л.Х. Механизм фазовых превращений карбида титана в условиях дозированного окисления // Известия СО АН СССР, Серия хим. наук. – 1982. – Вып. 6. – С. 53–57.
  3. Гегузин Я.Е., Когановский Ю.С. Диффузионные процессы на поверхности кристалла. – М.: Энергоатомиздат, 1984. – 124 с.
  4. Каур И., Густ В. Диффузия по границам зерен и фаз. – М.: Машиностроение, 1991. – 448 с.
  5. Дыбков В.И. Кинетика твердофазных химических реакций. – Киев: Наукова думка, 1992. – 128 с.
  6. Физическая химия. – М.: Высшая школа, 2003. – 527 с.
  7. Архаров В.И. Об уточнении понятий «энергия активации» и «элементарный акт» для твердого состояния вещества // Журн. техн. физики. – 1954. – Т. 24, № 3. – С. 375–387.
  8. Танцов Н.В. О направлении самопроизвольной кристаллизации и химических превращений // Журн. рос. физ.-хим. общества. – 1923. – Т. 55, Вып. 5/9. – С. 342–367.
  9. Танцов Н.В. Физико-химический процесс как совокупность элементарных превращений, протекающих адиабатным путем // Журн. рос. физ.-хим. общества. – 1926. – Т. 58, Вып. 8. – С. 947–956.
  10. Малкин В.И., Покидышев В.С. Термодинамические свойства карбида ванадия // Изв. АН СССР, Металлы. – 1969. – № 2. – С. 183–187.
  11. Storms E., Calkin B., Yenche A. The Vaporization Behavior of the Defect Carbides. Part I. The Nb–C System // High Temp. Sci. – 1969. – Vol. 1. – P. 430–455.
  12. Малкин В.И., Покидышев В.С. Изучение термодинамических свойств карбида титана в области гомогенности методом э.д.с. // Журн. физ. химии. – 1971. – Т. 45, № 8. – С. 2044–2046.
  13. Guillermet A.F. Analysis of thermochemical properties and phase stability in the zirconium – carbon system // J. Alloys and Compounds. – 1995. – Vol. 217. – P. 69–89.
  14. Свойства, получение и применение тугоплавких соединений. – М.: Металлургия, 1986. – 928 с.
  15. Голутвин Ю.М. Теплоты образования и типы химической связи в неорганических кристаллах. – М.: Изд-во АН СССР, 1962. – 96 с.
  16. Ивановский А.Л., Жуков В.П., Губанов В.А. Электронное строение тугоплавких карбидов и нитридов переходных металлов. – М.: Наука, 1990. – 224 с.
  17. Mizuno M., Tanaka I., Adachi H. Chemical bonding in titanium-metalloid compounds // Phys. Rev. B. – 1999. – Vol. 59, № 23. – P. 15033–15047.
  18. A theoretical study on the chemical bonding of 3d-transition metal carbides / Y. Zhang, J. Li, L. Zhou, S. Xiang // Solid State Communication. – 2002. – Vol. 121. – P. 411–416.
  19. First-principles calculations of mechanical properties of TiC and TiN / Y. Yang, H. Lu, C. Yu, J.M. Chen // J. Alloys and Compounds. – 2009. – Vol. 485. – P. 542–547.
  20. Electronic structure and elastic constants of TiCxN1–x, ZrxNb1–x and HfCxN1–x alloys: A first-principles study / W. Weng, S. Cui, H. Hu, G. Zhang, Z. Li // Physica B. – 2011. – Vol. 406. – P. 3631–3635.
  21. Yamashita J., Asano S. Electronic structure of CsCl-type transition metal alloys // Progr. Theor. Phys. B. – 1972. – Vol. 48. – P. 2119–2131.
  22. Алямовский С.И., Зайнулин Ю.Г., Швейкин Г.П. Оксикарбиды и оксинитриды металлов IVa и Va подгрупп. – М.: Наука, 1981. – 144 с.

Statistics

Views

Abstract - 22

PDF (Russian) - 57

Refbacks

  • There are currently no refbacks.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies