Density Matrix Renormalization Group of Gapless Systems

We investigate convergence of the density matrix renormalization group (DMRG) in the thermodynamic limit for gapless systems. Although the DMRG correlations always decay exponentially in the thermodynamic limit, the correlation length at the DMRG fixed-point scales as $\xi \sim m^{1.3}$, where $m$ is the number of kept states, indicating the existence of algebraic order for the exact system. The single-particle excitation spectrum is calculated, using a Bloch-wave ansatz, and we prove that the Bloch-wave ansatz leads to the symmetry $E(k)=E(\pi -k)$ for translationally invariant half-integer spin-systems with local interactions. Finally, we provide a method to compute overlaps between ground states obtained from different DMRG calculations.Comment: 11 pages, RevTex, 5 figure

The Dynamical Mean Field Theory phase space extension and critical properties of the finite temperature Mott transition

We consider the finite temperature metal-insulator transition in the half filled paramagnetic Hubbard model on the infinite dimensional Bethe lattice. A new method for calculating the Dynamical Mean Field Theory fixpoint surface in the phase diagram is presented and shown to be free from the convergence problems of standard forward recursion. The fixpoint equation is then analyzed using dynamical systems methods. On the fixpoint surface the eigenspectra of its Jacobian is used to characterize the hysteresis boundaries of the first order transition line and its second order critical end point. The critical point is shown to be a cusp catastrophe in the parameter space, opening a pitchfork bifurcation along the first order transition line, while the hysteresis boundaries are shown to be saddle-node bifurcations of two merging fixpoints. Using Landau theory the properties of the critical end point is determined and related to the critical eigenmode of the Jacobian. Our findings provide new insights into basic properties of this intensively studied transition.Comment: 11 pages, 12 figures, 1 tabl

Scale effects for strength, ductility, and toughness in "brittle” materials

Decreasing scales effectively increase nearly all important mechanical properties of at least some "brittle” materials below 100 nm. With an emphasis on silicon nanopillars, nanowires, and nanospheres, it is shown that strength, ductility, and toughness all increase roughly with the inverse radius of the appropriate dimension. This is shown experimentally as well as on a mechanistic basis using a proposed dislocation shielding model. Theoretically, this collects a reasonable array of semiconductors and ceramics onto the same field using fundamental physical parameters. This gives proportionality between fracture toughness and the other mechanical properties. Additionally, this leads to a fundamental concept of work per unit fracture area, which predicts the critical event for brittle fracture. In semibrittle materials such as silicon, this can occur at room temperature when the scale is sufficiently small. When the local stress associated with dislocation nucleation increases to that sufficient to break bonds, an instability occurs resulting in fractur

The Dynamical Mean Field Theory phase space extension and critical properties of the finite temperature Mott transition

We consider the finite temperature metal-insulator transition in the half filled paramagnetic Hubbard model on the infinite dimensional Bethe lattice. A new method for calculating the Dynamical Mean Field Theory fixpoint surface in the phase diagram is presented and shown to be free from the convergence problems of standard forward recursion. The fixpoint equation is then analyzed using dynamical systems methods. On the fixpoint surface the eigenspectra of its Jacobian is used to characterize the hysteresis boundaries of the first order transition line and its second order critical end point. The critical point is shown to be a cusp catastrophe in the parameter space, opening a pitchfork bifurcation along the first order transition line, while the hysteresis boundaries are shown to be saddle-node bifurcations of two merging fixpoints. Using Landau theory the properties of the critical end point is determined and related to the critical eigenmode of the Jacobian. Our findings provide new insights into basic properties of this intensively studied transition.Comment: 11 pages, 12 figures, 1 tabl

Clustering in mixing flows

We calculate the Lyapunov exponents for particles suspended in a random three-dimensional flow, concentrating on the limit where the viscous damping rate is small compared to the inverse correlation time. In this limit Lyapunov exponents are obtained as a power series in epsilon, a dimensionless measure of the particle inertia. Although the perturbation generates an asymptotic series, we obtain accurate results from a Pade-Borel summation. Our results prove that particles suspended in an incompressible random mixing flow can show pronounced clustering when the Stokes number is large and we characterise two distinct clustering effects which occur in that limit.Comment: 5 pages, 1 figur

A Two-dimensional Infinte System Density Matrix Renormalization Group Algorithm

It has proved difficult to extend the density matrix renormalization group technique to large two-dimensional systems. In this Communication I present a novel approach where the calculation is done directly in two dimensions. This makes it possible to use an infinite system method, and for the first time the fixed point in two dimensions is studied. By analyzing several related blocking schemes I find that there exists an algorithm for which the local energy decreases monotonically as the system size increases, thereby showing the potential feasibility of this method.Comment: 5 pages, 6 figure

The critical behaviour of the 2D Ising model in Transverse Field; a Density Matrix Renormalization calculation

We have adjusted the Density Matrix Renormalization method to handle two dimensional systems of limited width. The key ingredient for this extension is the incorporation of symmetries in the method. The advantage of our approach is that we can force certain symmetry properties to the resulting ground state wave function. Combining the results obtained for system sizes up-to $30 \times 6$ and finite size scaling, we derive the phase transition point and the critical exponent for the gap in the Ising model in a Transverse Field on a two dimensional square lattice.Comment: 9 pages, 8 figure

Thermodynamic limit of the density matrix renormalization for the spin-1 Heisenberg chain

The density matrix renormalization group (``DMRG'') discovered by White has shown to be a powerful method to understand the properties of many one dimensional quantum systems. In the case where renormalization eventually converges to a fixed point we show that quantum states in the thermodynamic limit with periodic boundary conditions can be simply represented by a special type of product ground state with a natural description of Bloch states of elementary excitations that are spin-1 solitons. We then observe that these states can be rederived through a simple variational ansatz making no reference to a renormalization construction. The method is tested on the spin-1 Heisenberg model.Comment: 13 pages uuencoded compressed postscript including figure

The Density Matrix Renormalization Group technique with periodic boundary conditions

The Density Matrix Renormalization Group (DMRG) method with periodic boundary conditions is introduced for two dimensional classical spin models. It is shown that this method is more suitable for derivation of the properties of infinite 2D systems than the DMRG with open boundary conditions despite the latter describes much better strips of finite width. For calculation at criticality, phenomenological renormalization at finite strips is used together with a criterion for optimum strip width for a given order of approximation. For this width the critical temperature of 2D Ising model is estimated with seven-digit accuracy for not too large order of approximation. Similar precision is reached for critical indices. These results exceed the accuracy of similar calculations for DMRG with open boundary conditions by several orders of magnitude.Comment: REVTeX format contains 8 pages and 6 figures, submitted to Phys. Rev.

Unmixing in Random Flows

We consider particles suspended in a randomly stirred or turbulent fluid. When effects of the inertia of the particles are significant, an initially uniform scatter of particles can cluster together. We analyse this 'unmixing' effect by calculating the Lyapunov exponents for dense particles suspended in such a random three-dimensional flow, concentrating on the limit where the viscous damping rate is small compared to the inverse correlation time of the random flow (that is, the regime of large Stokes number). In this limit Lyapunov exponents are obtained as a power series in a parameter which is a dimensionless measure of the inertia. We report results for the first seven orders. The perturbation series is divergent, but we obtain accurate results from a Pade-Borel summation. We deduce that particles can cluster onto a fractal set and show that its dimension is in satisfactory agreement with previously reported in simulations of turbulent Navier-Stokes flows. We also investigate the rate of formation of caustics in the particle flow.Comment: 39 pages, 8 figure