Прямые методы определения структуры: фазовая проблема

Когда монокристалл облучают рентгеном, детектор фиксирует только интенсивности дифракционных пятен, то есть квадраты модулей структурных факторов $|F_{hkl}|^2$. А чтобы построить карту электронной плотности и увидеть, где сидят атомы, нужны не модули, а полные комплексные структурные факторы - с фазами. Фазы при измерении теряются: это и есть знаменитая фазовая проблема кристаллографии. Прямые методы определения структуры решают её без всяких тяжёлых атомов и изоморфных производных, опираясь только на сами амплитуды и два физических ограничения на плотность. Ниже разберём логику метода по шагам, а собрать конкретную задачу по своему набору рефлексов можно сразу в форме.

Что такое фазовая проблема

Электронная плотность в элементарной ячейке раскладывается в ряд Фурье по структурным факторам:

$$\rho(\mathbf{r}) = \frac{1}{V}\sum_{hkl} |F_{hkl}|\,e^{i\varphi_{hkl}}\,e^{-2\pi i(\mathbf{h}\cdot\mathbf{r})}.$$

Модули $|F_{hkl}|$ мы извлекаем из интенсивностей напрямую, а фазы $\varphi_{hkl}$ - нет. Без них ряд Фурье не собрать: одни и те же амплитуды с разными наборами фаз дают совершенно разные карты плотности. Именно восстановление фаз и есть главная задача структурного анализа, и прямые методы определения структуры берутся за неё статистически, не прибегая к дополнительным экспериментам.

Рассчитать для своей задачи

Два физических ограничения на плотность

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

• Неотрицательность. Электронная плотность нигде не бывает отрицательной: $\rho(\mathbf{r}) \ge 0$. Это резко сужает множество допустимых наборов фаз.
• Атомарность. Плотность сосредоточена в компактных пиках на атомах, между которыми почти ноль. Кристалл - это набор примерно одинаковых разрешённых пиков, а не размазанное облако.

Эти два условия связывают фазы разных рефлексов между собой. Прямые методы переводят такую связь в вероятностные соотношения: зная фазы нескольких сильных отражений, можно с высокой вероятностью предсказать фазу четвёртого.

Уравнение Сейра и нормированные структурные факторы

Дэвид Сейр в 1952 году показал, что для структуры из одинаковых разрешённых атомов структурные факторы подчиняются точному тождеству. В терминах нормированных структурных факторов $E_{hkl}$ (амплитуд, очищенных от спада атомного фактора рассеяния с углом) оно записывается так:

$$E_{\mathbf{h}} = \frac{\theta}{V}\sum_{\mathbf{k}} E_{\mathbf{k}}\,E_{\mathbf{h}-\mathbf{k}}.$$

Нормированный фактор $E_{hkl}$ удобнее обычного $F_{hkl}$: его средний квадрат равен единице, $\langle |E|^2\rangle = 1$, а распределение модулей сразу подсказывает, центросимметрична структура или нет. Сильные рефлексы с большими $|E|$ - те самые, на которых строят фазовые соотношения, потому что для них вероятностные предсказания самые надёжными.

Рассчитать для своей задачи

Триплетное соотношение - сердце метода

Из уравнения Сейра следует ключевое правило. Для трёх рефлексов, индексы которых дают ноль в сумме ($\mathbf{h} + \mathbf{k} + \mathbf{l} = 0$), фазы почти наверняка связаны:

$$\varphi_{\mathbf{h}} + \varphi_{\mathbf{k}} + \varphi_{\mathbf{l}} \approx 0.$$

Это триплетное (или $\Sigma_2$) соотношение. В центросимметричной структуре, где фазы принимают только значения $0$ или $\pi$ (знаки $+$ и $-$), оно превращается в правило знаков: произведение трёх знаков положительно. Надёжность триплета тем выше, чем больше произведение модулей $|E_{\mathbf{h}} E_{\mathbf{k}} E_{\mathbf{l}}|$, и описывается параметром

$$\kappa = \frac{2}{\sqrt{N}}\,|E_{\mathbf{h}} E_{\mathbf{k}} E_{\mathbf{l}}|,$$

где $N$ - число атомов в ячейке. Чем сильнее отражения и чем меньше атомов, тем увереннее триплет. Похожую логику опоры на сильные сигналы используют и метод наименьших квадратов в физических измерениях, где вес результата растёт с надёжностью отсчёта.

Формула касательной и раскрутка фаз

Один триплет даёт оценку фазы с разбросом. Но к каждому рефлексу обычно ведёт много триплетов, и их нужно объединить. Это делает формула касательной Карле и Хауптмана:

$$\tan\varphi_{\mathbf{h}} = \frac{\sum_{\mathbf{k}} \kappa_{\mathbf{h},\mathbf{k}}\,\sin(\varphi_{\mathbf{k}} + \varphi_{\mathbf{h}-\mathbf{k}})}{\sum_{\mathbf{k}} \kappa_{\mathbf{h},\mathbf{k}}\,\cos(\varphi_{\mathbf{k}} + \varphi_{\mathbf{h}-\mathbf{k}})}.$$

Алгоритм запускают с нескольких стартовых фаз (часть назначают произвольно, выбирая начало координат), затем итеративно пересчитывают фазы всех сильных рефлексов по формуле касательной, пока набор не сойдётся. Из сошедшегося набора строят карту $E$-плотности; правильное решение узнают по чётким атомным пикам и по фигуре доверия (figure of merit). Именно за эту схему Джером Карле и Герберт Хауптман получили Нобелевскую премию по химии 1985 года.

Рассчитать для своей задачи

Где прямые методы работают, а где нет

Классические прямые методы определения структуры хороши для малых и средних молекул - примерно до сотни-двух атомов в асимметричной части ячейки. Им нужны данные высокого разрешения (лучше атомного, около 1.2 Å), чтобы атомы разрешались как отдельные пики. Для крупных белков прямой перебор фаз становится слишком ненадёжным, и там работают другие подходы: молекулярное замещение, изоморфное замещение, аномальное рассеяние. Но идея прямых методов живёт и в макромолекулярной кристаллографии - в виде dual-space алгоритмов (например, Shake-and-Bake), которые чередуют уточнение фаз и чистку карты в реальном пространстве.

Распределение Вильсона: центросимметрия по статистике

Перед запуском прямых методов структуру тестируют на наличие центра симметрии, и делают это по статистике распределения нормированных факторов $|E|$. Артур Вильсон показал, что эта статистика различается для двух классов структур. У центросимметричных кристаллов слабых и очень сильных рефлексов больше, среднее значение $\langle |E| \rangle$ около $0.798$, а доля отражений с $|E| > 2$ заметно выше. У нецентросимметричных распределение глаже: $\langle |E| \rangle \approx 0.886$, экстремальных значений меньше. Практический тест строится на среднем квадрате отклонения:

$$\langle |E^2 - 1| \rangle = \begin{cases} 0.968, & \text{центросимметрия}, \\ 0.736, & \text{без центра}. \end{cases}$$

Этот тест задаёт всю дальнейшую тактику: в центросимметричной группе работают со знаками $\pm$, в нецентросимметричной - с непрерывными фазами и полноценной формулой касательной. Ошибка на этом шаге уводит весь расчёт фаз в неверное русло, поэтому распределение $|E|$ всегда смотрят первым делом.

Частые ошибки

• Путать $F$ и $E$. Триплетные соотношения и параметр $\kappa$ работают с нормированными $|E|$, а не с сырыми $|F|$: без снятия спада атомного рассеяния статистика фаз ломается.
• Брать слабые рефлексы. Фазовые соотношения надёжны только для сильных отражений с большим $|E|$; включать в стартовый набор слабые - значит сеять ошибки на первой же итерации.
• Игнорировать симметрию. Выбор начала координат и приписывание стартовых фаз зависят от пространственной группы; неучтённая центросимметрия превращает фазы в знаки и меняет весь алгоритм.
• Ждать решения для белка. Прямые методы в чистом виде не тянут структуры с тысячами атомов и средним разрешением - это не их рабочий диапазон.
• Считать сошедшийся набор истиной. Алгоритм сходится и к ложным минимумам; решение проверяют по фигуре доверия и по виду карты плотности, а не по факту остановки итераций.

FAQ

Чем прямые методы отличаются от метода тяжёлого атома? Метод тяжёлого атома (Паттерсона) восстанавливает фазы по положению одного сильно рассеивающего атома, который надо сначала найти. Прямые методы не требуют тяжёлых атомов вообще: они извлекают фазы из статистических связей между амплитудами всех атомов, поэтому идеальны для лёгких органических молекул.

Почему нужны именно нормированные структурные факторы? Обычный $|F|$ падает с углом дифракции из-за форм-фактора атома, и сильные пики на малых углах маскируют информацию на больших. Нормировка $E = F / \sqrt{\langle F^2\rangle}$ выравнивает вклад всех углов, как если бы атомы были точечными, и только тогда триплетная статистика становится корректной.

Можно ли применять прямые методы к порошковым данным? Можно, но сложнее: в порошковой дифракции рефлексы перекрываются, и извлечь индивидуальные $|F_{hkl}|$ труднее. После аккуратного разделения интенсивностей прямые методы применимы, и для не слишком сложных структур они дают решение.

Коротко

Прямые методы определения структуры решают фазовую проблему рентгеноструктурного анализа, опираясь на неотрицательность и атомарность электронной плотности. Через нормированные структурные факторы $E$, уравнение Сейра и триплетные соотношения они вероятностно связывают фазы сильных рефлексов, а формула касательной итеративно раскручивает весь набор фаз до карты плотности с чёткими атомными пиками. Метод работает для малых и средних молекул высокого разрешения и лежит в основе современных dual-space алгоритмов.