Особенности плавления металлов при действии ультракоротких лазерных импульсов Текст научной статьи по специальности «Физика»

УДК 538.95

Е. Б. Яковлев, В. В. Свирина, О. Н. Сергаева

ОСОБЕННОСТИ ПЛАВЛЕНИЯ МЕТАЛЛОВ ПРИ ДЕЙСТВИИ УЛЬТРАКОРОТКИХ ЛАЗЕРНЫХ ИМПУЛЬСОВ

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

Ключевые слова: лазерное плавление, вакансионная модель плавления, ультракороткие лазерные импульсы, численное моделирование.

Введение. В последние годы произошел возврат интереса к изучению лазерно-индуцированных фазовых переходов (плавление, кристаллизация) в различных материалах. Было выяснено, что при быстром лазерном нагревании кинетика таких переходов существенно отличается от их протекания в условиях медленного нагревания [1].

Как показывают исследования, при действии ультракоротких лазерных импульсов на металлы плавление происходит при температурах, превышающих температуру плавления. Кроме того, если, в соответствии с существующими представлениями, фронт плавления распространяется с поверхности, то скорость его продвижения оказывается больше, чем допустимая скорость распространения возмущения в металлах — скорость звука [2].

Существует множество теорий плавления [см., например, 3], однако ни одна из них не дает возможности проанализировать кинетику плавления при действии ультракоротких лазерных импульсов, когда необходимо детально знать механизм перехода вещества из твердого состояния в жидкое. Например, предложен механизм гомогенной нуклеации, заключающийся в том, что при значительном перегреве в объеме материала могут возникнуть зародыши жидкой фазы, произойдет объемное плавление [2]. Против этого можно возразить: объем расплава большинства металлов больше, чем объем твердой фазы, поэтому при возникновении зародыша расплава в нем должно вырасти давление, что приведет к увеличению температуры плавления: вероятность образования зародыша резко снижается.

Теории плавления не могут также объяснить, как переносится энергия к фронту плавления при отсутствии градиента температуры в слое расплава, и когда начинается увеличение температуры расплава на поверхности, т. е. когда возникает градиент температуры в расплаве.

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

В работе [4] предложена модифицированная вакансионная модель плавления, которая может быть применена для описания процесса плавления металла при действии ультракоротких лазерных импульсов.

Основные положения, описывающие построение модели плавления, базирующейся на анализе изменения концентрации вакансий, предложены Френкелем [5, 6]. Эти положения можно сформулировать следующим образом:

— явление плавления связано с резким увеличением числа вакансий в кристалле при увеличении температуры, при этом необходимо учитывать уменьшение энергии образования вакансий при росте их концентрации;

— с увеличением концентрации вакансий связана скрытая теплота плавления и изменение объема материала при плавлении;

— концентрация вакансий увеличивается за счет их диффузии с поверхности, поэтому плавление происходит с поверхности.

Согласно этой модели нарушение дальнего порядка при плавлении происходит за счет деформации кристаллической решетки вакансиями (см. рис. 1, где а — единичная вакансия,

б — структура, искаженная взаимодействую-

а)

б) ^

щими вакансиями).

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

Состояние деформированного вакансиями кристалла описывается его свободной энергией Е :

Е = Еп-ап [1-ехр (-пу))-кТп 1п ,

где Е — энергия образования вакансии; п — концентрация вакансий; а — энергия взаимодействия вакансий; у — объем, в котором взаимодействуют вакансии; N — концентрация узлов кристаллической решетки; Т — температура; к — постоянная Больцмана; е — основание натурального логарифма.

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

пп — число вакансий во взаимодействующих парах); пп можно определить как произведение

общего числа вакансий п на вероятность того, что около вакансии в некотором объеме у , определяемом радиусом взаимодействия вакансий, окажется еще хотя бы одна вакансия: т.е. по

закону Пуассона пп = п [1-exp (-пу)) .

Уравнение диффузии вакансий при плавлении можно представить в следующем виде [4]:

^дп = д (р^дп^ & дх ^ дх /

где д = дЕ/ дп — химический потенциал вакансий, Е — коэффициент их диффузии; граничные условия сводятся к определению теплового потока на границе: q = Ег(п) дп(0)/дх (в(п) — энергия образования и растворения вакансии).

Основные результаты. Результаты численного решения задачи нагревания серебра лазерными импульсами микросекундной длительности приведены на рис. 2, где представлена зависимость температуры поверхности Т, глубины слоя расплава h и отражения R от времени (кривые 1—3 соответственно). Фрагмент экспериментальной зависимости R(t), рассмотренной в работе [7], показан на рис. 3.

Я, o.e.

Т, К

1800 Н 1600 1400 Н 1200 1000 -800 600 400

h, см 1,0-10-4 8,0-10-5 6,0-10-5 4,0-10-5 2,0-10~5 0,0

2,5

3,0 3,5

Рис. 2

4L

0,8 0,7 0,6 0,5 0,4 0,3

R, о.е.

0,5

4,0 t, сх10-

0,2

Рис. 3

0,4 t, мкс

При вычислениях приняты следующие значения параметров: теплопроводность

3

Х=4,29 Вт/(см-К); объемная теплоемкость с=2,467 Дж/(см-К); а=0,5756 эВ; Е=1,1 эВ;

#=5,856-1022 1/см3; энергия активации диффузии вакансии Ed =0,66 эВ; V=1,421-10 20 см3;

7 2

максимальная плотность мощности поглощенного лазерного излучения ^0 =7-10 Вт/см ; начальная температура образца, согласно работе [7], принята равной 573 К.

( ^—^шах )

20

Форма импульса определена в виде q0 ехР

t

где t =0,3 мкс, tmax =0,6 мкс.

Анализ рис. 2 и 3 показывает удовлетворительную согласованность расчета с экспериментом.

Темп нагревания при лазерном воздействии определяется большой скоростью подвода энергии: так, например, при действии пикосекундных импульсов скорость передачи энергии от электронов, нагретых до температур порядка нескольких десятков тысяч градусов Кельвина, к решетке может достигать примерно 1014 Вт/см3. Процесс плавления в при этом отличается от хорошо известного плавления в квазистационарных условиях.

Температуру металла при нагревании импульсом пикосекундной длительности можно определить из решения системы уравнений:

дТе

С

dt дх

T _д_(

1 dt dx

д (x дТе

дх

=a

, q0t ехР (_ap x_t/t0 ) /10 _ße1 (Te _T1);

- f ^ (Te _T1

х дт,(0,1), дт,(0,)_0

дх дх

где Те, ] , Хе, X,, се, с, — температура, теплопроводность и теплоемкость электронного газа (индекс „е") и кристаллической решетки (индекс „/") соответственно; Рег- — коэффициент теплообмена между электронным газом и решеткой; а pqot ехр (-а рх - ^^ ) / to — поглощенная плотность мощности лазерного излучения; а р — коэффициент поглощения; ^ — характерная длительность лазерного импульса.

0

3

При вычислениях принято: Хе =2,37 Вт/(см-К); = 1 Вт/(см-К); се =0,04035 Дж/(см -К);

3 -5 —1 22

сI =2,43 Дж/(см -К); ар =10 см ; £=1,1 эВ; концентрация атомов и0=5,856-10 ; ел=0,66 эВ; длительность воздействия т=4-10-11 с; д0=1,7-10п Вт/см2.

Результаты численного решения задачи нагревания серебра импульсами пикосекундной длительности графически представлены на рис. 4 и 5.

Т, К

И, нм

4000 3000 2000

1000 0

10 Рис. 4

и пс

600 400 200

0

10 15 Рис. 5

20

и, пс

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

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

Определим скорость подачи энергии, при которой плавление металла при температуре его плавления не происходит:

— при температуре плавления

(

ЛТ

ктт ехр

Е.

Л

Л

кТт

->-

Л т0(Е + Её ) — при температурах, превышающих Тт ,

п2

->-ехр

т )

ЛТ

Л %0(Еу + ЕЛ)

-Е± + К (Т-Тт ) кТ кТТт

Здесь Еу — энергия образования вакансий с учетом их взаимодействия, т0=10 14 с.

Зависимость температуры кристаллической решетки от времени характеризуется выражением

Т =-

дт

(

(

1- ехр

Х (1+сг1се ) У|

))

1 (се + с1 )И

Время достижения температуры Т поверхности определяется как

с Г ТИ(с + с )\

и =--, ^ . . 1п

1 ТИ ( се + сг) дт

Х (1+^1се )

Температура, при которой начнется плавление, определяется из уравнения

кТ 2

т0(Еу + ЕЛ )

ехр

-+ Е (Т-Тт ) кТ кТТт

с1сеИ

\_дт-ТИ(се + с )] .

5

5

В результате можно найти время проплавления

где 5 — постоянная кристаллической решетки.

Если температура начала плавления такова, что дальний порядок нарушается вследствие генерации френкелевских пар, то процесс плавления начинается одновременно по всему объему прогретого материала. Время плавления в этом случае определяется временем генерации достаточного числа дефектов решетки:

Если принять v таким же, как при вакансионном плавлении, то можно оценить энергию взаимодействия а^ френкелевских пар:

где EF — энергия образования френкелевских пар.

Заключение. Приведенные экспериментальные данные и численные расчеты подтверждают возможность применения вакансионной модели плавления к описанию процессов, происходящих при плавлении материалов ультракороткими лазерными импульсами. При нагревании металлов импульсами микросекундной длительности результаты численного расчета хорошо согласуются с экспериментальными данными. Используя вакансионную модель плавления, можно однозначно объяснить постоянство температуры поверхности при лазерном плавлении и наличие потока энергии к фронту плавления без градиента температуры в расплавленном слое. Такой режим возможен до тех пор, пока плотность потока мощности лазерного излучения на поверхности не превысит некоторого критического значения, определяемого скоростью генерации вакансий на поверхности расплава и величиной диффузионного потока вакансий, который определяется градиентом концентрации вакансий на фронте плавления. Поэтому глубина проплавления будет определяться плотностью потока лазерного излучения и свойствами материала.

Работа выполнена при поддержке Российского фонда фундаментальных исследований, гранты 09-02-00932-а и 09-02-01065-а, и государственного контракта № П1134.

1. Rethfeld B., Socolowski-Tinten K., von der Lin de D., Anisimov L. D. Ultrafast thermal melting of laser-excited solids by homogeneous nucleation // Phys. rev. B. 2002. Vol. 65, N 9. P. 1—4.

2. Ретфельд Б., Темнов В. В., Соколовски-Тинтен К. и др. Сверхбыстрое термическое плавление твердых тел под воздействием фемтосекундных лазерных импульсов // Опт. журн. 2004. Т. 71, № 6. С. 18—23.

3. Уббелоде А. Р. Расплавленное состояние вещества. М.: Металлургия, 1982. 256 с.

4. ЯковлевЕ. Б. Перегрев твердых тел при плавлении // Изв. АН СССР. Сер. физ. 1989. Т. 53, № 3. С. 591—594.

5. Френкель Я. И. Собрание избранных трудов. М. — Л.: Изд-во АН СССР, 1958. Т. 2.

6. ФренкельЯ.И. Статистическая физика. М. — Л.: Изд-во АН СССР, 1948.

7. Анисимов С. И., Имас Я. А., Романов, Г. С., Ходыко Ю. В. Действие излучения большой мощности на металлы. М.: Наука, 1970. 272 с.

a F

список литературы

Евгений Борисович Яковлев

Вера Владимировна Свирина

Ольга Николаевна Сергаева

Сведения об авторах д-р техн. наук, профессор; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра лазерных технологий и экологического приборостроения; E-mail: [email protected]

студентка; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра лазерных технологий и экологического приборостроения; E-mail: [email protected]

студентка; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра лазерных технологий и экологического приборостроения; E-mail: [email protected]

Рекомендована кафедрой лазерных технологий и экологического приборостроения

Поступила в редакцию 07.12.09 г.