14.09.2014 Views

pdf, 9 MiB - Infoscience - EPFL

pdf, 9 MiB - Infoscience - EPFL

pdf, 9 MiB - Infoscience - EPFL

SHOW MORE
SHOW LESS

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.

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

Saved successfully!

Ooh no, something went wrong!