Расчет сил локального трения при атомно-силовой микроскопии поверхности металлов
Поступила в редакцию: 12 января 1996 г.
Выставление онлайн: 20 января 1997 г.
*Введение В последние годы значительное внимание уделяется измерению горизонтальной составляющей силы, действующей на острие атомно-силового микроскопа (АСМ) со стороны исследуемой поверхности (см., например, [1-5]). Это связано прежде всего с тем, что трение само по себе является одним из наиболее общих и сложных физических явлений. Поэтому детальное изучение соответствующих процессов с помощью АСМ представляет несомненный интерес. Кроме того, силы трения, возникающие между зондирующим острием и исследуемой поверхностью, могут оказывать значительное влияние на пространственное разрешение АСМ и топографию получаемых с помощью АСМ изображений [4]. В то же время получение в процессе сканирования двух независимых рельефов - поверхности постоянной вертикальной составляющей силы и соответствующего рельефа для ее горизонтальной составляющей - дает дополнительную информацию об исследуемом образце. Технические трудности, связанные с одновременным измерением вертикальной и горизонтальной составляющих силы, действующей на зондирующее острие, к настоящему времени в значительной мере преодолены (см., например, [5,6]). Однако теоретическая интерпретация полученных результатов требует дополнительных исследований. В настоящей работе предлагается подход, позволяющий не только описать сканирование в режиме сил отталкивания, но и вычислить значения горизонтальных составляющих силы, действующей на острие АСМ. В качестве примера приводятся результаты расчетов для сканирования над поверхностью решетки плотной упаковки атомов, как невозмущенной, так и содержащей точечные дефекты различных типов. Подобные расчеты, как указывалось ранее [7,8], могут быть использованы при диагностике точечных дефектов на поверхности металлов методом АСМ. skip *Методика расчета В предположении о парном взаимодействии атомов острия и поверхности потенциальную энергию взаимодействия можно записать в виде U(x,y,z)=i=1nj=1m u(xi,yi,zi,x+xj,y+yj,y+yj,z+zj). (1) Здесь x, y, z - координаты кончика острия относительно исследуемой поверхности; индекс i нумерует атомы поверхности с координатами xi, yi, zi; индекс j нумерует атомы острия, имеющие координаты xj, yj, zj, в системе координат, начало которой совпадает с кончиком острия. Если пренебречь смещением атомов поверхности и острия в процессе сканирования (что во всяком случае, безусловно, справедливо при силе сканирования F0~10-9 H [9,10]), вертикальную и горизонтальную (вдоль оси сканирования x) составляющие силы, действующей на острие, то можно найти в результате дифференцирования потенциальной энергии (1) Fz(x,y,z)=-(d)/( d z)U(x,y,z), Fx(x,y,z)=-(d)/(d x)U(x,y,z). (2) При сканировании на постоянной высоте z=const сила Fx в (2) представляет собой собственно силу локального (атомного) трения. Сканирование в АСМ, однако, проводится вдоль поверхности постоянной силы Fz=const (точнее, постоянной вертикальной составляющей силы). Соответственно измеряемая при этом горизонтальная составляющая силы представляет собой функцию Fx(x,y)=Fx(x,y,z)|Fz(x,y,z)=const. (3) Для того чтобы найти величину (3), необходимо в каждой точке (x,y) решить относительно z уравнение -(d)/(d z)U(x,y,z)=F0, (4) где F0 - сила сканирования. [!tb] Двумерное изображение поверхности постоянной силы ( а) и профиля горизонтальной составляющей силы Fx ( б) при сканировании острия АСМ над решеткой плотной упаковки атомов. а - hmax=0.5a, hmin=0.367a; б - Fxmax=0.473F0, Fxmin=-0.473F0, При этом будет получена поверхность постоянной силы z=z(x,y). Затем в каждой точке этой поверхности следует вычислить силу Fx(x,y,z(x,y)). Результаты такого расчета, естественно, чувствительны к выбору потенциала парного взаимодействия, особенно в режиме сил отталкивания, к выбору модели атомного строения острия. Простейшая модель одноатомного острия (см., например, [7,9,10]) позволяет при любых начальных высотах (т. е. при любых силах сканирования) достаточно хорошо описать поверхность постоянной силы в области непосредственно над атомами исследуемой поверхности. Однако при достаточно малых начальных высотах найти непрерывную поверхность постоянной силы в области между атомами оказывается невозможным [10,12,13]. Для создания более сложных моделей кончика острия, позволяющих адекватно описывать экспериментальные результаты, используются два пути. Первый связан с моделированием той или иной кластерной структуры кончика острия (см., например, [12-15]). Второй путь состоит в представлении кончика остиря в виде некоторого квазиатома и подборе подходящего модельного потенциала взаимодействия этого квазиатома с атомами поверхности [16]. При этом, однако, вполне может оказаться, что модельный потенциал, подобранный так, чтобы обеспечить соответствие теоретических и экспериментальных результатов для поверхности постоянной силы, не обеспечивает адекватное описание горизонтальной составляющей силы, действующей на острие. Поэтому более естественным представляется моделирование кончика острия атомным кластером, соответствующим структуре материала (решетка плотной упаковки [12,13] либо тетрагональная структура [15]). [!tb] skip Двумерное изображение поверхности постоянной силы ( а) и профиля горизонтальной составляющей силы Fx при сканировании острия АСМ над вакансией в решетке плотной упаковки атомов: а - hmax=0.5a, hmin=0.284a; б - Fxmax=0.588F0, Fxmin=-0.456F0. [!b] Двумерное изображение поверхности постоянной силы ( a) и профилей горизонтальных составляющих силы Fx ( б) и Fy ( в) при сканировании острия АСМ над дивакансией в решетке плотной упаковки атомов. а - hmax=0.5a, hmin=0.34a; б - Fxmax=0.508F0, Fxmin=-0.532F0.; в - Fymax=0.585F0, Eymin=-0.585F0. В настоящей работе используется модель параболического острия [8,17], на кончике которого находится кластер из 7 атомов [12,13], один из которых опущен на расстояние h=0.22a относительно плоскости остальных шести, расположенных в виде шестиугольника со стороной b=0.85a (a - равновесное расстояние между атомами в решетке плотной упаковки). Такая модель острия позволяет непротиворечиво описать процесс сканирования на начальных высотах d~(0.5-0.6)a, когда силами притяжения можно пренебречь по сравнению с силами отталкивания. В этом случае парный потенциал в (1) можно взять, например, в виде [18] [b] u&(xi,yi,zi,x+xj,y+yj,z+zj) &= (alpha)/([(x+xj-xi)2+(y+yj-yi)2+ (z+zj-zi)2]6). (5) Для того чтобы избежать зависимости результатов от значения полуэмпирической постоянной alpha, удобно все расчеты проводить в относительных величинах, измеряя все расстояния в единицах межатомного равновесного расстояния a, а силы в единицах силы сканирования F0. Кроме того, поверхность постоянной силы, полученная при решении уравнения (1), и соответствующая горизонтальная составляющая силы зависят от ориентации острия относительно кристаллографических осей поверхности. При этом методика расчетов существенно зависит от того, как производится сканирование. Если при сканировании ориентация микрокластера на кончике острия остается строго фиксированной относительно поверхности, то необходимо проводить все расчеты именно для этой фиксированной ориентации. Если же в процессе сканирования ориентация микрокластера относительно поверхности случайным образом меняется, то необходимо проводить соответствующее усреднение [12,13]. Ниже приводятся результаты расчетов усредненных поверхностей постоянной силы и соответствующих средних значений горизонтальной составляющей силы для решетки плотной упаковки, моделирующей поверхность металлов. *Сканирование над решеткой плотной упаковки На рис. 1, а приведены результаты расчета усредненной поверхности постоянной силы при сканировании острия АСМ над невозмущенной решеткой плотной упаковки. Двумерные изображения здесь и далее представляют собой сечения соответствующих поверхностей, полученных при решении уравнения (4), десятью горизонтальными плоскостями, проведенными с шагом Delta z=(zmax-zmin)/10. Более светлый тон соответствует большей высоте, более темный - меньшей. Во всех рассматриваемых случаях zmax=0.5a. На рис. 1, б приведено двумерное изображение профиля горизонтальной силы (3). Белый тон на двумерных изображениях силы соответствует максимальному положительному значению Fxmax, когда горизонтальная сила действует на острие в направлении сканирования. Черный тон соответствует минимальному, т. е. максимальному по модулю отрицательному, значению Fxmin, когда горизонтальная сила действует на острие в направлении, противоположном направлению сканирования. Каждое изменение тона от более светлого к более темному соответствует уменьшению силы на Delta Fx=(Fxmax-Fxmin)/10. Как видно из сравнения рис. 1, а и б, горизонтальная составляющая силы обращается в нуль непосредственно над атомами исследуемой поверхности и в точках между ними, что согласуется с экспериментальными результатами (см.: например, [6]). Полученная глубина рельефа для горизонтальной составляющей (Fxmax-Fxmin)~ F0 также хорошо согласуется с экспериментальными данными [6]. *Вакансия и дивакансия в решетке плотной упаковки Рассмотрим теперь некоторые типы точечных дефектов в решетке плотной упаковки. На рис. 2 приведены результаты расчетов для дефекта типа вакансии: двумерное изображение усредненной поверхности постоянной силы (рис. 2, а) и соответствующего профиля горизонтальной составляющей силы Fx (рис. 2, б). Сравнение рис. 1, а, б с рис. 2, а, б соответственно наглядно демонстрирует возможность обнаружения на поверхности образца дефектов типа вакансии по данным АСМ. Для дивакансии двумерное изображение поверхности постоянной силы приведено на рис. 3, а. Заметим, что при сканировании острия АСМ вблизи дивакансии вдоль взаимно перпендикулярных осей x и y горизонтальные составляющие совершенно различны и уже не могут быть выражены друг через друга с использованием свойств симметрии решетки, как в случаях невозмущенной решетки и вакансии. Поэтому были подробно исследованы режимы сканирования вдоль обеих осей. Полученные профили горизонтальных компонент силы Fx и Fy приведены на рис. 3, б и в соответственно. Полученные результаты ясно показывают, что наличие на исследуемой поверхности точечных дефектов заметно изменяет не только характер поверхности постоянной силы, но и профиль горизонтальной составляющей силы, действующей на острие в процессе сканирвания в режиме сил отталкивания. Следовательно, измерение горизонтальной составляющей силы дает дополнительную информацию для диагностики точечных дефектов методом АСМ.
- Mate C.M., McClellard, Erlandsson R. // Phys. Rev. Lett. 1987. Vol. 59. N 4. P. 1942-1948
- Marti O., Colchero J., Mlynek J. // Nanotechnology. 1990. Vol. 1. N 1. P. 141--144
- Meyer G., Amer M.N. // Appl. Phys. Lett. 1990. Vol. 57. N 20. P. 2089--2091
- Моисеев Ю.Н., Панов В.И., Савинов С.В. // Письма в ЖТФ. 1991. Т. 17. Вып. 10. С. 24--28
- Marti O. // Physica Scripta. 1993. Vol. 149. N 1. P. 599--604
- Marti O. // Surface and Coating Technology. 1993. Vol. 62. P. 510--516
- Благов Е.В., Климчицкая Г.Л., Мостепаненко В.М. и др. // Письма в ЖТФ. 1993. Т. 19. Вып. 8. С. 73--78
- Благов Е.В., Климчицкая Г.Л., Мостепаненко В.М. и др. // Письма в ЖТФ. 1994. Т. 20. Вып. 1. С. 71--77
- Gould S.A.C., Burke K., Hansma P.K. // Phys. Rev. B. 1989. Vol. 40. N 8. P. 5363--5366
- Abraham F.F., Batra I.R. // Surf. Sci. 1989. Vol. 209. N 1/2. P. L125--L132
- Albrecht T.R., Quate C.F. // J. Appl. Phys. 1987. Vol. 62. N 7. P. 2599--2602
- Благов Е.В., Климчицкая Г.Л., Лобашев А.А., Мостепаненко В.М. // Письма в ЖТФ. 1995. Т. 21. Вып. 3. С. 73--80
- Blagov V., Klimchitskaya G., Lobashev A., Mostepanenko V. // Surf. Sci. In press. Preprint Friedmann Laboratory Publishing. N FL-030495. St. Petersburg, 1995
- Abraham F.F., Batra I.P., Ciraci S. // Phys. Rev. Lett. 1988. Vol. 60. N 13. P. 1314--1317
- Tang H., Joachim C., Devillers J. // Surf. Sci. 1993. Vol. 291. N 3. P. 439--450
- Overnay G., Tomanek D., Zhong W. et al. // J. Phys. C. 1994. Vol. 4. N 17. P. 4233--4249
- Bordag M., Klimchitskaya G., Mostepanenko V. // Surf. Sci. 1995. Vol. 328. N 1/2. P. 129--134
- Israelachvili J.N. Intermolecular and Surface Forces. New York; London; Academic Press, 1985
Подсчитывается количество просмотров абстрактов ("html" на диаграммах) и полных версий статей ("pdf"). Просмотры с одинаковых IP-адресов засчитываются, если происходят с интервалом не менее 2-х часов.
Дата начала обработки статистических данных - 27 января 2016 г.