Plaquette Renormalized Tensor Network States: Application to Frustrated Systems
|
|
- Oswald Lester
- 5 years ago
- Views:
Transcription
1 Workshop on QIS and QMP, Dec 20, 2009 Plaquette Renormalized Tensor Network States: Application to Frustrated Systems Ying-Jer Kao and Center for Quantum Science and Engineering Hsin-Chih Hsiao, Ji-Feng Yu(NTU) Anders W. Sandvik, Ling Wang (Boston)
2 Frustration AF or Kagome Pyrochlore Hyperkagome Garnet ZnCu3(OH)6Cl2 Tb2Ti2O7 Na4Ir3O8 Gd3Ga5O12 Large number of degenerate classical ground states Emergence of novel spin-disordered ground states due to quantum fluctuations Hard to study numerically: size limitation in Exact Diagonalization, sign problem in QMC (also fermion), dimension limit for DMRG (d=1)
3 Negative Sign Problem H = ij H ij = ij(s + i S j + S i S+ j ) flip α flip flip α p(c) = ( 1) 3 α H 12 H 23 H 31 α < 0 A p(c) = As p(c) s p(c) Negative matrix elements Fluctuations in signs:ensemble average diverges at low temperature
4 Simulation of Quantum Systems Is it possible to classically simulate faithfully a quantum system? To represent a quantum state: ψ = classical representation requires d n complex coefficients which grows exponentially with n Ground states of quantum many-body systems are usually entangled. d i 1,i 2,...,i n =1 c i1 i 2...i n i 1 i 2... i n
5 Entanglement Complementary viewpoints on entanglement: Quantum information theory: a crucial resource to process and send information in novel ways Quantum many-body physics: entanglement gives rise to exotic phases of matter Numerical simulation of strongly correlated quantum systems: source of difficulties!! What kind of superpositions appear in nature? symmetries, local interactions, little entanglement
6 Entanglement: Matrix Product State i α A [1]i 1 α 1 α 2 A [2]i 2 α 2 α 3 A [3]i 3 α 3 α 4 A [4]i 4 α 4 α 5 A [5]i 5 α 5 α 6 A [6]i 6 α 6 α 7 A [7]i 7 α 7 α 8 Between each nearest-neighbor lattice site, we introduce a D- dimensional maximally entangled state: φ >= 1 D αα > D α=1 At each lattice site there are two ends to these maximally entangled states Operator A: project D 2 to d (AKLT GS: D=2, d=3) ψ = d i 1,...,i N =1 A 1 (i 1 )A 2 (i 2 ) A N (i N ) i 1, i 2,..., i N We reexpress the 2 N coefficients of in terms of about 2D 2 N parameters (linear in N) Key ingredient behind the success of DMRG S. Ostlund and S. Rommer, Phys. Rev. Lett. 75, 3537(1995) Afflect, Kennedy, Lieb, Tasaki, Commun Math Phys 115 (1988)
7 Tensor product states Extension similar ideas to higher dimension: tensor product states, projected entangled paired state (PEPS) ψ = {s} ttr{a(s 1 )A(s 2 ) A(s N )} s 1, s 2,..., s N s =±1 i j A s k l =1...D Aijkl: rank-4 tensor A s 1 A s 4 A s 7 A s 2 A s 3 A s 5 A s6 A s 8 A s9 F. Verstraete and J. I. Cirac, arxiv:cond-mat/
8 Variational wave function Use tensor product state as trial wave function E = ψ H ψ ψ ψ ψ ψ = E 0 H = s 1,s 2,... i,j,k,l,... = ttr [T T T...] i H i H i = Ô0 i + Ô1 i Ô2 j +... A s 1 i,j,k,l A s 1 i,j,k,l As 2 k,m,n,o A s 1 k,m,n,o Double tensor T = s,s A s A s ) T
9 Variational Energy H = i H i H i = Ô0 i + Ô1 i Ô2 j +... ψ H ψ = ttr Ti 0 T T... +ttr Ti 1 Tj 2 T T a = s,s A s A s s Ôa s Computationally intensive: D N_bond Direct computation of contraction is impossible: approximation (RG schemes)
10 Tensor Contraction Contracting the internal indices, the four-leg tensor can be viewed as a single tensor. External link dimension becomes D 2 after one contraction; exponential growth as we keep contracting D 4,D 8... Computationally intensive; Impossible to store intermediate results. Need some RG scheme.
11 Tensor entanglement renormalization Singular Value Decomposition: rank-3 tensors S T αβµν T αβµν D cut γ =1 D cut γ =1 S 1 µνγ S3 αβγ S 2 ανγ S4 µβγ Introduce cut-off for the bond-contraction T γσλρ αβµν S 1 βαγs 2 µβσs 3 νµγs 4 ανρ H. C. Jiang, Z. Y. Weng, and T. Xiang, Phys. Rev. Lett. 101, (2008) Gu, Levin, Wen, Phys. Rev. B, 78, (2008)
12 Tranverse ising model H = ij σ z i σ z j h i σ x i D=2, Dcut=18, non-mf results is obtained Method is not variational Locally optimized SVD truncation Translational invariance explicitly broken Gu, Levin, Wen, Phys. Rev. B, 78, (2008)
13 Distribution of Singular Values Transverse Ising Model, 8x8, h=3.08 Fully Frustrated XX model, 8x8, h=1 i i i Frustrated Systems Remain Difficult i
14 Plaquette-Renormalization 8-index tensor: D 8 4-index tensor: Dcut 4 8-index dbl-tensor: D 16 4-index dbl-tensor: Dcut 8 e f d g c h b a k j l i d l d r e l e r c l c r b l b r k i f l a l f r a r g l h l l g r h r j a,b,...g,h=1...d T 0 T 0 i,j,k,l=1...dcut S 1 S 1 a,b,...g,h=1...d i,j,k,l=1...dcut 2 Renormalization of an 8-index plaquette tensor using auxiliary 3-index tensors S. T 0 T 0 T 0 T 0 S 1 T k 1 i S 1 T 0 T 0 j l Wang, Kao, Sandvik, arxiv:
15 Plaquette renormalized trial wave function We treat the elements in the A and S tensors as variational parameters, and minimize the total energy. Principal axis method: derivative-free, but computationally expensive.
16 Expectation values
17 Plaquette-Renormalization of TNS Effective reduced tensor network for a 8x8 lattice Summing over all unequivalent bonds and sites Method is variational Optimize T and R globally Size scaling: L 2 Log(L) R 1 T 1 T 0 T 2 R 2 Wang, Kao, Sandvik, arxiv:
18 Transverse Ising Model m z m x D=Dcut= Assume translational invariance: all initial T s are the same. Globally optimized T and R. hc=3.33 (3.04,QMC) mz~ (h-hc) β, β~0.40 h~hc, β~0.50 mean-field like. (Sanvik s talk) h/j Wang, Kao, Sandvik, arxiv:
19 Transverse Ising Model Spins at different lattice sites inside the plaquette have different environments. We use different tensors inside a plaquette.
20 J1-J2 Heisenberg Model H = J 1 S i S j + J 2 S i S j ij ij H. J. Schulz, T. A. Ziman and D. Poilblanc, J. Phys. I 6 (1996)
21 Error in energy Improvement in accuracy as D increase. D also sets a length scale. J2=0 A. W. Sandvik, Phys. Rev. B 56, (1997)
22 Building block: plaquette α 12 internal sums. Maximum computational effort: D 8. Each free index and summation contributes D. γ γ1 γ2 S 1 T 0 T 0 l α1 S 1 i j α2 T 0 k T 0 S 1 δ1 δ2 δ D 7 D 7 D 8 D 7 D 6 β2 S 1 β β1 T (1),αβ γδ = α 1,α 2,...,δ 1,δ 2,ijkl S (1),α α 1 α 2 S (1),β β 1 β 2 S (1),γ γ 1 γ 2 S (1),δ δ 1 δ 2 T (0),α 1l γ 2 i T (0),α 2j iδ 1 T (0),jβ 1 kδ 2 T (0),lβ 2 γ 1 k
23 Computational Costs Globally optimized T and R. Contraction calculation is highly parallelizable. IBM Blue Gene/L at BU. D=2, takes weeks to optimize. Bottleneck: plaquette internal contractions. Take advantage of GPU (Graphic Processing Unit). Wang, Kao, Sandvik, arxiv:
24 GPU supercomputing The GPU is specialized for compute-intensive, massively data parallel computation 1.4 billion transistors 1 Teraflop of processing power 240 processing cores nvidia Tesla 10-Series Processor NTU CQSE GPU cluster 16x Tesla S1070=64 Teraflops
25 S 1 GPU kernel Strategy Inputs: T (D 4 ), S (D 3 ) (unchanged during contraction) S 1 T 0 T 0 T 0 T 0 S 1 Load into texture memory (optimized for 2D read) Output: T (D 4 ) block geometry up to 3D (tx,ty,tz), grid of blocks geometry up to 2D (bx,by). Reshape T into a D 2 xd 2 matrix, S into a DxD 2 matrix. Use D 2 xd 2 thread blocks, each block has DxD threads. Each block computes an element of T. Use block indices (bx,by) to represent external indices. Host Kernel 1 Kernel 2 Block (1, 1) (0, 0) (0, 1) (0, 2) Device (1, 0) (1, 1) (1, 2) S 1 Grid 1 Block (0, 0) Block (0, 1) Grid 2 (2, 0) (2, 1) (2, 2) Block (1, 0) Block (1, 1) (3, 0) (3, 1) (3, 2) (4, 0) (4, 1) (4, 2) Block (2, 0) Block (2, 1) NVIDIA
26 Benchmark: One plaquette contraction CPU code: Intel Core 2 Duo E GHz D=6, 2.1 GFLOPS D=8, 3.4 GFLOPS GPU, Single Precision code (unoptimized): D=6, 24.9 GFLOPS (GTX8800), 49.7GFLOPS (GTX280) D=8, 92.2 GFLOPS (GTX280) 10x to 30x speedup. A factor of at least 8 performance hit for DP. GPU code requires optimization Bank conflict in shared memory R/W Low stream multiprocessor occupancy Current scheme D bound by shared memory size Global memory bandwidth hiding
27 Conclusion Tensor network states are promising candidates to understand frustrated quantum spin systems. In plaquette renormalized tensor network representation, no approximations are made when contracting the effective renormalized tensor network. Non-MFT results even with the smallest possible nontrivial tensors and truncation (D = 2) Larger internal bond dimension D is necessary to get the right physics. GPU can potentially speed up the computationally intensive part of the calculation. Use Monte Carlo sampling in the future.
Matrix-Product states: Properties and Extensions
New Development of Numerical Simulations in Low-Dimensional Quantum Systems: From Density Matrix Renormalization Group to Tensor Network Formulations October 27-29, 2010, Yukawa Institute for Theoretical
More informationQuantum many-body systems and tensor networks: simulation methods and applications
Quantum many-body systems and tensor networks: simulation methods and applications Román Orús School of Physical Sciences, University of Queensland, Brisbane (Australia) Department of Physics and Astronomy,
More informationTensor network methods in condensed matter physics. ISSP, University of Tokyo, Tsuyoshi Okubo
Tensor network methods in condensed matter physics ISSP, University of Tokyo, Tsuyoshi Okubo Contents Possible target of tensor network methods! Tensor network methods! Tensor network states as ground
More informationRenormalization of Tensor Network States
Renormalization of Tensor Network States I. Coarse Graining Tensor Renormalization Tao Xiang Institute of Physics Chinese Academy of Sciences txiang@iphy.ac.cn Numerical Renormalization Group brief introduction
More informationRenormalization of Tensor- Network States Tao Xiang
Renormalization of Tensor- Network States Tao Xiang Institute of Physics/Institute of Theoretical Physics Chinese Academy of Sciences txiang@iphy.ac.cn Physical Background: characteristic energy scales
More informationQuantum simulation with string-bond states: Joining PEPS and Monte Carlo
Quantum simulation with string-bond states: Joining PEPS and Monte Carlo N. Schuch 1, A. Sfondrini 1,2, F. Mezzacapo 1, J. Cerrillo 1,3, M. Wolf 1,4, F. Verstraete 5, I. Cirac 1 1 Max-Planck-Institute
More informationTime Evolving Block Decimation Algorithm
Time Evolving Block Decimation Algorithm Application to bosons on a lattice Jakub Zakrzewski Marian Smoluchowski Institute of Physics and Mark Kac Complex Systems Research Center, Jagiellonian University,
More informationAdvanced Computation for Complex Materials
Advanced Computation for Complex Materials Computational Progress is brainpower limited, not machine limited Algorithms Physics Major progress in algorithms Quantum Monte Carlo Density Matrix Renormalization
More informationClassical Monte Carlo Simulations
Classical Monte Carlo Simulations Hyejin Ju April 17, 2012 1 Introduction Why do we need numerics? One of the main goals of condensed matter is to compute expectation values O = 1 Z Tr{O e βĥ} (1) and
More informationScaling analysis of snapshot spectra in the world-line quantum Monte Carlo for the transverse-field Ising chain
TNSAA 2018-2019 Dec. 3-6, 2018, Kobe, Japan Scaling analysis of snapshot spectra in the world-line quantum Monte Carlo for the transverse-field Ising chain Kouichi Seki, Kouichi Okunishi Niigata University,
More informationParamagnetic phases of Kagome lattice quantum Ising models p.1/16
Paramagnetic phases of Kagome lattice quantum Ising models Predrag Nikolić In collaboration with T. Senthil Massachusetts Institute of Technology Paramagnetic phases of Kagome lattice quantum Ising models
More informationQuantum Hamiltonian Complexity. Itai Arad
1 18 / Quantum Hamiltonian Complexity Itai Arad Centre of Quantum Technologies National University of Singapore QIP 2015 2 18 / Quantum Hamiltonian Complexity condensed matter physics QHC complexity theory
More informationLoop optimization for tensor network renormalization
Yukawa Institute for Theoretical Physics, Kyoto University, Japan June, 6 Loop optimization for tensor network renormalization Shuo Yang! Perimeter Institute for Theoretical Physics, Waterloo, Canada Zheng-Cheng
More informationThe density matrix renormalization group and tensor network methods
The density matrix renormalization group and tensor network methods Outline Steve White Exploiting the low entanglement of ground states Matrix product states and DMRG 1D 2D Tensor network states Some
More informationIntroduction to numerical computations on the GPU
Introduction to numerical computations on the GPU Lucian Covaci http://lucian.covaci.org/cuda.pdf Tuesday 1 November 11 1 2 Outline: NVIDIA Tesla and Geforce video cards: architecture CUDA - C: programming
More informationarxiv: v1 [hep-lat] 7 Oct 2010
arxiv:.486v [hep-lat] 7 Oct 2 Nuno Cardoso CFTP, Instituto Superior Técnico E-mail: nunocardoso@cftp.ist.utl.pt Pedro Bicudo CFTP, Instituto Superior Técnico E-mail: bicudo@ist.utl.pt We discuss the CUDA
More informationIntroduction to tensor network state -- concept and algorithm. Z. Y. Xie ( 谢志远 ) ITP, Beijing
Introduction to tensor network state -- concept and algorithm Z. Y. Xie ( 谢志远 ) 2018.10.29 ITP, Beijing Outline Illusion of complexity of Hilbert space Matrix product state (MPS) as lowly-entangled state
More informationMachine Learning with Tensor Networks
Machine Learning with Tensor Networks E.M. Stoudenmire and David J. Schwab Advances in Neural Information Processing 29 arxiv:1605.05775 Beijing Jun 2017 Machine learning has physics in its DNA # " # #
More informationTensor network renormalization
Coogee'15 Sydney Quantum Information Theory Workshop Tensor network renormalization Guifre Vidal In collaboration with GLEN EVENBLY IQIM Caltech UC Irvine Quantum Mechanics 1920-1930 Niels Bohr Albert
More informationIPAM/UCLA, Sat 24 th Jan Numerical Approaches to Quantum Many-Body Systems. QS2009 tutorials. lecture: Tensor Networks.
IPAM/UCLA, Sat 24 th Jan 2009 umerical Approaches to Quantum Many-Body Systems QS2009 tutorials lecture: Tensor etworks Guifre Vidal Outline Tensor etworks Computation of expected values Optimization of
More informationTwo case studies of Monte Carlo simulation on GPU
Two case studies of Monte Carlo simulation on GPU National Institute for Computational Sciences University of Tennessee Seminar series on HPC, Feb. 27, 2014 Outline 1 Introduction 2 Discrete energy lattice
More informationQuantum Monte Carlo Simulations in the Valence Bond Basis
NUMERICAL APPROACHES TO QUANTUM MANY-BODY SYSTEMS, IPAM, January 29, 2009 Quantum Monte Carlo Simulations in the Valence Bond Basis Anders W. Sandvik, Boston University Collaborators Kevin Beach (U. of
More informationQuantum spin systems - models and computational methods
Summer School on Computational Statistical Physics August 4-11, 2010, NCCU, Taipei, Taiwan Quantum spin systems - models and computational methods Anders W. Sandvik, Boston University Lecture outline Introduction
More informationGPU Computing Activities in KISTI
International Advanced Research Workshop on High Performance Computing, Grids and Clouds 2010 June 21~June 25 2010, Cetraro, Italy HPC Infrastructure and GPU Computing Activities in KISTI Hongsuk Yi hsyi@kisti.re.kr
More informationarxiv: v1 [cond-mat.str-el] 1 Apr 2014
Replica exchange molecular dynamics optimization of tensor networ states for quantum many-body systems Wenyuan Liu, 1,2 Chao Wang, 1,2 Yanbin Li, 1 Yuyang Lao, 1 Yongjian Han, 1,2, Guang-Can Guo, 1,2 and
More informationQuantum Information and Quantum Many-body Systems
Quantum Information and Quantum Many-body Systems Lecture 1 Norbert Schuch California Institute of Technology Institute for Quantum Information Quantum Information and Quantum Many-Body Systems Aim: Understand
More informationRenormalization of Tensor Network States. Partial order and finite temperature phase transition in the Potts model on irregular lattice
Renormalization of Tensor Network States Partial order and finite temperature phase transition in the Potts model on irregular lattice Tao Xiang Institute of Physics Chinese Academy of Sciences Background:
More informationScale invariance on the lattice
Coogee'16 Sydney Quantum Information Theory Workshop Feb 2 nd - 5 th, 2016 Scale invariance on the lattice Guifre Vidal Coogee'16 Sydney Quantum Information Theory Workshop Feb 2 nd - 5 th, 2016 Scale
More informationPopulation annealing study of the frustrated Ising antiferromagnet on the stacked triangular lattice
Population annealing study of the frustrated Ising antiferromagnet on the stacked triangular lattice Michal Borovský Department of Theoretical Physics and Astrophysics, University of P. J. Šafárik in Košice,
More informationJournal Club: Brief Introduction to Tensor Network
Journal Club: Brief Introduction to Tensor Network Wei-Han Hsiao a a The University of Chicago E-mail: weihanhsiao@uchicago.edu Abstract: This note summarizes the talk given on March 8th 2016 which was
More informationDirect Self-Consistent Field Computations on GPU Clusters
Direct Self-Consistent Field Computations on GPU Clusters Guochun Shi, Volodymyr Kindratenko National Center for Supercomputing Applications University of Illinois at UrbanaChampaign Ivan Ufimtsev, Todd
More informationNeural Network Representation of Tensor Network and Chiral States
Neural Network Representation of Tensor Network and Chiral States Yichen Huang ( 黄溢辰 ) 1 and Joel E. Moore 2 1 Institute for Quantum Information and Matter California Institute of Technology 2 Department
More informationFermionic tensor networks
Fermionic tensor networks Philippe Corboz, Institute for Theoretical Physics, ETH Zurich Bosons vs Fermions P. Corboz and G. Vidal, Phys. Rev. B 80, 165129 (2009) : fermionic 2D MERA P. Corboz, R. Orus,
More informationQuantum Lattice Models & Introduction to Exact Diagonalization
Quantum Lattice Models & Introduction to Exact Diagonalization H! = E! Andreas Läuchli IRRMA EPF Lausanne ALPS User Workshop CSCS Manno, 28/9/2004 Outline of this lecture: Quantum Lattice Models Lattices
More informationGPU accelerated Monte Carlo simulations of lattice spin models
Available online at www.sciencedirect.com Physics Procedia 15 (2011) 92 96 GPU accelerated Monte Carlo simulations of lattice spin models M. Weigel, T. Yavors kii Institut für Physik, Johannes Gutenberg-Universität
More informationEfficient time evolution of one-dimensional quantum systems
Efficient time evolution of one-dimensional quantum systems Frank Pollmann Max-Planck-Institut für komplexer Systeme, Dresden, Germany Sep. 5, 2012 Hsinchu Problems we will address... Finding ground states
More informationFermionic topological quantum states as tensor networks
order in topological quantum states as Jens Eisert, Freie Universität Berlin Joint work with Carolin Wille and Oliver Buerschaper Symmetry, topology, and quantum phases of matter: From to physical realizations,
More informationMachine Learning with Quantum-Inspired Tensor Networks
Machine Learning with Quantum-Inspired Tensor Networks E.M. Stoudenmire and David J. Schwab Advances in Neural Information Processing 29 arxiv:1605.05775 RIKEN AICS - Mar 2017 Collaboration with David
More informationThe Density Matrix Renormalization Group: Introduction and Overview
The Density Matrix Renormalization Group: Introduction and Overview Introduction to DMRG as a low entanglement approximation Entanglement Matrix Product States Minimizing the energy and DMRG sweeping The
More information3 Symmetry Protected Topological Phase
Physics 3b Lecture 16 Caltech, 05/30/18 3 Symmetry Protected Topological Phase 3.1 Breakdown of noninteracting SPT phases with interaction Building on our previous discussion of the Majorana chain and
More informationT ensor N et works. I ztok Pizorn Frank Verstraete. University of Vienna M ichigan Quantum Summer School
T ensor N et works I ztok Pizorn Frank Verstraete University of Vienna 2010 M ichigan Quantum Summer School Matrix product states (MPS) Introduction to matrix product states Ground states of finite systems
More informationNon-magnetic states. The Néel states are product states; φ N a. , E ij = 3J ij /4 2 The Néel states have higher energy (expectations; not eigenstates)
Non-magnetic states Two spins, i and j, in isolation, H ij = J ijsi S j = J ij [Si z Sj z + 1 2 (S+ i S j + S i S+ j )] For Jij>0 the ground state is the singlet; φ s ij = i j i j, E ij = 3J ij /4 2 The
More informationQuantum s=1/2 antiferromagnet on the Bethe lattice at percolation I. Low-energy states, DMRG, and diagnostics
Quantum s=1/2 antiferromagnet on the Bethe lattice at percolation I. Low-energy states, DMRG, and diagnostics Hitesh J. Changlani, Shivam Ghosh, Sumiran Pujari, Christopher L. Henley Laboratory of Atomic
More informationPositive Tensor Network approach for simulating open quantum many-body systems
Positive Tensor Network approach for simulating open quantum many-body systems 19 / 9 / 2016 A. Werner, D. Jaschke, P. Silvi, M. Kliesch, T. Calarco, J. Eisert and S. Montangero PRL 116, 237201 (2016)
More information4 Matrix product states
Physics 3b Lecture 5 Caltech, 05//7 4 Matrix product states Matrix product state (MPS) is a highly useful tool in the study of interacting quantum systems in one dimension, both analytically and numerically.
More informationarxiv:quant-ph/ v2 24 Dec 2003
Quantum Entanglement in Heisenberg Antiferromagnets V. Subrahmanyam Department of Physics, Indian Institute of Technology, Kanpur, India. arxiv:quant-ph/0309004 v2 24 Dec 2003 Entanglement sharing among
More informationCLASSIFICATION OF SU(2)-SYMMETRIC TENSOR NETWORKS
LSSIFITION OF SU(2)-SYMMETRI TENSOR NETWORKS Matthieu Mambrini Laboratoire de Physique Théorique NRS & Université de Toulouse M.M., R. Orús, D. Poilblanc, Phys. Rev. 94, 2524 (26) D. Poilblanc, M.M., arxiv:72.595
More informationExcursion: MPS & DMRG
Excursion: MPS & DMRG Johannes.Schachenmayer@gmail.com Acronyms for: - Matrix product states - Density matrix renormalization group Numerical methods for simulations of time dynamics of large 1D quantum
More informationFrustration-free Ground States of Quantum Spin Systems 1
1 Davis, January 19, 2011 Frustration-free Ground States of Quantum Spin Systems 1 Bruno Nachtergaele (UC Davis) based on joint work with Sven Bachmann, Spyridon Michalakis, Robert Sims, and Reinhard Werner
More informationUni10 The Universal Tensor Network Library
Uni0 The Universal Tensor Network Library Ying-Jer Kao Department of Physics National Taiwan University National Center for Theoretical Sciences http://www.uni0.org TNQMP 06, ISSP Graphical Representation
More information2D tensor network study of the S=1 bilinear-biquadratic Heisenberg model
2D tensor network study of the S=1 bilinear-biquadratic Heisenberg model Philippe Corboz, Institute for Theoretical Physics, University of Amsterdam AF phase Haldane phase 3-SL 120 phase? ipeps 2D tensor
More informationFrustration without competition: the SU(N) model of quantum permutations on a lattice
Frustration without competition: the SU(N) model of quantum permutations on a lattice F. Mila Ecole Polytechnique Fédérale de Lausanne Switzerland Collaborators P. Corboz (Zürich), A. Läuchli (Innsbruck),
More informationEntanglement in Valence-Bond-Solid States on Symmetric Graphs
Entanglement in Valence-Bond-Solid States on Symmetric Graphs Shu Tanaka A, Hosho Katsura B, Naoki Kawashima C Anatol N. Kirillov D, and Vladimir E. Korepin E A. Kinki University B. Gakushuin University
More informationNumerical diagonalization studies of quantum spin chains
PY 502, Computational Physics, Fall 2016 Anders W. Sandvik, Boston University Numerical diagonalization studies of quantum spin chains Introduction to computational studies of spin chains Using basis states
More informationTensor network renormalization
Walter Burke Institute for Theoretical Physics INAUGURAL CELEBRATION AND SYMPOSIUM Caltech, Feb 23-24, 2015 Tensor network renormalization Guifre Vidal Sherman Fairchild Prize Postdoctoral Fellow (2003-2005)
More informationTensor network simulation of QED on infinite lattices: learning from (1 + 1)d, and prospects for (2 + 1)d
Tensor network simulation of QED on infinite lattices: learning from (1 + 1)d, and prospects for (2 + 1)d Román Orús University of Mainz (Germany) K. Zapp, RO, Phys. Rev. D 95, 114508 (2017) Goal of this
More informationCRYPTOGRAPHIC COMPUTING
CRYPTOGRAPHIC COMPUTING ON GPU Chen Mou Cheng Dept. Electrical Engineering g National Taiwan University January 16, 2009 COLLABORATORS Daniel Bernstein, UIC, USA Tien Ren Chen, Army Tanja Lange, TU Eindhoven,
More informationTensor network simulations of strongly correlated quantum systems
CENTRE FOR QUANTUM TECHNOLOGIES NATIONAL UNIVERSITY OF SINGAPORE AND CLARENDON LABORATORY UNIVERSITY OF OXFORD Tensor network simulations of strongly correlated quantum systems Stephen Clark LXXT[[[GSQPEFS\EGYOEGXMZMXMIWUYERXYQGSYVWI
More informationSolving the sign problem for a class of frustrated antiferromagnets
Solving the sign problem for a class of frustrated antiferromagnets Fabien Alet Laboratoire de Physique Théorique Toulouse with : Kedar Damle (TIFR Mumbai), Sumiran Pujari (Toulouse Kentucky TIFR Mumbai)
More informationSpin liquid phases in strongly correlated lattice models
Spin liquid phases in strongly correlated lattice models Sandro Sorella Wenjun Hu, F. Becca SISSA, IOM DEMOCRITOS, Trieste Seiji Yunoki, Y. Otsuka Riken, Kobe, Japan (K-computer) Williamsburg, 14 June
More informationGround State Projector QMC in the valence-bond basis
Quantum Monte Carlo Methods at Work for Novel Phases of Matter Trieste, Italy, Jan 23 - Feb 3, 2012 Ground State Projector QMC in the valence-bond basis Anders. Sandvik, Boston University Outline: The
More informationQuantum and classical annealing in spin glasses and quantum computing. Anders W Sandvik, Boston University
NATIONAL TAIWAN UNIVERSITY, COLLOQUIUM, MARCH 10, 2015 Quantum and classical annealing in spin glasses and quantum computing Anders W Sandvik, Boston University Cheng-Wei Liu (BU) Anatoli Polkovnikov (BU)
More information2. Spin liquids and valence bond solids
Outline 1. Coupled dimer antiferromagnets Landau-Ginzburg quantum criticality 2. Spin liquids and valence bond solids (a) Schwinger-boson mean-field theory - square lattice (b) Gauge theories of perturbative
More informationStochastic series expansion (SSE) and ground-state projection
Institute of Physics, Chinese Academy of Sciences, Beijing, October 31, 2014 Stochastic series expansion (SSE) and ground-state projection Anders W Sandvik, Boston University Review article on quantum
More informationLecture 3: Tensor Product Ansatz
Lecture 3: Tensor Product nsatz Graduate Lectures Dr Gunnar Möller Cavendish Laboratory, University of Cambridge slide credits: Philippe Corboz (ETH / msterdam) January 2014 Cavendish Laboratory Part I:
More informationEntanglement spectrum and Matrix Product States
Entanglement spectrum and Matrix Product States Frank Verstraete J. Haegeman, D. Draxler, B. Pirvu, V. Stojevic, V. Zauner, I. Pizorn I. Cirac (MPQ), T. Osborne (Hannover), N. Schuch (Aachen) Outline Valence
More informationEntanglement signatures of QED3 in the kagome spin liquid. William Witczak-Krempa
Entanglement signatures of QED3 in the kagome spin liquid William Witczak-Krempa Aspen, March 2018 Chronologically: X. Chen, KITP Santa Barbara T. Faulkner, UIUC E. Fradkin, UIUC S. Whitsitt, Harvard S.
More informationH ψ = E ψ. Introduction to Exact Diagonalization. Andreas Läuchli, New states of quantum matter MPI für Physik komplexer Systeme - Dresden
H ψ = E ψ Introduction to Exact Diagonalization Andreas Läuchli, New states of quantum matter MPI für Physik komplexer Systeme - Dresden http://www.pks.mpg.de/~aml laeuchli@comp-phys.org IPAM QS2009 Workshop,
More informationH ψ = E ψ. Introduction to Exact Diagonalization. Andreas Läuchli, New states of quantum matter MPI für Physik komplexer Systeme - Dresden
H ψ = E ψ Introduction to Exact Diagonalization Andreas Läuchli, New states of quantum matter MPI für Physik komplexer Systeme - Dresden http://www.pks.mpg.de/~aml laeuchli@comp-phys.org Simulations of
More informationc 2017 Society for Industrial and Applied Mathematics
MULTISCALE MODEL. SIMUL. Vol. 15, No. 4, pp. 1423 1447 c 2017 Society for Industrial and Applied Mathematics TENSOR NETWORK SKELETONIZATION LEXING YING Abstract. We introduce a new coarse-graining algorithm,
More informationCluster Density Matrix Embedding Theory for Quantum Spin Systems
Cluster Density Matrix Embedding Theory for Quantum Spin Systems Zhuo Fan 1, and Quan-lin Jie *, 1 1 Department of Physics, Wuhan University, Wuhan 43007, People s Republic of China School of Biomedical
More informationClassical and quantum simulation of dissipative quantum many-body systems
0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 Classical and quantum simulation of dissipative quantum many-body systems
More informationQuantum Monte Carlo Simulations in the Valence Bond Basis. Anders Sandvik, Boston University
Quantum Monte Carlo Simulations in the Valence Bond Basis Anders Sandvik, Boston University Outline The valence bond basis for S=1/2 spins Projector QMC in the valence bond basis Heisenberg model with
More informationJ. Phys.: Condens. Matter 10 (1998) L159 L165. Printed in the UK PII: S (98)90604-X
J. Phys.: Condens. Matter 10 (1998) L159 L165. Printed in the UK PII: S0953-8984(98)90604-X LETTER TO THE EDITOR Calculation of the susceptibility of the S = 1 antiferromagnetic Heisenberg chain with single-ion
More informationGPU-based computation of the Monte Carlo simulation of classical spin systems
Perspectives of GPU Computing in Physics and Astrophysics, Sapienza University of Rome, Rome, Italy, September 15-17, 2014 GPU-based computation of the Monte Carlo simulation of classical spin systems
More informationQuantum Annealing in spin glasses and quantum computing Anders W Sandvik, Boston University
PY502, Computational Physics, December 12, 2017 Quantum Annealing in spin glasses and quantum computing Anders W Sandvik, Boston University Advancing Research in Basic Science and Mathematics Example:
More informationQuantum computing with superconducting qubits Towards useful applications
Quantum computing with superconducting qubits Towards useful applications Stefan Filipp IBM Research Zurich Switzerland Forum Teratec 2018 June 20, 2018 Palaiseau, France Why Quantum Computing? Why now?
More informationSpinon magnetic resonance. Oleg Starykh, University of Utah
Spinon magnetic resonance Oleg Starykh, University of Utah May 17-19, 2018 Examples of current literature 200 cm -1 = 6 THz Spinons? 4 mev = 1 THz The big question(s) What is quantum spin liquid? No broken
More informationTensor Network Computations in Quantum Chemistry. Charles F. Van Loan Department of Computer Science Cornell University
Tensor Network Computations in Quantum Chemistry Charles F. Van Loan Department of Computer Science Cornell University Joint work with Garnet Chan, Department of Chemistry and Chemical Biology, Cornell
More informationQuantum Phase Transitions
1 Davis, September 19, 2011 Quantum Phase Transitions A VIGRE 1 Research Focus Group, Fall 2011 Spring 2012 Bruno Nachtergaele See the RFG webpage for more information: http://wwwmathucdavisedu/~bxn/rfg_quantum_
More informationarxiv: v1 [cond-mat.str-el] 19 Jan 2012
Perfect Sampling with Unitary Tensor Networks arxiv:1201.3974v1 [cond-mat.str-el] 19 Jan 2012 Andrew J. Ferris 1, 2 and Guifre Vidal 1, 3 1 The University of Queensland, School of Mathematics and Physics,
More informationMatrix Product States for Lattice Field Theories
a, K. Cichy bc, J. I. Cirac a, K. Jansen b and H. Saito bd a Max-Planck-Institut für Quantenoptik, Hans-Kopfermann-Str. 1, 85748 Garching, Germany b NIC, DESY Zeuthen, Platanenallee 6, 15738 Zeuthen, Germany
More informationSimulations of Quantum Dimer Models
Simulations of Quantum Dimer Models Didier Poilblanc Laboratoire de Physique Théorique CNRS & Université de Toulouse 1 A wide range of applications Disordered frustrated quantum magnets Correlated fermions
More informationFrustration and Area law
Frustration and Area law When the frustration goes odd S. M. Giampaolo Institut Ruder Bošković, Zagreb, Croatia Workshop: Exactly Solvable Quantum Chains Natal 18-29 June 2018 Coauthors F. Franchini Institut
More informationGPU acceleration of Newton s method for large systems of polynomial equations in double double and quad double arithmetic
GPU acceleration of Newton s method for large systems of polynomial equations in double double and quad double arithmetic Jan Verschelde joint work with Xiangcheng Yu University of Illinois at Chicago
More informationEfficient description of many body systems with prjected entangled pair states
Efficient description of many body systems with prjected entangled pair states J. IGNACIO CIRAC Workshop on Quantum Information Science and many-body systems, National Cheng Kung University, Tainan, Taiwan,
More informationTime-dependent DMRG:
The time-dependent DMRG and its applications Adrian Feiguin Time-dependent DMRG: ^ ^ ih Ψ( t) = 0 t t [ H ( t) E ] Ψ( )... In a truncated basis: t=3 τ t=4 τ t=5τ t=2 τ t= τ t=0 Hilbert space S.R.White
More informationUni10: an open-source library for tensor network algorithms
Journal of Physics: Conference Series PAPER OPEN ACCESS Uni10: an open-source library for tensor network algorithms To cite this article: Ying-Jer Kao et al 2015 J. Phys.: Conf. Ser. 640 012040 View the
More informationShortest Lattice Vector Enumeration on Graphics Cards
Shortest Lattice Vector Enumeration on Graphics Cards Jens Hermans 1 Michael Schneider 2 Fréderik Vercauteren 1 Johannes Buchmann 2 Bart Preneel 1 1 K.U.Leuven 2 TU Darmstadt SHARCS - 10 September 2009
More informationSpin liquids in frustrated magnets
May 20, 2010 Contents 1 Frustration 2 3 4 Exotic excitations 5 Frustration The presence of competing forces that cannot be simultaneously satisfied. Heisenberg-Hamiltonian H = 1 J ij S i S j 2 ij The ground
More informationQuantum annealing for problems with ground-state degeneracy
Proceedings of the International Workshop on Statistical-Mechanical Informatics September 14 17, 2008, Sendai, Japan Quantum annealing for problems with ground-state degeneracy Yoshiki Matsuda 1, Hidetoshi
More informationMatrix product states for the fractional quantum Hall effect
Matrix product states for the fractional quantum Hall effect Roger Mong (California Institute of Technology) University of Virginia Feb 24, 2014 Collaborators Michael Zaletel UC Berkeley (Stanford/Station
More informationIntroduction to Quantum Monte Carlo
Entanglement in Strongly Correlated Systems @ Benasque Feb. 6-17, 2017 Introduction to Quantum Monte Carlo Naoki KAWASHIMA (ISSP) 2017.02.06-07 Why bother? Estimating scaling dimension by TRG, TNR, etc
More informationFermionic tensor networks
Fermionic tensor networks Philippe Corboz, ETH Zurich / EPF Lausanne, Switzerland Collaborators: University of Queensland: Guifre Vidal, Roman Orus, Glen Evenbly, Jacob Jordan University of Vienna: Frank
More informationSpin liquids on ladders and in 2d
Spin liquids on ladders and in 2d MPA Fisher (with O. Motrunich) Minnesota, FTPI, 5/3/08 Interest: Quantum Spin liquid phases of 2d Mott insulators Background: Three classes of 2d Spin liquids a) Topological
More informationFrustration-free Ground States of Quantum Spin Systems 1
1 FRG2011, Harvard, May 19, 2011 Frustration-free Ground States of Quantum Spin Systems 1 Bruno Nachtergaele (UC Davis) based on joint work with Sven Bachmann, Spyridon Michalakis, Robert Sims, and Reinhard
More informationAccelerating Quantum Chromodynamics Calculations with GPUs
Accelerating Quantum Chromodynamics Calculations with GPUs Guochun Shi, Steven Gottlieb, Aaron Torok, Volodymyr Kindratenko NCSA & Indiana University National Center for Supercomputing Applications University
More informationFrom real materials to model Hamiltonians with density matrix downfolding
From real materials to model Hamiltonians with density matrix downfolding (Converting a many-electron problem to a fewer-electron one) Hitesh J. Changlani Florida State University Based on HJC, H. Zheng,
More informationA Quantum Chemistry Domain-Specific Language for Heterogeneous Clusters
A Quantum Chemistry Domain-Specific Language for Heterogeneous Clusters ANTONINO TUMEO, ORESTE VILLA Collaborators: Karol Kowalski, Sriram Krishnamoorthy, Wenjing Ma, Simone Secchi May 15, 2012 1 Outline!
More informationAntti-Pekka Hynninen, 5/10/2017, GTC2017, San Jose CA
S7255: CUTT: A HIGH- PERFORMANCE TENSOR TRANSPOSE LIBRARY FOR GPUS Antti-Pekka Hynninen, 5/10/2017, GTC2017, San Jose CA MOTIVATION Tensor contractions are the most computationally intensive part of quantum
More information