14.09.2014 Views

CASINO manual - Theory of Condensed Matter

CASINO manual - Theory of Condensed Matter

CASINO manual - Theory of Condensed Matter

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

[<br />

where the approximation<br />

ρ (1)T<br />

α<br />

(r)ρ (1)T<br />

β<br />

(r)<br />

] R<br />

≈ ρ<br />

(1)T R<br />

Another method, proposed here, is to define a modified TBDM,<br />

˜ρ (2)<br />

αβ (r 1, r 2 ; r ′ 1, r ′ 2) = N α (N β − δ αβ )<br />

and compute the condensate fraction as<br />

α<br />

∫<br />

|Ψ(r1 , r 2 )| 2 [ Ψ(r ′ 1 ,r′ 2 )<br />

(r)ρ (1)T R<br />

β<br />

(r) has been used.<br />

Ψ(r − Ψ(r′ 1 ,r2) Ψ(r 1,r ′ 2 )<br />

1,r 2) Ψ(r 1,r 2) Ψ(r 1,r 2)<br />

∫<br />

|Ψ(R)| 2 dR<br />

]<br />

dr 3 . . . dr N<br />

, (407)<br />

c = Ω2<br />

lim<br />

N α r→∞ ˜ρ(2)T R<br />

αβ<br />

(r) , (408)<br />

which also achieves the same purposes, but benefits from correlated sampling and is somewhat cheaper<br />

to evaluate, as the wave function updates required for the OBDM can be re-used in the evaluation <strong>of</strong><br />

the TBDM.<br />

The three condensate fraction estimators have been computed for a two-dimensional electron–hole<br />

bilayer (r s = 5, d = 1, N e = N h = 58), and are represented in the figure below. From top to bottom,<br />

the TBDM estimator [Eq. (403)], the TBDM-OBDM estimator [Eq. (406)], and the modified-TBDM<br />

estimator [Eq. (408)]. The advantages <strong>of</strong> Eq. (408) are evident in the short-range region, while the<br />

long-range region seems to display a slightly noisier behaviour than the other two.<br />

0.5<br />

(Ω 2 /N α<br />

) ρ αβ<br />

(2)TR<br />

(r)<br />

Condensate fraction estimators<br />

0<br />

0.5<br />

0<br />

0.5<br />

(Ω 2 /N α<br />

) [ρ αβ<br />

(2)TR<br />

(r)-ρα<br />

(1)TR<br />

(r)ρβ<br />

(1)TR<br />

(r)]<br />

(Ω 2 /N α<br />

) ρ’ αβ<br />

(2)TR<br />

(r)<br />

0<br />

0 1 2 3 4 5 6 7<br />

r/r s<br />

33.6 Momentum density<br />

K eyword: mom den<br />

The momentum density is the Fourier transform <strong>of</strong> the one-body density matrix, and is explicitly<br />

calculated as such. The k-vectors <strong>of</strong> the transformation are the reciprocal-lattice vectors <strong>of</strong> the simulation<br />

cell. Notice that for a homogeneous system these vectors are affected by keyword k <strong>of</strong>fset,<br />

hence running different accumulations with different k <strong>of</strong>fset values allows evaluating the momentum<br />

density on an arbitrarily fine grid.<br />

33.7 Localization tensor<br />

K eyword: loc tensor<br />

An early success <strong>of</strong> quantum mechanics was to explain the distinction between metal and non-metal<br />

using band theory. The system is metallic if the conduction and valence band overlap and more than<br />

195

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!