journals.iop.org home page electronic journals * User guide   * Site map   | Quick Search:Help  
Journal of Physics: Condensed Matter
Athens/Institutional login
IOP login: Password:   
Create account | Alerts | Contact us
Journals Home | Journals List | EJs Extra | This Journal | Search | Authors | Referees | Librarians | User Options | Help |

Arrow diagram theory for non-orthogonal electronic groups: the continued fractions method

Yu Wang et al 2009 J. Phys.: Condens. Matter 21 474204 (10pp)   doi: 10.1088/0953-8984/21/47/474204  Help

   PDF (318 KB) | References

Yu Wang and Lev Kantorovich
Department of Physics, King's College London, Strand, London WC2R 2LS, UK
E-mail: lev.kantorovitch@kcl.ac.uk

Abstract. The group function theory by Tolpygo and McWeeny is a useful tool in treating quantum systems that can be represented as a set of localized electronic groups (e.g. atoms, molecules or bonds). It provides a general means of taking into account intra-correlation effects inside the groups without assuming that the interaction between the groups is weak. For non-orthogonal group functions the arrow diagram (AD) technique provides a convenient procedure for calculating matrix elements \langle \Psi |\hat {O}|\Psi \rangle
of arbitrary symmetrical operators \hat {O} which are needed, for example, for calculating the total energy of the system or its electron density. The total wavefunction of the system \Psi=\hat {A}\prod_{{I}}\Phi
_{{I}} is represented as an antisymmetrized product of non-orthogonal electron group functions ΦI of each group I in the system. However, application of the AD theory to extended (e.g. infinite) systems (such as biological molecules or crystals) is not straightforward, since the calculation of the mean value of an operator requires that each term of the diagram expansion be divided by the normalization integral S = langΨ|Ψrang which is given by an AD expansion as well. In our previous work, we cast the mean value \langle \Psi |\hat {O}|\Psi \rangle of a symmetrical operator \hat {O} in the form of an AD expansion which is a linear combination of linked (connected) ADs multiplied by numerical pre-factors. To obtain the pre-factors, a method based on power series expansion with respect to overlap was developed and tested for a simple 1D Hartree–Fock (HF) ring model. In the present paper this method is first tested on a 2D HF model, and we find that the power series expansion for the pre-factors converges extremely slowly to the exact solution. Instead, we suggest another, more powerful, method based on a continued fraction expansion of the pre-factors that approaches the exact solution much faster. The method is illustrated on the calculation of the electron density for the 2D HF model. It provides a powerful technique for treating extended systems consisting of a large number of strongly localized electronic groups.

Print publication: Issue 47 (25 November 2009)
Received 26 February 2009
Published 5 November 2009

Bookmark and Share Post to CiteUlike | Post to Connotea | Post to Bibsonomy

 

Find related articles





Article options

Authors & Referees

PhysicsWorld, subscribe nowPW launch banner
 
Content finder
  Full Search
  Help


  
Setup information is available for Adobe Acrobat.
EndNote, ProCite ® and Reference Manager ® are registered trademarks of ISI Researchsoft.
Copyright © Institute of Physics and IOP Publishing Limited 2009.
Use of this service is subject to compliance with the Terms and Conditions of use. In particular, reselling and systematic downloading of files is prohibited.
Help: Cookies | Data Protection. Privacy policy Disclaimer
 
Bioinspiration and Biomimetics reasearch banner