Редактирование: Переход к тепловому равновесию в гармонической ГЦК решетке
Внимание! Вы не авторизовались на сайте. Ваш IP-адрес будет публично видимым, если вы будете вносить любые правки. Если вы войдёте или создадите учётную запись, правки вместо этого будут связаны с вашим именем пользователя, а также у вас появятся другие преимущества.
Правка может быть отменена. Пожалуйста, просмотрите сравнение версий, чтобы убедиться, что это именно те изменения, которые вас интересуют, и нажмите «Записать страницу», чтобы изменения вступили в силу.
Текущая версия | Ваш текст | ||
Строка 20: | Строка 20: | ||
===Вывод уравнений=== | ===Вывод уравнений=== | ||
− | Рассмотрим кристаллическую ГЦК решетку, состоящую из одинаковых частиц массой <math> m </math>, соединенных линейными пружинками жесткостью <math> c </math>. Уравнения движения частицы с радиус-вектором <math> \textbf{r} </math> имеют следующий вид: <br /> <math> | + | Рассмотрим кристаллическую ГЦК решетку, состоящую из одинаковых частиц массой <math> m </math>, соединенных линейными пружинками жесткостью <math> c </math>. Уравнения движения частицы с радиус-вектором <math> \textbf{r} </math> имеют следующий вид: <br /> <math> \ddot{\textbf{u}}(\textbf{r}) = \sum_\alpha \textbf{C}_\alpha \textbf{u}(\textbf{r}+\textbf{a}_\alpha)</math>, <br /> |
− | где <math> \textbf{u}(\textbf{r}) = (u_x, u_y, u_z)^\top </math> - вектор-столбец, состоящий из компонент вектора перемещения частицы с радиус-вектором <math> \textbf{r} </math>, <math> \textbf{a}_\alpha </math> - векторы, соединяющие частицу с радиус-вектором <math> \textbf{r} </math> с ближайшими соседями. <math> \textbf{C}_\alpha </math> - матрицы, коэффициенты которых определяют вклад частицы номер <math> \alpha </math> в суммарную силу, действующую на частицу с радиус-вектором <math> \textbf{r} </math>. <math> \alpha = \pm 1...\pm 6 </math>, <math> \textbf{C}_\alpha = c\textbf{n}_\alpha \textbf{n}_\alpha</math>. <br /> | + | где <math> \textbf{u}(\textbf{r}) = (u_x, u_y, u_z)^\top </math> - вектор-столбец, состоящий из компонент вектора перемещения частицы с радиус-вектором <math> \textbf{r} </math>, <math> \textbf{a}_\alpha </math> - векторы, соединяющие частицу с радиус-вектором <math> \textbf{r} </math> с ближайшими соседями. <math> \textbf{C}_\alpha </math> - матрицы, коэффициенты которых определяют вклад частицы номер <math> \alpha </math> в суммарную силу, действующую на частицу с радиус-вектором <math> \textbf{r} </math>. <math> \alpha = \pm 1...\pm 6 </math>, <math> \textbf{C}_\alpha = \frac{c}{2}\textbf{n}_\alpha \textbf{n}_\alpha</math>. <br /> |
Векторы <math> \textbf{n}_\alpha = \frac{\textbf{a}_\alpha}{|\textbf{a}_\alpha|}</math> в ГЦК решетке имеют следующий вид: <br /> | Векторы <math> \textbf{n}_\alpha = \frac{\textbf{a}_\alpha}{|\textbf{a}_\alpha|}</math> в ГЦК решетке имеют следующий вид: <br /> | ||
<math> \textbf{n}_{\pm1}=\pm\frac{(\textbf{e}_x+\textbf{e}_y)}{\sqrt{2}}, \textbf{n}_{\pm4} = \pm(\textbf{n}_3-\textbf{n}_2) </math><br /> | <math> \textbf{n}_{\pm1}=\pm\frac{(\textbf{e}_x+\textbf{e}_y)}{\sqrt{2}}, \textbf{n}_{\pm4} = \pm(\textbf{n}_3-\textbf{n}_2) </math><br /> | ||
Строка 27: | Строка 27: | ||
<math> \textbf{n}_{\pm3}=\pm\frac{(\textbf{e}_x+\textbf{e}_z)}{\sqrt{2}}, \textbf{n}_{\pm6} = \pm(\textbf{n}_1-\textbf{n}_2) </math>, <br /> | <math> \textbf{n}_{\pm3}=\pm\frac{(\textbf{e}_x+\textbf{e}_z)}{\sqrt{2}}, \textbf{n}_{\pm6} = \pm(\textbf{n}_1-\textbf{n}_2) </math>, <br /> | ||
где <math> \textbf{e}_x, \textbf{e}_y, \textbf{e}_z </math> - орты декартового базиса, направленные вдоль осей кубической симметрии. <br/> | где <math> \textbf{e}_x, \textbf{e}_y, \textbf{e}_z </math> - орты декартового базиса, направленные вдоль осей кубической симметрии. <br/> | ||
− | Сделаем следующую подстановку в | + | Сделаем следующую подстановку в уравнение движения для получения дисперсионного соотношения <math> \omega </math>: <br /> |
− | <math> \textbf{u}(\textbf{r}) = e^ | + | <math> \textbf{u}(\textbf{r}) = e^{i(\omega t + \textbf{k} \cdot \textbf{r})} </math>, |
<br /> | <br /> | ||
где <math> \textbf{k} </math> - волновой вектор, и получим следующее уравнение: <br /> | где <math> \textbf{k} </math> - волновой вектор, и получим следующее уравнение: <br /> | ||
− | <math> (\textbf{D}-\omega^2 \textbf{E})=0, \textbf{D} = -\frac{1}{m} \sum_\alpha \textbf{C}_\alpha e^ | + | <math> (\textbf{D}-\omega^2 \textbf{E})=0, \textbf{D} = -\frac{1}{m} \sum_\alpha \textbf{C}_\alpha e^{\textbf{k} \cdot \textbf{a}_\alpha} </math>. <br/> |
− | + | <math> \omega^2_j </math> - собственные числа динамической матрицы <math> \textbf{D} </math>. <br/> Формула для для кинетической температуры <math> T </math>: <br/> | |
<math> T = \frac{T_0}{2} + B_1 + B_2 + B_3, \quad B_j = \int_\textbf{k} (\textrm{cos} (2\omega_j t)) \textrm{d} \textbf{k} </math>, <br /> | <math> T = \frac{T_0}{2} + B_1 + B_2 + B_3, \quad B_j = \int_\textbf{k} (\textrm{cos} (2\omega_j t)) \textrm{d} \textbf{k} </math>, <br /> | ||
− | где <math> T_0 </math> - начальное значение кинетической температуры. | + | где <math> T_0 </math> - начальное значение кинетической температуры. |
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
===Результаты=== | ===Результаты=== |