Особенности локального атомного окружения и электронной структуры иона меди Cu 2+в водных растворах

Особенности локального атомного окружения и электронной структуры

иона меди Си2в водных растворах

1 12 М.А. Солдатов , М.А. Кременная , А.П. Чайников

1 Южный федеральный университет, Ростов-на-Дону 2Ростовский государственный университет путей сообщения

Аннотация: Исследование окружения иона меди Си2+ в водных растворах было проведено методами рентгеновской спектроскопии поглощения. Были измерены спектры рентгеновского поглощения для различных концентраций водных растворов СиС12. Была проведена отладка методики расчётов околопороговой тонкой структуры спектроврентгеновского поглощенияза ^-краем Си в координированных комплексах иона Си2+ методом конечных разностей для решения уравнения Шредингера.Проведены расчеты для различных электронных конфигураций иона меди Си2+. Предложена модель локального окружения иона меди Си2+в водных растворах.

Ключевые слова: спектроскопия рентгеновского поглощения, ионы 3dметаллов, водные растворы, локальное окружение, метод конченых разностей, координационные комплексы.

Введение

В настоящей статье речь пойдет об исследовании окружения иона меди Си2+ в водных растворах методами рентгеновской спектроскопии поглощения (ХАБ - Х-гауАЬвогрйопЗреСхоБСору). Данная задача была поставлена с целью отладки методики расчетов околопороговой тонкой структуры спектров рентгеновского поглощения (далее ХАКЕБ - X-rayЛbsorptionNearEdgeStructure) за ^-краем Си (8979 эВ для медной фольги) в координационных комплексах иона Си2+. Результаты этой работы, в частности, будут в дальнейшем использованы для структурного ХЛКЕБ-анализа центра связывания Си2+ в в-амилоиде [1,2]. Однако, необходимо отметить, что эта задача несомненно представляет и самостоятельный интерес. В настоящее время хорошо установлено, что для глубокого понимания многих химических свойств ионных растворов необходимо точное знание их координационной структуры. Данной тематике посвящено большое количество работ [3, 4], и для многих ионов металлов структура первой, а в некоторых случаях и второй гидратационных оболочек достаточно хорошо известна. В противоположность этому, геометрия

окружения иона Си2+ в водном растворе не вполне понятна и является предметом современных исследований. Повышенный интерес обусловлен также тем, что медь является одним из наиболее важных переходных металлов в биологических системах.

Перед тем, как начать обзор литературы по структурным исследованиям аквакомплексов иона Си2+ в водных растворах, необходимо сказать несколько общих слов о стереохимии соединений двухвалентной меди вообще. В природе в подавляющем числе случаев встречается три типа координации Си2+: Ян-Теллеровски искаженная октаэдрическая координация с 4 короткими (экваториальными) и 2 длинными (осевыми) междуатомными расстояниями,плоская квадратная координация и тетраэдрическая [5]. Плоское квадратное окружение можно рассматривать как предельный случай искаженного октаэдрического с бесконечно длинными осевыми расстояниями. Тетраэдрические же соединения малочисленны, поэтому основным типом координации Си2+ является искаженная октаэдрическая.

Необходимо отметить, что этот тип окружения является отличительным свойством меди, выделяющим ее в ряду двухвалентных ионов других 3б/-металлов, таких как Со, 7п, N1, обладающих регулярной октаэдрической координацией с шестью одинаковыми расстояниями. В случае же иона меди Си2+ в регулярном октаэдрическом поле дублет электронных состояний должен быть заполнен тремя электронами, что можно сделать двумя способами: 2 электрона на 3dz2и 1 на3dx2_y2, или наоборот. Таким образом, энергетический терм является вырожденным, и в результате эффекта Яна-Теллера происходит искажение окружения, понижающее симметрию до D4h и снимающее это вырождение.

В соответствии с вышесказанным на протяжении последних 50 лет предполагалось, что ионы Си2+ в водных растворах представлены в основном 6-координированными аквакомплексами Си (Н20)62+ с Ян-Теллеровски

искаженной октаэдрической симметрией. Однако, результаты, получаемые с помощью различных экспериментальных методов, такие как рентгеновская дифракции (XRD), дифракция нейтронов (ND) и спектроскопия рентгеновского поглощения (XAS) приводили к различным значениям длин связи. Так, например, расстояние до осевых атомов кислорода варьировалось от 2,1 до 2,6 Á [6], что с химической точки зрения является очень большой неопределенностью. Одной из причин, осложняющих определение структурных параметров, является тот факт, что структура аквакомплексов не статична и может иметь несколько различных конфигураций, в которых анионы растворяемых солей замещают одну или несколько молекул воды из гидратационной оболочки иона. Таким образом, при структурном анализе приходится учитывать все возможные варианты. В некоторых случаях, как например, в водных растворах хлорида меди CuCl2, ион хлора может

вытеснить молекулу воды из осевого узла, образуя хлорокупратный

2+

аквакомплекс Cu(H2Ü)5 сг. При этом расстояниеR(Cu-Clax) близко к расстояниям R(Cl-H2Ü) иR(H2Ü - H2Ü), что существенно затрудняет извлечение структурной информации из данных XRD [7] и ND [8]. Более того, существует возможность уменьшить влияние конфигураций, содержащих анионы солей, исследуя растворы маленькой концентрации. Однако, это также осложняет дифракционные исследования в связи с возрастающим соотношением сигнал/шум и не имеет смысла в тех случаях, когда представляетинтереснекоторая определенная концентрация.

В этом плане спектроскопия XAS является очень перспективным методом, так как она обладает элементной чувствительностью и позволяет исследовать локальную геометрическую вокруг интересующих атомов, а так же и электронную структуру. При этом с использованием синхротронного излучения становится возможным исследовать растворы в широком

И

диапазоне концентраций растворенного вещества, включая сильно разбавленные растворы с молярностью порядка 1 тМ [9].

Экспериментальная часть

Было приготовлено шесть образцов:

1. Водный раствор СиС12^2(Н20) с концентрацией 0,1 М;

2. Водный раствор СиС12^2(Н20) с концентрацией 1,0 М;

3. Водный раствор СиС12^2(Н20) с концентрацией 1,9 М;

4. Раствор СиС12^2(Н20) в ацетоне (С3Н60) с концентрацией ~ 0,1

М;

5. Раствор СиС12^2(Н20) в водном растворе этилового спирта с концентрацией 0,1 М;

6. Порошок 2-водного кристаллогидрата СиС12^2(Н20). Растворы 1-5 были получены путем растворения необходимого

количества 2-водного кристаллогидрата хлорида меди(11) СиС12^2(Н20) в соответствующих растворителях при комнатной температуре и нормальном атмосферном давлении.

Масса растворяемого порошка, необходимого для получения раствора заданной концентрации (молярности), рассчитывалась по формуле (1):

т(СиС12 • 2(Я20)) = = М (СиС12 • 2 (Я20)) • ^раствора • Ц (СиС12 • 2 (Я20)), (1) где М {СиС12 • 2 (Я20)) - концентрация раствора в (моль/л); ^раствора- объем раствора, задаваемый в нашем случае объемом мерного стаканчика, в котором приготовляется раствор (использовались два мерных стаканчика объемом 0,025 и 0,100 л);ц (СиС12 • 2 (Я20)) = 170,5 (—)- молярная масса

растворяемого вещества. В результате получаем т(СиС12 • 2(Я20)) в граммах.

Растворы были приготовлены следующим образом: половина необходимого объема растворителя помещалась в мерную колбу, на аналитических весах отбиралась необходимая масса растворяемого вещества и в мерной колбе смешивалась с растворителем до полного растворения. Далее мерная колба заполнялась растворителем до необходимого объема.

Для образов, описанных выше, были получены экспериментальные XANES спектры за К-краем меди. Измерения проводились на синхротроне национального исследовательского центра "Курчатовский институт" на станции "Структурное материаловедение" (СТМ).

Все спектры были измерены в диапазоне энергий рентгеновского фотона от 8750 эВ (~ 229 эВ до К-края медной фольги) до 9875 эВ (~ 896 эВ за К-краем медной фольги, длина волны фотоэлектрона ~ 15 А-1).

Измерения спектров растворов 2, 3 (водные растворы СиС12^2(Н20) с концентрациями 1,0 и 1,9 М соответственно) и порошка 2-водного кристаллогидрата СиС12^2(Н20) проводились в режиме "на прохождение". Растворы были помещены в узкие кюветы (шириной до 1 мм) заклеенные с двух сторон каптоновым скотчем. Для контроля энергии фотонов, отраженных от кристалла-монохроматора, при каждой регистрации спектра образца измерялся спектр стандарта (медной фольги). При обработке спектров положению СиК-края медной фольги приписывалось табличное значение 8979 эВ, что позволяло получить спектры в абсолютной шкале энергии фотонов. Для каждого образца было измерено по одному спектру (по одному "проходу").

Измерение спектров растворов 1, 4 и 5 (~ 0,1 М растворы СиС12^2(Н20) с различными растворителями: вода, ацетон и водный раствор этилового спирта соответственно) проводилось с помощью регистрации флуоресцентного излучения. Образцы находились в пробирках Эппендорфа. В связи с сильным поглощением рентгеновского излучения, проходящего

через образец, спектры стандартов измерить не удалось. Однако, обработка спектров растворов 2, 3 и образца 6, показала, что сдвиг энергетической шкалы при проведении серии измерений составляет не более 0,2 эВ, что позволяет представить результаты для образцов 1, 4 и 5 так же в абсолютной шкале энергии.

Энергия, эВ

Рис. 1. - Экспериментальные XANES спектры СиК-края водных растворов СиС122(Н20) различной концентрации: 0,1 М (черные линии), 1,0 М (синие

линии) и 1,9 М (красные линии).

На рис.1 представлены XANES спектры для водных растворов СиС12^2(Н20) различных концентраций: наблюдаются небольшие отличия в области "белой линии". Они более отчетливо видны на графике производных спектров (Рис. 2): с увеличением концентрации раствора от 0,1 до 1,0 М наблюдается перераспределение интенсивностей максимумов a и Ь (a -уменьшается, Ь - увеличивается, становится шире и сдвигается в сторону

И

меньших энергий). При концентрации 1,9 М говорить о наличии двух максимумов а и Ь сложно из-за относительно большого шума. Однако, в этом случае наблюдается особенность а' (эту ситуацию можно интерпретировать как переход а ^ а' + уширение и сдвижка в область меньших энергий максимума Ь).

9000 9010 9020 Энергия, эВ

Рис. 2. - Первая производная коэффициента рентгеновского поглощения экспериментальныхХЛКЕБ спектров за£-краемСи водных растворов СиС122(Н20) различной концентрации: 0,1 М (черные линии), 1,0 М (синие

линии) и 1,9 М (красные линии).

Отладка методики расчета спектров XANES и оптимизация неструктурных параметров

Расчеты спектров ХЛКЕБ Си^-края для аквакомплексов меди проводились в программном коде ЕБМЫЕ8[10].

Эффективный одноэлектронный потенциал (2), используемый в FDMNES, состоит из двух частей:

Vir) = VCoulomb{r) + Vxc

Кулоновской7С0иг0тЬ (г)и обменно-корреляционной7хс (г). Кулоновская часть находится из уравнения Пуассона для суперпозиции электронных плотностей изолированных атомов с контролируемыми числами заполнения орбиталей(З):

осс

&Vcouiomb

a,i

Для обменно-корреляционной части используется формулировка Хедина-Ландквиста[11]. Очень важной особенностью программного комплекса FDMNES является возможность проводить расчеты в полном потенциале, т.е. без использования маффин-тин приближения. Эта возможность реализуется в методе конечных разностей. Необходимо отметить, что расчеты методом конечных разностей очень трудоемкие и требуют значительно больше вычислительного времени.

Отладка методики расчета спектров состояла из "оптимизации" чисел заполнения атомных состояний, которые используются для построения электронной плотности (см. формулу (3)). Для этого были проведены расчеты для различных электронных конфигураций иона меди Cu2+. Изолированный атом меди имеет электронную конфигурацию [Ar]3d104s1,

2+ 9

ион Cu - [Ar]3d . Однако формальный заряд иона меди как в аквакомплексе, так и в других ионных соединениях, может сильно отличаться от реального. В случае ионных соединений 3 ^-переходных металлов разработчики программы FDMNES рекомендуют для имитирования переноса заряда использовать число заполнения 3d орбиталей такое же, как при формальном заряде (в нашем случаепСи,3^ = 9), но

сохранить нейтральность атома, разместив соответствующее число электронов на 4s- и 4^-орбитали, имеющие большие радиусы. Были проведенырасчёты следующих конфигураций:

1. [Аг^948°4р° (медь имеет заряд +2);

2. [Аг^94в04р2;

3. [Аг^^Чр1;

4. [Лг^94в24р0;

5. [Лr]3d94s04p0 с использованием процедуры самосогласования.

В качестве модели аквакомплекса использовалась структура, полученная в работе [12] - медь + первая гидратационная оболочка, представляющая Ян-Теллеровски искаженный кислородный октаэдр с

о о

расстояниями — 1,968 А и — 2,28 А и атомами водорода в так

называемой дипольной конфигурации. Молекулы воды ориентированы, как показано на рис.3. Точечная группа аквакомплекса -

2+

Рис. 3. - Модель аквакомплекса иона меди Си . Расчет сечений проводился с учетом электрических дипольных и квадрупольных переходов методом конечных разностей в полном потенциале с обменно-корреляционной частью в виде потенциала Хедина-Ландквиста. Потенциал остовной вакансии на 1^-уровне был учтен искусственным

И

перебросом остовного электрона на первую не полностью заполненную орбиталь (3^). Результаты расчетов представлены на рис. 4.

Энергия,

Рис. 4. - XANES спектры СиК-края аквакомплексов меди: черная линия -экспериментальный спектр водного раствора СиС12^2(Н20) концентрации 0,1 М, красная линия - теоретический спектр для 3d94s04p0 электронной конфигурациимеди, зеленая линия - теоретический спектр для 3d94s04p2 электронной конфигурации меди,синяя линия - теоретический спектр для

3d94s14p1 - электронной конфигурации меди, голубая линия -3d94s24p0электронной конфигурации меди, розовая линия - теоретический

расчет с учетом самосогласования На Рис.4 представлены XANES спектры СиК-края аквакомплексов меди: черная линия - экспериментальный спектр водного раствора СиС12^2(Н20) концентрации 0,1 М, цветные линии - теоретические спектры,

полученные при использовании различных электронных конфигураций иона меди (см. подписи на картинке). Теоретические спектры представлены без учёта экспериментального уширения, уширения внутреннего уровня и уширения связанного со временем жизни остовной вакансии.

Из качественного сравнительного анализа теоретических и экспериментальных спектров на рис.5 видно, что лучшее согласие с экспериментом достигается при использовании электронных конфигураций, соответствующих нейтральному атому меди, и при самосогласованном расчете. Чувствительность же спектров к различным электронным конфигурациям 45- и 4р-орбиталей проявляется главным образом в области основного максимума (белой линии, 8985-9000 эВ).

8986 8988 8990 8992 8994 8996 8998 9000 9002 9004

Энергия, эВ

Рис. 5. - ХЛКББ спектры Си^-края аквакомплексов меди: черная линия -экспериментальный спектр водного раствора СиС12^2(Н2О) концентрации 0,1

М, зеленая линия - теоретический спектр для 3d94s04p2электронной конфигурациимеди, синяя линия - теоретический спектр для 3d94s14p1-электронной конфигурации меди, голубая линия - 3d94s24p0электронной

конфигурации меди

Рис.5 показывает область белой линии XANES спектров Си К-края аквакомплексов меди: черная линия - экспериментальный спектр водного раствора СиС12^2(Н20) концентрации 0,1 М, цветные линии - теоретические спектры, полученные при использовании различных электронных конфигураций иона меди (см. подписи на картинке).

При "перебросе" электронов с 4s на 4р наблюдается увеличение интенсивностей пиков A' и А" и уменьшение пика А. Поскольку в эксперименте в области энергий, соответствующей пику A", наблюдается лишь слабовыраженное плечо S, значительно уступающее в интенсивности главному максимуму, то ясно, что расчетный спектр для электронной конфигурации [Ar]3d94s24p0 лучше остальных соответствует эксперименту.

Заключение

По итогам проведенного исследования, можно сделать вывод, что более правильным подходом является использование самосогласованного потенциала. Однако, процедура самосогласования, реализованная в коде FDMNES, использует МТ-приближение, не вполне обоснованное для рассматриваемой нами задачи. В связи с этим, а также с тем, что различие между спектрами, полученными для самосогласованного потенциала и для электронной конфигурации [Ar]3d94s24p0, незначительны, подход, основанный на искусственном перебросе электрона на 4^-орбиталь является вполне приемлемым, и дальнейшие расчеты будут проводитьсяв его рамках.

Благодарность

Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 14-02-31790 мол_а и Совета по грантам Президента Российской Федерации (стипендия для молодых учёных СП-2373.2013.5).

Авторы благодарны Я.В.Зубавичусу и коллективу станции Структурное материаловедение Курчатовского института за помощь в проведении экспериментов и В.Н.Дацюкуза выделение компьютерного времени вычислительного кластера ЮГИНФО для проведения расчётов.

Литература

1. V. A. Streltsov, S. J. Titmuss, V. C. Epa, K. J. Barnham, C. L. Masters, andJ. N. Varghese // Biophysical Journal.2008. № 95. p. 3447.

2. А. P. Chaynikov, M. A. Soldatov, V. Streltsov, and A. V. Soldatov// Journal of Physics: Conference Series.2013. № 430.p. 012042.

3. Кременная М.А., Солдатов М.А., Чайников А.П., Подковырина Ю.С., Бугаев А.Л., Ломаченко К.А., Кравцова А.Н. Рентгеноспектральное исследование и компьютерное моделирование локальной атомной структуры центра связывания иона меди в бета амилоиде // Инженерный вестник Дона, 2013, № 2, URL: ivdon.ru/ru/magazine/archive/n2y2013/1640.

4. Положенцев О.Е., Гуда А.А., Сафонова О.В., Ван Бокховен Д.А.и Солдатов А.В. Методика in-situ исследования катализаторов с помощью спектроскопии рентгеновского поглощения// Инженерный вестник Дона, 2013, № 2, URL:ivdon.ru/ru/magazine/archive/n2y2013/1649.

5. L. E. OrgelandJ. D. Dunitz // Nature. 1957. № 179. p. 462.

6. H. OhtakiandT. Radnai // ChemicalReviews. 1993. № 93. p. 1157.

7. M. Magini //The Journal of Chemical Physics. 1981. №74. p. 2523.

8. S. Ansell, R. H. Tromp, and G. W. Neilson // Journal of Physics: Condensed Matter. 1995. №7. p. 1513.

9. M. Benfatto, P. D'Angelo, S. Delia Longa, and N. Pavel // Physical Review B. 2002. № 65.p.174205.

10. O. Bunau and Y. Joly // Journal of Physics: Condensed Matter. 2009. № 21. p. 345501.

11. L. Hedin and B. I. Lundqvist //Journal of Physics C: Solid State Physics. 1971. № 4. p. 2064.

12. P. D'Angelo, E. Bottari, and M. R. Festa // The Journal of Chemical Physics. 1997. № 107. p. 2807.

References

1. V. A. Streltsov, S. J. Titmuss, V. C. Epa, K. J. Barnham, C. L. Masters, andJ. N. Varghese.Biophysical Journal. 2008. № 95. p. 3447.

2. А. P. Chaynikov, M. A. Soldatov, V. Streltsov, and A. V. Soldatov. Journal of Physics: Conference Series. 2013. № 430. p. 012042.

3. Kremennaja M.A., Soldatov M.A., Chajnikov A.P., PodkovyrinaJu.S., Bugaev A.L., Lomachenko K.A., Kravcova A.N. Inzenernyj vestnik Dona (Rus), 2013, № 2, URL: ivdon.ru/ru/magazine/archive/n2y2013/1640.

4. Polozhencev O.E., Guda A.A., Safonova O.V., Van Bokhoven D.A. and Soldatov A.V. Inzenernyj vestnik Dona (Rus), 2013, № 2, URL: ivdon.ru/ru/magazine/archive/n2y2013/1649.

5. L. E. OrgelandJ. D. Dunitz.Nature. 1957. № 179. p. 462.

6. H. OhtakiandT. Radnai.ChemicalReviews. 1993. № 93. p. 1157.

7. M. Magini.The Journal of Chemical Physics. 1981. №74. p. 2523.

8. S. Ansell, R. H. Tromp, and G. W. Neilson. Journal of Physics: Condensed Matter. 1995. №7. p. 1513.

9. M. Benfatto, P. D'Angelo, S. Della Longa, and N. Pavel.Physical Review B. 2002.№ 65.p. 174205.

10. O. Bunauand Y. Joly. Journal of Physics: Condensed Matter. 2009. № 21. p. 345501.

11. L. Hedin and B. I. Lundqvist. Journal of Physics C: Solid State Physics. 1971. № 4. p. 2064.

12. P. D'Angelo, E. Bottari, and M. R. Festa.The Journal of Chemical Physics. 1997. № 107. p. 2807.