Spin bipolaron in the framework of emery model for high-T(sub c) copper oxide superconductors

The high-T(sub c) oxide compounds discovered recently exhibit a number of interesting physical properties. Two-dimensional antiferromagnetic spin order has been observed in these materials at the oxygen deficiency. This fact can be explained by strong correlation of the spins, situated on Cu sites in the conducting planes of the oxide superconductors. The doping or the oxygen deficiency lead to the occurrence of holes, occupying the oxygen p-orbitals according to the Emery model. At the small hole concentration they can move along the antiferromagnetic lattice of spins, localized on Cu sites. Researchers consider the two holes situation and describe in what way their behavior depends on the antiferromagnetic exchange interation J. It is known that in the framework of Hubbard model with strong on-site Coulomb repulsion, a single hole can form a spin polaron of the large radius. It is reasonable to admit that two holes with parallel spins (triplet) form the spin bipolaron complex owing to the hole excitations' capability to polarize Cu spin surroundings. Such an excitation was considered in the phenomenological way. Here the problem is discussed on the basis of the microscopic approach in the framework of the variational principle. A special kind of wave function is used for such a purpose. The wave function is constructed by generalizing the trial functions proposed in over two holes excitation situation (triplet) and then the region of spin bipolaron existance in the framework of Emery model is studied. In this model the Hamiltonian can be easily rewritten by forming the oxygen states transforming as the irreducible representations of the group D(sub 4)

A Fast and Efficient Algorithm for Slater Determinant Updates in Quantum Monte Carlo Simulations

We present an efficient low-rank updating algorithm for updating the trial wavefunctions used in Quantum Monte Carlo (QMC) simulations. The algorithm is based on low-rank updating of the Slater determinants. In particular, the computational complexity of the algorithm is O(kN) during the k-th step compared with traditional algorithms that require O(N^2) computations, where N is the system size. For single determinant trial wavefunctions the new algorithm is faster than the traditional O(N^2) Sherman-Morrison algorithm for up to O(N) updates. For multideterminant configuration-interaction type trial wavefunctions of M+1 determinants, the new algorithm is significantly more efficient, saving both O(MN^2) work and O(MN^2) storage. The algorithm enables more accurate and significantly more efficient QMC calculations using configuration interaction type wavefunctions