Математические модели в инженерии
Моделирование термостабильности органических светоизлучающих устройств с использованием методов Ab Initio
Table of Contents
Введение в OLED и вызов термической стабильности
Органические светоизлучающие устройства (OLED) стали краеугольным камнем современной технологии отображения и освещения, ценимые за их гибкость, энергоэффективность и способность производить яркие цвета без необходимости подсветки. Несмотря на эти преимущества, широкое внедрение OLED в требовательных приложениях, таких как автомобильное освещение, дисплеи большой площади и вывески на открытом воздухе, по-прежнему ограничено их восприимчивостью к термической деградации. Когда OLED работает, нагревание джоуля повышает внутреннюю температуру органических слоев, ускоряя химические реакции, которые приводят к потере яркости, сдвигу цвета и возможному отказу устройства. Понимание и прогнозирование того, как OLED материалы ведут себя при тепловом напряжении, поэтому необходимо для проектирования следующего поколения высокопроизводительных устройств с длительным сроком службы.
Термическая стабильность в OLED-дисплеях не является единственным свойством, а представляет собой сложное взаимодействие молекулярной жесткости, прочности связи, геометрии упаковки и энергии возбужденных состояний. При повышенных температурах органические молекулы могут подвергаться фазовым переходам (например, кристаллизации аморфных пленок), расщеплению химических связей или диффузии допантов и примесей. Даже умеренное повышение температуры - от комнатной температуры до 80-100°C - может резко сократить срок службы устройства. Следовательно, исследователи обратились к методам вычислительных материалов, в частности ab initio , чтобы получить атомистическое понимание этих механизмов деградации и направлять рациональный дизайн более надежных органических полупроводников.
Методы Аб Инитио в моделировании материалов: от первых принципов до прогнозного дизайна
Методы Ab initio, полученные из квантовой механики, позволяют вычислять свойства материала, не полагаясь на эмпирические параметры или экспериментальную подгонку. Решая уравнение Шрёдингера с множеством электронов — или их приближение — эти методы обеспечивают строгое описание электронной структуры, энергии и динамики молекул и твердых тел. В контексте OLED методы ab initio используются для прогнозирования энергии диссоциации связей, потенциалов ионизации, сродства электронов и энергетических барьеров для конформационных изменений или химических реакций, которые происходят при тепловом напряжении. Наиболее широко принятой структурой является теория функционала плотности (DFT), но методы волновой функции на основе пост-Хартри-Фока также используются, когда требуется более высокая точность.
Функциональная теория плотности (DFT) и ее роль в стабильности OLED
DFT переформулирует многоэлектронную проблему с точки зрения плотности электронов, делая возможными крупномасштабные квантовые вычисления для систем, содержащих сотни атомов — типичный размер для молекул-излучателей OLED или комплексов-хозяев. Точность расчета DFT зависит от выбора функциональной корреляции обмена (]Mardirossian & Head-Gordon, 2017] Для органических молекул гибридные функции, такие как B3LYP, PBE0 или разделенный по диапазону ωB97X-D, являются общими, потому что они уравновешивают вычислительные затраты с надежным прогнозированием прочности связей и энергий возбуждения. При скрининге потенциальных OLED-материалов для термической стабильности DFT может вычислять энтальпии ключевых реакций разрыва связей (например, C-C, C-N или металлолигандные связи в фосфоресцирующих излучателях) при комнатной температуре и выше. Сравнивая эти энергии, исследователи могут идентифицировать молекулярные фрагменты, которые наиболее уязвимы для теплового
Помимо энергии наземного состояния, DFT также может моделировать влияние температуры посредством расчетов гармонической частоты. Вычисленный вибрационный спектр дает энергию нулевой точки и термодинамические поправки (энталпии, энтропии, свободной энергии) в качестве функции температуры, что позволяет строить профили реакций в условиях эксплуатации. Это особенно важно для OLED, потому что излучающее состояние часто является триплетным экситоном в фосфоресцирующих излучателях, а энергетический разрыв между триплетом и диссоциативной потенциальной энергетической поверхностью может быть небольшим. Поиски точки пересечения на основе минимальной энергии DFT показали, что термическая популяция низменных триплетных металлцентрических состояний является основным путем деградации для зеленых и синих фосфоресцирующих излучателей (]Hedley et al., 2019 .
За пределами DFT: волновая функция и пост-хартрей-фоковские методы
Хотя DFT эффективен, он иногда не подходит для систем с сильной статической корреляцией, таких как переходные состояния разрыва связи или промежуточные соединения с открытой оболочкой, которые распространены в термической деградации. В таких случаях методы волновой функции, такие как сопряженный кластер (CCSD(T)) или теория многоканальной возмущения (CASPT2, NEVPT2), обеспечивают энергии эталонного качества. Хотя эти методы являются вычислительно дорогостоящими и обычно ограничены небольшими модельными соединениями (20-50 атомов), они служат золотыми стандартами для проверки результатов DFT и обучения линейно-масштабным подходам. Например, энергии диссоциации связи общих OLED-связей (например, карбазол-триазиновые связи в материалах-хозяевах), рассчитанные на уровне CCSD (T) / CBS, могут использоваться для параметризации более дешевых методов DFT для высокопроизводительного скрининга (]Nayak et al., 2019 . Кроме того, зависящий от времени DFT (TD-DFT) и его аналог волновой
Молекулярная динамика Ab Initio и моделирование конечных температур
Статические расчеты электронной структуры не фиксируют полную сложность теплового движения. При рабочих температурах молекулярные вибрации могут вызывать растяжение связи, кручение вращения и даже изомеризацию, все из которых влияют на электронные свойства и стабильность. Молекулярная динамика Ab initio (AIMD) имитирует временную эволюцию атомов под квантовыми силами, обычно используя DFT в качестве двигателя электронной структуры. В AIMD используется приближение Борна-Оппенгеймера: для каждого ядерного шага электронное основное состояние вычисляется самосогласованно. Траектории 10-100 пс достаточны для наблюдения редких событий, таких как разрыв связи при применении усиленных методов отбора проб (например, метадинамика).
Молекулярная динамика Car-Parrinello, вариант, который интегрирует электронные и ядерные степени свободы, была использована для изучения термического разложения OLED-материалов, таких как Alq3 (tris(8-hydroxyquinolinato)aluminum), классический зеленый излучатель. Эти моделирования показали, что связь Al-O значительно ослабевает выше 100°C, что приводит к отслоению лиганда и образованию кластеров Al2O3, которые действуют как места утечки тока (]Pingel et al., 2007 . Для современных N-гетероциклических карбеновых (NHC) комплексов иридия AIMD показала, что связь Ir-C удивительно стабильна до 200°C, объясняя исключительную термическую устойчивость этих глубоко-синих излучателей.
Моделирование путей термического разложения: от облигаций к отказу устройства
Термическая деградация в OLEDs — это многоступенчатый процесс, который начинается на молекулярном уровне и в конечном итоге проявляется как макроскопическая потеря эффективности.Методы Ab initio позволяют исследователям связать эти масштабы, моделируя элементарные реакции, которые инициируют сбой.
Энергии диссоциации облигаций и радикальное формирование
Наиболее простой механизм термической деградации - гомолитическая или гетеролитическая расщепление связей. Используя DFT или CCSD(T), можно рассчитать энталпию диссоциации связей (BDE) для каждого ковалентного звена в органическом полупроводнике. Материалы с BDE ниже 60-80 кДж/моль обычно считаются нестабильными для длительной работы при 80°C. Например, кислородно-углеродные связи в эфирных материалах-хозяевах часто имеют BDE около 200 кДж/моль, но присутствие близлежащих электронно-выводящих групп может снизить это на 30-50 кДж/моль, делая их восприимчивыми к термическому гомолизу. Полученные радикалы могут затем реагировать с кислородом или с соседними молекулами, создавая люминесцентные гашения. Систематическое исследование DFT общих OLED функциональных групп (например, триазин, карбазол, трифениламин) предоставило базу данных стабильности, которая позволяет синтетическим химикам расставлять приоритеты для надежных строительных блоков ([FLT:
Возбужденная государственная динамика и тепловое закалывание
OLED-материалы не деградируют только от основного состояния; возбуждённое состояние (сингл или триплет) часто имеет различную потенциальную энергетическую поверхность с более низкими барьерами для диссоциации. Например, состояние триплета синего флуоресцентного излучателя может заполнять диссоциативное состояние посредством межсистемного пересечения или прямого спин-орбитального соединения. Расчеты Ab initio кривых потенциальной энергии по соответствующим координатам (например, растяжение связи) могут идентифицировать «горячие» режимы, которые активируются при повышенных температурах. Кроме того, тепловые эффекты могут увеличивать скорость нерадиационного распада (термическое закалывание) путем открытия вибрационных релаксационных каналов. TD-DFT и многоканальные методы могут использоваться для отображения топологии грунта и возбуждённых состояний популярных излучателей, таких как 4,4'-бис(]-карбазолил)-1,1'-бифенил (CBP) и трис(2-фенилпиридин)иридий (I
Интеграция вычислений с экспериментом для ускоренного открытия
В то время как методы ab initio дают фундаментальные идеи, их реальная сила заключается в руководстве экспериментальными усилиями. Вычислительный скрининг виртуальных молекулярных библиотек до синтеза может сэкономить месяцы пробной и ошибочной работы. Появились две дополнительные стратегии: высокопроизводительный виртуальный скрининг (HTVS) и ускоренное моделирование машинного обучения (ML).
Высокопроизводительный скрининг для термостабильных OLED-материалов
Рабочие процессы HTVS объединяют автоматизированные расчеты DFT с хемоинформатикой для оценки тысяч молекул-кандидатов на такие свойства, как уровни энергии HOMO/LUMO, энергия триплета, температура перехода стекла (T]g) и BDE слабейшей связи. Например, скрининг 20 000 производных карбазола с использованием DFT (B3LYP/6-31G*) выявил 2200 кандидатов как с высокой энергией триплета (>2,8 эВ), так и по меньшей мере одну связь с BDE > 300 кДж/моль (]Hachmann et al., 2019. Последующие тесты синтеза и термостабильности (термогравиметрический анализ и DSC) подтвердили, что верхние 10 соединений имели температуры разложения >400°C по сравнению с ~350°C для обычных материалов. Этот подход был расширен до фосфоресцирующих материалов-хозяев, где комбинация высоких Tg
Машинное обучение и суррогатные модели
Даже с дешевыми функциями DFT скрининг десятков тысяч молекул требует вычислительных усилий (несколько тысяч ядер-часов на партию). Чтобы преодолеть это, исследователи обучили модели машинного обучения на существующих данных DFT прогнозировать дескрипторы термостабильности с точностью почти до DFT при доле стоимости. Графовые нейронные сети, например, напрямую изучают взаимосвязь между структурой молекулярного графа и BDE, позволяя оценивать миллионы соединений за минуты. В недавнем исследовании использовался ансамбль из 100 нейронных сетей, обученных на 50 000 DFT-расчетных BDE, чтобы предсказать самую слабую связь в 1 миллионе гипотетических излучателей OLED. Модель ML воспроизвела рейтинг DFT со средней абсолютной ошибкой всего 8 кДж/моль, и 70% из 10 000 лучших кандидатов впоследствии были подтверждены как стабильные с помощью моделирования AIMD (]Choi et al., 2022 . Такие гибридные рабочие процессы DFT + ML теперь являются самыми современными для обнаружения материалов в области OLED.
Вызовы и направления будущего в моделировании термостабильности OLED
Несмотря на значительный прогресс, остаются несколько проблем. Во-первых, вычислительная стоимость точных методов ab initio (CCSD(T), CASPT2) является непомерной для систем с более чем ~50 атомами, но реальные OLED-материалы часто содержат > 100 атомов. Алгоритмы линейного масштабирования и методы на основе фрагментов (например, MFCC, GMBE) активно разрабатываются для расширения высокоуровневой точности на более крупные системы. Во-вторых, большинство симуляций предполагают идеальные условия газовой фазы или статические кристаллические структуры, тогда как рабочая среда OLED включает интерфейсы, аморфную морфологию и носители заряда. Многомасштабное моделирование, которое сочетает вычисления ab initio с классическими силовыми полями (QM/MM) или континуумными моделями, необходимо для захвата полного поведения устройства. В-третьих, длинные временные масштабы деградации (часы до лет) недоступны для стандартного AIMD. Усиленные временные масштабы деградации (часы до лет) являются недоступными для стандартных AIMD. Усиленные временные масштаб
Еще одним рубежом является включение экситонной динамики и электронно-фононной связи в предсказания термостабильности. Неадьябатическая связь между электронными состояниями может привести к теплоактивируемой передаче заряда, что в свою очередь может ускорить деградацию. Разработки в режиме реального времени TD-DFT и динамика поверхностного скачка начинают предоставлять эту информацию, но они остаются на стадии исследования. Наконец, интеграция данных ab initio в модели процессов осаждения пара и изготовления устройства поможет оптимизировать условия производства, чтобы минимизировать тепловое напряжение с самого начала.
Заключение
Методы Ab initio стали незаменимыми для понимания и улучшения термической стабильности органических светоизлучающих устройств. Обеспечивая строгое, первое описание прочности связей, путей реакции и динамического поведения при конечной температуре, эти вычислительные методы позволяют исследователям идентифицировать горячие точки деградации, экранные виртуальные библиотеки и проектировать новые материалы с исключительной надежностью. Продолжающаяся разработка более точных и эффективных методов квантовой химии в сочетании с ускорением машинного обучения и многомасштабной интеграцией обещает вступить в эпоху действительно прогностических материалов для OLED. По мере роста спроса на долговечные, высокосветовые дисплеи и освещение, синергия между моделированием ab initio и экспериментальной валидации будет ключом к преодолению узкого места термической стабильности и раскрытию полного потенциала технологии OLED.