arxiv:hep-th/ v1 3 Jul 2000

Size: px

Start display at page:

Download "arxiv:hep-th/ v1 3 Jul 2000"

Kerry Booth
5 years ago
Views:

1 TUTORIAL ON QUASI-SPARSE EIGENVECTOR DIAGONALIZATION DEAN LEE Uiversity of Massachusetts, Amherst, MA USA arxiv:hep-th/ v1 3 Jul 2000 We review several topics related to the diagoalizatio of quatum field Hamiltoias usig the quasi-sparse eigevector (QSE method. 1 Itroductio Quasi-sparse eigevector (QSE diagoalizatio is a ew computatioal method which fids the low-lyig eigevalues ad eigevectors for a geeral quatum field Hamiltoia 1. It is able to hadle the expoetial icrease i the size of Fock space for large systems by exploitig the sparsity of the Hamiltoia. QSE diagoalizatio ca eve be applied directly to ifiitedimesioal systems. The method is most effective whe the splittig betwee low-lyig eigevalues is ot too small compared to the size of the offdiagoal Hamiltoia matrix etries. I such cases the low-lyig eigevectors are quasi-sparse, meaig that the vector is domiated by a small fractio of its largest compoets. The QSE algorithm ca the be applied to the Hamiltoia H usig the followig steps: 1. Select a subset of orthoormal basis vectors {e i1,,e i } ad call the correspodig subspace S. 2. Diagoalize H restricted to S ad fid oe eigevector v. 3. Sort the basis compoets e ij v accordig to their magitude ad remove the least importat basis vectors. 4. Replace the discarded basis vectors by ew basis vectors. These are selected at radom from a pool of cadidate basis vectors which are coected to the old basis vectors through o-vaishig matrix elemets of H. 5. Redefie S as the subspace spaed by the updated set of basis vectors ad repeat steps 2 through 5. If the subset of basis vectors is sufficietly large, the exact eigevectors will be stable fixed poits of the update process. guagzhou2000: submitted to World Scietific o November 10,

2 The purpose of this short tutorial is to provide additioal iformatio for those iterested i writig or usig codes that implemet QSE methods. We discuss the calculatio of Hamiltoia matrix elemets i the mometum Fock space represetatio ad the geeratio ad selectio of ew basis states. We also iclude a simple program which applies the QSE method to fid the groud state of a iput matrix. Readers iterested i a itroductory overview of QSE diagoalizatio are ecouraged to cosult the refereces i 1. 2 Fock states ad the Hamiltoia matrix It is most effective to describe the details of a QSE calculatio i the cotext of a specific example. I this discussio we cosider φ 4 theory i dimesios, the same example as i 1. Some of the techiques described here have bee specifically optimized for scalar field theories. Other systems require a somewhat differet set of tools ad will be discussed i the future. The Hamiltoia desity for φ 4 theory i 1+1 dimesios has the form H = 1 2 ( φ t ( φ x 2 + µ 2 2 φ2 + λ 4! :φ4 :, where the ormal orderig is with respect to the mass µ. We cosider the system i a periodic box of legth 2L. We the expad i mometum modes ad reiterpret the problem as a equivalet Schrödiger equatio 2. We will use a mometum cutoff scheme where the modes q are restricted i mometum so that N max. The resultig Hamiltoia is where H = 1 2 q q λ 4!2L =0 ω (µ = ( ω 2 q 1 q 2 q 3 q 4, λb(µ (µ 4L ad b(µ is the coefficiet for the mass couterterm b(µ = q q (1 2 π 2 L 2 +µ 2 (2 1 2ω (µ. (3 It is coveiet to split the Hamiltoia ito free ad iteractig parts with respect to a arbitrary mass µ : guagzhou2000: submitted to World Scietific o November 10,

3 H free = 1 2 H = H free q q ( + λ 4!2L =0 ω 2 (µ q q, (4 µ 2 µ 2 λb(µ 4L q 1 q 2 q 3 q 4. q q (5 µ will be used to defie the basis states of our Fock space. Sice H is idepedet µ, it is useful to perform calculatios for differet µ to obtai a estimate of the error. It is also helpful to fid the rage of values for µ which maximizes the quasi-sparsity of the eigevectors. This ca be used to improve the accuracy of the QSE calculatio. Let us defie ladder operators with respect to µ, a (µ 1 = 2ω(µ a (µ 1 = 2ω(µ These satisfy the usual commutatio relatios, [ ] q ω (µ + q (6 [ ] q ω (µ q. (7 The Hamiltoia ca ow be rewritte as H = [ a,a ] m = δm (8 [a,a m ] = [ a,a m] = 0. (9 ω (µ a a (µ2 µ 2 λb(µ 4L + λ 192L =0 [ (a1 +a 1 ω1 (µ (a 2 +a 2 ω2 (µ (a +a (a +a ω (µ (10 (a 3 +a 3 ω3 (µ ] (a 4 +a 4. ω4 (µ I (10 we have omitted costats cotributig oly to the vacuum eergy. We ca represet ay mometum-space Fock state as a strig of occupatio umbers, where o Nmax,,o Nmax, (11 a a o Nmax,,o Nmax = o o Nmax,,o Nmax. (12 guagzhou2000: submitted to World Scietific o November 10,

4 From the usual ladder operator relatios, it is straightforward to calculate the matrix elemet of H betwee two arbitrary Fock states. I 1 auxiliary cutoffs were used to reder the dimesio of Fock space fiite. This is i fact a uecessary restrictio ad QSE diagoalizatio ecouters o difficulties dealig with the full ifiite-dimesioal space. 3 New basis vectors Aside from calculatig matrix elemets, the oly other fudametal operatio ivolvig basis vectors is the geeratio of ew basis vectors. As metioed before, the ew states should be coected to a old basisvectorthrough ovaishig matrix elemets of H. Let us refer to the old basis vector as e. For this example there are two types of terms i our iteractio Hamiltoia, a quartic iteractio (a 1 +a (a 1 2 +a (a 2 3 +a (a a, , 2, 3 (13 ad a quadratic iteractio ( a +a ( a +a. (14 Toproduceaewvectorfrom e wesimplychooseoeofthepossibleoperator moomials a 1 a 2 a 3 a 1 2 3, a 1 a 2 a 3 a 1 2 3,, (15 a a, a a, ad let it act upo e. Our experiece is that the iteractios ivolvig the small mometum modes are geerally more importat tha those for the large mometum modes, a sigal that the ultraviolet divergeces have bee properly reormalized. For this reaso it is best to arrage the selectio probabilities such that the smaller values of 1, 2, 3 ad are chose more ofte. We ote that the ew basis vector selectio will occasioally fail. Either the vector is zero due to a aihilatio operator actig o a uoccupied state or the vector is already i our subset of basis vectors ad therefore ot ew. I either case we simply select a ew basis vector agai. If the selectio process fails repeatedly the a differet old basis vector e is used. guagzhou2000: submitted to World Scietific o November 10,

5 4 Sample code We have cosidered two basic operatios, the calculatio of matrix elemets ad the geeratio of ew basis states. The subrouties which perform these tasks will deped o the form of the quatum field Hamiltoia. The mai structure of the QSE algorithm, however, is idepedet of the details of the Hamiltoia. We demostrate the essetial features of the algorithm with the followig short MATLAB program. I this example we fid the groud state of a fiite matrix. To keep the example as simple as possible, we avoid calculatig matrix elemets ad geeratig ew basis states by assumig that the etire Hamiltoia matrix ca be stored i memory. iter = 30; active = 100; retai = 80; H = loadsparse( samplematrix ; vecs = [1:active]; for iter = 1:iter [v,d] = eigs(h(vecs,vecs, SR,1; weights1 = v.^2; [sortelemets, sortorder] = sort(-weights1; vecs = vecs(sortorder(1:retai; wtsum1 = cumsum(weights1/sum(weights1; while (legth(vecs <active over1 = fid(rad*wtsum1(retai <wtsum1; seedvec = vecs(over1(1; cadidates = fid(h(:,seedvec; cadidates = cadidates(cadidates ~= seedvec; weights2 = abs(h(cadidates,seedvec; wtsum2 = cumsum(weights2/sum(weights2; over2 = fid(rad <wtsum2; ewvec = cadidates(over2(1; vecs = uique([vecs ewvec]; ed ed save results.mat vecs v d iter is the umber of iteratios of the algorithm. active is the umber of basis states used to form the subspace S, ad retai is the umber of basis states kept after removig the tail of the eigevector distributio. loadsparse is a program which reads i a datafile called samplematrix. This guagzhou2000: submitted to World Scietific o November 10,

6 data is stored i the matrix variable H. The remaiig code ca be broke dow i terms of the five basic steps of the QSE algorithm. Select a subset of orthoormal basis vectors {e i1,,e i } ad call the correspodig subspace S. vecs = [1:active]; Thevariablevecsisarowvectorwhichlabelsthebasisvectorsithesubspace S. Iourexamplevecsisiitializedtobethefirstactivebasisvectors. We are assumig that the groud state of the system has a o-zero overlap with at least oe of the first active basis vectors. A more geeral iitializatio would be to choose a radom subset of basis vectors. Diagoalize H restricted to S ad fid oe eigevector v. [v,d] = eigs(h(vecs,vecs, SR,1; The fuctio eigs is a sparse matrix Aroldi diagoalizatio routie. It is beig called with the parameters SR,1 which tells eigs to fid the eigevalue ad eigevector with the smallest real part. The Hamiltoia submatrix correspodig with basis vectors i vecs is diagoalized. The eigevector is stored i the colum vector v ad the eigevalue is stored as the variable d. Sort the basis compoets e ij v accordig to their magitude ad remove the least importat basis vectors. weights1 = v.^2; [sortelemets, sortorder] = sort(-weights1; vecs = vecs(sortorder(1:retai; weights1 is a colum vector which cotais the squares of the compoets of the eigevector v. sortelemets is a list of the elemets of -weights1 from smallest to greatest. sortorder is the correspodig list of locatios for these elemets. If for example weights1 = [.2.5.3], the sortorder = [2 3 1]. The list sortelemets is ot used. Oly the retai most importat basis vectors are kept i vecs. Replace the discarded basis vectors by ew basis vectors. These are selected at radom from a pool of cadidate basis vectors which are coected to the old basis vectors through o-vaishig matrix elemets of H. wtsum1 = cumsum(weights1/sum(weights1; while (legth(vecs <active over1 = fid(rad*wtsum1(retai <wtsum1; seedvec = vecs(over1(1; cadidates = fid(h(:,seedvec; cadidates = cadidates(cadidates ~= seedvec; weights2 = abs(h(cadidates,seedvec; guagzhou2000: submitted to World Scietific o November 10,

7 wtsum2 = cumsum(weights2/sum(weights2; over2 = fid(rad <wtsum2; ewvec = cadidates(over2(1; vecs = uique([vecs ewvec]; ed cumsum is a fuctio which takes the iput list ad creates aother list cosistig of partial sums. wtsum1 is therefore a ascedig strig of umbers betwee 0 ad 1 such that the differece betwee a etry ad its precedig etry is proportioal to the correspodig value of weights1. over1 is a locatio list of elemets of wtsum1 which are greater tha the product of a geerated radom umber betwee 0 ad 1 ad wtsum1(retai. over1(1 is the first etry of over1. We ote that over1(1 must be less tha or equal to retai. seedvec is the basis vector which correspods with over1(1. cadidates is a colum vector cosistig of all basis vectors which have a o-zero matrix elemet with seedvec. Sice the task is to fid ew basis vectors, seedvec is removed from cadidates. weights2 is a colum vector cotaiig the absolute values of the matrix trasitio elemets from seedvec. wtsum2 is a ascedig strig betwee 0 ad 1 such that the differece betwee a etry ad its precedig etry is proportioal to the value of weights2. over2 is a locatio list of elemets of wtsum2 greater tha a geerated radom umber. ewvec is the basis vector which correspods with over2(1. ewvec is appeded to the set vecs. The fuctio uique sorts the iput list ad removes ay repetitios. This loop is repeated util vecs has active elemets. Redefie S as the subspace spaed by the updated set of basis vectors ad repeat steps 2 through 5. for iter = 1:iter... ed save results.mat vecs v d The algorithm is iterated iter times ad the fial results are saved i the file results.mat. More examples of QSE diagoalizatio codes for various field theory systems writte i MATLAB, Fortra, ad/or C++ will be made available i the future. 5 Summary I this short tutorial we have listed some supplemetal material for those iterested i writig or usig codes that implemet QSE methods. We have guagzhou2000: submitted to World Scietific o November 10,

8 discussed the calculatio of Hamiltoia matrix elemets ad the selectio of ew basis states. We have also preseted a simple MATLAB program which applies the QSE method to fid the groud state of a matrix. The discussio here was limited to the basic QSE method. As preseted i the Guagzhou workshop there have also bee iterestig ew developmets i QSE diagoalizatio regardig the techique of stochastic error correctio. This techique will be preseted i a forthcomig work 4. Ackowledgmets I thak my collaborators o the papers cited here ad the orgaizers ad participats of the Iteratioal Workshop o Noperturbative Methods ad Lattice QCD i Guagzhou. Fiacial support was provided by the Natioal Sciece Foudatio. Refereces 1. D. J. Lee, N. Salwe, D. D. Lee, hep-th/ ; D. Lee, hep-th/ P. Marrero, E. Roura, D. Lee, Phys. Lett. B471 ( N. Salwe, D. Lee, Phys. Rev. D62 ( D. Lee, N. Salwe, M. Widoloski, work i progress. guagzhou2000: submitted to World Scietific o November 10,

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001. Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;