В настоящее время работая с моделированием молекулярной динамики, я хотел бы вычислить корреляции деформации сдвига в моей двумерной системе.
Как я делал вещи
Накопленная деформация сдвига в положении между делом и определяется как
Можно заметить, что
Этот метод работает, но, к сожалению, очень медленный, несмотря на все мои усилия по улучшению моего кода...
Как я хотел бы делать вещи
Имеется у B. Illing, S. Fritschi, D. Hajnal, C. Klix, P. Keim и M. Fuchs, Phys. Преподобный Летт. 117 , 208002 (2016) (доступен здесь с дополнительными материалами ) метод вычисления корреляций деформации сдвига на основе преобразования Фурье смещения.
Для этого вводят — без особых пояснений — поперечное и продольное «коллективное среднеквадратичное смещение» в пространстве Фурье соответственно. и , с волновой вектор, а затем утверждают, что (см. уравнение 10 в дополнительном материале)
Что я не понимаю
Прежде всего, я изо всех сил пытался понять значение и . Вдохновленный Ф. Леонфорте, Р. Буасьером, А. Танги, Дж. П. Виттмером и Ж.-Л. Баррат, физ. Rev. B 72 , 224206 (2005) (доступно здесь ), я использовал следующие определения
Однако я не могу посчитать и найти выражение корреляции деформации из этих определений. Отсутствие твердого математического доказательства также не позволяет узнать, забыл ли я некоторые факторы или полностью ошибся.
Если вы знаете это доказательство или правильные определения коллективных среднеквадратических перемещений и , или видел любой из них в другом месте, это мне очень поможет! Спасибо!
Поскольку никто не ответил на это, я попробую. Скорее всего, вы уже разобрались с этим, но другие могут столкнуться с этим вопросом, поэтому он может быть полезен.
Я считаю, что ваши определения и верны, и, надеюсь, я смогу пролить свет. Я думаю, что коллективный тензор среднеквадратичного смещения определяется
Различные коэффициенты в сложном выражении корреляции деформации, уравнение (10) в дополнительном материале к статье Иллинга, просто необходимы, чтобы повернуть матрицу назад от этой диагональной формы. к пространственно-фиксированной форме , и использовать его для расчета желаемой величины, связанной с деформацией а не просто смещение. Фиксированный в пространстве система, конечно, произвольная, но фиксированная; в то время как вы будете рассматривать широкий спектр векторы. Косинус и синус угла поворота между двумя системами координат просто связаны с компонентами единичного вектора, полученного из . Формула преобразования