skip to main content


Title: Comb Tensor Networks
In this paper we propose a special type of a tree tensor network that has the geometry of a comb—a one-dimensional (1D) backbone with finite 1D teeth projecting out from it. This tensor network is designed to provide an effective description of higher-dimensional objects with special limited interactions or, alternatively, one-dimensional systems composed of complicated zero-dimensional objects. We provide details on the best numerical procedures for the proposed network, including an algorithm for variational optimization of the wave function as a comb tensor network and the transformation of the comb into a matrix product state. We compare the complexity of using a comb versus alternative matrix product state representations using density matrix renormalization group algorithms. As an application, we study a spin-1 Heisenberg model system which has a comb geometry. In the case where the ends of the teeth are terminated by spin-1/2 spins, we find that Haldane edge states of the teeth along the backbone form a critical spin-1/2 chain, whose properties can be tuned by the coupling constant along the backbone. By adding next-nearest-neighbor interactions along the backbone, the comb can be brought into a gapped phase with a long-range dimerization along the backbone. The critical and dimerized phases are separated by a Kosterlitz-Thouless phase transition, the presence of which we confirm numerically. Finally, we show that when the teeth contain an odd number of spins and are not terminated by spin-1/2's, a special type of comb edge states emerge.  more » « less
Award ID(s):
1812558
NSF-PAR ID:
10155591
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Physical review
Volume:
99
ISSN:
1550-235X
Page Range / eLocation ID:
235426
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. In this paper we study the critical properties of the Heisenberg spin-1/2model on a comb lattice --- a 1D backbone decorated with finite 1D chains --the teeth. We address the problem numerically by a comb tensor network thatduplicates the geometry of a lattice. We observe a fundamental difference betweenthe states on a comb with even and odd number of sites per tooth, whichresembles an even-odd effect in spin-1/2 ladders. The comb with odd teeth isalways critical, not only along the teeth, but also along the backbone, whichleads to a competition between two critical regimes in orthogonal directions.In addition, we show that in a weak-backbone limit the excitation energy scales as1/(NL), and not as 1/N or 1/L typical for 1D systems. For even teeth in theweak backbone limit the system corresponds to a collection of decoupledcritical chains of length L, while in the strong backbone limit, one spin from eachtooth forms the backbone, so the effective length of a critical toothis one site shorter, L-1. Surprisingly, these two regimes are connected via astate where a critical chain spans over two nearest neighbor teeth, with an effectivelength 2L. 
    more » « less
  2. null (Ed.)
    One dimensional (1d) interacting systems with local Hamiltonianscan be studied with various well-developed analytical methods.Recently novel 1d physics was found numerically in systems witheither spatially nonlocal interactions, or at the 1d boundary of2d quantum critical points, and the critical fluctuation in thebulk also yields effective nonlocal interactions at the boundary.This work studies the edge states at the 1d boundary of 2dstrongly interacting symmetry protected topological (SPT) states,when the bulk is driven to a disorder-order phase transition. Wewill take the 2d Affleck-Kennedy-Lieb-Tasaki (AKLT) state as anexample, which is a SPT state protected by the SO(3) spinsymmetry and spatial translation. We found that the original(1+1)d boundary conformal field theory of the AKLT state isunstable due to coupling to the boundary avatar of the bulkquantum critical fluctuations. When the bulk is fixed at thequantum critical point, within the accuracy of our expansionmethod, we find that by tuning one parameter at the boundary,there is a generic direct transition between the long rangeantiferromagnetic Néel order and the valence bond solid (VBS)order. This transition is very similar to the Néel-VBStransition recently found in numerical simulation of a spin-1/2chain with nonlocal spatial interactions. Connections between ouranalytical studies and recent numerical results concerning theedge states of the 2d AKLT-like state at a bulk quantum phasetransition will also be discussed. 
    more » « less
  3. Abstract The $p$-tensor Ising model is a one-parameter discrete exponential family for modeling dependent binary data, where the sufficient statistic is a multi-linear form of degree $p \geqslant 2$. This is a natural generalization of the matrix Ising model that provides a convenient mathematical framework for capturing, not just pairwise, but higher-order dependencies in complex relational data. In this paper, we consider the problem of estimating the natural parameter of the $p$-tensor Ising model given a single sample from the distribution on $N$ nodes. Our estimate is based on the maximum pseudolikelihood (MPL) method, which provides a computationally efficient algorithm for estimating the parameter that avoids computing the intractable partition function. We derive general conditions under which the MPL estimate is $\sqrt N$-consistent, that is, it converges to the true parameter at rate $1/\sqrt N$. Our conditions are robust enough to handle a variety of commonly used tensor Ising models, including spin glass models with random interactions and models where the rate of estimation undergoes a phase transition. In particular, this includes results on $\sqrt N$-consistency of the MPL estimate in the well-known $p$-spin Sherrington–Kirkpatrick model, spin systems on general $p$-uniform hypergraphs and Ising models on the hypergraph stochastic block model (HSBM). In fact, for the HSBM we pin down the exact location of the phase transition threshold, which is determined by the positivity of a certain mean-field variational problem, such that above this threshold the MPL estimate is $\sqrt N$-consistent, whereas below the threshold no estimator is consistent. Finally, we derive the precise fluctuations of the MPL estimate in the special case of the $p$-tensor Curie–Weiss model, which is the Ising model on the complete $p$-uniform hypergraph. An interesting consequence of our results is that the MPL estimate in the Curie–Weiss model saturates the Cramer–Rao lower bound at all points above the estimation threshold, that is, the MPL estimate incurs no loss in asymptotic statistical efficiency in the estimability regime, even though it is obtained by minimizing only an approximation of the true likelihood function for computational tractability. 
    more » « less
  4. The boundary modes of topological insulators are protected by the symmetries of the nontrivial bulk electronic states. Unless these symmetries are broken, they can give rise to novel phenomena, such as the quantum spin Hall effect in one-dimensional (1D) topological edge states, where quasiparticle backscattering is suppressed by time-reversal symmetry (TRS). Here, we investigate the properties of the 1D topological edge state of bismuth in the absence of TRS, where backscattering is predicted to occur. Using spectroscopic imaging and spin-polarized measurements with a scanning tunneling microscope, we compared quasiparticle interference (QPI) occurring in the edge state of a pristine bismuth bilayer with that occurring in the edge state of a bilayer, which is terminated by ferromagnetic iron clusters that break TRS. Our experiments on the decorated bilayer edge reveal an additional QPI branch, which can be associated with spin-flip scattering across the Brioullin zone center between time-reversal band partners. The observed QPI characteristics exactly match with theoretical expectations for a topological edge state, having one Kramer’s pair of bands. Together, our results provide further evidence for the nontrivial nature of bismuth and in particular, demonstrate backscattering inside a helical topological edge state induced by broken TRS through local magnetism.

     
    more » « less
  5. Abstract

    Transition metal dichalcogenide materials are unique in the wide variety of structural and electronic phases they exhibit in the two-dimensional limit. Here we show how such polymorphic flexibility can be used to achieve topological states at highly ordered phase boundaries in a new quantum spin Hall insulator (QSHI), 1T′-WSe2. We observe edge states at the crystallographically aligned interface between a quantum spin Hall insulating domain of 1T′-WSe2and a semiconducting domain of 1H-WSe2in contiguous single layers. The QSHI nature of single-layer 1T′-WSe2is verified using angle-resolved photoemission spectroscopy to determine band inversion around a 120 meV energy gap, as well as scanning tunneling spectroscopy to directly image edge-state formation. Using this edge-state geometry we confirm the predicted penetration depth of one-dimensional interface states into the two-dimensional bulk of a QSHI for a well-specified crystallographic direction. These interfaces create opportunities for testing predictions of the microscopic behavior of topologically protected boundary states.

     
    more » « less