Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: 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 )
Abstract:

Answer: the density matrix renormalization group ( dmrg ) has become a powerful numerical method that can be applied to low - dimensional strongly correlated fermionic and bosonic systems . it allows for a very precise calculation of static , dynamical and thermodynamical properties . its field of applicability has now extended beyond condensed matter , and is successfully used in statistical mechanics and high energy physics as well . in this article , we briefly review the main aspects of the method . we also comment on some of the most relevant applications so as to give an overview on the scope and possibilities of dmrg and mention the most important extensions of the method such as the calculation of dynamical properties , the application to classical systems , inclusion of temperature , phonons and disorder , field theory , time - dependent properties and the _ ab initio _ calculation of electronic states in molecules .