Efficiency of symmetric targeting for finite-T DMRG

Two targeting schemes have been known for the density matrix renormalization group (DMRG) applied to non-Hermitian problems; one uses an asymmetric density matrix and the other uses symmetric density matrix. We compare the numerical efficiency of these two targeting schemes when they are used for the finite temperature DMRG.Comment: 4 pages, 3 Postscript figures, REVTe

Hyperbolic Deformation on Quantum Lattice Hamiltonians

A group of non-uniform quantum lattice Hamiltonians in one dimension is introduced, which is related to the hyperbolic $1 + 1$-dimensional space. The Hamiltonians contain only nearest neighbor interactions whose strength is proportional to $\cosh j \lambda$, where $j$ is the lattice index and where $\lambda \ge 0$ is a deformation parameter. In the limit $\lambda \to 0$ the Hamiltonians become uniform. Spacial translation of the deformed Hamiltonians is induced by the corner Hamiltonians. As a simple example, we investigate the ground state of the deformed $S = 1/2$ Heisenberg spin chain by use of the density matrix renormalization group (DMRG) method. It is shown that the ground state is dimerized when $\lambda$ is finite. Spin correlation function show exponential decay, and the boundary effect decreases with increasing $\lambda$.Comment: 5 pages, 4 figure

Wave-packet dynamics at the mobility edge in two- and three-dimensional systems

We study the time evolution of wave packets at the mobility edge of disordered non-interacting electrons in two and three spatial dimensions. The results of numerical calculations are found to agree with the predictions of scaling theory. In particular, we find that the $k$-th moment of the probability density $(t)$ scales like $t^{k/d}$ in $d$ dimensions. The return probability $P(r=0,t)$ scales like $t^{-D_2/d}$, with the generalized dimension of the participation ratio $D_2$. For long times and short distances the probability density of the wave packet shows power law scaling $P(r,t)\propto t^{-D_2/d}r^{D_2-d}$. The numerical calculations were performed on network models defined by a unitary time evolution operator providing an efficient model for the study of the wave packet dynamics.Comment: 4 pages, RevTeX, 4 figures included, published versio

Stochastic Light-Cone CTMRG: a new DMRG approach to stochastic models

We develop a new variant of the recently introduced stochastic transfer-matrix DMRG which we call stochastic light-cone corner-transfer-matrix DMRG (LCTMRG). It is a numerical method to compute dynamic properties of one-dimensional stochastic processes. As suggested by its name, the LCTMRG is a modification of the corner-transfer-matrix DMRG (CTMRG), adjusted by an additional causality argument. As an example, two reaction-diffusion models, the diffusion-annihilation process and the branch-fusion process, are studied and compared to exact data and Monte-Carlo simulations to estimate the capability and accuracy of the new method. The number of possible Trotter steps of more than 10^5 shows a considerable improvement to the old stochastic TMRG algorithm.Comment: 15 pages, uses IOP styl

Product Wave Function Renormalization Group: construction from the matrix product point of view

We present a construction of a matrix product state (MPS) that approximates the largest-eigenvalue eigenvector of a transfer matrix T, for the purpose of rapidly performing the infinite system density matrix renormalization group (DMRG) method applied to two-dimensional classical lattice models. We use the fact that the largest-eigenvalue eigenvector of T can be approximated by a state vector created from the upper or lower half of a finite size cluster. Decomposition of the obtained state vector into the MPS gives a way of extending the MPS, at the system size increment process in the infinite system DMRG algorithm. As a result, we successfully give the physical interpretation of the product wave function renormalization group (PWFRG) method, and obtain its appropriate initial condition.Comment: 8 pages, 8 figure

Path integral Monte Carlo simulations of silicates

We investigate the thermal expansion of crystalline SiO$_2$ in the $\beta$-- cristobalite and the $\beta$-quartz structure with path integral Monte Carlo (PIMC) techniques. This simulation method allows to treat low-temperature quantum effects properly. At temperatures below the Debye temperature, thermal properties obtained with PIMC agree better with experimental results than those obtained with classical Monte Carlo methods.Comment: 27 pages, 10 figures, Phys. Rev. B (in press

Application of the density matrix renormalization group method to finite temperatures and two-dimensional systems

The density matrix renormalization group (DMRG) method and its applications to finite temperatures and two-dimensional systems are reviewed. The basic idea of the original DMRG method, which allows precise study of the ground state properties and low-energy excitations, is presented for models which include long-range interactions. The DMRG scheme is then applied to the diagonalization of the quantum transfer matrix for one-dimensional systems, and a reliable algorithm at finite temperatures is formulated. Dynamic correlation functions at finite temperatures are calculated from the eigenvectors of the quantum transfer matrix with analytical continuation to the real frequency axis. An application of the DMRG method to two-dimensional quantum systems in a magnetic field is demonstrated and reliable results for quantum Hall systems are presented.Comment: 33 pages, 18 figures; corrected Eq.(117

Staggered dimer order in S=1/2 quantum spin ladder system with four spin exchange

We study the S=1/2 quantum spin ladder system with the four-spin exchange, using density matrix renormalization group method and an exact diagonalization method. Recently, the phase transition in this system and its universality class are studied. But there remain controversies whether the phase transition is second order type or the other type and the nature of order parameter. There are arguments that the massless phase appears. But this does not agree with our previous result. Analyzing DMRG data, we try a new approach in order to determine a phase which appears after the phase transition. We find that the edge state appears in the open boundary condition, investigating excitation energies of states with higher magnetizations.Comment: Submitted to Phys. Rev. B, (REVTeX4

The Seroepidemiology of Haemophilus influenzae Type B Prior to Introduction of an Immunization Programme in Kathmandu, Nepal.

Haemophilus influenzae type b (Hib) is now recognized as an important pathogen in Asia. To evaluate disease susceptibility, and as a marker of Hib transmission before routine immunization was introduced in Kathmandu, 71 participants aged 7 months-77 years were recruited and 15 cord blood samples were collected for analysis of anti-polyribosylribitol phosphate antibody levels by enzyme-linked immunosorbent assay. Only 20% of children under 5 years old had levels considered protective (>0.15 µg/ml), rising to 83% of 15-54 year-olds. Prior to introduction of Hib vaccine in Kathmandu, the majority of young children were susceptible to disease