Two-dimensional Quantum Lattice Research Articles

Two-dimensional quantum lattice models via mode optimized hybrid CPU-GPU density matrix renormalization group method

Two-dimensional quantum lattice models via mode optimized hybrid CPU-GPU density matrix renormalization group method

Dissipative quantum many-body dynamics in (1+1)D quantum cellular automata and quantum neural networks

Classical artificial neural networks, built from elementary units, possess enormous expressive power. Here we investigate a quantum neural network (QNN) architecture, which follows a similar paradigm. It is structurally equivalent to so-called (1+1)D quantum cellular automata, which are two-dimensional quantum lattice systems on which dynamics takes place in discrete time. Information transfer between consecutive time slices—or adjacent network layers—is governed by local quantum gates, which can be regarded as the quantum counterpart of the classical elementary units. Along the time-direction an effective dissipative evolution emerges on the level of the reduced state, and the nature of this dynamics is dictated by the structure of the elementary gates. We show how to construct the local unitary gates to yield a desired many-body dynamics, which in certain parameter regimes is governed by a Lindblad master equation. We study this for small system sizes through numerical simulations and demonstrate how collective effects within the quantum cellular automaton can be controlled parametrically. Our study constitutes a step towards the utilization of large-scale emergent phenomena in large QNNs for machine learning purposes.

TeNeS: Tensor network solver for quantum lattice systems

TeNeS (Tensor Network Solver) [1,2] is a free/libre open-source software program package for calculating two-dimensional many-body quantum states based on the tensor network method and the corner transfer matrix renormalization group (CTMRG) method. This package calculates ground-state wavefunctions for user-defined Hamiltonians and evaluates user-defined physical quantities such as magnetization and correlation functions. For certain predefined models and lattices, there is a tool that makes it easy to generate input files. TeNeS uses an OpenMP/MPI hybrid parallelized tensor operation library and thus can perform large-scale calculations using massively parallel machines. Program summaryProgram Title: TeNeSCPC Library link to program files:https://doi.org/10.17632/psm26xxbvd.1Code Ocean capsule:https://codeocean.com/capsule/1510058Licensing provisions: GNU General Public License version 3Programming language:C++ and python3External routines/libraries: mptensorNature of problem: TeNeS calculates the approximate ground states and their properties of user-defined two-dimensional quantum lattice models using user-friendly input files. Numerically exact solutions of such tasks generally require an exponentially diverging computational time, whereas the error in the output of TeNeS is well controlled and can be reduced with a polynomial cost.Solution method: TeNeS implements the tensor networks method based on a tensor-product-state (TPS) wavefunction and the corner transfer matrix renormalization group method. TeNeS also supports massively parallel computing using the library mptensor, which implements OpenMP/MPI hybrid parallelized tensor operations.

One- and two-dimensional quantum lattice algorithms for Maxwell equations in inhomogeneous scalar dielectric media I: theory

ABSTRACT A quantum lattice algorithm (QLA) is developed for Maxwell equations in scalar dielectric media using the Riemann–Silberstein representation on a Cartesian grid. For x-dependent and y-dependent dielectric inhomogeneities, the corresponding QLA requires a minimum of 8 qubits/spatial lattice site. This is because the corresponding Pauli spin matrices have off-diagonal components which permit the local collisional entanglement of these qubits. However, z-dependent inhomogeneities require a QLA with a minimum of 16 qubits/lattice site since the Pauli spin matrix is diagonal. For two-dimensional inhomogeneities, one can readily couple the 8–8 qubit schemes for x−y variations. z−x and y−z variations can be treated by either a 16–8 qubit scheme or a 16–16 qubit representation.

One- and two-dimensional quantum lattice algorithms for Maxwell equations in inhomogeneous scalar dielectric media. II: Simulations

Long-time quantum lattice algorithm (QLA) simulations are performed for the multiple reflection–transmission of an initial electromagnetic pulse propagating normally to a boundary layer region joining two media of different refractive index. For these one-dimensional (1D) simulations, there is excellent agreement between x-, y- and z-representations, as well as very good agreement with nearly all the standard plane wave boundary condition results for reflection and transmission off a dielectric discontinuity. In the QLA simulation, no boundary conditions are imposed at the continuous, but sharply increasing, dielectric boundary layers. Two-dimensional (2D) QLA scattering simulations in the x–z plane are performed for an electromagnetic pulse interacting with a conical dielectric obstacle for the 8–16 qubit model.

An algebraic approach to discrete time integrability

We propose the systematic construction of classical and quantum two-dimensional space-time lattices primarily based on algebraic considerations, i.e. on the existence of associated r-matrices and underlying spatial and temporal classical and quantum algebras. This is a novel construction that leads to the derivation of fully discrete integrable systems governed by sets of consistent integrable non-linear space-time difference equations. To illustrate the proposed methodology, we derive two versions of the fully discrete non-linear Schrödinger type system. The first one is based on the existence of a rational r-matrix, whereas the second one is the fully discrete Ablowitz–Ladik model and is associated to a trigonometric r-matrix. The Darboux-dressing method is also applied for the first discretization scheme, mostly as a consistency check, and solitonic as well as general solutions, in terms of solutions of the fully discrete heat equation, are also derived. The quantization of the fully discrete systems is then quite natural in this context and the two-dimensional quantum lattice is thus also examined.

Simulating excitation spectra with projected entangled-pair states

We develop and benchmark a technique for simulating excitation spectra of generic two-dimensional quantum lattice systems using the framework of projected entangled-pair states (PEPS). The technique relies on a variational ansatz for capturing quasiparticle excitations on top of a PEPS ground state. Our method perfectly captures the quasiparticle dispersion relation of the square-lattice transverse-field Ising model, and reproduces the spin-wave velocity and the spin-wave anomaly in the square-lattice Heisenberg model with high precision.

Tensor Network Annealing Algorithm for Two-Dimensional Thermal States.

Tensor network methods have become a powerful class of tools to capture strongly correlated matter, but methods to capture the experimentally ubiquitous family of models at finite temperature beyond one spatial dimension are largely lacking. We introduce a tensor network algorithm able to simulate thermal states of two-dimensional quantum lattice systems in the thermodynamic limit. The method develops instances of projected entangled pair states and projected entangled pair operators for this purpose. It is the key feature of this algorithm to resemble the cooling down of the system from an infinite temperature state until it reaches the desired finite-temperature regime. As a benchmark, we study the finite-temperature phase transition of the Ising model on an infinite square lattice, for which we obtain remarkable agreement with the exact solution. We then turn to study the finite-temperature Bose-Hubbard model in the limits of two (hard-core) and three bosonic modes per site. Our technique can be used to support the experimental study of actual effectively two-dimensional materials in the laboratory, as well as to benchmark optical lattice quantum simulators with ultracold atoms.

Nature of continuous phase transitions in interacting topological insulators

We revisit the effects of the Hubbard repulsion on quantum spin Hall effects (QSHE) in two-dimensional quantum lattice models. We present both unbiased exact diagonalization and density-matrix renormalization group simulations with numerical evidences for a continuous quantum phase transition (CQPT) separating QSHE from the topologically trivial antiferromagnetic phase. Our numerical results suggest that, the nature of CQPT exhibits distinct finite-size scaling behaviors, which may be consistent with either Ising or XY universality classes for different time-reversal symmetric QSHE systems.

Decay of Correlations in 2D Quantum Systems with Continuous Symmetry

We study a large class of models of two-dimensional quantum lattice systems with continuous symmetries, and we prove a general McBryan–Spencer–Koma–Tasaki theorem concerning algebraic decay of correlations. We present applications of our main result to the Heisenberg, Hubbard, and t-J models, and to certain models of random loops.

Geometric entanglement and quantum phase transitions in two-dimensional quantum lattice models

Geometric entanglement(GE), as a measure of multipartite entanglement, has been investigated as a universal tool to detect phase transitions in quantum many-body lattice models. We outline a systematic method to compute GE for two-dimensional (2D) quantum many-body lattice models based on the translational invariant structure of infinite projected entangled pair state (iPEPS) representations. By employing this method, the $q$-state quantum Potts model on the square lattice with $q \in \{2, 3 ,4, 5\}$ is investigated as a prototypical example. Further, we have explored three 2D Heisenberg models, such as the spin-$\frac{1}{2}$ XXX and antiferromagnetic anisotropic XYX models in an external magnetic field, and the spin-1 antiferromagnetic XXZ model. We find that continuous GE does not guarantee a continuous phase transition across a phase transition point. We observe and thus classify three different types of continuous GE across a phase transition point: (i) GE is continuous with maximum value at the transition point and the phase transition is continuous, (ii) GE is continuous with maximum value at the transition point but the phase transition is discontinuous, and (iii) GE is continuous with non-maximum value at the transition point and the phase transition is continuous. For the models under consideration we find that the second and the third types are related to a point of dual symmetry and a fully polarized phase, respectively.

Theory of network contractor dynamics for exploring thermodynamic properties of two-dimensional quantum lattice models

Based on the tensor network state representation, we develop a nonlinear dynamic theory, coined network contractor dynamics (NCD), to explore the thermodynamic properties of two-dimensional quantum lattice models. By invoking the rank-1 decomposition in the multilinear algebra, the NCD scheme makes the contraction of the tensor network of the partition function be realized through a contraction of a local tensor cluster with vectors on its boundary. An imaginary-time-sweep algorithm for implementation of the NCD method is proposed for practical numerical simulations. We benchmark the NCD scheme on the square Ising model, which shows great accuracy. Also, the results on the spin-1/2 Heisenberg antiferromagnet on a honeycomb lattice are disclosed to be in good agreement with the quantum Monte Carlo calculations. The quasientanglement entropy $S$, Lyapunov exponent ${I}^{\text{Lya}}$, and loop character ${I}^{\text{loop}}$ are introduced within the dynamic scheme, which are found to display ``nonlocality'' near the critical point, and can be applied to determine the thermodynamic phase transitions of both classical and quantum systems.

On the quantum L-operator for the two-dimensional lattice Toda model

We consider the two-dimensional quantum lattice Toda model for affine and simple Lie algebras of type A. For its known L-operator, the second-order correction in lattice parameter e is found. It is proved that the equation determining the third-order correction in e has no solutions. Bibliography: 9 titles.

Optimized decimation of tensor networks with super-orthogonalization for two-dimensional quantum lattice models

A novel algorithm based on the optimized decimation of tensor networks with super-orthogonalization (ODTNS) that can be applied to simulate efficiently and accurately not only the thermodynamic but also the ground state properties of two-dimensional (2D) quantum lattice models is proposed. By transforming the 2D quantum model into a three-dimensional (3D) closed tensor network (TN) comprised of the tensor product density operator and a 3D brick-wall TN, the free energy of the system can be calculated with the imaginary time evolution, in which the network Tucker decomposition is suggested for the first time to obtain the optimal lower-dimensional approximation on the bond space by transforming the TN into a super-orthogonal form. The efficiency and accuracy of this algorithm are testified, which are fairly comparable with the quantum Monte Carlo calculations. Besides, the present ODTNS scheme can also be applicable to the 2D frustrated quantum spin models with nice efficiency.

Studying Two-Dimensional Systems with the Density Matrix Renormalization Group

The density matrix renormalization group (DMRG) is one of the most powerful numerical methods for studying two-dimensional quantum lattice systems, despite a perception that it is only suitable for one dimension. Reviewing past applications of DMRG in 2D demonstrates its success in treating a wide variety of problems, although it remains underutilized in this context. We present techniques for performing cutting-edge 2D DMRG studies including methods for ensuring convergence, extrapolating finite-size data, and extracting gaps and excited states. Finally, we consider a selection of recently developed 2D tensor network methods and compare the performance of one of these to 2D DMRG.

Capturing long range correlations in two-dimensional quantum lattice systems using correlator product states

We study the suitability of correlator product states for describing ground-state properties of two-dimensional spin models. Our ansatz for the many-body wave function takes the form of either plaquette or bond correlator product states and the energy is optimized by varying the correlators using Monte Carlo minimization. For the Ising model we find that plaquette correlators are best for estimating the energy while bond correlators capture the expected long-range correlations and critical behavior of the system more faithfully. For the antiferromagnetic Heisenberg model, however, plaquettes outperform bond correlators at describing both local and long-range correlations because of the substantially larger number of local parameters they contain. These observations have quantitative implications for the application of correlator product states to other more complex systems, and give important heuristic insights: in particular the necessity of carefully tailoring the choice of correlators to the system considered, its interactions and symmetries.

Direct Extension of Density-Matrix Renormalization Group to Two-Dimensional Quantum Lattice Systems: Studies of Parallel Algorithm, Accuracy, and Performance

We parallelize the density-matrix renormalization group method to directly extend it to two-dimensional ( n -leg) quantum lattice models. The parallelization is performed mainly on the diagonalization for the superblock Hamiltonian since this step requires enormous memory space as the leg number n increases. The superblock Hamiltonian is divided into three parts, and the corresponding superblock vectors are transformed into matrices whose elements are uniformly distributed into processors. Parallel efficiency increases as the number of states kept, m , increases, and the obtained ground-state energy rapidly converges within a few sweeps. Furthermore, the present algorithm applied to doped 4-leg Hubbard ladders reaches their ground states, which satisfy the Lieb–Mattis theorem with m much smaller than those used in different indirect algorithms.

Entanglement Renormalization in Two Spatial Dimensions

We propose and test a scheme for entanglement renormalization capable of addressing large two-dimensional quantum lattice systems. In a translationally invariant system, the cost of simulations grows only as the logarithm of the lattice size; at a quantum critical point, the simulation cost becomes independent of the lattice size and infinite systems can be analyzed. We demonstrate the performance of the scheme by investigating the low energy properties of the 2D quantum Ising model on a square lattice of linear size L={6, 9, 18, 54, infinity} with periodic boundary conditions. We compute the ground state and evaluate local observables and two-point correlators. We also produce accurate estimates of the critical magnetic field and critical exponent beta. A calculation of the energy gap shows that it scales as 1/L at the critical point.

Classical Simulation of Infinite-Size Quantum Lattice Systems in Two Spatial Dimensions

We present an algorithm to simulate two-dimensional quantum lattice systems in the thermodynamic limit. Our approach builds on the projected entangled-pair state algorithm for finite lattice systems [F. Verstraete and J. I. Cirac, arxiv:cond-mat/0407066] and the infinite time-evolving block decimation algorithm for infinite one-dimensional lattice systems [G. Vidal, Phys. Rev. Lett. 98, 070201 (2007)10.1103/PhysRevLett.98.070201]. The present algorithm allows for the computation of the ground state and the simulation of time evolution in infinite two-dimensional systems that are invariant under translations. We demonstrate its performance by obtaining the ground state of the quantum Ising model and analyzing its second order quantum phase transition.

Topological phase separation in a two-dimensional quantum lattice Bose-Hubbard system away from half-filling

We suppose that the doping of the two-dimensional (2D) hard-core (hc) boson system away from half-filling may result in the formation of a multicenter topological inhomogeneity (defect), such as charge order (CO) bubble domain(s) with Bose superfluid (BS) and extra bosons both localized in domain wall(s), or else in a topological $\text{CO}+\text{BS}$ phase separation rather than a uniform mixed $\text{CO}+\text{BS}$ supersolid phase. Starting from the classical model, we predict the properties of the respective quantum system. The long-wavelength behavior of the system is believed to be reminiscent of that of granular superconductors, charge-density wave materials, Wigner crystals, and multi-Skyrmion systems akin to the quantum Hall ferromagnetic state of a 2D electron gas. To elucidate the role played by quantum effects and that of the lattice discreteness, we have addressed the simplest nanoscopic counterpart of the bubble domain in the checkerboard CO phase of a 2D hc Bose-Hubbard (BH) square lattice. It is shown that the relative magnitude and symmetry of a multicomponent order parameter are mainly determined by the sign of the nearest-neighbor and next-nearest-neighbor transfer integrals. In general, the topologically inhomogeneous phase of the hc-BH system away from the half-filling can exhibit the signatures of the $s$, $d$, and $p$ symmetries of the off-diagonal order.