Possible persistent current in a ring made of the perfect crystalline insulator

A mesoscopic conducting ring pierced by magnetic flux is known to support the persistent electron current. Here we propose possibility of the persistent current in the ring made of the perfect crystalline insulator. We consider a ring-shaped lattice of one-dimensional "atoms" with a single energy level. We express the Bloch states in the lattice as a linear combination of atomic orbitals. The discrete energy level splits into the energy band which serves as a simple model of the valence band. We show that the insulating ring (with the valence band fully filled by electrons) supports a nonzero persistent current, because each atomic orbital overlaps with its own tail when making one loop around the ring. In the tight-binding limit only the neighboring orbitals overlap. In that limit the persistent current at full filling becomes zero which is a standard result.Comment: Conference proceedings. Accepted for publication in Physica

Numerical study of the hard-core Bose-Hubbard Model on an Infinite Square Lattice

We present a study of the hard-core Bose-Hubbard model at zero temperature on an infinite square lattice using the infinite Projected Entangled Pair State algorithm [Jordan et al., Phys. Rev. Lett. 101, 250602 (2008)]. Throughout the whole phase diagram our values for the ground state energy, particle density and condensate fraction accurately reproduce those previously obtained by other methods. We also explore ground state entanglement, compute two-point correlators and conduct a fidelity-based analysis of the phase diagram. Furthermore, for illustrative purposes we simulate the response of the system when a perturbation is suddenly added to the Hamiltonian.Comment: 8 pages, 6 figure

Infinite boundary conditions for matrix product state calculations

We propose a formalism to study dynamical properties of a quantum many-body system in the thermodynamic limit by studying a finite system with infinite boundary conditions (IBC) where both finite size effects and boundary effects have been eliminated. For one-dimensional systems, infinite boundary conditions are obtained by attaching two boundary sites to a finite system, where each of these two sites effectively represents a semi-infinite extension of the system. One can then use standard finite-size matrix product state techniques to study a region of the system while avoiding many of the complications normally associated with finite-size calculations such as boundary Friedel oscillations. We illustrate the technique with an example of time evolution of a local perturbation applied to an infinite (translationally invariant) ground state, and use this to calculate the spectral function of the S=1 Heisenberg spin chain. This approach is more efficient and more accurate than conventional simulations based on finite-size matrix product state and density-matrix renormalization-group approaches.Comment: 10 page

Tricritical point of J1-J2 Ising model on hyperbolic lattice

A ferromagnetic-paramagnetic phase transition of the two-dimensional frustrated Ising model on a hyperbolic lattice is investigated by use of the corner transfer matrix renormalization group method. The model contains ferromagnetic nearest-neighbor interaction J_1 and the competing antiferromagnetic interaction J_2. A mean-field like second-order phase transition is observed when the ratio \kappa = J_2 / J_1 is less than 0.203. In the region 0.203 < \kappa < 1/4, the spontaneous magnetization is discontinuous at the transition temperature. Such tricritical behavior suggests that the phase transitions on hyperbolic lattices need not always be mean-field like.Comment: 7 pages, 13 figures, submitted to Phys. Rev.

Phase transition of clock models on hyperbolic lattice studied by corner transfer matrix renormalization group method

Two-dimensional ferromagnetic N-state clock models are studied on a hyperbolic lattice represented by tessellation of pentagons. The lattice lies on the hyperbolic plane with a constant negative scalar curvature. We observe the spontaneous magnetization, the internal energy, and the specific heat at the center of sufficiently large systems, where the fixed boundary conditions are imposed, for the cases N>=3 up to N=30. The model with N=3, which is equivalent to the 3-state Potts model on the hyperbolic lattice, exhibits the first order phase transition. A mean-field like phase transition of the second order is observed for the cases N>=4. When N>=5 we observe the Schottky type specific heat below the transition temperature, where its peak hight at low temperatures scales as N^{-2}. From these facts we conclude that the phase transition of classical XY-model deep inside the hyperbolic lattices is not of the Berezinskii-Kosterlitz-Thouless type.Comment: REVTeX style, 4 pages, 6 figures, submitted to Phys. Rev.

Simulation of two dimensional quantum systems on an infinite lattice revisited: corner transfer matrix for tensor contraction

An extension of the projected entangled-pair states (PEPS) algorithm to infinite systems, known as the iPEPS algorithm, was recently proposed to compute the ground state of quantum systems on an infinite two-dimensional lattice. Here we investigate a modification of the iPEPS algorithm, where the environment is computed using the corner transfer matrix renormalization group (CTMRG) method, instead of using one-dimensional transfer matrix methods as in the original proposal. We describe a variant of the CTMRG that addresses different directions of the lattice independently, and use it combined with imaginary time evolution to compute the ground state of the two-dimensional quantum Ising model. Near criticality, the modified iPEPS algorithm is seen to provide a better estimation of the order parameter and correlators.Comment: 5 pages, 5 figures, 1 table, revised version, accepted in PR

The large system asymptotics of persistent currents in mesoscopic quantum rings

We consider a one-dimensional mesoscopic quantum ring filled with spinless electrons and threaded by a magnetic flux, which carries a persistent current at zero temperature. The interplay of Coulomb interactions and a single on-site impurity yields a non-trivial dependence of the persistent current on the size of the ring. We determine numerically the asymptotic power law for systems up to 32000 sites for various impurity strengths and compare with predictions from Bethe Ansatz solutions combined with Bosonization. The numerical results are obtained using an improved functional renormalization group (fRG) method. We apply the density matrix renormalization group (DMRG) and exact diagonalization methods to benchmark the fRG calculations. We use DMRG to study the persistent current at low electron concentrations in order to extend the validity of our results to quasi-continuous systems. We briefly comment on the quality of calculated fRG ground state energies by comparison with exact DMRG data.Comment: REVTex, 12 pages, 12 figs, accepted in Phys. Rev.