The Density Matrix Renormalization Group
The Density Matrix Renormalization Group (DMRG) method is a widely-utilized sophisticated algorithm to obtain low-lying eigenvalues and eigenvectors of very large matrices such as those found in quantum many-body problems. DMRG possesses features that make it extremely powerful: it can treat systems with hundreds of quantum spins or electrons, provide extremely accurate ground-state energies, and compute gaps in low-dimensional systems. Together with quantum Monte Carlo methods, it dominates most of the numerical research in the field of strongly correlated electron systems.
The DMRG is particularly suitable for one-dimensional systems, where it is definitely the method of choice if one is only interested in ground state properties and low-lying eigenstates. The algorithm is much more efficient in systems with open boundary conditions, and quasi-2d systems (ladders) with cylindrical boundary conditions, where convergence is achieved with a relatively small number of states.
S. R. White
Density matrix formulation for quantum renormalization groups
Phys. Rev. Lett. 69, 2863 (1992).
S. R. White
Density-matrix algorithms for quantum renormalization groups
Phys. Rev. B 48, 10345 (1993).
The density-matrix renormalization group
Rev. Mod. Phys. 77, 259 (2005).
Number of eigenstates and energies to calculate. Default is 1, should be set to 2 to calculate gaps.
Number of DMRG sweeps. Each sweep involves a left-to-right half-sweep, and a right-to-left half-sweep.
Number of initial states to grow the DMRG blocks. If not specified, the algorithm will use default value used of 20.
Number of DMRG states kept on each half sweep. The user should specify either 2*SWEEPS different values of STATES or one MAXSTATES or NUMSTATES value.
Maximum number of DMRG states kept. The user should choose to specify either STATES values for each half-sweep, or a MAXSTATES or NUMSTATES that the program will use to grow the basis. The program will automatically determine how many states to use for each sweep, growing in steps of STATES/(2*SWEEPS) until reaching MAXSTATES.
Constant number of DMRG states kept for all sweeps.
TRUNCATION_ERROR The user can choose to set the tolerance for the simulation, instead of the number of states. The program will automagically determine how many states to keep in order to satisfy this tolerance. Care must be taken, since this could lead to an uncontrollable growth in the basis size, and a crash as a consequence. It is therefore advisable to also specify the maximum number of states as a constraint, using either MAXSTATES or NUMSTATES, as explained before.
Quantum numbers conserved by the model of interest. They will be used in the code in order to reduce matrices in block form. If no value is specified for a particular quantum number, the program will work in the grand canonical. For instance in spin chains if you don't specify Sz_total, the program will run using a Hilbert space with dim=2^N states. Running in the "canonical" (by setting Sz_total=0, for instance) will improve performance considerably by working in a subspace with a reduced dimension. For an example of how to do this, take a look at the parms file included with the dmrg code.
If set to an integer > 0, it will print extra output information, such as density-matrix eigenvalues. There are different verbose levels up to a maximum of 3, for debugging purposes, although the user shouldn't need a level larger than 1.
(Available after v1.3b6) Starting sweep for resuming a simulation that was interrupted, or extending it with a new set of states.
Starting direction for the resumed simulation. Can assume the values 0 or 1, for "left-to-right" or "right-to-left", respectively. It has effect only in the presence of START_SWEEP. Its default value is 0.
Starting iteration for the resumed simulation. It has effect only in the presence of START_SWEEP. Its default value is 1.