Сканирующая зондовая микроскопия (СЗМ), объединяющая широкий спектр современных методов исследования поверхности, насчитывает более двадцати лет своей истории - с момента создания Биннигом и Рорером сканирующего туннельного микроскопа (СТМ) [1,2]. За прошедшие годы применение зондовой микроскопии позволило достичь уникальных научных результатов в различных областях физики, химии и биологии. Наиболее яркими демонстрациями возможностей этого экспериментального направления при исследовании поверхностей твердых тел могут служить: результаты по прямой визуализации поверхностной реконструкции [3], манипуляция отдельными атомами для записи информации с рекордной плотностью, исследование локального влияния поверхностных дефектов на зонную структуру образца [4] и пр.
Новые возможности рассматриваемого направления в сравнении с традиционными методами исследования поверхности делают особенно перспективным применение зондовой микроскопии (в частности атомно-силовой микроскопии (АСМ) [5]) для изучения биологических и органических материалов. На этом пути в последние годы также был достигнут значительный прогресс. В частности, применительно к исследованиям нуклеиновых кислот, можно упомянуть о таких результатах, как визуализация отдельных молекул ДНК [6] и исследование их конформационного состояния в жидких средах, прямое измерение сил взаимодействия комплементарных нуклеотидов [7], визуализация в реальном масштабе времени процессов взаимодействия ДНК с белками [8].
В то же время, зондовая микроскопия биологических и органических объектов и структур, в том числе полимеров, остается более сложной задачей в сравнении с СЗМ-анализом поверхностей твердых тел. Действительно, прошло более десяти лет с момента возникновения СЗМ, прежде чем была убедительно показана ее адекватность для исследований биообъектов на примере молекулы ДНК [6]. Это связано с такими особенностями подобных объектов, как низкая проводимость (важно при СТМ-исследованиях) и невысокая механическая жесткость. При исследовании микрообъектов актуальной является проблема их иммобилизации на поверхностях твердых подложек в процессах приготовления образцов и при исследованиях (особенно в жидких средах). Важно, чтобы объекты были зафиксированы на подложке в таком состоянии, чтобы было возможно исследовать их интересующие структурные особенности. Возможным подходом к решению этой задачи может служить, например, модификация свойств подложки путем контролируемого осаждения на ее поверхность тонких органических пленок с заданными свойствами.
Весьма важным для адекватного применения зондовых микроскопов в широкомасштабных научных исследованиях является отслеживание и систематизация возможных механизмов возникновения артефактов, т.е. аппаратных эффектов, приводящих к наблюдению ложных или искаженных свойств исследуемого объекта, которое может быть обусловлено, например, воздействием на объект самого инструмента исследования и пр.
Действительно, сканирующий зондовый микроскоп представляет собой ``проектор'', проецирующий объекты и явления микромира на доступный нашему восприятию ``экран'' - в силу многих причин удобно, чтобы им служил экран монитора компьютера. В этом случае проекция становиться отчасти ``осязаемой'', поскольку допускает возможность дополнительного анализа с помощью соответствующего программного обеспечения. Однако подобное ``проецирование'' несет только частичную информацию об объекте, к тому же отчасти искаженную влиянием самого ``проектора''. Восстановление по проекции реальных свойств исследуемых объектов является типичной обратной задачей, требующей решения и для зондовой микроскопии.
Общей чертой всех сканирующих зондовых микроскопов (и определяющей их название) является наличие микроскопического зонда, который приводится в контакт (не всегда речь идет о механическом контакте) с исследуемой поверхностью и, в процессе сканирования, перемещается по некоторому участку поверхности заданного размера.
Контакт зонда и образца подразумевает их взаимодействие. Строго говоря, в общем случае это взаимодействие носит сложный характер. Чтобы осуществлять исследование с помощью конкретного прибора, из широкого спектра выбирается какое-либо одно рабочее взаимодействие. Природа этого выбранного взаимодействия и определяет принадлежность прибора к тому или иному типу в рамках семейства зондовых микроскопов. Информация о поверхности извлекается путем фиксации (при помощи системы обратной связи) или детектирования взаимодействия зонда и образца.
В туннельном микроскопе это взаимодействие проявляется в протекании постоянного тока в туннельном контакте. В основе атомно-силовой микроскопии лежит взаимодействие зонда и образца с силами притяжения или отталкивания. Можно упомянуть о таких разновидностях зондовых микроскопов, как магнитно-силовой микроскоп [9] (зонд и образец взаимодействуют с магнитными силами), микроскоп ближнего поля [10] (оптические свойства образца детектируются через миниатюрную диафрагму, находящуюся в ближней зоне источника фотонов), поляризационный силовой микроскоп [11] (с образцом взаимодействует проводящий заряженный зонд) и т.д.
Миниатюрные размеры зонда и высокая чувствительность детектирующей системы зондовых микроскопов позволяет достигать нано- и субнанометрового пространственного разрешения при детектировании поверхностных свойств (разрешающая способность прибора, как правило, тем выше, чем более короткодействующий характер имеет взаимодействие зонда и образца.)
Процесс сканирования осуществляется при помощи пьезокерамического манипулятора (или системы манипуляторов). Зонд движется последовательно, строка за строкой, вдоль поверхности (изменяются координаты X и Y). Для оцифровки данных участок сканирования разбивается на N строк, а каждая строка на M точек, таким образом, положение иглы в плоскости XY описывается двумя координатами Xi, Yj из множества {Xi, Yj} N×M точек (обычно выбирают N=M). Результатом работы сканирующего зондового микроскопа является установление соответствия между каждой парой координат из множества {Xi, Yj} и некоторым числовым значением (или рядом значений), характеризующим анализируемый параметр поверхности (или ряд параметров).
По способу движения иглы над поверхностью можно провести следующую дифференциацию работы СЗМ.
Таким образом, результатом СЗМ-исследования является получение функциональных зависимостей двух типов: по способу постоянной высоты: F|Z=const(X, Y) и по способу постоянного взаимодействия: Z|A=const(X,Y) (``топография''), плюс какая-либо дополнительная зависимость F|A=const(X,Y). С помощью компьютерного программного обеспечения можно проводить анализ полученных зависимостей (анализ характерных латеральных и вертикальных размеров поверхностных особенностей, построение сечений, фурье-анализ, оценка шероховатости и т.п.), отображать полученные зависимости на экране монитора и выводить их на принтер.
Следует учитывать отличие ``топографического'' массива, полученного в режиме постоянного взаимодействия: {Zij|A=const, Xi, Yj} от реальной топографии поверхности. В случае неоднородного распределения поверхностных свойств, определяющих интенсивность взаимодействия зонда и образца, для извлечения точной информации о топографии объекта необходимо в каждой точке проводить дополнительный анализ взаимодействия зонда и образца.
Например, в туннельной микроскопии реальная геометрия поверхности и карта поверхностного распределения электронных свойств могут быть полностью [12] реконструированы путем анализа трех измеряемых массивов: ``топографии'' в режиме постоянного взаимодействия, первой производной туннельного тока по напряжению смещения и первой производной туннельного тока по величине туннельного зазора.
В сканирующем туннельном микроскопе взаимодействие зонда и поверхности проявляется в протекании постоянного тока в туннельном зазоре между ними. Для плотности туннельного тока (в приближении плоских металлических электродов и вакуумного туннелирования) справедлива формула [13]:
| (1.1) |
| (1.2) |
Из анализа формулы (1.1) следует, что при изменении расстояния зонд-образец на один ангстрем величина туннельного тока изменяется на порядок. Поскольку величина взаимодействия зонд-образец столь существенно зависит от расстояния d, то это позволяет системе обратной связи поддерживать величину d постоянной в процессе сканирования с высокой точностью. Данное обстоятельство обуславливает высокое пространственное разрешение СТМ при определении ``топографической'' функции Z|It=const(X,Y).
Наряду этой зависимостью - ``топографией в режиме постоянного тока'' - в СТМ возможно получение зависимостей типа F|Z=const(X,Y) или F|A=const(X,Y). К первому типу относятся ``токовые'' изображения, полученные в режиме постоянной высоты: It|Z=const(X,Y). Ко второму типу относятся: d ln(It)/d z |It=const(X,Y) и d It/d U|It=const(X,Y), связанные с поверхностным распределением работы выхода j(X,Y) (1.1) и поверхностным распределением плотности электронных состояний r(X,Y,Ef ±eU). Последнее определяется формулой [14]:
| (1.3) |
Получение функциональных зависимостей d ln(It)/d z |It=const(X,Y) и d It/d U|It=const(X,Y) позволяет учесть их вклад в несовпадение Z|It=const(X,Y) с реальной топографией исследуемой поверхности.
В атомно-силовом микроскопе [5] взаимодействие A(X,Y,Z) является силовым взаимодействием зонда и образца.
Характер данного взаимодействия в общем случае достаточно сложен, поскольку определяется свойствами зонда, образца и среды, в которой проводится исследование. В случае исследований незаряженных поверхностей в естественной атмосфере (на воздухе) основной вклад в силовое взаимодействие зонда и образца дают: силы отталкивания, вызванные механическим контактом крайних атомов зонда и образца, силы Ван-дер-Ваальса, а также капиллярные силы, связанные с наличием пленки адсорбата (воды) на поверхности образца.
Зонды для атомно-силового микроскопа имеют форму конуса или пирамиды (рис. 1.1), кончик характеризуется радиусом кривизны R, лежащим, согласно данным фирм-производителей, в диапазоне 5ё40 нм, см. таблицу A.1.
Вычисляя дисперсионное взаимодействие зонда, изображенного на рис. 1.1а, и плоского образца в приближении аддитивности дисперсионного взаимодействия, а также:
| (1.4) |
| (1.5) |
| (1.6) |
Здесь r1 и r2 - плотности материалов зонда и образца, b - константа в законе Лондона [16] для энергии дисперсионного взаимодействия двух нейтральных атомов в вакууме:
| (1.7) |
Для данной константы справедлива оценка b ~ a1a2 Eатом, где a1 и a2 - поляризуемости взаимодействующих атомов, Eатом ~ me4/(h/2p)2 - ``атомная'' энергия (m и e - масса и заряд электрона). Воспользовавшись тем, что a×r - величина безразмерная и для многих элементов имеет порядок 1/10, получим оценку для константы Гамакера: A ~ Eатом/10 ~ 10-19 Дж. Экспериментально измеренные значения константы Гамакера для различных материалов по порядку величины совпадают с выписанной оценкой [17]: Aэксп ~ (0,4ё4)×10-19 Дж.
Проводя оценки для типичных условий АСМ-эксперимента в режиме контакта, получаем для величины ван-дер-ваальсового притяжения: FВдВ ~ 10-8ё10-9 Н.
Стоит отметить, что предположение об аддитивности дисперсионного взаимодействия, сделанное Гамакером при выводе формулы (1.5), вообще говоря, не обосновано. Другие существующие теории (например, подход Лифшица [18], связывающий интенсивность дисперсионного взаимодействия с диэлектрическими проницаемостями материалов взаимодействующих тел и среды исследования) позволяют объяснить экспериментально наблюдаемые в ряде случаев силы дисперсионного отталкивания [19]. Однако, несмотря на то, что природа дисперсионного взаимодействия достаточно понятна, разработать точный подход к описанию явления оказалось трудно, и поэтому современные теории являются приблизительными или модельными; вследствие этого совпадение в реальном исследовании расчетных и экспериментальных данных по порядку величины считается вполне удовлетворительным [20].
Ввиду наличия адсорбированной пленки воды на исследуемой поверхности (в случае ее достаточной гидрофильности), при проведении АСМ-экспериментов на воздухе между зондом и образцом формируется мениск, что приводит к возникновению дополнительной силы притяжения между контактирующими поверхностями [21]:
| (1.8) |
Приведенные значения позволяют по порядку величины оценить величину капиллярных сил в режиме контакта АСМ-исследований: Fкап ~ 10-8 Н.
Сущность этого эффекта - в возникновении интенсивных сил притяжения между двумя гидрофобными поверхностями, находящимися в водной среде. В настоящее время считается, что этот эффект носит сложный характер и при различных условиях может быть обусловлен различными механизмами, проявляющимися, тем не менее, одинаково. Согласно имеющимся представлениям, вклад в эффект могут давать энтропийные процессы, связанные со структурным упорядочиванием молекул воды вблизи гидрофобной поверхности [22], термодинамические процессы, связанные с возрастанием роли субкритических флуктуаций в зазоре [23] и пр.
Деление АСМ по способу измерения и фиксации силового взаимодействия зонда и образца позволяет выделить два основных случая: контактная атомно-силовая микроскопия и АСМ прерывистого контакта.
Для измерения величины силового взаимодействия в контактном режиме используется следующая схема, включающая в качестве миниатюрного динамометра упругую консоль (называемую левером или кантилевером), имеющую на свободном конце зонд (другой конец левера заделан в держателе). При сканировании баланс сил взаимодействия зонда и образца приводит к изгибу левера (т.е. сумма всех сил уравновешивается упругой силой изогнутого левера, см. рис. 1.1б); величина изгиба детектируется прецизионным датчиком. В большинстве атомно-силовых микроскопов для этого используют оптические датчики, реализованные по следующей схеме: луч лазерного диода падает под углом на поверхность левера и отражается в центр четырехсекционного фотодиода (рис. 1.2). При изгибе левера в нормальном направлении или при его кручении возникает разница в сигналах соответствующих участков фотодиода: верхние сегменты/нижние сегменты или правые сегменты/левые сегменты. Первый сигнал несет информацию о балансе сил притяжения и отталкивания, а второй - о латеральных силах взаимодействия зонда и образца (рис. 1.3).
|
|
В процессе сканирования система обратной связи поддерживает на заданном уровне величину изгиба кантилевера (следовательно, и силы воздействия зонда на образец) посредством вариации Z-координаты точки закрепления кантилевера. Сигнал обратной связи несет, таким образом, информацию о топографии поверхности Z|F=const(X,Y).
Упругая сила изгиба кантилевера, действующая на зонд, может быть направлена как в сторону образца, так и в обратную (прямой и обратный изгиб кантилевера, см. рис. 1.1б); в первом случае она увеличивает давление зонда на образец, во втором уменьшает. Как правило, при сканировании стремятся уменьшить воздействие на образец, для чего выбирают минимальную величину прямого изгиба кантилевера, или, по мере возможности, максимальную величину обратного изгиба кантилевера, при которой еще сохраняется механический контакт. Сканирование при обратном изгибе может быть нестабильным. Если силы притяжения (капиллярные, дисперсионные) действуют на зонд неодинаково на всем участке сканирования, то в тех местах, где они меньше, зонд может ``оторваться'' от поверхности (выйти из контакта), если величина упругих сил обратного изгиба превысит силы притяжения. Это обстоятельство является основным препятствием при минимизации силы воздействия зонда.
Разница сигналов правых и левых сегментов фотодиода отображает величину сил трения, действующих на зонд при сканировании, что позволяет исследовать распределение локальных фрикционных свойств поверхности. Информацию о градиенте к исследуемой поверхности несет сигнал отклонений Fоткл (X,Y) детектируемый при сканировании как разностный сигнал верхних и нижних сегментов фотодиода, см. рис. 1.3. Оказывается, что в эксперименте зависимости Fтр(X,Y) и Fоткл(X,Y) часто характеризуются большей латеральной разрешающей способностью, чем топографический сигнал Z|F=const(X,Y), в силу чего оказывается возможным детектирование более мелких деталей поверхности.
Некоторые ограничения. Используемая схема измерения силового взаимодействия в АСМ проста и удобна, но в некоторых случаях силы трения между зондом и образцом могут искажать ``топографический'' сигнал, см. рис. 1.4.
|
В случае, когда направление сканирования совпадает с осью левера, карта поверхностного распределения сил трения будет напрямую накладываться на карту измеряемой ``топографии'' поверхности, поскольку действие сил трения будет варьировать разностный сигнал верхних и нижних сегмента фотодиода (именно этот сигнал фиксируется системой обратной связи).
В случае, когда направление сканирования перпендикулярно оси левера, действие сил трения сведется к вариации крутильного изгиба левера и проявится в разностном сигнале правых и левых сегментов фотодиоа, что позволяет различать вклад нормальных и латеральных сил. Однако в реальности и в этом случае фрикционный сигнал может накладываться на топографический, что может быть обусловлено асимметрией фокусировки лазерного луча относительно оси левера или асимметрией ориентации левера (и зонда) относительно подложки. Эти артефакты могут искажать результаты АСМ-измерения высот объектов. Отследить их возникновение можно, анализируя зависимость измеряемых высот от направления сканирования.
Разрешающая способность АСМ.
При сканировании обратная связь фиксирует разностный сигнал верхних и нижних сегментов фотодиода, нормированный на величину суммарного сигнала всех сегментов фотодиода. Это исключает влияние шумов лазерного диода на точность измерения изгиба кантилевера. Влияние сейсмических шумов в достаточной степени исключается использованием простейших антисейсмических фильтров: например, демпфирующей каучуковой прокладкой под гранитным основанием, на котором устанавливается прибор. Поэтому разрешающая способность атомно-силового микроскопа по нормали (в направлении Z) ограничена другими шумами: пьезоманипулятора, кантилевера и электронного блока (предусилителя, цепи обратной связи и высоковольтных усилителей, задающих сигналы на электродах пьезоманипулятора, см. рис. 1.2).
Критерием разрешающей способности по нормали может служить минимальное изменение Z-координаты иглы при сканировании, детектируемое на уровне шумов. Последний существенно зависит от параметров сканирования (скорости, параметров пропорционального и интегрального звеньев цепи обратной связи, размера кадра), а также от вязкоупругих свойств исследуемого образца. Обычно предел разрешения по нормали составляет доли- единицы ангстрем (в зависимости от параметров эксперимента).
Процедура определения разрешающей способности в латеральном направлении не устоялась. Представляется возможным определить ее следующим образом (рис. 1.5).
Пусть зондирующее острие характеризуется радиусом кривизны R, а разрешаемые особенности поверхности - r (рис. 1.5). Тогда возможность латерального разрешения поверхностных особенностей будет связана с пределом разрешения по нормали Dz: критерием разрешения является условие возможности детектирования разницы в значениях вертикальной координаты иглы над объектами и между ними. Геометрический анализ рис. 1.5 позволяет получить соотношение для минимального расстояния между разрешаемыми поверхностными особенностями, при котором ``провал'' между ними на АСМ-изображении еще может быть детектирован (т.е., когда он равен пределу Dz):
| (1.9) |
|
Поскольку достижимое пространственное разрешение должно являться инвариантной характеристикой прибора (не зависящей от объекта исследования), то его следует определить, рассматривая условие детектирования двух точечных объектов (r=0). Тогда соотношение (1.9) примет вид:
| (1.10) |
Стоит отметить, что введенная процедура определения латерального пространственного разрешения упрощена и не учитывает, например, влияния зонда на структуру объекта (за счет контактных деформаций). Введенная процедура, по-видимому, неприменима при определении разрешающей способности в исследованиях молекулярной (атомной) структуры поверхности (см. раздел 2.1.5).
Отдельно следует подчеркнуть отличие предложенного подхода от процедуры определения разрешающей способности оптических приборов, когда используют критерий Релея, соответствующий тому, что освещенность между изображениями детектируемых точек составляет 74% максимальной. То есть для оптического прибора берут в расчет относительную величину провала между изображениями объектов. В случае АСМ важна абсолютная глубина провала - она должна быть детектируемой на уровне шумов. Это различие является следствием того, что оптические изображения формируются путем сложения интенсивностей изображений отдельных объектов (принцип и), напротив, АСМ-изображение формируется путем исключения вклада в изображение более ``низких'' поверхностных особенностей (принцип или).
Для измерения и фиксации при сканировании интенсивности силового взаимодействия зонда и образца в АСМ прерывистого контакта (tapping mode) используется резонансная схема. Дополнительный пьезоэлемент возбуждает вынужденные колебания левера на его резонансной частоте (вдали от поверхности образца). При сближении зонда и образца возникновение дополнительного градиента сил их взаимодействия приводит к сдвигу резонансной частоты (изменению эффективной жесткости) и частичному выходу системы из резонанса. Наряду с этим, при соударениях зонда и образца увеличивается демпфирование колебаний за счет неупругих процессов. Следствием обоих механизмов является уменьшение амплитуды колебаний.
При сканировании АСМ в режиме прерывистого контакта система обратной связи поддерживает на заданном уровне именно величину амплитуды колебаний: Dz=const. В силу высокой чувствительности амплитуды колебаний к среднему значению расстояния между зондом и образцом, можно получать информацию о топографии поверхности (Z|Dz=const(X,Y)) с достаточно высоким пространственным разрешением.
Информация о тонкой структуре (и локальных вязкоупругих свойствах) исследуемой поверхности может быть получена из измерений зависимости: Dj|Dz=const(X,Y), где Dj - сдвиг фаз между колебаниями левера и внешней вынуждающей силы. Вдали от поверхности колебательная система настраивается в резонанс (j = p/2), однако при сближении зонда и образца она, частично, выходит из резонанса и вклад в сдвиг фаз будут давать упругие (изменение резонансной частоты) и диссипативные (увеличение декремента затухания) механизмы. Оказывается, что в ряде случаев сигнал Dj|Dz=const(X,Y) характеризуется большей латеральной разрешающей способностью, чем топографический Z|Dz=const(X,Y), позволяя разрешить более мелкие детали поверхности.
Стоит отметить, что в настоящее время не существует законченной теории АСМ прерывистого контакта, которая позволяла бы количественно связать параметры эксперимента (величину амплитуды и сдвига фаз колебаний левера) с интенсивностью силового воздействия зонда на образец и с локальными вязкоупругими свойствами образца. Это обстоятельство осложняет количественную интерпретацию фазовых АСМ-изображений (Dj|Dz=const(X,Y)) микроскопии прерывистого контакта.
Метод модуляции силы рассмотрим применительно к атомно-силовой микроскопии контакта. В этом случае общая схема измерений такова. Как и в микроскопии прерывистого контакта в экспериментальную схему введен дополнительный пьезоэлемент (биморф), который возбуждает колебания левера, но в этом случае не на его резонансной частоте, а на резонансной частоте самого биморфа (существенно более массивного). Данное обстоятельство приводит к тому, что, в отличие от АСМ прерывистого контакта, в этом случае силовое взаимодействие зонда и образца не влияет на собственную резонансную частоту системы (система не выходит из резонанса в случае, когда проявляется взаимодействие зонда и образца).
Как и в контактной АСМ в процессе сканирования система обратной связи поддерживает на заданном уровне изгиб кантилевера. Но в этом случае фиксируется среднее за период значение изгиба, поскольку кантилевер совершает вынужденные колебания. Величина амплитуды колебаний будет зависеть от упругих свойств поверхности: более жесткие участки поверхности будут ``продавливаться'' в меньшей степени и там амплитуда колебаний будет меньше.
Как и для микроскопии прерывистого контакта в настоящее время не существует удовлетворительной теории, позволяющей проводить количественный анализ упругих свойств поверхностей на основании анализа АСМ-изображений, полученных в режиме модуляции силы.
Поэтому в литературе по АСМ модуляции силы при интерпретации экспериментальных результатов ограничиваются лишь качественным анализом, привлекая некоторые модельные представления.
В настоящее время в арсенале исследователей, изучающих биологические и органические структуры, имеется широкий спектр различных физических методов исследования. Для того чтобы прояснить возможности и перспективы метода сканирующей зондовой микроскопии и обозначить ее место в ряду других методов исследования, ниже приведена краткая характеристика основных экспериментальных подходов к исследованию биологических и органических структур.
К дифракционным относятся методы реконструкции реальной структуры объекта путем анализа дифракции взаимодействующего с ним пучка электронов, нейтронов или рентгеновский лучей: математическая обработка дифракционной картины позволяет решить обратную задачу и восстановить структуру объекта. Принципиальная трудность метода заключается в невозможности получения непосредственно из эксперимента информации о фазах дифрагировавших лучей. Один из возможных путей решения данной проблемы - использование меток - атомов тяжелых металлов, присоединяемых к определенным группам исследуемого объекта. Так, последовательное мечение различных специфических групп гемоглобина позволило получить необходимую информацию для расчета структуры этого белка [24] по данным рентгеноструктурного анализа его кристаллов. Этот метод нашел применение и при исследовании некристаллических (например, волокнистых) структур. Именно по косвенным рентгеноструктурным данным Уотсон и Крик впервые определили для молекулы ДНК диаметр и параметры спирали (наличие 10 эквивалентных групп в одном обороте, равном 3,4 нм). На основании этих результатов ими и была впервые предложена общепринятая ныне модель строения молекулы ДНК [25].
Дифракционные методы успешно используются и для анализа свойств поверхности, позволяя достигать высокой точности в определении реконструируемых пространственных параметров (при условии достаточной упорядоченности исследуемой структуры). Особенно широко используется метод электронографии, поскольку, из-за сильного взаимодействия электронов с веществом, определяющим в дифракционной картине является вклад именно ближайших к поверхности атомных слоев. В случае рентгеноструктурных исследований поверхности приходится использовать ``скользящую'' геометрию падающего луча и достаточно мощный источник излучения.
Метод рентгеноструктурного анализа применим и для исследований систем частиц, расположение которых не характеризуется ближним или дальним порядком: в этом случае анализируют рассеяние рентгеновских лучей под малыми углами. Этот подход позволяет определять либо форму, либо размеры микрообъектов системы (для решения задачи какая-то информация должна быть известна [26]). Если объекты исследуемого ансамбля характеризуются некоторым распределением по параметру, через который выражается искомый, то точность его восстановления ограничена. Например, при восстановлении формы частиц, разброс их по размерам резко снижает точность решения обратной задачи (если оно вообще может быть найдено) [26], и наоборот. Дополнительная информация, необходимая для корректной постановки обратной задачи, может быть получена с помощью других методов исследования - в том числе сканирующей зондовой микроскопии.
Метод электронной микроскопии (ЭМ), за годы своего существования, был дополнен множеством вспомогательных методик препарирования образцов (при исследовании биообъектов) [27]. Проблема в том, что биообъекты состоят из веществ с малыми атомными номерами, следствием чего является малое число рассеянных электронов и, в итоге, отсутствие разрешения. Сегодня наиболее часто используемые способы контрастирования биообъектов для ЭМ-исследований - это приготовление реплик, оттенение запылением металлами, ``окрашивание'' отдельных частиц и молекул ионами тяжелых металлов. Высокое разрешение дает предложенный Холлом [28] метод негативного контрастирования, принцип которого состоит в помещении исследуемых образцов в пленку аморфного вещества высокой плотности (т.е. хорошо рассеивающего электроны).
При исследовании нуклеиновых кислот методом ЭМ важной проблемой является нанесение исследуемых структур на подложку. При высушивании капли разбавленного раствора молекул на поверхности наиболее часто используемых в ЭМ подложек (коллодий, формвар, углерод) происходит агрегация исследуемых структур. Нанесение молекул ДНК с помощью пульверизатора приводит к их деградации. Поэтому большое распространение получила методика Кляйншмидта и Цана [29], основанная на взаимодействии нуклеиновой кислоты с основными белками и способности этих белков образовывать монослой на поверхности раздела фаз воздух/водный раствор. Молекулы нуклеиновых кислот, связываясь с белковой пленкой, переходят из состояния трехмерного клубка в двумерное расправленное состояние на поверхности водной субфазы. Белковую пленку с молекулами нуклеиновой кислоты переносят на сетку, после чего применяют стандартные методы контрастирования и осуществляют исследование [30]. Недостатком метода является то, что белковая пленка, покрывающая молекулу, увеличивает ее ``толщину'' до 5- 8 нм и затрудняет решение ряда задач (например, исследования взаимодействия нуклеиновых кислот с белками).
Существенной проблемой электронной микроскопии является большая вероятность возникновения артефактов, связанных с необходимостью использования специфических методик контрастирования и фиксирования биообъектов на подложке.
Зондовая микроскопия по общей схеме эксперимента близка ЭМ и для многих объектов позволяет достигать латерального разрешения, близкого к разрешению ЭМ. В то же время, латеральное разрешение СЗМ по нормали выше (метод запыления под углом в ЭМ не дает высокой точности в определении высот объектов из-за конечного размера зерна запыления), что позволяет успешно решать специфические задачи, весьма сложные для ЭМ: например, по увеличению высоты локального участка АСМ-изображения молекулы ДНК, детектировать местоположение триплекса в составе молекулы [31].
Оптические методы исследования получили в последнее время большое распространение (особенно с применением лазерного излучения). Среди наиболее часто используемых методик - спектрофотометрия, основанная на изменении спектра поглощения молекулы при димеризации (в общем случае - полимеризации) за счет взаимодействия соседних мономеров (в основном играет роль диполь-дипольное взаимодействие). Спектрофотометрия позволяет определять концентрации нуклеиновых кислот, определять степень спиральности и исследовать переход клубок ® глобула. Однако результаты исследования чувствительны к наличию в растворе ионов, способных взаимодействовать с исследуемыми молекулами (в том числе протонов, определяющих pH среды). Это связано с возможностью влияния подобного взаимодействия на контролируемые в процессе исследования спектральные изменения. Для выяснения природы спектральных изменений необходимы вспомогательные эксперименты.
Наряду со спектрофотометрией в исследованиях биообъектов широко применяются инфракрасная спектроскопия, спектрополяриметрия, спектрально-флуоресцентные исследования и т.д. Общей чертой данных методик является их косвенный, непрямой (относительный) характер - результатом исследования является усредненный по ансамблю отклик объектов. Проблемой является и необходимость решения обратной задачи для извлечения информации об объекте, что, в ряде случаев, может привести либо к значительной погрешности в результатах, либо к неадекватности решения.
Отдельно следует упомянуть метод флуоресцентной микроскопии, который нашел применение при визуализации одиночных макромолекул в нативных условиях и исследовании изменений их морфологии в процессе взаимодействия с различными лигандами и пр. К недостаткам метода следует отнести невысокую пространственную разрешающую способность (ограниченную пределом оптического микроскопа l/2, поэтому успешно можно исследовать лишь достаточно большие молекулы), и необходимость комплексообразования макромолекулы с молекулами флуоресцентного ``красителя''.
Широко используемые гидродинамические методы исследования биополимеров (вискозиметрия, седиментация, двойное лучепреломление в потоке) позволяют определить приближенную форму, размеры, молекулярный вес и деформируемость исследуемых молекул. Информация, получаемая с помощью гидродинамических методов, также является неабсолютной (усредненной по ансамблю), что подразумевает необходимость сравнения с результатами других методов исследования. Актуальной является проблема корректной постановки и решения обратной задачи для определения параметров объектов.
В этой связи методы сканирующей зондовой микроскопии призваны занять важное место в ряду традиционных исследовательских подходов к изучению биообъектов. Поскольку приборы зондовой микроскопии предназначены для исследования свойств поверхности, то стандартная схема СЗМ-исследования биообъектов включает использование достаточно ровной подложки (проводящей в случае СТМ), на которую осуществляется адсорбция изучаемых структур. Наиболее распространенные в сканирующей зондовой микроскопии подложки - атомарно гладкие сколы слюды, пирографита, дисульфида молибдена, а также эпитаксиально выращенные пленки металлов.
Необходимость использования подложки сближает методы СЗМ с электронной микроскопией. Основные отличия и преимущества сканирующей зондовой микроскопии в сравнении с ЭМ - возможность проведения исследований в различных средах (в том числе в жидкости), отсутствие необходимости какого-либо дополнительного контрастирования исследуемых структур, большая разрешающая способность по нормали, неразрушающий характер исследования (обусловленный невысокой интенсивностью взаимодействия зонд-образец). Т.о. процедура приготовления образцов для СЗМ проще и исследования можно проводить в условиях, близких к нативным.
Все это свидетельствует о том, что зондовая микроскопия является мощным современным методом исследования поверхности, призванным обеспечить дальнейшее развитие традиционным подходам. Результаты исследований последних лет убедительно продемонстрировали, что для решения многих задач метод зондовой микроскопии является самодостаточным (т.е. нет необходимости контрольного сравнения с результатами других исследований); к таковым задачам относятся, например, физическое картирование молекул нуклеиновых кислот [32], исследование изменения их морфологии в жидких средах, анализ их взаимодействия с белками [33,34], анализ структуры тонкопленочных покрытий (см. главу 5) и пр.
Тем не менее, в процессе исследования большое внимание следует уделять возможным механизмам возникновения артефактов. Достаточно очевидным фактом является необходимость учитывать конечность размеров зонда, приводящую к завышению латеральных размеров исследуемых структур (раздел 2.2), а также эффект занижения высот, связанный с контактными деформациями (раздел 2.1). Кроме того, при исследовании многокомпонентных образцов, надлежит проводить дополнительный анализ каждой из компонент, с целью определения их парциального влияния на результат исследования, что позволит избежать ошибок в идентификации комплексов, полученных в результате взаимодействия исходных объектов.
Постановка адекватных контрольных экспериментов, безусловно, способствует уменьшению вероятности ошибки при интерпретации результатов СЗМ-исследований; избежать ошибки позволяет и сравнение с результатами других методов исследования (рассмотренных выше).
Среди множества биологических объектов, к исследованию которых применялись методы зондовой микроскопии, первой была молекула дезоксирибонуклеиновой кислоты. Однако на начальном этапе этих исследований отмечались сложности в интерпретации получаемых изображений. Так, авторы работ [35,36] представили данные сканирующей туннельной микроскопии молекул ДНК, адсорбированных на подложку из высокоориентированного пиролитического графита (пирографита). Выбор данного материала в качестве подложки был обусловлен его проводящими свойствами, а также простотой приготовления атомарно гладких участков поверхности значительной площади. На СТМ-изображениях наблюдались протяженные цепочки, обладающие продольной периодичностью (период около 2-3 нм), на основании чего авторы делали вывод о визуализации витков двойной спирали.
Однако позднее были опубликованы данные [37], свидетельствующие о возможности наблюдения сходных ``ДНК-подобных'' структур при исследовании методом СТМ чистой, свежесколотой поверхности пирографита. Возникновение таких структур объясняется чувствительностью СТМ к электронным свойствам поверхности - происходит визуализация доменных стенок пирографита [38]. Эти результаты, а также слабая воспроизводимость и отсутствие безусловно необходимых контрольных экспериментов в первых работах по СТМ-визуализации ДНК вызывали сомнения в применимости методов зондовой микроскопии для исследования макромолекул.
Существенное изменение ситуации произошло в 1992 г., когда были опубликованы первые надежные и воспроизводимые результаты исследования ДНК методом атомно-силовой микроскопии (АСМ) [6,39]. Таким образом, прошло более пяти лет с момента создания атомно-силового микроскопа (1986 г.) и более десяти лет существования зондовой микроскопии, прежде чем был достигнут определенный успех в применении метода к исследованию биологических объектов. В то же время методы зондовой микроскопии весьма успешно применялись в течение этих десяти лет для исследований структуры поверхности твердых (кристаллических) тел. За эти годы возникло понимание, что потенциал возможностей зондовой микроскопии в исследованиях биологических объектов может быть реализован лишь при решении в каждом конкретном случае основной экспериментальной задачи: определения адекватной методики препарирования образцов.
Универсальной методики приготовления образцов для решения широкого спектра задач зондовой микроскопии пока не существует, поэтому используемые методики определяются спецификой задачи, т.е. являются методиками ``на конкретный случай''. Применительно к исследованиям нуклеиновых кислот во многих случаях бывает важно, чтобы макромолекулы адсорбировались на подложку в расправленном состоянии. Для успешного исследования необходимо также, чтобы исследуемые молекулы достаточно прочно фиксировались на подложке. В противном случае, при исследованиях в жидкой среде, молекулы будут увлекаться зондом в процессе сканирования; при исследовании на воздухе, силы, действующие на молекулу при высыхании капли препарата на подложке, также будут увлекать адсорбируемые структуры, приводя к неоднородному (с образованием агрегатов) поверхностному распределению исследуемых объектов. Ниже рассмотрим основные методики, позволяющие зафиксировать молекулы нуклеиновых кислот на поверхности в расправленном состоянии.
Первые надежные результаты исследования ДНК методом атомно-силовой микроскопии были получены в том случае, когда изучаемые структуры наносили из капли рабочего раствора (раствора буфера, водно-спиртовой среды и т.п.) на поверхность слюды, модифицированной ионами двух- (и более) валентных металлов [6,39]. Эти ионы, по-видимому, служат связующими мостиками между отрицательно заряженной слюдой (в водных растворах) и отрицательно заряженными фосфатными группами молекулы ДНК. Модификация слюды может предшествовать процессу адсорбции макромолекул - в этом случае свежесколотую слюду помещают для предварительной обработки на некоторое время (минуты, часы) в раствор, содержащий катионы металлов, затем промывают (водой, водным раствором этанола и пр.) и высушивают. После этого на модифицированную поверхность наносят рабочий раствор с исследуемыми структурами. Другой подход - добавление солей металлов непосредственно к рабочему раствору - позволяет исключить этап предварительной обработки слюды (препарат наносят непосредственно на поверхность свежего скола). После высыхания капли рабочего раствора образцы, как правило, подвергают дополнительной промывке.
В работе [40] методом АСМ проводили сравнительный анализ процесса адсорбции макромолекул ДНК на слюду из раствора, как с катионами металлов, так и без них. Было показано, что во втором случае адсорбированные молекулы нестабильны в процессе промывки образцов и при сканировании. Кроме того, АСМ-изображения макромолекул характеризуются большим количеством запутанных и перекрученных участков, тогда как при присутствии в растворе катионов металлов и применении предшествующей сканированию процедуры промывки образцов молекулы ДНК фиксируются на поверхности подложки в расправленном состоянии.
Авторы работы [41] исследовали методом АСМ в жидкостной ячейке в реальном масштабе времени обратимое осаждение молекул ДНК на слюду при различных параметрах рабочего раствора. Было показано, что адсорбция объектов на поверхность подложки и величина силы адгезии существенно зависят от концентрации связующих катионов (оптимальное значение около 2 мМ для Zn 2+) и pH среды (оптимальное значение около 7,5).
На сегодняшний день использование катионов металлов для связывания ДНК с подложкой является, пожалуй, наиболее распространенной методикой приготовления образцов при АСМ-исследованиях, что объясняется ее простотой и высокой воспроизводимостью результатов.
На рис. 1.6 представлены АСМ-изображения молекул ДНК, адсорбированных на поверхность слюды в присутствии катионов Ni 2+ (а) или уранилацетата (б). Из рисунка видно, что молекулы зафиксированы на подложке в достаточно расправленном состоянии.
Приводимые изображения молекул характеризуются завышенными латеральными размерами - 10-15 нм - и заниженными вертикальными - около 0,7 нм, в то время как диаметр одиночной молекулы ДНК около 2 нм. Завышение латеральных размеров полученных изображений обусловлено конечной величиной радиуса кривизны зонда (количественная методика учета данного артефакта и определения реальной ширины объекта разработана в разделе 2.2). Причиной занижения вертикальных размеров может являться деформация молекулы под действием зонда (рассмотрение теории контактных деформаций и применение ее для описания эффекта занижения высот АСМ-профиля объекта приведено в разделе 2.1). Следует подчеркнуть, что эти артефакты являются типичными при АСМ-исследованиях молекул нуклеиновых кислот [42].
Другая широко используемая методика химической модификации подложки для стабилизации молекул нуклеиновых кислот на поверхности - это силанизация слюды [43]. Процесс силанизации существенно не изменяет топографических особенностей поверхности слюды - она остается достаточно гладкой, в то же время улучшается связь макромолекулы с модифицированной подложкой. В работе [44] при проведении АСМ-исследований нуклеиновых кислот (ДНК, двунитевой рибонуклеиновой кислоты (РНК)) в качестве подложки использовали поверхность слюды, химически модифицированную 3-аминопропилтриэтоксисиланом (APTES). Модификация включала ковалентную привязку аминогрупп молекул APTES к слюде, степень модификации контролировали с помощью АСМ, и выбирали таким образом, чтобы шероховатость получаемой поверхности не затрудняла идентификацию исследуемых структур. Это обстоятельство (необходимость контроля степени шероховатости модифицированной подложки) незначительно усложняет методику.
Было показано [45], что этот подход препарирования образцов достаточно эффективен в рутинных исследованиях распределения молекул нуклеиновых кислот по длинам. Представленные авторами результаты совпадали с данными электронной микроскопии (исследовались молекулы двунитевой РНК ретровируса), при этом было отмечено, что при равенстве разрешающих способностей важным преимуществом АСМ является существенно менее сложная методика приготовления образцов.
Была продемонстрирована возможность использования модифицированной APTES слюды при проведении исследований как на воздухе, так и в водной среде - в обоих случаях величина адгезии молекул нуклеиновых кислот к подложке позволяла осуществлять сканирование в контактном режиме. В работе [46] осуществляли сканирование одного и того же участка поверхности с адсорбированными макромолекулами на воздухе и в воде, при этом во втором случае отмечалось увеличение разрешающей способности (примерно в три раза), что может быть связано с исключением негативного влияния капиллярных сил (см. раздел 1.3.1). Капиллярные силы увеличивают силовое воздействие зонда на образец, что вызывает дополнительные деформации макромолекулы, ее разрушение или нестабильность в процессе сканирования, и приводит к уменьшению разрешающей способности.
Другие исследователи также сообщают, что, если АСМ-исследования проводятся либо в сухой газовой атмосфере (иногда для уменьшения относительной влажности температуру рабочей атмосферы повышают до 60-100°С [47]), либо в жидкостной ячейке [48,49], то влияние капиллярных сил уменьшается, что увеличивает достигаемое пространственное разрешение. Применение режима прерывистого контакта (tapping mode, см. раздел 1.3.2) [50,51] также позволяет исключить влияние капиллярных сил; микроскопия прерывистого контакта позволяет сегодня достичь разрешения (при проведении исследований в жидких средах), необходимого для визуализации витков двойной спирали ДНК [52].
На пути подбора адекватного подхода к приготовлению образцов может оказаться полезным накопленный за десятилетия исследований багаж методик препарирования образцов электронной микроскопии. В упрощенном варианте многие из этих методик могут применяться и для решения задач зондовой микроскопии. Так, модифицированный метод Кляйншмидта [29,30] (метод белковой пленки, см. раздел 1.4) дает хорошие результаты при использовании гидрофобных подложек [53]. Недостатком метода может служить то, что комплексообразование молекул нуклеиновых кислот с белковой пленкой (часто используют пленку цитохромов C) затрудняет исследование взаимодействия этих молекул с другими белками, что ограничивает круг решаемых задач.
Для фиксации развернутых молекул ДНК на поверхности подложки применяют также бензилдиметилалкиламмоний хлорид (BAC) [54,55] - широко используемый в электронной микроскопии реагент, представляющий собой катионное поверхностно-активное вещество (ПАВ). BAC в малых концентрациях (около 5×10-5%) добавляют непосредственно в раствор, содержащий макромолекулы, перед нанесением препарата на поверхность слюды. Вследствие свойств ПАВ, на поверхности капельки, наносимой на подложку, формируется стабильная пленка молекул BAC (со связанными с ней молекулами ДНК), которая распределяется по поверхности слюды и стабилизируется за счет сил электростатического взаимодействия. Использование других ПАВ (в малых концентрациях) также способствует стабилизации молекул нуклеиновых кислот на поверхности слюды. Так, в работе [56] для разворачивания ДНК на поверхности, успешно применили 2,4,6-трис(диметиламинометил)фенол (DMP-30) и хлорид цетилпиридиния (CP), представляющие собой, соответственно, неионное и катионное ПАВ. На рис. 4.1 представлен результат АСМ-исследования ДНК бактериофага Т4, адсорбированной на поверхность слюды из раствора, содержащего малые концентрации цетилтриметиламмоний бромида (ЦТАБ).
К недостаткам метода с применением ПАВ можно отнести некоторую поверхностную неоднородность получаемых образцов, что обусловлено неконтролируемостью процессов формирования пленки и комплексообразования с исследуемыми макромолекулами. Удачным дополнением метода может служить использование методики Ленгмюра-Блоджетт: формирование на границе раздела фаз жидкость/газ пленки амфифильных молекул и последующее контролируемое перенесение ее на поверхность твердой подложки (см. главу 5).
В работе [57] для иммобилизации молекул ДНК применили другую методику, сходную с традиционной для электронной микроскопии. Макромолекулы адсорбировали на поверхность слюды, сверху напыляли слой углерода. После этого на углерод приклеивали стальную пластинку, слюду удаляли, и открытой для сканирования оказывалась внутренняя (обращенная прежде к слюде) поверхность углеродной пленки с жестко закрепленными в ней макромолекулами. Шероховатость поверхности углеродной пленки составляла, по оценкам авторов, единицы ангстрем и допускала однозначную идентификацию отдельных молекул ДНК. Авторы работы отмечают высокую стабильность приготовленных образцов - макромолекулы не разрушались при больших силах воздействия зонда, результаты не зависели от среды исследования (эксперименты проводились как на воздухе, так и в жидкости), длительное хранение не влияло на качество образцов. Недостаток этого метода состоит в его относительной экспериментальной трудоемкости.
Успехи применения атомно-силовой микроскопии к исследованию молекул нуклеиновых кислот позволили достичь прогресса и при проведении подобных экспериментов методами сканирующей туннельной микроскопии. Хорошие результаты были получены при СТМ-исследовании молекул ДНК, адсорбированных на химически модифицированных поверхностях металлов [58]. Химическая модификация включала ковалентную связь химически поляризуемых групп тиолов с чистой поверхностью металла. Адсорбция ДНК осуществлялась за счет кулоновского взаимодействия молекулы с плотно упакованной мономолекулярной пленкой тиолов, ориентированной положительными функциональными группами к молекуле ДНК. Характерной особенностью представленных в работе [59] изображений молекул ДНК (приготовленных по описываемой методике) является отрицательный контраст (при движении над молекулой туннельная игла опускается ниже, чем при движении над поверхностью подложки), связанный, по-видимому, со слабой проводимость макромолекулы.
Другой перспективный подход к исследованию макромолекул методами СТМ - применение низкотоковой туннельной микроскопии [60] с рабочим диапазоном туннельного тока менее 1 пА. В этом случае экспериментальные результаты [61] свидетельствуют о возможности визуализации молекул ДНК, адсорбированных на поверхности диэлектрической подложки (слюды). Транспорт заряда в системе зонд-образец-подложка осуществляется за счет проводимости ультратонкой пленки воды, покрывающей поверхности образца и подложки при проведении исследований в условиях контролируемой относительной влажности. Оптимальное значение последней для подобных экспериментов составляет 60- 65%.
Т.о. спектр применяемых методик препарирования образцов нуклеиновых кислот для СЗМ весьма широк. В то же время стоит отметить, что описанные методики разработаны для исследования свободных молекул ДНК, поэтому априори неизвестно адекватны ли они для изучения, например, взаимодействия нуклеиновых кислот с белками или поверхностно-активными веществами. Упомянем также, что в литературе по СЗМ описаны, в основном, лишь исследования двунитевых молекул нуклеиновых кислот: РНК и ДНК, имеющих двойную спираль, и практически отсутствуют данные об исследовании однонитевых молекул (которые сложнее зафиксировать на подложке в расправленном состоянии), см. главу 3.
При подборе адекватной методики приготовления образцов сканирующий зондовый микроскоп представляет собой удобный и надежный прибор для исследования свойств биологических структур на молекулярном уровне. Убедительным доказательством этого утверждения могут служить результаты работ, посвященных исследованию конформационных свойств молекул нуклеиновых кислот и их комплексов с белками, ПАВ и т.п.
Авторы работы [33] исследовали влияние специфических лигандов (дистамицина и микрогонотропена-6b) на конформационные свойства нескольких типов молекул ДНК методом АСМ прерывистого контакта. Результаты исследований не только обеспечили визуальное доказательство влияния данных лигандов на конформацию ДНК, но и позволили сделать определенные количественные оценки, основанные на анализе представленных авторами гистограмм распределения расстояния между концами ДНК.
В работе [34] исследовали специфическое взаимодействие РНК-полимеразы из Escherihia coli с ДНК, содержащей lPL промотор. Образованные комплексы были визуализированы методом АСМ, на полученных изображениях наблюдались молекулы ДНК с РНК-полимеразой, локализованной на расстоянии 4/9 от одного из концов. Проводили анализ структуры комплекса РНК-полимеразы и ДНК как в открытых комплексах, связанных с промотором (OPCs), так и в стабильных элонгационных комплексах (С15) с образованным 15-нуклеотидным транскриптом. Было показано, что в области прикрепления РНК-полимеразы происходит изгиб молекулы ДНК, причем средний угол изгиба для OPCs и С15 составлял 54±31° и 92±37°, соответственно. На основании обнаруженного различия в значениях средних углов изгиба для OPCs и C15 авторы делают вывод, что процесс транскрипции сопровождается изменением структуры комплекса РНК-полимеразы с ДНК, обусловленной конформационными изменениями полимеразы.
Возможность АСМ-визуализации молекулярных процессов продемонстрирована в работе [8] на основе последовательного анализа образования неспецифичных комплексов молекулы ДНК и РНК-полимеразы в реальном масштабе времени. Схема эксперимента включала в себя адсорбцию макромолекул ДНК на слюду из 10 мM HEPES буфера в присутствии ионов магния (1-10 мM MgCl2), затем следовали промывка и высушивание образцов в эксикаторе. После этого образцы помещали в жидкостную ячейку микроскопа, начинали процесс сканирования, а затем в раствор вводили РНК-полимеразу. Было показано, что образование комплексов белок-ДНК наблюдается уже через несколько секунд после добавления полимеразы, что свидетельствует о сохранении нативной конформации РНК-полимеразы и ДНК в процессе приготовления образцов и проведения исследований. Эти результаты демонстрируют возможность использования зондовой микроскопии для исследования процессов, ответственных за распознавание и сборку макромолекулярных комплексов в физиологических условиях.
Большое прикладное значение имеет возможность решения методом АСМ задачи физического картирования ДНК: специфические участки макромолекулы помечаются определенными маркерами, анализ местоположения этих маркеров на получаемых АСМ-изображениях позволяет составлять карты относительного расположения специфических последовательностей нуклеотидов. Так, в работе [42] продемонстрирована возможность визуализации местоположения биотина, ковалентно привязанного к первому нуклеотиду праймера. Биотин помечался белковым комплексом стрептавидин-стафиллококковый белок А (стрептавидин имеет высокую способность связывания с биотином, а белок А увеличивает размеры маркера для однозначной его идентификации на АСМ-изображениях).
В работах [62,32] был использован метод картирования клонированных в плазмидный вектор последовательностей LTR (длинные концевые повторы) с помощью специфических маркеров, имеющих характерную форму - R-петель. R-петли формировались последовательностями РНК из 345 и 380 нуклеотидных оснований, комплементарными к U3 и U5 областям LTR человеческого эндогенного ретровируса K-10 (HERV-K10). В процессе образования петли РНК формировала двойную спираль с комплементарным участком ДНК, вытесняя вторую нить ДНК, которая коллапсировала в результате воздействия ионов Mg 2+; характерная форма образованной структуры (петля) позволяла уверенной идентифицировать ее на АСМ-изображениях. Различная длина зондов позволяла уже из одной гистограммы определить как положение, так и ориентацию LTR. Полученные результаты позволяют предположить, что в будущем разрешение приборов СЗМ может быть достаточным для определения последовательности нуклеотидов ДНК.
Уникальную возможность зондового микроскопа как прибора, позволяющего проводить прецизионные исследования локальных свойств поверхности, продемонстрировали авторы работы [7]. Они проводили прямые исследования силового взаимодействия, ответственного за формирование витков молекулы ДНК, по следующей схеме.
На поверхностях подложки и кремниевого микрозонда создавали два типа покрытия, представляющего собой слой ковалентно привязанных за один из концов комплементарных олигомеров - в одном случае (АЦТГ)5, в другом (ЦАГТ)5. В процессе взаимодействия пары данных олигонуклеотидов (длиной в 20 нуклеотидов) возможно образование комплексов с 20, 16, 12, 8 и 4 парами взаимодействующих оснований.
В ходе эксперимента многократно измеряли кривую силового взаимодействия F(z) между зондом и подложкой. На основании этих результатов определяли силу адгезии и строили гистограммы ее распределения. В случае отсутствия специфического взаимодействия между комплементарными участками наблюдалось бы однородное гауссово распределение для силы адгезии. Однако на полученных гистограммах четко прослеживались отдельные пики (1,52±0,19; 1,1±0,13; 0,83±0,11 нН), отображающие, очевидно, специфическое взаимодействие 20, 16 и 12 пар комплементарных оснований (авторы указывают, что комплексы с 8 и 4 взаимодействующими парами оснований термодинамически нестабильны при 27°С - температуре проведения исследований).
Сходную экспериментальную схему - с образованием между взаимодействующими поверхностями ``мостиков'' в виде отдельных нитей ДНК длиной в 160 оснований - использовали для анализа внутримолекулярных упругих свойств. Таким образом, авторам удалось продемонстрировать возможность применения АСМ для прямого количественного анализа межмолекулярного и внутримолекулярного взаимодействия в комплексах биологических и синтетических макромолекул, что открывает новые перспективы для обнаружения и локализации специфических последовательностей нуклеотидов с ангстремным разрешением.
Несмотря на возможность достижения высокого пространственного разрешения, информация, получаемая методами зондовой микроскопии (в частности - АСМ), может неадекватно отображать реальные особенности поверхности, что является следствием влияния инструмента исследования на объект и приводит к наблюдению артефактов. Эти артефакты, как правило, легко учитываются на качественном уровне при интерпретации АСМ-результатов, однако специфика ряда задач может потребовать количественных оценок и методов восстановления реальной геометрии объектов.
Ниже проанализированы два основных артефакта АСМ, влияние которых существенно при проведении исследований отдельных микрообъектов, адсорбированных на поверхность твердой подложки: эффекта уширения профиля и эффекта занижения высот АСМ-изображения объектов исследования. Построены количественные методики учета влияния рассматриваемых эффектов на результаты исследования АСМ. В свете результатов теории контактных деформаций рассмотрен механизм визуализации атомной (или молекулярной) структуры поверхности при исследованиях АСМ.
С первых работ по АСМ-визуализации молекул нуклеиновых кислот [6,39] отмечалось, что высоты АСМ-изображений ДНК существенно занижены в сравнении с модельными представлениями о структуре молекулы. В то же время, для ряда других объектов (с близкими физическими свойствами, но отличными геометрическими размерами) эффект занижения высот проявляется не столь выраженно. Так было обнаружено (см. раздел 3.1), что эффект занижения высот для вирусных частиц несущественен (при минимизации силы воздействия зонда), в то время, как высота АСМ-изображений молекул нуклеиновых кислот в тех же условиях эксперимента занижена более чем на 50%, несмотря на то, что и те и другие объекты визуализированы на одном кадре. Применение излагаемой ниже методики позволило объяснить данную закономерность и связать ее с различием радиусов частиц ВТМ ( ~ 10 нм) и нуклеиновых кислот ( Ј 1 нм).
Следуя [63] предположим, что эффект занижения высот АСМ-изображений объектов связан с контактными деформациями. Действительно, в процессе сканирования, как в режиме контакта, так и прерывистого контакта, зонд воздействует на образец с некоторой силой (см. раздел 1.3), величина которой обычно лежит в диапазоне (1ё100)×10-9 Н. Действие этой силы, в связи с чрезвычайно малым радиусом кривизны зондирующего острия (обычно около 10 нм), приводит к возникновению значительного контактного давления, которое должно вызывать контактные деформации.
Впервые задача о контактных деформациях двух тел была рассмотрена Г. Герцем в 1882 г. [64]. Остановимся на основных соотношениях решения контактной задачи Герца [65]. Рассмотрим два тела, которые соприкасаются друг с другом в некоторой точке их поверхностей, см. рис. 2.1 а, на рисунке изображен разрез соприкасающихся поверхностей, O - точка контакта.
|
Вблизи точки касания O уравнения поверхностей могут быть записаны в виде:
| (2.1) |
Если теперь эти два тела сдавливаются некоторой силой F, то они будут деформироваться и сблизятся на некоторое расстояние h (рис. 2.1 б), при этом областью соприкосновения будет уже не одна точка, а некоторый участок конечной площади S.
Дальнейший анализ, проведенный Герцем, включает рассмотрение тензора cab + cўab, главные значения которого A и B могут быть выражены через R1, R2, Rў1, Rў2, соответствующие формулы для общего случая приведены в [65]:
| (2.2) |
Известные соотношения Герца показывают, что формой области контакта является эллипс с полуосями a и b, причем (при условии малости деформаций в сравнении с соответствующими радиусами кривизны):
| (2.3) |
| (2.4) |
Воспроизведем, также, формулу, описывающую распределения давления в области контакта:
| (2.5) |
Соотношения Герца (2.3) связывают искомые параметры задачи (полуоси a, b эллиптической области контакта и величину сближения h) с известными параметрами (сдавливающей силой F, модулями Юнга зонда и образца, параметрами A и B) через интегральные зависимости, сводимые к эллиптическим интегралам. Однако в силу того, что (2.3) являются системой нелинейных уравнений с неявными зависимостями от искомых параметров a и b, то, для возможности их применения при интерпретации экспериментальных результатов, необходимы либо реализация общего численного решения, либо дополнительный анализ с привлечением упрощающих предпосылок.
Известно, что соотношения Герца (2.3) легко упрощаются в случае, когда контактирующие тела имеют сферическую геометрию (рассмотрение этого случая входит в стандартный курс теории упругости, см., например, [65]).
Такая модель геометрии контакта может найти широкое применение в задачах, связанных с интерпретацией результатов СЗМ - при рассмотрении контактных деформаций, возникающих в процессе сканирования микрообъектов, форма которых может быть аппроксимирована сферой (например, молекул ряда белков и пр.), а также плоских образцов, например, тонких пленок.
Воспроизведем известные формулы [65] и проведем оценки. Если зонд и образец вблизи точки контакта описываются сферическими поверхностями и характеризуются радиусами кривизны R и Rў, то, согласно (2.2),
|
| (2.6) |
Для величины h - сближения зонда и образца за счет контактных деформаций - в этом случае справедлива формула:
| (2.7) |
И приведем формулу для вычисления давления на образец:
| (2.8) |
Рассмотрим контакт кремниевого зонда (E = 1,5×1011 Па) радиусом кривизны R = 10 нм (типичное значение, см. таблицу A.1) c плоским образцом, и проведем оценки для двух значений величины сдавливающей силы: F = 5×10-9 Н и F = 5×10-8 Н (таблица 2.1).
Модуль упругости образца, Па | a, нм - радиус области контакта | h, нм - сближение за счет деформации | P, ГПа - контактное давление | |||
108 | 7,2 | 16 | 5,2 | 24 | 0,03 | 0,07 |
109 | 3,4 | 7,2 | 1,1 | 5,2 | 0,14 | 0,3 |
1010 | 1,6 | 3,4 | 0,25 | 1,1 | 0,63 | 1,4 |
1011 | 0,9 | 1,8 | 0,07 | 0,3 | 2,2 | 4,7 |
1012 | 0,7 | 1,4 | 0,04 | 0,2 | 3,7 | 7,9 |
при величине сближающей силы (нН): | ||||||
5 | 50 | 5 | 50 | 5 | 50 |
Таблица перекрывает диапазон значений модулей упругости образцов от 108 Па (значение модуля упругости пенопласта) до 1012 Па (значение модуля упругости алмаза). Для многих материалов значения модуля Юнга лежат в диапазоне 1010ё1011 Па; для биополимеров - как правило, в диапазоне 109ё1010 Па. Один из существенных выводов анализа таблицы тот, что величина контактного давления составляет значительную величину. Из анализа таблицы можно предположить, например, что АСМ-исследования очень твердых материалов (типа алмаза) при использовании обычных кремниевых зондов могут приводить к разрушению последних (предел прочности кремния: sв = 0,7×109 Па), если не уделять внимание минимизации силы взаимодействия зонда и образца.
Тот результат, что для более жестких образцов выше значение контактного давления при АСМ-исследованиях, служит объяснением следующей экспериментально наблюдаемой закономерности. Поверхность слюды можно локально разрушить при сканировании участка небольшой площади (100×100 нм2), если увеличить силу воздействия зонда до величины ~ 10-7 нН. В то же время, в тождественных условиях эксперимента поверхность пирографита разрушить не удается. Это объясняется меньшим значением контактного давления во втором случае в силу меньшей жесткости пирографита (значение микротвердости пирографита Hm ~ 1010 Па, слюды Hm ~ 2×1011 Па).
С другой стороны, согласно таблице 2.1, типичные значения контактного давления в АСМ-эксперименте часто превышают пределы прочности исследуемых образцов и при меньших силах воздействия зонда (пределы прочности (по сжатию) слюды и графита составляют sв.сж @ 0,5×109 Па и sв.сж @ 0,03×109 Па соответственно [66]) .
В чем же тогда причина возможности проведения методом АСМ неразрушающих исследований широкого спектра материалов?
Вот аргументы, объясняющие возможность проведения неразрушающих исследований методом АСМ:
Действительно, во-первых, значения пределов прочности определяют приложением разрушающего давления к некоторой поверхности образца в заданном направлении (одном). В нашем случае мы имеем дело с приложением к образцу локального давления, которое перераспределяется по трем направлениям.
Во-вторых, при сканировании зонд оказывает воздействие на локальный участок образца в течение промежутка времени t ~ a/Vскан ~ a/Lf, где a - латеральный размер области контакта (2.6), L - длина скана (строки сканирования), f - частота строчной развертки, Vскан - скорость сканирования. Приведем оценку границ диапазона характерного времени взаимодействия. Согласно таблице 2.1, типичные значения a ~ 1 нм. Для частоты сканирования 10 Гц при размере кадра 15 мкм имеем для времени взаимодействия оценку t ~ 10-5 сек, а при размере кадра 100 нм - t ~ 10-3 сек. Полученные оценки и определяют границы диапазона типичных значений времен взаимодействия зонда и локального участка образца при сканировании.
Поскольку известен экспериментальный факт, что в результате замедления скорости сканирования (увеличения t от значения 10-5 сек до 10-3 сек при достаточной величине силы воздействия зонда), ряд исследуемых поверхностей может разрушаться, то можно предположить, что характерные времена процессов разрушения этих поверхностей под локальным воздействием зонда попадают в указанный диапазон.
Эффект разрушения исследуемых поверхностей при замедлении скорости сканирования.
Замедление скорости сканирования (при фиксированной величине воздействующей силы) может приводить к разрушению исследуемой поверхности. Например, уменьшение размера кадра (при той же частоте строчной развертки приводит к уменьшению скорости сканирования (перемещения зонда по поверхности) и может вызывать разрушение поверхности, успешно сканируемой при больших размерах кадра при том же значении силы воздействия зонда.
На этом принципе основана методика формирования в тонких органических пленках искусственных дефектов: участок поверхности заданной площади сканируют при медленной скорости, что приводит к локальному разрушению пленки зондом и удалению ее материала с этого участка; затем, при сканировании кадра большего размера (при увеличенной скорости сканирования), можно визуализировать искусственный дефект в пленке, размеры и форма которого совпадают с участком предварительного сканирования (см. рис. 5.10). С другой стороны, стоит отметить, что именно разрушение поверхности при уменьшении размера кадра (замедлении скорости перемещения зонда по поверхности) является основным препятствием при исследованиях молекулярной упаковки тонких пленок, см. раздел 5.3.
Можно было бы предположить, что механизм эффекта следующий. При большой скорости сканирования объект (например, тонкая пленка) не ``успевает'' деформироваться и эффективная глубина проникновения зонда меньше, чем определяемая формулой (2.7). При замедлении скорости перемещения зонд глубже ``вязнет'' в образце и, если глубина проникновения h становиться сравнимой с толщиной пленки, то имеет место разрушение.
Однако, по-видимому, это объяснение не может быть принято как универсальное. Действительно, согласно оценками раздела 2.1.4, характерное время установления упругих контактных деформаций составляет величину ~ 10-6 сек, что меньше или много меньше, чем время взаимодействия зонда и локального участка образца (10-3ё10-5 сек, см. выше). Поэтому применение статических формул теории контактных деформаций (типа (2.7)) для анализа результатов АСМ является оправданным, в силу чего глубину проникновения зонда в образец за счет упругих деформаций следует считать не зависящей от скорости сканирования. Т.о. проводимые в разделе 2.1.4 оценки позволяют сделать вывод, что наблюдаемый эффект: разрушения поверхности при замедлении скорости сканирования не связан с динамикой процесса упругих деформаций и должен быть обусловлен более медленными, неупругими процессами. Т.е. разрушение поверхности имеет место в случае, когда характерное время взаимодействия зонда и локального участка образца становиться сравнимым со временами неупругих процессов, связанных с локальным воздействием зонда (при достаточной величине силы воздействия).
Модель цилиндрического образца может найти применение при анализе деформаций микрочастиц цилиндрической формы (вирусных частиц, линейных макромолекул и пр.). Однако в этом случае (контакт сферического зонда с боковой поверхностью цилиндра) соотношения Герца (2.3) не упрощаются. Ниже изложена реализация общего численного решения системы (2.3), а также проанализированы два частных случая (близких и различающихся главных значений суммарного тензора кривизны контактирующих поверхностей), для которых можно получить аналитические формулы, выражающие искомые величины явно через известные параметры задачи. Из формул (2.3) следует:
| (2.9) |
Соотношения (2.3) являются системой нелинейных уравнений, где искомые величины выражены неявно через эллиптические интегралы. Это обстоятельство осложняет реализацию численного решения, основанного непосредственно на (2.3). Поэтому проведем дополнительные аналитические преобразования, основываясь на соотношениях (2.9).
Воспользовавшись формулой [67]:
|
| (2.10) |
Соотношение (2.10) является нелинейным уравнением с одним неизвестным - отношением a/b, что позволяет реализовать его решение. Далее, из (2.9) получаем зависимость:
| (2.11) |
Таким образом мы свели систему нелинейных уравнений к независимым уравнениям: из соотношения (2.10) можно численно определить a/b, затем, из соотношения (2.11), значение b. Воспользовавшись первым из уравнений (2.9) определяем h. Так получаем все искомые параметры: h, a и b, т.е. задача решена численно в общем виде без использования упрощающих предположений. После этого давление (среднее) в области контакта также можно определить: P=F/pab. Текст программы (написанной на языке C++), реализующей данное численное решение, приведен в Приложении A.1.
Реализованное решение справедливо для произвольной геометрии зонда и образца. Необходимые для решения значения A и B определяются именно этой геометрией: для случая контакта сферы радиуса R и боковой поверхности цилиндра радиуса Rў они выражаются из соотношений (2.2) следующим образом:
| (2.12) |
Рассмотрим теперь дополнительные упрощающие предположения частных случаев.
Случай различающихся главных значений суммарного тензора кривизны контактирующих поверхностей.
В случае контакта зонда и боковой поверхности цилиндра, при условии, что радиус цилиндра меньше радиуса зонда, из формулы (2.12) следует, что главные значения суммарного тензора кривизны поверхностей различаются: A > B. Из формулы (2.10) (или непосредственно из формул (2.3)) следует, что в этом случае a < b. Если различие составляет достаточную величину, то можно воспользоваться асимптотикой полного эллиптического интеграла, справедливой при условии a2 << b2, что, очевидно, не является жестким условием [67]:
| (2.13) |
| (2.14) |
|
| (2.15) |
Формулы для параметров эллиптической области контакта a и b несколько громоздки:
| (2.16) |
| (2.17) |
| (2.18) |
Таким образом, полученные нами формулы (2.17), (2.16) и (2.18) решают поставленную задачу. Действительно, искомые параметры выражены явно через известные величины (F, D, A, B) и параметр C, который можно определить из соотношения (2.15) или воспользоваться оценкой.
Для удобства сравнения результатов точного решения (по соотношениям (2.10), (2.11) и первым из (2.9)) и приближенного (по (2.17), (2.16) и (2.18)), последнее также включено в программу решения задачи о контактных деформациях, см. A.1.
Случай близких главных значений суммарного тензора кривизны контактирующих поверхностей.
Случай близких значений величин A и B реализуется для задачи контакта сферического зонда и боковой поверхности цилиндра при условии, что радиус цилиндра много больше радиуса зонда. Тогда, в силу соотношений (2.12), действительно A ~ B и, в силу (2.10) (или непосредственно (2.3)), имеем: a ~ b. В этом случае асимптотика (2.13) теряет применимость.
При выполнении условия a ~ b следует воспользоваться другой асимптотикой полного эллиптического интеграла [67]:
|
В этом случае для параметров области контакта a и b с привлечением соотношений (2.12) и (2.9) можно вывести следующие зависимости:
| (2.19) |
| (2.20) |
Для сближения зонда и образца за счет деформации получим:
| (2.21) |
И для величины давления на образец имеем выражение:
| (2.22) |
Итак, и для рассматриваемого случая получены аналитические формулы (2.21, 2.20), решающие задачу. Для удобства сравнительного анализа с точным численным решением они также добавлены в текст программы A.1.
Выше мы рассмотрели контактные деформации в области соприкосновения зонда и образца. Однако общая деформация, определяющая занижение высоты АСМ-профиля, включает еще и вклад деформаций в области контакта образца и подложки (имеется в виду случай, когда сверху на образец давит зонд). Для этого случая остаются справедливыми приведенные выше решения, нужно лишь соответствующим образом переопределить параметры A и B.
Геометрия контакта образца радиуса Rў (на который сверху давит зонд радиуса R - образец следует рассматривать в этом случае как изогнутый цилиндр с радиусом изгиба поверхности, контактирующей с подложкой: R + 2Rў) и плоской подложки позволяет определить значения A и B по формулам (2.2):
| (2.23) |
Алгоритм анализа этого случая идейно не отличается от проведенного выше для значений A и B, определяемых формулой (2.12), его результаты вошли в таблицу 2.2.
Вообще, вклад в формулу для занижения высот АСМ-профиля объектов исследования за счет контактных деформаций при сканировании дают три слагаемых:
| (2.24) |
Однако при анализе ситуации с ``мягкой'' подложкой третье слагаемое формулы (2.24) следует учитывать. Так, из таблиц 2.1 и 2.2 следует, что занижение высот цилиндрического образца (радиус 1 нм, модуль Юнга 1010 Па, сила воздействия зонда 5 нН) в случае слюдяной подложки (E ~ 1011 Па) будет в ~ 1,4 раза больше, чем в случае пирографитовой (E ~ 1010 Па) подложки. Этот вывод согласуется с экспериментом: действительно, при использовании в качестве подложки слюды занижение высот АСМ-изображений молекул нуклеиновых кислот больше, чем при использовании пирографитовой подложки.
В качестве иллюстрации применим рассмотренный алгоритм для вычисления контактных деформаций в модельных случаях цилиндрического образца с радиусом 1 нм и радиусом 10 нм и сферического образца с теми же значениями радиуса. Результаты приведены в таблице 2.2, где для удобства сравнительного анализа для всех случаев используются одинаковые параметры: значения модулей Юнга образцов (Eў1 = Eў2 = 1010 Па), зонда (E = 1011 Па), силы воздействия зонда F = 5×10-9 Н и радиуса кривизны кончика зонда R = 10 нм. Выбор значений модулей Юнга образов (10 ГПа) обусловлен тем, что это значение является верхней границей модулей упругости биополимеров. Поэтому значения, приведенные в таблице, должны рассматриваться в смысле ``как минимум'' при АСМ-исследованиях биообъектов; в таком случае из анализа таблицы следует вывод, что роль контактных деформаций весьма значительна и должна учитываться при интерпретации результатов подобных экспериментов.
Rў, нм | Область контакта | a и b, нм | P, Па | h, нм | e, % |
сферический образец | |||||
1 нм | зонд/образец | 0,7 нм | 3×109 Па | 0,57 нм | 29% |
образец/подложка | 0,7 нм | 3×109 Па | 0,55 нм | 28% | |
суммарная деформация: | 1,12 нм | 56% | |||
10 нм | зонд/образец | 1,3 нм | 1×109 Па | 0,32 нм | 1,6% |
образец/подложка | 1,6 нм | 0,6×109 Па | 0,26 нм | 1,3% | |
суммарная деформация: | 0,58 нм | 2,9% | |||
цилиндрический образец | |||||
1 нм | зонд/образец | 0,46 и 2,2 нм | 1,6×109 Па | 0,36 нм | 18% |
образец/подложка | 0,47 и 2,4 нм | 1,5×109 Па | 0,34 нм | 17% | |
суммарная деформация: | 0,7 нм | 35% | |||
10 нм | зонд/образец | 1,1 и 1,8 нм | 0,8×109 Па | 0,29 нм | 1,4% |
образец/подложка | 1,3 и 2,7 нм | 0,4×109 Па | 0,21 нм | 1% | |
суммарная деформация: | 0,5 нм | 2,5% |
В таблице для цилиндрического образца приведены результаты численных расчетов точной методики, реализованной по соотношениям (2.9), (2.10) и (2.11). Расчеты по приближенным методикам дают незначительные различия с точным решением: для случая Rў = 1 нм решение по формулам (2.17), (2.16) дает отличие в значениях a и b около 2%, и в значении h - 0,5%; для случая Rў = 10 нм решение по формулам (2.21), (2.20) дает отличие от точного решения для a и b около 10%, для h около1%.
Из таблицы следует, что при прочих равных условиях относительные деформации объектов с меньшим радиусом кривизны существенно выше. Это подтверждается экспериментом: так относительные деформации молекул нуклеиновых кислот (радиус кривизны около 1 нм) превышают относительные деформации вирусных частиц табачной мозаики, даже в том случае, когда объекты визуализованы на одном кадре (что определяет абсолютную тождественность условий эксперимента).
Итак, применив теорию контактных деформаций, можно объяснить основные закономерности эффекта занижения высот АСМ-изображений при проведении АСМ-исследований.
Можно также показать, что зависимость величины деформации от параметров задачи имеет в общем случае (см. формулы (2.7), (2.17) и (2.21)) вид:
| (2.25) |
С целью экспериментальной иллюстрации закона h ~ F2/3 ниже приводятся результаты исследования тест-объектов - вирусных частиц табачной мозаики и молекул ДНК (имеющих цилиндрическую форму) - при различных значениях нагружающей силы при сканировании.
Для вирусных частиц табачной мозаики наблюдается хорошее совпадение эксперимента с теорией (т.е. с законом ``две третьих'' (2.25)), см. рис. 2.2
|
Экспериментальную зависимость измеряли путем анализа высот АСМ-изображений вирусных частиц, полученных при различной величине силы воздействия зонда при сканировании. Теоретическая зависимость рассчитана по изложенной выше методике (точное решение), где использовались значения радиуса зонда R = 25 нм (значение определено из эксперимента для конкретного зонда) и модуля упругости объекта Eў = 4×109 Па (значение определено аппроксимацией экспериментальных точек).
Можно предложить методику измерения модуля упругости отдельного микрообъекта, адсорбированного на поверхность твердой подложки. Для этого необходимо измерить высоту АСМ-профиля объекта при различных значениях силы воздействия зонда при сканировании, а затем, из аппроксимирующей зависимости (закона ``две третьих''), определить искомое значение D (связанное с модулем упругости объекта по (2.4)). Использование данной процедуры для определения модуля Юнга отдельных частиц ВТМ показало, что определяемые значения лежат в диапазоне 3ё4 ×109 Па. Эта методика определения упругих параметров адсорбированных объектов может успешно применяться в случае, если относительные деформации объектов невелики (иначе теряет применимость теория контактных деформаций).
Погрешности методики связаны с погрешностью калибровки сканера, погрешностью определения параметра геометрии иглы (R), определения уровня ``нулевого'' изгиба кантилевера (минимального воздействия зонда) и определения жесткости кантилевера фирмой-производителем (что приводит к погрешностям в определении значений силы воздействия зонда). Однако потенциально жесткость кантилевера можно определять с большей точностью, используя калибровку по резонансной частоте колебаний. Это может повысить точность определения упругих параметров отдельных микрообъектов до уровня, определяемого величиной набранной статистики высот объектов и дисперсией этих значений.
Из рис. 2.2 следует, что закон ``две третьих'' (2.25) справедлив для исследуемого случая в широком диапазоне сил, за исключением области минимальных воздействий. Это связано со следующим: за нулевое значение в эксперименте выбирали минимально возможное значение нагружающей силы, при котором удавалось осуществлять сканирование. Однако в силу того, что данный эксперимент проводили в естественной атмосфере, наличие капиллярных сил (капиллярного мостика), очевидно, не позволяло минимизировать силу воздействия до величины, меньшей, чем несколько наноньютонов.
Для анализа контактных деформаций молекул нуклеиновых кислот проводили сходный эксперимент, соответствующая экспериментальная зависимость изображена на рис. 2.3.
|
Однако в рассматриваемом случае деформация молекул ДНК уже не описывается законом ``две третьих'' (2.25). Вместо этого наблюдаем обычную линейную зависимость, т.е. закон Гука. Но это и естественно, поскольку теория контактных деформаций справедлива при условии малости деформаций объекта в сравнении с его радиусом кривизны. Последнее не выполняется для случая АСМ-исследований молекул ДНК, поскольку их относительные деформации велики в стандартных экспериментальных условиях, см. таблицу 2.2. Кроме того, вообще применимость континуальной теории и использование макроскопических параметров (типа модуля Юнга) для описания деформации отдельной молекулы ДНК, вызывает сомнения. Тем не менее, смысл проведенного анализа в том, что он позволяет качественно объяснить, почему относительные деформации под воздействием зонда тем выше, чем меньше радиус кривизны объектов, см. таблицу 2.2.
То, что для случая молекул ДНК экспериментально измеренные относительные деформации велики даже при малых силах воздействия зонда, объясняется, опять же, присутствием капиллярных сил (капиллярного мостика), не позволяющих минимизировать силу воздействия зонда на образец при сканировании на воздухе до значений меньших, чем несколько наноньютонов; из-за малости радиуса кривизны молекулы ДНК этого значения силы уже достаточно для того, чтобы значения относительной деформации были велики, см. таблицу 2.2.
Покажем, что характерное время установления контактных деформаций меньше типичных значений времени взаимодействия зонда и локального участка образца при сканировании. Этим обоснуем применимость статических формул для контактных деформаций в задачах, связанных с АСМ-исследованиями.
Рассмотрим контакт зонда и плоского образца. При приложении силы F в статическом случае сближение за счет деформации будет определяться формулой (2.7). Однако при мгновенном приложении силы к касающимся телам это значение h будет достигнуто не сразу, а в результате переходного процесса динамического изменения h во времени. Уравнение, описывающее рассматриваемый переходный процесс для h(t), может быть получено из элементарных рассмотрений баланса сил взаимодействия (прилагаемой извне F и обусловленной упругим откликом ЦRh3/2/D, см. (2.7)):
| (2.26) |
|
График численного решения уравнения (2.26) приведен на рисунке 2.4. При решении использовали значения таблицы A.1 для оценки масс зондов и кантилеверов (лежат в диапазоне 10-10ё10-11 кг). Для других параметров задачи выбраны следующие типичные значения: F = 5×10-9 Н, D = 10-10 Па-1, R = 10×10-9 нм. Из решения следует, что характерные времена деформации имеют порядок t ~ 10-6 сек (та же оценка может быть получена из соотношения t ~ Ц{mh/F}, которое следует из (2.26), и где h оценивается по (2.7)).
Итак, характерные времена установления контактных деформаций (т.е. характерные времена упругого отклика системы) для случая АСМ меньше или много меньше характерных времен взаимодействия зонда и локального участка образца (10-3ё10-5 сек, см. 2.1.2). За эти времена взаимодействия переходный процесс установления контактных деформаций заведомо закончится, и поэтому использование при анализе задач АСМ формул статической теории контактных деформаций вполне оправдано.
Соотношение, следующее из формулы (2.6):
| (2.27) |
Рассчитанная по формуле (2.27) оценка, что при типичных условиях АСМ-эксперимента a ~ 1 нм (см. также таблицу 2.1), означает, что площадь области контакта pa2 имеет значение S ~ 3 нм2. Если, например, сравнить с типичной величиной площади, приходящейся на одну молекулу в плотноупакованном ЛБ слое (0,2 нм2, см. раздел 5.3), то очевидно, что при сканировании в каждый момент времени имеет место контакт зонда не с одной, а с десятком и более молекул. Почему же, в таком случае, АСМ позволяет получать молекулярное (или атомное) разрешение при исследовании широкого спектра поверхностей кристаллических материалов и тонких пленок?
Авторы работы [63] для того, чтобы совместить результаты теории контактных деформаций (значительную величину области контакта) с возможностью наблюдения атомной и молекулярной структуры поверхности методом АСМ прибегают к предположению об игле, имеющей некоторую ``особенность'' ангстремных размеров. Подобное представление, что при визуализации атомной структуры с поверхностью контактирует лишь некоторая ``особенность'' (один крайний атом) зонда, широко распространено в силу своей наглядности (в этом случае говорят об ``истинном'' атомном разрешении АСМ [68]). Однако имеется также представление, что визуализация двумерной периодической структуры возможна и в случае, когда игла контактирует с исследуемой поверхностью несколькими атомами своего кончика [69] (в этом случае говорят о ``ложном'' атомном разрешении).
Анализ влияния геометрии кончика иглы (один или несколько атомов, контактирующих с поверхностью) на формирование АСМ-изображений атомной упаковки исследуемых поверхностей рассматривался в литературе. Например, авторы работ [70,71] методом компьютерного моделирования исследовали ряд моделей: взаимодействие зонда с двумерной периодической структурой (атомной структурой поверхности) в случае, если зонд контактирует с поверхностью одним, тремя, четырьмя и девятью крайними атомами (при различных расстояниях между атомами). Был сделан вывод, что ``истинное'' атомное разрешение возможно лишь при наличии единственного контактирующего с поверхностью атома иглы. Если контактирующих атомов несколько, то было показано, что и в этом случае возможна визуализация двумерной периодической структуры, характеризующейся теми же параметрами элементарной ячейки, что и реальная поверхностная решетка. Однако структура самой ячейки отображается неадекватно. При определенных условиях возможна инверсия контраста, т.е. наблюдение минимумов АСМ-изображения над атомами исследуемой поверхности и максимумов между ними. Точечный дефект (пропуск одного атома решетки) также неадекватно отображается на АСМ-изображении: имеет место перераспределение его вклада по некоторой области, и, при определенных условиях, возможна визуализация ``ложного'' атома на месте дефекта. Т.о. речь идет о достижении ``ложного'' атомного разрешения.
Авторы отмечают корреляцию полученных ими результатов с экспериментальными наблюдениями и делают вывод, что для достижения ``истинного'' атомного разрешения необходимо использовать иглу с единственным атомом на кончике (предлагая методику тестирования геометрии кончика иглы путем исследования точечных дефектов в атомной структуре поверхности тест-объекта).
В то же время, авторы рассматриваемых работ интерпретируют этот параметр - ``количество атомов на кончике иглы'' - как случайный, не допускающий контролируемого изменения со стороны экспериментатора. Напротив, ниже предлагается подход, основанный на результатах теории контактных деформаций, который связывает наблюдаемые экспериментальные особенности АСМ-визуализации атомной (молекулярной) упаковки поверхности с контролируемыми параметрами эксперимента.
Рассмотрим при каких условиях возможно достижение ``истинного'' (используя терминологию авторов [68,69]) атомного (молекулярного) разрешения. Согласно формуле (2.27) размер контактной площадки pa2 будет равен 0,2 нм2 (контакт с единственным атомом при расстоянии между атомами ~ 0,5 нм), при радиусе кривизны зонда 10 нм и достаточно жестком образце (модуль Юнга 1011 Па), лишь при минимизации величины силы воздействия зонда при сканировании до значения 0,1 нН, что требует специальных экспериментальных условий (исследования в жидких средах или в вакууме).
В случае проведения АСМ-исследований на воздухе минимизация силы воздействия на образец возможна, по-видимому, лишь до единиц наноньютонов. Поэтому при исследованиях на воздухе (особенно ``мягких'' образцов, например, тонких пленок ЛБ, см. раздел 5.3) получаемая в результате АСМ-исследования двумерная картина является ``ложным'' молекулярным разрешением.
Ниже покажем, что теория контактных деформаций допускает возможность достижения ``ложного'' атомного (молекулярного) разрешения, т.е. визуализации двумерной периодической структуры, характеризующейся теми же параметрами элементарной ячейки, что и решетка исследуемой поверхности.
В силу неоднородного распределения давления в области контакта (см. формулу (2.5)) вклад каждого атома поверхности в силовое взаимодействие зонда и образца будет определяться его положением относительно центра области контакта; количество атомов, дающих вклад во взаимодействие, определяется радиусом области контакта (2.27).
Рассмотрим следующую модель: зонд сканирует поверхность плоского образца (радиус области контакта определяется формулой (2.27)), тогда, согласно формуле (2.5), распределение давления в области контакта неоднородно, что позволяет ввести аппаратную функцию вида:
| (2.28) |
|
Введенная аппаратная функция позволяет связать получаемое АСМ-изображение (учитывая вклад каждого атома в суммарное силовое взаимодействия зонда и образца) f(x,y) с реальной геометрией атомной структуры объекта исследования j(x,y):
| (2.29) |
Физический смысл уравнения (2.29) таков: степень вклада каждого атома в силовое взаимодействие с иглой определяется тем, насколько он далеко находится от центра области контакта (это и описывается аппаратной функцией A(x-xў,y-yў)).
Применим разработанный подход для простейшей модели: вычисления АСМ-изображения в случае, когда атомная геометрия объекта описывается гексагональной решеткой с параметром ячейки d, а в качестве атомных функций взяты d-функции, см. рис. 2.6.
|
Вычислим функцию f(x,y), определяемую (2.29), для двух случаев, изображенных на рис. 2.6:
|
|
Проводя вычисления с использованием (2.28) и (2.29), обозначая для двух случаев вычисляемые функции как f1 и f2, воспользовавшись свойством d-функций, получим (пренебрегая числовыми множителями):
| (2.30) |
| (2.31) |
| (2.32) |
В силу изложенного в разделе 1.3.2 (для разрешающей способности АСМ) отличие АСМ от оптических приборов заключается в том, что для характеристики разрешающей способности и контраста изображений следует использовать абсолютную (а не относительную) разницу между максимумами и минимумами АСМ-изображения. Поэтому мы будем использовать именно разностную функцию D АСМ-изображения, а не функцию видности S (которую традиционно используют для характеристики оптических приборов).
Проведем расчет разностной функции (2.32) для случая d = 0,52 нм (параметр ячейки слюды). Результаты представлены в виде графика на рис. 2.7, где приведены зависимости разностной функции и ее модуля от радиуса области контакта a.
|
Как следует из рисунка, значение разностной функции существенно варьирует в зависимости от параметра a, определяемого соотношением (2.27). Более того, при некоторых значениях a может изменяться знак разностной функции, что должно приводить к инвертированному АСМ-изображению (инверсия контраста). Но особенно следует подчеркнуть, что по мере увеличения a (и вовлечения в контакт все новых атомов) не наблюдается тенденция к снижению максимумов разностной функции, т.е. сохраняется возможность достижения ``атомного'' разрешения. Т.о. ``ложное'' атомное разрешение может быть достигнуто и при значительной величине области контакта a (т.е., например, при ``тупой'' игле, ``мягком'' образце или значительной силе воздействия зонда F).
Эти выводы, несмотря на использование упрощенной модели, качественно согласуются с наблюдаемыми экспериментальными закономерностями. Действительно, контраст и качество АСМ-изображений молекулярной (или атомной) упаковки исследуемой поверхности существенно зависит от величины приложенной силы воздействия зонда F, определяющей, согласно (2.27), величину области контакта a.
Стоит отметить, что, несмотря на искажения реальной картины атомной или молекулярной упаковки при АСМ-визуализации (инверсия контраста, ``потеря'' одного атома в базисе для случая пирографита, различная спектральная плотность сечений вдоль эквивалентных направлений и т.п.), метод АСМ применим для определения параметров решетки поверхности. Действительно, визуализируемая двумерная структура характеризуется теми же векторами трансляции, что и реальная исследуемая поверхность (преобразование (2.29) не может исказить значений пространственных частот двумерной периодической функции исследуемой поверхности).
Важной отличием изложенного подхода в сравнении с описанными в литературе является то, что он позволяет однозначно связать особенности (контраст, качество и пр.) получаемого АСМ-изображения атомной структуры поверхности с реальными параметрами эксперимента (силой воздействия зонда, модулями упругости зонда и образца, радиусом кривизны иглы и степенью ее асимметрии), а не с абстрактным ``количеством атомов на кончике иглы''.
Эффект уширения проявляется в том, что отдельные микрообъекты, адсорбированные на плоскую подложку и визуализированные АСМ, имеют завышенные значения ширины профиля. Этот эффект может быть полезен. Так, например, при АСМ исследованиях нуклеиновых кислот он облегчает идентификацию молекул: ``уширенные'' молекулы (ширина профиля молекулы ДНК завышается в 5-10 раз) легче обнаружить на кадре значительной площади, что облегчает набор статистики. В силу этого, при исследовании нуклеиновых кислот эффект уширения позволяет обходится без дополнительного контрастирования макромолекул (комплексообразованием с белками и пр.).
Эффект уширения связан с тем, что зондирующее острие микроскопа имеет конечный радиус кривизны кончика. Эту аппаратную погрешность трудно преодолеть, поскольку уменьшение радиуса кривизны кончика зонда (использование более острых зондов) приводит к увеличению давления в области контакта (при том же значении величины контактных сил, см. формулу (2.8)); значительная величина контактного давления может привести, в свою очередь, к частичному разрушению зонда при сканировании. Уменьшить контактное давление можно при наблюдении поверхности образцов в жидкостях, поскольку в этом случае можно поддерживать контактные силы на существенно более низком уровне [48].
Подобные эффекты могут весьма существенно (особенно при значительном радиусе кривизны кончика зондирующего острия) искажать морфологию поверхности при ее визуализации методом АСМ. Наглядное сравнение результатов АСМ и электронной микроскопии при исследовании одних и тех же поверхностей проводят авторы работы [72]. Приведенные изображения показывают, что при исследовании поверхностей, характеризующихся сложной морфологией, конечность радиуса кривизны иглы может привести к потере информации на получаемых изображениях АСМ.
При исследованиях биологических и органических объектов часто возникает задача определения влияния эффекта уширения на АСМ-профиль объекта, адсорбированного на поверхность твердой подложки. Искажающее влияние эффекта уширения и контактных деформаций приводит к тому, что для восстановления реальной геометрической формы объекта исследования по его АСМ-изображению необходим дополнительный математический анализ с использованием определенных модельных представлений о геометрии зонда (конус со сферическим кончиком, параболоид и пр.) и априорных представлений о форме объекта исследования.
В работе [73] предложена универсальная компьютерная методика деконволюции АСМ-изображений, включающая два этапа: определение геометрии используемого острия с помощью тест-объектов и свертку инвертированной геометрии острия с измеренным АСМ-профилем; эта процедура позволяет во многих случаях восстановить исходный профиль объекта с высокой точностью.
Мы проводили тестирование данной методики, решая задачу восстановления геометрии объектов, адсорбированных на поверхность плоской подложки. Анализ показал, что восстановленные по данной методике латеральные размеры объекта (ширина на полувысоте) существенно завышены при условии, что радиус кривизны объекта меньше радиуса кривизны кончика зонда (ошибка тем выше, чем больше разница радиусов). Это обстоятельство осложняет применимость рассматриваемой методики при исследованиях ряда биообъектов (отдельных макромолекул, их комплексов и пр.) в силу малости размеров последних в сравнении с радиусом кривизны кончика зонда АСМ. Так, при радиусе кривизны зонда 15 нм и радиусе объекта 2,5 нм (типичные значения для исследований биообъектов) погрешность (завышение) восстановленных значений ширины составляет для тестируемой методики более 50%.
В работе [74] предложена методика восстановления объема исследуемых частиц (по АСМ-профилю), не включающая стадию предварительного тестирования зонда: как геометрия зонда, так и геометрия объектов исследования, определяемые соответствующими радиусами кривизны, может быть восстановлена путем анализа одного и того же АСМ-изображения. Однако данная методика, включающая использование для описания геометрии поверхностных структур значение радиуса кривизны, сводит задачу к рассмотрению геометрии контакта сферического зонда и сферического образца.
Возможность получения в этом случае дополнительной информации (т.е. не только параметров объекта, но и зонда из одного АСМ-изображения) может быть проиллюстрирована следующим. Ниже приводится решение задачи определения истинной ширины (2a) эллипсоидального в сечении объекта по известным параметрам: высоте (2b), ширине АСМ-профиля на полувысоте (2d) и радиусу зонда R. При этом возможным оказывается получить дополнительные соотношения (2.38, 2.39), связывающие указанные параметры. Модель сферического образца позволяет описать его геометрию одним параметром - радиусом кривизны (без использования предположения о связи этого параметра с высотой объекта над поверхностью подложки). Подобное уменьшение числа параметров задачи позволяет извлечь, путем анализа геометрии контакта, дополнительную информацию.
Авторы работы [75], отмечая, что при достаточно ``остром'' профиле исследуемых поверхностных структур на АСМ-изображении будет отображаться лишь профиль более ``тупой'' иглы, при сравнительном анализе также оперируют параметрами соответствующих радиусов кривизны (зонда и исследуемых поверхностных структур).
Но, в силу имеющихся представлений о существенной роли контактных деформаций (раздел 2.1) в исследованиях АСМ, можно предположить, что модель сферической геометрии образца и использование для ее описания соответствующего радиуса кривизны вряд ли адекватна при решении задачи восстановления геометрии адсорбированных биообъектов, характеризующихся, как известно, невысокими значениями модуля упругости. Более общей является модель, позволяющая учесть эффект уширения при контакте иглы с деформированной (сплюснутой) частицей, имеющей эллипсоидальное сечение. Решение данной задачи позволяет определить объем частицы, выразив его через параметры эллиптического сечения: найденного значения a и связанного с высотой частицы значения b = h/2. Ниже приведен алгоритм решения данной задачи.
Применим для описания эффекта уширения геометрическую модель (рис. 2.8), учитывающую взаимодействие объекта только с кончиком зонда (предполагается, что нет контакта со стенками пирамиды). Это оправдано в том случае, когда высота исследуемых структур над подложкой не превышает радиуса кривизны кончика иглы.
|
Кончик зонда будем аппроксимировать либо полусферой радиуса R, либо параболоидом вращения (z = kr2, где k - коэффициент параболы). При этом результаты применения двух методик фактически тождественны: количественное расхождение в ответах составляет менее 10%. Стоит отметить, что аппроксимация иглы с помощью полусферы более наглядна и более широко используется в литературе.
Исследуемую частицу опишем моделью эллипсоида. Т.е. мы исходим из априорных представлений о контактной деформации образца под действием зонда, см. раздел 2.1. В сечении объекта исследования - эллипс с полуосями a и b. Ставится задача: по заданным значениям высоты объекта (h = 2b), параметра геометрии зонда R (или k для параболической аппроксимации) и ширины профиля АСМ-изображения частицы на ее полувысоте (2d) найти значение истинной ширины объекта 2a.
Запишем систему уравнений для эллипса (сечение профиля исследуемой частицы) и окружности (сечение профиля иглы):
|
Для решения задачи запишем условия для точки контакта с координатами (x0, y0): равенство касательных к эллипсу и к окружности (параболы) и удовлетворение координат точки контакта уравнениям эллипса и окружности (параболы):
| (2.33) |
Данная система аналитически сводится к одному уравнению с одним неизвестным (x0), которое можно решать численно. Для сферической иглы эти уравнения выглядят следующим образом:
| (2.34) |
| (2.35) |
Аналогично, для случая параболической иглы имеем уравнение относительно x0:
| (2.36) |
| (2.37) |
Аналитически можно показать, что полученное уравнение (2.34) для сферической иглы не имеет решения (на интересующем нас интервале) только в том случае, когда выполняется система неравенств:
| (2.38) |
Аналогично, уравнение для параболической иглы (2.36) не имеет решения на интересующем нас интервале в случае, когда выполняется условие:
| (2.39) |
Смысл данных ограничений очевиден: если измеренный АСМ-профиль объекта достаточно ``острый'', то игла, которой он был получен, не могла быть ``тупой''. Полученные соотношения (2.38, 2.39) являются, по-видимому, частным случаем общего подхода к оценке геометрии иглы, изложенного в работе [76]. Автор данной работы теоретически обосновывает возможность оценки ``остроты'' иглы (определение верхнего предела) путем анализа морфологии особенностей получаемого АСМ-изображения при произвольной, но достаточно неоднородной структуре исследуемой поверхности. Однако полученные соотношения (2.38, 2.39) более просты и наглядны, чем общий подход, изложенный в работе [76], и могут быть весьма полезны при интерпретации экспериментальных результатов.
Уравнения (2.34-2.37) использованы для построения численного решения (см. текст программы A.2): уравнения (2.34) и (2.36) решаются методом вилки [77] относительно x0, далее по уравнениям (2.35) и (2.37) определяется искомое значение a.
Изложенный метод является достаточно простым, кроме того, он имеет важную особенность: проверка выполнения условий (2.38) и (2.39) позволяет заранее определить случаи отсутствия решений и получить информацию о допустимом значении параметров геометрии зонда R или k.
Определение точного значения радиуса R (или k) для конкретного зонда требует его тестирования непосредственно перед использованием (по результатам исследования тест-объектов, например, вируса табачной мозаики - удобного природного ``наноразмерного'' объекта с фиксированными геометрическими параметрами, см. раздел 3.1). Однако и в этом случае существует вероятность того, что в процессе сканирования форма иглы претерпит изменения в результате взаимодействия с объектом. В этой связи чрезвычайно полезным является получение информации о форме зонда непосредственно из АСМ-изображений объекта исследования.
При анализе АСМ-изображений объектов, для которых необходимо восстановить истинную ширину, измеряется пара значений: высота и ширина на полувысоте; результатом анализа является набранная статистика пар этих значений. Анализируя выполнимость соотношений (2.38, 2.39) при различных параметрах геометрии зонда, можно определить предельное значение R (или k) выше (для k - ниже) которого наблюдается рост числа случаев отсутствия решения (для набранной статистики пар значений параметров объекта b и d). Это и будет верхнее граничное значение для оцениваемого радиуса кривизны иглы. Типичный пример зависимости числа случаев отсутствия решения для набранной статистики пар b и h от радиуса кривизны иглы представлен на рис. 4.5.
Следует подчеркнуть, что предложенная методика не учитывает возможного дополнительного вклада в уширение, связанного с частичным увлечением образца зондом при сканировании. Эффект увлечения обусловлен латеральными силами взаимодействия зонда и образца, которые характеризуются достаточно большой интенсивностью при проведении исследований в контактном режиме на воздухе. По нашим оценкам данный эффект, в основном, может проявляться при исследовании объектов, имеющих небольшую величину площади сечения (например, одиночных молекул ДНК), и приводить к 1,5-2,5-кратному завышению значения ширины объекта 2d и, как следствие, восстанавливаемого значения a.
Ослабить влияние эффекта увлечения образца зондом можно путем уменьшения интенсивности силового латерального взаимодействия зонда и образца. Это достигается при применении режима прерывистого контакта, особенно в случае проведения исследований в жидких средах.
Основная погрешность восстановленных значений a вносится погрешностью измерения полуширины АСМ-профиля объекта d. В меньшей степени погрешность решения (Da) зависит от погрешностей параметра геометрии зонда (DR или Dk) и погрешности измеренной высоты объекта (Db).
Погрешность измерения d может быть обусловлена описанным выше эффектом частичного увлечения объектов зондом при сканировании. Значительная асимметрия иглы также увеличивает эту погрешность, поскольку значения ширины в этом случае зависят от ориентации объектов. В силу этого стандартное отклонение восстановленных значений истинной ширины объектов исследования по результатам анализа их АСМ-изображений может превышать 20%. Увеличить точность результата исследования можно путем уменьшения эффекта увлечения (исследования в режиме прерывистого контакта в жидкости), а также (потенциально) путем определения и учета степени асимметрии иглы. Кроме того, величина ошибки может быть уменьшена при наборе достаточно большой статистики анализируемых параметров АСМ-изображений объектов.
Особенностями изложенного алгоритма восстановления реальной ширины объектов, адсорбированных на поверхность подложки, являются:
Этот алгоритм может служить прототипом для общего подхода восстановления геометрии объектов с учетом трехмерной геометрии иглы. Действительно, возможно обобщение задачи на трехмерный случай, что позволит использовать модель асимметричной иглы.
Вирусные частицы являются простейшими природными объектами, способными к репликации. Они состоят из носителя генетической информации (молекулы нуклеиновой кислоты) и оболочки, образованной белковыми молекулами. Репликация вирусных частиц осуществляется только внутри клеток живых организмов, при этом ее процесс регулируется спецификой взаимодействия молекулы нуклеиновой кислоты с белками.
Вирус табачной мозаики представляет собой полый цилиндр, имеющий длину около 300 нм, внутренний диаметр 4 нм и внешний диаметр 18 нм (согласно данным электронной микроскопии [27]). Каждая вирусная частица состоит из 2200 белковых субъединиц, организованных в правовинтовую спираль с шагом 16 1/3 субъединиц на виток. Вирусная РНК, состоящая из 6400ё6600 нуклеотидов, также образует спираль, каждый виток которой располагается между двумя витками спирали белковых субъединиц, см. рис. 3.1.
|
Вирус табачной мозаики является удобной модельной системой для апробирования методик приготовления образцов, позволяющих проводить АСМ-исследования взаимодействия вирусной РНК с белками.
В настоящее время разработаны различные методы фиксации макромолекул ДНК и РНК на поверхностях различных подложек для исследований СЗМ, см. раздел 1.5.1. Так, авторы работы [78] использовали традиционную методику иммобилизации нуклеиновых кислот на поверхности слюды катионами Mg2+, добавляемыми в рабочий раствор. В качестве объектов исследования они выбирали, в частности, короткие (25 и 50 оснований) молекулы однонитевой ДНК, а также РНК-гомополимера поли-A. Авторы обнаружили, что по используемой методике препарирования образцов только молекулы поли-А возможно иммобилизовать на подложке в расправленном состоянии - в отличие от адсорбированных молекул однонитевой ДНК, которые характеризовались компактной глобулярной структурой. Это наблюдение показывает, что основным препятствием при расправлении на подложке однонитевых молекул нуклеиновых кислот является образование внутримолекулярных пар оснований (по механизму Уотсона-Крика [25]). Действительно, молекулы ДНК, используемые авторами работы [78], могли формировать подобные внутримолекулярные пары (судя по приводимой авторами последовательности нуклеотидов этих молекул), в отличие от молекул поли-A, являющегося гомополимером, образованным лишь одним типом нуклеотидов.
Традиционная методика приготовления РНК в расправленном состоянии для электронной микроскопии включает образование комплекса РНК с монослойной пленкой белков (цитохром С или диизопропилфосфаттрипсин) на поверхности раздела фаз воздух/вода и последующий перенос пленки с молекулами РНК на поверхность пленки-подложки [30]. Данная методика позволяет получать ЭМ-изображения молекул РНК, находящихся в двумерном расправленном состоянии [79], однако в этом случае тот факт, что РНК находится в комплексе с белками (в белковой ``шубе''), может ограничить решение задач, связанных с тонкой структурой РНК и с исследованиями РНК-белкового взаимодействия.
Наблюдение свободной РНК позволяет проводить как физическое картирование элементов вторичной структуры РНК, так и исследовать специфическое взаимодействие с ней белков, узнающих определенные (например, регуляторные) последовательности нуклеотидов. В свою очередь, исследование взаимодействия вирусной РНК и белков может дать вклад в понимание механизмов процесса сборки и размножения вирусов, структурной организации и функционирования рибонуклеопротеидов.
Ниже изложена методика приготовления образцов для АСМ-исследования вирусной РНК, апробированная совместно с д. х. н. Ю. Ф. Дрыгиным (НИИ ФХБ им. А. Н. Белозерского МГУ).
Для приготовления всех растворов использовали деионизованную бидистиллированную воду. Препарат вируса (ВТМ, штамм U1), полученный методом дифференциального центрифугирования и центрифугирования в градиенте плотности раствора сахарозы, был любезно предоставлен Е. Н. Добровым (НИИ ФХБ МГУ).
В экспериментах по депротеинизации использовали суспензию ВТМ в воде при концентрации 2,6 мг/мл. Депротеинизацию вирусных частиц производили:
Для разделения отдельных молекул белковой оболочки ВТМ, свободной РНК и частично разрушенных вирусных частиц использовали метод высокоэффективной жидкостной хроматографии с гелем Toyo-Pearl HW-55. В трех используемых методиках колонки (150 мкл) предварительно уравновешивали водными растворами 36% ДМСО, 1,6 М мочевины или глицина (концентрация 0,1 М при pH 10,3), соответственно.
В качестве подложек использовали как чистые (свежий скол), так и модифицированные поверхности высокоориентированного пиролитического графита (пирографита) и слюды. Поверхность подложки модифицировали или катионным ПАВ - бензилдиметилалкиламмоний хлоридом (BAC) (1% раствор в формамиде, экспозиция 45 мин) [83] или 100 мМ MgCl2 (см. раздел 1.5.1), экспозиция 75 мин. Процесс модификации завершали промывкой водой (в 0,2 мл четыре раза по 1-2 мин), излишек воды удаляли фильтровальной бумагой и подложки высушивали в струе фильтрованного сжатого воздуха.
Молекулы РНК, высвобожденные из белковой оболочки ВТМ, или частично разрушенные вирусные частицы посредством обработки мочевиной или диметилсульфоксидом наносили на подложку из 1,6 М раствора мочевины или 36% диметилсульфоксида соответственно, разделив предварительно эти структуры на фракции по размерам методом высокоэффективной жидкостной хроматографии, см. выше.
Схема эксперимента без использования хроматографии включала следующие этапы. Вирусные частицы ВТМ (0,138 мг/мл) в 2 мМ MgCl2, 200 мМ KCl экспонировали 50% ДМСО и 3,7% формальдегидом в течении 5 мин при 45°С для денатурации РНК. Затем к смеси добавили 1% BAC и адсорбировали каплю препарата на поверхность слюды.
Молекулы РНК вируса энцефаломиокардита мышей, выделенные из вирусных частиц, иммобилизовали на подложке по следующей схеме. Макромолекулы (концентрация 60 мкг/мл) в 0,01 М HEPES буфере (pH 7,4) и 7,5% формальдегиде экспонировали в течении 15 мин при 55°C для денатурации. Затем или к препарату добавляли равный объем 0,1% раствора BAC в формамиде (до конечной концентрации BAC 0,05%) или к 3 частям препарата добавляли 27 частей того же раствора BAC (до конечной концентрации BAC 0,01%).
Адсорбцию исследуемых структур осуществляли (во всех случаях) из капли препарата (3-15 мкл), которую наносили на поверхность подложки и экспонировали во влажной среде под колпаком. Длительность процесса экспозиции варьировали от нескольких секунд до часа и более. После завершения экспозиции поверхность подложки также промывали водой и высушивали в струе фильтрованного сжатого воздуха.
АСМ-исследования проводили на воздухе в режимах постоянного или прерывистого контактов на приборе ``Nanoscope-IIIa'' (``Digital Instruments'', США). При проведении исследований в контактном режиме использовали кантилеверы жесткостью 0,06 Н/м из нитрида кремния. Для исследований в режиме прерывистого контакта применяли кремниевые кантилеверы с резонансной частотой в диапазоне 350-380 кГц, (параметры кантилеверов см. в таблице A.1).
Частицы вируса табачной мозаики были исследованы в режимах постоянного и прерывистого контактов на поверхностях пирографита и слюды. Было обнаружено, что скорость сорбирования вирусных частиц на свежий скол поверхности существенно выше при использовании пирографитовой подложки, чем при использовании подложки из слюды. В первом случае при равных условиях сорбирования (температура, объемная концентрация вирусов в препарате, время промывки и пр.) поверхностная плотность вирусных частиц существенно выше. В то же время на пирографитовой подложке при сканировании в режиме постоянного контакта вирусные частицы нестабильны и увлекаются иглой даже при минимизации воздействующей силы до величины менее 1 нН. Стабильные АСМ-изображения вирусных частиц на пирографите удается получить только в режиме прерывистого контакта. В случае же использования слюдяной подложки сила адгезии вирусных частиц достаточна для проведения исследований и в режиме постоянного контакта с величиной воздействующей силы в десятки нН.
Большая скорость сорбции вирусных частиц на поверхность пирографита в сравнении со слюдой может объясняться влиянием гидрофобного эффекта, см. раздел 1.3.1. В пользу данного утверждения говорит тот факт, что в случае гидрофобизации поверхности слюды посредством обработки BAC скорость сорбции вирусных частиц увеличивается.
Нестабильность адсорбированных на пирографит вирусных частиц в процессе сканирования в режиме постоянного контакта объясняется слабой адгезией частицы, обусловленной инертностью подложки (гидрофобный эффект не дает вклад в силу адгезии после высушивания образцов, поскольку проявляется только в водной среде).
АСМ-изображения частиц вируса табачной мозаики приведены на рис. 3.2. Согласно анализу полученных данных, для вирусных частиц среднее значение длины составляет 310 нм, при этом среднеквадратичное отклонение от среднего составляет 30%. Столь широкое распределение объясняется, во-первых, частичным разрушением вирусных структур при приготовлении препарата и, во-вторых, агрегацией вирусных частиц торец в торец (это, в частности, проявляется в наличии пиков на удвоенной, утроенной и т.д. длине вирусной частицы). С началом процесса обработки вирусных структур с целью высвобождения вирусной РНК максимум распределения вирусных частиц по длинам смещается в более короткую область.
| ,
Figure |
Анализ высот АСМ-изображений вирусных частиц над поверхностью подложки дает значение 18±2 нм (независимо от режима исследования), что согласуется с результатами электронной микроскопии [27]. Как показали дополнительные исследования, высота АСМ-профиля вирусной частицы над поверхностью подложки уменьшается при увеличении силы воздействия зонда, что связано с контактными деформациями, см. раздел 2.1; результаты экспериментального исследования эффекта уменьшения высоты АСМ-изображений частиц ВТМ при увеличении силы воздействия изображены на рис. 2.2.
Частицы ВТМ, являясь природными ``наноразмерными'' объектами с фиксированными геометрическими параметрами, могут служить удобным тест-объектом для анализа качества зондирующего острия (определения параметров его геометрии). При минимизации силы воздействия зонда величина деформации вирусных частиц незначительна (см. экспериментальную зависимость рис. 2.2 или рассчитанные значения деформации для модельных случаев, приведенные в таблице 2.2). В силу этого можно считать, что частица вируса имеет в своем сечении окружность. Оценить эффект уширения (рассмотрению которого посвящен раздел 2.2) для частицы, имеющей в сечении окружность, существенно проще, чем для частицы, имеющей эллипсоидальную форму (этот случай подробно рассматривается в разделе 2.2). Согласно элементарным геометрическим построениям для острия с радиусом кривизны кончика R значение реально измеряемой ширины на полувысоте АСМ-изображений вирусных частиц радиуса r должно составлять величину:
| (3.1) |
Экспериментально определенные значения ширины АСМ-изображений вирусных частиц для режимов прерывистого и постоянного контактов были равны 30±5 нм и 32±6 нм соответственно. Воспользовавшись формулой (3.1), можно оценить радиус кривизны для кончиков зондов используемых типов кантилеверов: R ~ 5ё15 нм.
В мягких условиях депротеинизации наблюдали рибонуклеопротеиды (РНП) с выходящими из них (с одного конца или с обоих) нитями РНК, см. рис. 3.3 и 3.4. Мы протестировали три известные методики частичной депротеинизации ВТМ и высвобождения РНК, описанные в литературе [80,81,82]. Сравнительный анализ показал, что наиболее перспективным является метод с применением ДМСО (рис. 3.4). Метод депротеинизации в растворах со щелочным значением рН характеризуется большим фоном получаемых АСМ-изображений. Кроме того, молекулы РНК в этом случае склонны к агрегации и в значительной степени сохраняют элементы вторичной структуры. Метод с использованием мочевины позволяет получать АСМ-изображения с незначительным фоном, но при этом наблюдается агрегация вирусных РНП (рис. 3.3).
| ,
Figure |
| ,
Figure |
АСМ-изображения рисунков 3.3 и 3.4 получены в топографическом режиме Z(X,Y) и в режиме измерения сил трения Fтр(X,Y), (см. раздел 1.3); последние изображения, в силу значительной разности высот частиц ВТМ и молекул РНК, характеризуются лучшим контрастом в сравнении с топографическими, позволяя различать на АСМ-изображении оба типа структур. Белковые остовы могут рассматриваться как маркеры молекул РНК, позволяя уверенно отождествлять наблюдаемые протяженные нити, выходящие из разрушенных вирусных структур, именно с молекулами вирусной РНК.
Наблюдение вирусных остовов с двумя нитями РНК, выходящими из одного конца вирусной частицы, может объясняться лишь процессом частичного их восстановления после того, как экспозиция с целью разрушения была завершена. Действительно, согласно работе [84], только в результате самосборки вирусной частицы возможно наблюдение двух нитей РНК, выходящих из вирусной частицы с одного конца, что связано с особенностью процесса самоорганизации вирусной частицы. Процесс самосборки вирусных частиц мог протекать только до разделения компонент препарата методом высокоэффективной жидкостной хроматографии, поскольку в процессе этой обработки молекулы белковой оболочки, молекулы свободной РНК и частично разрушенные вирусные структуры разделялись на отдельные фракции, что исключало возможность восстановления разрушенных вирусных частиц.
На заключительной стадии обработки вирусных частиц белковая оболочка разрушается, и макромолекула РНК полностью высвобождается. На рис. 3.5 приведены типичные результаты АСМ-исследования свободных РНК, адсорбированных на подложку из слюды.
| ,
Figure |
Молекулы, изображенные на рис. 3.5а и б, отделены от частично разрушенных вирусных частиц и молекул белковой оболочки методом высокоэффективной жидкостной хроматографии и адсорбированы на свежий скол слюды из 36% ДМСО; именно присутствие денатурирующего ДМСО способствует расправлению молекул.
Мы протестировали ряд методик модификации слюды до нанесения препарата: обработка в течении 75 мин 100 мМ MgCl2, 150 мМ KCl, 100 мМ MgCl2 + 150 мМ KCl; обработка 1% раствором BAC в формамиде или APTES (2 часа при 105° C), см. раздел 1.5.1. Было обнаружено, что методика с использованием MgCl2 приводит к компактизации молекул ДНК, адсорбируемых на подложке. Оптимальной, по-видимому, является методика с применением BAC; при этом хорошие результаты были получены и в случае, когда использовали свежий скол слюды, а BAC (конечная концентрация 0,01-1%) добавляли непосредственно к препарату, содержащему исследуемые макромолекулы. Таким образом, именно присутствие BAC способствует иммобилизации молекул на подложке.
Молекулы, изображенные на рис. 3.5в, также иммобилизованы на поверхности свежего скола слюды в присутствии 1% ВАС; методика приготовления препарата не включала в этом случае разделение его на фракции, поэтому для предотвращения самосборки частиц ВТМ молекулы РНК денатурировали обработкой 3,7% формальдегидом (см. описание экспериментальных методик 3.1).
Экспериментальный анализ ширины макромолекул РНК, адсорбированных на слюду, дает значения 10ё15 и 20ё30 нм для режимов прерывистого и постоянного контактов соответственно. Анализ значений высот макромолекул показал, что эти величины достаточно сильно варьируют как для различных молекул, так и для различных участков одной и той же молекулы, и заключены в диапазоне 0,3ё1,5 нм (для исследований как в режиме постоянного, так и прерывистого контактов). Данный разброс может объясняться сохранением макромолекулой элементов вторичной структуры (неполное разворачивание на поверхности подложки, наличие скрученных участков).
О частичном сохранении элементов вторичной структуры макромолекул говорит также наличие петель и ``шпилек''. Расправлению макромолекул может способствовать влияние потока при высыхании капли препарата (см. рис. 3.5). Это влияние обусловлено силами, действующими со стороны высыхающей капли, на молекулы, прикрепленные к подложке. Природа этих сил может быть гидродинамической или электростатической: действительно, капля препарата, находящаяся на поверхности слюды, обогащена ионами K+, перешедшими в раствор при частичной диссоциации слюды, поверхность молекулы РНК заряжена отрицательно. При высыхании капли концентрация катионов K+ увеличивается, что может вызвать электростатическое взаимодействие капли с молекулой РНК и объяснять их наблюдаемую ориентацию.
При этом на АСМ-изображении четко идентифицируются участки прикрепления макромолекулы к подложке. Видно, что молекула, будучи закрепленной в данной точке, ориентируется, расправляясь по потоку. Стабилизация макромолекул на поверхности подложки в расправленном состоянии позволяет провести анализ распределения макромолекул РНК по длинам. В то же время, если силы воздействия на молекулы значительны, то это может приводить к их разрыву, что способно затруднить интерпретацию получаемых распределений РНК по длинам, см. ниже.
На рисунке 3.6 представлена гистограмма распределения по длинам молекул РНК, высвобожденных из белковой оболочки посредством обработки ДМСО. Характер распределения одинаков как для молекул, адсорбированных из 36% ДМСО после разделения на фракции методом хроматографии, так и для молекул, адсорбированных в присутствии 1% BAC и денатурированных 3,7% формальдегидом (в обоих случаях адсорбция осуществлялась на свежий скол слюды). Из рисунка видно, что гистограмма распределения может быть описана суперпозицией трех гауссовых распределений (пунктирными линиями показаны соответствующие аппроксимирующие распределения, сплошной линией их сумма), характеризующихся значениями средней длины молекулы [`l]: 0,7 мкм, 1,4 мкм и 2,0 мкм. Значения ширины соответствующих распределений: 0,5 мкм, 0,3 мкм и 0,7 мкм соответственно.
|
Исследования зависимости длины высвобожденной РНК от длины частично разрушенного вирусного остова позволили сделать вывод, что полностью расправленным молекулам РНК соответствует пик со значением максимума 2,0 мкм. Поскольку РНК исследуемого штамма U1 включает около 6400 нуклеотидов, можно оценить величину межнуклеотидной связи в молекуле РНК, обработанной и иммобилизованной на подложку слюды по описанной методике: lн/н=0,31±0,02 нм.
С целью проверки универсальности апробированной при визуализации РНК ВТМ методики иммобилизации молекул на подложке с применением BAC, были проведены сходные исследования РНК вируса энцефаломиокардита (ВЭМК) мышей. Молекулы РНК, выделенные из препарата вируса, были денатурированы обработкой 7,5% формальдегидом в течение 15 мин при 55°C и адсорбированы на свежий скол слюды в присутствии 0,01 и 0,05% BAC, описание экспериментальных методик см. 3.1. Применение данного алгоритма позволило визуализировать молекулы исследуемой РНК в расправленном состоянии, см. рис. 3.7.
Этот результат позволяет сделать вывод о возможность иммобилизации однонитевой высокомолекулярной вирусной РНК (при условии ее денатурации) на поверхности слюды в расправленном состоянии по методике с использованием BAC; причем разработанный подход является, по-видимому, универсальным.
Иммобилизация РНК ВЭМК мышей на подложке в расправленном состоянии позволила провести анализ распределения их по длинам, см. рис. 3.8. При этом обнаружено, что характер распределения не зависит от используемой концентрации BAC.
|
Обнаружено, что распределение РНК ВЭМК мышей по длинам, изображенное на рис. 3.8 характеризуется несколькими максимумами, как и распределение по длинам молекул РНК ВТМ, приготовленных по сходной методике, см. рис. 3.6. Поэтому можно предположить, что специфический характер распределений обусловлен именно используемой методикой приготовления образцов.
Разработка методик иммобилизации вирусной РНК на поверхности подложки имела своей целью отработку рутинного подхода к исследованию специфического РНК/белкового взаимодействия; к картированию как данного взаимодействия, так и элементов вторичной структуры РНК с целью определения местоположения специфических последовательностей нуклеотидов. Тот факт, что молекулы РНК, иммобилизованные на подложке по апробированной методике, характеризуются несколькими максимумами в распределении по длинам, может быть учтен при решении задачи картирования следующим образом. Если построенная карта будет характеризоваться также несколькими максимумами, то это будет свидетельствовать о разной длине межнуклеотидной связи молекул, адсорбированных на подложке по используемой методике. В этом случае для определения адекватного местоположения исследуемой специфической последовательности нуклеотидов каждый максимум для расстояния от конца молекулы до специфического участка должен быть нормирован на соответствующий максимум для полной длины молекулы. Если же построенная карта будет характеризоваться единственным максимумом, то его следует нормировать на величину родительского максимума распределения молекул по длинам.
Молекула дезоксирибонуклеиновой кислоты является носителем генетической информации живых организмов (совокупности данных, необходимых для синтеза белков), которая ``записана'' в структуре молекулы в ``четвертичном'' (по числу нуклеотидов) коде. Согласно модели Уотсона и Крика [25] молекула ДНК представляет собой две нити, которые формируют двойную спираль, стабилизированную взаимодействием пар комплементарных нуклеодтидов (вторая нить образована нуклеотидами, комплементарными нуклеотидам первой нити).
В живой природе (внутри вирусов и клеток) молекула ДНК находится в существенно ``компактизованном'' состоянии, занимая объем на несколько порядков (до 106) меньший, чем объем молекулы ДНК в растворе (в хорошем растворителе). Таким образом, именно компактная конформация (конформация ``глобулы'') является ``естественной'' для молекулы ДНК, хотя в лабораторных исследованиях чаще имеют дело с ДНК в растворах, т.е. расправленной (конформация ``клубка'').
Перевод молекулы ДНК из развернутого в компактное (конденсированное) состояние может быть осуществлен посредством разнообразных конденсирующих агентов и пр.; структуру, образованную в результате конденсации ДНК в присутствии полимеров и солей, иногда называют psi-ДНК (psi расшифровывается как Polymer and Salt Induced) или y-ДНК [85]. Начало исследованию процессов компактизации ДНК было положено Лерманом [86], который методом седиментации обнаружил, что в разбавленных растворах ДНК, при увеличении концентрации нейтрального полимера (полиэтилен гликоля) и соли одновалентного металла (NaCl) выше некоторого порогового значения, формируются частицы высокой плотности.
В настоящее время установлено, что конденсацию молекул ДНК, находящихся в водных растворах, способны вызвать различные мультивалентные катионы: спермидин3+ [87], спермин4+ [88,89], катионные полипептиды: полилизин, гистоны [90], катионные ПАВ [91,92]. Двух- (и более) валентные металлические катионы способствуют конденсации молекул ДНК, находящихся в водно-спиртовых (этанол, метанол, изопропанол) средах [93].
Структура компактизованной ДНК исследовалась различными методами (краткий обзор основных результатов приводится в работе Блюмфельда [94]) при этом непосредственная визуализация объектов исследования осуществлялась методами электронной [90,95,96], атомно-силовой [97,98] и (с меньшим разрешением) флуоресцентной [91,92,99] микроскопии.
Результаты исследований показали [94], что двумя основными формами y-компактизованной молекулы являются тороидальная и стержневая. Для обеих форм сечение компактизованных структур, имеет, по-видимому, одно и то же значение (если структуры сформированы одинаковыми молекулами). Достаточно интригующим представляется то, что размер (радиус) тороидальной структуры не зависит от длины молекулы (если она лежит внутри диапазона 400- 40000 п.о.) и определяется, по-видимому, внутренними механизмами процесса конденсации; типичные значения внешнего радиуса тороидальной частицы Rout лежат в диапазоне 35-90 нм [94]. Длинные молекулы ДНК образуют компактные структуры из одной молекулы, при этом удавалось наблюдать тороидальную геометрию компактной частицы для молекул, длина которых не превышала 40000 п.о. Короткие ДНК формируют структуры, образованные в результате самоорганизации нескольких молекул, при этом совсем короткие молекулы (менее 150 п.о.), по-видимому, не могут быть компактизованы [100].
Метод электронной микроскопии, позволяя достигать высокого пространственного разрешения в латеральном направлении, дает значительную ошибку при определении высоты объекта, если она сравнима с размером зерна напыления; это обстоятельство не позволяет с большой точностью определять объем исследуемых частиц (необходимый для анализа их количественного молекулярного состава) по результатам исследований ЭМ. Более того, при расчете объемов частиц компактизованной ДНК часто не задаются целью определить высоту структур над подложкой, полагая, что сечение объекта представляет собой окружность [94]; нам представляется, что такое предположение не всегда обосновано (см. раздел 4.1.2).
Стоит также отметить, что специфика препарирования образцов для исследований ЭМ может искажать структуру объекта исследования.
Этого недостатка лишен метод флуоресцентной микроскопии, позволяющий проводить в реальном масштабе времени исследования компактизации молекул ДНК, находящихся в нативных условиях (в растворах с подходящим pH). Однако разрешающая способность метода ограничена пределом оптического микроскопа l/2. Поэтому в качестве объектов исследования приходится выбирать достаточно большие молекулы ДНК (например, гигантскую ДНК Т4, которая включает 166 тыс. п.о.) [91,92,99].
Применение метода АСМ для исследований компактизации молекул ДНК представляется весьма перспективным, поскольку оказывается возможным проводить исследования в нативных или близких к нативным условиях, при этом, как было показано в разделе 2.2, геометрия исследуемых объектов может быть восстановлена путем анализа измеренных параметров их АСМ-изображений (высоты над подложкой и ширины на полувысоте). В то же время, разрешающая способность АСМ существенно выше, чем у метода флуоресцентной микроскопии и близка разрешающей способности ЭМ. Однако в настоящее время в литературе имеется лишь незначительное число публикаций [97,98], посвященных применению СЗМ для исследований процессов компактизации ДНК.
Авторы работы [97] исследовали начальные стадии конденсации плазмидной поли-А ДНК в результате ее взаимодействия со спермидином. Используемые концентрации спермидина были малыми (верхняя граница 150 мкМ) поскольку, как сообщают авторы, при больших концентрациях спермидин негативно влияет на возможность АСМ-визуализации исследуемых структур.
Эксперименты проводились на воздухе и адсорбция осуществлялась из капли препарата на поверхность свежего скола подложки (слюды). Авторы обнаружили, что при концентрации спермидина в диапазоне 7,5-15 мкМ отдельные участки ДНК образуют петли, радиус которых составляет величину около 20 нм. При больших концентрациях спермидина имеет место межмолекулярная агрегация ДНК, причем, в исследуемом диапазоне, с увеличением концентрации спермидина возрастают размеры агрегатов; специфическая форма этого класса структур оправдывает название ``цветка'', которое используют авторы для описания морфологии комплекса.
Авторы полагают, что наблюдаемые образования могут являться промежуточной стадией компактизации молекул в тороидальные структуры. Однако обращает на себя внимание специфическая плоская геометрия визуализованных образований: они характеризуются микронными значениями латерального размера при высоте всего 3-7 нм. Представляется, что в исследуемой системе весьма важен вклад влияния подложки на геометрию компактизованных комплексов, поэтому обнаруженная специфическая морфология может быть связана именно с тем, что изучаемая компактизация является, по сути, компактизацией на поверхности подложки; в этом случае естественно ожидать иной геометрии компактных структур, в сравнении с компактизацией в объеме.
Наиболее широко используемым способом иммобилизации молекул ДНК на подложке в расправленном состоянии 1.5.1 является метод добавления катионов металлов в рабочий раствор, содержащий нуклеиновые кислоты, с последующей адсорбцией капли препарата на поверхность слюды (свежий скол).
Типичные концентрации компонент в препарате, наносимом на поверхность слюды, обычно выбирают в следующем диапазоне: концентрация ДНК около 0,5 мкМ×п.о., концентрация MgCl2 или NiCl2 около 1 мМ. Каплю препарата объемом 0,5-5 мкл наносят на свежий скол слюды. На рис. 4.1а приведены результаты АСМ-исследования молекул ДНК бактериофага ФX-174 (5386 п.о.), которые были адсорбированы на слюду по описанной методике в присутствии NiCl2.
Другой метод фиксации макромолекул на подложке, это добавление к рабочему раствору, содержащему ДНК, молекул катионных ПАВ в малых концентрациях (менее 0,1%), см. рис. 4.1б. Типичные концентрации ДНК в препарате выбирают в том же диапазоне, что и выше. По-видимому, в этом случае имеет место комплексообразование молекул ДНК с молекулами ПАВ, образующими тонкую пленку, покрывающую снаружи каплю препарата, наносимую на подложку. Этот экспериментальный подход эквивалентен изложенному в разделе 3.1, где для иммобилизации вирусной РНК применяли метод комплексообразования молекулы нуклеиновой кислоты с ВАС, также являющимся катионным ПАВ.
Компактизация ДНК при взаимодействии с катионными ПАВ и прохождение комплексов ДНК-ПАВ через межфазную границу вода/малополярный органический растворитель может рассматриваться как отдаленная модель транспорта генетической информации в клетку [101].
Ранее процесс компактизации ДНК при взаимодействии с катионными ПАВ исследовался, в основном, лишь методом флуоресцентной микроскопии (ФМ) [91,92,99]. Это объясняется тем, что данный метод позволяет визуализировать только те объекты, которые ``окрашены'' молекулами флуоресцентного красителя, поэтому вклад во флуоресцентное изображение дает лишь молекула ДНК (именно с ней связываются молекулы красителя), а присутствие в растворе молекул ПАВ не влияет на получаемое изображение. В этом отличие и преимущество данного метода в сравнении с электронной и атомно-силовой микроскопией, для которых вопрос об адекватном распознавании и идентификации наблюдаемых структур может являться весьма нетривиальным. Этот аспект особенно существенен при исследованиях, связанных с молекулами ПАВ, в силу их способности к самоорганизации, в результате которой возможно формирование структур с различной геометрией [102].
Согласно результатам исследований авторов работы [59] свободная молекула ДНК характеризуется невысокой проводимостью при исследовании СТМ, что выражается в визуализации молекулы по способу отрицательного контраста (т.е. при движении над макромолекулой, в силу диэлектрических свойств последней, сканирующая игла опускается ниже, чем при движении над подложкой). Несмотря на это, мы обнаружили при исследовании комплексов ДНК-ПАВ, что последние могут быть успешно визуализованы методом СТМ (по способу прямого контраста).
Ниже изложена методика приготовления образцов для СЗМ-исследования комплексов ДНК-ПАВ, апробированная совместно с к.х.н. О. А. Пышкиной и к.х.н. В. Г. Сергеевым (Химический факультет МГУ).
В эксперименте использовали ДНК из эритроцитов цыплят (``Союзхимреактив'', Россия) 1000-3000 п.о. и дистеарилдиметиламмоний хлорид (ДСДАХ) (``Tokyo Kasei Kogyo Co.'', Япония).
Эксперименты проводили на СТМ СКАН-8 (``Центр перспективных технологий'', Россия) [103], в режиме постоянного тока (рабочие параметры: Ut ~ 0,3 В, It ~ 0,5 нА). АСМ-эксперименты проводили на приборе ``Nanoscope-IIIa'' (``Digital Instruments'', США) в режиме контакта с использованием кантилеверов из нитрида кремния (см. таблицу A.1).
Раствор ДНК (10-5 М) и ПАВ (9×10-5 М) в 0,5-М ТВЕ буфере объемом 2 мкл помещали на подложку (свежий скол пирографита или слюды) и высушивали в течение 3 мин. При нанесении комплексов из хлороформа их первоначально формировали в водной фазе при взаимодействии ДНК с ПАВ, затем комплексы высушивали и растворяли в хлороформе (в отличие от свободных молекул ДНК их комплексы с ПАВ растворимы в хлороформе [104]).
Было обнаружено, что комплексы ДНК-ПАВ могут быть визуализированы методом СТМ, см. рис. 4.2 а и б, где исследуемые объекты нанесены на подложку пирографита из водной фазы (а) и хлороформа (б). Для исключения эффекта влияния подложки и проверки адекватности применения метода СТМ к исследованию диэлектрических структур проводили контрольные измерения с помощью атомно-силового микроскопа. При этом были визуализованы объекты с идентичной морфологией (рис. 4.2 в).
Было обнаружено, что частицы комплекса ДНК-ПАВ хорошо адсорбируются на поверхность гидрофобной пирографитовой подложки из водной фазы, что может объясняться гидрофобным эффектом, см. 1.3.1 (поверхность исследуемых комплексов характеризуется гидрофобными свойствами [101]). Методика нанесения частиц комплекса из хлороформа позволяет избежать агрегации исследуемых структур в процессе приготовления образцов. Это может быть связанно с тем, что в хлороформе гидрофобное взаимодействие отсутствует и агрегировавшие образования реассоциируют.
Можно предположить, что возможность визуализации таких объектов, как комплексы ДНК-ПАВ, методом СТМ связана с поверхностной проводимостью слоя молекул ПАВ, покрывающих комплекс снаружи. Это предположение согласуется с обнаруженной нами ранее (см. работу [105]) возможностью СТМ-визуализации бислойных липидных образований - липосом. Действительно, исследованные липосомы были образованы молекулами липидов, свойства и структура которых близки молекулам используемого ПАВ (ДСДАХ).
В силу того, что комплексы ДНК-ПАВ нерастворимы в водной фазе, но растворимы в органических растворителях (хлороформе, гептане, циклогексане) [104], можно провести исследование процесса их переноса через границу раздела фаз в органический растворитель и исследовать морфологию перешедших комплексов. Подобная система может служить модельной при изучении механизма трансмембранного переноса ДНК [101].
Знание геометрии компактизованных комплексов является весьма важным, например, для определения состава комплекса (количества образующих его молекул). Точность определения геометрии органического объекта по результатам исследований электронной микроскопии ограничена (трудно получить точную информацию о высоте объекта: методика запыления под углом дает значительную ошибку при малых высотах объекта).
Определение геометрии комплексов ДНК-ПАВ по результатам относительных методов исследования (см. раздел 1.4) осложнено в силу вариации комплексов по размерам.
Метод АСМ позволяет непосредственно определять форму и размеры отдельных объектов, адсорбированных на поверхность подложки, но для восстановления их реальной геометрии необходим дополнительный математический анализ (учет эффекта уширения, см. раздел 2.2). Здесь мы применили изложенную в разделе 2.2 методику для определения формы комплексов ДНК-ПАВ, что позволило сделать выводы о количестве молекул ДНК, входящих в комплекс.
Ниже изложена методика приготовления образцов для АСМ-исследования комплексов ДНК-ПАВ, апробированная совместно с к.х.н. О. А. Пышкиной и к.х.н. В. Г. Сергеевым (Химический факультет МГУ).
Процесс переноса комплексов осуществляли из водной среды в хлороформ. В качестве ПАВ нами был выбран дистеарилдиметиламмоний хлорид (ДСДАХ) ``Tokio Kasei Kogyo Co.'', Япония. В работе использовали ДНК из спермы лосося (``ГОСНИИОХТ'', Россия) 300-500 п. о.
Свободные молекулы ДНК нерастворимы в хлороформе и могут быть переведены в органическую фазу только в виде комплексов с ПАВ (которые, как упоминалось, в хлороформе растворяются). Контроль процесса переноса комплексов ДНК-ПАВ через границу раздела фаз проводили методом УФ-спектрофотометрии на спектрофотометре ``Specord M40'' (Германия). Было показано, что при добавлении 1,5-кратного избытка ДСДАХ к двухфазной системе, состоящей из водного раствора ДНК и фазы хлороформа, наблюдается практически 100% перенос молекул ДНК в органическую фазу.
Адсорбцию образцов на слюду осуществляли стандартным образом: каплю препарата (хлороформ с комплексами после проведения 70% переноса) объемом 3 мкл помещали на свежий скол слюды, выдерживали в течение 3-5 минут, промывали бидистиллированной водой и высушивали на воздухе, после чего осуществляли сканирование.
АСМ исследования осуществляли на приборе ``Nanoscope-IIIa'' (Digital Instruments, США) в режиме контакта с кантилеверами из нитрида кремния жесткостью 0,06 Н/м.
Методика нанесения комплексов ДНК-ПАВ на поверхность подложки из хлороформа, по-видимому, имеет следующие преимущества в сравнении с нанесением из водных растворов. Комплексы ДНК-ПАВ являются гидрофобными, и в водной среде могут взаимодействовать с силами притяжения за счет гидрофобного эффекта, что может привести к агрегации комплексов при приготовлении образцов и осложнить их АСМ-визуализацию. В хлороформе гидрофобные взаимодействия не проявляются, и перешедшие туда комплексы стабилизированы лишь внутренними механизмами конденсации ДНК, в силу чего каждый комплекс, по-видимому, состоит из молекул ДНК, покрытых снаружи слоем молекул ПАВ. Капля хлороформа, содержащая исследуемые комплексы, при нанесении на подложку высыхает достаточно быстро, что уменьшает вероятность образования регулярных структур вследствие механизмов агрегации или самоорганизации ПАВ.
С помощью АСМ было обнаружено, что при нанесении на подложку структур, оставшихся в водной фазы после проведения 70%-переноса, на АСМ-изображениях присутствуют агрегаты молекул ДНК в виде нитей длины 250-2000 нм. Образование длинных тяжей происходит, по всей видимости, вследствие межмолекулярной агрегации отдельных молекул ДНК при воздействии ПАВ. Измерение значений высоты и ширины на полувысоте АСМ-изображений данных структур дает значение 0,3-0,9 нм и 13-25 нм соответственно.
На рис. 4.3 приведены результаты АСМ-исследования комплексов ДНК-ПАВ, перешедших в хлороформ, адсорбированных затем на свежий скол слюды и исследованных на воздухе. С целью анализа морфологии комплексов была проведена статистическая обработка измеренных параметров АСМ-изображений тороподобных структур: D - диаметра тора, d - диаметра сечения тора и h - высоты тороидальной структуры над подложкой.
|
Гистограммы распределения параметров тороидальных комплексов представлены на рис. 4.4, статистика набрана в результате измерения 180 АСМ-профилей комплексов, сплошной линией указана аппроксимация полученных распределений по закону распределения Гаусса.
|
Средние значения для параметров тора, определенные по аппроксимации распределением Гаусса, составляют следующие величины (в качестве погрешности указано стандартное отклонение):
D = 100±30 нм, d = 25±9 нм и h = 5±2 нм.
Т.о. параметры частиц (внешний радиус тора Rout = (D+d)/2 = 62±16 нм) хорошо согласуются с типичными значениями для молекул ДНК, имеющих различную длину и сконденсированных различными агентами [94].
В силу того, что длина молекул лежит в диапазоне 100-170 нм, а средняя длина окружности тороидальной частицы pD составляет величину 300 нм, очевидно, что в образовании комплекса принимает участие, по крайней мере, две-три молекулы ДНК. Более подробный анализ требует восстановления реальной геометрии комплекса.
Занижение значений вертикальных размеров исследуемых структур обусловлено контактными деформациями под действием зонда, особенно существенными при исследовании объектов, радиус кривизны которых меньше радиуса кривизны кончика иглы, см. раздел 2.1. Предполагая, что деформация объекта не сопровождается изменением его объема, мы делаем вывод о том, что уменьшение высоты АСМ-изображения (занижение параметра h = 2b) будет скомпенсировано соответствующим увеличением его ширины (параметра d). Таким образом, мы делаем априорное предположение, что тороидальные комплексы ДНК-ПАВ имеют ``сплюснутую'' форму.
Завышение значений ширины АСМ-изображений является в основном следствием того, что кончик зондирующего острия имеет конечный радиус кривизны, см. раздел 2.2. Рассчитать вклад в уширение можно на основании анализа геометрии контакта исследуемой структуры и кончика зонда с помощью программы A.1.
Результаты применения методики учета эффекта уширения (раздел 2.2) к восстановлению ширины сечения исследуемых объектов.
Геометрию кончика иглы (радиус кривизны) можно определить при исследовании объектов, имеющих в сечении окружность известного радиуса, см. 3.1. Согласно результатам исследования тест-объектов (частиц вируса табачной мозаики, см. раздел 3.1) радиус кривизны кончика иглы кантилеверов используемого типа имеет величину от 5 до 15 нм (по данным фирмы-производителя до 40 нм, см. таблицу A.1). Определение точного значения радиуса R для конкретного зонда требует его тестирования непосредственно перед использованием. Однако и в этом случае существует вероятность того, что в процессе сканирования форма иглы претерпит изменения в результате взаимодействия с объектом. В этой связи чрезвычайно полезным является получение информации о форме зонда непосредственно из АСМ-изображений объекта исследования. Эта информация может быть получена в методике учета эффекта уширения (раздел 2.2).
Так, на рис. 4.5 представлен график зависимости числа случаев отсутствия решения (для статистики пар значений d и h) от радиуса аппроксимации зонда, построенный путем тестирования выполнения ``реперных'' условий (2.38) или (2.39). На основании рис. 4.5 можно сделать вывод, что верхняя граница значения R, характеризующего используемый зонд, составляет 12 нм - выше этого значения наблюдается линейный рост числа случаев отсутствия решения. Нижняя граница для радиуса кривизны зонда определяется контактными деформациями.
|
Восстановление истинной ширины тороидальных структур дает среднее значение параметра a: 10-11 нм (относительное отклонение e = 0,5) в моделии сферической иглы. Результаты применения параболической модели иглы дают значения 11-12 нм. Таким образом, действительно, формой комплекса ДНК-ПАВ является сплюснутый тор.
Расчет числа молекул ДНК, образующих комплекс с ПАВ.
Предполагая плотную (гексагональную) упаковку молекул ДНК в комплексе, получим, что объем, приходящийся на одну молекулу в составе комплекса, имеет величину 4 R2 sin60° L, где R - радиус молекулы ДНК, а L - ее длина.
Объем тороидального комплекса будет определяться формулой p2 abD, где a - найденное численное решение, b = h/2, D - диаметр тора. С учетом наличия слоя ПАВ толщиной l получим формулу для числа молекул:
|
На рис. 4.6 приведены гистограммы распределения найденных значений числа молекул при двух значениях параметра R (радиуса кривизны иглы) - 6 и 12 нм.
|
Анализ гистограмм показывает, что, согласно полученному решению, 80% комплексов имеет в своем составе менее 15 молекул ДНК (при R=12 нм).
В разделе 4.1 мы рассматривали конформационные изменения, которые претерпевают молекулы ДНК при взаимодействии с молекулами ПАВ в водных растворах. Методом АСМ было подтверждено, что в этих условиях может происходить компактизация молекул ДНК, при этом была визуализирована тороидальная структура компактных комплексов.
Стоит отметить, что система ДНК + ПАВ достаточно сложна для исследований прямыми методами визуализации геометрии микрообъектов, к которым относится и АСМ. Действительно, присутствие ПАВ в качестве компонента системы требует проведения трудоемких контрольных экспериментов для определения методики приготовления образцов, при использовании которой структуры, образованные только ПАВ (без ДНК), не приводили бы к ошибкам в интерпретации АСМ-результатов. Это связано со способностью ПАВ к самоорганизации [106], которая, например, может иметь место в процессе высыхания капли препарата, нанесенной на поверхность подложки для приготовления образцов (при высыхании капли будет повышаться концентрация ПАВ). Поэтому необходимо определить такие условия процессов нанесения образцов, высушивания капли препарата и промывки, при которых наличие на подложке объектов, связанных лишь с самоорганизацией ПАВ (без участия ДНК), минимально. Или, если они и присутствуют, то их геометрия принципиально отличается от исследуемых комплексов ДНК-ПАВ: например, наличие на подложке самоорганизованных плоских моно- или мультислойных пленок ПАВ не осложняет идентификацию искомых ДНК-ПАВ комплексов, но, напротив, часто способствует иммоблизации исследуемых объектов в процессе сканирования.
Поэтому достаточно показательным было бы подтвердить адекватность применения СЗМ для исследования процессов конденсации ДНК на примере более простых систем. Таковой системой, на наш взгляд, может считаться молекула ДНК, находящаяся в водно-спиртовой среде (этанол, изопропанол). При достаточно высоких концентрациях спирт становится для молекул ДНК плохим растворителем (в присутствии катионов Na+ - при работе с натриевыми солями ДНК), что приводит к переходу макромолекул из состояния развернутого клубка в состояние компактной глобулы [107].
Ниже изложены методики приготовления образцов для АСМ-исследования конденсации молекул ДНК, апробированные совместно с к.х.н. О. А. Пышкиной и к.х.н. В. Г. Сергеевым (Химический факультет МГУ).
Эксперименты в естественной атмосфере осуществляли по следующей схеме. Молекулы ДНК конденсировали в водно-спиртовой среде при концентрациях спирта, приводящих к переходу клубок ® глобула, и наносили каплю препарата на подложку, затем образцы высушивали, промывали и снова высушивали, т.е. схема эксперимента была стандартной для исследований на воздухе и идентичной случаю анализа конденсации комплексов ДНК-ПАВ.
Для экспериментов выбрали ДНК из эритроцитов цыплят (1000-3000 п. о., ``Союзхимреактив'', Россия ), в качестве спирта использовали изопропанол и этанол (х. ч.). Контрольные эксперименты по исследованию молекул ДНК, приготовленных нанесением из водного раствора буфера (без спирта), демонстрируют отсутствие компактизации макромолекул в данных условиях - рис. 4.7. Здесь мы использовали в качестве подложки немодифицированную слюду, поскольку такую же подложку выбирали и для конденсированных молекул. Использование немодифицированной (т.е. отрицательно заряженной) слюды в качестве подложки объясняет тот факт, что адсорбированные молекулы недостаточно расправлены и частично проагрегировали.
|
Добавление к исходному раствору 50% и более спирта (изопропанола или этанола) существенно влияет на конформацию ДНК. Результаты экспериментов показывают, что большая часть молекул находится в одной из двух типов конформаций - тороидальной или стержневой (rod-like). Следует подчеркнуть, что при высокой концентрации молекул ДНК в исходном растворе тороидальные и стержневые частицы агрегируют в единую структуру. На рис. 4.8 представлен результат АСМ визуализации стержневой структуры ДНК, полученной при адсорбции макромолекул из 66% водного раствора этанола (концентрация ДНК составляла 5×10-5 М).
|
Анализ ширины профиля стержневых структур показывает, что они образованы, по крайней мере, несколькими молекулами ДНК (``стянутыми в жгут''); высокая степень агрегации стержневых структур объясняется значительной концентрацией ДНК в рабочем растворе.
Несколько отличные результаты были получены при использовании изопропанола. На рис. 4.9 представлены результаты АСМ-визуализации тороидальных структур, образованных молекулами ДНК при нанесении на поверхность слюды 4×10-7 М раствора ДНК в 50% изопропаноле (в присутствии 0,05 М ТВЕ буфера). Результаты свидетельствуют, что в этом случае (в сравнении с этанолом) более часто наблюдаются тороподобные структуры, нежели стержневые. Причина может заключаться в различных значениях диэлектрических проницаемостей изопропанола (e = 18,3) и этанола (e = 24,3), что обеспечивает лучшую компенсацию заряда поверхности макромолекулы противоионами Na+ в случае изопропанола.
|
Стоит подчеркнуть, что при проведении АСМ-исследований на воздухе необходимой стадией приготовления образцов является процедура высушивания капли препарата на подложке. В процессе высушивания концентрации компонент препарата в капле неравновесно возрастают, что может неконтролируемо влиять на морфологию адсорбируемых на подложку структур. Исключить процедуру высушивания возможно при проведении АСМ-исследований в жидких средах.
В этой серии экспериментов проводили исследования конденсации молекул ДНК непосредственно в жидкостной ячейке АСМ, что позволяло исключить возможное негативное влияние процесса высушивания капли препарата в стандартной схеме экспериментов на воздухе.
Использовали ДНК бактериофага Т4 и водные растворы изопропанола и этанола. Ранее было показано [107], что гигантская ДНК Т4 в водно-спиртовой смеси при концентрации спирта более 50% (для изопропанола) претерпевает конформационный переход клубок ® глобула. При уменьшении концентрации спирта до 40% наблюдалось частичное разворачивание глобул. Динамика перехода исследовалась методом флуоресцентной микроскопии. Однако, в силу ограничения разрешающей способности метода (предел оптического микроскопа), оказалось трудным определить точную геометрию компактных глобул и структуру частично компактизованных образований.
Представлялось показательным провести АСМ-исследования процесса компактизации именно для этого объекта в тех же условиях.
В качестве объектов исследования использовали ДНК бактериофага Т4 (``Nippone Gen'', Япония), состоящей из 166 тыс. п.о. При проведении экспериментов использовали бидистиллированную воду, очищенную с помощью системы ``Milli-Q'' (``Millipore'', США). Водно-спиртовые среды готовили с использованием изопропанола (х. ч.).
Исследования конденсации молекул ДНК проводили методом атомно-силовой микроскопии в водно-спиртовой смеси с использованием подложек слюды. Свежий скол слюды и кантилевер (с зондом) помещали в жидкостную ячейку АСМ, которую последовательно заполняли рабочими растворами. Препарат, содержащий молекулы ДНК Т4 в водно-спиртовой смеси при 80% концентрации спирта, получали, смешивая один объем раствора ДНК в 0,5-ТВЕ-буфере с четырьмя объемами изопропанола. Препарат молекул ДНК в 40% изопропаноле получали, смешивая два объема раствора ДНК в буфере с тремя объемами изопропанола. В обоих случаях конечная концентрация ДНК в смеси составляла 1×10-6 М.
Эксперименты проводили в жидкостной ячейке АСМ ``Nanoscope-III'' (Digital Instruments, США) в режиме прерывистого контакта с использованием заостренных кантилеверов из нитрида кремния жесткостью 0,32 Н/м, см. таблицу A.1.
Определение адекватной методики проведения АСМ-эксперимента включало анализ наиболее подходящих концентраций водно-спиртовой смеси, а также последовательность заполнения жидкостной ячейки рабочими растворами (тестировали различные варианты).
В результате тестирования остановились на следующей методике:
Препарат, содержащий молекулы ДНК Т4 в водно-спиртовой смеси при 80% концентрации изопропанола, вводили в жидкостную ячейку, заполненную предварительно водно-спиртовой смесью с той же концентрацией спирта. При этих условиях макромолекулы ДНК осаждались на поверхность слюды в виде компактных глобул (рис. 4.10а и в).
| , Figure |
Применение методики восстановления реальных размеров исследуемых объектов по измеренной высоте и ширине профиля АСМ-изображения (учет эффекта уширения), которая изложена в разделе 2.2, позволяет определить объем глобулярных структур. Проводили расчеты эффекта уширения для двух значений радиуса кривизны иглы: 5 и 10 нм, результаты приведены в таблице 4.1. Согласно оценкам (с применением тест-объектов) радиус кривизны иглы, используемой при записи приведенных изображений, лежит именно в диапазоне 5ё10 нм. Значения параметров глобулы, определенные с использованием двух граничных значений для радиуса кривизны иглы фактически идентичны. Восстановленные значения параметров полуосей a и b позволяют сделать вывод, что геометрической формой глобулы является сплюснутый и слегка вытянутый эллипсоид.
Стандартные отклонения величин, приведенных в таблице 4.1, велики, что объясняется статистическим разбросом измеряемых параметров для различных глобул. Более наглядно это видно из анализа гистограмм распределения анализируемых параметров.
R, нм | a, нм | b, нм | c, нм | V, ×105 нм3 |
5 | 41±12 | 55±17 | 23±9 | 2,5±1,7 |
10 | 37±12 | 52±17 | 23±9 | 2,1±1,5 |
| ,
Figure |
На рисунке 4.11 представлены гистограммы распределения значений объема глобулярных частиц. Наличие в набранной статистике глобул, характеризующихся малыми значениями объема (меньше, чем объем одной молекулы ДНК) объясняется, по-видимому, присутствием в используемом препарате фрагментов молекул ДНК. Из анализа рисунка 4.11 (ср. с таблицей 4.1) следует, что объем глобулы, образованной цельной молекулой ДНК (не фрагментом) составляет величину порядка 3-5×105 нм3, что превышает величину объема молекулы ДНК Т4: 1,7×105 нм3.
Т.о. образом можно предположить, что каждая глобула образована одной молекулой ДНК, находящейся в достаточно плотно упакованном состоянии.
Частично декомпактизованные структуры макромолекул ДНК Т4, получаемые при разбавлении 80% спиртового раствора в воде до 40-50%, исследовали следующим образом. В жидкостную ячейку, содержащую ДНК Т4 в виде глобул в 80% изопропаноле (20% воды), добавляли смесь изопропанола с водой (40% изопропанола + 60% воды), с целью понизить концентрацию спирта в ячейке.
Поскольку процесс компактизации ДНК в спирте обратим, то при понижении концентрации спирта в жидкостной ячейке глобулы претерпевают конформационные изменения, связанные с началом процесса декомпактизации. На рисунках 4.10 б и г изображены те же участки поверхности, что и на рисунках 4.10 а и в, но при понижении концентрации спирта в ячейке до значения 40-50%. Из рисунков видно, что в процессе декомпактизации изменились очертания глобулярных структур и их высота (ср. а и б).
Некоторые выводы о структуре компактизованной глобулы можно сделать, анализируя стадии динамического процесса декомпактизации, изображенные на рис. 4.10 в и г. Анализ рисунка показывает, что первоначально компактная глобула (рис. 4.10 в, отметка ``Гл'' ) в процессе декомпактизации уменьшается в высоте и позволяет визуализировать центральную полость, см. рис. 4.10 г и врезку.
Однако оказалось, что дальнейший процесс разворачивания глобул не происходит. Это можно объяснить наличием сил взаимодействия между молекулой и подложкой, которые и препятствуют разворачиванию при уменьшении концентрации спирта. (Оказалось, что в условиях водно-спиртовой среды исследуемые объекты достаточно стабильны на подложке и не увлекаются зондом при сканировании).
Поэтому промежуточные состояния компактизации ДНК (между глобулой и клубком) получали следующим образом. Молекулы ДНК, находящиеся в 40% изопропаноле, вводили в жидкостную ячейку, а затем концентрацию спирта в ячейке повышали прокачиванием 80% изопропанола. При повышении концентрации спирта молекулы выпадали на поверхность подложки в частично компактизованном состоянии. Дальнейшая компактизация ДНК, адсорбированных на подложку, не наблюдалась, что, по-видимому, опять же объясняется взаимодействием молекул ДНК со слюдой.
АСМ-изображения частично компактизованных структур приведены на рис. 4.12. Было обнаружено, что начальным процессом компактизации глобул является закручивание отдельных участков макромолекулы ДНК в тороидальную структуру, эти компактные участки, по-видимому, и являются центрами дальнейшей компактизации.
| ,
Figure |
На основании наблюдений можно предположить, что и компактные глобулы, визуализованные на рисунке 4.10 а и в являются продуктом именно тороидальной компактизации макромолекул ДНК и представляют собой плотные ``клубки'', в которых отдельные участки молекулы закручены вокруг тороидальных структур.
Следует отметить, что макромолекулы ДНК при повышении концентрации спирта в ячейке (от 40%) образуют не только тороидальные, но и (реже) стержнеобразные структуры, см. рисунок 4.13.
| |
Стрежневые структуры, изображенные на рисунке 4.13, имеют ту же природу, что и приведенные на рисунке 4.8, где изображены результаты исследований, проводимых на воздухе. Видно, что при исследованиях в жидкости достижимо большее пространственное разрешение, что позволяет визуализировать отдельные нити молекулы ДНК, образующей исследуемую структуру. Из рисунка 4.13 б видно, что при сканировании под воздействием зонда микроскопа стержневая структура частично ``расплелась''.
Наблюдение в ряде случаев компактизованных структур, имеющих морфологию, отличную от тороидальной, может объясняться тем, что в эксперименте исследовались промежуточные стадии процесса компактизации, в которых некоторое количество молекул может находиться в нестационарных морфологических состояниях
Описанные подходы проведения АСМ-эксперимента в жидких средах (т.е., потенциально, в условиях, близких нативным) могут быть полезны при исследованиях модельных систем транспорта генетического материала внутрь живых клеток. Подобные исследования могут также дать вклад в понимание механизмов процессов компактизации молекулы ДНК под воздействием различных компактизующих реагентов.
Ленгмюровская пленка [108] может быть определена как упорядоченный нерастворимый мономолекулярный слой, сформированный на границе фаз жидкость - газ. Пленкообразующее вещество должно быть амфифильным, т.е. состоять из компактной полярной головки (гидрофильной) и гидрофобного фрагмента (длинного алкильного радикала). Процесс формирования ленгмюровской пленки состоит в следующем. Вначале, при нанесении на поверхность раздела фаз пленкообразующего вещества в подходящем растворителе, площадь поверхности пленки сохраняется достаточно большой, чтобы после испарения растворителя с поверхности жидкости образовалась двумерная газообразная пленка. При сокращении площади, занимаемой пленкой, с помощью подвижного барьера происходит ориентация и упорядочивание молекул в пленке и она последовательно проходит различные двумерные фазовые состояния: от ``газообразного'' до ``твердого'' [20], что сопровождается изменением характера зависимости поверхностного давления от площади, приходящейся на молекулу p = p(A); эти зависимости являются аналогом изотерм p = p(V) для трехмерного случая [102]. Однако фазовые диаграммы двумерного случая носят более сложный характер, так, например, существуют двумерные фазовые состояния ``жидко-растянутая'', ``жидко-конденсированная пленка'' [109] (которые различаются по характеру зависимости p = p(A) на этих участках), не имеющие аналогов на изотермах p = p(V) трехмерного случая.
Пленками Ленгмюра-Блоджетт (ЛБ) называют совокупность ленгмюровских пленок на твердой поверхности [110]. Перенос пленки осуществляют обычно вертикальным способом (метод Ленгмюра-Блоджетт), пропуская подложку сквозь монослой, причем давление в пленке поддерживают постоянным путем сокращения площади пленки на водной поверхности подвижным барьером в процессе переноса монослоя. Структура сформированного покрытия определяется способом выделения (рис. 5.1).
|
Технология ЛБ позволяет достаточно просто изменять свойства поверхности и формировать качественные пленочные покрытия за счет точного контроля толщины пленки в процессе выделения (количества наносимых слоев), однородности покрытия, низкой шероховатости и высокой (при подборе адекватных условий) адгезии пленки к поверхности. Свойства пленок также можно легко варьировать, изменяя структуру полярной головки амфифильной молекулы, состав монослоя (двух- и многокомпоненные смеси молекул), а также варьируя условия выделения - состав субфазы и поверхностное давление. Это объясняет интерес к потенциальному использованию моно- и мультислойных пленок ЛБ в высокотехнологичных отраслях: микро- и наноэлектронике, оптоэлектронике, при разработке сенсоров и биосенсоров, био- и ультрафильтрационных мембран, создании нелинейнооптических материалов и оптических элементов (волноводы, люминесцентные устройства) с заданными свойствами, и т.д. Однако для широкомасштабного применения пленок данного типа необходимо прояснить основные механизмы, определяющие структурные и физические свойства сформированных тонкопленочных покрытий.
Морфология пленки на твердой поверхности зависит от совокупности факторов - условий формирования монослоя на межфазной поверхности, способа и условий переноса пленки на твердую поверхность, природы подложки и т.д. Весьма критичным для качественного переноса пленки является структура поверхностно-активного вещества и, в частности, природа полярной головки, которая в значительной мере определяет свойства монослоя. Действительно, многие амфифильные молекулы образуют стабильные монослои на водной поверхности, однако не все материалы могут быть успешно выделены в виде пленок ЛБ традиционным (вертикальным) методом выделения.
Выделение монослоев на твердую поверхность и формирование мультимолекулярных пленок вертикальным методом ЛБ происходит при циклическом погружении и подъеме подложки сквозь монослой, причем давление в пленке поддерживают постоянным в процессе переноса (рис. 5.1а). В этом случае формируется пленка (называемая пленкой Y-типа), состоящая из центросимметричных бислоев. В зависимости от строения, некоторые молекулы могут выделяться в виде монослоя только при погружении (Х-тип выделения) или только при извлечении подложки из субфазы (Z-тип). При этом последовательно нанесенные монослои не обязательно обладают фиксированной ориентацией. Известно [20] (по данным рентгеноструктурного анализа), что многие пленки X-типа в результате переорганизации имеют внутреннюю ориентацию, совпадающую с пленками Y-типа (считается, что Y-пленки более стабильны). Чтобы избежать явлений, затрудняющих перенос монослоя обычным вертикальным методом выделения, таких, как опрокидывание молекул в процессе переноса, кристаллизация вещества в мениске жидкости и т.д., используют метод горизонтального осаждения (ГО), см. рис. 5.1б.
Существенное влияние на свойства монослоев оказывает состав субфазы и природа противоиона. Наиболее ярко это проявляется на примере поверхностно-активных кислот. Перенос на твердую подложку монослоев жирных высших кислот, сформированных на чистой водной поверхности, весьма проблематичен. При использовании подложки из слюды это связано с отрицательным зарядом, приобретаемым слюдой в водных растворах (карбоксильная группа также заряжается отрицательно), т.е. возникает та же проблема, что и при иммобилизации на поверхности подложки молекул нуклеиновых кислот. Поэтому, обычно, монослои жирных кислот переносят на подложку методом ЛБ в виде их солей с ионами различных металлов [111]. Однако методом горизонтального осаждения удается переносить качественную монослойную пленку бегеновой кислоты и с водной поверхности на слюду, см. рис. 5.2а. На этом рисунке отчетливо видно сосуществование двух типов областей пленки, характеризующихся различной толщиной слоя, а, следовательно, и различной ориентацией углеводородных хвостов (различный угол с направлением нормали к поверхности). По-видимому, наблюдаемые области являются участками различного двумерного фазового состояния пленки (``твердая'' и ``жидко-конденсированная'' фазы). При использовании традиционного вертикального метода ЛБ формирования покрытия детектировать подобные особенности организации монослоя на твердой поверхности (сосуществование фаз различной структуры) не удается.
Актуальной задачей наноэлектроники является создание гибридных органическо-неорганических наноструктур на твердой поверхности. Полупроводниковые кластеры представляют собой один из наиболее интересных и перспективных объектов исследования; их уникальные свойства обусловлены квантоворазмерными эффектами и проявляются, например, во взаимодействии с лазерным излучением: процессах рассеяния и поглощения, генерации второй гармоники и пр. Отличительной особенностью нанокомплексов на основе халькогенидов Cd, стабилизированных органическими лигандами, является наличие реакционноспособных так называемых Е-центров, центров свечения, (например, атомов Te в CdTe кластерах) на внешней поверхности нанокластеров. Особенности кластеров таковы, что большинство поверхностных атомов Cd связано с органическими лигандами, в то время как Е-центры остаются несвязанными и влияют на процесс взаимодействия кластеров с электромагнитным излучением. Помимо этого, поскольку Е-центры свободны, существует возможность их химической модификации и создания новых структур [112].
Была изучена возможность введения CdTe кластеров в монослой с последующим переносом наночастиц на поверхность на примере Cd54Te32(SCH2CH2OH)x » 50. Исследованные нанокластеры имеют предположительно тетраэдрическое строение, причем тиоглицерол-стабилизирующие лиганды располагаются на углах и ребрах, оставляя открытыми для химического контакта атомы Te на плоских гранях тетраэдра [112]. С помощью методов кварцевого микровзвешивания, ИК- и УФ-спектроскопии было установлено, что тиоглицерол-стабилизированные кластеры можно успешно внедрять в монослой бегеновой кислоты и переносить на твердую подложку с поверхности их водных растворов.
Ниже изложены методики формирования тонкопленочных покрытий на твердой подложке для АСМ-исследований, апробированные совместно с к.х.н. Г. К. Жавнерко (Институт новых материалов Академии Наук Белоруси).
Образцы пленок были приготовлены вертикальным методом, см. рис. 5.1а, и методом горизонтального осаждения (рис. 5.1б), который в ряде случаев позволяет получать более качественные монослойные покрытия в сравнении с традиционным вертикальным способом. Действительно, сформированный на водной поверхности монослой осаждается методом ГО ``как есть'' на предварительно погруженную под водную поверхность подложку в процессе медленного (около 0,01 мм в минуту) понижения уровня жидкости (в результате формируется структура Z-типа, рис. 5.1б). Пленки, полученные методом ЛБ при пропускании подложки сквозь монослой вертикально со скоростью 4 мм/мин, имели, как правило, структуру Y-типа.
Тиоглицерол-стабилизированные CdTe кластеры, использованные для внедрения в пленки, были приготовлены по методике, описанной в работе [112]. Монослой бегеновой кислоты формировали на поверхности водных растворов CdTe кластеров с концентрацией 0,1-1×10-4 М. Монослой Z-типа или 1-Y бислой бегеновой кислоты выделяли с водной поверхности или с субфазы, содержащей CdTe кластеры на поверхность подложки (слюды или пирографита), очищенную стандартным методом межслойного скалывания непосредственно перед нанесением монослоя.
АСМ исследования проводили на приборе Nanoscope-IIIa (Digital Instruments, США) с использованием кантилеверов из нитрида кремния жесткостью 0,32 Н/м с заостренными зондами в режиме контакта на воздухе. Было обнаружено, что использование данного контактного кантилевера (см. таблицу A.1) предпочтительнее в сопоставлении с другими того же типа, поставлямыми фирмой-производителем (жесткостью 0,06 Н/м, 0,12 Н/м или 0,58 Н/м), при исследованиях тонких пленок, поскольку разрушающее воздействие зонда на образец в этом случае минимально. Это может быть связано с тем, что данный кантилевер обладает наименьшей массой, и, как следствие, наименее инерционен. Силу воздействия при сканировании варьировали в широком диапазоне: от единиц до десятков и сотен наноньютонов, при этом в ряде случаев возможность проведения неразрушающих исследований требовала минимизации силы воздействия зонда. Частоту строчной развертки выбирали в диапазоне 5-60 Гц; при получении молекулярного разрешения для минимизации теплового дрейфа использовали максимальное значение (60 Гц).
Было обнаружено, что осажденная с CdTe субфазы на поверхность подложки (метод ГО) мономолекулярная пленка бегеновой кислоты Z-типа содержит спонтанно сформированные кристаллические образования, имеющие ламелярную структуру. Близкую морфологию имела пленка, сформированная по процедуре получения 1-Y бислоя (метод ЛБ). Анализ сечения АСМ-изображений этих образований показывает, что их высота над поверхностью слюды близка утроенной длине молекулы бегеновой кислоты в максимально вытянутой конформации (около 3 нм), см. рис 5.3. В случае использования пирографитовой подложки и следования процедуре формирования 1-Y бислоя (метод ЛБ) структура покрытия, фактически, идентична: также наблюдаются включения аналогичных кристаллических структур, см. рис. 5.4, но морфология покрытия, окружающего кристаллиты, менее однородна. Согласно результатам АСМ-исследования высота монослоя над поверхностью подложки (слюды) варьируется в диапазоне 2,2-2,7 нм (занижение высот обусловлено контактными деформациями, см. раздел 2.1). Высота кристаллитов над поверхностью монослоя (для слюды) лежит в интервале 4,2-6,2 нм. Высота кристаллитов над поверхностью пирографитовой подложки приблизительно равна 8 нм (т.е. соотносится со значениями для подложки-слюды). Поскольку, как показали контрольные эксперименты, подобные кристаллиты не формируются при нанесении аналогичного монослоя с поверхности воды или CdCl2 субфазы, мы полагаем, что частичная кристаллизация бегеновой пленки при нанесении с поверхности CdTe субфазы обусловлена влиянием именно кластеров, которые могут являться центрами начальной кристаллизации.
| ,
Figure, |
| ,
Figure |
| ,
Figure |
Предположение о том, что CdTe кластеры находятся именно внутри кристаллитов, подтверждается тем фактом, что спектроскопические методы исследования детектируют присутствие кластеров в приготовленном покрытии; причем их детектируемое количество для конкретного образца напрямую коррелирует с величиной поверхностной концентрации кристаллических структур, визуализованных АСМ.
Было установлено, что ориентация ламелей кристаллических структур подчиняется определенным закономерностям. Так, углы между ламелями кратны 60°, как при использовании подложки из слюды, так и графита (рис. 5.5).
Гистограммы показывают, что среднее значение величины угла между ламелями составляет 60±7°, как в случае подложки из слюды (рис. 5.5а), так и пирографита (рис. 5.5б). Возможную погрешность измерения, связанную с наличием теплового дрейфа, мы оценивали путем контрольного измерения углов между кристаллографическими осями подложки (угол 60°). Было показано, что в условиях данного эксперимента анализируемая погрешность составляет величину меньшую, чем 3°.
Для установления ориентации ламелей относительно кристаллографических осей подложки использовали следующую процедуру. Определяли ориентацию ламелей. Затем в монослойном покрытии формировали искусственный дефект путем воздействия на эту область зонда микроскопа с увеличенной силой (до 100 нН и более) при медленной скорости сканирования (частота строчной развертки менее 3 Гц). В результате материал пленки удаляли зондом микроскопа в пределах выбранного окна сканирования. В этом окне получали АСМ-изображение атомной решетки подложки и определяли направление основных поверхностных кристаллографических осей.
Как показали эксперименты по АСМ-визуализации молекулярной структуры поверхности кристаллитов (см. ниже раздел 5.3) вектор a поверхностной решетки (ЦП-тип) сонаправлен с ориентацией ламели. Оказалось, что для подложки из слюды направление конкретной ламели кристаллической структуры образует с одной из основных осей решетки подложки ([10], [01] или [11]) угол близкий 20°. На рис. 5.6 а представлена гистограмма распределения значений угла между направлением ламели кристаллической структуры и одной из кристаллографических осей слюды (ближайшей, к направлению ламели).
Из гистограммы следует, что среднее значение угла между кристаллографической осью поверхностной решетки слюды и направлением ориентации ламелей кристаллической структуры составляет 20 ±4°. Это позволяет предположить, что ламели предпочтительно ориентированы вдоль направлений [2[`(1)]], [32], [13] или [31], [23], [[`(1)]2] решетки слюды. Действительно, углы между направлениями [31] и [10], [23] и [11], [[`(1)]2] и [01] имеют величину 19,1°, а углы между [2[`(1)]] и [10], [32] и [11], [13] и [01] составляют ту же величину, но со знаком минус.
Те же измерения проводили для пленок, сформированных на подложке из пирографита. Измерение угла между направлениями кристаллических ламелей и кристаллографическими осями графита дает значение 30±7°, гистограмма распределения результатов измерения приведена на рисунке 5.6б. Полученное значение позволяет предположить, что предпочтительным направлением ориентации ламелей в этом случае является направление [12] поверхностной решетки графита.
Несмотря на локальную устойчивость кристаллических образований к силовому воздействию зонда на небольшой участок поверхности, именно они (а не монослой) в первую очередь разрушаются в случае, когда величина воздействующей силы увеличена (до 100 нН), а размер кадра составляет значительную величину (до 15×15 мкм2). После нескольких сканирований с увеличенной силой воздействия зонда на АСМ-изображении наблюдается следующая картина: участки пленки, соответствовавшие монослойным участкам, характеризуются лишь увеличенным количеством пор и мелких дефектов, в то время как кристаллические образования ``исчезают''. Анализ высот АСМ-изображения показывает, что в тех местах, где прежде были кристаллические структуры, остается лишь чистая поверхность слюды. Визуально это наблюдается в виде дефектов (дыр) в монослое, контуры которых в точности повторяют контуры кристаллических структур, присутствовавших на этих местах прежде, см. рис. 5.3б.
Данный эффект может свидетельствовать о склонности кристаллических структур взаимодействовать с иглой микроскопа. Причем, при последующих сканированиях с большим размером кадра и минимизованной силой, на краях области сканирования наблюдалось лишь незначительное количество удаленного материала, что позволяет сделать вывод о его адсорбции на зонде. Действительно, площадь боковых граней зонда составляет несколько микрон (см. таблицу A.1), поэтому значительное количество удаленного материала пленки может при сканировании некоторое время ``облеплять'' боковые стенки зонда, несущественно влияя на качество получаемых АСМ-изображений.
Если после разрушения кристаллических структур увеличенной силой продолжать сканирование, уменьшив силу воздействия зонда, то через несколько последовательных сканирований материал вновь появится в тех участках, где прежде были кристаллические структуры, см. рис. 5.3 (в и г). При этом восстанавливается и высота кристаллических образований, что свидетельствует об их значительной устойчивости. Сходный эффект переноса материала пленки с подложки на зонд и обратно был описан в работе [62].
Т.о. визуализирован динамический и обратимый процесс разрушения кристаллических структур, что позволяет сделать ряд заключений.
При исследовании пленки бегеновой кислоты на пирографите наблюдали, что материал пленки организовывался в упорядоченные протяженные структуры (см. рис. 5.7), похожие на ``грядки''. Ориентация и организация грядок подчиняется определенным закономерностям: углы между ними кратны 60°, а расстояние между грядками составляет величину 6,3±0,4 нм. Угол между направлениями грядок и основными кристаллографическими осями поверхностной решетки пирографита составляет 30° (рис. 5.8), т.о. грядки ориентированы вдоль того же направления подложки пирографита, что и описанные выше кристаллические ламели.
| ,
Figure, |
| ,
Figure |
|
Было обнаружено, что сформированное покрытие на графите характеризуется нестабильностью и разрушается иглой при сканировании при незначительном увеличении воздействующей силы (или увеличении характерного времени воздействия t при уменьшении размеров кадра). При исследовании динамики разрушения покрытия был обнаружен эффект временного формирования сходных упорядоченных структур (``грядок''), характеризующихся той же фиксированной взаимоориентацией (угол между этими структурами также близок 60°) и тем же фиксированным расстоянием между грядками (около 6,3 нм). Было обнаружено, что ориентация описываемых структур также составляет угол близкий 30° с направлениями кристаллографических осей графита, что свидетельствует о выделенности данного направления для эффекта взаимодействия молекул бегеновой кислоты с пирографитовой подложкой. Следовательно, можно сделать вывод, что в обоих случаях формирование грядок обусловлено одним механизмом.
| ,
Figure, |
| ,
Figure |
Наблюдаемые структуры, по-видимому, являются полуцилиндрическими мицеллами, формируемыми на поверхности гидрофобной пирографитовой подложки молекулами бегеновой кислоты. Авторы работ [113,114] обнаружили сходные протяженные образования, сформированные молекулами ПАВ на поверхностях слюды и пирографита. Так в работе [114] исследовали процесс самоорганизации молекул катионного ПАВ: тетрадецилтриметиламмоний бромида (C14TAB), которые, как показали авторы, формировали на поверхности подложки цилиндрические (в случае подложки из слюды) и полуцилиндрические (при использовании пирографита) образования. В случае полуцилиндров они были образованы молекулами ПАВ, стабилизированными гидрофобным взаимодействием с подложкой и друг с другом (поэтому полярные группы молекул формировали внешнюю поверхность полуцилиндра). При этом направления ориентации ``грядок'' на поверхности пирографита, как и в нашем случае, составляла угол в 30° с направлением основных кристаллографических осей подложки. Этот эффект может объясняться тем, что в случае такой ориентации полуцилиндрических мицелл отдельные молекулы ПАВ, лежащие на подложке, ориентированы вдоль кристаллографических направлений графита, причем периодичность решетки графита вдоль этих направлений (0,246 нм) и период углеводородной цепи молекулы (0,254 нм) имеют близкие значения.
Период упаковки полуцилиндрических мицелл, обнаруженный авторами рассматриваемой работы, составил 4,7±0,3 нм. Длина углеводородных цепей молекул бегеновой кислоты (C22), объекта исследования в нашем случае, в 1,5 раза больше, чем молекул C14TAB, исследованных в [114]. Поэтому наблюдаемый период ``грядок'', формируемых молекулами бегеновой кислоты (6,3±0,4 нм) вполне соотносится со значением, полученным в упомянутой работе. Эти закономерности позволяют сделать вывод, что особенности самоорганизации молекул бегеновой кислоты на поверхности пирографита объясняются именно формированием полуцилиндрических мицелл.
Иммобилизация белковых молекул на твердой поверхности достаточно актуальна как для практических задач (создание био- и имунносенсоров, например), так и для фундаментальных исследований строения белка [20]. При решении этих задач существует проблема стабилизации белков на поверхности подложки [27]. Ввиду большого разнообразия белковых молекул, чрезвычайно сложно найти универсальную методику их жесткой фиксации на подложке. Основные белки можно рассматривать как амфифильные макромолекулы [115], поскольку они имеют гидрофильные и гидрофобные зоны. Следовательно, можно использовать все достоинства технологии ЛБ для создания однородных монослойных пленок белков на водной поверхности и переноса их на твердую подложку. На примере ряда белков (гемоглобин лошади, бычий сывороточный альбумин (БСА) и др.) была исследована возможность формирования на различных поверхностях как индивидуальных, так и композиционных пленок.
Ниже изложена методика формирования белковых покрытий на твердой подложке для АСМ-исследований, апробированная совместно с к.х.н. Г. К. Жавнерко (Институт новых материалов Академии Наук Белоруси).
Монослои белковых молекул на водной поверхности формировали, используя различные варианты доставки вещества на поверхность. Во-первых, использовали тот факт, что исследуемые белки самопроизвольно концентрируются на водной поверхности после введения их растворов в водную субфазу (рН 5,8). В этом случае сжимали монослой спустя примерно через 30 минут после нанесения вещества на поверхность. Во-вторых, белок непосредственно на водную поверхность наносили из спиртового раствора, который готовили, растворяя 1-2 капли водного раствора белка, содержащего необходимую навеску, в 0,5 мл этилового спирта. В третьих, на водную поверхность наносили раствор белка, приготовленный в смеси хлороформ: ацетон: спирт в объемном соотношении 1:1:1 (ХАС раствор).
Особое внимание уделили гемоглобину. Гемоглобин, как основной компонент эритроцитов крови, осуществляющий доставку кислорода от легких к тканям, достаточно доступен. При растворении гемоглобина в ХАС растворе спустя 2-3 часа белок выделял гем (окрашенный в розовый цвет раствор), а в осадок выпадал белок белого цвета. Оказалось, что осадок опять возможно перевести в водный раствор (рН 5,1-5,3), как и при растворении исходного гемоглобина. Как белок, выпавший в осадок (глобин), так и геминовая часть обладали поверхностно-активными свойствами, что позволило формировать ЛБ покрытия как из белков без геминовой группы, так и, наоборот, только из геминовых групп.
Пленки (монослои) были перенесены на поверхности подложек (слюда, кремний, пирографит) с поверхности водной субфазы методом ГО при давлении выделения p в диапазоне 20ё25 мН/м.
АСМ-измерения проводили на приборе ``Nanoscope-IIIa'' в режимах контакта и модуляции силы (см. раздел 1.3.2) с использованием стандартных кантилеверов (таблица A.1) жесткостью 0,32 Н/м.
Сопоставление морфологии поверхности пленки гемоглобина, ковалентно сшитой (с помощью глутарового альдегида) с гидрофильной поверхностью, обработанной аминопропилтриэтоксисиланом, и монослойной пленки гемоглобина, осажденной с водной поверхности методом ГО на подложку, показывает преимущества второго способа формирования пленки. Если в первом случае имела место хаотичная агрегация белка из раствора на поверхности, то в случае осаждения монослоя белка с водной субфазы образуется однородное монослойное покрытие (рис. 5.10). Однородную пленку молекул гемоглобина удается сформировать на водной поверхности и перенести на поверхность подложки при различном способе доставки вещества на водную субфазу: как из спиртового раствора (рис. 5.10а), б)), так и из ХАС раствора (рис. 5.10в), г)).
Оказалось возможным сформировать пленки из разных фрагментов гемоглобина. Так, и на поверхность графита (рис. 5.10з), и)), и на поверхность слюды (рис. 5.10к)) методом ГО переносится качественная пленка из нанесенного на водную поверхность окрашенного (геминовой группой) ХАС раствора (после удаления из этого раствора выпавшего белого осадка белка). Вещество, оставшееся в ХАС растворе, во-первых, оказалось поверхностно-активным, что позволило сформировать пленку на водной поверхности, во-вторых, эта пленка, перенесенная на подложку, оказалась значительно более устойчива на поверхности графита по сравнению с поверхностью слюды. Отметим также, что после расщепления гемоглобина белковая часть молекулы также остается поверхностно-активной, поскольку и из нее можно формировать твердую пленку на водной поверхности и успешно переносить ее на гидрофильный кремний (рис. 5.10л)) или слюду (рис. 5.10м)).
Гемоглобин растворяется в достаточно узком интервале рН, поэтому можно было предположить, что его переход на водную поверхность из объема обусловлен изменением растворимости при изменении рН. В то же время бычий сывороточный альбумин растворим в достаточно широком интервале рН. Тем не менее, также удалось сформировать твердую пленку БСА на водной поверхности при доставке белка на субфазу из водного раствора: пленка успешно была перенесена на твердую подложку (рис. 5.10ж)). Аналогичная однородная монослойная пленка формируется и при использовании процедуры доставки на водную поверхность альбумина из спиртового раствора (рис. 5.10д), е)).
Искусственные дефекты в белковых пленках. Высокое качество пленки может вызвать неоднозначность в интерпретации результата АСМ-исследования, поскольку отсутствие каких-либо характерных дефектов (пор и т.п.) не позволяет оценить толщину пленки, и вообще сделать вывод о присутствии пленки на подложке. Действительно, однородная поверхность может наблюдаться и в том случае, если пленку не удалось успешно перенести с границы раздела фаз на твердую подложку (для ряда материалов достаточно сложно определить адекватные условия переноса: состав субфазы, поверхностное давление и пр.).
Полученные белковые пленки характеризовались высокой однородностью поверхности, поэтому проводили их дополнительный анализ, путем воздействия сканирующего зонда на некоторый выбранный участок при увеличенной величине силы (до 100 нН) и медленной скорости сканирования (частота строчной развертки около 1 Гц). Данная процедура, в случае наличия на подложке органической пленки, позволяла формировать в ней искусственный дефект квадратной формы, в котором материал пленки был удален зондом при сканировании. При указанных условиях эта процедура не приводит еще к разрушению подложки и позволяет проводить как контроль наличия пленки на поверхности, так и измерение ее толщины (следует, однако, учитывать занижение результатов АСМ-измерения этого параметра в силу контактных деформаций, см. раздел 2.1). Типичные АСМ-изображения подобных искусственно сформированных дефектов в белковых пленках изображены на рис. 5.10.
Пары изображений а) и б), в) и г), д) и е), з) и и) рис. 5.10 получены в режиме модуляции силы (см. раздел 1.3.2), первое изображение пары является обычным топографическим, а второе отображает величину амплитуды колебаний левера, находящегося в контакте с поверхностью и испытывающего со стороны пьезоманипулятора периодическое воздействие. Согласно изложенному в разделе 1.3.2, более светлые участки этого амплитудного изображения соответствуют менее жестким участкам исследуемой поверхности. Белковые пленки были выбраны в качестве объекта этих измерений именно в силу того, что они достоверно характеризуется меньшим значением модуля упругости в сравнении с подложкой.
Метод формирования искусственных дефектов в тонких пленках может найти универсальное применение. Формирование искусственных дефектов позволяет не только детектировать присутствие пленки в случае ее высокого качества, но и определять толщину покрытия, а также, например, при исследовании молекулярной упаковки, анализировать взаимоориентацию решеток пленки и подложки. Более того, путем определения предельного силового воздействия зонда, необходимого для того, чтобы искусственный дефект был сформирован, можно проводить сравнительный анализ механической прочности тонкопленочных покрытий.
Основным преимуществом метода сканирующей зондовой микроскопии при исследовании молекулярной структуры тонких пленок является возможность определения параметров элементарной ячейки для локально-неупорядоченных, островковых и доменных пленок. Действительно, в настоящее время для исследования структуры поверхности широко используются дифракционные методы: дифракция электронов, рентгеновских лучей, нейтронов, ионов и нейтральных атомов [116,117]. Эти методы позволяют определять структуру поверхности с высокой степенью точности лишь при условии, что исследуемая поверхность характеризуется определенной степенью упорядоченности во всей области взаимодействия с зондирующим пучком. Наличие в пленке разориентированных доменов, островков, кристаллических включений и пр. осложняет дифракционную картину и создает сложности при интерпретации результатов и определении искомых параметров.
В этой связи метод атомно-силовой микроскопии представляется весьма перспективным для исследования структуры молекулярной упаковки тонких органических пленок, в том числе в сравнении и с микроскопическими методами. Действительно метод АСМ характеризуется невысокой интенсивностью воздействия на образец (в отличие от электронной и ионной микроскопии), позволяет проводить исследования диэлектрических структур (в отличие от СТМ) и характеризуется высоким пространственным разрешением, поскольку позволяет определять параметры молекулярной или атомной упаковки поверхности.
Влияние принципа плотной упаковки углеводородных цепей на параметры ячейки н-парафинов было рассмотрено Китайгородским [118,119] на основе геометрического анализа. (Нормальные парафины (CnH2n+2) являются простейшими полимерами, и выяснение их строения имеет большое значение для понимания общих закономерностей органической кристаллохимии). Исходным положением служило рассмотрение плотного контакта молекул плоской зигзагообразной формы, оси которых параллельны. Данный анализ позволил определить возможные упаковки молекул в н-парафинах, см. таблицу 5.1. Известные кристаллические структуры н-парафинов относятся к небольшому числу предсказанных в результате рассматриваемого анализа. Как оказалось, принцип плотной упаковки углеводородных цепей может определять и молекулярную упаковку амфифильных молекул, имеющих полярную группу и длинный углеводородный фрагмент (молекул жирных кислот и т.п.).
Слой | Параметры ячейки слоя, нм | ja | jb | jc* |
H [0, 0] | a=0,48, g = 120° | 90° | 90° | 0° |
R [0, 0] | a=0,742, b=0,496, g = 90° | 90° | 90° | 0° |
R [0, ±2] | a=0,900, b=0,496, g = 90° | 90° | 55°, 124° | 34° |
R [±1, 0] | a=0,742, b=0,557, g = 90° | 63°, 117° | 90° | 27° |
R [0, ±1] | a=0,785, b=0,496, g = 90° | 90° | 71°, 109° | 19° |
R [±1, ±1] | a=0,785, b=0,557, g = 81° | 63°, 117° | 71°, 109° | 31° |
R [-+1, ±1] | a=0,785, b=0,557, g = 98° | 63°, 117° | 71°, 109° | 31° |
M [0, 0] | a=0,42, b=0,44, g = 111° | 90° | 90° | 0° |
M [±1, 0] | a=0,49, b=0,44, g = 107° | 59°, 121° | 90° | 32° |
M [0, ±1] | a=0,42, b=0,51, g = 107° | 90° | 60°, 120° | 32° |
T [±1/2, 0] | a=0,43, b=0,45, g = 103° | 73°, 107° | 90° | 18° |
T [±1/2, 1] | a=0,43, b=0,52, g = 109° | 73°, 107° | 120°, 60° | 36° |
В работе [120] авторы проводили сравнительный анализ применения методов АСМ и дифракции электронов к исследованию молекулярной упаковки твердых растворов н-парафинов и многокомпонентных восков. Методом атомно-силовой микроскопии (результаты находились в согласии с данными дифракции электронов) были определены параметры элементарной ячейки (прямоугольная) как для поверхности кристалла чистого н-C36H74 (a = 0,47(5) нм, b = 0,78(5) нм), так и поверхностей твердых растворов: н-C35H72/ н-C36H74 (a = 0,46(6) нм, b = 0,76(1) нм) и C32H66/ н-C36H74 (a = 0,59(6) нм, b = 0,78(5) нм). Возможность наблюдения с помощью АСМ двумерной упорядоченной структуры в последнем случае достаточно неожиданна, поскольку строение молекул компонентов твердого раствора различается на четыре метиленовых группы. Можно предположить, что визуализация молекулярной упаковки была возможна на участках локальной однородности пленок.
Работа [121] посвящена АСМ-исследованию пленок транс-6-октадеценовой кислоты, перенесенных по методике ЛБ (25 или 5 слоев) на поверхность стекла с субфазы двух типов: CdCl2 или TbCl3 (концентрацией 10-4 М). Авторы обнаружили, что тип упаковки пленки существенно зависит от природы противоина. Так, пленки, перенесенные на подложку с субфазы, содержащей ионы Tb 3+, обнаруживали в фурье-спектре несколько кратных рефлексов, причем их значения позволяют предположить, что упаковка пленки характеризуется переорганизованной в сверхструктуру центрированной прямоугольной (ЦП) решеткой с параметрами a = 0,46±0,02 нм и b = 0,89±0,03 нм. Площадь, приходящаяся на молекулу, в этом случае составляет 0,205±0,010 нм2, что хорошо согласуется со значением предельной молекулярной площади 0,210±0,004 нм2, которое было определено авторами по изотерме сжатия пленки на поверхности субфазы. В то же время пленки, перенесенные с CdCl2 субфазы характеризовались косоугольной решеткой с параметрами: a1 = 0,48±0,03 нм, a2 = 0,51±0,03 нм и g = 68±3°, что соответствует площади на молекулу 0,23±0,02 нм2.
Авторы работы [122] исследовали методом АСМ молекулярную упаковку гидрофобной поверхности верхнего слоя 3-, 5- и 7-слойного покрытия арахидата кадмия (CdA2). Авторы определили параметры ячейки для данной поверхности (тип: прямоугольная с двумя молекулами в базисе), которые составили: a = 0,48±0,03 нм, b = 0,74±0,04 нм. Авторы отмечают, что в ряде измерений параметры ячейки были искажены и связывают это с влиянием зонда микроскопа. Однако представляется, что локальные деформации или локальные смещения молекул под воздействием зонда в процессе сканирования не могут сказаться на результатах определения параметров элементарной ячейки при измерениях АСМ, поскольку воздействие зонда не может приводить к ``переупаковке'' молекул на всем сканируемом участке. По всей видимости, наблюдаемые искажения могут быть объяснены влиянием температурного дрейфа. Стоит отметить, что обнаруженные авторами параметры решетки близки параметрам плотной упаковки слоев углеводородных цепей R[0,0], см. таблицу 5.1.
Сосуществование трех различных типов молекулярных упаковок в трехслойных ЛБ пленках арахидата бария (BaA2) было обнаружено авторами работы [123]. Это демонстрирует высокую чувствительность организации пленок к структуре противоиона (который определяет площадь полярной группы на поверхности субфазы): замена иона Cd 2+ ионом Ba 2+ существенно изменяет параметры молекулярной упаковки. Три обнаруженных типа упаковки для исследуемых пленок характеризовались близкими значениями площади, приходящейся на молекулу. При этом каждый тип упаковки различался ориентацией молекул относительно направления, нормального к поверхности подложки (что определялось по различающейся высоте соответствующих участков АСМ-изображения пленки над поверхностью подложки).
Параметры первого из обнаруженных типов элементарной ячейки (триклинная ячейка, a1 = 0,42±0,01 нм, a2 = 0,45±0,01 нм, g = 76±1°) позволили авторам сделать вывод, что в этом случае структура пленки определяется принципом плотной упаковки углеводородных цепей с образованием слоев, характеризующихся ``Т''-подъячейкой (см. таблицу 5.1). Второй тип упаковки идентифицировался авторами как структура с ``R''-подъячейкой (R[0, ±1], см. таблицу 5.1), также обусловленной принципом плотной упаковки углеводородных цепей (ЦП ячейка, a = 0,51±0,01 нм, b = 0,79±0,01 нм). Третья обнаруженная структура упаковки была менее упорядочена в сравнении с первыми двумя, что не позволило определить параметры ячейки. Угол наклона молекул (измерялся по толщине слоев) составлял 26° и 19° соответственно для первых двух типов, третий тип упаковки был образован, по-видимому, вертикально ориентированными молекулами. Весьма интересно, что все типы упаковки характеризовались близкими значениями площади, приходящейся на молекулу, причем при исследованиях конкретных пленок в различных участках их поверхности было возможно обнаружить все три типа упаковки (авторы определили, что первый тип упаковки занимает 70% поверхности пленки, второй и третий -5 и 25% соответственно).
Как показано в работе [124], весьма специфическое влияние на структуру пленки арахидатов оказывает использование в качестве противоиона Ca 2+. Авторы исследовали трехслойную ЛБ пленку арахидата кальция (CaA2) и обнаружили, что структура молекулярной упаковки пленки не обладает зеркальной симметрией. Пленка имеет доменную структуру с двумя типами решетки, которые являются зеркальным отражением друг друга. Элементарная ячейка включает восемь молекул, с углом наклона их к подложке в 26°, параметры ячейки: a1 = 1,86±0,01 нм, a2 = 0,894±0,006 нм, g = 71,3±0,5°). Авторы связывают наблюдаемую специфику структуры пленки с конкуренцией требования плотной упаковки углеводородных цепей и значением площади полярной группы, определяемым радиусом ее экранировки противоионом кальция.
Обобщающее рассмотрение влияния природы противоина, присутствующего в субфазе, на структуру пленок солей жирных кислот приводится в работе [125]. Авторы указывают на корреляцию между значениями площади, приходящейся на молекулу слоя, и электроотрицательностью противоина, см. таблицу 5.2. По мере уменьшения электроотрицательности противоина возрастает степень ионности связи и, как оказалось, возрастает площадь, приходящаяся на молекулу в пленке. Так, Cd, Mn и Pb образуют ковалентную связь, а щелочноземельные (Ba, Ca и Mg) взаимодействуют с ионами жирных кислот электростатически. Очевидно, именно тип связи определяет площадь, приходящуюся на отдельную полярную группу.
Противоион | тип элементарной ячейки | наклон молекул | молекулярная площадь, нм2 | электроотрицательность |
Pb | ЦП | 0° | 0,179 | 1,8 |
Cd | ЦП | 0° | 0,18 | 1,7 |
Mg | ЦП | 0° | 0,187 | 1,3 |
Mg | триклинная | 0° | 0,187 | 1,3 |
Mn | ЦП | 19° | 0,195 | 1,5 |
Ca | триклинная | 26° | 0,196 | 1,0 |
Ba | ЦП | 19° | 0,202 | 0,9 |
Ba | триклинная | 26° | 0,204 | 0,9 |
Zn | гексагональная | 33° | 0,222 | 1,6 |
При этом требование плотной упаковки углеводородных частей молекул также влияет на структуру пленок. Если площадь, приходящаяся на полярные группы молекул больше, чем площадь поперечного сечения плотноупакованных углеводородных цепей (около 0,18 нм2), то удовлетворить требованию плотной упаковки ``хвостов'' можно лишь допустив возможность наклона осей молекул к подложке. Как показали авторы, значение угла наклона (и параметров ячейки), как правило, не произвольно, но подчиняется закономерностям, предсказанным Китайгородским [118] для структуры упаковки н-парафинов (см. таблицу 5.1).
В то же время, как отмечают авторы рассматриваемой работы, неожиданными были результаты исследования молекулярной упаковки пленки (7 слоев) арахидата цинка. Они обнаружили, что структура пленки описывается гексагональной решеткой с параметром 0,47-0,48 нм. Высокая симметрия упаковки пленки соответствует, т.н. ``газокристаллическому'' состоянию (характеризующемуся невысокой плотностью упаковки), когда молекулы сохраняют некоторую ротационную свободу (это обеспечивает аксиальную симметрию области, занимаемой каждой молекулой) [119].
В работе рассматривается также эффект влияния свойств подложки на структуру молекулярной упаковки сформированных на ее поверхности пленок слой жирных кислот. Это влияние может проявляться весьма различным образом для разных противоинов. Отметим, в частности, факт уменьшения молекулярной площади в пленках стеарата свинца (с 0,21 до 0,18 нм2), сформированных на поверхности слюды, при увеличении числа слоев в пленке (1-3-5-7). Молекулярная упаковка верхнего слоя 7-слойного покрытия релаксирует, по-видимому, к невозмущенному влиянием подложки состоянию, поскольку те же параметры упаковки характеризуют 200-слойное покрытие (по данным рентгеноструктурного анализа). При использовании в качестве подложки поверхности кремния (покрытой оксидом и существенно более инертной, чем слюда) молекулярная упаковка верхнего слоя уже трехслойного покрытия характеризуется теми же невозмущенными параметрами.
Авторы работы рассматривают также механизмы формирования сверхструктур в молекулярной упаковке поверхности, связывая их возникновение с присутствием регулярных дефектов упаковки, обусловленных необходимостью удовлетворения требованиям плотной упаковки углеводородных частей и полярных групп молекул.
Сходное обобщение АСМ-результатов исследования упаковки молекул жирных кислот в зависимости от условий формирования пленок (влияние субфазы и природы подложки) проводится в обзоре [126]. Подчеркивается, что параметры элементарной ячейки, симметрия и площадь на молекулу не зависят от длины алкильной цепи жирных кислот (С16-С22), а контролируются особенностями взаимодействия противоиона с карбоксильной группой (СООН). Отмечен факт [126] невозможности достигнуть молекулярного разрешения на монослойных пленках солей жирных кислот (за исключением стеарата свинца), в отличие от мультислойных покрытий.
Комплексное исследование пленок жирных кислот, перенесенных на поверхности слюды, пирографита или кремния, осуществили авторы работы [127]. Применение метода АСМ позволило авторам определить параметры элементарной ячейки для кадмиевых солей стеариновой, арахиновой, бегеновой и лигноцериновой кислот, см. табл. 5.3. Как отмечают авторы, молекулярное разрешение могло быть достигнуто лишь в случае, когда толщина покрытия превышала три монослоя (даже для пленок, осажденных горизонтальным методом). Значения параметров элементарной ячейки, приведенные в двух нижних строках табл. 5.3, могут быть обусловлены принципом плотной упаковки углеводородных цепей: ср. со значениями параметров a = 0,496 нм и b = 0,785 нм упаковки слоев углеводородных цепей R[0, ±1] таблицы 5.1.
a, нм | b, нм | |
стеарат (C18) | 0,45-0,47 | 0,74-0,76 |
арахидат (C20) | 0,46-0,48 | 0,73-0,76 |
бегенат (C22) | 0,47-0,49 | 0,73-0,78 |
лигноцеринат (C24) | 0,46-0,49 | 0,79-0,80 |
Авторы работы [128] сообщают о том, что ими впервые методом АСМ была визуализована молекулярная упаковка монослойных пленок жирных кислот (лигноцериновой и стеариновой), перенесенных на слюду с водной поверхности. В схеме эксперимента использовалось последовательное многоступенчатое сжатие монослоя на поверхности раздела фаз, что, по-видимому, и позволило перенести качественный монослой на поверхность подложки по традиционной вертикальной методике ЛБ. Авторы определяют молекулярную упаковку поверхностей монослоев как гексагональную с параметром ячейки a = 0,43 нм для лигноцериновой и a = 0,42 нм для стеариновой кислот. Однако эти значения соответствуют 0,15 и 0,16 нм2 площади, приходящейся на молекулу в слое, что, по-видимому, маловероятно, поскольку минимальное значение данного параметра для кристаллов рассматриваемых жирных кислот составляет около 0,18 нм2, см. [118]. Наблюдаемое расхождение объясняется, по-видимому, ошибкой в калибровке пьезоманипулятора, допущенной авторами: при калибровке использовалось значение 0,46 нм для параметра ячейки слюды, в то время как этот параметр составляет 0,52 нм. Если воспользоваться данным значением и пересчитать результаты измерений авторов, то для параметров ячейки лигноцериновой и стеариновой кислот можно получить a = 0,49 нм и a = 0,47 нм соответственно, что, в свете результатов рассмотренных выше работ, представляется вполне разумными значениями (и соответствует параметрам ``газокристаллического'' состояния).
Метод АСМ позволяет исследовать молекулярную упаковку некоторого выбранного локального участка поверхности заданной площади. Эту площадь участка в экспериментах выбирают исходя из следующих двух требований.
Основным фактором, ограничивающим точность метода АСМ при исследованиях молекулярной упаковки поверхности твердых тел, является температурный дрейф зонда микроскопа относительно образца. Если скорость температурного дрейфа vдрейф в направлении оси сканирования сравнима со скоростью перемещения левера в этом направлении, то возникают значительные искажения при измерении расстояния между двумя точками поверхности. В зондовой микроскопии различают направления быстрого сканирования (направление сканов) и медленного сканирования (перпендикулярное направление, в котором перемещается левер при переходе от скана к скану). В направлении оси медленного сканирования левер перемещается в среднем со скоростью vскан = Lf/N, а в направлении оси быстрого сканирования -со скоростью vскан = 2Lf; здесь L -размер АСМ-изображения (считаем, что Lx = Ly), f -частота строчной развертки сканирования, N -число строк в АСМ-изображении. Пусть температурный дрейф направлен вдоль оси медленного сканирования и имеет скорость vдрейф, и пусть необходимо измерить отрезок длины a, параллельный оси медленного сканирования, тогда ошибка измерения, обусловленная температурным дрейфом, может быть определена следующим образом:
| (5.1) |
| (5.2) |
Было обнаружено, что в стандартных условиях эксперимента (без применения термостатирования) искажающее влияние температурного дрейфа может проявляться и при максимально возможных параметрах f и L (fmax = 60 Гц, Lmax ~ 70ё80 нм), что позволяет оценить по формуле (5.1) величину скорости теплового дрейфа: vдрейф <~0,8 нм/с (N=512, e = 0,1, f=60 Гц, L=70 нм).
Оказывается, что температурный дрейф характеризуется постоянным направлением. Поскольку он особенно чувствителен в том случае, когда компонента vдрейф, направленная вдоль оси медленного сканирования, принимает значительную величину, то полезно при сканировании записывать АСМ-изображения парами -при двух взаимнопротивоположных направлениях медленного сканирования. В таком случае, при усреднении измеренных параметров элементарной ячейки для первого и второго изображения имеет место взаимокомпенсация ошибки, связанной с влиянием температурного дрейфа в направлении оси медленного сканирования. Для уменьшения влияния температурного дрейфа в направлении перпендикулярной оси быстрого сканирования следует определять параметры молекулярной упаковки путем усреднения значительного числа измерений АСМ-изображений, полученных при разных направления сканирования.
Усреднение параметров элементарной ячейки на основании нескольких АСМ-изображений с различной ориентацией осей сканирования относительно образца позволяет свести вносимую дрейфом погрешность в измерениях к ошибке среднего арифметического, и достигнуть точности в определении искомых параметров на уровне нескольких процентов, см. ниже.
Образцы пленок были приготовлены вертикальным методом (Y-пленки) и методом горизонтального осаждения (Z-пленки) при комнатной температуре (t ~ 25°С) и поверхностном давлении p в диапазоне 20-40 мН/м. Монослой, сформированный на поверхности субфазы, осаждали методом ГО на предварительно погруженную под водную поверхность подложку в процессе медленного (около 0,01 мм в минуту) понижения уровня жидкости (рис. 5.1б). Пленки, полученные вертикальным методом ЛБ, формировали при пропускании подложки сквозь монослой вертикально со скоростью 4 мм/мин. Пленкообразующее вещество наносили на поверхность субфазы из раствора в хлороформе (х. ч.). При приготовлении пленок пчелиного воска использовали препарат, очищенный разделением на фракции.
Калибровка прибора. Все эксперименты по определению молекулярной упаковки проводили на приборе ``Nanoscope-IIIa'' (Digital Instruments, США), оборудованном ``D''-сканером (максимальный размер кадра 15 мкм), в контактном режиме с использованием кантилевера жесткостью 0,32 Н/м, см. таблицу A.1. Для калибровки сканера использовали поверхности тест-объектов: слюды и пирографита, для которых измеряли параметры решетки. Полученные значения усредняли по результатам обработки около 100 АСМ-изображений, сравнивали с известными параметрами гексагональной упаковки используемых тест-объектов (0,519 нм для слюды и 0,246 нм для графита), и определяли коэффициент калибровки. Полученные значения для коэффициента калибровки совпали для двух типов используемых тест-объектов с точностью до 2%. Измеренные значения характеризовались среднеквадратичной ошибкой меньше 2%. Таким образом, можно оценить, что точность достигнутой калибровки не хуже, чем 2%.
Основным препятствием при АСМ-исследованиях молекулярной упаковки тонких пленок является локальное разрушение покрытия под воздействием зонда микроскопа при уменьшении размера кадра до 70-80 нм (уменьшение размера кадра приводит к уменьшению скорости перемещения зонда и, следовательно, увеличению характерного времени взаимодействия зонда и образца, см. раздел 2.1.2).
Тем не менее, для ряда тонкопленочных покрытий (устойчивых к воздействию зонда) оказалось возможным достичь с помощью АСМ молекулярного разрешения и определить параметры элементарной ячейки. Примеры АСМ-визуализации молекулярных упаковок тонких пленок приведены на рис. 5.11, рис. 5.13 и рис. 5.15.
| ,
Figure |
В таблице 5.4 приведены значения средних арифметических для измеряемых параметров, полученные путем усреднения результатов обработки нескольких десятков АСМ-изображений. В качестве погрешности измерения приведена величина стандартного отклонения среднего арифметического.
No, тип | подложка | a, нм | b, нм | b/a | S, нм2 | N |
1 | слюда | 0,48±0,03 | 0,77±0,06 | 1,59±0,16 | 0,185±0,017 | 60 |
ЦП | e = 6% | e = 8% | e = 10% | e = 9% | ||
2 | графит | 0,50±0,02 | 0,78±0,06 | 1,56±0,14 | 0,194±0,016 | 70 |
ЦП | e = 4% | e = 8% | e = 9% | e = 8% | ||
3 | слюда | 0,47±0,02 | 0,88±0,08 | 1,9±0,2 | 0,210±0,018 | 70 |
ЦП | e = 4% | e = 9% | e = 11% | e = 9% | ||
4 | слюда | 0,53±0,03 | 0,77±0,04 | 1,46±0,13 | 0,203±0,007 | 10 |
ЦП | e = 5% | e = 5% | e = 9% | e = 3% | ||
5 | слюда | 0,477±0,014 | 0,197±0,007 | 46 | ||
Г | e = 3% | e = 4% | ||||
6 | слюда | 0,476±0,018 | 0,196±0,008 | 30 | ||
Г | e = 4% | e = 4% | ||||
7 | слюда | 0,477±0,011 | 0,197±0,006 | 50 | ||
Г* | e = 2% | e = 3% | ||||
8 | слюда | 0,566±0,015 | 0,80±0,03 | 1,42±0,08 | 0,227±0,011 | 50 |
ЦП | e = 3% | e = 4% | e = 6% | e = 5% | ||
9 | слюда | 0,482±0,010 | 0,202±0,005 | 15 | ||
Г | e = 2% | e = 5% | ||||
10 | слюда | 0,548±0,009 | 1,56±0,08 | 2,85±0,09 | 0,427±0,014 | 20 |
ЦП* | e = 2% | e = 5% | e = 3% | e = 3% | ||
11а | слюда | 0,548±0,009 | 1,31±0,06 | 2,40±0,11 | 0,360±0,018 | 67 |
ЦП | e = 2% | e = 5% | e = 5% | e = 5% | ||
11б | слюда | 0,550±0,012 | 1,30±0,07 | 2,37±0,014 | 0,36±0,02 | 45 |
ЦП | e = 2% | e = 6% | e = 6% | e = 5% | ||
12 | слюда | 0,508±0,007 | 1,37±0,08 | 2,71±0,017 | 0,35±0,02 | 40 |
ЦП | e = 1% | e = 5% | e = 6% | e = 5% | ||
13 | графит | 0,50±0,04 | 0,79±0,07 | 1,6±0,2 | 0,195±0,015 | 15 |
ЦП | e = 7% | e = 9% | e = 13% | e = 8% |
No | материал пленки | давление p (мН/м) | подложка и способ нанесения | субфаза | сформированное покрытие |
1 | Бегеновая кислота | 30 | слюда, ГО | CdTe | ламели |
2 | 30 | графит, ЛБ | CdTe | ламели | |
3 | 30 | слюда, ЛБ | CdTe | монослой | |
4 | 30 | слюда, ГО | H2O | монослой | |
5 | N,N'-диоктадецилпропандиамид | 30 | слюда, ГО | H2O | трислой |
6 | N-гексадецилацетамид | 30 | слюда, ГО | H2O | трислой |
7 | N,N'-дигексадецилпропандиамид | 25-30 | слюда, ГО | H2O | трислой |
8 | N-гексадецил-N'-(2-нафтил)пропандиамид | 30 | слюда, ГО | H2O | монослой |
9 | N-гексадецил-N'-(6-хинолин)пропандиамид | 32-34 | слюда, ГО | H2O | монослой |
10 | 5-октадецил-2,4,6 | 30 | слюда, ГО | H2O | монослой |
11а | (1Н, 3Н, 5Н) пири- | 30 | слюда, ГО | H2O | ламели |
11б | мидинтрион | 35 | слюда, ГО | H2O | ламели |
12 | 35 | слюда, ГО | H2O | островки | |
13 | Пчелиный воск | 20 | графит, ГО | H2O | монослой |
Бегеновая кислота. Структурная формула: С21H43COOH.
Молекулярное разрешение удавалось стабильно получать на поверхностях кристаллитов (3 слоя), спонтанно сформированных в монослоях бегеновой кислоты, переносимой на поверхность подложки (слюда и пирографит) с субфазы, содержащей CdTe-кластеры, см. раздел 5.1.
Было установлено, что молекулярная упаковка верхнего слоя кристаллической структуры на слюде, может быть описана как прямоугольная центрированная (No 1 таблицы 5.4). При этом направление вектора a совпадает с направлением ламелей кристаллических структур. Т.е. для направления вектора a справедливы все результаты, изложенные выше относительно направления ламелей (об углах с кристаллическими осями слюды и о взаимной ориентации ламелей, относящихся к различным кристаллическим структурам, см. гистограммы рис. 5.5а и рис. 5.6а).
Молекулярная упаковка верхнего слоя кристаллической структуры, образованной на поверхности пирографита, также описывается ЦП-решеткой с близкими значениями параметров ячейки (No 2 таблицы 5.4). При этом направление вектора a также совпадает с направлением ламели (об ориентации ламелей на пирографите см. гистограммы рис. 5.5б и рис. 5.6б).
В обоих случаях параметры решетки близки параметрам для R[0,±1]-слоев плотноупакованных молекул н-парафинов (см. таблицу 5.1). Было обнаружено, что те же параметры характеризуют упаковку молекул монослоя пчелиного воска, перенесенного на поверхность пирографита (No 13 таблицы 5.4). Близкие параметры зафиксировали методом АСМ авторы работы [120] для поверхностей кристаллов н-парафинов, а также авторы работы [127] для тонких пленок (более 3-х слоев) бегената и лигноцерината кадмия (см. таблицу 5.3). Это позволяет предположить, что молекулярная упаковка поверхности кристаллитов, сформированных в монослоях бегеновой кислоты на CdTe-субфазе, в незначительной степени подвержена влиянию подложки (действительно, параметры решетки кристаллитов близки в обоих случаях: пирографитовой подложки и подложки из слюды), но определяется принципом плотной упаковки углеводородных цепей.
По-видимому, в большей степени влияние подложки проявляется при формировании структуры монослойных пленок бегеновой кислоты, перенесенных на поверхность слюды (No 3 и No 4 таблицы 5.4). В отличие от авторов работы [128], также исследовавших монослои жирных кислот, перенесенных на подложку с водной субфазы, мы обнаружили не гексагональную (``газокристаллическая фаза''), а типичную ЦП-решетку молекулярной упаковки монослоев в этих случаях.
Значение параметра b = 0,77±0,04 нм молекулярной упаковки монослоя бегеновой кислоты, сформированного методом ГО (No 4), близко значению b = 0,785 нм, определяемому принципом плотной упаковки углеводородных частей молекул для слоев с R-подъячекой, см. таблицу 5.1. При этом значение параметра a = 0,53±0,03 нм может объясняться влиянием подложки.
С другой стороны, параметры решетки (a = 0,47±0,02 нм и b = 0,88±0,08 нм) монослоя, перенесенного на слюду вертикальным методом ЛБ (No 3 таблицы 5.4) близки параметрам, обнаруженным авторами работы [121] для верхнего слоя мультислойных пленок транс-6-октадеценовой кислоты, перенесенных на стекло с TbCl3-субфазы, когда влияние подложки вряд ли может иметь место.
Пленки амфифильных кетоамидов. В таблице 5.4 приведены результаты исследования пленок, сформированных на основе: N,N'-диоктадецилпропандиамида (No 5) C18H37NH-CO-CH2-CO-NHC18H37, N-гексадецилацетамида (No 6) C16H33NH-CO-CH2-CO-CH3, N,N'-дигексадецилпропандиамида (No 7) C16H33NH-CO-CH2-CO-NHC16H33, N-гексадецил-N'-(2-нафтил)пропандиамида (No 8), N-гексадецил-N'-(6-хинолин)пропандиамида(No 9)
Изображения структуры исследуемых пленок приведены на рис. 5.12. Соединения No 8 и No 9 позволяют формировать на своей основе качественные монослои с незначительным количеством дефектов, см. рис. 5.12г и д. В то же время из рисунка видно, что для образцов No 5-7, имеет место структурная перестройка молекул пленки с формированием ими образований, высота которых превышает высоту монослоя и составляет величину 5,5-6,0 нм (см. рис. 5.12а-в). Можно предположить, что эти структуры образованы тремя монослоями молекул (высоты тонких пленок по результатам АСМ-исследований занижены в силу контактных деформаций, см. раздел 2.1). С учетом того обстоятельства, что молекулы образцов No 5 и No 7 имеют по две углеводородных цепи равной длины можно предположить одинаковую ориентацию этих цепей относительно полярной группы молекулы. В пользу этого заключения говорит тот факт, что структура пленки образца No 6 сходна с образцами No5 и No 7, хотя в первом случае молекула имеет лишь один длинный углеводородный фрагмент.
Было обнаружено (рис. 5.13), что за одним исключением (No 8) АСМ-изображения пленок исследованных соединений характеризуются гексагональной упаковкой, с параметром ячейки близким значению 0,48 нм, которое характеризует упаковку слоев H[0, 0] (таблица 5.1) углеводородных цепей с гексагональной подъячейкой. При этом упаковка АСМ-изображения образца No 7 характеризуется сверхструктурой: в фурье-спектре наблюдается максимум на пространственной частоте 0,826 нм (a×Ц3) вдоль направления [21], см. рис. 5.13в. В прямом пространстве это проявляется в виде чередования яркости АСМ-изображения рядов (через один ряд), перпендикулярных направлению [21].
Можно предположить, что в этих случаях АСМ-изображение двумерной периодической структуры является отображением упаковки углеводородных цепей. Молекулы образцов No 5 и No 7 имеют по два углеводородных фрагмента равной длины. Поэтому можно предположить, что каждой молекуле этих соединений на АСМ-изображении структур исследуемых пленок соответствуют две пучности. Т.о. мы имеем дело с гексагональной упаковкой углеводородных цепей, причем именно тот факт, что эти цепи попарно объединены в молекулы, возможно, и проявляется в наблюдаемой сверхструктуре (для образца No 7).
Факт упаковки углеводородных цепей в гексагональную решетку для молекул, имеющих по два углеводородных ``хвоста'' является несколько неожиданным. Действительно, согласно имеющимся представлениям [119] гексагональная упаковка углеводородных цепей в слой H[0, 0] не характеризуется высокой плотностью: молекулы в слое ориентированы вертикально, но в то же время молекулярная площадь составляет около 0,2 нм2, тогда как более плотные упаковки вертикально ориентированных молекул могут характеризоваться значением площади 0,18 нм2. Считается, что молекулы в слое H[0, 0] сохраняют некоторую свободу вращения вокруг вертикальной оси, что приводит к аксиальной симметрии области, занимаемой молекулой, и обуславливает высокую симметрию молекулярной упаковки. Это состояние органического вещества называют [119] ротационно-кристаллическим или ``газокристаллическим'' состоянием. Поэтому для объяснения высокой симметрии упаковки в нашем случае следует допустить, что два углеводородных хвоста молекул исследуемых структур сохраняют возможность независимого вращения вокруг своей оси.
Отдельно стоит отметить, что молекулярная упаковка пленки на основе соединения No 9 характеризуется относительной неоднородностью и направления векторов трансляции различны в пределах различных участков пленки, причем размеры участков фиксированной ориентации составляют лишь сотни квадратных нанометров. При этом четко визуализировать границы доменов не удавалось, что, по-видимому, обусловлено изложенными в разделе 2.1.5 выводами относительно ``ложности'' молекулярного разрешения в АСМ при типичных параметрах эксперимента на воздухе.
Сосуществование различных типов упаковок в структуре пленки 5-октадецил-2,4,6(1Н, 3Н, 5Н)пиримидинтриона. Для этого соединения были обнаружены три различных типа упаковки, см. таблицу 5.4: для монослоя (No 10, таблица 5.4), для кристаллитов (No 11) и для островковых образований (No 12).
Структура пленки визуализирована методом АСМ на рис. 5.14, где изображены области различной структурной организации сформированных покрытий: монослой, островки и ламели. Анализ высот показывает, что островки сформированы, по-видимому, двумя слоями пленки (высота этих образований около 4 нм, высота монослоя 1,5-2 нм), в то время как кристаллиты (ламели) образованы 3, 6 или 9 слоями (высоты около 5,3-5,7 нм, а также удвоенные и утроенные значения.)
Молекулярная упаковка островковых и кристаллических образований характеризовалась центрированной прямоугольной ячейкой с отличающимися параметрами a и b, но близкими значениями площади, приходящейся на молекулу, см. таблицу 5.4.
В то же время АСМ-изображение молекулярной упаковки монослоя характеризовалось, во-первых, сверхструктурой (рис. 5.15а), а во-вторых, большей площадью на молекулу. Последнее может объясняться влиянием подложки. Действительно, было обнаружено, что направление вектора b для монослойного покрытия в этом случае совпадает с одной из кристаллографических осей слюды. Периодичность упаковки пленки вдоль этого направления характеризовалась значением 1,56 нм, в то время как значение этого параметра для подложки 0,519 нм, что в три раза меньше. Таким образом имеет место формирование монослоем решетки, один из параметров которой соизмерим с параметром решетки подложки (второй параметр пленки имеет то же значение, что и для кристаллических ламелярных структур).
Успешная визуализация молекулярной структуры пленки возможна лишь при условии достаточно прочной связи молекул пленки друг с другом и с подложкой (что позволяет проводить неразрушающие исследования). Эта связь может возникнуть в результате следующих причин. Во-первых, за счет гидрофобного взаимодействия алкильных участков молекул и гидрофобной подложки. Это позволило наблюдать молекулярную упаковку пленок пчелиного воска (No 13 таблицы 5.4). Во-вторых, за счет процессов спонтанной кристаллизации или самоорганизации молекул пленки, приводящих к формированию плотноупакованных структур, характеризующихся высокой связанностью молекул (No 1, No 2, No 5-7, No 11 и No 12 таблицы 5.4).
В отсутствие этих факторов достижение молекулярного разрешения методом АСМ остается сложной задачей: структура пленки разрушается под воздействием зонда микроскопа при уменьшении размера кадра (при уменьшении размера кадра уменьшается скорость сканирования и увеличивается характерное взаимодействие зонда с локальным участком образца) несмотря на минимизацию силового взаимодействия зонда и образца. Возможно, решить эту проблему удастся при проведении исследований в жидких средах.
Тем не менее, в тех случаях, когда пленка характеризовалась достаточной прочностью, что позволяло проводить неразрушающие исследования малого локального участка поверхности на воздухе, оказалось возможным определить молекулярную структуру ряда монослойных тонкопленочных покрытий (No 3, No 4, No 8-10).
Было обнаружено, что для многих пленок параметры ячейки подчиняются принципу плотной упаковки углеводородных цепей молекул (ср. с таблицей 5.1). Так, образцы No 1, No 2, No 13 характеризуются ЦП-решеткой с параметрами, близкими параметрам слоя R[0, ±1] плотноупакованных молекул н-парафинов (таблица 5.1). В свою очередь молекулярные пленки образцов No 5-7 и No 9 характеризуются гексагональной упаковкой и параметр ячейки близок значению для слоя H[0,0] (таблица 5.1) ``газокристаллической'' фазы н-парафинов.
С другой стороны молекулярная структура пленок No 10-12 не может быть обусловлена принципом плотной упаковки углеводородных цепей. Площадь, приходящаяся на молекулу, в этом случае (0,35-0,36 нм2) значительно превышает площадь поперечного сечения углеводородной цепи (около 0,18 нм2). По-видимому, в этом случае параметры решетки определяются значением площади, занимаемой полярной группой на поверхности субфазы.
Наряду с этим отличие, наблюдаемое для ряда пленок, параметров упаковки от значений, реализующих плотную упаковку углеводородных групп, может быть обусловлено и влиянием подложки. В частности, обнаружено формирование молекулами монослоя решетки с параметрами, соизмеримыми с параметрами подложки (No 10).
Т.о. параметры молекулярной упаковки тонких пленок определяются несколькими конкурирующими факторами: принципом плотной упаковки углеводородных цепей, значением площади полярной группы и влиянием подложки.
Ниже приведен текст программы, написанной на языке C++.
#include <math.h>
#include <stdlib.h>
#include <io.h>
#include <iostream.h>
double F,E1,E2,D,A,B,a,b,a0,b0,R1,R2,h,b_a_2,k,a_b,P,C;
double K_ell(double k_1) //first elleptical integral
{
double K_count = 1.0, q_count = 1.0, m;
m = (1.0 - k_1)/(1.0 + k_1);
for (int i = 1; i < 1000; i += 2)
{
q_count = q_count*i*i/(i + 1)/(i + 1)*m*m;
K_count = K_count + q_count;
};
return (M_PI/2.0*(1.0 + m)*K_count);
}
double E_ell(double k_1) //second elleptical integral
{
double E_count = 1.0, q_count = 1.0, m;
m = (1.0 - k_1)/(1.0 + k_1);
for (int i = 1; i < 1000; i += 2)
{
q_count = q_count*(i-2)*(i-2)/(i+1)/(i+1)*m*m;
E_count = E_count + q_count;
}
return (M_PI/2.0/(1.0 + m)*E_count);
}
main()
{
//Main parameters input
cout << "F = ? (nN)" << endl;
cin >> F;
F = F*1.0e-9;
cout << "E_sample = ? (GPa)" << endl;
cin >> E1;
E1 = E1*1.0e9;
cout << "E_probe = ? (GPa)" << endl;
cin >> E2;
E2 = E2*1.0e9;
D = 3.0/4.0*(1.0/E1 + 1.0/E2);
cout << "R_sample (cil, nm) = ?" << endl;
cin >> R1;
R1 = R1*1.0e-9;
cout << "R_probe (sph, nm) = ?" << endl;
cin >> R2;
R2 = R2*1.0e-9;
// probe/sample
A = 0.5*(1.0/R1 + 1.0/R2);
B = 0.5/R2;
// sample/substrate
// A = 0.5/R1;
// B = 0.5/(R2+2*R1);
//Exact solution
a = 0.1;
b = 1.0;
b_a_2 = b*b/a/a;
for (int i = 1; i < 1000; i++)
{
a_b = sqrt(1.0/b_a_2);
b_a_2 = A/B*K_ell(a_b)/E_ell(a_b);
b_a_2 = b_a_2 - A/B + K_ell(a_b)/E_ell(a_b);
};
b = 2.0*F*D/M_PI/B/(1.0 - a_b*a_b);
b = b*(K_ell(a_b) - E_ell(a_b));
b = pow(b, 1.0/3.0);
a = a_b*b;
h = 2.0*F*D/M_PI/b*K_ell(a_b);
P = F/M_PI/a/b;
cout << " Test: exactly" << endl;
a = a*1e9, b = b*1e9, h = h*1e9, P = P*1e-9;
cout << " a = " << a << " nm";
cout << " b = " << b << " nm" << endl;
cout << " h = " << h << " nm" << endl;
cout << " P = " << P << " MPa" << endl;
//Approximating solutions
// first: a sim b
cout << " Test: a sim b" << endl;
a = pow(2.0*F*D*(1.0/A + 1.0/B), 1.0/3.0);
a = a/(1.0 + A/B);
b = A*a/B;
h = 2.0*F*D/(a + b);
P = F/M_PI/a/b;
a = a*1e9, b = b*1e9, h = h*1e9, P = P*1e-9;
cout << " a = " << a << " nm";
cout << " b = " << b << " nm" << endl;
cout << " h = " << h << " nm" << endl;
cout << " P = " << P << " MPa" << endl;
//second: b>a
cout << " Test: b > a " << endl;
a = 0.1;
b = 1.0;
b_a_2 = b*b/a/a;
for (int i=1; i < 1000; i++)
{
C = log(4.0*sqrt(b_a_2)) - 1.0;
b_a_2 = A/B*C;
};
b = pow (2.0*C*F*D/M_PI/B, 1.0/3.0);
a = pow (4.0*B*F*F*D*D/M_PI/M_PI/C/A/A/A, 1.0/6.0);
P = F/M_PI/a/b;
h = 2*F*D*(C + 1.0)/M_PI/b;
a = a*1e9, b = b*1e9, h = h*1e9, P = P*1e-9;
cout << " a = " << a << " nm";
cout << " b = " << b << " nm" << endl;
cout << " h = " << h << " nm" << endl;
cout << " P = " << P << " MPa" << endl;
}
\tth_fileclosefebb02a4
Программа представляется достаточно простой: в качестве функций K_ell и E_ell реализованы процедуры вычисления полных эллиптических интегралов по асимпотикам, приведенным в [67]. Основная часть начинается с запроса параметров задачи, затем определяются параметры A и B (по формуле (2.12), если интересна область контакта образца и подложки, то следует использовать формулу (2.23)), после чего реализовано точное решение (по формулам (2.10), (2.11), (2.9). Для приближенного решения реализованы два случая: a ~ b (с использованием формул (2.21), (2.22) и a < b по формулам (2.17), (2.18). В качестве результата работы программа выдает найденные значения параметров области контакта (a и b), величину деформации h и контактное давление P.
Ниже приведен текст программы, написанной на языке С++.
#include <conio.h>
#include <stdio.h>
#include <math.h>
#include <fstream.h>
#include <stdlib.h>
#include <io.h>
#include <iostream.h>
double d_vtm, k, a, b, R, d, x, x1, x2;
int i, N, N1, N2, N_sph_bad, N_par_bad;
double sq(double z)
{
return (sqrt(fabs(R*R-(z-d)*(z-d))));
}
double func_sph(double z)
{
double temp = (b*b-(R-sq(z))*(R-sq(z)))*sq(z);
return (temp+z*(z-d)*(R-sq(z)));
}
double func_par(double z)
{
return (k*k*(z-d)*(z-d)*(z-d)*(z+d)+b*b);
}
double a_sph()
{
if (d > R) x1 = d-R;
else
{
if (b > R-sqrt(R*R-d*d)) return 0.0;
else x1 = 1e-6;
};
x2 = d;
for (int i = 1; i < N; i++)
{
x = (x1+x2)/2.0;
if (func_sph(x1)*func_sph(x) > 0.0) x1=x;
else x2=x;
};
if ((R-sq(x))*(R-sq(x))/b/b < 1.0)
{
return (x/sqrt(1.0-(R-sq(x))*(R-sq(x))/b/b));
}
else return 0.0; // flag: no solution
};
double a_par()
{
x1 = 1e-6;
x2 = d;
if (b > d*d*k) return 0.0;
else
{
for (int i = 1; i < N; i++)
{
x = (x1+x2)/2.0;
if (func_par(x1)*func_par(x)>0.0) x1 = x;
else x2 = x;
};
};
return sqrt (0.5*x*(x+d));
}
main()
{
ofstream outFile;
ifstream inFile;
outFile.open("RESULT.DAT");
inFile.open ("DATA.DAT");
cout <<" data file lenght?" << endl;
cin >> N1;
cout << " testing phase?" << endl;
cout << " 1 - yes" << endl;
cout << " 0 - no" << endl;
cin >> N2;
if (N2 == 1)
{
cout << " VTM width at half-height?" << endl;
cin >> d_vtm;
d_vtm = d_vtm/2.0;
R = (d_vtm*d_vtm-81.0)/18.0;
x = 0.5*(sqrt (d_vtm*d_vtm+648.0)-d_vtm);
k = sqrt (81.0-x*x)/(x-d_vtm)/(x-d_vtm);
cout << " for sphere, R=" << R << endl;
cout << " for parabola, k=" << k << endl;
cout << " Press any key..." << endl;
N2 = getch();
}
else
{
cout << " R (sphere)?" << endl;
cin >> R;
cout << " k (parabola)?" << endl;
cin >> k;
};
N = 1000;
cout << " sphere parabola" << endl;
for (int j = 1; j < N1; j++)
{
inFile >> d;
d = d/2.0;
inFile >> b;
b = b/2.0;
a = a_sph();
if (a != 0.0)
{
outFile << a << " ";
cout <<" a = "<< a << " ";
}
else
{
N_sph_bad++;
cout << " no" << " ";
outFile << "-" << " ";
};
a = a_par();
if (a != 0.0)
{
outFile << a << endl;
cout << " a = " << a << endl;
}
else
{
N_par_bad++;
cout << " no" << endl;
outFile << "-" << endl;
};
};
cout <<" Sphere wrong:" << N_sph_bad << endl;
cout <<" Parabola wrong:" << N_par_bad << endl;
};
\tth_fileclosefebb0304
Программа реализована по следующей схеме. Как показал анализ, погрешность определяемых значений существенно зависить от погрешности входных параметров (особенно измеренной ширины профиля на полувысоте d). Для уменьшения этой погрешности необходим набор достаточной статистики для значений входных параметров - h (высоты) и d. Поэтому в программе ввод входных данных осуществляется в виде файла ("data.dat") состоящего из двух столбцов значений для h и для d соответственно. На этапе вычисления программа запрашивает о величине набранной статистики (количестве строк в файле "data.dat"). Вывод значений также осуществляется в файл (result.dat) в две колонки (для аппроксимации иглы сферой и параболой соответственно).
Функции func_sph и func_par построены по формулам (2.34, 2.36), далее, функции a_sph и a_par реализуют решение уравнений (2.34, 2.36) относительно x0 и, по найденному x0, вычисление искомого значения a по формулам (2.35, 2.37). Одновременно, эти функции осуществляют проверку выполнения условий (2.38, 2.39) отсутствия решений; подсчет числа случаев отсутствия решений ведут счетчики N_sph_bad и N_par_bad для сферической и параболической аппроксимаций соответственно.
Работа программы по желанию может быть предварена тестирующей фазой, в которой по измеренным параметрам АСМ-профиля тест-объекта (вируса табачной мозаики) определяются аппроксимирующие параметры иглы (R и k), минуя тестирующую фазу эти параметры можно задать вручную.
Основные параметры зондов промышленных кантилеверов (Nanoprobe, Digital Instruments, США) для атомно-силовых микроскопов приведены в таблице A.1. Эти кантилеверы использовались в настоящей работе.
Модель | NP-S | FESP(W) | TESP(W) |
Oxide Sharpened Silicon Nitride | Force Modulation Etched Silicon | Tapping Mode Etched Silicon | |
Жесткость | 0,06 Н/м, 0,12 Н/м, 0,38 Н/м, 0,58 Н/м | 1-5 Н/м | 20-100 Н/м |
Резонансная частота | 5-50 КГц | 60-100 КГц | 200-400 КГц |
Радиус кривизны кончика иглы | 5-40 нм | 5-10 нм | 5-10 нм |
Высота иглы | 3 мкм | 10-15 мкм | 10-15 мкм |
Угол раствора конуса иглы | 35° | 18° Ъ25° Ъ10° | 18° Ъ25° Ъ10° |
Отражающее покрытие | Золото | нет | нет |
Длина кантилевера | 120 мкм, 200 мкм | 225 мкм | 125 мкм |
Ширина кантилевера | 15 мкм, 40 мкм | 30 мкм | 30-40 мкм |
Толщина кантилевера | 0,4-0,6 мкм | 2-3 мкм | 3,5-5 мкм |
Форма кантилевера | V-образная | прямая | прямая |