CC BY
69
(1-4) важно иметь в виду следующие два обстоятельства. Во-первых, каждый из четырех рассмотренных механизмов может реализоваться при более, чем одном наборе соответствующих одноэле-ктронных состояний. Это не означает, что вероятность перехода будет равна сумме вероятностей из-за наличия интерференционных слагаемых,
х м а
= Х |маТ • Последние могут либо суще-
ственно уменьшить результирующую вероятность, если частные значения оценочных амплитуд М™
имеют противоположные знаки, либо наоборот -увеличить, если амплитуды одного знака. Во-вторых, в совокупности вкладов в полную амплитуду могут оказаться слагаемые, отвечающие разным из рассмотренных модельных механизмов. И в этом случае интерференция различных парциальных каналов может существенно уменьшить или увеличить вероятность перехода.
В обоих случаях уже нельзя, строго говоря, приписывать переходу тот или иной механизм в качестве доминирующего - требуется учет всех корреляционных взаимодействий.
2
СПИСОК ЛИТЕРАТУРЫ
1. Килин В.А. Методика расчета амплитуд и вероятностей переходов в атомах с учетом корреляций в рамках многочастичной нестационарной теории возмущений // Известия Томского политехнического университета. —2004. —Т. 307. —№ 6. — C. 6-11.
2. Купляускис З.Й., Купляускене Ф.В., Тутлис В.И. Об изучении возбужденных состояний атомов с использованием неортогональных радиальных орбиталей // Известия вузов. Физика. — 1981. — № 3. —С. 7—11.
3. Ландау Л.Д., Лифшиц Е.М. Квантовая механика. — М.: Наука, 1974. —752 с.
4. Килин В.А., Ли И.С. Двойной Оже-распад в рамках МТВ // Известия вузов. Физика. —1989. — № 7. — С. 78—82.
5. Carlson T.A., Krause M.O. Experimental Evidence for Double Electron Emission in an Auger Process // Phys. Rev. Lett. —1965. — V. 14. — № 11. —P. 390—392.
6. Kilin VA., Ehresmann A., Vollveiler F., Schartner K.-H., Schmoran-zer H. Perturbation theory study of triple photoionization. I. Two-step approximations in triple photoionization of Kr in the exciting photon energy region of KrI 3d9 np resonances (90 eV-100 eV) // J. Phys. B.:
Atom. Mol. Opt. Phys. -1997. -V. 30. -№ 24. - P. 5715-5727.
7. Kilin V.A., Lee I.S. Participator-spectator-vacancy satellites in Auger spectra. Probabilities and angular distribution // Proc. of XXII EGAS, Uppsala, Sweden, 1990. - P. 629-631.
8. Amusia M.Ya., Kilin V.A., Ehresmann A., Schmoranzer H., Schartner K.-H. Double-autoionization decay of resonantly excited single-electron state // J. Phys. B: Atom. Mol. Opt. Phys. -1993. -V. 26. - № 26. -P. 1281-1300.
9. Kilin V.A., Ehresmann A., Schmoranzer H., Schartner K.-H. Indirect observation of new three electron Auger transitions by PIFS // Abstracts IV ECAMP, Riga, Latvia, 1992. - P. 167.
10. Kilin V.A., Kharlova A.N., Ehresmann A., Schmoranzer H., Schartner K.-H. Competition between non-correlative visible and correlative fluorescence transitions in KrIII // J. Phys. B.: Atom. Mol. Opt. Phys. - 1995 -V. 28. - № 22. - P. 4723-4732.
11. Ehresmann A., Kilin V.A., Chernysheva L.V., Schmoranzer H., Amusia M.Ya., Schartner K.-H. Three-electron radiative transitions // J. Phys. B: Atom. Mol. Opt. Phys. -1993. -V. 26. -№ 5. - P. L97-L102.
12. Carter S. L., Kelly H.P. Double photoionization of neon and argon// Phys. Rev. A. -1977. -V. 16. - № 4. - P. 1525-1534.
УДК 621.314
ДИНАМИКА ПАРАМЕТРА ПОРЯДКА В СТРУКТУРНОНЕУСТОЙЧИВОМ КРИСТАЛЛЕ
Е.Е. Слядников
Отдел проблем информатизации ТНЦ СО РАН. г. Томск E-mail: [email protected]
Теоретически показано, что мягкая мода в окрестности структурного перехода исходная структура — предпереходное состояние — конечная структура, вынужденного как изменением температуры, так и внешней силы, переторможена и, следовательно, динамика параметра порядка имеет релаксационный характер.
1. Введение
Экспериментально обнаружено, что кристалл, неустойчивый относительно структурного перехода исходная - конечная структура, вызванного как изменением температуры, так и внешней силы, в окрестности структурного перехода находится в пред-переходном состоянии [1, 2]. Для теоретического описания структурного перехода и предпереходно-го состояния сформулирована микроскопическая модель [3], описывающая структурнонеустойчивый
кристалл как квантовую систему псевдоспинов. Причина возникновения предпереходного состояния в кристалле в том, что изменение внешнего воздействия (температуры, механической силы), стимулирующее структурный переход исходная -конечная структура, уменьшает площадь горба, разделяющего минимумы двухямного кристаллического потенциала атома. Это приводит к существенному увеличению квантового туннелирования атома и уменьшению асимметрии двухямного кристаллического потенциала. Возникает неустойчивость
состояния исходной кристаллической решетки с асимметричным двухямным потенциалом относительно возникновения предпереходного состояния решетки с симметричным двухямным потенциалом. В результате на верхней границе неустойчивости (при температуре T+) возбуждается мягкая коллективная мода - псевдоспиновая волна, переводящая кристалл из исходной структуры в предпе-реходное состояние. При дальнейшем уменьшении внешнего воздействия на нижней границе неустойчивости (при температуре T-) происходит конденсация мягкой псевдоспиновой волны, которая переводит кристалл из предпереходного состояния в конечную структуру. В рамках приближения молекулярного поля (ПМП) [3] можно изучить только термодинамические (стационарные) свойства структурного перехода и предпереходного состояния. Поэтому для изучения динамических свойств структурнонеустойчивого кристалла, например, частоты мягкой псевдоспиновой волны, необходимо выйти за рамки ПМП и использовать приближение хаотических фаз (ПХФ) [4].
В настоящей работе с помощью ПХФ показано, что при почти непрерывном переходе исходная структура - предпереходное состояние - конечная структура собственная коллективная мода системы псевдоспинов смягчается и спонтанно возникает разница заполнения левой и правой ям потенциального рельефа, то есть происходит Бозе-конденсация псевдоспиновой волны. Если не учитывать релаксацию системы псевдоспинов в окрестности структурного перехода, то это приводит к динамике с незаторможенной мягкой модой. Однако, хорошо известно, что вблизи структурного перехода динамика параметра порядка является колебательно-релаксационной и описывается уравнением типа Ландау-Халатникова [5], причем мягкая мода, как правило, переторможена [6]. В этой связи возникает задача учета релаксации системы псевдоспинов в окрестности структурного перехода исходная структура -предпереходное состояние - конечная структура.
2. Гамильтониан системы и основные уравнения
Поскольку нас интересует отклик системы псевдоспинов на компоненту механического поля напряжений Q;(t), изменяющую асимметрию двухямно-го потенциального рельефа и зависящую от времени и пространственных координат, запишем гамильтониан квантовой системы псевдоспинов в виде [3]
H = -X {-ЙюД' + [(1/2) X j +
• i
+(1/3) X hIiJ-sm + йц (t)] }. (l)
Здесь Sf, Sf, Sf - операторы Паули для спина Ню0 - расщепление энергий четного и нечетного состояний атома в двухямном потенциале, J (hlljm) - константа двухчастичного (трехчастичного) взаимодействия псевдоспинов. Вследствие временной зависимости H средние значения спино-
вых переменных также зависят от времени, поэтому будем обозначать их S"(0=Sf>;.
Для исследования динамики системы псевдоспинов, описываемой гамильтонианом (1), используем стандартную методику построения теории хаотических фаз [4]. Гейзенберговские уравнения движения для средних значений спиновых операторов (в системе единиц, где Й=1) имеют вид
dS. (t) / dt = -i < [S., H] > . (2)
В приближении хаотических фаз, являющемся прямым обобщением ПМП на случай задач с временной зависимостью уравнения движения (2) приобретают вид
dSt (t)/ dt = Si (t) x h (t), (3)
где зависящее от времени молекулярное поле h;(t) дается выражением
И.(t) = -[д < H >t /д§ (t)]. (4)
Эти уравнения эквивалентны уравнениям для свободной прецессии псевдоспина вокруг мгновенного значения молекулярного поля в данном узле. Однако в ур. (3, 4) не содержится релаксация системы псевдоспинов, которая, как было сказано выше, играет принципиальную роль в окрестности структурного перехода. Остается, таким образом, дополнить (3, 4) членами, описывающими релаксацию. Вначале заметим, что (3) обладает хорошо известным интегралом движения
Sf = (SX )2 + (sy )2 + S )2 = (S01 )2, (5)
где S0x - термодинамически равновесное значение Sx в предпереходном состоянии кристалла.
Структурный переход и динамика параметра порядка в его окрестности - это коллективные явления. Отсюда можно предположить, что эффекты смягчения собственной коллективной моды и релаксации в окрестности перехода также имеют коллективный (когерентный) характер и поэтому не нарушают закон сохранения (5), как это имеет место в теории ферромагнетизма [7]. Следуя и дальше отмеченной аналогии, перепишем (3) с учетом релаксации, опуская индекс узла i, в виде
dS / dt = S x И + R, (6)
R = -X[S x[S xh ]] = -X[(Sh)S - (S0 z )2 h ]], (7)
где R - релаксационный член, а Х- феноменологический параметр. В [8] было показано, что в низкочастотном приближении механизмы нелинейности, с одной стороны, дисперсии и диссипации -с другой, могут быть учтены аддитивным образом. Используем низкочастотное приближение, выражающееся в нашем случае неравенством
ю < min(fflc,т_1). (8)
Здесь ю - характерная частота динамического процесса, юс - частота мягкой коллективной моды, т - время коллективной релаксации параметра по-
рядка Sz. Принимая во внимание замечание о возможности аддитивного учета нелинейности, с одной стороны, дисперсии и диссипации - с другой, в приближении (8) линеаризуем (7) по 31, 3У, 3х, П. Введем эффективную статическую восприимчивость %4, определяемую соотношением
я0 = Хе/Ь°, (9)
где £0 и 6 0 - соответствующие равновесные значения 3 и к . Из (3, 4, 9) при й3/сИ=0 следует, что
Xе =-Я0х /а0 = а0-1Л(Йа0/2квТ), (10)
где кв - постоянная Больцмана. Здесь использовано выражение для термодинамически равновесного значения 3х в предпереходном состоянии кристалла [3]. Тогда после линеаризации Я из (6, 10) находим векторное уравнение
/& = Яхк-у(Я- %е/к), (11)
у = -Хю0Я02 и Аа0 / и.
(12)
Компоненты вектора молекулярного поля, согласно (4), имеют вид
к = (-а0,0,3 0 Я2 + /0( Я2 )2 +П (/)),
(13)
где А=?4г
1 1,т
3. Уравнения динамики системы псевдоспинов без учета релаксации
Поскольку 3у=0 (11-13), как в предпереходном состоянии, так в исходной и конечной структурах, то уравнения движения дл я компонент вектора флуктуации псевдоспина 83( без учета релаксации приобретают вид
= (X\JjS0 +Х ^я* ]у>3? + +(Х [3,1+ 2Х ^ 8 +Ц }8ЯУ, (14)
,а8Я? = (а -Я?XV» + 2X1 »Я:^}-
-(X[3Я2 +Х1 »тЯЯ2]}8Я* --Я,0Х8Ц -(X[V + 2ХIтЯт+Ц}8Я*, (15)
¡а8Я, = -а08Яу.
(16)
чений операторов флуктуации псевдоспина являются искомыми коллективными переменными:
8Я, = X8 Я, ехр(-,|Д), 8^ = XЦ ехр(-,|?). (17)
/ /
Рассмотрим свободную прецессию псевдоспина в предпереходном состоянии кристалла, достигнутом за счет изменения температуры, где 3;0г=0. Используя координаты (17), получаем для прецессионного движения с волновым вектором д систему нелинейных уравнений Блоха
,а8ЯI = + 80,-}8Я?, (18)
1а8Я? = (а0 - 3,Я0х }8Я| -
-Я Х8П, -{3?8Я, +8П4}8Я1. , а8Я, =-а08Я,.
(19)
(20)
Однородная система линеаризованных уравнений (18-20) при сЮ.д=0 имеет нетривиальное решение только в случае, если определитель этой системы тождественно равен нулю. Решая секулярное уравнение для собственных частот в духе [4], получим
а,® = 0, (21)
а23(|) =а0(а0 -3,Я0х), Я0х = (1/2)Ш(а0/2квТ). (22)
Решение а1(д )=0 (21) соответствует продольной моде, то есть флуктуации в направлении молекулярного поля, параллельного 3х в предпереходном состоянии кристалла. Два других решения а2,3( 6) (22) отвечают поперечным флуктуациям и описывают свободную прецессию псевдоспинов относительно молекулярного поля. Однако на нижней и верхней границе интервала (Т-,Т+) частота прецессии становится малой величиной вследствие ослабления двухчастичного взаимодействия псевдоспинов когда
2а0 -/2квТ) — 0, (23)
то есть при Т—^Тс. Разумно предположить, что, по крайней мере, в пределе ТС—ТМА (23), структурный переход из предпереходного состояния в исходную (конечную) фазу связан с конденсацией мягкой псевдоспиновой волны. В окрестности Тс можно разложить а23(д0) по степеням Т-Тс
2(,0) = (а/Ы/дТ)т =Т(Т- Тс) = = А(,0)(Т - Тс)/Тс,
(24)
Теперь (14-16) - это система 3N нелинейных уравнений, связывающих амплитуды флуктуаций псевдоспинов относительно молекулярного поля с амплитудой зависящей от времени компоненты механического поля напряжений. Собственные частоты системы псевдоспинов можно вычислить, полагая 0;=0 и решая однородную систему линейных уравнений (14-16) относительно а. Представленная выше система 3N уравнений может быть сведена к N системам трех уравнений для заданного волнового вектора путем преобразования к нормальным координатам. Фурье-образы средних зна-
А(,0) = (®о 3 (,0)/4 квТ2)сЪ-\а,/2 квТс), (25)
в то время как а232(д^д 0)Т=Т^0, где д 0 волновой вектор мягкой (критической) моды. В предпереход-ном состоянии при Т— Тш квадрат частоты мягкой псевдоспиновой волны, как и в случае мягких фо-нонов, является линейной функцией Т—Тс (24, 25) в окрестности точки Тс.
Исходная (конечная) структура определяется "замороженными" смещениями в левую (правую) яму двухямного потенциального рельефа, отвечающими мягкой моде. Параметром порядка, который
представляет собой замороженную координату, соответствующую мягкой нормальной моде, в данном случае служит однородная спонтанная разность заполнения (поляризация) левого и правого минимума двухямного потенциального рельефа. Положительное значение разности заполнения соответствует исходной структуре, а отрицательное значение - конечной фазе. В предпереходном состоянии разность заполнения левого и правого минимума двухямного потенциального рельефа равна нулю.
В состоянии исходной (конечной) структуры величина Я01, так же как и ^0г, отлична от нуля и молекулярное поле имеет в плоскости направление, в общем случае не совпадающее с направлением осей (11-13). Используя ту же процедуру, получим нелинейные уравнения движения, соответствующие волновому вектору д
= (У +102 f}SSy +
+{(У-+21, S02 +8П, , (26)
iкёSy = К -(У + 21-S02)S0x}8S¡ --{У 0 s0 2 +10^0 2 )2 }8sq -
-S0ёО. - (У + 21-S02 )ёSz¡ + ёО. ё/, (27) iшss; =-ШoёSy. (28)
Однородная система линеаризованных уравнений (26-28) при й0^=0 имеет нетривиальное решение только в случае, если определитель этой системы тождественно равен нулю. Решая секулярное уравнение аналогичным образом, получим в пред-переходном состоянии
К(-) = 0, (29)
а два других имеют вид
<(-) = [ У oS 0 2 +1 - ^ 02 )2]2 +
+®оК - У-S0х - 21-S02S0х). (30)
Поскольку для исходной (конечной) структуры при Т—^ Тс выполняется ^ 0х—ю0/Л, то моды псевдоспиновых волн к=0, а к,3 является критической только в длинноволновом пределе д—0 (29, 30). В этом предельном случае квадрат частоты мягкой моды определяется спонтанной разностью заполнения
К- = [JoS0z +1-)2]2-(21-/ЛК^02. (31)
При Т—Тс и Т—Тш параметр порядка ^0г стремится к нулю [3], следовательно, собственные частоты ю2,3(0) также стремятся к нулю (31). В то время как в предпереходном состоянии мода псевдоспиновых волн может стать мягкой модой при любом значении критического волнового вектора д0 (из области 0<д0<дгзБ), для которого имеет максимум, в исходной (конечной) структуре мода становится мягкой только при д —0. Если критическое значение волнового вектора в предпереходном состоянии равно д—0, тогда собственные векторы флуктуаций псевдоспинов одинаковы как в предпереходном состоянии, так и в исходной (конечной) фазе.
4. Уравнения динамики системы псевдоспинов с учетом релаксации
Перепишем для удобства векторное ур. (6) в виде системы трех уравнений для компонент псевдоспина S
Sx = [J0Sz +10(Sz )2 + Q]Sy - - S0 ], (32)
S y = -[ J о Sz + 1o(Sz )2 + Q]S - o^S - 7V , (33)
Sz = a0Sy-7{Sz + [J0Sz + 10(Sz)2 +Q]S01 a0-1}. (34)
Покажем, что из (32-34) в низкочастотном режиме вытекает уравнение типа Ландау-Халатнико-ва. Дифференцируя по времени и используя простые алгебраические преобразования, сведем (32-34) к системе двух уравнений:
Sz +у[2 + a,1 S0x ( J0 + 210 SQz )]S +
+K01S0xQ + [a02 + Y2](1 + a01S0xJ0)Sz + +K0\S0 % S0z [a2 +72](Sz )2 + +a0[J0Sz + 10(Sz)2 + Q]SSx = -a0S01 [a02 + 72]Q, (35)
SSx = -ySSx + a0-1[ J 0 Sz +10(Sz )2 +Q] x
x(£z +7{Sz + S0 xa(-1[ J 0 Sz + l0(Sz )2 + 0]}), (36)
где SS^S-S01. Разрешая систему ур. (35, 36) относительно Q как функции параметра порядка Sz и его производных методом последовательных приближений по дисперсии, диссипации и нелинейности, получим во втором порядке:
Sz +YSz +cÙ2cSz - ( 10/ J0)[a02 +y2]X
x(Sz )2 + (1/ 2)( Sz )3 = a2xn, (37)
X = X /(1 - J0Xef ) = {J0(1 + a-1J0 S01 )}-1,
К = K2 +Y2](1+C01 J0 S0x ). (38)
С другой стороны, согласно уравнению Ландау-Халатникова [5], можно записать
Sz +ySz = -vSF/SSz, (39)
где v - коэффициент, F -свободная энергия системы псевдоспинов в механическом поле, которая вблизи Tc имеет вид [3]
F = F0 + hJ0{(1 / 2)[1 - qth(ha0 / 2kBT)] x x(Sz)2 - (1/3)( 10/ J0)(5e )3 +
+(1/4)(J02/2a2)[1-n](Sz )4}-hQ.Sz, (40)
П = (1/ 2)(q - qln(q +1/q-1),
+(1/4)(J2/2a2)[1-n](Sz)4}-¿QSz. (41)
Из (37-41) видно, что статическая восприимчивость x имеет полюс при T=Tc, как того требуют общетеоретические положения. Подставляя (40) в (39) и приравнивая друг другу правые части (37) и (39) получим
J2 =®02[п/(1 -П)], * = 4q[n /(1 -П)]1/2,
v= [a^/ Й^(1 -п)], (42)
= ®о [1 -п]-1(! + <' Л 30"),5 =- /„• (43)
Таким образом, постулировав выполнение закона сохранения (5), и при наличии релаксации, нам удалось получить выражение (42) для диссипа-тивного параметра у в виде функций от ю0,Л-
Как следует из (42, 43), условие незаторможенности мягкой моды может быть записано в виде
П < [1 - ЧЩНа0/2квТ))1 (44)
Анализ (44) показывает, что данное условие выполняется лишь в узком интервале значений параметра д«1. Отсюда приходим к выводу о том, что в окрестности структурного перехода исходная структура - предпереходное состояние - конечная структура, где д>>1 (слабое, среднее туннелирова-ние [3]), мягкая мода переторможена и динамика параметра порядка становится чисто релаксационной. Общее исследование низкочастотной динамики параметра порядка в окрестности структурного перехода может быть проведено на основе ур. (37). Однако соответствующий анализ является достаточно громоздким. Случай же д>>1 представляется более простым и физически наглядным. В пределе д>>1, ур. (37), при 0=0 принимает вид уравнения Колмогорова-Петровского-Пискунова в пространственно-одномерном случае:
32 =-у[1 -9Ш( Пю012 квТ)] +
+У(10 /Л)()2 - у-\¿2 /2)(з2)3 + В52/522, (45)
где В&уа2&а0а2&10ио-1 (10-8 см=10-3 см2-е-1 - эффективный коэффициент диффузии, а - межатомное расстояние в кристаллической решетке, ухуш= =(ю0^)1/2«105 см-с-1 -скорость звука в металлах. Как отмечалось выше, в низкочастотном длинноволновом приближении слагаемое, ответственное за пространственную дисперсию, может быть учтено аддитивным образом, простым добавлением в (37).
Полученные результаты для перехода исходная фаза - предпереходное состояние - конечная фаза, вызванного изменением температуры, могут быть пролонгированы на случай структурного перехода, стимулированного изменением внешней силы. Согласно результатам [3] в окрестности структурного перехода, вынужденного изменением внешней силы, поведение параметра порядка имеет релаксационный характер, что приводит к связи параметра порядка с компонентой механического поля напряжений 0;=ЙА-1(ст)^?. Поэтому в окрестности структурного перехода, стимулированного деформацией, учет внешней силы сводится к перенормировке константы двухчастичного взаимодействия /0=/0+ЙА-1(ст), а динамика параметра порядка подчиняется ур. (45).
5. Обсуждение результатов и выводы
Одним из главных свойств рассматриваемого структурного перехода исходная структура - предпе-реходное состояние - конечная структура является то, что спонтанное появление отличного от нуля параметра порядка в исходной (конечной) фазе нарушает внутреннюю симметрию системы. Частота моды, восстанавливающей симметрию, должна непрерывно уменьшаться при Т—Т+(Т-). В случае структурного перехода почти второго рода при Т=Т+(Т-) частота этой мягкой моды практически обращается в нуль. С другой стороны, при подходе к Т+(Т-) со стороны предпереходного состояния частота мягкой моды понижающей симметрию предпереходного состояния также уменьшается при Т—Т+(Т-).
Восстанавливающая сила для смещений атомов, отвечающих такой моде, стремится к нулю до тех пор, пока псевдоспиновая волна не конденсируется на границе устойчивости. Следовательно, статические смещения атомов при почти непрерывном переходе из предпереходного состояния в исходную (конечную) структуру представляют собой замороженные смещения мягкой псевдоспиновой волны. Параметром порядка при таком переходе является статическая компонента собственного вектора мягкой псевдоспиновой волны. Так как исходная (конечная) структура характеризуется макроскопической спонтанной положительной (отрицательной) разницей заполнения левого и правого минимума одночастич-ного потенциального рельефа, одинаковой на каждом узле, то мягкая псевдоспиновая волна является длинноволновой (6 —^0). Следовательно, псевдоспиновая волна является мягкой модой, с одной стороны, понижающей симметрию предпереходного состояния, а, с другой стороны, восстанавливающей симметрию исходной (конечной) структуры.
Феноменологический учет релаксации псевдоспинов в окрестности структурного перехода, основанный на выполнении закона сохранения квадрата длины вектора Блоха, который подчеркивает коллективный, когерентный характер релаксационных процессов, приводит к следующим результатам. Во-первых, удалось подтвердить такие известные эффекты, как смягчение коллективной моды, критическое замедление. Во-вторых, из предложенного подхода естественным образом вытекает важный вывод о том, что мягкая мода в окрестности структурного перехода исходная структура -предпереходное состояние -конечная структура переторможена и, следовательно, динамика параметра порядка имеет релаксационный характер. В-третьих, поскольку релаксационные процессы приводят к связи компоненты механического поля напряжений 0.а с параметром порядка [3], динамика параметра порядка в деформируемом кристалле также имеет релаксационный характер.
СПИСОК ЛИТЕРАТУРЫ
1. Пушин В.Г., Кондратьев В.В., Хачин В.Н. Предпереходные явления и мартенситные превращения. - Екатеринбург: УрО РАН, 1998. -367 с.
2. Панин В.Е., Лихачев В..А., Гриняев Ю.В. Структурные уровни деформации твердых тел. -Новосибирск: Наука, 1985. -229 с.
3. Слядников Е.Е. Предпереходное состояние и структурный переход в деформированном кристалле // Физика твердого тела. -2004. -Т. 46. - В. 6. -С. 1065-1071.
4. Блинц Р., Жекш Б. Сегнетоэлектрики и антисегнетоэлектрики. -М.: Мир, 1975. -270 с.
5. Лифшиц Е.М., Питаевский Л.П. Физическая кинетика. — М.: Наука, 1974. —522 с.
6. Брус А., Каули Р. Структурные фазовые переходы. — М.: Мир, 1984. —408 с.
7. Барьяхтар В.Г. Интегрируемость и кинетические уравнения для солитонов. — Киев: Наукова думка, 1990. — 386 с.
8. Беленов Э.М., Назаркин А.В. Низкочастотное приближение в нелинейных уравнениях // Письма в ЖЭТФ. —1990. —Т. 51. — № 5. —С. 252—256.
УДК 533.75
ОПИСАНИЕ ПАРАМЕТРОВ СОСТОЯНИЯ ВЕЩЕСТВА В ЗАДАЧАХ МОДЕЛИРОВАНИЯ ДИНАМИКИ СИСТЕМЫ "КОНЦЕНТРИРОВАННЫЙ ПОТОК ЭНЕРГИИ - МЕТАЛЛ"
В.И. Бойко, Ю.В. Данейкин, К.В. Юшицин
Томский политехнический университет E-mail: [email protected]
Предложена модель описания параметров состояния металлов, испытывающих мощное энергетическое воздействие. Показана работоспособность модели в широком диапазоне изменения плотностей и температур. Фундаментальность физических принципов положенных в основу модели, позволила получить новые результаты при рассмотрении классических задач гидродинамики.
Развитие систем генерации сильноточных электронных пучков, мощных импульсных пучков ионов, лазерного излучения создало уникальную возможность получения новых видов воздействия концентрированных потоков энергии (КПЭ) на материалы. Процессы, сопровождающие взаимодействие КПЭ с веществом, могут использоваться для решения большого числа научных и технологических задач, в том числе: получение сверхплотной плазмы инерциального термоядерного синтеза; модификация прочностных, триботехнических, антикоррозийных свойств металлов и сплавов; синтез новых композиционных материалов, включающих метастабильные фазы и соединения; получение тонких пленок и покрытий посредством осаждения абляционной плазмы; разрушение и резка твердых материалов. Детально области применения КПЭ рассмотрены в монографии [1] и обзоре [2], а также в библиографии к ним.
Общим признаком всех областей использования КПЭ является возбуждение необходимого процесса в веществе, что реализуется определенным выбором интенсивности и времени облучения, а также видом КПЭ и подбором оптимальных амплитудно-временных параметров импульса. Современные системы генерации КПЭ позволяют получать потоки в широком диапазоне интенсивности воздействия Ж=105...1013 Вт/см2. При длительности импульса 10—8...10—6 с материал может испытывать воздействие с флюенсом 10—1.103 Дж/см2/на импульс. Такой широкий диапазон энергетического воздействия КПЭ определяет возможность генерации в объеме мишени большого числа разнородных физических явлений и, как следствие, многообразие результирующих эффектов, возникающих как при протекании отдельных процессов, так и при суперпозиции их совокупности.
Используя в качестве критерия возникновение фазовых превращений в объеме мишени при воздействии КПЭ, можно условно выделить три поддиапазона плотности мощности импульсного потока:
1. Для систем ионной имплантации и прямого нанесения на поверхность наноразмерных пленок ионным пучком используются низкоэнергетические пучки (до десятков кэВ) и соответствующими плотностями тока, обеспечивающими плотность мощности до 106 Вт/см2 [1]. Выделившейся энергии, недостаточно для развития фазовых превращений, а модифицирующий эффект обусловлен формированием дефектной структуры в приповерхностной области, не превышающей глубины проникновения ионов (десятки нм).
2. Используемые в экспериментах по инерциально-му термоядерному синтезу КПЭ предельно высоких плотностей мощности (более 1011 Вт/см2) при взаимодействии с веществом вызывают объемную сублимацию материала мишени. Рассмотрение динамики поведения мишени проводится в газовых и плазменных приближениях [2].
3. При рассмотрении КПЭ промежуточных плотностей мощности 107...1010 Вт/см2 в области энерговыделения и прилегающих к ней слоях мишени возможно формирование локальных зон различного фазового состава, в том числе с образованием слоев бинарных состояний "жидкость-твердое тело" и "жидкость-пар".
Наиболее многопараметрической в диапазоне промежуточных плотностей мощности является система "КПЭ-металл". На рис. 1 представлено характерное распределение зон различного фазового состава в приповерхностной области металлической мишени на момент окончания импульса ионного тока (10—7 с) и флюенсом около 30 Дж/см2. Несмотря на относи-
