Почему расчёт энергии отрыва электрона ошибался и как это исправили ✛

Оригинал: "Third-order Algebraic Diagrammatic Construction Renormalized by Particle-Particle Ladders"
· Arno Förster
arXiv:2608.28424v1 · 2026-08-28 · CC BY 4.0 · 1 мин
Простая поправка на «тесноту» электронов в молекуле делает расчёт почти в четыре раза точнее.

Чтобы вырвать электрон из молекулы, нужно затратить энергию — её называют потенциалом ионизации. Старый метод расчёта ADC(3) считал, что электрон уходит по пустой лестнице, не встречая других. Но на самом деле на ступеньках уже стоят другие электроны, и они отталкивают уходящего — как пассажиры в переполненном вагоне.

Новый метод EN-ADC(3) добавляет в расчёт эту «тесноту». Каждая ступенька лестницы получает поправку: если на ней уже есть электрон, подъём стоит дороже. Эта поправка — не одноразовая, а бесконечная цепочка взаимодействий, но считать её можно без дополнительных вычислений. В итоге средняя ошибка для 23 молекул падает с 0.326 эВ до 0.09 эВ — лучше, чем у более дорогого метода EOM-CCSD.

Поправка работает как плата за соседство: два электрона на одной орбитали расталкиваются, и энергия отрыва меняется. Для молекулы BN одна такая поправка сдвигала результат на 0.15 эВ.

Не для всех молекул улучшение одинаково. У молекул с открытыми оболочками, где электроны не спарены, ошибка снижается меньше. А для длинных цепочек метод может даже ухудшить результат. Значит, лестницу нужно достраивать дальше — но для большинства обычных молекул она уже готова.

Сдвиг Эпштейна–Несбета можно интерпретировать как добавление энергии отталкивания двух электронов, оказавшихся в одной и той же орбитали. Для пары виртуальных орбиталей этот сдвиг всегда положителен, что делает знаменатель более отрицательным и предотвращает расходимости при вырождении орбитальных энергий. В расчётах для молекулы BN второй потенциал ионизации менялся на 0.15 эВ только из-за включения или исключения одного члена сдвига — настолько чувствительной оказалась система.

\tilde{t}^{ab(1)}_{ij} = \frac{\langle ab||ij\rangle}{\epsilon_i + \epsilon_j - \epsilon_a - \epsilon_b - \langle ij||ij\rangle - \langle ab||ab\rangle}
Амплитуда перехода из занятых орбиталей i,j в виртуальные a,b с учётом отталкивания двух дырок и двух частиц
L_{pq,pq}(\omega) = \frac{(1-n_p)(1-n_q)}{\omega - \epsilon_p - \epsilon_q - \langle pq||pq\rangle + i\eta} - \frac{n_p n_q}{\omega - \epsilon_p - \epsilon_q + \langle pq||pq\rangle - i\eta}
Пропагатор пары частиц или дырок с учётом их взаимодействия в лестничном приближении
Учёные
Edward WittenJuan Maldacena
Понятия
самосборка, движимая энтропией электронная структура атомные единицы автоэнкодер