Szを指定したTPQ計算について

[ultimatile]

いつもお世話になっております. model="Spin"でのTPQ計算の物理量が正しく計算されていないように見えるのですが確認していただけないでしょうか.

StdFace.defの例(Sz=2)

a0W=4
a0L=1
a1W=-2
a1L=4
model="Spin"
method="TPQ"
lattice="triangular Lattice"
J=1.0
2Sz=4
InitialVecType="R"
Lanczos_max=10
LargeValue=2
NumAve=1

出力例としてはFlct_rand0.datは

# inv_temp, N, N^2, D, D^2, Sz, Sz^2, step_i
0.0055537298224793 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 2.0000000000000000 0.0000000000000000 1
0.0111059427906962 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 2.0000000000000000 0.0000000000000000 2
0.0166566471027138 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 2.0000000000000000 0.0000000000000000 3
0.0222058510680096 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 2.0000000000000000 0.0000000000000000 4
0.0277535631060710 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 2.0000000000000000 0.0000000000000000 5

のようになります.

[k-yoshimi]

@ultimatile ご指摘ありがとうございます。 スピンのSz保存系ではflctの計算が自明ですので、効率化のため保存量の値を用いて計算するようにしていました。ところが関数を確認したところ、Szのみ代入され他の物理量は0が格納されたままでした。

上記修正に関しては次のバージョンで取り入れると思いますが、早急な対応が必要な場合には、上のコメントを参照にexpec_energy_flct.c関数を書き換えることで正常に稼働しますので、取り入れていただければと思います。 よろしくお願いします。

[ultimatile]

ご返信,ご対応ありがとうございます. Flctは自明なので今回ご対応いただいた物理量そのものが必要になるわけではないのですが, 万が一ですがenergy関係量も変なことになっているのではないかと不安になったためご連絡させていただきました. ありがとうございました.