VMPS++
Loading...
Searching...
No Matches
Static Structure Factor for unit cells
VUMPS-Strukturfaktor mit Einheitszellen

VUMPS-Strukturfaktor mit Einheitszellen

\AtEveryCite\renewbibmacro

in: \defbibenvironmentbibliography

  1. \bibhang\bibitemsep\bibparsep
\DeclareCiteCommand

] \usebibmacroprenote \usebibmacrociteindex\printtext[bibhyperref]\usebibmacrocite \multicitedelim \usebibmacropostnote \DeclareCiteCommand*] \usebibmacroprenote \usebibmacrociteindex\printtext[bibhyperref]\usebibmacrociteyear \multicitedelim \usebibmacropostnote \DeclareCiteCommand\parencite[\mkbibparens] \usebibmacroprenote \usebibmacrociteindex\printtext[bibhyperref]\usebibmacrocite \multicitedelim \usebibmacropostnote \DeclareCiteCommand*\parencite[\mkbibparens] \usebibmacroprenote \usebibmacrociteindex\printtext[bibhyperref]\usebibmacrociteyear \multicitedelim \usebibmacropostnote \DeclareCiteCommand\footcite[\mkbibfootnote] \usebibmacroprenote \usebibmacrociteindex\printtext[bibhyperref] \usebibmacrocite \multicitedelim \usebibmacropostnote \DeclareCiteCommand\footcitetext[\mkbibfootnotetext] \usebibmacroprenote \usebibmacrociteindex\printtext[bibhyperref]\usebibmacrocite \multicitedelim \usebibmacropostnote \DeclareCiteCommand\textcite \usebibmacrociteindex\printtext[bibhyperref]\usebibmacrotextcite \multicitedelim \usebibmacrotextcite:postnote

  • m,n geht über verschiedene Einheitszellen

  • j,j geht innerhalb der Zelle von 0 bis Lc-1

  • Lc ist die Zellengröße

sαβ(q):=1||m,nj,jcelle-iq(nLc-mLc+j-j)<On,jβOm,jα>=nj,jcelle-iq(nLc+j-j)<On,jβO0,jα>=j,jcelle-iq(j-j)(ne-iqnLc<On,jβO0,jα>)=:j,jcelle-iq(j-j)sjjαβ(q) (1)

Strukturfaktor zwischen Einheitszellen:

sjjαβ(q):=ne-iqnLc<On,jβO0+jα> (2)

Zerfällt in n=0, n>0, n<0:

sjjαβ(q)=<O0+jβO0,jα>+ne-iqnLc<On,jβO0,jα>+ne+iqnLc<O-n,jβO0,jα>=<OnLc+jβO0,jα>+ne-iqnLc<O0,jαOn,jβ>+ne+iqnLc<O0,jβOn,jα>=sjj,cellαβ+sjj,interαβ(q) (3)

Es wird sjj,interαβ(q) von intercellSF (mehrere q-Werte), intercellSFpoint (ein q-Wert) berechnet. Die volle Transformation geschieht mit SF und SFpoint. Der Faktor Lc in der Exponentialfunktion hat vorher gefehlt.

Außerdem scheint die Annahme (Oα,β)=Oα,β, bzw. OαOβ=OβOα getroffen zu werden, da nur das Produkt vorkommt und man die konjugierten Operatoren nicht übergibt. In “Tangent-space methods for uniform matrix product states” (Kap. 2.5) verschwinden dann einfach die Kreuze…

Für j=j muss eigentlich nur ein Term berechnet werden:

sjj,interαβ(q)=n(e-iqnLc+e+iqnLc)<OjOj2Rene-iqnLc<OjOj> (4)

Die Gleichheit der Terme sieht man numerisch, wird aber nicht ausgenutzt.