Теоретическое моделирование электронных и магнитных свойств эндоэдральных фуллеренов и их производныхНИР

Theoretical modelling of electronic and magnetic properties of endohedral fullerenes and their derivatives

Источник финансирования НИР

грант РНФ

Этапы НИР

# Сроки Название
1 8 февраля 2021 г.-29 декабря 2021 г. Исследование региохимических тенденций последовательной функционализации рассматриваемых эндоэдральных молекул.
Результаты этапа: 1) С помощью метода функционала плотности рассмотрено присоединение CF2-группы к Sc3N@C78. Выяснено, что энергии изомеров CF2-производных с M3N@C78 существенно зависят от места присоединения заместителя, что согласуется с имеющимися литературными данными по близким реакциям. В среднем энергия присоединения оказывается существенно ниже, чем в M3N@C80 и С60. Обнаружено, что реакция происходит в условиях кинетического контроля, причем образуется сравнительно термодинамически невыгодный изомер Sc3N@C78(CF2), т.е. известный эмпирический принцип корреляции термодинамических и кинетических параметров сходных реакций нарушается. Механизм реакции функционализации является не карбеновым, как предполагалось ранее, а нуклеофильным. Его первая безбарьерная стадия состоит в присоединении к углеродному каркасу аниона CF2Cl-, вслед за которым происходит внутримолекулярное вытеснение хлорид-аниона с замыканием цикла. Интересным аспектом строения CF2-производных фуллеренов является возможность существования открытой или закрытой конформации CF2-групп, переключаемых при измерении зарядового состояния. Так, экспериментально образующийся изомер Sc3N@C78(CF2) является единственным, у которого в нейтральном состоянии эта конформация является закрытой, тогда как в анионном состоянии образуется дополнительный открытый минимум. Выяснено, что производные Sc3N@C78 близки к Sc3N@C80 по энергии нижних синглетных и триплетных возбужденных состояний, но при этом они обладают меньшими на несколько десятых эВ сродством к электрону и энергией ионизации, вследствие чего они могут оказаться более реакционноспособны. 2) Для производных молекул M3N@C80 с кластерами Sc2DyN и ScDy2N проведены исследования применимости различных эффективных остовных потенциалов к определению набора ориентаций кластера внутри углеродного каркаса с помощью метода функционала плотности. Сопоставлены стандартный остовный потенциал для диспрозия со включением в остов 28 электронов и часто используемый для упрощения расчетов псевдопотенциал для трехвалентного иона с большим остовом, заключающим в себя f-оболочку. Рассмотрены изомеры M3N@C80 с двумя типами внешнесферных заместителей - CF2-группой и пирролидиновым циклом. Найдено, что при использовании псевдопотенциалов с большим остовом разница в ориентациях кластеров Sc2DyN и ScDy2N связанна с не слишком сильным изменением энергии. Однако переход к малому остовному потенциалу приводит к неожиданным изменениям, и в некоторых случаях разница между различными ориентациями кластера начинает превышать 1 эВ. При этом иногда наблюдается сходимость к достаточно низким по энергии решениям с нетипичными спиновыми заселенностями атомов диспрозия, отличными от пяти неспаренных электронов. Полученные результаты свидетельствуют о необходимости проверки применимости используемых остовных потенциалов и других параметров расчета в целом по данным расчетных методов более высокого уровня. 3) Исследовано ориентационное распределение кластеров Sc2DyN и ScDy2N внутри углеродного каркаса С80 с CF2- и пирролидиновым заместителями. Полученные данные показывают, что для наиболее устойчивых структур с ориентацией атома диспрозия против внешнесферного заместителя существует несколько сопоставимых по энергии минимумов, отвечающих поворотам кластера около оси N-Dy. Структуры с ориентацией против внешнесферного заместителя атома скандия несколько менее выгодны, и при этом также существует целый ряд локальных минимумов, для которых в сторону заместителя не направлен ни один из атомов металла. Однако эти структуры характеризуются существенно более высокой энергией (приблизительно на 100 кДж/моль и более даже при использовании псевдопотенциалов с большим остовом), и поэтому являются лишь метастабильными. В ближайшее время предполагается исследование характеристик ориентационных переходов между полученными конформациями и их электрических и магнитных свойств. 4) Проведены многоконфигурационные расчеты скандий-диспрозиевых систем вида Sc2DyN@C80(CF2) с различными активными пространствами методом МКССП в полном активном пространстве (CASSCF) и квазивырожденной многоконфигурационной теории возмущений XMCQDPT2. В данных расчетах использован, конечно, базисный набор с малым остовом. В качестве базовой составляющей активного пространства используется f-оболочка диспрозия. Исследовано добавление к ней других возможно значимых орбиталей, таких как неподеленная пара атома азота, вакантные s- и p-уровни диспрозия и др. Обнаружено, что орбитали азота и скандия не вносят заметный вклад в статическую электронную корреляцию в M3N@C80(CF2), тогда как ситуация с другими орбиталями менее ясна. Данные XMCQDPT2 показывают, что имеет место весьма тесное вырождение 11 из 21 состояния, связанного с различными схемами заселения f-оболочки девятью электронами, причем полученная картина достаточно близка к картине МКССП. Анализ полученных нижних состояний показывает, что они в существенной степени многоконфигурационны, и поэтому представляется необходимым прямое сравнение результатов метода функционала плотности и реперных результатов теории возмущений для тестовых наборов избранных структур.
2 10 января 2022 г.-30 декабря 2022 г. Исследование ориентаций эндоэдрального кластера и многоконфигурационное описание эндоэдральных фуллеренов
Результаты этапа: 1. Найдено, что использование двухэкспонентных базисных наборов для углеродного каркаса с внешнесферными заместителями и трехэкспонентных базисных наборов для эндоэдрального кластера является оптимальным для расчетного предсказания геометрических, электронных и магнитных свойств функционализированных эндоэдральных систем как методом функционала плотности, так и методом многоконфигурационной теории возмущений XMCQDPT2. 2. С помощью сравнения результатов метода функционала плотности и метода теории возмущений показано, что использование в методе функционала плотности эффективных остовных потенциалов с большим остовом, заключающим в себя f-оболочку, позволяет удовлетворительно предсказывать оптимальные конфигурации эндоэдрального кластера внутри углеродного каркаса, обеспечивая при этом хорошую сходимость процедуры самосогласования. 3. Получены полные наборы устойчивых ориентаций эндоэдральных кластеров Sc2DyN и ScDy2N внутри ряда функционализированных углеродных каркасов C80. Найдено, что электронные свойства, такие как энергии граничных состояний, а также степень переноса заряда между эндоэдральным кластером и углеродным каркасом, достаточно слабо зависят от ориентации кластера внутри каркаса. Выявлены наборы наиболее устойчивых конфигураций кластера, подтверждены предположения о предпочтительности ориентации одного из атомов металла против внешнесферного заместителя. 4. Предложен способ расчета наборов нижних электронных состояний рассматриваемых эндоэдральных систем, учитывающий спин-орбитальные взаимодействия и эффекты динамической корреляции, где последния вводятся в виде получаемых с помощью теории возмущений поправок к диагональным элементам спин-орбитального гамильтониана, построенного с помощью метода МКССП в полном активном пространстве в базисе собственных состояний спина. 5. Показано, что многоконфигурационное описание нижних электронных состояний систем с эндоэдральным кластером Sc2DyN требует включения в активное пространство только f-оболочки диспрозия и рассмотрения только порождаемых ею секстетных состояний, тогда как влияние остальных состояний может быть рассмотрено пертурбативно, т.е. взаимодействие с состояниями, отвечающими другим степеням переноса заряда между кластером и каркасом несущественно. При этом обнаружено, что величины эффектов динамической корреляции для нижних электронных состояний достаточно близки друг к другу, и их имеет смысл учитывать только при рассмотрении спин-орбитальных взаимодействий для предсказания магнитных свойств, где разница между поправками по теории возмущений для различных состояний уже может быть значимой. 6.Показано, что эффекты ориентации иона диспрозия относительно внешнесферных заместителей могут оказывать существенное влияние на магнитное поведение молекул, благоприятствуя либо низкоспиновым, либо высокоспиновым состояниям в зависимости от того, с какой областью углеродного каркаса соседствует атом диспрозия. Таким образом, может играть ключевую роль подбор оптимальных ориентирующих внешнесферных заместителей. 7. Найдено, что в системах с эндоэдральным кластером ScDy2N взаимодействие между ионами диспрозия очень слабо, и сложение спиновых моментов их секстетных состояний, приводит к очень большой группе вырожденных состояний различной мультиплетности, описание спин-орбитального взаимодействия внутри которой сопряжено с существенными техническими сложностями.
3 1 января 2023 г.-31 декабря 2023 г. Настройка магнитных свойств эндоэдральных фуллеренов под влиянием внешнесферных заместителей
Результаты этапа: 1. Найден способ описания в рамках спин-орбитального КВ на полном активном пространстве f-состояний эндоэдральных систем с двумя слабо взаимодействующими атомами диспрозия, порождающими широкий набор близких по энергии состояний с различными мультиплетностями. Осуществлена предварительная характеристика электронных уровней модельного шестизарядного кластера ScDy2N, показывающая квазивырождение 121 нижнего электронного состояния для каждой из мультиплетностей от 1 до 11. Показано, что влияние других электронных состояний, расстояние до которых превышает характерные значения матричных элементов спин-орбитального взаимодействия, несущественны. 2. Охарактеризован типичный уровень ошибок метода функционала плотности при использовании базисов с большим остовом, включающим f-оболочки, в задачах анализа конформационного разнообразия эндоэдральных молекул. Верификация по широкому набору конформеров на основании сравнения с данными многоконфигурационной теории возмущений XMCQDPT2, явно рассматривающей эффекты нединамической электронной корреляции в f-оболочке и эффекты динамической корреляции в рамках углеродного каркаса, показала, что максимальный уровень расхождения относительных энергий не превышает 15 кДж/моль, но в большинстве случаев оказывается не выше 10 кДж/моль, что следует признать хорошим совпадением. 3. Изучены особенности полифункционализации диспрозиевых эндоэдральных соединений CF2-группами. Показано, что для двух групп предпочтительной является ориентация атомов металла по внешнесферным заместителям, и эндоэдральный кластер оказывает, таким образом, ориентирующий эффект на региохимию их присоединения, который в энергетическом выражении может превышать 100 кДж/моль. С другой стороны, результатом взаимного влияния внешнесферных заместителей и эндоэдрального кластера может оказываться деформация последнего с изменением расстояния между атомами металла, т.е. особенностей взаимодействия между ними. Это может служить способом настройки электронных свойств. 4. Исследована функционализация диспрозиевых эндоэдральных соединений CH2- и N-метилпирролидиновыми аддендами. Показано сходство конформационных особенностей для продуктов CH2- и СF2-функционализации, но существенное отличие таковых в случае пирролидинового адденда, где не проявляется явное предпочтение ориентации одного из атомов металла по внешнесферному заместителю. При этом пирролидиновые производные демонстрируют существенное сужение зазора между НВМО и ВЗМО. По сравнению с нефункционализированными эндоэдральными молекулами, пирролидиновые производные демонстрируют тенденцию к сохранению магнитных свойств соединений, что позволяет использовать такой тип функционализации для создания новых материалов. 5. Начато рассмотрение влияния димеризации энндоэдральных соединений через одинарную связь на их электронные свойства, определяемые обратным переходом электрона на эндоэдральный кластер под действием образования межкаркасной связи. Показана возможность стабилизации димерных соединений в полярных средах.
4 1 января 2024 г.-31 декабря 2024 г. Поиск наиболее перспективных эндоэдральных магнитных систем
Результаты этапа: 1. Построена схема нижних электронных состояний и особенностей расщепления в нулевом поле для эндоэдральных систем с двумя атомами диспрозия на основе Dy2ScN@C80. Отмечено, что ионы диспрозия с электронной конфигурацией 4f9 взаимодействуют между собой достаточно слабо, и различные схемы сложения их спиновых моментов 5/2 дают близкие к вырождению электронные состояния. Нижние состояния имеют тенденцию демонстрировать высокие проекции спина, т.е. свойства молекулярных магнетиков, однако конкретный результат может зависеть от непосредственной ориентации эндоэдрального кластера, в существенной степени задаваемой внешнесферным заместителем, ориентирующим атомы металла, преимущественно - диспрозия, на себя. 2. Многоконфигурационные расчеты политрифторметилированных производных соединений типа DySc2N@C80, основанные на мотивах присоединения трифторметильных групп, наблюдавшихся в эксперименте при функционализации Sc3N@C80 показали, что при умеренном числе аддендов электронной конфигурацией диспрозия остается 4f9, т.е. Dy3+, и степень переноса заряда с эндоэдрального кластера на углеродный каркас сохраняется. При этом для нижних электронных состояния продолжают наблюдаться тенденции расщепления в нулевом поле, свидетельствующие о сохранении свойств молекулярных магнетиков. Такие же результаты наблюдаются для случаев функционализации C80 группами CH2 и конденсированными циклами, однако конкретные величины расщепления в нулевом поле для нижних крамерсовских дублетов могут варьироваться в пределах 300-700 см-1. Для родственных политрифторметилированных углеродных каркасов обычных фуллеренов исследованы взаимосязи между мотивом трифторметилирования и акцепторными свойствами. 3. Найдено, что для широкого набора различных внешнесферных производных эндоэдральных фуллеренов эффективные остовные потенциалы для ионов Dy3+ с большим остовом, в который заключена 4f-оболочка диспрозия, обеспечивают надежное воспроизведение энергий эндоэдрального кластера ScDy2N и Sc2DyN внутри углеродного каркаса методом функционала плотности. Качество результатов оказывается сопоставимым с расчетами с малым 28-электронным остовом, причем данные метода функционала плотности в варианте PBE0 достаточно мало отличаются от данных многоконфигурационной теории возмущений второго порядка, не меняя качественной картины. 4. Методом функционала плотности были оценены энергии активации для реориентации диспрозиевых эндоэдральных кластеров внутри функционализированного каркаса С80. Оказалось, что для соседствующих устойчивых конформаций эти энергии могут оказываться ниже, чем в исходном С80, уменьшаясь до величин порядка 10 кДж/моль, т.е. пренебрежимо малы для комнатной и пониженных температур. Показано, что устойчивые конформации эндоэдральных кластеров демонстрируют достаточно широкий набор форм, однако это мало влияет на магнитные взаимодействия эндоэдральных атомов. Найдено, что энергии активации функционализации эндоэдральных молекул CF2-группами могут несколько повышаться в случае диспрозиевых эндоэдральных кластеров по сравнению со скандиевыми до величин порядка 160-180 кДж/моль, т.е. при введении диспрозия эти реакции могут требовать несколько более жестуких условий. 5. Показана необходимость тщательного отбора и исследования устойчивости результата при применении различных схем усреднения по состояниям к диспрозиевым эндоэдральным фуллеренам при расчете относительных энергий состояний с различными заселенностями 4f-оболочки, т.е. различными распределениями электронов между эндоэдральным кластером и углеродным каркасом. Показана необходимость тщательного подбора усреднения для сбалансированного описания формы f-орбиталей в зависимости от заселенностей. 6. Найдено, что непосредственные взаимодействия между эндоэдральными кластерами в ковалентно связанных друг с другом эндоэдральных производных С80, не проявляются явным и очевидным образом, и можно предпочтительно говорить о взаимодействии металлических кластеров с сайтами внешнесферной функционализации.

Прикрепленные к НИР результаты

Для прикрепления результата сначала выберете тип результата (статьи, книги, ...). После чего введите несколько символов в поле поиска прикрепляемого результата, затем выберете один из предложенных и нажмите кнопку "Добавить".