Математические модели квантовых эффектов магнитоплазмоникиНИР

Mathematical models of quantum effects of magnetoplasmonics

Соисполнители НИР

МГУ имени М.В.Ломоносова Координатор

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

грант РНФ

Этапы НИР

# Сроки Название
1 1 января 2022 г.-31 декабря 2022 г. Математические модели квантовых эффектов магнитоплазмоники
Результаты этапа: 1. Разработана и реализована компьютерная модель магнитоплазмонной наночастицы с цилиндрической симметрией с учетом пространственной дисперсии материала золотой оболочки. Разработана концепция вычисления сечения поглощения электромагнитной энергии частицей в аналитическом виде через амплитуды дискретных источников, позволяющая избежать численного интегрирования поля по внешней поверхности слоя. Показано, что за счет уменьшения толщины металлического слоя до 2нм и использования вещества ядра с большим индексом рефракции Fe2O3 вместо Fe3O4 удается обеспечить сдвиг максимума сечения поглощения в область прозрачности человеческих тканей 750-900нм. Это обстоятельство позволило при исследовании пространственных структур (слоистая сфера, сфероид) ограничиться минимальными размерами ядра 16нм и толщины оболочки в 2нм, что существенно сэкономило ресурсы и время поиска оптимальных параметров частицы. В результате моделирования определено, что деформация формы частицы существенно меняет как амплитуду сечения поглощения, так и расположение его максимума при различных направлениях внешнего возбуждения. Установлено, что учет пространственной дисперсии в металлической оболочке приводит к снижению интенсивности поглощения и сдвигу положения максимума в коротковолновую область. На основе метода Дискретных источников построена математическая модель 3х мерной слоистой наночастицы FenOm@Au, учитывающая пространственную дисперсию в золотой оболочке на основе теории обобщенного нелокального отклика (GNOR). Проведен анализ дополнительных граничных условий на границах раздела сред слоистой частицы, обеспечивающих однозначную разрешимость граничной задачи дифракции. Исследовано влияние пространственной дисперсии на коэффициент усиления (КУ) поля на внешней границе слоистой частицы, расположенную в воде и состоящей из магнитного ядра и золотой оболочки FenOm@Au. В результате проведенных исследований было установлено, что: 1. Для сферических частиц Fe3O4@Au при увеличении диаметра ядра или уменьшении толщины золотой оболочки (до 2нм) наблюдается увеличение КУ со сдвигом максимального значения в ближнюю инфракрасную область. Учет пространственной дисперсии в золотой оболочке приводит к существенному снижению величины КУ (почти на 50%) и небольшому смещению в область коротких волн на 10нм по сравнению с отсутствием пространственной дисперсии. 2. Использование материала Fe2O3 в ядре приводит к дальнейшему смещению в инфракрасную область и увеличению КУ на 50% по сравнению с Fe3O4 . Это является следствием более высокой вещественной части показателя преломления у Fe2O3 при практическом отсутствии мнимой части. 3. При асимметричном расположении ядра по отношению к оболочке со средней толщиной 2нм (1.3нм и 2.7нм), было проведено осреднение КУ по углам падения и поляризациям внешнего возбуждения. В этом случае учет пространственной дисперсии приводит к снижению амплитуды КУ в два раза при незначительном смещении максимума (до 12нм) в сторону коротких волн. Метод Дискретных источников был адаптирован для исследования задачи поглощения электромагнитной энергии гибридной несферической частицей, составленной из магнитного ядра FenOm@Au и золотой оболочки с учетом пространственной дисперсии золота в рамках теории (GNOR). При этом сечение поглощения вычислялось без использования процедуры интегрирования по поверхности, посредством оптической теоремы с аналитическим вычислением сечения рассеяния аналогично случаю цилиндрической частицы. Убедительно показано, что увеличение вытянутости частицы параллельно направлению вектора внешнего электрического поля при сохранении объема ядра и толщины оболочки позволяет существенно увеличить величину поглощенной энергии и одновременно сдвинуть максимум поглощения в область прозрачности биологических тканей 750-900нм. Учет же пространственной дисперсии приводит к снижению величины поглощенной энергии до 30% при сдвиге в область коротких волн на существенную величину в 25nm. В результате моделирования выяснилось, что осреднение сечения поглощения по направлениям внешнего возбуждения и поляризациям влечет за собой снижение поглощенной энергии в три раза. Данное обстоятельство демонстрирует важность ориентирования магнитоплазмонных частиц перпендикулярно направлению внешнего возбуждения и позволяет, при использовании P-поляризованнго источника излучения, оптимизировать процесс фототермического нагрева на умеренных оптических мощностях, способствуя сбережению здоровых тканей и их быстрому восстановлению. Резюмируя комплекс проведенных исследований, отмечаем, что учет пространственной дисперсии является необходимым условием для разработки оптимальной стратегии и практического применения магнитоплазмонных слоистых наночастиц при осуществлении как визуализации опухолей, так и фототермического нагрева онкологических образований с целью использования умеренных оптических мощностей и сохранения окружающих здоровых тканей.
2 1 января 2023 г.-31 декабря 2023 г. Математические модели квантовых эффектов магнитоплазмоники
Результаты этапа: В отчетный период продолжались работы по исследованию однозначной разрешимости исходной граничной задачи дифракции для системы уравнений Максвелла с учетом пространственной дисперсии в рамках теории обобщенного нелокального оптического отклика (GNOR), с условием непротекания тока проводимости на поверхности интерфейса металл-диэлектрик. Для рассеивателя, ограниченного произвольной гладкой цилиндрической поверхностью было получено базовое энергетическое соотношение для обобщенной системы Максвелла. Была доказана теорема единственности для решения обобщенной граничной задачи дифракции при дополнительных условиях: • немагнитной среды; • отличной от нуля мнимой части комплексной диэлектрической проницаемости. Следует отметить, что эти условия выполняются во всем оптическом диапазоне длин волн для благородных металлов таких, как золото и серебро. Тем самым установлена единственная разрешимость исходной граничной задачи дифракции для нелокального случая. Проведено изучение влияния зависимости от длины волны внешнего излучения коэффициента диффузии в рамках теории GNOR на сечение поглощения и коэффициент усиления поля магнитоплазмонной частицы, состоящей из магнитного ядра и золотой оболочки. В результате моделирования установлено, что использование динамического коэффициента диффузии приводит к росту сечения поглощения и коэффициента усиления поля, по сравнению с постоянным коэффициентом диффузии, сохраняя при этом положение плазмонного резонанса (ПР) в области длин волн. Как удалось установить, это различие может превышать 50% от амплитуды. Оно особенно заметно при оценке коэффициента усиления поля на поверхности плазмонной оболочки. Было исследовано влияние деформации магнитоплазмонной слоистой частицы на оптические характеристики, с учетом теории GNOR и динамического коэффициента диффузии. Рассматривались деформации двух видов: смещение центра ядра по отношению к оболочке и деформация формы частицы в сфероидальную (нанорисинки). При конкретном рассмотрении смешения центра ядра на 50% от толщины пленки, было установлено, что подобное изменение толщины золотой пленки лишь незначительно влияет на уровень поглощенной энергии. С другой стороны оказалось, что подобная асимметрия приводит к тому, что коэффициент усиления поля на внешней границе пленки возрастает в 1.5 раза. При деформации же формы сферической слоистой частицы в сфероидальную коэффициент усиления поля возрастает в три раза по сравнению со сферически симметричным случаем. Таким образом было установлено, что деформации формы магнитоплазмонных частиц не оказывают негативного влияния на амплитуду ПР, а также заметного влияния на его положение в ближней инфракрасной оптической области. Был проведен анализ влияния окружающей среды на сечение поглощения магнитоплазмонной наночастицы в случае различных тканей человеческого организма: печень, легкие, сердце, мозг, желудок и женская грудь. Было установлено, что с увеличением плотности внешней среды максимум сечения поглощения и коэффициента усиления поля немного сдвигается в инфракрасную область, а амплитуда незначительно возрастает. Все исследования проведены в диапазоне окна прозрачности биологических тканей 750-900нм с учетом пространственной дисперсии золота. Показано, что все перечисленные выше среды не выводят сечение поглощения за пределы окна прозрачности. Рассмотрена задачи дифракции плоской электромагнитной волны на кластере двух слоистых наноцилиндров, состоящих из магнитоплазмонного ядра и золотой оболочки. Учет эффектов пространственной дисперсии в золотой оболочке осуществлялся в рамках теории GNOR. Используя оптическую теорему, устанавливающую связь между сечением рассеяния, сечением поглощения в золотой пленке и сечением экстинкции, удалось получить представление для сечения поглощения в аналитическом виде, не прибегая к численному интегрированию по внешней поверхности частиц кластера. Это обстоятельство позволило существенно упростить и ускорить моделирования происходящих процессов в частотной области. На основе модифицированной схемы метода дискретных источников проведено исследование влияния взаимного расположения частиц, их деформации и эффекта пространственной дисперсии в золотой оболочке на поведение сечения поглощения в золотой оболочке и коэффициента усиления ближнего поля на внешней поверхности пленки. Показано, что положение максимума сечения поглощения можно сдвигать в окно прозрачности биологических тканей за счет вариации расстояния между частицами и выбора материала ядер. Установлено, что учет пространственной дисперсии в золотой оболочке приводит к снижению амплитуды ПР и небольшому сдвигу его положения в коротковолновую область (blue shift), не выводя, однако, за пределы окна прозрачности биологических тканей. Был подробно изучен вопрос о влиянии расстояния между частицами кластера на оптические характеристики. В результате моделирования было определено, что с увеличением зазора между частицами кластера амплитуда ПР уменьшается, а его положение незначительно смещается в коротковолновую область, все еще оставаясь в окне прозрачности. Было также установлено, что интенсивность в промежутке между частицами существенно возрастает при их сближении. Также продемонстрировано, что учет эффекта пространственной дисперсии приводит к незначительному снижению интенсивности ПР и его сдвигу в область коротких волн. Было изучено влияние асимметрии геометрии частиц, составляющих кластер, на сечения поглощения и коэффициента усиления поля в случае сдвига центра ядра относительно геометрического центра слоистой частицы вплоть 50% толщины золотой оболочки. Оказалось, что амплитуда ПР в этом случае может смещается в длинноволновую область на величины до 35нм, а значение максимума уменьшается на 20%. Однако, подобное смещение не является критическим, так как ПР все еще располагается в окне прозрачности даже при учете эффекта пространственной дисперсии. Основные выводы проведенного исследования в2023г.: 1.Эффект пространственной дисперсии оказывает заметное влияние на положение максимумов как сечения поглощения, так и коэффициента усиления поля на внешней границе золотой оболочки в частотной области. 2.Учет динамического коэффициента диффузии, лишь незначительно снижает амплитуды плазмонных резонансов и осуществляет сдвиг с область коротких волн - blue shift. Магнитоплазмонные наночастицы синтезированные таким образом, что их ПР попадают в середину окна прозрачности тканей человека ~820нм, будут сохранять положения ПР в окне прозрачности биотканей даже при деформациях. 3.При использовании кластера частиц нет необходимости для дополнительного покрытия золотой пленки какими либо материалами, чтобы избежать слипания частиц, так как это может приводить лишь к снижению эффективности осуществления фототермической терапии. 4. Синтезированные магнитоплазмонные наночастицы могут использоваться при терапии любых органов человека, так как увеличение плотности окружающей среды только увеличивает поглощение электромагнитной энергии и усиливает ближнее поле. Публикации. 1. Еремин Ю.А., Лопушенко В.В. О влиянии динамического коэффициента диффузии с параметром Фийбельмана на квантовый эффект нелокальности гибридных плазмонных наночастиц//Математическое моделирование. 2024. Т. 36. №1. с.3-8. 2. Еремин Ю.А., Пензарь А.С. Влияние пространственной дисперсии в металле на оптические характеристики кластера магнетоплазмонных слоистых наночастиц//Вестник МГУ, сер.3. Физика. Астрономия. 2023. Т.78. №5. 2350104. По результатам исследований подготовлены и опубликованы 2 статьи в журналах, индексируемых в Scopus и 4 доклада на конференции. Информационный ресурс, посвященный проекту: https://istina.msu.ru/projects/436597335/

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

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