Численные методы моделирования свойств нанокристаллов. Нагорнов Ю.С. - 33 стр.

UptoLike

Составители: 

33
изменения температуры. Длительность моделирования с шагом по времени
Dt=5×10
-15
с составляла для 30000 временных шагов 1.5×10
-10
с.
Если плавление при данной температуре происходило, то весь
процесс повторяли для более низкой температуры, до тех пор, пока не
переставали фиксировать плавление системы.
Для исключения эффектов, связанных с явлением метастабильности
кристаллической фазы при плавлении, проводился дополнительный
эксперимент по долговременному (50000¸120000 временных шагов или
(2,5¸12)×10
-10
с) моделированию эволюции кристаллита при более низких
по сравнению с плавлением температурах (примерно на 100-200 К. Если
под действием поверхностных флуктуаций плотности кристаллическая
структура разрушалась, то эксперимент повторяли, снизив температуру
моделирования. Та температура, при которой кристаллит еще плавился
принималась за температуру плавления. Погрешность определения
температуры плавления составила ~100 K.
3.4. Моделирование разупорядочения кристаллов
несехиометрического диоксида урана. Радиальные функции
распределения ионов урана и кислорода
Разупорядочение в системе контролировали с помощью анализа
парных радиальных функций распределения, рассчитываемых для
различных температур и визуального наблюдения за видом атомных
плоскостей моделируемого кристаллита.
На рисунке 3.4 приведен вид смежных плоскостей для кристаллита
UO
1.65
из 6144 частиц при температурах 300K и 2300 К. Как следует из
рисунка 3.4, при комнатной температуре кристалл UO
1.65
находится в
суперионном состоянии с «расплавленной» подрешеткой кислорода, а при
2300 Кположения ионов соответствуют состоянию предшествующему
плавлению.
На рисунке 3.5 приведены зависимости РФР для ионов кислорода
относительно ионов кислорода (внутренняя зона рис. 3.5 а и
поверхность рис. 3.5 б) для температур Т=300 K и Т=2000 K для
кристаллита UO
1.65
из 6144 частиц. Вид радиальной функции подтверждает
результаты визуального наблюдения кристаллита.