Корреляции деформации сдвига по преобразованию Фурье смещения

В настоящее время работая с моделированием молекулярной динамики, я хотел бы вычислить корреляции деформации сдвига в моей двумерной системе.

Как я делал вещи

Накопленная деформация сдвига в положении р между делом т и т + Δ т определяется как

ε Икс у ( р , т , т + Δ т ) "=" 1 2 ( Икс ты у ( р , т , т + Δ т ) + Икс ты у ( р , т , т + Δ т ) )
с ты ( р , т , т + Δ т ) "=" ( ты Икс ( р , т , т + Δ т ) ты у ( р , т , т + Δ т ) ) смещение частицы первоначально в положении р вовремя т между делом т и т + Δ т . Следовательно, автокорреляционная функция деформации сдвига
С ε Икс у ε Икс у ( Δ р , Δ т ) "=" г т г 2 р   ε Икс у ( р , т , т + Δ т ) ε Икс у ( р + Δ р , т , т + Δ т ) г т г 2 р   ε Икс у ( р , т , т + Δ т ) 2
который я хочу вычислить.

Можно заметить, что

г 2 р   ε Икс у ( р , т , т + Δ т ) ε Икс у ( р + Δ р , т , т + Δ т ) "=" Ф 1 { Ф { ε Икс у } * × Ф { ε Икс у } } ( Δ р , т , т + Δ т )
с Ф оператор преобразования Фурье. Говоря вычислительным языком, это тождество очень полезно для быстрой оценки корреляций. До сих пор я следовал следующему методу:

  1. Деформация сдвига крупных зерен в позициях, линейно распределенных на сетке от положений частиц между моментами времени т и т + Δ т , вслед за J. Chattoraj и A. Lemaître, Phys. Преподобный Летт. 111 , 066001 (2013) (доступно здесь ) и Goldhirsch, I. & Goldenberg, C. Eur. физ. J. E (2002) 9: 245 (доступно здесь ).
  2. Вычислите корреляции деформации сдвига, используя быстрое преобразование Фурье (БПФ), а затем обратное БПФ из полученной сетки.

Этот метод работает, но, к сожалению, очень медленный, несмотря на все мои усилия по улучшению моего кода...

Как я хотел бы делать вещи

Имеется у B. Illing, S. Fritschi, D. Hajnal, C. Klix, P. Keim и M. Fuchs, Phys. Преподобный Летт. 117 , 208002 (2016) (доступен здесь с дополнительными материалами ) метод вычисления корреляций деформации сдвига на основе преобразования Фурье смещения.

Для этого вводят — без особых пояснений — поперечное и продольное «коллективное среднеквадратичное смещение» в пространстве Фурье соответственно. С ( д , Δ т ) и С | | ( д , Δ т ) , с д "=" ( д Икс д у ) волновой вектор, а затем утверждают, что (см. уравнение 10 в дополнительном материале)

С ε Икс у ε Икс у ( Δ р , Δ т ) "=" Ф 1 { ( С ( д , Δ т ) С | | ( д , Δ т ) ) д Икс 2 д у 2 д 2 + С ( д , Δ т ) д Икс 2 + д у 2 4 } ( Δ р , Δ т )

Что я не понимаю

Прежде всего, я изо всех сил пытался понять значение С и С | | . Вдохновленный Ф. Леонфорте, Р. Буасьером, А. Танги, Дж. П. Виттмером и Ж.-Л. Баррат, физ. Rev. B 72 , 224206 (2005) (доступно здесь ), я использовал следующие определения

С ( д , Δ т ) "=" 1 д 2 | | д Ф { ты } ( д , т , т + Δ т ) | | 2 С ( д , Δ т ) "=" 1 д 2 | | д Ф { ты } ( д , т , т + Δ т ) | | 2
где обозначает среднее значение по времени т . Использование этих определений работает — почти — отлично, а вычисление деформации сдвига теперь происходит невероятно быстро.

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

Если вы знаете это доказательство или правильные определения коллективных среднеквадратических перемещений С ( д , Δ т ) и С | | ( д , Δ т ) , или видел любой из них в другом месте, это мне очень поможет! Спасибо!

Ответы (1)

Поскольку никто не ответил на это, я попробую. Скорее всего, вы уже разобрались с этим, но другие могут столкнуться с этим вопросом, поэтому он может быть полезен.

Я считаю, что ваши определения С ( д ) и С ( д ) верны, и, надеюсь, я смогу пролить свет. Я думаю, что коллективный тензор среднеквадратичного смещения определяется

С "=" ты ( д ) * ты ( д )
т.е. как диадический продукт, 2 × 2 матрица (в 2D). Я опускаю аргумент(ы) времени для ясности. Кроме того, мы обычно используем шляпу или тильду для обозначения переменных, преобразованных Фурье, но я также опускаю это. Теперь для ненулевого д , это не изотропный тензор, хотя материал (стекло) считается изотропным. Однако ясно (по симметрии), что в системе координат, основанной на единичных векторах ( е , е ) , определенный так, что д "=" д е , и е перпендикулярно д , тензор будет диагональным. Будет продольная составляющая ты , параллельно д , и поперечная составляющая, перпендикулярная д :
С "=" ( | ты ( д ) | 2 0 0 | ты ( д ) | 2 ) ( С ( д ) 0 0 С ( д ) )
Вся физика заключается в этих двух функциях С ( д ) и С ( д ) ; перекрестного термина нет. Определения этих функций, данные здесь, такие же (я полагаю), что и те, которые вы взяли из статьи группы Баррата.

Различные коэффициенты в сложном выражении корреляции деформации, уравнение (10) в дополнительном материале к статье Иллинга, просто необходимы, чтобы повернуть матрицу назад от этой диагональной формы. С к пространственно-фиксированной форме С , и использовать его для расчета желаемой величины, связанной с деформацией С ε Икс у ε Икс у а не просто смещение. Фиксированный в пространстве Икс у система, конечно, произвольная, но фиксированная; в то время как вы будете рассматривать широкий спектр д векторы. Косинус и синус угла поворота ф между двумя системами координат просто связаны с компонентами единичного вектора, полученного из д . Формула преобразования

( ты Икс ты у ) "=" ( потому что ф грех ф грех ф потому что ф ) ( ты ты ) "=" ( д Икс / д д у / д д у / д д Икс / д ) ( ты ты )
Деформация представляет собой симметризованный градиент смещения, поэтому соответствующий термин в пространстве Фурье равен (где я "=" 1 )
ε Икс у ( д ) "=" 1 2 [ я д Икс ты у ( д ) + я д у ты Икс ( д ) ] | ε Икс у | 2 "=" 1 4 [ д Икс 2 | ты у | 2 + д у 2 | ты Икс | 2 + д Икс д у ты Икс * ты у + д Икс д у ты у * ты Икс ]
Замена для ( ты Икс , ты у ) использование приведенной выше формулы вращения утомительно, но просто, и есть некоторое упрощение, поскольку все перекрестные члены исчезают:
С ε Икс у ε Икс у "=" | ε Икс у | 2 "=" д Икс 2 д у 2 д 2 | ты | 2 + д Икс 4 + д у 4 2 д Икс 2 д у 2 4 д 2 | ты | 2 "=" д Икс 2 д у 2 д 2 С + д Икс 4 + д у 4 2 д Икс 2 д у 2 4 д 2 С
Имея в виду, что д Икс 2 + д у 2 "=" д 2 , это идентично формуле, о которой вы беспокоились.