the basics of the density matrix renormalization group were developed by s. white in 1992@xcite and since then dmrg has proved to be a very powerful method for low dimensional interacting systems . its remarkable accuracy can be seen for example in the spin-1 heisenberg chain : for a system of hundreds of sites a precision of @xmath0 for the ground state energy can be achieved . since then it has been applied to a great variety of systems and problems including , among others , spin chains and ladders , fermionic and bosonic systems , disordered models , impurities and molecules and 2d electrons in high magnetic fields .
it has also been improved substantially in several directions like two ( and three ) dimensional ( 2d ) classical systems , stochastic models , the presence of phonons , quantum chemistry , field theory , the inclusion of temperature and the calculation of dynamical and time - dependent properties .
some calculations have also been performed in 2d quantum systems .
all these topics are treated in detail and in a pedagogical way in the book @xcite , where the reader can find an extensive review on dmrg . in this article
we will attempt to cover the different areas where it has been applied without entering into details but in a few cases , where we have chosen some representative contributions .
we suggest the interested reader to look for further information in the referenced work .
our aim here is to give the reader a general overview on the subject .
one of the most important limitations of numerical calculations in finite systems is the great amount of states that have to be considered and its exponential growth with system size .
several methods have been introduced in order to reduce the size of the hilbert space to be able to reach larger systems , such as monte carlo , renormalization group ( rg ) and dmrg .
each method considers a particular criterion of keeping the relevant information .
the dmrg was originally developed to overcome the problems that arise in interacting systems in 1d when standard rg procedures were applied .
consider a block b ( a block is a collection of sites ) where the hamiltonian @xmath1 and end - operators are defined .
these traditional methods consist in putting together two or more blocks ( e.g. b - b , which we will call the superblock ) , connected using end - operators , in a basis that is a direct product of the basis of each block , forming @xmath2 .
this hamiltonian is then diagonalized , the superblock is replaced by a new effective block @xmath3 formed by a certain number @xmath4 of lowest - lying eigenstates of @xmath2 and the iteration is continued ( see ref .
although it has been used successfully in certain cases , this procedure , or similar versions of it , has been applied to several interacting systems with poor performance .
for example , it has been applied to the 1d hubbard model keeping @xmath5 states . for 16 sites , an error of 5 - 10%
was obtained @xcite .
other results@xcite were also discouraging .
a better performance was obtained @xcite by adding a single site at a time rather than doubling the block size .
however , there is one case where a similar version of this method applies very well : the kondo ( and anderson ) model .
wilson@xcite mapped the one - impurity problem onto a one - dimensional lattice with exponentially descreasing hoppings .
the difference with the method explained above is that in this case , one site ( equivalent to an  onion shell " ) is added at each step and , due to the exponential decrease of the hopping , very accurate results can be obtained .
returning to the problem of putting several blocks together , the main source of error comes from the election of eigenstates of @xmath2 as representative states of a superblock .
since @xmath2 has no connection to the rest of the lattice , its eigenstates may have unwanted features ( like nodes ) at the ends of the block and this ca nt be improved by increasing the number of states kept .
based on this consideration , noack and white@xcite tried including different boundary conditions and boundary strengths .
this turned out to work well for single particle and anderson localization problems but , however , it did not improve the results significantly for interacting systems .
these considerations led to the idea of taking a larger superblock that includes the blocks @xmath6 , diagonalize the hamiltonian in this large superblock and then somehow project the most favorable states onto @xmath6 .
then @xmath6 is replaced by @xmath3 . in this way ,
awkward features in the boundary would vanish and a better representation of the states in the infinite system would be achieved .
white@xcite proposed the density matrix as the optimal way of projecting the best states onto part of the system and this will be discussed in the next section .
some considerations concerning the effect of boundary conditions in the nature of the states kept ( and an analogy to the physics of black holes ) is given in @xcite .
the justification of using the density matrix is given in detail in ref.@xcite .
a very easy and pedagogical way of understanding the basic functioning of dmrg is applying it to the calculation of simple quantum problems like one particle in a tight binding chain @xcite .    in the following section we will briefly describe the standard method ; in sect .
3 we will mention some of the most important applications ; in sect .
4 we review the most relevant extensions to the method and finally in sect .
5 we concentrate on the way dynamical calculations can be performed within dmrg .
the dmrg allows for a systematic truncation of the hilbert space by keeping the most probable states describing a wave function ( _ e.g .  _
the ground state ) instead of the lowest energy states usually kept in previous real space renormalization techniques .
the basic idea consists in starting from a small system ( _ e.g _ with @xmath7 sites ) and then gradually increase its size ( to @xmath8 , @xmath9 , ... ) until the desired length is reached .
let us call the collection of @xmath7 sites the _ universe _ and divide it into two parts : the _ system _ and the _ environment _ ( see fig .
[ figsuperblock ] ) .
the hamiltonian is constructed in the _ universe _ and its ground state @xmath10 is obtained .
this is considered as the state of the _ universe _ and called the _
target state_. it has components on the _ system _ and the _ environment_. we want to obtain the most relevant states of the _ system _ , i.e. , the states of the _ system _ that have largest weight in @xmath11 . to obtain this ,
the _ environment _ is considered as a statistical bath and the density matrix@xcite is used to obtain the desired information on the _
system_. so instead of keeping eigenstates of the hamiltonian in the block ( _ system _ ) , we keep eigenstates of the density matrix . we will be more explicit below .
= 2.0 in = 1.0 in    [ figsuperblock ]    let s define block [ * b * ] as a finite chain with @xmath12 sites , having an associated hilbert space with , @xmath4 states where operators are defined ( in particular the hamiltonian in this finite chain , @xmath1 and the operators at the ends of the block , useful for linking it to other chains or added sites ) .
except for the first iteration , the basis in this block is nt explicitly known due to previous basis rotations and reductions .
the operators in this basis are matrices and the basis states are characterized by quantum numbers ( like @xmath13 , charge or number of particles , etc ) .
we also define an added block or site as [ * a * ] having @xmath14 states .
a general iteration of the method is described below :    \i ) define the hamiltonian @xmath2 for the superblock ( the _ universe _ ) formed by putting together two blocks [ * b * ] and [ * b * ] and two added sites [ * a * ] and [ * a * ] in this way : [ * b a a b * ] ( the primes are only to indicate additional blocks , but the primed blocks have the same structure as the non - primed ones ; this can vary , see the finite - size algorithm below ) . in general , blocks
[ * b * ] and [ * b * ] come from the previous iteration .
the total hilbert space of this superblock is the direct product of the individual spaces corresponding to each block and the added sites . in practice
a quantum number of the superblock can be fixed ( in a spin chain for example one can look at the total @xmath15 subspace ) , so the total number of states in the superblock is much smaller than @xmath16 . in some cases ,
as the quantum number of the superblock consists of the sum of the quantum numbers of the individual blocks , each block must contain several subspaces ( several values of @xmath13 for example ) .
here periodic boundary conditions can be attached to the ends and a different block layout should be considered ( e.g. [ * b a b a * ] ) to avoid connecting blocks [ * b * ] and [ * b * ] which takes longer to converge .
the boundary conditions are between [ * a * ] and [ * b * ] . for closed chains
the performance is poorer than for open boundary conditions @xcite .
\ii ) diagonalize the hamiltonian @xmath2 to obtain the ground state @xmath11 ( target state ) using lanczos@xcite or davidson@xcite algorithms .
other states could also be kept , such as the first excited ones : they are all called _ target states_.
\iii ) construct the density matrix : @xmath17 on block [ * b a * ] , where @xmath18 , the states @xmath19 belonging to the hilbert space of the block [ * b a * ] and the states @xmath20 to the block [ * b a * ] .
the density matrix considers the part [ * b a * ] as a system and [ * b a * ] , as a statistical bath .
the eigenstates of @xmath21 with the highest eigenvalues correspond to the most probable states ( or equivalently the states with highest weight ) of block [ * b a * ] in the ground state of the whole superblock .
these states are kept up to a certain cutoff , keeping a total of @xmath4 states per block .
the density matrix eigenvalues sum up to unity and the truncation error , defined as the sum of the density matrix eigenvalues corresponding to discarded eigenvectors , gives a qualitative indication of the accuracy of the calculation .
\iv ) with these @xmath4 states a rectangular matrix @xmath22 is formed and it is used to change basis and reduce all operators defined in [ * b a * ] .
this block [ * b a * ] is then renamed as block [ * b@xmath23 * ] or simply [ * b * ] ( for example , the hamiltonian in block [ * b a * ] , @xmath24 , is transformed into @xmath25 as @xmath26 ) .
\v ) a new block [ * a * ] is added ( one site in our case ) and the new superblock [ * b a a b * ] is formed as the direct product of the states of all the blocks .
\vi ) this iteration continues until the desired length is achieved . at each step
the length is @xmath27 ( if [ * a * ] consists of one site ) .    when more than one target state is used , _
i.e _ more than one state is wished to be well described , the density matrix is defined as : @xmath28 where @xmath29 defines the probability of finding the system in the target state @xmath30 ( not necessarily eigenstates of the hamiltonian ) .
the method described above is usually called the _ infinite - system algorithm _ since the system size increases at each iteration .
there is a way to increase precision at each length @xmath7 called the _ finite - system algorithm_. it consists of fixing the lattice size and zipping a couple of times until convergence is reached . in this case and for the block configuration [ * b a a b * ] , @xmath31 where @xmath12 and @xmath32 are the number of sites in @xmath33 and @xmath34 respectively . in this step
the density matrix is used to project onto the left @xmath35 sites . in order to keep @xmath7
fixed , in the next block configuration , the right block @xmath34 should be defined in @xmath36 sites such that @xmath37 .
the operators in this smaller block should be kept from previous iterations ( in some cases from the iterations for the system size with @xmath38)@xcite .
the calculation of static properties like correlation functions is easily done by keeping the operators in question at each step and performing the corresponding basis change and reduction , in a similar manner as done with the hamiltonian in each block@xcite .
the energy and measurements are calculated in the superblock . in ref.@xcite
an interpretation of the correlation functions of systems at criticality is given in terms of wave function entanglement , conjecturing a modification of dmrg for these cases that preserves the entanglement .
a faster convergence of lanczos or davidson algorithm is achieved by choosing a good trial vector@xcite .
an interesting analysis on dmrg accuracy is done in ref .
fixed points of the dmrg and their relation to matrix product wave functions were studied in @xcite and an analytic formulation combining the block renormalization group with variational and fokker - planck methods in @xcite .
the connection of the method with quantum groups and conformal field theory is treated in @xcite .
there are also interesting connections between the density matrix spectra and integrable models@xcite via corner transfer matrices .
these articles give a deep insight into the essence of the dmrg method .
since its development , the number of papers using dmrg has grown enormously and other improvements to the method have been performed .
we would like to mention some applications where this method has proved to be useful .
other applications related to further developments of the dmrg will be mentioned in sect .
a very impressive result with unprecedented accuracy was obtained by white and huse @xcite when calculating the spin gap in a @xmath39 heisenberg chain obtaining @xmath40 .
they also calculated very accurate spin correlation functions and excitation energies for one and several magnon states and performed a very detailed analysis of the excitations for different momenta .
they obtained a spin correlation length of 6.03 lattice spacings .
simultaneously srensen and affleck@xcite also calculated the structure factor and spin gap for this system up to length 100 with very high accuracy , comparing their results with the nonlinear @xmath41 model .
in a subsequent paper@xcite they applied the dmrg to the anisotropic @xmath39 chain , obtaining the values for the haldane gap .
they also performed a detailed study of the @xmath42 end excitations in an open chain .
thermodynamical properties in open @xmath39 chains such as specific heat , electron paramagnetic resonance ( epr ) and magnetic susceptibility calculated using dmrg gave an excellent fit to experimental data , confirming the existence of free spins 1/2 at the boundaries@xcite .
a related problem , _
i.e. _ the effect of non - magnetic impurities in spin systems ( dimerized , ladders and 2d ) was studied in @xcite .
in addition , the study of magnon interactions and magnetization of @xmath39 chains was done in @xcite , supersymmetric spin chains modelling plateau transitions in the integer quantum hall effect in @xcite and esr studies in these systems was considered in @xcite . for larger integer spins
there have also been some studies .
nishiyama and coworkers@xcite calculated the low energy spectrum and correlation functions of the @xmath43 antiferromagnetic heisenberg open chain .
they found @xmath39 end excitations ( in agreement with the valence bond theory ) .
edge excitations for other values of @xmath44 have been studied in ref .
@xcite . almost at the same time
schollwck and jolicoeur@xcite calculated the spin gap in the same system , up to 350 sites , ( @xmath45 ) , correlation functions that showed topological order and a spin correlation length of 49 lattice spacings .
more recent accurate studies of @xmath43 chains are found in @xcite and of @xmath39 chains in staggered magnetic fields @xcite including a detalied comparison to the non - linear sigma model in @xcite . in ref .
@xcite the dispersion of the single magnon band and other properties of the @xmath43 antiferromagnetic heisenberg chains were calculated .
concerning @xmath42 systems , dmrg has been crucial for obtaining the logarithmic corrections to the @xmath46 dependence of the spin - spin correlation functions in the isotropic heisenberg model @xcite . for this , very accurate values for the energy and correlation functions were needed . for @xmath47 sites
an error of @xmath48 was achieved keeping @xmath49 states per block , comparing with the exact finite - size bethe ansatz results . for this model
it was found that the data for the correlation function has a very accurate scaling behaviour and advantage of this was taken to obtain the logarithmic corrections in the thermodynamic limit .
other calculations of the spin correlations have been performed for the isotropic @xcite and anisotropic cases @xcite .
luttinger liquid behaviour with magnetic fields have been studied in @xcite , field - induced gaps in @xcite , anisotropic systems in @xcite and the heisenberg model with a weak link in @xcite .
an analysis of quantum critical points and critical behaviour in spin chains by combining dmrg with finite - size scaling was done in @xcite .
similar calculations have been performed for the @xmath50 heisenberg chain @xcite . in this case
a stronger logarithmic correction to the spin correlation function was found .
for this model there was interest in obtaining the central charge @xmath51 to elucidate whether this model corresponds to the same universality class as the @xmath42 case , where the central charge can be obtained from the finite - size scaling of the energy .
although there have been previous attempts@xcite , these calculations presented difficulties since they involved also a term @xmath52 . with the dmrg
the value @xmath53 was clearly obtained .    in ref .
@xcite , dmrg was applied to an effective spin hamiltonian obtained from an su(4 ) spin - orbit critical state in 1d .
other applications were done to enlarged symmetry cases with su(4 ) symmetry in order to study coherence in arrays of quantum dots@xcite , to obtain the phase diagram for 1d spin orbital models@xcite and dynamical properties in a magnetic field@xcite .
dimerization and frustration have been considered in refs .
@xcite and alternating spin chains in @xcite .
the case of several coupled spin chains ( ladder models ) have been investigated in @xcite , spin ladders with cyclic four - spin exchanges in @xcite and kagome antiferromagnets in @xcite .
zigzag spin chains have been considered in @xcite and spin chains of coupled triangles in @xcite . as the dmrg s performance is optimal in open systems , an interesting analysis of the boundary effect on correlation functions is done in @xcite .
magnetization properties and plateaus for quantum spin ladder systems@xcite have also been studied .
an interesting review on the applications to some exact and analytical techniques for quantum magnetism in low dimension , including dmrg , is presented in @xcite .
there has been a great amount of applications to fermionic systems such as 1d hubbard and t - j models @xcite , luttinger liquids with boundaries@xcite , the falicov - kimball model @xcite , the quasiperiodic aubry - andre chain@xcite and fibonacci - hubbard models@xcite .
it has also been applied to field theory@xcite .
the method has been very successful for several band hubbard models@xcite , hubbard ladders @xcite and t - j ladders@xcite .
also several coupled chains at different dopings have been considered @xcite as well as flux phases in these systems @xcite .
time reversal symmetry - broken fermionic ladders have been studied in @xcite and power laws in spinless - fermion ladders in @xcite .
long - range coulomb interactions in the 1d electron gas and the formation of a wigner crystal was studied in @xcite .
several phases including the wigner crystal , incompressible and compressible liquid states , stripe and pairing phases , have been found using dmrg for 2d electrons in high magnetic fields considering different landau levels@xcite .
persistent currents in mesoscopic systems have been considered in @xcite .
quite large quasi-2d systems can be reached , for example in @xcite where a 4x20 lattice was considered to study ferromagnetism in the infinite - u hubbard model ; the ground state of a 4-leg t - j ladder in @xcite ; the one and two hole ground state in 9x9 and 10x7 t - j lattices in @xcite ; a doped 3-leg t - j ladder in @xcite ; the study of striped phases in @xcite ; domain walls in 19x8 t - j systems in @xcite ; the 2d t - j model in @xcite and the magnetic polaron in a 9x9 t - j lattice in @xcite . also big @xmath54 spin-1/2 lattices reaching 24x11 sites@xcite have been studied .
there have been some recent attempts to implement dmrg in two and higher dimensions @xcite but the performance is still poorer than in 1d .
a recent extension using a two - step dmrg algorithm for highly anisotropic spin systems has shown promising results@xcite .
impurity problems have been studied for example in one- @xcite and two - impurity @xcite kondo systems , in spin chains @xcite and in luttinger liquids @xcite .
there have also been applications to kondo and anderson lattices @xcite , kondo lattices with localized @xmath55 configurations @xcite , the two - channel kondo lattice on a ladder @xcite , a t - j chain coupled to localized kondo spins@xcite and ferromagnetic kondo models for manganites @xcite .
there have been several extensions to dmrg like the inclusion of symmetries to the method such as spin and parity@xcite .
total spin conservation and continuous symmetries have been treated in @xcite and in interaction - round a face hamiltonians@xcite , a formulation that can be applied to rotational - invariant sytems like @xmath39 and 2 chains@xcite .
a momentum representation of this technique @xcite that allows for a diagonalization in a fixed momentum subspace has been developed as well as applications in dimension higher than one@xcite and bethe lattices@xcite .
the inclusion of symmetries is essential to the method since it allows to consider a smaller number of states , enhance precision and obtain eigenstates with definite quantum numbers .
other recent applications have been in nuclear shell model calculations where a two level pairing model has been considered@xcite and in the study of ultrasmall superconducting grains , in this case , using the particle ( hole ) states around the fermi level as the system ( environment ) block@xcite .
a very interesting and successful application is a recent work in high energy physics@xcite . here
the dmrg is used in an asymptotically free model with bound states , a toy model for quantum chromodynamics , namely the two dimensional delta - function potential . for this case
an algorithm similar to the momentum space dmrg@xcite was used where the block and environment consist of low and high energy states respectively .
the results obtained here are much more accurate than with the similarity renormalization group@xcite and a generalization to field - theoretical models is proposed based on the discreet light - cone quantization in momentum space@xcite . below we briefly mention other important extensions , leaving the calculation of dynamical properties for the next section .
the dmrg has been very successfully extended to study classical systems . for a detailed description
we refer the reader to ref . @xcite .
since 1d quantum systems are related to 2d classical systems@xcite , it is natural to adapt dmrg to the classical 2d case .
this method is based on the renormalization group transformation for the transfer matrix @xmath56 ( tmrg ) .
it is a variational method that maximizes the partition function using a limited number of degrees of freedom , where the variational state is written as a product of local matrices@xcite . for 2d classical systems ,
this algorithm is superior to the classical monte carlo method in accuracy , speed and in the possibility of treating much larger systems .
a recent improvement of this method considering periodic boundary conditions is given in @xcite and a detailed comparison between symmetric and asymmetric targetting is done in @xcite .
tmrg has also been successfully used to renormalize stochastic transfer matrices in a study of cellular automatons@xcite .
the calculation of thermodynamical properties of 3d classical statistical systems has been proposed@xcite where the eigenstate of the transfer matrix with maximum eigenvalue is represented by the product of local tensors optimized using dmrg .
a further improvement to this method is based on the corner transfer matrix@xcite , the ctmrg@xcite and can be generalized to any dimension@xcite .
it was first applied to the ising model@xcite and also to the potts model@xcite , where very accurate density profiles and critical indices were calculated .
further applications have included non - hermitian problems in equilibrium and non - equilibrium physics . in the first case
, transfer matrices may be non - hermitian and several situations have been considered : a model for the quantum hall effect@xcite , the @xmath57-symmetric heisenberg chain related to the conformal series of critical models@xcite and the anisotropic triangular nearest and next - nearest neighbour ising models@xcite . in the second case , the adaptation of the dmrg to non - equilibrium physics like the asymmetric exclusion problem@xcite and reaction - diffusion problems @xcite has shown to be very successful .
it has also been applied to stochastic lattice models like in @xcite and to the 2d xy model @xcite .
the adaptation of the dmrg method for classical systems paved the way for the study of 1d quantum systems at non zero temperature , by using the trotter - suzuki method @xcite . in this case
the system is infinite and the finiteness is in the level of the trotter approximation .
standard dmrg usually produces its best results for the ground state energy and less accurate results for higher excitations .
a different situation occurs here : the lower the temperature , the less accurate the results .
very nice results have been obtained for the dimerized , @xmath42 , xy model , where the specific heat was calculated involving an extremely small basis set@xcite ( @xmath58 ) , the agreement with the exact solution being much better in the case where the system has a substantial gap .
it has also been used to calculate thermodynamical properties of the anisotropic @xmath42 heisenberg model , with relative errors for the spin susceptibility of less than @xmath59 down to temperatures of the order of @xmath60 keeping @xmath61 states@xcite
. a complete study of thermodynamical properties like magnetization , susceptibility , specific heat and temperature dependent correlation functions for the @xmath42 and 3/2 heisenberg models was done in @xcite .
other applications have been the calculation of the temperature dependence of the charge and spin gap in the kondo insulator@xcite , the calculation of thermodynamical properties of ferrimagnetic chains@xcite and spin ladders@xcite , the study of impurity properties in spin chains@xcite , frustrated quantum spin chains@xcite , t - j@xcite and spin ladders@xcite and dimerized frustrated heisenberg chains@xcite .
an alternative way of incorporating temperature into the dmrg procedure was developed by moukouri and caron@xcite .
they considered the standard dmrg taking into account several low - lying target states ( see eq .
[ eq : pl ] ) to construct the density matrix , weighted with the boltzmann factor ( @xmath62 is the inverse temperature ) : @xmath63 with this method they performed reliable calculations of the magnetic susceptibility of quantum spin chains with @xmath42 and @xmath64 , showing excellent agreement with bethe ansatz exact results .
they also calculated low temperature thermodynamical properties of the 1d kondo lattice model@xcite and of organic conductors @xcite .
zhang et al.@xcite applied the same method in the study of a magnetic impurity embedded in a quantum spin chain .
a significant limitation to the dmrg method is that it requires a finite basis and calculations in problems with infinite degrees of freedom per site require a large truncation of the basis states@xcite .
however , jeckelmann and white developed a way of including phonons in dmrg calculations by transforming each boson site into several artificial interacting two - state pseudo - sites and then applying dmrg to this interacting system@xcite ( called the  pseudo - site system " ) .
the idea is based on the fact that dmrg is much better able to handle several few - states sites than few many - state sites@xcite .
the key idea is to substitute each boson site with @xmath65 states into @xmath7 pseudo - sites with 2 states@xcite .
they applied this method to the holstein model for several hundred sites ( keeping more than a hundred states per phonon mode ) obtaining negligible error .
in addition , up to date , this method is the most accurate one to determine the ground state energy of the polaron problem ( holstein model with a single electron ) .
an alternative method ( the  optimal phonon basis")@xcite is a procedure for generating a controlled truncation of a large hilbert space , which allows the use of a very small optimal basis without significant loss of accuracy .
the system here consists of only one site and the environment has several sites , both having electronic and phononic degrees of freedom .
the density matrix is used to trace out the degrees of freedom of the environment and extract the most relevant states of the site in question . in following steps , more bare phonons
are included to the optimal basis obtained in this way .
this method was successfully applied to study the interactions induced by quantum fluctuations in quantum strings , as an application to cuprate stripes@xcite and the dissipative two - state system @xcite .
a variant of this scheme is the  four block method " , as described in @xcite .
they obtain very accurately the luttinger liquid - cdw insulator transition in the 1d holstein model for spinless fermions .
the method has also been applied to pure bosonic systems such as the disordered bosonic hubbard model@xcite , where gaps , correlation functions and superfluid density are obtained .
the phase diagram for the non - disordered bose - hubbard model , showing a reentrance of the superfluid phase into the insulating phase was calculated in ref .
it has also been used to study a chain of oscillators with optical phonon spectrum @xcite and optical phonons in the quarter - filled hubbard model for organic conductors@xcite .
the dmrg has also been generalized to 1d random and disordered systems , and applied to the random antiferromagnetic and ferromagnetic heisenberg chains@xcite , including quasiperiodic exchange modulation@xcite and a detailed study of the haldane@xcite and griffiths phase@xcite in these systems .
strongly disordered spin ladders have been considered in @xcite .
it has also been used in disordered fermi systems such as the spinless model@xcite . in particular , the transition from the fermi glass to the mott insulator and the strong enhancement of persistent currents in the transition was studied in correlated one - dimensional disordered rings@xcite .
disorder - induced crossover effects at quantum critical points were studied in @xcite
.      there have been several applications to molecules and polymers , such as the pariser - parr - pople ( ppp ) hamiltonian for a cyclic polyene@xcite ( where long - range interactions are included ) , magnetic keplerate molecules @xcite , molecular iron rings@xcite and polyacenes considering long range interactions @xcite .
it has also been applied to conjugated organic systems ( polymers ) , adapting the dmrg to take into account the most important symmetries in order to obtain the desired excited states@xcite . also conjugated one - dimensional semiconductors @xcite have been studied , in which the standard approach can be extended to complex 1d oligomers where the fundamental repeat is not just one or two atoms , but a complex molecular building block .
relatively new fields of application are the calculation of dynamical properties in the rubinstein - duke model for reptons @xcite and excitons in dendrimer molecules@xcite .
recent attempts to apply dmrg to the _ ab initio _ calculation of electronic states in molecules have been successful@xcite . here , dmrg is applied within the conventional quantum chemical framework of a finite basis set with non - orthogonal basis functions centered on each atom .
after the standard hartree - fock ( hf ) calculation in which a hamiltonian is produced within the orthogonal hf basis , dmrg is used to include correlations beyond hf , where each orbital is treated as a  site " in a 1d lattice .
one important difference with standard dmrg is that , as the interactions are long - ranged , several operators must be kept , making the calculation somewhat cumbersome .
however , very accurate results have been obtained in a check performed in a water molecule ( keeping up to 25 orbitals and @xmath66 states per block ) , obtaining an offset of 0.00024hartrees with respect to the exact ground state energy@xcite , a better performance than any other approximate method@xcite .    in order to avoid the non - locality introduced in the treatment explained above , white introduced the concept of _ orthlets _ , local , orthogonal and compact wave functions that allow prior knowledge about singularities to be incorporated into the basis and an adequate resolution for the cores@xcite .
the most relevant functions in this basis are chosen via the density matrix .
an application based on the combination with the momentum version of dmrg is used in @xcite to calculate the ground state of several molecules .
the dmrg was originally developed to calculate static ground state properties and low - lying energies .
however , it can also be useful to calculate dynamical response functions .
these are of great interest in condensed matter physics in connection with experiments such as nuclear magnetic resonance ( nmr ) , neutron scattering , optical absorption , photoemission , etc .
we will describe three different methods in this section .
a recent development for calculating response functions in single impurity systems in the presence of a magnetic field was done in @xcite by using the dmrg within wilson s nrg to obtain the green s function .
an interesting extension of dmrg to tackle time - dependent quantum many - body systems out of equilibrium was considered by cazalilla and marston@xcite .
an effective way of extending the basic ideas of this method to the calculation of dynamical quantities is described in ref.@xcite .
it is important to notice here that due to the particular real - space construction , it is not possible to fix the momentum as a quantum number .
however , we will show that by keeping the appropriate target states , a good value of momentum can be obtained .
we want to calculate the following dynamical correlation function at @xmath67 : @xmath68 where @xmath69 is the hermitean conjugate of the operator @xmath70 , @xmath71 is the heisenberg representation of @xmath70 , and @xmath72 is the ground state of the system .
its fourier transform is : @xmath73 where the summation is taken over all the eigenstates @xmath74 of the hamiltonian @xmath75 with energy @xmath76 , and @xmath77 is the ground state energy .    defining the green s function @xmath78
the correlation function @xmath79 can be obtained as @xmath80    the function @xmath81 can be written in the form of a continued fraction : @xmath82 the coefficients @xmath83 and @xmath84 can be obtained using the following recursion equations @xcite : @xmath85 where @xmath86    for finite systems the green s function @xmath87 has a finite number of poles so only a certain number of coefficients @xmath83 and @xmath84 have to be calculated .
the dmrg technique presents a good framework to calculate such quantities . with it , the ground state , hamiltonian and the operator @xmath70 required for the evaluation of @xmath79
are obtained .
an important requirement is that the reduced hilbert space should also describe with great precision the relevant excited states @xmath88 .
this is achieved by choosing the appropriate target states . for most systems
it is enough to consider as target states the ground state @xmath11 and the first few @xmath89 with @xmath90 and @xmath91 as described above . in doing so
, states in the reduced hilbert space relevant to the excited states connected to the ground state via the operator of interest @xmath70 are included .
the fact that @xmath92 is an excellent trial state , in particular , for the lowest triplet excitations of the two - dimensional antiferromagnet was shown in ref .
of course , if the number @xmath4 of states kept per block is fixed , the more target states considered , the less precisely each one of them is described . an optimal number of target states and @xmath4
have to be found for each case . due to this reduction ,
the algorithm can be applied up to certain lengths , depending on the states involved . for longer chains , the higher
energy excitations will become inaccurate .
proper sum rules have to be calculated to determine the errors in each case .    as an application of the method
we calculate @xmath93 for the 1d isotropic heisenberg model with spin @xmath42 .
the spin dynamics of this model has been extensively studied .
the lowest excited states in the thermodynamic limit are the des  cloiseaux - pearson triplets @xcite , having total spin @xmath94 .
the dispersion of this spin - wave branch is @xmath95 . above this lower boundary
there exists a two - parameter continuum of excited triplet states that have been calculated using the bethe ansatz approach @xcite with an upper boundary given by @xmath96 .
it has been shown @xcite , however , that there are excitations above this upper boundary due to higher order scattering processes , with a weight that is at least one order of magnitude lower than the spin - wave continuum .    in fig .
2 we show the spectrum for @xmath97 and @xmath98 for different values of @xmath4 , where exact results are available for comparison .
the delta peaks of eq .
( [ eq : szzn ] ) are broadened by a lorentzian for visualizing purposes . as expected
, increasing @xmath4 gives more precise results for the higher excitations .
this spectra has been obtained using the infinite - system method and more precise results are expected using the finite - system method , as described later .
= 3.3 in = 2.75 in    in fig .
3 we show the spectrum for two systems lengths and @xmath97 and @xmath99 keeping @xmath100 states and periodic boundary conditions . for this case
it was enough to take 3 target states , _
i. e.  _ @xmath11 , @xmath101 and @xmath102 .
here we have used @xmath103 pairs of coefficients @xmath83 and @xmath84 , but we noticed that if we considered only the first ( @xmath104 ) coefficients @xmath83 and @xmath84 , the spectrum at low energies remains essentially unchanged .
minor differences arise at @xmath105 .
this is another indication that only the first @xmath106 are relevant for the low energy dynamical properties for finite systems .    in the inset of fig .
3 the spectrum for @xmath99 and @xmath107 is shown .
for this case we considered 5 target states _
i. e.  _ @xmath11 , @xmath108 , @xmath109 and @xmath100 . here ,
and for all the cases considered , we have verified that the results are very weakly dependent on the weights @xmath29 of the target states ( see eq.([eq : pl ] ) ) as long as the appropriate target states are chosen . for lengths where
this value of @xmath57 is not defined we took the nearest value .
= 3.3 in = 2.75 in    even though we are including states with a given momentum as target states , due to the particular real - space construction of the reduced hilbert space , this translational symmetry is not fulfilled and the momentum is not fixed .
to check how the reduction on the hilbert space influences the momentum @xmath57 of the target state @xmath110 , we calculated the expectation values @xmath111 for all @xmath112 .
if the momenta of the states were well defined , this value is proportional to @xmath113 if @xmath114 . for @xmath115 , @xmath116 .    the momentum distribution for @xmath97 is shown in fig . 4 in a semilogarithmic scale where
the @xmath117-axis has been shifted by .003 so as to have well - defined logarithms .
we can see here that the momentum is better defined , even for much larger systems , but , as expected , more weight on other @xmath112 values arises for larger @xmath7 .    as a check of the approximation we calculated the sum rule @xmath118 for @xmath107 , 5 target states and @xmath100 .
we obtain a relative error of 0.86% .    recently , important improvements to this method have been published @xcite : by considering the finite system method in open chains , khner and white obtained a higher precision in dynamical responses of spin chains . in order to define a momentum in an open chain and to avoid end effects ,
they introduce a filter function with weight centered in the middle of the chain and zero at the boundaries .
recent applications of this method include the calculation of excitations in spin - orbital models ( su(4 ) ) in a magnetic field@xcite , spin dynamics in models for cuprate spin ladders including cyclic spin exchange@xcite , optical conductivity of the ionic hubbard model@xcite , excitations in the one - dimensional bose - hubbard model@xcite and the optical response in 1d mott insulators@xcite .    in this section
we have presented a method of calculating dynamical responses with dmrg .
although the basis truncation is big , this method keeps only the most relevant states and , for example , even by considering a @xmath119 of the total hilbert space ( for @xmath107 only @xmath120 40000 states are kept ) a reasonable description of the low energy excitations is obtained .
we show that it is also possible to obtain states with well defined momenta if the appropriate target states are used .      introduced in ref .
@xcite in the dmrg context and improved in ref .
@xcite , this method focuses on a particular energy or energy window , allowing a more precise description in that range and the possibility of calculating spectra for higher energies . instead of using the tridiagonalization of the hamiltonian , but in a similar spirit regarding the important target states to be kept , the spectrum can be calculated for a given @xmath121 by using a correction vector ( related to the operator @xmath70 that can depend on momentum @xmath57 ) .
following ( [ eq : din ] ) , the ( complex ) correction vector @xmath122 can be defined as : @xmath123 so the green s function can be calculated as @xmath124    separating the correction vector in real and imaginary parts @xmath125 we obtain @xmath126 and @xmath127 the former equation is solved using the conjugate gradient method . in order to keep the information of the excitations at this particular energy the following states are targeted in the dmrg iterations : the ground state @xmath128 , the first lanczos vector @xmath129 and the correction vector @xmath122 . even though only a certain energy is focused on
, dmrg gives the correct excitations for an energy range surrounding this particular point so that by running several times for nearby frequencies , an approximate spectrum can be obtained for a wider region @xcite .
a variational formulation of the correction vector technique which leads to more accurate excited energies and spectral weights has been developed in @xcite .
it has been successfully applied to calculate the optical conductivity of mott insulators@xcite .
this method@xcite relies on a moment expansion of the dynamical correlations using sum rules that depend only on static correlation functions which can be calculated with dmrg . with these moments
, the green s functions can be calculated using the maximum entropy method .
the first step is the calculation of sum rules . as an example , and following @xcite , the spin - spin correlation function @xmath130 of the heisenberg model is calculated where the operator @xmath70 of eq .
( [ eq : ca ] ) is @xmath131 and the sum rules are@xcite : @xmath132,s^z(-q ) ] \rangle \nonumber \\ & = & 2[1-\cos(q)]\sum_i\langle s^+_{i}s^-_{i+1}+s^-_{i}s^+_{i+1}\rangle\end{aligned}\ ] ] where @xmath133 is the static susceptibility .
these sum rules can be easily generalized to higher moments : @xmath134 ... ],s^z(-q)]\rangle \end{aligned}\ ] ] for @xmath12 odd .
a similar expression is obtained for @xmath12 even , where the outer square bracket is replaced by an anticommutator and the total sign is changed . here
@xmath75 appears in the commutator @xmath135 times .
apart from the first moment which is given by the static susceptibility , all the other moments can be expressed as equal time correlations ( using a symbolic manipulator ) .
the static susceptibility @xmath136 is calculated by applying a small field @xmath137 and calculating the density response @xmath138 with dmrg .
then @xmath139 for @xmath140 .
these moments are calculated for several chain lengths and extrapolated to the infinite system .
once the moments are calculated , the final spectra is constructed via the maximum entropy method ( me ) , which has become a standard way to extract maximum information from incomplete data ( for details see ref . @xcite and references therein ) .
reasonable spectra are obtained for the xy and isotropic models , although information about the exact position of the gaps has to be included .
otherwise , the spectra are only qualitatively correct .
this method requires the calculation of a large amount of moments in order to get good results : the more information given to the me equations , the better the result .      in order to include temperature in the calculation of dynamical quantities , the transfer matrix rg described above ( tmrg@xcite )
was extended to obtain imaginary time correlation functions@xcite .
after fourier transformation in the imaginary time axis , analytic continuation from imaginary to real frequencies is done using maximum entropy ( me ) .
the combination of the tmrg and me is free from statistical errors and the negative sign problem of monte carlo methods .
since we are dealing with the transfer matrix , the thermodynamic limit can be discussed directly without extrapolations . however , in the present scheme , only local quantities can be calculated .
a systematic investigation of local spectral functions is done in ref .
@xcite for the anisotropic heisenberg antiferromagnetic chain .
the authors obtain good qualitative results especially for high temperatures but a quantitative description of peaks and gaps are beyond the method , due to the severe intrinsic limitation of the analytic continuation .
this method was also applied with great success to the 1d kondo insulator@xcite .
the temperature dependence of the local density of states and local dynamic spin and charge correlation functions were calculated .
we have presented here a very brief description of the density matrix renormalization group technique , its applications and extensions .
the aim of this article is to give the unexperienced reader an idea of the possibilities and scope of this powerful , though relatively simple method .
the experienced reader can find here an extensive ( however incomplete ) list of references covering most applications to date using dmrg in a great variety of fields such as condensed matter , statistical mechanics and high energy physics .
the author acknowledges hospitality at the centre de recherches mathematiques , university of montreal and at the physics department of the university of buenos aires , argentina , where this work has been performed .
we thank s. white for a critical reading of the manuscript and all those authors that have updated references and sent instructive comments .
k. h. is a fellow of conicet , argentina .
grants : pict 03 - 00121 - 02153 and pict 03 - 00000 - 00651 .
999 s. white , phys .
lett . * 69 * , 2863 ( 1992 ) , edited by i. peschel , x. wang , m. kaulke and k. hallberg ( series : lecture notes in physics , springer , berlin , 1999 ) s. white , phys .
b * 48 * , 10345 ( 1993 ) j. bray and s. chui , phys . rev .
b * 19 * , 4876 ( 1979 ) c. pan and x. chen , phys .
b * 36 * , 8600 ( 1987 ) ; m. kovarik , phys . rev .
b * 41 * , 6889 ( 1990 ) t. xiang and g. gehring , phys .
b * 48 * , 303 ( 1993 ) k. wilson , _ rev . mod . phys . _
* 47 * , 773 ( 1975 ) s. white and r. noack , phys .
. lett . * 68 * , 3487 , ( 1992 ) ; r. noack and s. white , phys .
b * 47 * , 9243 ( 1993 ) j. gaite , mod .
a*16 * , 1109 ( 2001 ) r. noack and s. white in ref .
@xcite , chap .
2(i ) m. martn - delgado , g. sierra and r. noack , j. phys .
a : math . gen . * 32 * , 6079 ( 1999 ) r. feynman , _ statistical mechanics : a set of lectures _ , ( benjamin , reading , ma , 1972 ) s. r. white , i. affleck and d. scalapino , phys .
b * 65 * , 165122 ( 2002 ) see e. dagotto , rev .
mod . phys . * 66 * , 763 ( 1994 ) e. r. davidson , j. comput . phys . * 17 * , 87 ( 1975 ) ; e.r . davidson , _ computers in physics _ * 7 * , no . 5 , 519 ( 1993 ) .
t. osborne and m. nielsen , quantum information processing , volume 1 , issue 1 - 2 , 45 ( 2002 ) ( available on - line : http://www.kluweronline.com/issn/1570-0755/ ) s. white , phys .
lett . * 77 * , 3633 ( 1996 ) t. nishino and k. okunishi , phys .
64 * , 4084 ( 1995 ) ; u. schollwck , phys . rev .
b * 58 * , 8194 ( 1998 ) and phys . rev .
b * 59 * , 3917 ( 1999 ) ( erratum ) .
. legeza and g. fth , phys .
b * 53 * , 14349 ( 1996 ) ; m - b lepetit and g. pastor , phys . rev .
b * 58 * , 12691 ( 1998 ) s. stlund and s. rommer , phys . rev .
lett . * 75 * , 3537 ( 1995 ) ; phys .
b * 55 * , 2164 ( 1997 ) ; m. andersson , m. boman and s. stlund , phys .
b * 59 * , 10493 ( 1999 ) ; h. takasaki , t. hikihara and t. nishino , j. phys . soc
* 68 * , 1537 ( 1999 ) ; k. okunishi , y. hieida and y.akutsu , phys .
e * 59 * , r6227 ( 1999 ) m. a. martn - delgado and g. sierra , int
a , * 11 * , 3145 ( 1996 ) g. sierra and m. a. martn - delgado , in `` the exact renormalization group '' by a krasnitz , r potting , y a kubyshin and p. s. de sa ( eds . ) world scientific pub co ; isbn : 9810239394 , ( 1999 ) , ( cond - mat/9811170 ) i. peschel , m. kaulke and . legeza , annalen der physik * 8 * , 153 ( 1999),(cond - mat/9810174 ) s. r. white and d. huse , phys . rev .
b * 48 * , 3844 ( 1993 ) e. s. srensen and i. affleck , phys .
b * 49 * , 13235 ( 1994 ) e. s. srensen and i. affleck , phys . rev .
b * 49 * , 15771 ( 1994 ) ; e. polizzi , f. mila and e. srensen , phys . rev .
b * 58 * , 2407 ( 1998 ) c. batista , k. hallberg and a. aligia , phys . rev .
b * 58 * , 9248 ( 1998 ) ; phys . rev .
b * 60 * , 12553 ( 1999 ) ; physica b , * 259 * , 1017 ( 1999 ) ; e. jannod , c. payen , k. schoumacker , c. batista , k. hallberg and a. aligia , phys .
rev b , * 62 * , 2998 ( 2000 ) m. laukamp et al .
b * 5 * , 10755 ( 1998 ) t - k .
ng , j. lou and z. su , phys .
b * 61 * , 11487 ( 2000 ) j. lou , s. qin , t - k .
ng , z. su and i. affleck , phys .
b * 62 * , 3786 ( 2000 ) s - w tsai and j. b. marston , ann . phys .
( leipzig ) 8 , special issue , 261 ( 1999 ) m. sieling , u. lw , b. wolf , s. schmidt , s. zvyagin and b.lthi , phys .
b * 61 * , 88 ( 2000 ) y. nishiyama , k. totsuka , n. hatano and m. suzuki , j. phys . soc .
jpn * 64 * , 414 ( 1995 ) s. qin , t. ng and z - b su , phys .
b * 52 * , 12844 ( 1995 ) u. schollwck and t. jolicoeur , europhys
. lett . * 30 * , 493 ( 1995 ) x. wang , s. qin and lu yu , phys .
b. * 60 * , 14529 ( 1999 ) w. tatsuaki , phys .
e * 61 * , 3199 ( 2000 ) ; t. wada and t. nishino , cond - mat/0103508 ( proceedings of the conference on computational physics 2000 ( ccp2000 ) , gold coast , queensland , australia , 3 - 8 december 2000 ) m. capone and s. caprara , phys .
b * 64 * , 184418 ( 2001 ) j. lou , x. dai , s. qin , z. su and l. yu , phys . rev .
b * 60 * , 52 ( 1999 ) e. ercolessi , g. morandi , p. pieri and m. roncaglia , europhys .
lett . * 49 * , 434 ( 2000 ) s. qin , x. wang and lu yu , phys .
b * 56 * , r14251 ( 1997 ) k. hallberg , p. horsch and g. martnez , phys . rev .
b , * 52 * , r719 ( 1995 ) h. e. boos , v. e. korepin , y. nishiyama and m. shiroishi , j. phys .
a*35 * , 4443 ( 2002 ) ; v. e. korepin , s. lukyanov , y. nishiyama and m. shiroishi , cond - mat/0210140 m. shiroishi , m. takahashi and y. nishiyama , j. phys .
* 70 * , 3535 ( 2001 ) t. hikihara and a. furusaki , phys .
b * 58 * , r583 ( 1998 ) g. fath , cond - mat/0208580 j. lou , s. qin , c. chen , z. su and l. yu , phys . rev .
b * 65 * , 064420 ( 2002 ) f. capraro and c. gros , cond - mat/0207279 y. hieida , k. okunishi and y. akutsu , phys . rev .
b * 64 * , 224422 ( 2001 ) t. byrnes , r. bursill , h - p .
eckle , c. hamer and a. sandvik , cond - mat/0205140 s - w . tsai and j. b. marston , phys .
b * 62 * , 5546 ( 2000 ) k. hallberg , x. wang , p. horsch and a. moreo , phys . rev . lett . * 76 * , 4955 ( 1996 ) a. moreo , phys . rev .
b * 35 * , 8562 ( 1987 ) ; t. ziman and h. schulz , phys .
lett * 59 * , 140 ( 1987 ) y. yamashita , n. shibata and k. ueda , phys .
b * 58 * , 9114 ( 1998 ) ; j. phys .
vol.69 , 242 - 247 ( 2000 ) ; phys .
b * 61 * , 4012 ( 2000 ) a. onufriev and b. marston , phys . rev .
b * 59 * , 12573 ( 1999 ) c. itoi , s. qin and i. affleck , phys . rev .
b * 61 * , 6747 ( 2000 ) w. yu and s. haas , phys . rev .
b * 63 * , 024423 ( 2001 ) r.  j.  bursill , t.  xiang and g.  a.  gehring , j.phys .
a * 28 * 2109 ( 1994 ) r.  j.  bursill , g.  a.  gehring , d.  j.  j.  farnell , j.  b.  parkinson , t. xiang and c. zeng , j.  phys .
c * 7 * 8605 ( 1995 ) u.  schollwck , th .
jolicoeur and t.  garel , phys .
b * 53 * , 3304 ( 1996 ) r.  chitra , s.  pati , h.  r.  krishnamurthy , d.  sen and s.  ramasesha , phys .
b * 52 * , 6581 ( 1995 ) ; s.  pati , r.  chitra , d.  sen , h.  r.  krishnamurthy and s.  ramasesha , europhys .
* 33 * , 707 ( 1996 ) ; j. malek , s. drechsler , g. paasch and k. hallberg , phys .
b * 56 * , r8467 ( 1997 ) ; e. srensen et al in ref.@xcite , chap.1.2 ( part ii ) and references therein ; d. augier , e. srensen , j. riera and d. poilblanc , phys . rev .
b * 60 * , 1075 ( 1999 ) y.  kato and a.  tanaka , j.  phys .
jpn .  * 63 * , 1277 ( 1994 ) s.  r.  white and i.  affleck , phys .
b * 54 * , 9862 ( 1996 ) g. bouzerar , a. kampf and g. japaridze , phys . rev .
b * 58 * , 3117 ( 1998 ) ; g. bouzerar , a. kampf and f. schnfeld , cond - mat/9701176 , unpublished ; m - b .
lepetit and g. pastor , phys .
b * 56 * , 4447 ( 1997 ) m. kaburagi , h. kawamura and t. hikihara , j. phys . soc . jpn .
* 68*. 3185 ( 1999 ) t. hikihara , j. phys .
. jpn . * 71 * , 319 ( 2002 ) ; t. hikihara , m. kaburagi , h. kawamura and t. tonegawa , phys .
b * 63 * , 174430 ( 2001 ) s. pati , s. ramasesha and d. sen , phys . rev .
b * 55 * , 8894 ( 1996 ) ; j. phys . cond
* 9 * , 8707 ( 1997 ) ; t. tonegawa _ et al _ , j. mag . * 177 - 181 * , 647 ( 1998 ) , ( cond - mat/9712298 ) m.  azzouz , l.  chen and s.  moukouri , phys .
b * 50 * 6223 ( 1994 ) ; s.  r.  white , r.  m.  noack and d.  j.  scalapino , phys .
lett .  * 73 * 886 ( 1994 ) ; k.  hida , j.  phys .
jpn .  * 64 * 4896 ( 1995 ) ; t.  narushima , t.  nakamura and s.  takada , j.  phys .
* 64 * 4322 ( 1995 ) ; u.  schollwck and d.  ko , phys .
b * 53 * 240 ( 1996 ) ; g. sierra , m. a. martn - delgado , s. white and j. dukelsky , phys .
b * 59 * , 7973 ( 1999 ) ; s.
white , phys . rev .
b * 53 * , 52 ( 1996 ) a. kawaguchi , a. koga , k. okunishi and n. kawakami , cond - mat/0205635 a. trumper and c. gazza , phys . rev .
b * 64 * , 134408 ( 2001 ) m. roger , phys . rev .
b * 63 * , 144433 ( 2001 ) t. hikihara and a. furusaki , phys . rev .
b * 63 * , 134438 ( 2001 ) a. laeuchli , g. schmida and m. troyer , cond - mat/0206153 t. hikihara , t. momoi and x. hu , cond - mat/0206102 t. nunner , p. brune , t. kopp , m. windt and m. grninger , cond - mat/0203472 y. honda and t. horiguchi , cond - mat/0106426 s. k. pati and r. r. p. singh , phys . rev .
b * 60 * , 7695 ( 1999 ) ; s. r. white and r. singh , phys .
lett . * 85 * , 3330 ( 2000 ) n. maeshima , m. hagiwara , y. narumi , k. kindo , t. c. kobayashi and k. okunishi , cond - mat/0208373 ; n. maeshima and k. okunishi , phys .
b * 62 * , 934 - 939 ( 2000 ) c. itoi and s. qin , phys .
b * 63 * , 224423 ( 2001 ) j. lou , j. dai , s. qin , z. su and l. yu , phys . rev .
b * 62 * , 8600 ( 2000 ) c. raghu , i. rudra , s. ramasesha and d. sen , phys . rev .
b * 62 * , 9484 ( 2000 ) a. honecker , m. kaulke and k.d .
schotte , eur .
j. b * 15 * , 423 ( 2000 ) a. gendiar and a. surda , phys .
b * 62 * , 3960 ( 2000 ) k. tandon et al , phys . rev .
b * 59 * , 396 ( 1999 ) ; r. citro , e. orignac , n. andrei , c. itoi and s. qin , j. phys.:cond . mat .
* 12 * , 3041 ( 2000 ) j. lou , c. chen and s. qin , phys . rev .
b * 64 * , 144403 ( 2001 ) x. wang and l. yu , phys
lett . * 84 * , 5399 ( 2000 ) s. k. pati , s. ramasesha and d. sen , in _ magnetism : molecules to materials iv _ ,
j.s . miller and m. drillon ( wiley - vch weinheim , 2002 ) , chapter 4 , ( cond - mat/0106621 ) l.  chen and s.  moukouri , phys .
b * 53 * 1866 ( 1996 ) ; s.  j.  qin , s.  d.  liang , z.  b.  su and l.  yu , phys .
b * 52 * r5475 ( 1995 ) ; r. noack in ref.@xcite , chap.1.3 ( part ii ) ; s. daul and r. noack , phys . rev .
b * 61 * , 12361 ( 2000 ) ; m. vojta , r. hetzel and r. noack , phys .
b * 60 * , r8417 ( 1999 ) k. penc , k. hallberg , f. mila and h. shiba , phys .
* 77 * , 1390 ( 1996 ) ; phys .
b , * 55 * , 15475 ( 1997 ) n. bulut , cond - mat/0207186 , to appear in advances in physics , * 51 * , no . 6 ( 2002 )
a. malvezzi , t. paiva and r. dos santos , cond - mat/0205266 e. jeckelmann , cond - mat/0204244 y. zhang , c. wu and h. q.lin , phys .
b * 65 * , 115101 ( 2002 ) c. aebischer , d. baeriswyl and r. m. noack , phys . rev . lett .
* 86 * , 468 ( 2001 ) r. arita and h. aoki , phys .
b * 61 * , 12261 ( 2000 ) ; r. arita , y. shimoi , k. kuroki and h. aoki , phys . rev .
b * 57 * , 10609 ( 1998 ) s. daul , cond - mat/9911361 s. qin , j. lou , t. xiang , g - s .
tian and z. su , cond - mat/0004162 a. aligia , k. hallberg , c. batista and g. ortiz , phys .
b * 61 * , 7883 ( 2000 ) v. meden , w. metzner , u. schollwoeck , o. schneider , t. stauber and k. schoenhammer , , europhys .
j b * 16 * , 631 ( 2000 ) p. farkasovsky , phys . rev .
b * 65 * , 081102 ( 2002 ) c. schuster , r. roemer and m. schreiber , phys . rev .
b * 65 * , 115114 ( 2002 ) k. hida , phys .
lett . * 86 * , 1331 ( 2001 ) w. lay and j. rudnick , phys .
lett . * 88 * , 057203 ( 2002 ) ; y. nishiyama , j.phys .
a*34 * , 11215 ( 2001 ) ; s. g. chung , phys .
e * 62 * , 3262 ( 2000 ) b. srinivasan and m. b. lepetit , phys .
b * 66 * , 024421 ( 2002 ) ; h. sakamoto , t. momoi and k. kubo , phys .
b * 65 * , 224403 ( 2002 ) k. hamacher , c. gros and w. wenzel , phys . rev
. lett . * 88 * , 217203 ( 2002 ) m. vojta , a. huebsch and r. m. noack , phys .
b * 63 * , 045105 ( 2001 ) ; z. weihong , j. oitmaa , c.j . hamer and r.j .
bursill , j. phys .
c * 13 * , 433 ( 2001 ) youngho park , s. liang and t. k. lee , phys .
b * 59 * , 2587 ( 1999 ) s. rommer , s. r. white and d. j. scalapino , phys . rev .
b * 61 * , 13424 ( 2000 ) r. m. noack , s. r. white and d. j. scalapino , phys
lett . * 73 * , 882 ( 1994 ) ; s. r. white , r. m. noack and d. j. scalapino , j. low temp . phys . * 99 * , 593 ( 1995 ) ; r. m. noack , s. r. white and d. j. scalapino , europhys . lett . * 30 * , 163 ( 1995 ) ; c. a. hayward , d. poilblanc , r. m. noack , d. j. scalapino and w. hanke , phys .
lett . * 75 * , 926 ( 1995 ) ; s. white and d. scalapino , phys . rev . lett . * 81 * , 3227 ( 1998 ) ; e. jeckelmann , d. scalapino and s. white , phys . rev .
b * 58 * , 9492 ( 1998 ) s. nishimoto , e. jeckelmann and d. scalapino , cond - mat/0208189 j. b. marston , j. o. fjaerestad and a. sudbo , phys . rev . lett . * 89 * , 056404 ( 2002 ) u.schollwoeck , sudip chakravarty , j. o. fjaerestad , j. b. marston and m. troyer , cond - mat/0209444 l. g. caron and c. bourbonnais , phys . rev .
b * 66 * , 045101 ( 2002 ) g. fano , f. ortolani , a. parola and l. ziosi , phys . rev .
b * 60 * , 15654 ( 1999 ) n. shibata and d. yoshioka , phys .
86 * , 5755 ( 2001 ) ; physica e * 12 * , 43 ( 2002 ) d. yoshioka and n. shibata , j. phys .
. jpn . * 70 * , 3690 ( 2001 ) v. meden and u. schollwck , cond - mat/0209588 s. liang and h. pang , europhys .
lett . * 32 * , 173 ( 1995 ) s. white and d. scalapino , phys . rev .
b * 55 * , 14701 ( 1997 ) s. white and d. scalapino , phys . rev .
b * 55 * , 6504 ( 1997 ) s. white and d. scalapino , phys . rev .
b * 57 * , 3031 ( 1998 ) e. arrigoni , a. p. harju , w. hanke , b. brendel and s.a .
kivelson , phys .
b * 65 * , 134503 ( 2002 ) ; j. bonca , j. e. gubernatis , m. guerrero , e. jeckelmann and s. r. white , phys . rev .
b * 61 * , 3251 ( 2000 ) ; a. l. chernyshev , s. white and a. h. castro neto , phys . rev .
b , * 65 * , 214527 ( 2002 ) ; s. white and d. scalapino , phys . rev .
b * 61 * , 6320 ( 2000 ) s. white and d. scalapino , phys . rev . lett . *
80 * , 1272 ( 1998 ) ( see also cond - mat/9907375 ) i. p. mcculloch , a. r. bishop and m. gulacsi , phil . mag .
b * 81 * , 1603 ( 2001 ) s. r. white and i. affleck , phys . rev .
b * 64 * , 024411 ( 2001 ) t. xiang , j. lou and z. su , phys . rev .
b * 64 * , 104414 ( 2001 ) ; p. henelius , phys . rev .
b * 60 * , 9561 ( 1999 ) s. nishimoto , e. jeckelmann , f. gebhard and r. noack , phys .
b * 65 * , 165114 ( 2002 ) n. maeshima , y. hieida , y. akutsu , t. nishino and k. okunishi , phys . rev .
e * 64 * ( 2001 ) 016705 m.a .
martin - delgado , j. rodriguez - laguna and g. sierra , nucl .
b * 601 * , 569 ( 2001 ) p. henelius , phys .
b * 60 * , 9561 ( 1999 ) s. moukouri and l. g. caron , cond - mat/0210668 t.  a.  costi , p.  schmitteckert , j.  kroha and p.  wlfle , phys .
lett .  * 73 * , 1275 ( 1994 ) ; s.  eggert and i.  affleck , phys .  rev .  lett .  * 75 * , 934 ( 1995 ) ; e.  s.  srensen and i. affleck , phys .
b * 51*. 16115 ( 1995 ) ; x.  q.  wang and s.  mallwitz , phys .  rev .
b * 53 * , r492 ( 1996 ) ; w.  wang , s.  j.  qin , z.  y.  lu , l. yu and z.  b.  su , phys .
b * 53 * , 40 ( 1996 ) ; c.  c.  yu and m.  guerrero , phys . rev .
b * 54 * , 15917 ( 1996 ) ; a. furusaki and t. hikihara , phys . rev .
b * 58 * , 5529 ( 1998 ) k. hallberg and r. egger , phys . rev .
b * 55 * , 8646 ( 1997 ) c. schuster and u. eckern , ann . phys .
( leipzig ) , * 11 * , 901 ( 2002 ) ( cond - mat/0201390 ) ; w. zhang , j. igarashi and p. fulde , phys .
b * 56 * , 654 ( 1997 ) v. meden , w. metzner , u. schollwoeck and k. schoenhammer , phys . rev .
b * 65 * , 045318 ( 2002 ) c.  c.  yu and s.  r.  white , phys .  rev .  lett .
* 71 * , 3866 ( 1993 ) ; c.  c.  yu and s.  r.  white , physica b * 199 * , 454 ( 1994 ) s.  moukouri and l.  g.  caron , phys .  rev .
b * 52 * , 15723 ( 1995 ) ; n.  shibata , t.  nishino , k.  ueda and c.  ishii , phys .
b * 53 * , r8828 ( 1996 ) ; m.  guerrero and r.  m.  noack , phys .
b * 53 * , 3707 ( 1996 ) h.  otsuka and t.  nishino , phys .
b * 52 * , 15066 ( 1995 ) ; s.  moukouri , l.  g.  caron , c.  bourbonnais and l.  hubert , phys .
b * 51 * , 15920 ( 1995 ) j. c. xavier , e. novais and e. miranda , phys .
b * 65 * , 214406 ( 2002 ) t. yamamoto , r. manago and y. mori , cond - mat/0204588 i. p. mcculloch , a. juozapavicius , a. rosengren and m. gulacsi , phys . rev .
b * 65 * , 052410 ( 2002 ) m. guerrero and r.m .
noack , phys . rev .
b * 63 * , 144423 ( 2001 ) s. watanabe , j. phys .
. jpn . * 69 * , 2947 ( 2000 ) i. p. mcculloch , m. gulacsi , s. caprara and a. juozapavicius , j. low temp .
117 , 323 ( 1999 ) j. xavier , r. pereira , e. miranda and i. affleck , cond - mat/0209623 s. watanabe , y. kuramoto , t. nishino and n. shibata , j. phys . soc .
jpn * 68 * , 159 ( 1999 ) j. moreno , s. qin , p. coleman and l. yu , phys .
b * 64 * , 085116 ( 2001 ) s.  moukouri , l.  chen and l.  g.  caron , phys .  rev .
b * 53 * , r488 ( 1996 ) j. riera , k. hallberg and e. dagotto , phys . rev
. lett . * 79 * , 713 ( 1997 ) ; e. dagotto et al . ,
b * 58 * , 6414 ( 1998 ) d. garca , k. hallberg , c. batista , m. avignon and b. alascio , phys .
rev lett . , * 85 * , 3720 ( 2000 ) ; d. j. garca , k. hallberg , c. d. batista , s. capponi , d. poilblanc , m. avignon and b. alascio , phys . rev .
b. * 65 * , 134444 ( 2002 ) b. ammon and m. imada , j. phys .
. jpn . * 70 * , 547 ( 2001 ) s. ramasesha et al , phys .
b * 54 * , 7598 ( 1996 ) ; synth .
metals * 85 * , 1019 ( 1997 ) e.  s.  srensen and i. affleck , phys .
b * 51 * , 16115 ( 1995 ) i. mcculloch and m. gulacsi , europhys .
lett . * 57 * , 852 ( 2002 ) g. sierra and t. nishino , nucl .
b * 495 * , 505 ( 1997 ) ( cond - mat/9610221 ) t. xiang , phys .
b * 53 * , r10445 ( 1996 ) m.s.l .
du croo de jongh and j.m.j .
van leeuwen , phys .
b * 57 * , 8494 ( 1998 ) m - b .
lepetit , m. cousy and g. pastor , eur .
j. b * 13 * , 421 ( 2000 ) ; h. otsuka , phys . rev .
b * 53 * , 14004 ( 1996 ) j. dukelsky and g. dussel , phys . rev .
b * 59 * , r3005 ( 1999 ) and references therein .
j. dukelsky and g. sierra , phys .
lett . * 83 * , 172 ( 1999 ) m. a. martn - delgado and g. sierra , phys .
lett . * 83 * , 1514 ( 1999 ) s. glazek and k. wilson , phys .
d * 48 * , 5863 ( 1993 ) ; * 49 * , 4214 ( 1994 ) t. eller , h .- c .
pauli and s. brodsky , phys .
d * 35 * , 1493 ( 1987 ) t. nishino , j. phys .
. jpn . * 64 * , 3598 ( 1995 ) ; see also t. nishino in ref .
@xcite , chap .
5(i ) ; t. nishino and k. okunishi in _ strongly correlated magnetic and superconducting systems _ , ed . g. sierra and m. a. martn - delgado ( springer , berlin , 1997 ) h. trotter , proc .
* 10 * , 545 ( 1959 ) ; m. suzuki , prog .
. phys . * 56 * , 1454 ( 1976 ) ; r. feynman and a. hibbs _ quantum mechanics and path integrals _ ( mcgraw - hill , 1965 ) a. gendiar and a. surda , phys .
b * 63 * , 014401 ( 2001 ) t. nishino and n. shibata , j. phys . soc . jpn . *
68 * , 3501 ( 1999 ) ; t. enss and u. schollwoeck , j. phys .
a * 34 * , 7769 ( 2001 ) a. kemper , a. schadschneider and j. zittartz , j. phys .
a : math . gen . * 34 * , l279 ( 2001 ) ; t. enss and u. schollwck , j phys a : math . gen . *
34 * , 7769 ( 2001 ) r. baxter , j. math . phys . * 9 * , 650 ( 1968 ) ; j. stat phys . * 19 * , 461 ( 1978 ) t. nishino and k. okunishi , j. phys . soc
. jpn . * 65*,891 ( 1996 ) ; ibid . * 66 * , 3040 ( 1997 ) ; t. nishino , k. okunishi and m. kikuchi , phys .
a * 213 * , 69 ( 1996 ) h. takasaki , t. nishino and y. hieida , j. phys .
. jpn . * 70 * , 1429 ( 2001 ) t. nishino , y. hieida , k. okunishi , n. maeshima , y. akutsu and a. gendiar , prog . theor . phys . * 105 * , 409 ( 2001 ) c. ritter and g. von gehlen , in `` quantization , gauge theory and strings '' , ed .
et al . , vol.i , p.563 - 578 , scientific world publ.co .
( 2001 ) ( cond - mat/0009255 ) t. nishino and k. okunishi , j. phys .
. jpn . * 67 * , 3066 ( 1998 ) ; t. nishino , k. okunishi , y. hieida , n. maeshima and y. akutsu , nucl . phys .
b * 575 * , 504 ( 2000 ) e. carlon and a. drzewiski , phys . rev . lett . *
79 * , 1591 ( 1997 ) ; phys .
e * 57 * , 2626 ( 1998 ) ; e. carlon , a. drzewiski and j. rogiers , phys .
b * 58 * , 5070 ( 1998 ) ; a. drzewiski , a. ciach and a. maciolek , eur .
j. b * 5 * , 825 ( 1998 ) ; phys .
e * 60 * , 2887 ( 1999 ) a. drzewinski , a. maciolek and r. evans , phys .
lett . * 85 * , 3079 ( 2000 ) ; a. drzewinski , phys . rev .
e * 62 * , 4378 ( 2000 ) m .- c .
chung , m. kaulke , i. peschel , m. pleimling and w. selke , eur .
j. b * 18 * , 655 ( 2000 ) e. carlon and f. igli , phys .
b * 57 * , 7877 ( 1998 ) ; phys . rev .
b * 59 * , 3783 ( 1999 ) ; e. carlon , c. chatelain and b. berche , phys . rev .
b * 60 * , 12974 ( 1999 ) j. kondev and j. marston , nucl .
b * 497 * , 639 ( 1997 ) ; t. senthil , b. marston and m. fisher , phys .
b * 60 * , 4245 ( 1999 ) ; j. marston and s. tsai , phys .
lett . * 82 * , 4906 ( 1999 ) ; s. tsai and j. marston , ann .
( leipzig ) * 8 * , special issue , 261(1999 ) m. kaulke and i. peschel , eur .
j. b * 5 * , 727 ( 1998 ) y. hieida , j. phys .
. jpn . * 67 * , 369 ( 1998 ) i. peschel and m. kaulke , in ref .
@xcite , chap .
3.1(ii ) e. carlon , m. henkel and u. schollwck , eur . phys .
j. b * 12 * , 99(1999 ) ; e. carlon , m. henkel and u. schollwoeck , phys .
e * 63 * , 036101 ( 2001 ) j. hooyberghs , e. carlon and c. vanderzande , phys .
e * 62 * , 036124 ( 2001 ) s.g .
chung , phys .
b * 60 * , 11761(1999 ) r. bursill , t. xiang and g. gehring , j. phys .
c * 8 * , l583 ( 1996 ) h. trotter , proc .
* 10 * , 545 ( 1959 ) ; m. suzuki , prog .
. phys . * 56 * , 1454 ( 1976 ) x. wang and t. xiang , phys .
b * 56 * , 5061 ( 1997 ) n. shibata and k. ueda , j. phys . soc
vol.70 , 3690 ( 2001 ) n. shibata , j. phys .
. jpn . * 66 * , 2221 ( 1997 ) t. xiang , phys .
b * 58 * , 9142 ( 1998 ) n. shibata , b. ammon , t. troyer , m. sigrist and k. ueda , j. phys .
. jpn . * 67 * , 1086 ( 1998 ) k. maisinger , u. schollwck , s. brehmer , h - j .
mikeska and s. yamamoto , phys .
b * 58 * , r5908 ( 1998 ) s. rommer and s. eggert , phys . rev .
b * 59 * , 6301 ( 1999 ) i. maruyama , n. shibata and k. ueda , phys . rev .
b * 65 * , 174421 ( 2002 ) k. maisinger and u. schollwck , phys . rev . lett .
* 81 * , 445 ( 1999 ) b. ammon , m. troyer , t. rice and n. shibata , phys .
lett . * 82 * , 3855 ( 1999 ) ; n. shibata and h. tsunetsugu , j. phys .
. jpn . * 68 * , 3138 ( 1999 ) f. naef and x. wang , phys .
lett . * 84 * , 1320 ( 2000 ) a. klmper , r. raupach and f. schnfeld , phys .
b * 59 * , 3612 ( 1999 ) s. moukouri and l. caron , phys . rev . lett . *
77 * , 4640 ( 1996 ) .
s. moukouri and l. caron , see ref.@xcite , chap .
4.5(ii ) s. moukouri , cond - mat/0011169 w. zhang , j. igarashi and p. fulde , j. phys .
. jpn . * 66 * , 1912 ( 1997 ) l. caron and s. moukouri , phys .
lett . * 76 * , 4050 ( 1996 ) ; phys .
b * 56 * , r8471 ( 1997 ) e. jeckelmann and s. white , phys . rev .
b * 57 * , 6376 ( 1998 ) r. noack , s. white and d. scalapino in _ computer simulations in condensed matter physics vii _ , edited by d. landau , k .- k . mon and h .- b .
schttler ( springer verlag , heidelberg and berlin , 1994 ) e. jeckelmann , c. zhang and s. white in ref .
@xcite , chap . 5.1
( ii ) c. zhang , e. jeckelmann and s. white , phys .
lett . * 80 * , 2661 ( 1998 ) ; e. jeckelmann , c. zhang and s. r. white , physical review b * 60 * , 7950 ( 1999 ) y. nishiyama , cond - mat/0102123 eur .
j. b * 12 * , 547 ( 1999 ) r. bursill , y. mckenzie and c. hammer , phys .
lett . * 80 * , 5607 ( 1998 ) ; * 83 * , 408 ( 1999 ) ; r. bursill , phys . rev .
b * 60 * , 1643 ( 1999 ) r. pai , r. pandit , h. krishnamurthy and s. ramasesha , phys .
lett . * 76 * , 2937 ( 1996 ) ( see also the comment by n. v. prokofev and b. v. svistunov , phys . rev . lett . * 80 * , 4355 ( 1998 ) ) ; s. rapsch , u. schollwck and w. zwerger , europhys . lett . *
46 * , 559 ( 1999 ) t. khner and h. monien , phys .
b * 58 * , r14741 ( 1998 ) i. peschel and m - c . chung , j. phys .
gen . * 32 * , 8419 ( 1999 ) p. maurel and m - b .
lepetit , phys .
b * 62 * , 10744 ( 2000 ) ; p. maurel , m - b .
lepetit and d. poilblanc , eur .
j. * b 21 * , 481 ( 2001 ) k. hida , j. phys .
. jpn . * 65 * , 895 ( 1996 ) and 3412 ( 1996 ) ( erratum ) ; j. phys .
. jpn . * 66 * , 330 ( 1997 ) ; j. phys .
. jpn . * 66 * , 3237 ( 1997 ) ; prog .
. suppl . * 145 * , 320 ( 2002 ) k. hida , j. phys .
. jpn . * 68 * , 3177 ( 1999 ) k. hida , phys .
lett . * 83 * , 3297 ( 1999 ) f. igloi , r. juhasz and p. lajko , phys .
lett . * 86 * , 1343 ( 2001 ) r. mlin , y - c .
lin , p. lajk , h. rieger and f. igli , phys .
b * 65 * , 104415 ( 2002 ) p. schmitteckert , t. schulze , c. schuster , p. schwab and u. eckern , phys . rev . lett .
* 80 * , 560 ( 1998 ) ; p. schmitteckert and u. eckern , phys .
b * 53 * , 15397 ( 1996 ) d. weinmann , p. schmitteckert , r. jalabert and j. pichard , eur
j. b * 19 * , 139 - 156 ( 2001 ) p. schmitteckert , r. jalabert , d. weinmann and j. l. pichard , phys .
lett . * 81 * , 2308 ( 1998 ) e. carlon , p. lajko and f. igli , phys .
lett . * 87 * , 277201 ( 2001 ) g. fano , f. ortolani and l. ziosi , j. chem . phys . * 108 * , 9246 ( 1998),(cond - mat/9803071 ) ; r. bursill and w. barford , phys . rev . lett .
* 82 * , 1514 ( 1999 ) m. exler and j. schnack , cond - mat/0205068 b. normand , x. wang , x. zotos and daniel loss , phys .
b * 63 * , 184409 ( 2001 ) c. raghu , y. anusooya pati and s. ramasesha , journal - ref : j. phys .
a * 34 * , 11215 ( 2001 ) w. barford and r. bursill , chem .
. lett . * 268 * , 535(1997 ) ; w. barford , r. bursill and m. lavrentiev , j. phys : cond .
matt , * 10 * , 6429 ( 1998 ) ; w. barford in ref.@xcite , chap.2.3 ( part ii ) and references therein ; m. lavrentiev , w. barford , s. martin , h. daly , r. bursill , physical review b * 59 * , 9987 ( 1999 ) e. carlon , a. drzewinski and j. van leeuwen , j. chem . phys . * 117 * , 2425 ( 2002 ) ; phys . rev .
e * 64 * , 010801(r ) ( 2001 ) ; m. paessens and g. schtz , phys .
e * 66 * , 021806 ( 2002 ) ; m. a. martn - delgado , j. rodriguez - laguna and g. sierra , phys .
b * 65 * , 155116 ( 2002 ) s. white and r. martin , j. chem . phys .
* 110 * , 4127 ( 1999 ) ; see also s. white in @xcite , chap .
s. white in ref .
@xcite , chap . 2.1 .
. legeza , j. rder and b.a .
hess , cond - mat/0208187 s. daul , i. ciofini , c. daul and s. r. white , int .
j. of quantum chem .
, * 79 * , 331(2000 ) ( cond - mat/9912348 ) c. bauschlicher and p. taylor , j. chem . phys .
* 85 * , 2779 ( 1986 ) . legeza , j. roder and b. a. hess , cond - mat/0204602 w. hofstetter , phys .
lett . * 85 * , 1508 ( 2000 ) m. a. cazalilla and j. b. marston , phys . rev .
lett . * 88 * , 256403 ( 2002 ) k. hallberg , phys .
b * 52 * , 9827 ( 1995 ) .
e. r. gagliano and c. a. balseiro , phys .
lett . * 59 * , 2999 ( 1987 ) .
g. grosso and g. partori parravicini , in _ memory function approaches to stochastic problems in condensed matter _ , adv . in chemical physics ,
* 62 * , 133 ( wiley , n. y. , 1985 ) p. horsch and w. von der linden , z. phys .
b * 72 * 181 ( 1981 ) j. des cloiseaux and j. j. pearson , phys . rev . * 128 * , 2131 ( 1962 ) t. yamada , prog .
. jpn . * 41 * , 880 ( 1969 ) ; l. d. fadeev and l. a. takhtajan , phys . lett . * 85 a * , 375 ( 1981 ) g. mller , h. thomas , h. beck and j. bonner , phys .
b * 24 * , 1429 ( 1981 ) and references therein .
s. haas , j. riera and e. dagotto , phys .
b * 48 * , 3281 ( 1993 ) t. khner and s. white , phys . rev .
b * 60 * , 335(1999 ) p. brune , g. japaridze and a. p. kampf , cond - mat/0106007 ( unpublished ) t. d. khner , s. r. white and h. monien , phys .
b * 61 * , 12474 ( 2000 ) s. s. kancharla and c. j. bolech , phys . rev .
b * 64 * , 085119 ( 2001 ) y. anusooya , s. pati and s. ramasesha , j. chem . phys . *
106 * , 1 ( 1997 ) ; s. ramasesha , k. tandon , y. anusooya and s. pati , proc . of spie ,
* 3145 * , 282 ( 1997 ) ; s. ramasesha , z. shuai and j. brdas , chem .
. lett . * 245 * , 224 ( 1995 ) e. jeckelmann , phys .
b * 66 * , 045114 ( 2002 ) e. jeckelmann and f. gebhard , phys . rev .
lett . * 85 * , 3910 ( 2000 ) ; f. essler , f. gebhard and e. jeckelmann , phys . rev .
b * 64 * , 125119 ( 2001 ) ; e. jeckelmann , cond - mat/0208480 h.  b.  pang , h.  akhlaghpour and m.  jarrell , phys.rev .
b * 53 * 5086 ( 1996 ) p. hohenberg and w. brinkman , phys .
b * 10 * , 128 ( 1974 ) x. wang , k. hallberg and f. naef in ref.@xcite , chap.7(i ) t. mutou , n. shibata and k. ueda , phys . rev .
lett * 81 * , 4939 ( 1998 ) ( erratum : phys .
lett . * 82 * , 3727 ( 1999 ) ) f. naef , x. wang , x. zotos and w. von der linden , phys .
b * 60 * , 359 ( 1999 )