4.6. DLA の出力ファイル

4.6.1. フォーマット

DLA は計算結果を行区切りのプレーンテキストファイルで出力します. 行頭の文字はその行の意味を表します.

P <name> = <value>

入力パラメータファイルと格子ファイルから読み取ったパラメータ.

R <name> = <mean> <error>

計算で求められた物理量. <mean> は平均値を, <error> は標準誤差を示します.

I <text> = <value>

その他計算中に得られた情報.

C <text>

コメント.

以下にサンプル(反強磁性ハイゼンベルグ鎖)を示します.

C This is DSQSS ver.v2.0.0

I N_PROC = 1
P D       =            1
P L       =            8
P BETA    =      10.0000000000000000
P NSET    =           10
P NPRE    =         1000
P NTHERM  =         1000
P NDECOR  =         1000
P NMCS    =         1000
P SEED    =    198212240
P NSEGMAX =        10000
P NVERMAX =        10000
P BETA    = 10.000000000000
P NTAU    =           10
P NCYC    =            7
P ALGFILE = algorithm.xml
P LATFILE = lattice.xml
P WVFILE = wv.xml
P DISPFILE  = disp.xml
P OUTFILE    = sample.log
P CFOUTFILE  = cf.dat
P SFOUTFILE  = sf.dat
P CKOUTFILE  = ck.dat
P SIMULATIONTIME   =     0.000000
R sign = 1.00000000e+00 0.00000000e+00
R anv = 3.03805000e+00 1.25395375e-02
R ene = -4.55991910e-01 1.20267537e-03
R spe = -1.76672204e-02 4.09064489e-02
R som = -1.76672204e-01 4.09064489e-01
R len = 1.20014021e+01 4.78403202e-02
R xmx = 3.00035053e-01 1.19600800e-03
R amzu = -2.00000000e-04 1.08972474e-04
R bmzu = -2.00000000e-04 1.08972474e-04
R smzu = 1.32382500e-03 1.40792745e-04
R xmzu = 1.32382500e-02 1.40792745e-03
R ds1 = -1.32954309e-03 7.87178338e-04
R wi2 = 2.31040000e+01 3.83762890e-01
R rhos = 1.44400000e-01 2.39851806e-03
R rhof = inf nan
R comp = 2.43165481e+35 1.71412709e+35
R amzs0 = -2.00000000e-04 1.08972474e-04
R bmzs0 = 1.65625000e-04 1.76161818e-05
R smzs0 = 1.32382500e-03 1.40792745e-04
R xmzs0 = 1.32382500e-02 1.40792745e-03
R amzs1 = -9.25000000e-04 4.02247160e-03
R bmzs1 = 1.09209375e-01 1.12051866e-03
R smzs1 = 8.72503175e-01 8.93939492e-03
R xmzs1 = 3.00500011e+00 2.99056535e-02
R time = 1.03679300e-07 1.22794234e-09
I [the maximum number of segments]          = 123
I [the maximum number of vertices]          = 66
I [the maximum number of reg. vertex info.] = 3

以下の物理量定義に現れる記号の意味を示します.

Ns

サイト数.

Q(k)

格子点 i 上で定義される任意の演算子 Qi のフーリエ変換.

Q(k)1NsNsiQieirik

Q(τ)

虚時間 τ における演算子.

Q(τ)exp[τH]Q(τ=0)exp[τH]

˜Q

任意の演算子 Q について, 虚時間方向の平均 1ββ0dτQ(τ)

Mz

局所自由度の量子化軸方向成分. たとえばスピン系では局在スピン演算子の z 成分 Sz で, ボース粒子系では数演算子 n です.

M±

Mz の昇降演算子. スピン系では M±S± , ボース粒子系では生成消滅演算子 M+b および Mb .

Mx

非対角秩序変数. スピン系では Mx(S++S)/2 , ボース粒子系では Mx(b+b) .

T

温度.

β

逆温度.

h

Mz に共役な外場. スピン系では縦磁場, ボース粒子系では化学ポテンシャル.

Q

任意の演算子 Q のグランドカノニカル平均.

4.6.2. メイン出力

メイン出力ファイルは, 入力パラメータファイルの outfile キーワードで指定した名前で出力されます.

sign

重みの符号.

iWii|Wi|, ここで i はモンテカルロサンプルの番号.

anv

平均バーテックス数.

NvNs

ene

エネルギー密度.

ϵ1Ns(E0TNv)

spe

比熱.

CVϵT

som

比熱と温度の比.

γCVT=βCV

len

平均ワーム長さ.

xmx

横感受率.

amzu

「磁化」(uniform, τ=0).

mz1NsNsiMzi としたときの mz .

bmzu

「磁化」(uniform, τ 平均). ˜mz .

smzu

構造因子(uniform).

Szz(k=0)1Nsi,jeik(rirj)[MziMzjMziMzj]|k=0=Ns[(mz)2mz2]

xmzu

縦感受率(uniform).

χzz(k=0,ω=0)˜mzh=βNs[(˜mz)2˜mz2]

amzsK

「磁化」("staggered", τ=0)

mzK1NsNsiMzicos(kri) としたときの mzs . K は波数ベクトルXMLファイルで定義された波数 k のインデックス.

bmzuK

「磁化」("staggered", τ 平均). ˜mzK .

smzsK

構造因子 ("staggered").

Szz(k)=Ns[(mzK)2mzK2]

xmzsK

縦感受率 ("staggered").

χzz(k,ω=0)=βNs[(˜mzK)2˜mzK2]

wi2

ワインディングナンバー.

W2=Dd=1L2dW2d

rhos

超流動密度.

ρs=W22DVβ

rhof

超流動比 (superfluid fraction).

ρsmz

comp

圧縮率 (compressibility).

χzz(k=0,ω=0)˜mz2

time

モンテカルロスイープあたりの計算時間(秒).

4.6.3. 構造因子出力ファイル

構造因子出力ファイルは, 入力パラメータファイルの sfoutfile キーワードで指定した名前で出力されます. このファイルには虚時間構造因子

Szz(k,τ)Mz(k,τ)Mz(k,0)Mz(k,τ)Mz(k,0)

が出力されます. 波数 k や虚時間 τ の値は, 物理量名を用いて

R C0t0 = 1.32500000e-03 1.40929454e-04
R C0t1 = 1.32500000e-03 1.40929454e-04
R C1t0 = 7.35281032e-02 3.18028565e-04

のように C<k>t<t> という形で区別されます. ここで <k> は波数ベクトル XMLファイルの kindex (KR タグの最終要素) で指定される波数のインデックスで, <t> は離散化した虚時間のインデックス.

4.6.4. 実空間表示温度グリーン関数出力ファイル

実空間表示温度グリーン関数出力ファイルは, 入力パラメータファイルの cfoutfile キーワードで指定した名前で出力されます. このファイルには温度グリーン関数

G(rij,τ)M+i(τ)Mj

が出力されます. 変位 rij や虚時間 τ の値は構造因子と同様に, C<k>t<t> という形で物理量名によって区別されます. ここで <k> は変位XML ファイルの kind (R タグの第一要素) で指定される変位のインデックスで, <t> は離散化した虚時間のインデックス.

4.6.5. 波数表示温度グリーン関数出力ファイル

波数表示温度グリーン関数出力ファイルは, 入力パラメータファイルの ckoutfile キーワードで指定した名前で出力されます. このファイルには温度グリーン関数

G(k,τ)M+(k,τ)M(k,0)

が出力されます. 波数 k や虚時間 τ の値は構造因子と同様に, C<k>t<t> という形で物理量名によって区別されます. ここで <k> は波数ベクトルXMLファイルの kindex (RK タグの最終要素) で指定される波数のインデックスで, <t> は離散化した虚時間のインデックス.