Формализм мультипольной плотности - Multipole density formalism

Изоповерхность статической электронной плотности для доксициклин в результате уточнения мультипольной модели на ${ Displaystyle 1.5 mathrm { AA} / е ^ {3}}$ уровень.

Эквивалентная изоповерхность, полученная в результате независимой модели атома. Обратите внимание на отсутствие электронов ковалентных связей.

В Формализм мультипольной плотности (также называемый Формализм Хансена-Коппенса) является Рентгеновская кристаллография метод электронная плотность моделирование, предложенное Нильс К. Хансен и Филип Коппенс в 1978 году. В отличие от обычно используемых Независимая модель атома формализм Хансена-Коппенса представляет асферический подход, позволяющий моделировать электрон распространение вокруг ядро отдельно в разных направлениях и поэтому описывают многочисленные химические особенности молекула внутри ячейка обследованного кристалл в деталях.

Теория

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

Независимая модель атома

Независимая модель атома (сокращенно IAM), на которой основана многополюсная модель, является методом плотность заряда моделирование. Он основан на предположении, что распределение электронов вокруг атом является изотропный, и поэтому плотность заряда зависит только от расстояния до ядра. Выбор радиальная функция Используемая для описания этой электронной плотности является произвольной, при условии, что ее значение в начале координат конечно. На практике либо Гауссовский или же Slater-type 1s-орбитальный используются функции.^[1]

Благодаря упрощенному подходу, этот метод обеспечивает простую модель, не требующую дополнительных параметров (кроме позиционных и Факторы Дебая – Валлера ) для уточнения. Это позволяет IAM удовлетворительно работать при относительно небольшом объеме данных из дифракционный эксперимент доступен. Однако фиксированная форма сингулярного базисная функция предотвращает любое подробное описание асферических атомных элементов.

Каппа формализм

Чтобы скорректировать некоторые валентной оболочки параметрам предложен формализм Каппа.^[2] Он вводит два дополнительных уточняемых параметра: совокупность внешней оболочки (обозначается как ${ displaystyle P_ {val}}$) и его расширение / сжатие (${ displaystyle kappa}$). Следовательно, электронная плотность формулируется как:

${ displaystyle rho _ {atom} = rho _ {core} + rho '_ {valence} ( kappa r) = rho _ {core} + P_ {val} kappa ^ {3} rho _ {val} ( kappa r)}$

Пока ${ displaystyle P_ {val}}$, отвечающая за расходную часть, является линейно связанный с частичная оплата, то нормализованный ${ displaystyle kappa}$ масштаб параметров радиальная координата ${ displaystyle r}$. Поэтому снижение ${ displaystyle kappa}$ Параметр приводит к расширению внешней оболочки и, наоборот, ее поднятие вызывает сжатие.^[3] Хотя формализм Каппа все еще, строго говоря, сферический метод, это важный шаг к пониманию современных подходов, поскольку он позволяет различать химически разные атомы одного и того же элемент.

Многополюсное описание

В описании мультипольной модели плотность заряда вокруг ядра задается следующим уравнением:

${ displaystyle rho _ {atom} (r, theta, phi) = P_ {core} rho _ {core} (r) + P_ {val} kappa ^ {3} rho _ {val} ( kappa r) + sum _ {l = 0} ^ {l_ {max}} kappa '^ {3} R_ {l} ( kappa' r) sum _ {m = -l} ^ {l} P_ {lm pm} Y_ {lm pm} ( theta, phi)}$

Сферическая часть остается практически неотличимой от формализма Каппа, с той лишь разницей, что один параметр соответствует заселенности внутренняя оболочка. Настоящая сила формализма Хансена-Коппенса заключается в правой, деформационной части уравнения. Здесь ${ displaystyle kappa '}$ выполняет роль, аналогичную ${ displaystyle kappa}$ в формализме Каппа (расширение / сжатие асферической части), тогда как отдельные ${ displaystyle R_ {l}}$ фиксированные сферические функции, аналогичные ${ displaystyle rho}$. Сферические гармоники ${ displaystyle Y_ {lm}}$ (каждый со своим популяционным параметром ${ displaystyle P_ {lm}}$), однако, вводятся для моделирования электрически анизотропного распределения заряда.^[4]

В этом подходе фиксированный система координат для каждого атома нужно применить. Хотя на первый взгляд кажется практичным произвольно и без разбора сделать его зависимым от элементарной ячейки для всех присутствующих атомов, гораздо выгоднее назначить каждому атому его собственную местные координаты, что позволяет сосредоточиться на гибридизация -специфический взаимодействия. Хотя единственное сигма-облигация из водород можно хорошо описать с помощью некоторой z-параллели псевдоорбитали, ориентированная в плоскости xy многополюсники с 3-кратным вращательная симметрия окажется более выгодным для квартиры ароматный конструкции.^[5]

Приложения

Карта изоповерхности электронной плотности вокруг ковалентной связи, смоделированная с помощью мультипольной модели, с популяционными параметрами, взятыми из базы данных ELMAM2. Обратите внимание на удлиненную область высокой плотности рядом с атомом водорода, указывающую в направлении кислорода.

Карта изоповерхности электронной плотности вокруг ковалентной связи смоделирована с помощью модели независимого атома в том же масштабе. Ядроцентрированные функции снижают плотность заряда на пути связи.

Основным преимуществом формализма Хансена-Коппенса является его способность освобождать модель от сферических ограничений и гораздо более точно описывать окружение ядра. Таким образом, становится возможным исследовать некоторые молекулярные особенности, которые обычно можно только приблизительно оценить или полностью игнорировать.

Позиционирование водорода

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

Однако атомы водорода обладают уникальной среди всех элементов особенностью - они обладают ровно одним электроном, который, кроме того, находится на их валентной оболочке и, следовательно, участвует в создании сильных ковалентные связи с атомами различных других элементов. Во время формирования связи максимум функции электронной плотности значительно удаляется от ядра в сторону другого атома. Это предотвращает правильное определение положения водорода при любом сферическом приближении. Поэтому обычно положение водорода оценивается на основе данных нейтронной кристаллографии для подобных молекул,^[6] или вообще не моделируется в случае дифракционных данных низкого качества.

Вполне возможно (хотя и спорный) на позицию свободно Усовершенственных атомы водорода с использованием формализма Hansen-Коппенс, после освобождения длины связей от каких-либо ограничений, полученных из нейтронных измерений.^[7] Орбиталь скрепления, смоделированная с помощью адекватных мультиполей, точно описывает распределение плотности, сохраняя при этом правдоподобную длину скрепления. Возможно, стоит приблизить атомы водорода » параметры анизотропного смещения, например с помощью ОТТЕНОК, прежде чем ввести формализм и, возможно, отказаться от ограничений расстояния связи.^[8]

Моделирование склеивания

Чтобы проанализировать длину и силу различных взаимодействий внутри молекулы, Ричарда Бадера "Атомы в молекулах Теорема может быть применена. Благодаря сложному описанию электронного поля, обеспечиваемому этой асферической моделью, становится возможным установить реалистичное пути облигаций между взаимодействующими атомами, а также найти и охарактеризовать их критические точки. Более глубокое понимание этих данных дает полезную информацию о прочность сцепления, тип, полярность или эллиптичность, а по сравнению с другими молекулы позволяет лучше понять реальную электронную структуру исследуемого соединения.^[9]

Поток заряда

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

Недостатки и ограничения

Хотя многополюсный формализм является простым и понятным альтернативным средством уточнения структуры, он определенно не безупречен. Хотя обычно для каждого атома требуется уточнение трех или девяти параметров, в зависимости от того, учитывается ли анизотропное смещение или нет, полное мультипольное описание тяжелых атомов, принадлежащих четвертому и последующим периоды (Такие как хлор, утюг или же бром ) требует уточнения до 37 параметров.^[10] Это оказывается проблематичным для любого кристаллы обладая большим асимметричные блоки (особенно макромолекулярный соединений) и делает уточнение с использованием формализма Хансена-Коппенса недостижимым для данных низкого качества с неудовлетворительным соотношением независимых отражений к уточненным параметрам.

Следует соблюдать осторожность при одновременном уточнении некоторых параметров (т. Е. ${ displaystyle kappa}$ или же ${ displaystyle kappa '}$, многополюсные населенности и тепловые параметры), поскольку они могут сильно коррелировать, что приводит к нестабильному уточнению или нефизический значения параметров. Применение дополнительных ограничений в результате локальная симметрия для каждого атома в молекуле (что уменьшает количество уточненных мультиполей)^[1] или импорт популяционных параметров из существующих баз данных^[11][12] также может быть необходимо для достижения приемлемой модели. С другой стороны, вышеупомянутые подходы значительно сокращают объем информации, необходимой для экспериментов, сохраняя при этом некоторый уровень детализации, касающийся асферического распределения заряда.^[5] Поэтому даже макромолекулярный структуры с удовлетворительными данными дифракции рентгеновских лучей могут быть смоделированы асферически аналогичным образом.^[13]

Несмотря на их сходство, отдельные мультиполи не соответствуют атомным проекциям молекулярных орбиталей волновой функции, как результат квантовые вычисления. Тем не менее, как блестяще резюмировал Стюарт, «структура модельной кристаллической плотности как суперпозиция псевдоатомов [...] действительно имеет количественные характеристики, которые близки ко многим результатам, основанным на квантово-химических расчетах».^[14] Если перекрытие между атомными волновыми функциями достаточно мало, как это происходит, например, в комплексах переходных металлов, атомные мультиполи могут быть коррелированы с атомными валентными орбиталями, а мультиполярные коэффициенты могут быть коррелированы с популяциями металлических d-орбиталей.^[15]Более сильная корреляция между измеренными рентгеновскими лучами дифрагированными интенсивностями и квантово-механическими волновыми функциями возможна с использованием методов, основанных на волновых функциях. ^[16] из Квантовая кристаллография, как, например, рентгеновская модель атомной орбиты,^[17] так называемая экспериментальная волновая функция ^[18] или Hirshfeld Atom Refinement.^[19]

Рекомендации