Scherdehnungskorrelationen aus der Verschiebungs-Fourier-Transformation

Ich arbeite derzeit mit molekulardynamischen Simulationen und möchte Scherspannungskorrelationen in meinem 2-dimensionalen System berechnen.

Wie ich es früher gemacht habe

Kumulierte Scherdehnung an Position R zwischen den Zeiten T Und T + Δ T ist definiert als

ε X j ( R , T , T + Δ T ) = 1 2 ( X u j ( R , T , T + Δ T ) + X u j ( R , T , T + Δ T ) )
mit u ( R , T , T + Δ T ) = ( u X ( R , T , T + Δ T ) u j ( R , T , T + Δ T ) ) die Verschiebung des Partikels anfänglich an Position R zum Zeitpunkt T zwischen den Zeiten T Und T + Δ T . Daher die Autokorrelationsfunktion der Scherdehnung
C ε X j ε X j ( Δ R , Δ T ) = D T D 2 R   ε X j ( R , T , T + Δ T ) ε X j ( R + Δ R , T , T + Δ T ) D T D 2 R   ε X j ( R , T , T + Δ T ) 2
die ich berechnen möchte.

Das merkt man

D 2 R   ε X j ( R , T , T + Δ T ) ε X j ( R + Δ R , T , T + Δ T ) = F 1 { F { ε X j } × F { ε X j } } ( Δ R , T , T + Δ T )
mit F der Fourier-Transformationsoperator. Rechnerisch ist diese Identität sehr nützlich, um Korrelationen schnell auszuwerten. Bisher bin ich dann nach folgender Methode vorgegangen:

  1. Grobkörnige Scherdehnung an Positionen, die linear auf einem Gitter von den Positionen der Partikel zwischen den Zeiten verteilt sind T Und T + Δ T , im Anschluss an J. Chattoraj und A. Lemaître, Phys. Rev. Lett. 111 , 066001 (2013) ( hier erhältlich ) und Goldhirsch, I. & Goldenberg, C. Eur. Phys. J. E. (2002) 9: 245 ( hier verfügbar ).
  2. Berechnen Sie die Scherspannungskorrelationen mithilfe der schnellen Fourier-Transformation (FFT) und dann der inversen FFT aus dem erhaltenen Gitter.

Diese Methode funktioniert, ist aber leider sehr langsam, trotz meiner Bemühungen, meinen Code zu verbessern ...

Wie ich Dinge tun möchte

Es gibt in B. Illing, S. Fritschi, D. Hajnal, C. Klix, P. Keim und M. Fuchs, Phys. Rev. Lett. 117 , 208002 (2016) ( hier verfügbar mit ergänzendem Material ) eine Methode zur Berechnung von Scherdehnungskorrelationen aus der Verschiebungs-Fourier-Transformation.

Dafür führen sie – ohne viele Erklärungen – die transversale bzw. longitudinale „kollektive Mean-Square-Verschiebung“ im Fourier-Raum ein C ( Q , Δ T ) Und C | | ( Q , Δ T ) , mit Q = ( Q X Q j ) den Wellenvektor, und behaupten Sie dann, dass (siehe Gleichung 10 im Zusatzmaterial)

C ε X j ε X j ( Δ R , Δ T ) = F 1 { ( C ( Q , Δ T ) C | | ( Q , Δ T ) ) Q X 2 Q j 2 Q 2 + C ( Q , Δ T ) Q X 2 + Q j 2 4 } ( Δ R , Δ T )

Was ich nicht verstehe

Zunächst einmal habe ich mich bemüht, die Bedeutung von zu verstehen C Und C | | . Inspiriert von F. Leonforte, R. Boissière, A. Tanguy, JP Wittmer und J.-L. Barrat, Phys. Rev. B 72 , 224206 (2005) ( hier verfügbar ), habe ich die folgenden Definitionen verwendet

C ( Q , Δ T ) = 1 Q 2 | | Q F { u } ( Q , T , T + Δ T ) | | 2 C ( Q , Δ T ) = 1 Q 2 | | Q F { u } ( Q , T , T + Δ T ) | | 2
Wo bezeichnet einen zeitlichen Durchschnitt T . Die Verwendung dieser Definitionen funktioniert – fast – gut, und die Berechnung der Scherspannung ist jetzt unglaublich schneller.

Ich bin jedoch nicht in der Lage, die Berechnung durchzuführen und den Dehnungskorrelationsausdruck aus diesen Definitionen zu finden. Wenn ich keinen soliden mathematischen Beweis habe, weiß ich auch nicht, ob ich einige Faktoren vergessen habe oder ob ich mich völlig irre.

Wenn Sie diesen Beweis oder die korrekten Definitionen der kollektiven mittleren quadratischen Verschiebungen kennen C ( Q , Δ T ) Und C | | ( Q , Δ T ) , oder woanders gesehen haben, das würde mir sehr helfen! Danke schön!

Antworten (1)

Da niemand darauf geantwortet hat, werde ich es versuchen. Höchstwahrscheinlich haben Sie das inzwischen herausgefunden, aber andere könnten auf diese Frage stoßen, daher kann sie hilfreich sein.

Ich glaube, dass Ihre Definitionen von C ( Q ) Und C ( Q ) sind richtig, und ich hoffe, ich kann etwas Licht ins Dunkel bringen. Ich denke, dass der kollektive Mean-Square-Verschiebungstensor definiert ist

C = u ( Q ) u ( Q )
dh als dyadisches Produkt, a 2 × 2 Matrix (in 2D). Ich lasse die Zeitargumente aus Gründen der Klarheit durchgehend weg. Außerdem würden wir normalerweise einen Hut oder eine Tilde verwenden, um Fourier-transformierte Variablen anzuzeigen, aber ich lasse das auch weg. Nun, für Nicht-Null Q , dies ist kein isotroper Tensor, obwohl das Material (ein Glas) als isotrop angenommen wird. Es ist jedoch (durch Symmetrie) klar, dass in einem Koordinatensystem, das auf Einheitsvektoren basiert ( e , e ) , so dass definiert Q = Q e , Und e steht senkrecht dazu Q , wird der Tensor diagonal sein. Es wird eine Längskomponente von geben u , neben Q , und eine Querkomponente, senkrecht zu Q :
C ' = ( | u ( Q ) | 2 0 0 | u ( Q ) | 2 ) ( C ( Q ) 0 0 C ( Q ) )
Die ganze Physik liegt in diesen beiden Funktionen C ( Q ) Und C ( Q ) ; es gibt keinen Kreuzbegriff. Die hier gegebenen Definitionen dieser Funktionen sind (glaube ich) die gleichen wie die, die Sie dem Papier der Barrat-Gruppe entnommen haben.

Die verschiedenen Koeffizienten in dem komplizierten Dehnungskorrelationsausdruck, Gl. (10) im Zusatzmaterial zum Illing-Artikel, sind einfach das, was benötigt wird, um die Matrix von dieser diagonalen Form zurück zu drehen C ' zur raumfesten Form C zu berechnen und daraus die gewünschte Menge bezogen auf den Stamm zu berechnen C ε X j ε X j nicht nur die Verdrängung. Der Raum fixiert X j System ist natürlich willkürlich, aber festgelegt; während Sie eine Vielzahl von in Betracht ziehen Q Vektoren. Kosinus und Sinus des Rotationswinkels ϕ zwischen den beiden Koordinatensystemen werden einfach auf die Komponenten des daraus abgeleiteten Einheitsvektors bezogen Q . Die Umrechnungsformel lautet

( u X u j ) = ( cos ϕ Sünde ϕ Sünde ϕ cos ϕ ) ( u u ) = ( Q X / Q Q j / Q Q j / Q Q X / Q ) ( u u )
Die Dehnung ist der symmetrisierte Gradient der Verschiebung, daher ist der geeignete Term im Fourier-Raum (wobei ich = 1 )
ε X j ( Q ) = 1 2 [ ich Q X u j ( Q ) + ich Q j u X ( Q ) ] | ε X j | 2 = 1 4 [ Q X 2 | u j | 2 + Q j 2 | u X | 2 + Q X Q j u X u j + Q X Q j u j u X ]
Ersatz für ( u X , u j ) Die Verwendung der obigen Rotationsformel ist mühsam, aber unkompliziert, und es gibt eine gewisse Vereinfachung, da alle Kreuzbegriffe verschwinden:
C ε X j ε X j = | ε X j | 2 = Q X 2 Q j 2 Q 2 | u | 2 + Q X 4 + Q j 4 2 Q X 2 Q j 2 4 Q 2 | u | 2 = Q X 2 Q j 2 Q 2 C + Q X 4 + Q j 4 2 Q X 2 Q j 2 4 Q 2 C
Unter Berücksichtigung, dass Q X 2 + Q j 2 = Q 2 , das ist identisch mit der Formel, um die Sie sich Sorgen gemacht haben.