CASINO manual - Theory of Condensed Matter
CASINO manual - Theory of Condensed Matter
CASINO manual - Theory of Condensed Matter
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