pdf, 9 MiB - Infoscience - EPFL
pdf, 9 MiB - Infoscience - EPFL
pdf, 9 MiB - Infoscience - EPFL
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
4.3. MODEL AND METHODS 87<br />
comb lattice with N =72andN = 144 sites. The N = 72 cluster is defined<br />
by the translation vector T 1 = nb 1 + nb 2 (b i are defined in Fig. 5.2) and the<br />
vector T 2 = nb 1 + mb 2 orthogonal to T 1 . We have used anti-periodic boundary<br />
conditions along T 1 and periodic ones along T 2 . This allows to reduce the finite<br />
size effect. The N = 144 cluster uses 2T 1 and 2T 2 .<br />
4.3.2 Variational wavefunctions<br />
Next we introduce the variational subspace that we use to study the ground state<br />
of the t−J model. Our variational wavefunctions are built from the ground-states<br />
|ψ MF 〉 of the mean-field like Hamiltonian 1.23. The set of variational parameters<br />
is given by χ ij allowing anisotropic nearest neighbor hopping, when χ ij = |χ ij |e iθ ij<br />
is complex, it is breaking time-reversal symmetry leading to a so-called flux phase,<br />
∆ ij the singlet BCS pairing, and h i responsible for SDW instabilities. In general<br />
|ψ MF 〉 is not a state of fixed number of particles due to the presence of the BCS<br />
pairing. In VMC simulations it is however necessary to work with states of fixed<br />
particle number and we therefore apply the projector P N which projects on the<br />
subspace of N electrons. Further, since the t − J model allows at most one<br />
fermion per site, we discard also all configurations with doubly occupied sites by<br />
applying the complete Gutzwiller projector P G = ∏ i (1 − n i↑n i↓ ). To summarize,<br />
our variational wavefunction is written as<br />
|ψ var 〉 = P G P N |ψ MF 〉 . (4.1)<br />
Now we give a list of all the variational states which we found to be relevant in<br />
the honeycomb lattice.<br />
• FS: The Gutzwiller projected Fermi sea is our reference state, i.e. t ij ≡ 1,<br />
all other parameters equal zero, and thus no minimization is necessary.<br />
• AF: A staggered antiferromagnetic order, i.e. ∆ ij ≡ 0andh i ≡ (−1) i h.<br />
• RVB: A singlet superconducting phase 1 , i.e. h i ≡ 0, and<br />
∆ ij ≡ ∆e iφα (4.2)<br />
where the singlet order parameter has a uniform amplitude but each nearest<br />
neighbor bond α =1, 2, 3isallowedtohaveitsownphase.<br />
• RVB/AF: A state mixing superconductivity and antiferromagnetism.<br />
• F: A ferromagnetic state with partial or full polarization, i.e. ∆ ij ≡ 0and<br />
h i ≡ h.<br />
1 We have also checked the possibility of triplet pairing, however the minimum of the energy<br />
was always found for singlet pairing trial functions.