Arxiv:2008.11206V2 [Cond-Mat.Str-El] 30 Apr 2021 Topological Paramagnets, Protected by Z2 and Z2 Symmetries, 47, 49, 52–55]

Home , Numerical sign problem

AAAB83icbVDLSsNAFL2pr1pfVZduBovgqiSloO6KblxWsK3YhDKZTtqhk0mYh1BCf8ONC0Xc+jPu/BsnbRbaemDgcM693DMnTDlT2nW/ndLa+sbmVnm7srO7t39QPTzqqsRIQjsk4Yl8CLGinAna0Uxz+pBKiuOQ0144ucn93hOViiXiXk9TGsR4JFjECNZW8v0Y63EYZo+zQWNQrbl1dw60SryC1KBAe1D98ocJMTEVmnCsVN9zUx1kWGpGOJ1VfKNoiskEj2jfUoFjqoJsnnmGzqwyRFEi7RMazdXfGxmOlZrGoZ3MM6plLxf/8/pGR5dBxkRqNBVkcSgyHOkE5QWgIZOUaD61BBPJbFZExlhiom1NFVuCt/zlVdJt1L1m/equWWtdF3WU4QRO4Rw8uIAW3EIbOkAghWd4hTfHOC/Ou/OxGC05xc4x/IHz+QP7v5Gr

AAAB9XicbVDLSsNAFL2pr1pfVZduBovgqiQ1oMuiG5cV7APbtEymk3boZBJmJkoJ/Q83LhRx67+482+ctFlo64GBwzn3cs8cP+ZMadv+tgpr6xubW8Xt0s7u3v5B+fCopaJEEtokEY9kx8eKciZoUzPNaSeWFIc+p21/cpP57UcqFYvEvZ7G1AvxSLCAEayN1O+FWI99P32YDWr9i0G5YlftOdAqcXJSgRyNQfmrN4xIElKhCcdKdR071l6KpWaE01mplygaYzLBI9o1VOCQKi+dp56hM6MMURBJ84RGc/X3RopDpaahbyazlGrZy8T/vG6igysvZSJONBVkcShIONIRyipAQyYp0XxqCCaSmayIjLHERJuiSqYEZ/nLq6RVqzpu1b1zK/XrvI4inMApnIMDl1CHW2hAEwhIeIZXeLOerBfr3fpYjBasfOcY/sD6/AEptZJL

AAAB7XicbVDLSgNBEOz1GeMr6tHLYBA8hV0J6DHoxWME84BkCb2TSTJmdmaZmRXCkn/w4kERr/6PN//GSbIHTSxoKKq66e6KEsGN9f1vb219Y3Nru7BT3N3bPzgsHR03jUo1ZQ2qhNLtCA0TXLKG5VawdqIZxpFgrWh8O/NbT0wbruSDnSQsjHEo+YBTtE5qdlEkI+yVyn7Fn4OskiAnZchR75W+un1F05hJSwUa0wn8xIYZasupYNNiNzUsQTrGIes4KjFmJszm107JuVP6ZKC0K2nJXP09kWFszCSOXGeMdmSWvZn4n9dJ7eA6zLhMUsskXSwapIJYRWavkz7XjFoxcQSp5u5WQkeokVoXUNGFECy/vEqal5WgWqneV8u1mzyOApzCGVxAAFdQgzuoQwMoPMIzvMKbp7wX7937WLSuefnMCfyB9/kDjcePHw==

AAAB6XicbVBNS8NAEJ3Ur1q/qh69LBbBiyWRgh6LXjxWsR/QhrLZbtqlm03YnQgl9B948aCIV/+RN/+N2zYHbX0w8Hhvhpl5QSKFQdf9dgpr6xubW8Xt0s7u3v5B+fCoZeJUM95ksYx1J6CGS6F4EwVK3kk0p1EgeTsY38789hPXRsTqEScJ9yM6VCIUjKKVHi68frniVt05yCrxclKBHI1++as3iFkacYVMUmO6npugn1GNgkk+LfVSwxPKxnTIu5YqGnHjZ/NLp+TMKgMSxtqWQjJXf09kNDJmEgW2M6I4MsveTPzP66YYXvuZUEmKXLHFojCVBGMye5sMhOYM5cQSyrSwtxI2opoytOGUbAje8surpHVZ9WrV2n2tUr/J4yjCCZzCOXhwBXW4gwY0gUEIz/AKb87YeXHenY9Fa8HJZ47hD5zPH+YojPM=

AAAB63icbVBNS8NAEJ34WetX1aOXxSJ4MSRSUW9FLx4r2A9oQ9lsN+3S3U3Y3Qgl9C948aCIV/+QN/+NmzYHbX0w8Hhvhpl5YcKZNp737aysrq1vbJa2yts7u3v7lYPDlo5TRWiTxDxWnRBrypmkTcMMp51EUSxCTtvh+C73209UaRbLRzNJaCDwULKIEWxy6dxzL/uVqud6M6Bl4hekCgUa/cpXbxCTVFBpCMdad30vMUGGlWGE02m5l2qaYDLGQ9q1VGJBdZDNbp2iU6sMUBQrW9Kgmfp7IsNC64kIbafAZqQXvVz8z+umJroOMiaT1FBJ5ouilCMTo/xxNGCKEsMnlmCimL0VkRFWmBgbT9mG4C++vExaF65fc28eatX6bRFHCY7hBM7Ahyuowz00oAkERvAMr/DmCOfFeXc+5q0rTjFzBH/gfP4Axw+Nbg==

AAAB6HicbVBNS8NAEJ3Ur1q/qh69LBbBU0mkoN6KXjy2YGuhDWWznbRrN5uwuxFK6C/w4kERr/4kb/4bt20O2vpg4PHeDDPzgkRwbVz32ymsrW9sbhW3Szu7e/sH5cOjto5TxbDFYhGrTkA1Ci6xZbgR2EkU0igQ+BCMb2f+wxMqzWN5byYJ+hEdSh5yRo2Vmm6/XHGr7hxklXg5qUCORr/81RvELI1QGiao1l3PTYyfUWU4Ezgt9VKNCWVjOsSupZJGqP1sfuiUnFllQMJY2ZKGzNXfExmNtJ5Ege2MqBnpZW8m/ud1UxNe+RmXSWpQssWiMBXExGT2NRlwhcyIiSWUKW5vJWxEFWXGZlOyIXjLL6+S9kXVq1Wvm7VK/SaPowgncArn4MEl1OEOGtACBgjP8ApvzqPz4rw7H4vWgpPPHMMfOJ8/fPmMwA==

AAAB63icbVBNS8NAEJ34WetX1aOXxSIIQkikot6KXjxWsB/QhrLZbtqlu5uwuxFK6F/w4kERr/4hb/4bN20O2vpg4PHeDDPzwoQzbTzv21lZXVvf2Cxtlbd3dvf2KweHLR2nitAmiXmsOiHWlDNJm4YZTjuJoliEnLbD8V3ut5+o0iyWj2aS0EDgoWQRI9jk0rnnXvYrVc/1ZkDLxC9IFQo0+pWv3iAmqaDSEI617vpeYoIMK8MIp9NyL9U0wWSMh7RrqcSC6iCb3TpFp1YZoChWtqRBM/X3RIaF1hMR2k6BzUgvern4n9dNTXQdZEwmqaGSzBdFKUcmRvnjaMAUJYZPLMFEMXsrIiOsMDE2nrINwV98eZm0Lly/5t481Kr12yKOEhzDCZyBD1dQh3toQBMIjOAZXuHNEc6L8+58zFtXnGLmCP7A+fwBxAGNbA==

AAAB6XicbVBNS8NAEJ3Ur1q/qh69LBZBEEoiBT0WvXisYj+gDWWz3bRLN5uwOxFK6D/w4kERr/4jb/4bt20O2vpg4PHeDDPzgkQKg6777RTW1jc2t4rbpZ3dvf2D8uFRy8SpZrzJYhnrTkANl0LxJgqUvJNoTqNA8nYwvp357SeujYjVI04S7kd0qEQoGEUrPVx4/XLFrbpzkFXi5aQCORr98ldvELM04gqZpMZ0PTdBP6MaBZN8WuqlhieUjemQdy1VNOLGz+aXTsmZVQYkjLUthWSu/p7IaGTMJApsZ0RxZJa9mfif100xvPYzoZIUuWKLRWEqCcZk9jYZCM0ZyokllGlhbyVsRDVlaMMp2RC85ZdXSeuy6tWqtftapX6Tx1GEEziFc/DgCupwBw1oAoMQnuEV3pyx8+K8Ox+L1oKTzxzDHzifP+MejPE=

AAAB9HicdVDLSgMxFM3UV62vqks3wSK4GmbaYqe7oiBuhAr2Ae1YMmnahiaZMckUytDvcONCEbd+jDv/xkxbQUUPBA7n3Ms9OUHEqNKO82FlVlbX1jeym7mt7Z3dvfz+QVOFscSkgUMWynaAFGFUkIammpF2JAniASOtYHyR+q0JkYqG4lZPI+JzNBR0QDHSRvK7HOmR5Mnl9eyu1MsXHNvzPLdUgY5dqlbO3GJK3Krh0LWdOQpgiXov/97thzjmRGjMkFId14m0nyCpKWZkluvGikQIj9GQdAwViBPlJ/PQM3hilD4chNI8oeFc/b6RIK7UlAdmMg2pfnup+JfXifXA8xMqolgTgReHBjGDOoRpA7BPJcGaTQ1BWFKTFeIRkghr01POlPD1U/g/aRZtt2xXb8qF2vmyjiw4AsfgFLigAmrgCtRBA2BwDx7AE3i2Jtaj9WK9LkYz1nLnEPyA9fYJLrCSaA==

AAAB7nicbVBNS8NAEJ3Ur1q/qh69LBbBU0lKQY9FLx4r9ENoQ9lsN+3SzSbsTgol9Ed48aCIV3+PN/+N2zYHbX0w8Hhvhpl5QSKFQdf9dgpb2zu7e8X90sHh0fFJ+fSsY+JUM95msYz1U0ANl0LxNgqU/CnRnEaB5N1gcr/wu1OujYhVC2cJ9yM6UiIUjKKVui0tpoLKQbniVt0lyCbxclKBHM1B+as/jFkacYVMUmN6npugn1GNgkk+L/VTwxPKJnTEe5YqGnHjZ8tz5+TKKkMSxtqWQrJUf09kNDJmFgW2M6I4NuveQvzP66UY3vqZUEmKXLHVojCVBGOy+J0MheYM5cwSyrSwtxI2ppoytAmVbAje+subpFOrevVq/bFWadzlcRThAi7hGjy4gQY8QBPawGACz/AKb07ivDjvzseqteDkM+fwB87nD2wXj6A=

AAAB83icdVDLSgMxFM34rPVVdekmWARXw0wrdrorduOygn1AO5RMmmlDk8yQh1CG/oYbF4q49Wfc+Tdm2goqeuDC4Zx7k3tPlDKqtOd9OGvrG5tb24Wd4u7e/sFh6ei4oxIjMWnjhCWyFyFFGBWkralmpJdKgnjESDeaNnO/e0+koom407OUhByNBY0pRtpKg2bCORHKSKTJsFT23CAI/GoNem61XrvyKznx65ZD3/UWKIMVWsPS+2CUYGMf0Jghpfq+l+owQ1JTzMi8ODCKpAhP0Zj0LRWIExVmi53n8NwqIxgn0pbQcKF+n8gQV2rGI9vJkZ6o314u/uX1jY6DMKMiNZoIvPwoNgzqBOYBwBGVBGs2swRhSe2uEE+QRFjbmIo2hK9L4f+kU3H9S7d+Wyk3rldxFMApOAMXwAc10AA3oAXaAIMUPIAn8OwY59F5cV6XrWvOauYE/IDz9gnuspJO

AAAB6nicdVDJSgNBEK1xjXGLevTSGARPQ88wmhyDXjxGskIyhJ5OT9KkZ6G7RwhDPsGLB0W8+kXe/Bs7i6CiDwoe71VRVS9IBVca4w9rbX1jc2u7sFPc3ds/OCwdHbdVkknKWjQRiewGRDHBY9bSXAvWTSUjUSBYJ5jczP3OPZOKJ3FTT1PmR2QU85BToo3UaNSbg1IZ227FwV4VYRtf4ivXNcTDVcf1kGPjBcqwQn1Qeu8PE5pFLNZUEKV6Dk61nxOpORVsVuxniqWETsiI9QyNScSUny9OnaFzowxRmEhTsUYL9ftETiKlplFgOiOix+q3Nxf/8nqZDqt+zuM00yymy0VhJpBO0PxvNOSSUS2mhhAqubkV0TGRhGqTTtGE8PUp+p+0XdvxbO/OLdeuV3EU4BTO4AIcqEANbqEOLaAwggd4gmdLWI/Wi/W6bF2zVjMn8APW2ydEU43M

AAAB9XicdVDLSgMxFM3UV62vqks3wSK4GmZasdNd0Y3uKlhbaMeSSTNtaB5DklFK6X+4caGIW//FnX9jpq2gogcuHM65N7n3RAmj2njeh5NbWl5ZXcuvFzY2t7Z3irt7N1qmCpMmlkyqdoQ0YVSQpqGGkXaiCOIRI61odJ75rTuiNJXi2owTEnI0EDSmGBkr3V4KLDknQqcKGdIrljw3CAK/UoWeW6lVT/1yRvya5dB3vRlKYIFGr/je7Uuc2gcMZkjrju8lJpwgZShmZFroppokCI/QgHQsFYgTHU5mW0/hkVX6MJbKljBwpn6fmCCu9ZhHtpMjM9S/vUz8y+ukJg7CCRVJaojA84/ilEEjYRYB7FNFsGFjSxBW1O4K8RAphI0NqmBD+LoU/k9uyq5/4tauyqX62SKOPDgAh+AY+KAK6uACNEATYKDAA3gCz8698+i8OK/z1pyzmNkHP+C8fQKI35M5

AAAB9HicdVDLSgMxFM3UV62vqks3wSK6GjLDaOuu6MZlxb6gHUomTdvQZGZMMoUy9DvcuFDErR/jzr8x01ZQ0QOBwzn3ck9OEHOmNEIfVm5ldW19I79Z2Nre2d0r7h80VZRIQhsk4pFsB1hRzkLa0Exz2o4lxSLgtBWMrzO/NaFSsSis62lMfYGHIRswgrWR/K7AeiRFelerz057xRKy3bKDvApENjpHF65riIcqjutBx0ZzlMAStV7xvduPSCJoqAnHSnUcFGs/xVIzwums0E0UjTEZ4yHtGBpiQZWfzkPP4IlR+nAQSfNCDefq940UC6WmIjCTWUj128vEv7xOogcVP2VhnGgaksWhQcKhjmDWAOwzSYnmU0MwkcxkhWSEJSba9FQwJXz9FP5Pmq7tePblrVeqXi3ryIMjcAzOgAPKoApuQA00AAH34AE8gWdrYj1aL9brYjRnLXcOwQ9Yb5/tC5I7

AAAB7XicdVBNS8NAEJ3Ur1q/qh69LBbBU0iq2PRW9OKxgv2ANpTNdtOu3WzC7kYoof/BiwdFvPp/vPlv3LQVVPTBwOO9GWbmBQlnSjvOh1VYWV1b3yhulra2d3b3yvsHbRWnktAWiXksuwFWlDNBW5ppTruJpDgKOO0Ek6vc79xTqVgsbvU0oX6ER4KFjGBtpLbSkiV0UK44tud57lkNOfZZvXbhVnPi1g1Hru3MUYElmoPye38YkzSiQhOOleq5TqL9DEvNCKezUj9VNMFkgke0Z6jAEVV+Nr92hk6MMkRhLE0Jjebq94kMR0pNo8B0RliP1W8vF//yeqkOPT9jIkk1FWSxKEw50jHKX0dDJinRfGoIJpKZWxEZY4mJNgGVTAhfn6L/Sbtqu+d2/aZaaVwu4yjCERzDKbhQgwZcQxNaQOAOHuAJnq3YerRerNdFa8FazhzCD1hvn2mRj7Y=

AAACDXicdVC7SgNBFJ31GeMramkzGIXYLLtJMJsuaGMZwTwgG8PsZJIMmX04c1cIS37Axl+xsVDE1t7Ov3HyUFT0wMC559zLnXu8SHAFlvVuLCwuLa+sptbS6xubW9uZnd26CmNJWY2GIpRNjygmeMBqwEGwZiQZ8T3BGt7wbOI3bphUPAwuYRSxtk/6Ae9xSkBLncyh6/G+yLnqWkJSGLvAfaa+Ku3J46tCJ5O1TMdx7EIJW2ahXDqx8xNilzXHtmlNkUVzVDuZN7cb0thnAVBBlGrZVgTthEjgVLBx2o0Viwgdkj5raRoQvbSdTK8Z4yOtdHEvlPoFgKfq94mE+EqNfE93+gQG6rc3Ef/yWjH0nHbCgygGFtDZol4sMIR4Eg3ucskoiJEmhEqu/4rpgEhCQQeY1iF8Xor/J/W8aRfN8kUxWzmdx5FC++gA5ZCNSqiCzlEV1RBFt+gePaIn4854MJ6Nl1nrgjGf2UM/YLx+AHSpnHo=

3 From trivial to topological paramagnets: The case of Z2 and Z2 symmetries in two dimensions

Maxime Dupont,1, 2 Snir Gazit,3 and Thomas Scaffidi4 1Department of Physics, University of California, Berkeley, California 94720, USA 2Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA 3Racah Institute of Physics and the Fritz Haber Center for Molecular Dynamics, The Hebrew University, Jerusalem 91904, Israel 4Department of Physics, University of Toronto, Toronto, Ontario, M5S 1A7, Canada Using quantum Monte Carlo simulations, we map out the phase diagram of Hamiltonians interpolating between Z Z3 trivial and nontrivial bosonic symmetry-protected topological phases, protected by 2 and 2 symmetries, in two dimensions. In all cases, we find that the trivial and the topological phases are separated by an intermediate phase in which the protecting symmetry is spontaneously broken. Depending√ on√ the model, we identify a variety of magnetic orders on the triangular lattice, including ferromagnetism, 3 × 3 order, and stripe orders (both commensurate and incommensurate). Critical properties are determined through a finite-size scaling analysis. Possible scenarios regarding the nature of the phase transitions are discussed.

I. INTRODUCTION p p 3 3 (a) SPT0 3 3 Trivial FM SPT AAACz3icjVHLTsJAFD3UF+ILdemmkZjghhQ00SXRjUtIBEyAmGkZoKGvtFMNIRi3/oBb/SvjH+hfeGccEpUYnabtmXPvOTP3Xjvy3ERY1mvGWFhcWl7JrubW1jc2t/LbO80kTGOHN5zQC+MrmyXccwPeEK7w+FUUc+bbHm/Zo3MZb93wOHHD4FKMI9712SBw+67DBFGdYsdnYmj3J2x6eJ0vWCVLLXMelDUoQK9amH9BBz2EcJDCB0cAQdgDQ0JPG2VYiIjrYkJcTMhVcY4pcqRNKYtTBiN2RN8B7dqaDWgvPROldugUj96YlCYOSBNSXkxYnmaqeKqcJfub90R5yruN6W9rL59YgSGxf+lmmf/VyVoE+jhVNbhUU6QYWZ2jXVLVFXlz80tVghwi4iTuUTwm7CjlrM+m0iSqdtlbpuJvKlOycu/o3BTv8pY04PLPcc6DZqVUPipV6seF6pkedRZ72EeR5nmCKi5QQ4O8IzziCc9G3bg17oz7z1QjozW7+LaMhw+kuJPv ⇥ 3 ↵ The Landau-Ginzburg-Wilson (LGW) paradigm for classi- Z2 1 0.5 0 +0.5 +1 fying phases of matter, based on their symmetry and its sub- sequent breaking [1,2], has been challenged by the discovery Commensurate Incommensurate of symmetry-protected topological (SPT) phases [3–9]. Of stripe Trivial stripe (b) SPT0 SPT particular interest are SPT phases that arise due to strong cor- AAACz3icjVHLTsJAFD3UF+ILdemmkZjghhQ00SXRjUtIBEyAmLYM0NBXplMNIRi3/oBb/SvjH+hfeGccEpUYnabtmXPvOTP3Xif2vURY1mvGWFhcWl7JrubW1jc2t/LbO80kSrnLGm7kR/zKsRPmeyFrCE/47CrmzA4cn7Wc0bmMt24YT7wovBTjmHUDexB6fc+1BVGdYiewxdDpT5zp4XW+YJUstcx5UNagAL1qUf4FHfQQwUWKAAwhBGEfNhJ62ijDQkxcFxPiOCFPxRmmyJE2pSxGGTaxI/oOaNfWbEh76ZkotUun+PRyUpo4IE1EeZywPM1U8VQ5S/Y374nylHcb09/RXgGxAkNi/9LNMv+rk7UI9HGqavCoplgxsjpXu6SqK/Lm5peqBDnExEncozgn7CrlrM+m0iSqdtlbW8XfVKZk5d7VuSne5S1pwOWf45wHzUqpfFSq1I8L1TM96iz2sI8izfMEVVyghgZ5x3jEE56NunFr3Bn3n6lGRmt28W0ZDx+nGpPw ↵ relations, generalizing the free fermion band-structure classifi- Z2 1 0.5 0 +0.5 +1 cation of topological insulators to generic interacting systems. SPT phases sharing the same protecting symmetry display identical physical properties in the bulk and are indistinguish- FIG. 1. Quantum phase diagram for the Hamiltonians (2.3). For Z3 Z able by symmetry-based probes such as local order param- each symmetry ( 2 and 2), the Hamiltonian interpolates between eters. Instead, the distinction between trivial and nontrivial the parent Hamiltonian of a paramagnet at 훼 = 0, the Hamiltonian of a nontrivial SPT at 훼 = 1, and the Hamiltonian of another nontrivial SPT phases is more subtle and manifests itself after gauging SPT (called SPT’) for 훼 = −1. The only difference between SPT and the protecting symmetry, or through properties like string or- SPT’ is a minus sign in front of the parent Hamiltonian, which confers der parameters, edge states, entanglement spectrum, strange a nontrivial weak 0D index to SPT’. In all cases, the transition from correlators, etc. [3–15]. Well known examples include the the trivial to the SPT phases happens via an intermediate symmetry- Haldane phase of odd-integer Heisenberg chains [8,9, 16–22] breaking phase, where the protecting symmetry is spontaneously Z3 and the bosonic integer quantum Hall phases in two dimen- broken. The 2 model at positive and negative 훼 as well as the sions [23–30], among others. Z2 model for 훼 < 0 are studied in this paper. The other case of the While the topological classification of SPT phases is by Z2 model with 훼 > 0 was thoroughly studied in Ref. 59. See TableI for a more quantitative summary of our findings. now fairly well established [3–9], our understanding of quantum phase transitions involving them is still lacking. Such transitions are expected to give rise to novel quantum critical behavior, going beyond the LGW predictions. In that respect, in many cases, quantum Monte Carlo techniques are plagued phase transitions between SPT phases and the more familiar by the numerical “sign problem” rendering the statistical er- symmetry-broken, ordered states [31–40] as well as transitions rors uncontrolled [62]. Therefore, identifying concrete lattice between trivial and nontrivial SPT phases [27, 30, 41–52] have models that exhibit transitions involving SPTs and that are both attracted tremendous attention. Previous works have amenable to an unbiased numerical solution is desirable. mostly considered transitions between SPT phases that are In this paper, we numerically investigate the phase diagram protected by continuous symmetries and uncovered remark- of two different models that interpolate between trivial and able relations with deconfined quantum criticality [30, 45– 3 arXiv:2008.11206v2 [cond-mat.str-el] 30 Apr 2021 topological paramagnets, protected by Z2 and Z2 symmetries, 47, 49, 52–55]. However, the study of microscopic mod- see also Ref. 59. Crucially, the interpolation does not explic- els with discrete symmetry groups has remained relatively itly break the protecting symmetry, and since it connects two scarce [37, 42, 48, 56–58]. distinct phases of matter, we expect a quantum phase transition In the absence of an overarching theoretical framework of to occur along the way. This can happen either through a single SPT criticality, one must resort to exact numerical methods transition point separating the two phases or a two-step tran- to determine their properties in specific cases. In that regard, sition via an intermediate phase, which spontaneously breaks in dimensions 퐷 > 1, exact diagonalization and density ma- the protecting symmetry. In all cases studied in this work, we trix renormalization group (DMRG) techniques [60, 61] are find that the latter scenario is realized, giving rise to an inter- restricted to small system sizes, which typically do not al- mediate magnetically ordered phase, see Fig.1. Interestingly, low studies of long-wavelength universal properties. Also, in certain cases, magnetic order is accompanied by additional AAAB6HicbVDLTgJBEOzFF+IL9ehlIjHxRHYNRo9ELx4hkUcCGzI79MLI7OxmZtZICF/gxYPGePWTvPk3DrAHBSvppFLVne6uIBFcG9f9dnJr6xubW/ntws7u3v5B8fCoqeNUMWywWMSqHVCNgktsGG4EthOFNAoEtoLR7cxvPaLSPJb3ZpygH9GB5CFn1Fip/tQrltyyOwdZJV5GSpCh1it+dfsxSyOUhgmqdcdzE+NPqDKcCZwWuqnGhLIRHWDHUkkj1P5kfuiUnFmlT8JY2ZKGzNXfExMaaT2OAtsZUTPUy95M/M/rpCa89idcJqlByRaLwlQQE5PZ16TPFTIjxpZQpri9lbAhVZQZm03BhuAtv7xKmhdlr1K+rFdK1ZssjjycwCmcgwdXUIU7qEEDGCA8wyu8OQ/Oi/PufCxac042cwx/4Hz+AOjRjQQ=

AAAB8XicdVDLSgMxFL1TX7W+qi7dBIvgqszU0dpdsRuXFewD26Fk0kwbmskMSUYoQ//CjQtF3Po37vwb03YEFT0QOJxzLzn3+DFnStv2h5VbWV1b38hvFra2d3b3ivsHbRUlktAWiXgkuz5WlDNBW5ppTruxpDj0Oe34k8bc79xTqVgkbvU0pl6IR4IFjGBtpLt+iPVYBWljNiiW7LJTdd1qFdllewFDzi7O7UoNOZlSggzNQfG9P4xIElKhCcdK9Rw71l6KpWaE01mhnygaYzLBI9ozVOCQKi9dJJ6hE6MMURBJ84RGC/X7RopDpaahbyYXCX97c/Evr5fo4NJLmYgTTQVZfhQkHOkIzc9HQyYp0XxqCCaSmayIjLHERJuSCqaEr0vR/6RdKTtuuXbjlupXWR15OIJjOAUHqlCHa2hCCwgIeIAneLaU9Wi9WK/L0ZyV7RzCD1hvnxpYkTs=

AAAB8XicdVDLSgMxFM34rPVVdekmWARXQ2Y6tHVX6sZlBfvAdiiZNNOGZjJDkhHK0L9w40IRt/6NO//GdFpBRQ8EDufcS849QcKZ0gh9WGvrG5tb24Wd4u7e/sFh6ei4o+JUEtomMY9lL8CKciZoWzPNaS+RFEcBp91gerXwu/dUKhaLWz1LqB/hsWAhI1gb6W4QYT1RYdacD0tlZFc9p+JWIbIdt1J3XUM8VKugKnRslKMMVmgNS++DUUzSiApNOFaq76BE+xmWmhFO58VBqmiCyRSPad9QgSOq/CxPPIfnRhnBMJbmCQ1z9ftGhiOlZlFgJvOEv72F+JfXT3VY9zMmklRTQZYfhSmHOoaL8+GISUo0nxmCiWQmKyQTLDHRpqSiKeHrUvg/6bi249mXN1650VzVUQCn4AxcAAfUQANcgxZoAwIEeABP4NlS1qP1Yr0uR9es1c4J+AHr7RMgBJE/

AAAB8XicdVDLSgMxFL1TX7W+qi7dBIvgakhf2O6qblxWsK3YlpJJM21oJjMkGaEM/Qs3LhRx69+4829MpxVU9EDgcM695NzjRYJrg/GHk1lZXVvfyG7mtrZ3dvfy+wdtHcaKshYNRahuPaKZ4JK1DDeC3UaKkcATrONNLud+554pzUN5Y6YR6wdkJLnPKTFWuusFxIy1n5zPBvkCdus1jMslhN1qFddL2JJytVKuYVR0cYoCLNEc5N97w5DGAZOGCqJ1t4gj00+IMpwKNsv1Ys0iQidkxLqWShIw3U/SxDN0YpUh8kNlnzQoVb9vJCTQehp4djJN+Nubi3953dj4tX7CZRQbJuniIz8WyIRofj4acsWoEVNLCFXcZkV0TBShxpaUsyV8XYr+J+2SW6y49etKoXGxrCMLR3AMp1CEM2jAFTShBRQkPMATPDvaeXRenNfFaMZZ7hzCDzhvnydCkUQ=

AAAB6HicbVBNS8NAEJ3Ur1q/qh69LBbBU0mkoseiF48t2FpoQ9lsJ+3azSbsboQS+gu8eFDEqz/Jm//GbZuDtj4YeLw3w8y8IBFcG9f9dgpr6xubW8Xt0s7u3v5B+fCoreNUMWyxWMSqE1CNgktsGW4EdhKFNAoEPgTj25n/8IRK81jem0mCfkSHkoecUWOl5qRfrrhVdw6ySrycVCBHo1/+6g1ilkYoDRNU667nJsbPqDKcCZyWeqnGhLIxHWLXUkkj1H42P3RKzqwyIGGsbElD5urviYxGWk+iwHZG1Iz0sjcT//O6qQmv/YzLJDUo2WJRmApiYjL7mgy4QmbExBLKFLe3EjaiijJjsynZELzll1dJ+6Lq1aqXzVqlfpPHUYQTOIVz8OAK6nAHDWgBA4RneIU359F5cd6dj0VrwclnjuEPnM8f6lWNBQ==

AAAB+nicbVDLSgNBEJz1GeNro0cvg0HwFHaDoN6CXjxGyAuSJcxOepMhsw9mepWw5lO8eFDEq1/izb9xkuxBEwsaiqrume7yEyk0Os63tba+sbm1Xdgp7u7tHxzapaOWjlPFocljGauOzzRIEUETBUroJApY6Eto++Pbmd9+AKVFHDVwkoAXsmEkAsEZGqlvlxojBUB16kuGKDjovl12Ks4cdJW4OSmTHPW+/dUbxDwNIUIumdZd10nQy5gyz0mYFnuphoTxMRtC19CIhaC9bL76lJ4ZZUCDWJmKkM7V3xMZC7WehL7pDBmO9LI3E//zuikGV14moiRFiPjioyCVFGM6y4EOhAKOcmII40qYXSkfMcU4mrSKJgR3+eRV0qpW3IvK9X21XLvJ4yiQE3JKzolLLkmN3JE6aRJOHskzeSVv1pP1Yr1bH4vWNSufOSZ/YH3+AGMblBo=

AAAB9XicbVDLTgJBEJzFF+IL9ehlIjHxRHYJiR6JXjxiIo8EkMwOvTBhdnad6dWQDf/hxYPGePVfvPk3DrAHBSvppFLVne4uP5bCoOt+O7m19Y3Nrfx2YWd3b/+geHjUNFGiOTR4JCPd9pkBKRQ0UKCEdqyBhb6Elj++nvmtR9BGROoOJzH0QjZUIhCcoZXu65I9JIAIFEGH/WLJLbtz0FXiZaREMtT7xa/uIOJJCAq5ZMZ0PDfGXso0Ci5hWugmBmLGx2wIHUsVC8H00vnVU3pmlQENIm1LIZ2rvydSFhozCX3bGTIcmWVvJv7ndRIMLnupUHGCoPhiUZBIihGdRUAHQgNHObGEcS3srZSPmGbcZmAKNgRv+eVV0qyUvWq5elsp1a6yOPLkhJySc+KRC1IjN6ROGoQTTZ7JK3lznpwX5935WLTmnGzmmPyB8/kDoGiSmA==

AAAB83icbVBNS8NAEN34WetX1aOXxSJ4sSSloMeiF29WsB/QhrLZTtqlm03YnQgl9G948aCIV/+MN/+N2zYHbX0w8Hhvhpl5QSKFQdf9dtbWNza3tgs7xd29/YPD0tFxy8Sp5tDksYx1J2AGpFDQRIESOokGFgUS2sH4dua3n0AbEatHnCTgR2yoRCg4Qyv17tWlEQgUQUf9UtmtuHPQVeLlpExyNPqlr94g5mkECrlkxnQ9N0E/YxoFlzAt9lIDCeNjNoSupYpFYPxsfvOUnltlQMNY21JI5+rviYxFxkyiwHZGDEdm2ZuJ/3ndFMNrPxMqSREUXywKU0kxprMA6EBo4CgnljCuhb2V8hHTjNsMTNGG4C2/vEpa1YpXq9QequX6TR5HgZySM3JBPHJF6uSONEiTcJKQZ/JK3pzUeXHenY9F65qTz5yQP3A+fwCkk5Fu

broken symmetries, such as lattice translations or point group (c) Three sublattices On-site term (b) Plaquette term AAACz3icjVHLTsJAFD3UF+ILdemmkZjghhQ00SXRjUtIBEyAmHYYoKGvtFMNIRi3/oBb/SvjH+hfeGccEpUYnabtmXPvOTP3Xify3ERY1mvGWFhcWl7JrubW1jc2t/LbO80kTGPGGyz0wvjKsRPuuQFvCFd4/CqKue07Hm85o3MZb93wOHHD4FKMI9717UHg9l1mC6I6xY5vi6HTn7Dp4XW+YJUstcx5UNagAL1qYf4FHfQQgiGFD44AgrAHGwk9bZRhISKuiwlxMSFXxTmmyJE2pSxOGTaxI/oOaNfWbEB76ZkoNaNTPHpjUpo4IE1IeTFheZqp4qlyluxv3hPlKe82pr+jvXxiBYbE/qWbZf5XJ2sR6ONU1eBSTZFiZHVMu6SqK/Lm5peqBDlExEnco3hMmCnlrM+m0iSqdtlbW8XfVKZk5Z7p3BTv8pY04PLPcc6DZqVUPipV6seF6pkedRZ72EeR5nmCKi5QQ4O8IzziCc9G3bg17oz7z1QjozW7+LaMhw+pfJPx (a) AAACz3icjVHLTsJAFD3UF+ILdemmkZjghhQ00SXRjUtIBEyAmLYM0NBXplMNIRi3/oBb/SvjH+hfeGccEpUYnabtmXPvOTP3Xif2vURY1mvGWFhcWl7JrubW1jc2t/LbO80kSrnLGm7kR/zKsRPmeyFrCE/47CrmzA4cn7Wc0bmMt24YT7wovBTjmHUDexB6fc+1BVGdYiewxdDpT5zp4XW+YJUstcx5UNagAL1qUf4FHfQQwUWKAAwhBGEfNhJ62ijDQkxcFxPiOCFPxRmmyJE2pSxGGTaxI/oOaNfWbEh76ZkotUun+PRyUpo4IE1EeZywPM1U8VQ5S/Y374nylHcb09/RXgGxAkNi/9LNMv+rk7UI9HGqavCoplgxsjpXu6SqK/Lm5peqBDnExEncozgn7CrlrM+m0iSqdtlbW8XfVKZk5d7VuSne5S1pwOWf45wHzUqpfFSq1I8L1TM96iz2sI8izfMEVVyghgZ5x3jEE56NunFr3Bn3n6lGRmt28W0ZDx+nGpPw

AAACz3icjVHLTsJAFD3UF+ILdemmkZjghhQ00SXRjUtIBEyAmGkZoKGvtFMNIRi3/oBb/SvjH+hfeGccEpUYnabtmXPvOTP3Xjvy3ERY1mvGWFhcWl7JrubW1jc2t/LbO80kTGOHN5zQC+MrmyXccwPeEK7w+FUUc+bbHm/Zo3MZb93wOHHD4FKMI9712SBw+67DBFGdYsdnYmj3J2x6eJ0vWCVLLXMelDUoQK9amH9BBz2EcJDCB0cAQdgDQ0JPG2VYiIjrYkJcTMhVcY4pcqRNKYtTBiN2RN8B7dqaDWgvPROldugUj96YlCYOSBNSXkxYnmaqeKqcJfub90R5yruN6W9rL59YgSGxf+lmmf/VyVoE+jhVNbhUU6QYWZ2jXVLVFXlz80tVghwi4iTuUTwm7CjlrM+m0iSqdtlbpuJvKlOycu/o3BTv8pY04PLPcc6DZqVUPipV6seF6pkedRZ72EeR5nmCKi5QQ4O8IzziCc9G3bg17oz7z1QjozW7+LaMhw+kuJPv symmetries. C

AAACxHicjVHLSsNAFD2Nr1pfVZdugkVwVZIq6LIoiMsW7ANqkWQ6rdHJg8lEKEV/wK1+m/gH+hfeGVNQi+iEJGfOvefM3Hv9RASpcpzXgjU3v7C4VFwurayurW+UN7faaZxJxlssFrHs+l7KRRDxlgqU4N1Eci/0Be/4t6c63rnjMg3i6EKNE94PvVEUDAPmKaKaN1flilN1zLJngZuDCvLViMsvuMQAMRgyhOCIoAgLeEjp6cGFg4S4PibESUKBiXPco0TajLI4ZXjE3tJ3RLtezka0156pUTM6RdArSWljjzQx5UnC+jTbxDPjrNnfvCfGU99tTH8/9wqJVbgm9i/dNPO/Ol2LwhDHpoaAakoMo6tjuUtmuqJvbn+pSpFDQpzGA4pLwswop322jSY1teveeib+ZjI1q/csz83wrm9JA3Z/jnMWtGtV96Baax5W6if5qIvYwS72aZ5HqOMcDbSM9yOe8GydWcJKrewz1Srkmm18W9bDB0sqj3M= j kAAACxHicjVHLSsNAFD2Nr1pfVZdugkVwVZIq6LIoiMsW7ANqkWQ6rUMnD5KJUIr+gFv9NvEP9C+8M6agFtEJSc6ce8+Zuff6sRSpcpzXgrWwuLS8Ulwtra1vbG6Vt3faaZQljLdYJKOk63splyLkLSWU5N044V7gS97xx+c63rnjSSqi8EpNYt4PvFEohoJ5iqjm+KZccaqOWfY8cHNQQb4aUfkF1xggAkOGABwhFGEJDyk9PbhwEBPXx5S4hJAwcY57lEibURanDI/YMX1HtOvlbEh77ZkaNaNTJL0JKW0ckCaivISwPs028cw4a/Y376nx1Heb0N/PvQJiFW6J/Us3y/yvTteiMMSpqUFQTbFhdHUsd8lMV/TN7S9VKXKIidN4QPGEMDPKWZ9to0lN7bq3nom/mUzN6j3LczO861vSgN2f45wH7VrVParWmseV+lk+6iL2sI9DmucJ6rhEAy3j/YgnPFsXlrRSK/tMtQq5ZhfflvXwAU2Kj3Q= The rest of the paper is organized as follows: In Sec.II, AAACxHicjVHLSsNAFD2Nr1pfVZdugkVwVZIq6LIoiMsW7ANqkWQ6rdHJg8lEKEV/wK1+m/gH+hfeGVNQi+iEJGfOvefM3Hv9RASpcpzXgjU3v7C4VFwurayurW+UN7faaZxJxlssFrHs+l7KRRDxlgqU4N1Eci/0Be/4t6c63rnjMg3i6EKNE94PvVEUDAPmKaKaN1flilN1zLJngZuDCvLViMsvuMQAMRgyhOCIoAgLeEjp6cGFg4S4PibESUKBiXPco0TajLI4ZXjE3tJ3RLtezka0156pUTM6RdArSWljjzQx5UnC+jTbxDPjrNnfvCfGU99tTH8/9wqJVbgm9i/dNPO/Ol2LwhDHpoaAakoMo6tjuUtmuqJvbn+pSpFDQpzGA4pLwswop322jSY1teveeib+ZjI1q/csz83wrm9JA3Z/jnMWtGtV96Baax5W6if5qIvYwS72aZ5HqOMcDbSM9yOe8GydWcJKrewz1Srkmm18W9bDB0sqj3M= j building on several exactly solvable models of SPT phases B lAAACxHicjVHLSsNAFD2Nr1pfVZdugkVwVZIq6LIoiMsW7ANqkWQ6raGTB5OJUIr+gFv9NvEP9C+8M6agFtEJSc6ce8+Zuff6iQhS5TivBWthcWl5pbhaWlvf2Nwqb++00ziTjLdYLGLZ9b2UiyDiLRUowbuJ5F7oC97xx+c63rnjMg3i6EpNEt4PvVEUDAPmKaKa4qZccaqOWfY8cHNQQb4acfkF1xggBkOGEBwRFGEBDyk9PbhwkBDXx5Q4SSgwcY57lEibURanDI/YMX1HtOvlbER77ZkaNaNTBL2SlDYOSBNTniSsT7NNPDPOmv3Ne2o89d0m9Pdzr5BYhVti/9LNMv+r07UoDHFqagiopsQwujqWu2SmK/rm9peqFDkkxGk8oLgkzIxy1mfbaFJTu+6tZ+JvJlOzes/y3Azv+pY0YPfnOOdBu1Z1j6q15nGlfpaPuog97OOQ5nmCOi7RQMt4P+IJz9aFJazUyj5TrUKu2cW3ZT18AE/qj3U= [63, 64], we construct a family of one-parameter Hamiltoni- y A ans, which interpolate between trivial and topological SPT x phases. In Sec. III we present a sign-problem free quantum Monte Carlo method used to numerically study these models FIG. 2. The Hamiltonian (2.3) defined on the triangular lattice is and discuss the physical observables used to probe the vari- made of two distinct parts. (a) The first one (2.1) describes a trivial ous emergent phases and phase transitions. We present our gapped Ising paramagnet with onsite terms. (b) The second part (2.2) numerical results and discuss their physical interpretation in describes a topological gapped Ising paramagnet protected by the Sec.IV. Lastly, we summarize our findings and highlight fu- Z2 Ising spin-flip symmetry, with plaquette terms involving the six ture research directions in Sec.V. nearest neighbors of a given lattice site 푗. See Eqs. (2.4) and (2.6) for the microscopic definitions of the plaquette operator considered in Z3 this work. (c) Different sublattices for the 2 model on the triangular II. MODELS lattice geometry considered in this work.

In this section we present two single-parameter Hamiltoni- Z3 3 A. 2 model ans, admitting Z2 and Z2 symmetries, that interpolate between trivial and symmetry-protected topological phases. The de- 3 grees of freedom of our models are Ising spins residing on We first consider an SPT phase protected by a Z2 sym- the sites of a two-dimensional triangular lattice, see Fig.2.A metry corresponding to flipping all the spins belonging to trivial Ising paramagnet is simply defined by the Hamiltonian, each one of the three sublattices of the triangular lattice, see Fig.2 (c). The associated generators of this symmetry are ∑︁ 푥 Î 푥 Htri = − 휎 , (2.1) 퐺 / / = 휎 . 푗 푗 A B C 푖∈A/B/C 푖 The plaquette operator is defined as [4, 64, 65], 푥,푦,푧 푧 where 휎 are the standard Pauli matrices (휎 = ±1) de- 3 푗 푗 (Z ) Ö 1 ( − 푧 )( − 푧 ) 2 4 1 휎푘 1 휎푙 fined on site 푗, see Fig.2 (a). The above Hamiltonian has 휗푗 = (−1) . (2.4) 4 푗푘푙 a unique gapped ground state, given by the product state Î 푥 푥 To evaluate the above product, one counts the number of |Ψgsi = 푗 |휎푗 = +1i in the 휎 basis. On the same lattice, we can also define nontrivial SPTs with nearest-neighbor spin pairs belonging to the plaquette sur- − parent Hamiltonians taking the form [3–5, 63–65]:, rounding 푗 and both taking the value 1. If the number of such pairs is odd the product equals −1, and otherwise it ∑︁ 푥 (Z3) Htopo = − 휎 휗푗 . (2.2) 2 푗 푗 equals +1. With the above definition for 휗푗 , we can re- ( 3) Z2 Here, 휗푗 is a plaquette operator centered around site 푗 and in- late Htopo to the trivial paramagnet, Htri, by the following volving all its neighbors, as sketched in Fig.2 (b). It is diagonal unitary transformation, 푧 푧 푧 in the computational basis |휎 , 휎 , . . . 휎 i with eigenvalues 3 † 1 2 푁 (Z ) ( 3) ( 3) 2 Z2 Z2 ±1 . In a geometry with periodic boundary conditions and with Htopo = U Htri U ,. (2.5) appropriate choices of 휗푗 , as we detail in the following sub- (Z3) 푁 sections, the above Hamiltonian has a unique gapped ground Here, U 2 = (−1) 4−−− is a diagonal (in the computational state realizing a symmetry protected topological phase. basis) unitary operator, where 푁4−−− counts the number of Since the ground states of (2.1) and (2.2) describe a trivial triangles with three −1 spins in a given basis configuration. (Z3) and nontrivial SPT phases, respectively, then a quantum phase H 2 3 Importantly, topo commutes with the Z2 symmetry. The transition is expected to occur along a symmetry-preserving resulting SPT phase corresponds to a type-III cocycle and path in parameter space that interpolates between the two. therefore couples nontrivially to the three different sublattice Note that the ground states of Htopo and −Htopo have the same Ising symmetries. Gauging the Z3 symmetry gives rise to a strong SPT index, but a different weak zero-dimensional SPT 2 non-Abelian 퐷4 quantum double phase [64, 65, 67–70]. index [66]. (Z3) (Z3) † Since U 2 = (U 2 ) , we find that for 훼 ≥ 0, H(훼) is This enables us to study two distinct transitions for each related by a unitary transformation to H(1 − 훼). Furthermore, symmetry group, with the following Hamiltonian, since −Htopo is related to Htri by a unitary transformation with ( 3) Z2 Î 푧 H(훼) = 1 − |훼| Htri + 훼Htopo, (2.3) unitary operator U 푖 휎푖 , we also have a duality relating H(훼) to H (−1 − 훼) for 훼 ≤ 0. In other words, the phase with 훼 ∈ [−1, +1]. The transition for 훼 > 0 involves changing diagram is symmetric about 훼 = 0.5 for 훼 > 0 and around only the strong SPT index. The transition for 훼 < 0 involves 훼 = −0.5 for 훼 < 0. an additional change in the weak zero-dimensional SPT index. The above relations are a key property of our model since, In the following we give explicit expressions for the plaquette as we show below, in the computational 휎푧 basis, the Hamil- operators 휗푗 of Eq. (2.2). tonian is sign-problem free only within the parameter range AAAB9XicbVBNS8NAEN34WetX1aOXxSJ4KokU9Fj04rGC/YAmlsl2ky7d3YTdjVJi/4cXD4p49b9489+4bXPQ1gcDj/dmmJkXppxp47rfzsrq2vrGZmmrvL2zu7dfOThs6yRThLZIwhPVDUFTziRtGWY47aaKggg57YSj66nfeaBKs0TemXFKAwGxZBEjYKx0/+RrFgvwFciY036l6tbcGfAy8QpSRQWa/cqXP0hIJqg0hIPWPc9NTZCDMoxwOin7maYpkBHEtGepBEF1kM+unuBTqwxwlChb0uCZ+nsiB6H1WIS2U4AZ6kVvKv7n9TITXQY5k2lmqCTzRVHGsUnwNAI8YIoSw8eWAFHM3orJEBQQY4Mq2xC8xZeXSfu85tVr9dt6tXFVxFFCx+gEnSEPXaAGukFN1EIEKfSMXtGb8+i8OO/Ox7x1xSlmjtAfOJ8/7RySyw==

AAAB83icbVDLSgNBEOyNrxhfUY9eBoPgKeyGgB6DXjxGMImQLGF20psMmX0w0yuEkN/w4kERr/6MN//GSbIHTSxoKKq6Z7orSJU05LrfTmFjc2t7p7hb2ts/ODwqH5+0TZJpgS2RqEQ/BtygkjG2SJLCx1QjjwKFnWB8O/c7T6iNTOIHmqToR3wYy1AKTlbqKU4kBTIjCfvlilt1F2DrxMtJBXI0++Wv3iARWYQxCcWN6XpuSv6Ua/ukwlmplxlMuRjzIXYtjXmExp8udp6xC6sMWJhoWzGxhfp7YsojYyZRYDsjTiOz6s3F/7xuRuG1P5VxmhHGYvlRmClGCZsHwAZSoyA1sYQLLe2uTIy45oJsTCUbgrd68jpp16pevVq/r1UaN3kcRTiDc7gED66gAXfQhBYISOEZXuHNyZwX5935WLYWnHzmFP7A+fwBDpeRsg==

AAAB6nicbVBNS8NAEJ3Ur1q/qh69LBbBU0hKQY9FL96saGuhDWWz3bRLN5uwOxFK6U/w4kERr/4ib/4bt20O2vpg4PHeDDPzwlQKg5737RTW1jc2t4rbpZ3dvf2D8uFRyySZZrzJEpnodkgNl0LxJgqUvJ1qTuNQ8sdwdD3zH5+4NiJRDzhOeRDTgRKRYBStdH+bur1yxXO9Ocgq8XNSgRyNXvmr209YFnOFTFJjOr6XYjChGgWTfFrqZoanlI3ogHcsVTTmJpjMT52SM6v0SZRoWwrJXP09MaGxMeM4tJ0xxaFZ9mbif14nw+gymAiVZsgVWyyKMkkwIbO/SV9ozlCOLaFMC3srYUOqKUObTsmG4C+/vEpaVdevubW7aqV+lcdRhBM4hXPw4QLqcAMNaAKDATzDK7w50nlx3p2PRWvByWeO4Q+czx/kGo2K

AAAB63icbVBNS8NAEJ3Ur1q/qh69LBbBU0lKQY9FLx4r2g9oQ9lsN+3S3U3Y3Qgl9C948aCIV/+QN/+NmzQHbX0w8Hhvhpl5QcyZNq777ZQ2Nre2d8q7lb39g8Oj6vFJV0eJIrRDIh6pfoA15UzSjmGG036sKBYBp71gdpv5vSeqNIvko5nH1Bd4IlnICDaZ9BAzOarW3LqbA60TryA1KNAeVb+G44gkgkpDONZ64Lmx8VOsDCOcLirDRNMYkxme0IGlEguq/TS/dYEurDJGYaRsSYNy9fdEioXWcxHYToHNVK96mfifN0hMeO2nTMaJoZIsF4UJRyZC2eNozBQlhs8twUQxeysiU6wwMTaeig3BW315nXQbda9Zb943aq2bIo4ynME5XIIHV9CCO2hDBwhM4Rle4c0Rzovz7nwsW0tOMXMKf+B8/gAQiI5B

훼 ∈ [−0.5, 0.5]. We can, therefore, determine all physical p Spin properties within this range by means of quantum Monte Carlo M simulations and treat the rest of the phase diagram using the −1 above duality relations. M − 1 +1 Op. … AAAB+3icdVDLSsNAFJ34rPUV69LNYCu4Kklpqu6KbtxZwT6gDWUyvWmHTiZhZiItpb/ixoUibv0Rd/6N04egomd1OOde7rknSDhT2nE+rJXVtfWNzcxWdntnd2/fPsg1VJxKCnUa81i2AqKAMwF1zTSHViKBRAGHZjC8mvnNe5CKxeJOjxPwI9IXLGSUaCN17dxNApLoWGImejDChaTQtfNO0XMct1LBc+I5niHlC69ScrFbdObIoyVqXfu904tpGoHQlBOl2q6TaH9CpGaUwzTbSRUkhA5JH9qGChKB8ifz7FN8YpQeDk2AMBYaz9XvGxMSKTWOAjMZET1Qv72Z+JfXTnV47k+YSFINgi4OhSnHOsazInCPSaCajw0hVDKTFdMBkYRqU1fWlPD1Kf6fNEpFt1ws35by1ctlHRl0hI7RKXLRGaqia1RDdUTRCD2gJ/RsTa1H68V6XYyuWMudQ/QD1tsnU6+T/A== Operator index X B. Z model 2 I 1 The second model we consider is an SPT phase protected by a Z2 symmetry corresponding to a global Ising spin flip, 0 푧 푧 lattice site | i 휎푖 → −휎푖 ∀푖. In this case, the plaquette operator takes the form [3–5, 63], FIG. 3. Example of a quantum Monte Carlo configuration for the (Z ) Ö 1 (− − 푧 푧 ) 2 2 1 휎푘 휎푙 model (3.1) in the SSE formulation. It consists of a basis state 휗푗 = 푖 . (2.6) 4 푗푘푙 |휎i ≡ |휎1휎2 . . . 휎푁 i (at index 푝 = 0) and a sequence 푆푀 of operators of fixed length 푀. There are only two types of operators in our The above product gives a minus sign if the number of nearest case: Identity 퐼 and spin-flip 푋. Here, the identity operators have neighbor spin pairs belonging to the plaquette surrounding 푗 been assigned a lattice site. The algorithm aims at sampling the and pointing in the same direction equals to 2 or 6, and gives a configuration space {|휎i, 푆푀 }. plus sign otherwise. Such pairs can only come in even numbers (0, 2, 4, 6), such that 휗(Z2) is a Hermitian operator despite the 푗 Within the stochastic series expansion (SSE) formulation of (Z2) presence of the imaginary number 푖 in its definition. Htopo is quantum Monte Carlo, the partition function of the system at related to the trivial paramagnet Htri by the following unitary inverse temperature 훽 = 1/푇 reads [71, 72], transformation, 푀 ∑︁ ∑︁ (−훽)푛 (푀 − 푛)! D Ö E † ( ) Z = 휎 푂(푝) 휎 , (3.2) H Z2 = U (Z2) H U (Z2) . (2.7) 푀! topo tri {|휎i} {푆푀 } 푝=0 where the configuration space is defined by all possible combi- Here, U (Z2) = (−1) 푁DW , where 푁 counts the number of DW nations of basis states {|휎i ≡ |휎 , 휎 , . . . 휎 i} and sequence domain walls in a given spin configuration. The number of 1 2 푁 of operators {푆 } of fixed length 푀. For a given sequence of domain walls is well-defined on the triangular lattice since each 푀 operators 푆 , the operators are denoted by 푂(푝) with position spin configuration defines a unique configuration of closed and 푀 index 푝. They can either be the identity 퐼(푝) or 푋(푝), as de- nonintersecting domain walls on the dual honeycomb lattice. fined in Eq. (3.1), see Fig.3. Spin-flip operators 푋(푝) act on We note that U (Z2) is diagonal in the computational basis. lattice sites 푗 ∈ [1, 푁] but its index can be omitted because it Gauging the Ising symmetry in this model realizes the double is actually redundant with the sequence 푆 considered, unless semion phase [63]. 푀 specified otherwise. The integer 푛 is the number of nonidentity By the same argument as in the Z3 case, we find a duality 2 operators in the sequence 푆 . Note that 푀 should be taken relating 훼 and (1 − 훼) for 훼 ≥ 0, and relating 훼 and (−1 − 훼) 푀 large enough such that 푛 < 푀 is ensured in practice. One can for 훼 ≤ 0. As before, owing to this symmetry, and we can rewrite the partition function (3.2) as limit ourselves to calculations in the range 훼 ∈ [−0.5, 0.5] using a sign problem free quantum Monte Carlo simulations ∑︁ ∑︁ Z = 푊 |휎i, 푆 , (3.3) and obtain results for the rest of the phase diagram using the 푀 {| i} { } above duality relations. 휎 푆푀 where 푊 (|휎i, 푆푀 ) is the weight of a configuration with a probability 푃(|휎i, 푆푀 ) = 푊 (|휎i, 푆푀 )/Z. We want to sam- III. NUMERICAL METHODS: QUANTUM MONTE ple these configurations in a Monte Carlo fashion, which sup- CARLO AND OBSERVABLES poses that 푊 (|휎i, 푆푀 ) ≥ 0 for all configurations, otherwise we end up with the infamous sign problem [62]. Because the A. Stochastic Series expansion number 푛 of spin-flip operators (3.1) is even to respect the periodicity along the “operator index axis” (see Fig.3), this condition is fulfilled for all the models considered in the range For convenience, and without loss of generality, we first 훼 ∈ [−0.5, +0.5]. rewrite the Hamiltonian (2.3) as There are no known efficient loop or cluster-type updates [73–76] for the models (3.1), and we can only rely on ∑︁푁 푥 H = 푋 푗 with 푋 푗 = −휎 1 − |훼| + 훼휗푗 , (3.1) 푗=1 푗 local moves in the configuration space [77]. This limits the system sizes one can access in practice to a few hundred lat- with 푁 being the total number of lattice sites. In what follows, tice sites. Assuming some valid configuration is defined by we always impose periodic boundary conditions. {|휎i, 푆푀 }, the updates that we propose involve changes in the 4 sequence of operators 푆푀 that will indirectly involve changes as the Z2 model for 훼 < 0. The remaining case of the Z2 model in the basis state |휎i. There cannot be two operators with with 훼 > 0 has already been thoroughly studied by the same the same index 푝, and there can only be an even number of authors [59]. For completeness, we provide a brief account of 푋 operators at a given lattice site 푗; otherwise, the two states its physical properties. |휎i sandwiching the product of operators in (3.2) would be We begin our analysis by probing the equal-time (휏 = 0) two- different. Full details on the algorithm implementation are point correlation function (3.4), depicted in Figs.4 (a, b, c) for discussed in AppendixA. 12 × 12 lattices, at representative values of 훼, at which order Because of the constraints in the different models as 훼 → sets in. In all cases, we observe a clear maximum of intensity ±1/2 (some spin flips are strictly prohibited at 훼 = ±1/2, at specific wave vectors 풒order indicating the presence of long- as explained below), we have found that the SSE algorithm range magnetic order, as reported in TableI. The maximum with local updates gives incorrect results when one gets very of intensity serves as the definition for the order parameter close to 훼 = ±1/2 (by comparing to exact diagonalization on (squared) that can be systematically analyzed versus system small system sizes). We believe this is an ergodicity issue in size and 훼. For a continuous phase transition, the following the SSE configuration space due to the nature of our updates. finite-size critical scaling is expected [2, 81], Note that this problem does not concern the data shown in = / 2훽/휈 h 1/휈 i this work since they are relatively far away from 훼 1 2. 퐶 풒order, 0 × 퐿 = G 퐿 (훼 − 훼c) , (4.1) However, a study of the Z2 model at 훼 = 1/2 was necessary in Ref. 59. We have, therefore, developed a complementary with 훽 being the order parameter critical exponent, 휈 the cor- algorithm specifically to study that case. This algorithm is relation length critical exponent, 퐿 the linear system size, 훼c based on projective quantum Monte Carlo [78] and does not the critical point, and G a universal scaling function. Based suffer from ergodicity issues, see AppendixB. on Eq. (4.1) and as explained in AppendixC, we determine the position of the critical points 훼c and critical exponents 2훽/휈. The results of this analysis are shown in Fig.4 (d, e, f). Indeed, B. Physical observables we find that after rescaling curves corresponding to different systems sizes cross at a single point. For all the three cases To determine the different quantum phases and phase tran- considered we estimate 2훽/휈 ' 1. sition in each model, we focus primarily on the spin structure We now turn our attention to estimating the gap Δ between factor associated with the imaginary-time two-point correla- the two Z2 Ising symmetry sectors by a numerical fit of the tion function, imaginary-time correlation function to Eq. (3.4) at 풒order close to criticality 훼 ' 훼c. The relevant data set is shown in Ap- 1 ∑︁ −푖풒·풓 D 푧 푧 E 퐶 풒, 휏 = e 휎풓 (휏)휎0 (0) . (3.4) pendixD. For a continuous transition, the expected finite-size 푁 풓 scaling of the gap follows the form [2, 81], This quantity is readily computed in SSE simulation, since it −푧 is diagonal in the computational basis [71, 79, 80]. The equal Δ 퐿, 훼c ∼ 퐿 , (4.2) time correlation function (휏 = 0) probes potential spontaneous magnetic order marked by a peak at an ordering wave vector with 푧 being the dynamical critical exponent. In Fig.5, we show the gap versus 1/퐿, which displays a linear scaling, 풒 ≡ 풒order. The gap Δ between the two Z2 Ising symmetric sectors is also accessible by examining the long imaginary compatible with 푧 = 1 for each of the models. time asymptotic decay In principle, we can independently extract the correlation length exponent 휈 by rescaling the 푥 axis of Fig.4 (d, e, f) → 1/휈 ( − ) 퐶(풒order, 휏) ∝ exp(−Δ휏) (3.5) using the scaling argument 훼 퐿 훼 훼c . Doing so is expected to result in a curve collapse associated with the scal- with 휏 ∈ [0, 훽/2]. Periodic systems of finite size 푁 = 퐿 × 퐿 ing function G in Eq. (4.1). However, the limited system sizes are considered, with the lattice geometry of Fig.2. We set the accessible numerically do not allow for a reliable estimation inverse temperature of the SSE algorithm at 훽 = 4퐿, which of 휈. Thus, we cannot confidently determine the precise uni- we found to be sufficiently low to probe the ground state of the versality classes of each phase transition. In the following, models studied. we discuss several possible scenarios describing the various phase transitions, based on our numerical observations and symmetry arguments. IV. NUMERICAL RESULTS

Z3 There are a total of four cases to be investigated since there A. The 2 model with 훼 > 0 3 are two diﬀerent symmetry groups (Z2 and Z2) and 훼 ≶ 0.A quantum phase transition is expected in all cases with either a The peaks at 풒order in the structure factor, see TableI, are single transition point separating the two phases or a two-step consistent with a ferromagnetic phase within each sublattice transition via an intermediate phase breaking the protecting A, B, and C as displayed in Fig.6 (see also AppendixE). One symmetry. We ﬁrst turn our attention to the nature of the way to understand the emergence of this phase is to study the 3 / transition for the Z2 model at positive and negative 훼 as well Hamiltonian at 훼 = 1 2. At that value, it is easy to see that AAAB6nicbVDLTgJBEOzFF+IL9ehlIjHxRHYNRo9ELx4xyiMBQmaHXpgwO7vOzBrJhk/w4kFjvPpF3vwbB9iDgpV0UqnqTneXHwuujet+O7mV1bX1jfxmYWt7Z3evuH/Q0FGiGNZZJCLV8qlGwSXWDTcCW7FCGvoCm/7oeuo3H1FpHsl7M46xG9KB5AFn1Fjp7qH31CuW3LI7A1kmXkZKkKHWK351+hFLQpSGCap123Nj002pMpwJnBQ6icaYshEdYNtSSUPU3XR26oScWKVPgkjZkobM1N8TKQ21Hoe+7QypGepFbyr+57UTE1x2Uy7jxKBk80VBIoiJyPRv0ucKmRFjSyhT3N5K2JAqyoxNp2BD8BZfXiaNs7JXKZ/fVkrVqyyOPBzBMZyCBxdQhRuoQR0YDOAZXuHNEc6L8+58zFtzTjZzCH/gfP4AbzqN6A==

AAAB6nicbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0WPRi8eK9gPaUDbbTbt0s4m7E6GE/gQvHhTx6i/y5r9x2+agrQ8GHu/NMDMvSKQw6LrfTmFldW19o7hZ2tre2d0r7x80TZxqxhsslrFuB9RwKRRvoEDJ24nmNAokbwWjm6nfeuLaiFg94DjhfkQHSoSCUbTS/WNv3CtX3Ko7A1kmXk4qkKPeK391+zFLI66QSWpMx3MT9DOqUTDJJ6VuanhC2YgOeMdSRSNu/Gx26oScWKVPwljbUkhm6u+JjEbGjKPAdkYUh2bRm4r/eZ0Uwys/EypJkSs2XxSmkmBMpn+TvtCcoRxbQpkW9lbChlRThjadkg3BW3x5mTTPqt559eLuvFK7zuMowhEcwyl4cAk1uIU6NIDBAJ7hFd4c6bw4787HvLXg5DOH8AfO5w9wvo3p

AAAB6nicbVDLSgNBEOyNrxhfUY9eBoPgKexKRC9C0IvHiOYByRJmJ51kyOzsMjMrxCWf4MWDIl79Im/+jZNkD5pY0FBUddPdFcSCa+O6305uZXVtfSO/Wdja3tndK+4fNHSUKIZ1FolItQKqUXCJdcONwFaskIaBwGYwupn6zUdUmkfywYxj9EM6kLzPGTVWun+68rrFklt2ZyDLxMtICTLUusWvTi9iSYjSMEG1bntubPyUKsOZwEmhk2iMKRvRAbYtlTRE7aezUyfkxCo90o+ULWnITP09kdJQ63EY2M6QmqFe9Kbif147Mf1LP+UyTgxKNl/UTwQxEZn+TXpcITNibAllittbCRtSRZmx6RRsCN7iy8ukcVb2KuXzu0qpep3FkYcjOIZT8OACqnALNagDgwE8wyu8OcJ5cd6dj3lrzslmDuEPnM8f3ZuNiA==

AAAB8nicbVBNS8NAEJ3Ur1q/qh69LBbBU0lKQY9FLx4r2FZIQ9lsN+3SzSbsToQS+jO8eFDEq7/Gm//GbZuDtj4YeLw3w8y8MJXCoOt+O6WNza3tnfJuZW//4PCoenzSNUmmGe+wRCb6MaSGS6F4BwVK/phqTuNQ8l44uZ37vSeujUjUA05THsR0pEQkGEUr+V6jjyLmhniNQbXm1t0FyDrxClKDAu1B9as/TFgWc4VMUmN8z00xyKlGwSSfVfqZ4SllEzrivqWK2j1Bvjh5Ri6sMiRRom0pJAv190ROY2OmcWg7Y4pjs+rNxf88P8PoOsiFSjPkii0XRZkkmJD5/2QoNGcop5ZQpoW9lbAx1ZShTaliQ/BWX14n3Ubda9ab981a66aIowxncA6X4MEVtOAO2tABBgk8wyu8Oei8OO/Ox7K15BQzp/AHzucP2XiQUw==

AAACBnicbVDLSsNAFJ34rPUVdSlCsAiuSlIruiy6cVnBPrCNYTKZtENnMmFmIpQQN278FTcuFHHrN7jzb5ymWWjrgYEz59zLvff4MSVS2fa3sbC4tLyyWlorr29sbm2bO7ttyROBcAtxykXXhxJTEuGWIoribiwwZD7FHX90OfE791hIwqMbNY6xy+AgIiFBUGnJMw/6DKqh76e3mVe7O3nIv4KljAeYZp5Zsat2DmueOAWpgAJNz/zqBxwlDEcKUShlz7Fj5aZQKIIozsr9ROIYohEc4J6mEWRYuml+RmYdaSWwQi70i5SVq787UsikHDNfV062lLPeRPzP6yUqPHdTEsWJwhGaDgoTailuTTKxAiIwUnSsCUSC6F0tNIQCIqWTK+sQnNmT50m7VnXq1dPreqVxUcRRAvvgEBwDB5yBBrgCTdACCDyCZ/AK3own48V4Nz6mpQtG0bMH/sD4/AEG6pl/

AAAB8nicbVDLSgMxFL1TX7W+qi7dBIvgqs5IRZdFNy5cVLAPmA4lk2ba0EwSkoxQSj/DjQtF3Po17vwb03YWWj0QOJxzL7nnxIozY33/yyusrK6tbxQ3S1vbO7t75f2DlpGZJrRJJJe6E2NDORO0aZnltKM0xWnMaTse3cz89iPVhknxYMeKRikeCJYwgq2Twq7SUlmJgrO7XrniV/050F8S5KQCORq98me3L0mWUmEJx8aEga9sNMHaMsLptNTNDFWYjPCAho4KnFITTeYnT9GJU/ookdo9YdFc/bkxwakx4zR2kym2Q7PszcT/vDCzyVU0YUJllgqy+CjJOHIhZ/lRn2lKLB87golm7lZEhlhjYl1LJVdCsBz5L2mdV4Na9eK+Vqlf53UU4QiO4RQCuIQ63EIDmkBAwhO8wKtnvWfvzXtfjBa8fOcQfsH7+AZ7fpC8

AAACAHicbVDLSsNAFJ3UV62vqAsXbgaLUDclKQUVXBTduKxgH9iEMJlO2qGTSZiZCCXEhb/ixoUibv0Md/6NkzYLrR64cDjnXu69x48ZlcqyvozS0vLK6lp5vbKxubW9Y+7udWWUCEw6OGKR6PtIEkY56SiqGOnHgqDQZ6TnT65yv3dPhKQRv1XTmLghGnEaUIyUljzzwAmRGvt+epd5jYeag1g8RhfWiWdWrbo1A/xL7IJUQYG2Z346wwgnIeEKMyTlwLZi5aZIKIoZySpOIkmM8ASNyEBTjkIi3XT2QAaPtTKEQSR0cQVn6s+JFIVSTkNfd+bnykUvF//zBokKztyU8jhRhOP5oiBhUEUwTwMOqSBYsakmCAuqb4V4jATCSmdW0SHYiy//Jd1G3W7Wz2+a1dZlEUcZHIIjUAM2OAUtcA3aoAMwyMATeAGvxqPxbLwZ7/PWklHM7INfMD6+Ac1xlec=

AAAB8XicbVBNS8NAEJ34WetX1aOXxSJ4CokU1INQ9OKxgv3ANpTJdtMu3WzC7kYopf/CiwdFvPpvvPlv3LY5aOuDgcd7M8zMC1PBtfG8b2dldW19Y7OwVdze2d3bLx0cNnSSKcrqNBGJaoWomeCS1Q03grVSxTAOBWuGw9up33xiSvNEPphRyoIY+5JHnKKx0mMHRTrAa8+tdEtlz/VmIMvEz0kZctS6pa9OL6FZzKShArVu+15qgjEqw6lgk2In0yxFOsQ+a1sqMWY6GM8unpBTq/RIlChb0pCZ+ntijLHWozi0nTGagV70puJ/Xjsz0WUw5jLNDJN0vijKBDEJmb5PelwxasTIEqSK21sJHaBCamxIRRuCv/jyMmmcu37FvbqvlKs3eRwFOIYTOAMfLqAKd1CDOlCQ8Ayv8OZo58V5dz7mrStOPnMEf+B8/gBnzJAb

AAAB83icbVBNSwMxEM3Wr1q/qh69BIvgxWVXCtaDUPTisYL9gO5SZtNsG5rNhiQrlNK/4cWDIl79M978N6btHrT1wcDjvRlm5kWSM20879sprK1vbG4Vt0s7u3v7B+XDo5ZOM0Vok6Q8VZ0INOVM0KZhhtOOVBSSiNN2NLqb+e0nqjRLxaMZSxomMBAsZgSMlYIAuBzCzYXnVmu9csVzvTnwKvFzUkE5Gr3yV9BPSZZQYQgHrbu+J004AWUY4XRaCjJNJZARDGjXUgEJ1eFkfvMUn1mlj+NU2RIGz9XfExNItB4nke1MwAz1sjcT//O6mYlr4YQJmRkqyGJRnHFsUjwLAPeZosTwsSVAFLO3YjIEBcTYmEo2BH/55VXSunT9qnv9UK3Ub/M4iugEnaJz5KMrVEf3qIGaiCCJntErenMy58V5dz4WrQUnnzlGf+B8/gBRYZCU

AAAB8nicbVBNS8NAEN34WetX1aOXxSJ4MSRSUA9C0YvHCvYD0lAm2027dLMbdjdCCf0ZXjwo4tVf481/47bNQVsfDDzem2FmXpRypo3nfTsrq2vrG5ulrfL2zu7efuXgsKVlpghtEsml6kSgKWeCNg0znHZSRSGJOG1Ho7up336iSjMpHs04pWECA8FiRsBYKegCT4dwc+65tV6l6rneDHiZ+AWpogKNXuWr25ckS6gwhIPWge+lJsxBGUY4nZS7maYpkBEMaGCpgITqMJ+dPMGnVunjWCpbwuCZ+nsih0TrcRLZzgTMUC96U/E/L8hMfBXmTKSZoYLMF8UZx0bi6f+4zxQlho8tAaKYvRWTISggxqZUtiH4iy8vk9aF69fc64datX5bxFFCx+gEnSEfXaI6ukcN1EQESfSMXtGbY5wX5935mLeuOMXMEfoD5/MH1ByQUg==

AAACBHicbVDLSsNAFJ34rPUVddnNYBFclaQU1F3RjcsK9oFtCJPJpB06k4SZiVBCBDf+ihsXirj1I9z5N07SLLT1wMCZc+7l3nu8mFGpLOvbWFldW9/YrGxVt3d29/bNg8OejBKBSRdHLBIDD0nCaEi6iipGBrEgiHuM9L3pVe7374mQNApv1SwmDkfjkAYUI6Ul16yNOFITz0vvMrf5UHwET3nkE5a5Zt1qWAXgMrFLUgclOq75NfIjnHASKsyQlEPbipWTIqEoZiSrjhJJYoSnaEyGmoaIE+mkxREZPNGKD4NI6BcqWKi/O1LEpZxxT1fmW8pFLxf/84aJCs6dlIZxokiI54OChEEVwTwR6FNBsGIzTRAWVO8K8QQJhJXOrapDsBdPXia9ZsNuNS5uWvX2ZRlHBdTAMTgFNjgDbXANOqALMHgEz+AVvBlPxovxbnzMS1eMsucI/IHx+QPPqZje

AAACBHicbVDLSsNAFJ34rPUVddnNYBHqwpKUYnUhFN24rGAf0IQwmU7aoTNJmJkIJXThxl9x40IRt36EO//GSZuFth64cDjnXu69x48Zlcqyvo2V1bX1jc3CVnF7Z3dv3zw47MgoEZi0ccQi0fORJIyGpK2oYqQXC4K4z0jXH99kfveBCEmj8F5NYuJyNAxpQDFSWvLMkoNYPEKew5EaCZ7i6dWZVa2fNyq1U88sW1VrBrhM7JyUQY6WZ345gwgnnIQKMyRl37Zi5aZIKIoZmRadRJIY4TEakr6mIeJEuunsiSk80coABpHQFSo4U39PpIhLOeG+7sxulYteJv7n9RMVXLgpDeNEkRDPFwUJgyqCWSJwQAXBik00QVhQfSvEIyQQVjq3og7BXnx5mXRqVbtevbyrl5vXeRwFUALHoAJs0ABNcAtaoA0weATP4BW8GU/Gi/FufMxbV4x85gj8gfH5A3yPlrI=

AAACA3icbVDLSsNAFJ34rPUVdaebYBHqJiTaUl0IRTcuK9gHNCFMptN26MwkzEyEEgpu/BU3LhRx60+482+ctFlo64ELh3Pu5d57wpgSqRzn21haXlldWy9sFDe3tnd2zb39lowSgXATRTQSnRBKTAnHTUUUxZ1YYMhCitvh6Cbz2w9YSBLxezWOsc/ggJM+QVBpKTAPPUjjIQw8BtVQsBRNrhz7vFYtV08Ds+TYzhTWInFzUgI5GoH55fUilDDMFaJQyq7rxMpPoVAEUTwpeonEMUQjOMBdTTlkWPrp9IeJdaKVntWPhC6urKn6eyKFTMoxC3Vndqqc9zLxP6+bqP6FnxIeJwpzNFvUT6ilIisLxOoRgZGiY00gEkTfaqEhFBApHVtRh+DOv7xIWme2W7Ev7yql+nUeRwEcgWNQBi6ogTq4BQ3QBAg8gmfwCt6MJ+PFeDc+Zq1LRj5zAP7A+PwBDACWfA==

AAAB+3icbVBNS8NAEJ3Ur1q/Yj16CRahXmISCupBKHrxWMF+QBPKZrtpl242YXcjltK/4sWDIl79I978N27bHLT1wcDjvRlm5oUpo1I5zrdRWFvf2Nwqbpd2dvf2D8zDcksmmcCkiROWiE6IJGGUk6aiipFOKgiKQ0ba4eh25rcfiZA04Q9qnJIgRgNOI4qR0lLPLHt+SBQ693l27dqOV/XOembFsZ05rFXi5qQCORo988vvJziLCVeYISm7rpOqYIKEopiRacnPJEkRHqEB6WrKUUxkMJnfPrVOtdK3okTo4sqaq78nJiiWchyHujNGaiiXvZn4n9fNVHQZTChPM0U4XiyKMmapxJoFYfWpIFixsSYIC6pvtfAQCYSVjqukQ3CXX14lLc92a/bVfa1Sv8njKMIxnEAVXLiAOtxBA5qA4Qme4RXejKnxYrwbH4vWgpHPHMEfGJ8/DICShA==

AAAB+nicbVBNS8NAEN3Ur1q/Uj16CRahXmJSCtqDUPTisYL9gCaUzXbaLt1swu5GKbE/xYsHRbz6S7z5b9y2OWjrg4HHezPMzAtiRqVynG8jt7a+sbmV3y7s7O7tH5jFw5aMEkGgSSIWiU6AJTDKoamoYtCJBeAwYNAOxjczv/0AQtKI36tJDH6Ih5wOKMFKSz2zWPECUPjc48mVY9fK7lnPLDm2M4e1StyMlFCGRs/88voRSULgijAsZdd1YuWnWChKGEwLXiIhxmSMh9DVlOMQpJ/OT59ap1rpW4NI6OLKmqu/J1IcSjkJA90ZYjWSy95M/M/rJmpw6aeUx4kCThaLBgmzVGTNcrD6VABRbKIJJoLqWy0ywgITpdMq6BDc5ZdXSatiu1W7dlct1a+zOPLoGJ2gMnLRBaqjW9RATUTQI3pGr+jNeDJejHfjY9GaM7KZI/QHxucPoT2STw==

AAAB+3icbZDLSsNAFIYnXmu9xbp0EyxC3dSkFKwLoejGZQV7gSaUyXTSDp1MwsyJWEJfxY0LRdz6Iu58GydtFtr6w8DHf87hnPn9mDMFtv1trK1vbG5tF3aKu3v7B4fmUamjokQS2iYRj2TPx4pyJmgbGHDaiyXFoc9p15/cZvXuI5WKReIBpjH1QjwSLGAEg7YGZqnm+hTwhSuSa6dq25XG+cAs25oyWavg5FBGuVoD88sdRiQJqQDCsVJ9x47BS7EERjidFd1E0RiTCR7RvkaBQ6q8dH77zDrTztAKIqmfAGvu/p5IcajUNPR1Z4hhrJZrmflfrZ9A0PBSJuIEqCCLRUHCLYisLAhryCQlwKcaMJFM32qRMZaYgI6rqENwlr+8Cp1a1alXr+7r5eZNHkcBnaBTVEEOukRNdIdaqI0IekLP6BW9GTPjxXg3Phata0Y+c4z+yPj8ARKQkog=

AAACBHicbVDLSsNAFJ34rPUVddnNYBHqwpLUQnUhFN24rGAf0IQwmU7aoTNJmJkIJXThxl9x40IRt36EO//GSZuFth64cDjnXu69x48Zlcqyvo2V1bX1jc3CVnF7Z3dv3zw47MgoEZi0ccQi0fORJIyGpK2oYqQXC4K4z0jXH99kfveBCEmj8F5NYuJyNAxpQDFSWvLMkoNYPEKew5EaCZ7i6dWZVT1v1Cu1U88sW1VrBrhM7JyUQY6WZ345gwgnnIQKMyRl37Zi5aZIKIoZmRadRJIY4TEakr6mIeJEuunsiSk80coABpHQFSo4U39PpIhLOeG+7sxulYteJv7n9RMVXLgpDeNEkRDPFwUJgyqCWSJwQAXBik00QVhQfSvEIyQQVjq3og7BXnx5mXRqVbtevbyrl5vXeRwFUALHoAJs0ABNcAtaoA0weATP4BW8GU/Gi/FufMxbV4x85gj8gfH5A3f5lq8=

AAAB9HicbVBNS8NAEN3Ur1q/qh69LBbBU0hsS70IRS8eK9gPaEOZbDft0s0m7m4KJfR3ePGgiFd/jDf/jds2B219MPB4b4aZeX7MmdKO823lNja3tnfyu4W9/YPDo+LxSUtFiSS0SSIeyY4PinImaFMzzWknlhRCn9O2P76b++0JlYpF4lFPY+qFMBQsYAS0kbwe8HgEN45drpWr/WLJsZ0F8DpxM1JCGRr94ldvEJEkpEITDkp1XSfWXgpSM8LprNBLFI2BjGFIu4YKCKny0sXRM3xhlAEOImlKaLxQf0+kECo1DX3TGYIeqVVvLv7ndRMdXHspE3GiqSDLRUHCsY7wPAE8YJISzaeGAJHM3IrJCCQQbXIqmBDc1ZfXSevKdit29aFSqt9mceTRGTpHl8hFNVRH96iBmoigJ/SMXtGbNbFerHfrY9mas7KZU/QH1ucPzySQ0w==

AAAB9HicbVDJSgNBEK1xjXGLevTSGAQvDjMhLhch6MVjBLNAMoSaTk/SpGexuycQhnyHFw+KePVjvPk3dpaDJj4oeLxXRVU9PxFcacf5tlZW19Y3NnNb+e2d3b39wsFhXcWppKxGYxHLpo+KCR6xmuZasGYiGYa+YA1/cDfxG0MmFY+jRz1KmBdiL+IBp6iN5LVRJH28OXfs8mWpUyg6tjMFWSbunBRhjmqn8NXuxjQNWaSpQKVarpNoL0OpORVsnG+niiVIB9hjLUMjDJnysunRY3JqlC4JYmkq0mSq/p7IMFRqFPqmM0TdV4veRPzPa6U6uPYyHiWpZhGdLQpSQXRMJgmQLpeMajEyBKnk5lZC+yiRapNT3oTgLr68TOol2y3bFw/lYuV2HkcOjuEEzsCFK6jAPVShBhSe4Ble4c0aWi/Wu/Uxa12x5jNH8AfW5w/BZ5DK

AAAB9XicbVBNS8NAEJ3Ur1q/qh69LBbBiyGxhepBKHrxWMF+QBvLZrtpl242YXejlND/4cWDIl79L978N27bHLT1wcDjvRlm5vkxZ0o7zreVW1ldW9/Ibxa2tnd294r7B00VJZLQBol4JNs+VpQzQRuaaU7bsaQ49Dlt+aObqd96pFKxSNzrcUy9EA8ECxjB2kgPXczjIb46c+xy1S33iiXHdmZAy8TNSAky1HvFr24/IklIhSYcK9VxnVh7KZaaEU4nhW6iaIzJCA9ox1CBQ6q8dHb1BJ0YpY+CSJoSGs3U3xMpDpUah77pDLEeqkVvKv7ndRIdXHgpE3GiqSDzRUHCkY7QNALUZ5ISzceGYCKZuRWRIZaYaBNUwYTgLr68TJrntluxL+8qpdp1FkcejuAYTsGFKtTgFurQAAISnuEV3qwn68V6tz7mrTkrmzmEP7A+fwA3RpEK

AAAB83icbVBNS8NAEJ3Ur1q/qh69BItQLyUpgnorevFYwX5AE8pmu2mXbjZhdyKW0L/hxYMiXv0z3vw3btsctPXBwOO9GWbmBYngGh3n2yqsrW9sbhW3Szu7e/sH5cOjto5TRVmLxiJW3YBoJrhkLeQoWDdRjESBYJ1gfDvzO49MaR7LB5wkzI/IUPKQU4JG8qoesicMwmwwPe+XK07NmcNeJW5OKpCj2S9/eYOYphGTSAXRuuc6CfoZUcipYNOSl2qWEDomQ9YzVJKIaT+b3zy1z4wysMNYmZJoz9XfExmJtJ5EgemMCI70sjcT//N6KYZXfsZlkiKTdLEoTIWNsT0LwB5wxSiKiSGEKm5utemIKELRxFQyIbjLL6+Sdr3mXtSu7+uVxk0eRxFO4BSq4MIlNOAOmtACCgk8wyu8Wan1Yr1bH4vWgpXPHMMfWJ8/1meRkg== AAAB83icbVBNS8NAEN3Ur1q/qh69BItQLyUpgnorevFYwX5AE8pmO2mXbjZhdyKW0L/hxYMiXv0z3vw3btsctPXBwOO9GWbmBYngGh3n2yqsrW9sbhW3Szu7e/sH5cOjto5TxaDFYhGrbkA1CC6hhRwFdBMFNAoEdILx7czvPILSPJYPOEnAj+hQ8pAzikbyqh7CEwZhBtPzfrni1Jw57FXi5qRCcjT75S9vELM0AolMUK17rpOgn1GFnAmYlrxUQ0LZmA6hZ6ikEWg/m988tc+MMrDDWJmSaM/V3xMZjbSeRIHpjCiO9LI3E//zeimGV37GZZIiSLZYFKbCxtieBWAPuAKGYmIIZYqbW202oooyNDGVTAju8surpF2vuRe16/t6pXGTx1EkJ+SUVIlLLkmD3JEmaRFGEvJMXsmblVov1rv1sWgtWPnMMfkD6/MH1+2Rkw== AAAB83icbVBNS8NAEJ3Ur1q/qh69BItQLyUpgnorevFYwX5AE8pmu2mXbjZhdyKW0L/hxYMiXv0z3vw3btsctPXBwOO9GWbmBYngGh3n2yqsrW9sbhW3Szu7e/sH5cOjto5TRVmLxiJW3YBoJrhkLeQoWDdRjESBYJ1gfDvzO49MaR7LB5wkzI/IUPKQU4JG8qoesicMwiycnvfLFafmzGGvEjcnFcjR7Je/vEFM04hJpIJo3XOdBP2MKORUsGnJSzVLCB2TIesZKknEtJ/Nb57aZ0YZ2GGsTEm05+rviYxEWk+iwHRGBEd62ZuJ/3m9FMMrP+MySZFJulgUpsLG2J4FYA+4YhTFxBBCFTe32nREFKFoYiqZENzll1dJu15zL2rX9/VK4yaPowgncApVcOESGnAHTWgBhQSe4RXerNR6sd6tj0VrwcpnjuEPrM8f2XORlA==

AAAB+HicbZDLSsNAFIZP6q3WS6Mu3QwWoYKUpBR0WXTjsoK9QBvKZDpph04mYWYi1NAnceNCEbc+ijvfxkmbhbb+MPDxn3M4Z34/5kxpx/m2ChubW9s7xd3S3v7BYdk+Ou6oKJGEtknEI9nzsaKcCdrWTHPaiyXFoc9p15/eZvXuI5WKReJBz2LqhXgsWMAI1sYa2uWBz8a86lw6GciLoV1xas5CaB3cHCqQqzW0vwajiCQhFZpwrFTfdWLtpVhqRjidlwaJojEmUzymfYMCh1R56eLwOTo3zggFkTRPaLRwf0+kOFRqFvqmM8R6olZrmflfrZ/o4NpLmYgTTQVZLgoSjnSEshTQiElKNJ8ZwEQycysiEywx0SarkgnBXf3yOnTqNbdRa9w3Ks2bPI4inMIZVMGFK2jCHbSgDQQSeIZXeLOerBfr3fpYthasfOYE/sj6/AEokpIf

AAAB/nicbVDLSsNAFL3xWesrKq7cDBahgtREA7osunFZwT6gCWUynbRDJw9mJkIJBX/FjQtF3Pod7vwbJ20W2nrgwplz7mXuPX7CmVSW9W0sLa+srq2XNsqbW9s7u+befkvGqSC0SWIei46PJeUsok3FFKedRFAc+py2/dFt7rcfqZAsjh7UOKFeiAcRCxjBSks989D12YBXHTdh55dnyMqf4rRnVqyaNQVaJHZBKlCg0TO/3H5M0pBGinAsZde2EuVlWChGOJ2U3VTSBJMRHtCuphEOqfSy6foTdKKVPgpioStSaKr+nshwKOU49HVniNVQznu5+J/XTVVw7WUsSlJFIzL7KEg5UjHKs0B9JihRfKwJJoLpXREZYoGJ0omVdQj2/MmLpHVRs52ac+9U6jdFHCU4gmOogg1XUIc7aEATCGTwDK/wZjwZL8a78TFrXTKKmQP4A+PzB8hDlBY=

AAACCHicbVDLSsNAFJ3UV62vqEsXDhahgtTEBnRZdOOygn1AE8pkOmmHTpJxZiKU0KUbf8WNC0Xc+gnu/BsnbRbaeuDC4Zx7ufcenzMqlWV9G4Wl5ZXVteJ6aWNza3vH3N1ryTgRmDRxzGLR8ZEkjEakqahipMMFQaHPSNsfXWd++4EISePoTo058UI0iGhAMVJa6pmHrk8HrGKdujx0XE7Paq68FyqtTTJDnPTMslW1poCLxM5JGeRo9Mwvtx/jJCSRwgxJ2bUtrrwUCUUxI5OSm0jCER6hAelqGqGQSC+dPjKBx1rpwyAWuiIFp+rviRSFUo5DX3eGSA3lvJeJ/3ndRAWXXkojnigS4dmiIGFQxTBLBfapIFixsSYIC6pvhXiIBMJKZ1fSIdjzLy+S1nnVdqrOrVOuX+VxFMEBOAIVYIMLUAc3oAGaAINH8AxewZvxZLwY78bHrLVg5DP74A+Mzx9KE5jk

AAACEnicbVDLSsNAFJ3UV62vqEs3g0VoQWrSFnRZdOOygn1AE8pkOmmHTpJxZiKU0G9w46+4caGIW1fu/BsnbUBtPTBwOOdc7tzjcUalsqwvI7eyura+kd8sbG3v7O6Z+wdtGcUCkxaOWCS6HpKE0ZC0FFWMdLkgKPAY6Xjjq9Tv3BMhaRTeqgknboCGIfUpRkpLfbPseHTISg4PYNXh9Kx2+kMdeSdUUpumEVHum0WrYs0Al4mdkSLI0Oybn84gwnFAQoUZkrJnW1y5CRKKYkamBSeWhCM8RkPS0zREAZFuMjtpCk+0MoB+JPQLFZypvycSFEg5CTydDJAayUUvFf/zerHyL9yEhjxWJMTzRX7MoIpg2g8cUEGwYhNNEBZU/xXiERIIK91iQZdgL568TNrVil2v1G/qxcZlVkceHIFjUAI2OAcNcA2aoAUweABP4AW8Go/Gs/FmvM+jOSObOQR/YHx8A9WUnFg=

AAACEnicbVDLSgMxFM3UV62vUZdugkVoQepMW9Bl0Y3LCvYBnaFk0rQNTWZikhHK0G9w46+4caGIW1fu/Bsz7Sy0eiBwOOdcbu4JBKNKO86XlVtZXVvfyG8WtrZ3dvfs/YO2imKJSQtHLJLdACnCaEhammpGukISxANGOsHkKvU790QqGoW3eiqIz9EopEOKkTZS3y57AR2xkic4rHqCntVOPS4y6qk7qZPaLI3Ict8uOhVnDviXuBkpggzNvv3pDSIccxJqzJBSPdcR2k+Q1BQzMit4sSIC4QkakZ6hIeJE+cn8pBk8McoADiNpXqjhXP05kSCu1JQHJsmRHqtlLxX/83qxHl74CQ1FrEmIF4uGMYM6gmk/cEAlwZpNDUFYUvNXiMdIIqxNiwVTgrt88l/SrlbceqV+Uy82LrM68uAIHIMScME5aIBr0AQtgMEDeAIv4NV6tJ6tN+t9Ec1Z2cwh+AXr4xvVkZxY

AAACBHicbVBNS8NAEJ3Ur1q/oh57CRahgtSkVvRY9OKxgv2AJpTNdtMu3Wzi7kYooQcv/hUvHhTx6o/w5r9x0/ag1QcDj/dmmJnnx4xKZdtfRm5peWV1Lb9e2Njc2t4xd/daMkoEJk0csUh0fCQJo5w0FVWMdGJBUOgz0vZHV5nfvidC0ojfqnFMvBANOA0oRkpLPbPo+nTAyvZx1Y3piSvvhEpPJ5kojnpmya7YU1h/iTMnJZij0TM/3X6Ek5BwhRmSsuvYsfJSJBTFjEwKbiJJjPAIDUhXU45CIr10+sTEOtRK3woioYsra6r+nEhRKOU49HVniNRQLnqZ+J/XTVRw4aWUx4kiHM8WBQmzVGRliVh9KghWbKwJwoLqWy08RAJhpXMr6BCcxZf/kla14tQqZze1Uv1yHkceinAAZXDgHOpwDQ1oAoYHeIIXeDUejWfjzXifteaM+cw+/ILx8Q1ifpdP

AAACCHicbVDLSgMxFM3UV62vUZcuHCxCBakzWtFl0Y3LCvYBnVIyaaYNTTIxyQhl6NKNv+LGhSJu/QR3/o2ZdhbaeuByD+fcS3JPIChR2nW/rdzC4tLySn61sLa+sbllb+80VBRLhOsoopFsBVBhSjiua6IpbgmJIQsobgbD69RvPmCpSMTv9EjgDoN9TkKCoDZS1973A9KnJV+QY18w0058dS91cjZODXnUtYtu2Z3AmSdeRoogQ61rf/m9CMUMc40oVKrtuUJ3Eig1QRSPC36ssIBoCPu4bSiHDKtOMjlk7BwapeeEkTTFtTNRf28kkCk1YoGZZFAP1KyXiv957ViHl52EcBFrzNH0oTCmjo6cNBWnRyRGmo4MgUgS81cHDaCESJvsCiYEb/bkedI4LXuV8vltpVi9yuLIgz1wAErAAxegCm5ADdQBAo/gGbyCN+vJerHerY/paM7KdnbBH1ifP0ZOmYM=

AAACCnicbVDLSsNAFJ34rPUVdelmtAgVpCZtQZdFNy4r2Ac0oUymk3boZBJnJkIJXbvxV9y4UMStX+DOv3GSZqGtBy73cM69zNzjRYxKZVnfxtLyyuraemGjuLm1vbNr7u23ZRgLTFo4ZKHoekgSRjlpKaoY6UaCoMBjpOONr1O/80CEpCG/U5OIuAEacupTjJSW+uaR49EhK1ediJ7XzmDWHXkvVFKbppY47Zslq2JlgIvEzkkJ5Gj2zS9nEOI4IFxhhqTs2Vak3AQJRTEj06ITSxIhPEZD0tOUo4BIN8lOmcITrQygHwpdXMFM/b2RoEDKSeDpyQCpkZz3UvE/rxcr/9JNKI9iRTiePeTHDKoQprnAARUEKzbRBGFB9V8hHiGBsNLpFXUI9vzJi6Rdrdj1Sv22Xmpc5XEUwCE4BmVggwvQADegCVoAg0fwDF7Bm/FkvBjvxsdsdMnIdw7AHxifPxQEmUM=

AAAB63icbVDLSgMxFL3xWeur6tJNsAiuyoxUdFl047KCfUA7lEyaaUOTzJBkhDL0F9y4UMStP+TOvzHTzkJbD4QczrmXe+8JE8GN9bxvtLa+sbm1Xdop7+7tHxxWjo7bJk41ZS0ai1h3Q2KY4Iq1LLeCdRPNiAwF64STu9zvPDFteKwe7TRhgSQjxSNOic0lr+ZdDSpV982BV4lfkCoUaA4qX/1hTFPJlKWCGNPzvcQGGdGWU8Fm5X5qWELohIxYz1FFJDNBNt91hs+dMsRRrN1TFs/V3x0ZkcZMZegqJbFjs+zl4n9eL7XRTZBxlaSWKboYFKUC2xjnh+Mh14xaMXWEUM3drpiOiSbUunjKLgR/+eRV0r6s+fVa/aFebdwWcZTgFM7gAny4hgbcQxNaQGEMz/AKb0iiF/SOPhala6joOYE/QJ8/ygiNbA==

AAAB7HicbVBNTwIxEJ3FL8Qv1KOXRmLiiewSjB6JXjxi4gIJbEi3zEJDt7tpuyaE8Bu8eNAYr/4gb/4bC+xBwZc0fXlvJjPzwlRwbVz32ylsbG5t7xR3S3v7B4dH5eOTlk4yxdBniUhUJ6QaBZfoG24EdlKFNA4FtsPx3dxvP6HSPJGPZpJiENOh5BFn1FjJd6tu7apfrth/AbJOvJxUIEezX/7qDRKWxSgNE1TrruemJphSZTgTOCv1Mo0pZWM6xK6lksaog+li2Rm5sMqARImyTxqyUH93TGms9SQObWVMzUivenPxP6+bmegmmHKZZgYlWw6KMkFMQuaXkwFXyIyYWEKZ4nZXwkZUUWZsPiUbgrd68jpp1apevVp/qFcat3kcRTiDc7gED66hAffQBB8YcHiGV3hzpPPivDsfy9KCk/ecwh84nz87iY2o

AAAB6nicbVBNSwMxEJ2tX7V+VT16CRbBU9mVQj0WvXisaD+gXUo2nW1Ds9klyQpl6U/w4kERr/4ib/4b03YP2vog4fHeDDPzgkRwbVz32ylsbG5t7xR3S3v7B4dH5eOTto5TxbDFYhGrbkA1Ci6xZbgR2E0U0igQ2Akmt3O/84RK81g+mmmCfkRHkoecUWOlB7fqDsoV+y9A1omXkwrkaA7KX/1hzNIIpWGCat3z3MT4GVWGM4GzUj/VmFA2oSPsWSpphNrPFqvOyIVVhiSMlX3SkIX6uyOjkdbTKLCVETVjverNxf+8XmrCaz/jMkkNSrYcFKaCmJjM7yZDrpAZMbWEMsXtroSNqaLM2HRKNgRv9eR10r6qerVq7b5WadzkcRThDM7hEjyoQwPuoAktYDCCZ3iFN0c4L86787EsLTh5zyn8gfP5A1Q+jS0=

3 3 (↵ > 0) (↵ < 0) 2 (↵ < 0)

AAAC43icjVHLSsNAFD2Nr1pfVZeCBFuhbkraLnThoujGZQX7QKtlko42NC+SiVBKXblzJ279Abf6L+If6F94Z0zBB6ITkpw5954zc+81A8eOhGG8pLSJyanpmfRsZm5+YXEpu7zSiPw4tHjd8h0/bJks4o7t8bqwhcNbQciZazq8afb3Zbx5ycPI9r0jMQj4qcsuPPvctpggqpNdz7ddJnqmOTwenVU65fxVId9mTtBj+q6R3+pkc0bRUEv/CUoJyCFZNT/7jDa68GEhhgsOD4KwA4aInhOUYCAg7hRD4kJCtopzjJAhbUxZnDIYsX36XtDuJGE92kvPSKktOsWhNySljk3S+JQXEpan6SoeK2fJ/uY9VJ7ybgP6m4mXS6xAj9i/dOPM/+pkLQLn2FE12FRToBhZnZW4xKor8ub6p6oEOQTESdyleEjYUspxn3WliVTtsrdMxV9VpmTl3kpyY7zJW9KAS9/H+RM0ysVSpVg+LOeqe8mo01jDBgo0z21UcYAa6uR9jQc84knj2o12q919pGqpRLOKL0u7fwdu65nl ZAAAC43icjVHLSsNAFD2Nr1pfVZeCBFuhbkraLnQlRTcuK9gHWi2TdLSheZFMhFLqyp07cesPuNV/Ef9A/8I7Ywo+EJ2Q5My595yZe68ZOHYkDOMlpU1MTk3PpGczc/MLi0vZ5ZVG5MehxeuW7/hhy2QRd2yP14UtHN4KQs5c0+FNs78v481LHka27x2JQcBPXXbh2ee2xQRRnex6vu0y0TPN4fHorNIp568K+TZzgh7Td438ViebM4qGWvpPUEpADsmq+dlntNGFDwsxXHB4EIQdMET0nKAEAwFxpxgSFxKyVZxjhAxpY8rilMGI7dP3gnYnCevRXnpGSm3RKQ69ISl1bJLGp7yQsDxNV/FYOUv2N++h8pR3G9DfTLxcYgV6xP6lG2f+VydrETjHjqrBppoCxcjqrMQlVl2RN9c/VSXIISBO4i7FQ8KWUo77rCtNpGqXvWUq/qoyJSv3VpIb403ekgZc+j7On6BRLpYqxfJhOVfdS0adxho2UKB5bqOKA9RQJ+9rPOARTxrXbrRb7e4jVUslmlV8Wdr9O3Oxmec= 2 Z2 (b) Z (a) AAACz3icjVHLTsJAFD3UF+ILdemmkZjghhQ00SXRjUtIBEyAmLYM0NBXplMNIRi3/oBb/SvjH+hfeGccEpUYnabtmXPvOTP3Xif2vURY1mvGWFhcWl7JrubW1jc2t/LbO80kSrnLGm7kR/zKsRPmeyFrCE/47CrmzA4cn7Wc0bmMt24YT7wovBTjmHUDexB6fc+1BVGdYiewxdDpT5zp4XW+YJUstcx5UNagAL1qUf4FHfQQwUWKAAwhBGEfNhJ62ijDQkxcFxPiOCFPxRmmyJE2pSxGGTaxI/oOaNfWbEh76ZkotUun+PRyUpo4IE1EeZywPM1U8VQ5S/Y374nylHcb09/RXgGxAkNi/9LNMv+rk7UI9HGqavCoplgxsjpXu6SqK/Lm5peqBDnExEncozgn7CrlrM+m0iSqdtlbW8XfVKZk5d7VuSne5S1pwOWf45wHzUqpfFSq1I8L1TM96iz2sI8izfMEVVyghgZ5x3jEE56NunFr3Bn3n6lGRmt28W0ZDx+nGpPw AAACz3icjVHLTsJAFD3UF+ILdemmkZjghhQ00SXRjUtIBEyAmGkZoKGvtFMNIRi3/oBb/SvjH+hfeGccEpUYnabtmXPvOTP3Xjvy3ERY1mvGWFhcWl7JrubW1jc2t/LbO80kTGOHN5zQC+MrmyXccwPeEK7w+FUUc+bbHm/Zo3MZb93wOHHD4FKMI9712SBw+67DBFGdYsdnYmj3J2x6eJ0vWCVLLXMelDUoQK9amH9BBz2EcJDCB0cAQdgDQ0JPG2VYiIjrYkJcTMhVcY4pcqRNKYtTBiN2RN8B7dqaDWgvPROldugUj96YlCYOSBNSXkxYnmaqeKqcJfub90R5yruN6W9rL59YgSGxf+lmmf/VyVoE+jhVNbhUU6QYWZ2jXVLVFXlz80tVghwi4iTuUTwm7CjlrM+m0iSqdtlbpuJvKlOycu/o3BTv8pY04PLPcc6DZqVUPipV6seF6pkedRZ72EeR5nmCKi5QQ4O8IzziCc9G3bg17oz7z1QjozW7+LaMhw+kuJPv (c)

AAACz3icjVHLTsJAFD3UF+ILdemmkZjghhQ00SXRjUtIBEyAmHYYoKGvtFMNIRi3/oBb/SvjH+hfeGccEpUYnabtmXPvOTP3Xify3ERY1mvGWFhcWl7JrubW1jc2t/LbO80kTGPGGyz0wvjKsRPuuQFvCFd4/CqKue07Hm85o3MZb93wOHHD4FKMI9717UHg9l1mC6I6xY5vi6HTn7Dp4XW+YJUstcx5UNagAL1qYf4FHfQQgiGFD44AgrAHGwk9bZRhISKuiwlxMSFXxTmmyJE2pSxOGTaxI/oOaNfWbEB76ZkoNaNTPHpjUpo4IE1IeTFheZqp4qlyluxv3hPlKe82pr+jvXxiBYbE/qWbZf5XJ2sR6ONU1eBSTZFiZHVMu6SqK/Lm5peqBDlExEnco3hMmCnlrM+m0iSqdtlbW8XfVKZk5Z7p3BTv8pY04PLPcc6DZqVUPipV6seF6pkedRZ72EeR5nmCKi5QQ4O8IzziCc9G3bg17oz7z1QjozW7+LaMhw+pfJPx

2⇡/3, 2⇡/p3 ↵ =0.4 0, 4⇡/3p3 ↵ = 0.48 0, 2⇡/p3 ↵ = 0.4 ± 12 12 12 12 12 12 ⇥ ⇥ ⇥ 05 .

0.050 05 05 . . 0 0 0, 0 0 0.05 04 . 0 p , ⇡, ⇡/p3 0.05 04 2⇡/3, 2⇡/3 3 04 . . q 0 0 03 . ± ± 0 ± 0 0 03 03 .

0.025 . 0 0 , 02 AAACBnicbVDLSgMxFM3UV62vUZciBItQQcqMFHRZ7MZlBfuAzlAymbQNzSRjkhHK0JUbf8WNC0Xc+g3u/Bsz7Sy09UDI4Zx7ufeeIGZUacf5tgorq2vrG8XN0tb2zu6evX/QViKRmLSwYEJ2A6QIo5y0NNWMdGNJUBQw0gnGjczvPBCpqOB3ehITP0JDTgcUI22kvn3c8AI6ZBUvECxUk8h86f303MlUeda3y07VmQEuEzcnZZCj2be/vFDgJCJcY4aU6rlOrP0USU0xI9OSlygSIzxGQ9IzlKOIKD+dnTGFp0YJ4UBI87iGM/V3R4oila1oKiOkR2rRy8T/vF6iB1d+SnmcaMLxfNAgYVALmGUCQyoJ1mxiCMKSml0hHiGJsDbJlUwI7uLJy6R9UXVr1dptrVy/zuMogiNwAirABZegDm5AE7QABo/gGbyCN+vJerHerY95acHKew7BH1ifP9dUmLo= , C . 4⇡/3, 0 0 q q 0.025 q 02 02

y . .

0 0.025 0 01 . 0 AAACBnicbVDLSgMxFM3UV62vUZciBItQQcqMFHRZ7MZlBfuAzlAymbQNzSRjkhHK0JUbf8WNC0Xc+g3u/Bsz7Sy09UDI4Zx7ufeeIGZUacf5tgorq2vrG8XN0tb2zu6evX/QViKRmLSwYEJ2A6QIo5y0NNWMdGNJUBQw0gnGjczvPBCpqOB3ehITP0JDTgcUI22kvn3c8AI6ZBUvECxUk8h86f303MlUeda3y07VmQEuEzcnZZCj2be/vFDgJCJcY4aU6rlOrP0USU0xI9OSlygSIzxGQ9IzlKOIKD+dnTGFp0YJ4UBI87iGM/V3R4oila1oKiOkR2rRy8T/vF6iB1d+SnmcaMLxfNAgYVALmGUCQyoJ1mxiCMKSml0hHiGJsDbJlUwI7uLJy6R9UXVr1dptrVy/zuMogiNwAirABZegDm5AE7QABo/gGbyCN+vJerHerY95acHKew7BH1ifP9dUmLo= C AAACBnicbVDLSgMxFM3UV62vUZciBItQQcqMFHRZ7MZlBfuAzlAymbQNzSRjkhHK0JUbf8WNC0Xc+g3u/Bsz7Sy09UDI4Zx7ufeeIGZUacf5tgorq2vrG8XN0tb2zu6evX/QViKRmLSwYEJ2A6QIo5y0NNWMdGNJUBQw0gnGjczvPBCpqOB3ehITP0JDTgcUI22kvn3c8AI6ZBUvECxUk8h86f303MlUeda3y07VmQEuEzcnZZCj2be/vFDgJCJcY4aU6rlOrP0USU0xI9OSlygSIzxGQ9IzlKOIKD+dnTGFp0YJ4UBI87iGM/V3R4oila1oKiOkR2rRy8T/vF6iB1d+SnmcaMLxfNAgYVALmGUCQyoJ1mxiCMKSml0hHiGJsDbJlUwI7uLJy6R9UXVr1dptrVy/zuMogiNwAirABZegDm5AE7QABo/gGbyCN+vJerHerY95acHKew7BH1ifP9dUmLo= C 01 01

. p . 0 2⇡/3, 2⇡/3 3 0 00 . qx ± ⌥ 0.00 00 00 . . 0 0.0 0.00

0.6 1.0 0.8 3 3 ⇥ ↵c =0.375(5) ↵c = 0.467(2) 3 3 4 4 6 6 ⇥ ⇥ ⇥ 0.5 6 6 6 6 ⌫ ⇥ 0.00 0.01 0.02 0.03 0.04 0.05 0.00 0.01 0.029 9 0.03 0.04 0.05 0.00 0.01 0.02 0.03 0.04 0.905 9 0.8 ⇥ / ⇥ ⇥ 8 8 12 12 12 12 ⇥ 2 ⇥ 10 10 0.6 15 15 ⇥ ⇥ L ⇥ 15 15 12 12 · ⇥ ⇥ 0.4 0.6 18 18 ⇥ 0

, 2/⌫ =1.02(2) 0.4 2/⌫ =0.9(1) (d) 0.3 (e) 0.4 (f) order q AAACKXicbVBNSysxFM349bR+9enybYJFUJA6Uwq+ZdGNCxc+eFWhU0smc6cGM8mY3HlQhvk7bvwrbhQUdesfeZnahV8XQg7nnEtyTpRJYdH3n72p6ZnZuR/zC7XFpeWV1frPtROrc8Ohy7XU5ixiFqRQ0EWBEs4yAyyNJJxGlweVfvoPjBVa/cVRBv2UDZVIBGfoqEG9cxBGYii3wkjL2I5SdxVX5SBMGV5ESaFNDKbc8SuT2Q55rJEenRetMAJku6HKy0G94Tf98dCvIJiABpnM8aB+H8aa5yko5JJZ2wv8DPsFMyi4hLIW5hYyxi/ZEHoOKpaC7RfjpCXddExME23cUUjH7PuNgqW2SuGcVQL7WavI77RejsnvfiFUliMo/vZQkkuKmla10VgY4ChHDjBuhPsr5RfMMI6u3JorIfgc+Ss4aTWDdrP9p93o7E/qmCe/yAbZIgHZIx1ySI5Jl3ByTW7JA3n0brw778l7ebNOeZOddfJhvNf/EimnuQ== C 0.2 2/⌫ =1.00(8) 0.2 0.2 ↵ = 0.374(2) c 0.0 0.1 0.0 0.2 0.3 0.4 0.5 0.50 0.45 0.40 0.5 0.4 0.3 0.2 ↵AAACyXicjVHLTsJAFD3UF+ILdemmkZi4Ii2a6JLoxsQNJgImQMx0GKDSl+3UiMSVP+BWf8z4B/oX3hlLohKj07Q9c+49Z+be60Sem0jLes0ZM7Nz8wv5xcLS8srqWnF9o5GEacxFnYdeGF84LBGeG4i6dKUnLqJYMN/xRNMZHqt480bEiRsG53IUiY7P+oHbczmTRDXazIsG7LJYssqWXuY0sDNQQrZqYfEFbXQRgiOFD4EAkrAHhoSeFmxYiIjrYExcTMjVcYF7FEibUpagDEbskL592rUyNqC98ky0mtMpHr0xKU3skCakvJiwOs3U8VQ7K/Y377H2VHcb0d/JvHxiJQbE/qWbZP5Xp2qR6OFQ1+BSTZFmVHU8c0l1V9TNzS9VSXKIiFO4S/GYMNfKSZ9NrUl07aq3TMffdKZi1Z5nuSne1S1pwPbPcU6DRqVs75UrZ/ul6lE26jy2sI1dmucBqjhBDXXyvsIjnvBsnBrXxq1x95lq5DLNJr4t4+EDuyWRnQ== ↵AAACyXicjVHLTsJAFD3UF+ILdemmkZi4Ii2a6JLoxsQNJgImQMx0GKDSl+3UiMSVP+BWf8z4B/oX3hlLohKj07Q9c+49Z+be60Sem0jLes0ZM7Nz8wv5xcLS8srqWnF9o5GEacxFnYdeGF84LBGeG4i6dKUnLqJYMN/xRNMZHqt480bEiRsG53IUiY7P+oHbczmTRDXazIsG7LJYssqWXuY0sDNQQrZqYfEFbXQRgiOFD4EAkrAHhoSeFmxYiIjrYExcTMjVcYF7FEibUpagDEbskL592rUyNqC98ky0mtMpHr0xKU3skCakvJiwOs3U8VQ7K/Y377H2VHcb0d/JvHxiJQbE/qWbZP5Xp2qR6OFQ1+BSTZFmVHU8c0l1V9TNzS9VSXKIiFO4S/GYMNfKSZ9NrUl07aq3TMffdKZi1Z5nuSne1S1pwPbPcU6DRqVs75UrZ/ul6lE26jy2sI1dmucBqjhBDXXyvsIjnvBsnBrXxq1x95lq5DLNJr4t4+EDuyWRnQ== ↵AAACyXicjVHLTsJAFD3UF+ILdemmkZi4Ii2a6JLoxsQNJgImQMx0GKDSl+3UiMSVP+BWf8z4B/oX3hlLohKj07Q9c+49Z+be60Sem0jLes0ZM7Nz8wv5xcLS8srqWnF9o5GEacxFnYdeGF84LBGeG4i6dKUnLqJYMN/xRNMZHqt480bEiRsG53IUiY7P+oHbczmTRDXazIsG7LJYssqWXuY0sDNQQrZqYfEFbXQRgiOFD4EAkrAHhoSeFmxYiIjrYExcTMjVcYF7FEibUpagDEbskL592rUyNqC98ky0mtMpHr0xKU3skCakvJiwOs3U8VQ7K/Y377H2VHcb0d/JvHxiJQbE/qWbZP5Xp2qR6OFQ1+BSTZFmVHU8c0l1V9TNzS9VSXKIiFO4S/GYMNfKSZ9NrUl07aq3TMffdKZi1Z5nuSne1S1pwPbPcU6DRqVs75UrZ/ul6lE26jy2sI1dmucBqjhBDXXyvsIjnvBsnBrXxq1x95lq5DLNJr4t4+EDuyWRnQ==

= Z3 Z3 FIG. 4. (a)(b)(c) Equal-time (휏 0) two-point correlation (3.4) for the three models considered, (a) the 2 model with 훼 > 0,(b) the 2 model with 훼 < 0 and (a) the Z2 model with 훼 < 0, with 푁 = 12 × 12 and inverse temperature 훽 = 48, represented as an intensity map in the Brillouin zone of the triangular lattice. A representative value of 훼 has been chosen for each model. Intensity peaks are visible at specific wave vector 풒order reported in TableI, and 퐶(풒order, 0) serves as the definition for a potential order parameter (squared), studied versus 훼.(d)(e)(f) ( ) Z3 Z3 Z Finite-size scaling analysis of 퐶 풒order, 0 for (d) the 2 model with 훼 > 0,(e) the 2 model with 훼 < 0 and (f) the 2 model with 훼 < 0. For a second order phase transition it is expected to follow the finite-size critical scaling of Eq. (4.1). In the three different cases, we find a single crossing point of the different system sizes with 2훽/휈 ' 1, as reported in TableI (see AppendixC for information on how we estimated 2훽/휈 and 훼c). The data corresponds to the average over the different symmetry-related 풒order wave vectors.

1.0

3 ↵ =0.3735 Z2 model ↵ = 0.462 0.8 model ↵ = 0.3713 Z2 0.6 z =1 1/L / 0.4 FIG. 6. Magnetization for one of the ground states of the FM3 phase, Z3 0.2 which appears for the 2 model for 훼 > 0. The other ground states are obtained by flipping all the spins on one given sublattice (each sublattice is shown in a different color). 0.0 0.0 0.1 0.2 0.3 0.4 1/L certain spin flips become strictly disallowed, thereby creating a kinetic constraint. The disallowed spin flips are the ones that

FIG. 5. For the three models considered, Δ gap versus inverse linear would change the parity of 푁4−−− (the number of triangles with system size 퐿 at 훼 ' 훼c, extracted by an exponential fit of the three −1 spins). The magnetic order close to 훼 = 1/2 should, imaginary-time correlation function (3.4). For a second order phase therefore, try to maximize the number of flippable spins. It is transition, it is expected to scale according to Eq. (4.2). For the three easy to see that configurations for which each sublattice forms models, the linear scaling observed versus 1/퐿 is compatible with a a perfect ferromagnet have all of their spins flippable. dynamical exponent 푧 = 1, when excluding the smallest system size. 3 Due to the Z2 symmetry, there is no direct coupling (e.g., 푧 푧 휎A 휎B ) between the FM order parameters on each sublattice: They are only coupled through their energy density (like in 6

change the parity of the number of domain walls. It is easy Model Z3 ( 훼 > 0) Z3 ( 훼 < 0) Z ( 훼 < 0) 2 2 2 to see that a stripe phase of a period 2 has all of its spins in (0, 0) ± 2휋 , ± 2√휋 a flippable configuration, which is expected from an energetic 3 3 3 2휋 4휋 0, √ , 0 2휋 2√휋 3 point of view. 풒order 3 ± 3 , ∓ 3 3 휋, ± √휋 Based on symmetry alone, a simultaneous breaking of both 2휋 2√휋 4휋 3 3 , 0, ± √ 3 3 3 Ising and the rotational Z3 symmetry must occur via a first- √ √ 3 Commensurate order transition [90, 91]. In order to reconcile this scenario Order type FM3 × 3 3 stripe with our evidence for a second order phase transition, one would need to invoke a weakly first order transition, with a 훼 0.375(5) −0.467(2) −0.374(2) c correlation length (finite at the transition) that is larger than 2훽/휈 1.00(8) 0.9(1) 1.02(2) the system sizes numerically available. Another possibility is to break Z and Z through two consecutive transitions, with 푧 1 1 1 2 3 an intermediate nematic phase [90, 91]. In that case, the Z2 breaking transition would be expected to be in the Ising class TABLE I. The first row indicates the wave vector 풒order at which and would correspond to the transition we observe. Further magnetic order settles in for the different cases, see Fig.4. The work would be needed to distinguish these two scenarios. next second row labels the different ordered phases. The third row corresponds to our estimates of the position of the transition point 훼c for each model. The next row reports on the value of the combined critical exponents 2훽/휈 obtained from the finite-size scaling analysis D. The Z2 model with 훼 > 0 of the numerical data. See AppendixC for information on how we estimated 2훽/휈 and 훼c. The last row shows the dynamical exponent The Z2 model with 훼 > 0 was thoroughly investigated in 푧 obtained by the closing of the finite-size gap at criticality. Ref. 59, using the same quantum Monte Carlo methods (SSE and projective) that have been developed in this paper. We 2 provide a brief overview of the main results for completeness. the classical Ashkin-Teller model [82] in the case of Z2). If this energy density coupling were irrelevant, we would expect Similarly to what we have uncovered in this work, there three decoupled Ising critical points. However, renormaliza- is also an intermediate phase that spontaneously breaks the tion group calculations show that this coupling is relevant, and protecting Ising Z2 symmetry, and which displays stripe or- drives the transition either to the 푂(푁 = 3) universality if der around the wavevector |풒order| ' 2휋/5. A jump in the order parameter at 훼c ≈ −0.48 suggests a first order transition, 푁 < 푁푐, or to the case of 푂(푁 = 3) with cubic anisotropy if in agreement with a symmetry-based Ginzburg-Landau anal- 푁 > 푁푐 [83–85]. One would need to access larger system sizes to verify this scenario by extracting the critical exponents ysis prohibiting a continuous transition for the corresponding 휈 and 훽 with high precision. 풒order [87]. Remarkably, while one might have expected the intermediate phase to be gapped and confined, it was found to be gapless and dual to a deconfined U(1) gauge theory due to B. The Z3 model with 훼 < 0 the incommensurability of the stripe pattern, providing one of 2 the first observations of the “Cantor deconfinement” scenario in a microscopic model [92–96]. The peaks at 풒order√ in√ the structure factor, see TableI, are consistent with a ( 3 × 3) order within each sublattice [86– 88], see Fig.7. This phase breaks translation and the Z3 rotation symmetry of the triangular lattice on top of the Ising V. CONCLUSION AND PERSPECTIVES spin flip symmetries. Assuming we can neglect the coupling between the sublattices, we would find three independent tran- Employing numerical simulations based on stochastic series sitions in the XY class with Z6 anisotropy [87]. Interestingly, expansion quantum Monte Carlo, we have studied the quantum the Z6 anisotropy is predicted to be dangerously irrelevant, phase diagram of two Hamiltonians interpolating between triv- leading to an emergent U(1) symmetry in the ordered phase 3 ial and nontrivial paramagnets, protected by Z2 and Z2 sym- below a length scale which diverges as a power of the cor- metries, respectively. In all cases, we find that the transition relation length [89]. However, the different sublattices are, happens via an intermediate symmetry-breaking phase, where in reality, coupled and the impact of this coupling should be the protecting symmetry is spontaneously broken, displaying studied in future work. long-range magnetic order. By performing a finite-size scaling analysis of the order parameter, we precisely determined the location of the critical points. The phase diagram of the C. The Z2 model with 훼 < 0 various models that were investigated in this work are summa- rized in Fig.1. Moreover, we computed the gap Δ between the The ordered phase is a stripe phase that breaks translation two Z2 Ising symmetry sectors at criticality, and find that it and the Z3 rotation symmetry of the triangular lattice on top scales as the inverse linear system size of the system, compat- of the Z2 Ising spin flip symmetry, see Fig.8. At 훼 = −1/2, ible with a dynamical exponent 푧 = 1. We also discussed the certain spin flips become disallowed: the ones that would not different possible scenarios describing the nature of the phase 7

FIG. 8. Magnetization for one of the ground states with commensurate stripe order, which appears for the Z2 model for 훼 < 0.

compared to the transition from a trivial paramagnet to the √ √ same symmetry-breaking phase. Previous studies on such FIG. 7. Magnetization for one of the ground states of the ( 3 × 3)3 gapless SPT order, also called symmetry-enriched criticality, Z3 phase, which appears for the 2 model with 훼 < 0. The other ground states are obtained by flipping all the spins on one given sublattice has been mostly limited to one dimension [31–40] (except for (each sublattice is shown in a different color). Each sublattice is itself Refs. [32, 34, 106] which include higher dimensional cases). subdivided into three subsublattices (e.g., sublattice A is divided into The four models presented here provide an ideal platform to AA, AB, and AC), and the magnetization takes a pattern of the type study these phenomena in higher dimensions. 1 1 (1, − 2 , − 2 ) for these three subsublattices (1 is represented as a big Finally, the models presented here also provide a way of 1 arrow while − 2 is represented as a small one). For the particular studying transitions between discrete gauge theories. Whereas ground state shown in the figure, the magnetization has a simple the trivial paramagnet is dual to the toric code, the nontrivial stripe pattern, but this is not the case for all ground states. SPT phases are dual to the double semion (in the case of Z2) 3 and the non-Abelian 퐷4 quantum double (in the case of Z2). This gauge description was particularly useful to study the transitions, which we were not able to single out numerically transition between toric code and double semion in the case of in the present study. 훼 > 0 [59]. A generalization of this gauge theory description 3 Despite the fact that we have developed sign-problem-free to the Z2 case for 훼 < 0 and to the Z2 case is left for future algorithms for the models considered, there is no known ef- work. ficient update for sampling the configuration space such as loop updates or cluster-type updates [73–76]. Therefore, we can only rely on local updates in the configuration space [77], ACKNOWLEDGMENTS limiting the system sizes one can simulate. Accessing larger system sizes is paramount in identifying the exact nature of We are grateful to F. Alet, N. Bultinck, X. Cao, S. Capponi, the transitions taking place in these models, calling for the de- E. Fradkin, B. Kang, A. Paramekanti, D. Poilblanc, F. Poll- velopment of a better-suited quantum Monte Carlo algorithm. mann, P.Pujol, A. W.Sandvik, R. Vasseur, C. Xu and L. Zou for One could get inspired by the recent progress made for quan- interesting discussions. M.D. was supported by the U.S. De- tum dimer models on the square lattice, also displaying strong partment of Energy, Office of Science, Office of Basic Energy geometrical restrictions [97]. Sciences, Materials Sciences and Engineering Division under An important follow-up to this work would be to add terms Contract No. DE-AC02-05-CH11231 through the Scientific which frustrate the different magnetic orders, in order to reach Discovery through Advanced Computing (SciDAC) program multicritical points at which trivial and topological paramag- (KC23DAC Topological and Correlated Matter via Tensor Net- netic phases could potentially have a direct transition. For works and Quantum Monte Carlo). S.G. acknowledges support example, in the case of stripe phases, it might be possible to from the Israel Science Foundation, Grant No. 1686/18. T.S. reach a quantum Lifshitz point at which the stripe wavevector acknowledges the financial support of the Natural Sciences goes continuously to zero [93, 98–104]. Another possibility and Engineering Research Council of Canada (NSERC), in is to reach an instance of deconfined quantum critical points, particular the Discovery Grant [RGPIN-2020-05842], the Ac- which have been predicted to occur at the transition between celerator Supplement [RGPAS-2020-00060], and the Discov- different SPT orders in the presence of continuous symme- ery Launch Supplement [DGECR-2020-00222]. This research tries [30, 45–47, 49, 52–55, 105]. used the Lawrencium computational cluster resource provided In this work, we have focused on bulk properties of the by the IT Division at the Lawrence Berkeley National Labo- system and thus used periodic boundary conditions. However, ratory (Supported by the Director, Office of Science, Office it is important to note that the symmetry of the phase diagram of Basic Energy Sciences, of the U.S. Department of Energy around 훼 = 1/2 for 훼 > 0 and around 훼 = −1/2 for 훼 < 0 under Contract No. DE-AC02-05CH11231). This research applies to bulk properties, but not to edge properties. In fact, also used resources of the National Energy Research Scientific the transitions from a topological paramagnet to a symmetry- Computing Center (NERSC), a U.S. Department of Energy breaking phase is expected to have anomalous edge properties Office of Science User Facility operated under Contract No. 8

DE-AC02-05CH11231. T.S. contributed to this work prior to 1. Select a random site 푗 ∈ [1, 푁]. If the number of joining Amazon. operators in the sequence attached to the selected site 푗 is smaller than two, cancel the update. Otherwise, continue. Appendix A: Practical information regarding the SSE quantum Monte Carlo algorithm 2. In the list of operators attached to the site 푗, select one of them at random. We note it 푂(푝1). 1. Monte Carlo updates 3. Get the number of operators 푁ops acting on site 푗 between 푝1 and the first operator 푋 푗 encountered (the oper- We first discuss the two types of Monte Carlo moves that we ator at 푝1 is excluded from the count and the operator 푋 have implemented in order to sample the configuration space. included). If no operator 푋 is encountered before going They are called the identity and spin-flip updates. back to 푝1, the count runs up to the previous operator to 푂(푝1).

a. Identity update 4. Select at random with probability 1/푁ops an operator with 푝 > 푝1 acting on site 푗. The position of this For convenience, we also assign a real space position index 푗 operator is noted 푝2. to identity operators. The first update is a change of real space 5. To fulfill detailed balance in the selection of 푝 and position 푗 of such an operator to another real space position 푖: 1 푝 , the probability to select them in the configuration 퐼 (푝) → 퐼 (푝). It can be performed as follows, 2 푗 푖 before and after the update should be the same. This 1. Run a loop over each operator of the sequence 푆푀 of is the case in the selection of 푝1 but the probability to the current configuration. select 푝2 depends on the nature of the operator at 푝2. Consider this: If it is an 푋 operator then the probability 2. If the operator 푂(푝) is not an identity operator, we move to select it is 1/푁 . After an update changing 푋 to 퐼, to the next 푝 index. ops the probability to select the identity operator at 푝2 will be 3. If the operator at position 푝 is an identity, we get the site modified. In the current scheme, the probability to select 0 0 푗 on which it acts on. it would be 1/(푁ops + 푁ops) with 푁ops the number of operators between 퐼 푗 (푝2) and the next operator 푋 acting 4. We then select at random a site 푖 ∈ [1, 푁]. on site 푗. To correct this imbalance in the selection of 5. We change with probability one the site on which 퐼(푝) 푝2 if 푂 푗 (푝2) ≡ 푋 푗 (푝2), we cancel this selection with 0 is acting from 푗 to 푖. probability, 푃cancel = 1 − 푁ops/(푁ops + 푁ops). Basically, this move should always be accepted since it does 6. If the selection is not canceled, we suggest the update not change the configuration. Indeed, in the definition (3.2) of according to Eq. (A2). the partition function, the identity operators are not specifically associated to a lattice site. We only assign them a lattice site The probability to accept such an update involves the ratio in the algorithm because it makes it much easier to deal with of the weights of the configurations after “a” and before “b”, = ( ) them, especially in regards to the other update. i.e., 푃accept min 푅b→a, 1 , with,

a 푊 |휎ai, 푆푀 b. Spin-flip update 푅b→a = . (A3) b 푊 |휎bi, 푆푀 The second type of update involves two operators at a time, Specifically, by defining the ratio of matrix elements on different positions 푝1 and 푝2 in the sequence but at the same position 푗 in real space, Î h휎a| a 푂|휎ai 푂∈푆푀 h i h i 푟 = Î ≥ 0, (A4) 푋 (푝 ), 푋 (푝 ) ←→ 퐼 (푝 ), 퐼 (푝 ) , (A1) h휎b| b 푂|휎bi 푗 1 푗 2 푗 1 푗 2 푂∈푆푀 and one finds that the acceptance probability of the updates (A1) h i h i and (A2) is given by 퐼 푗 (푝1), 푋 푗 (푝2) ←→ 푋 푗 (푝1), 퐼 푗 (푝2) . (A2) h i 2 These updates change the configuration and should be accepted 푅푋 푋→퐼 퐼 = 푟 (푀 − 푛 + 2)(푀 − 푛 + 1) / 푁 훽 , or refused fulfilling detailed balance. In between 푝 and 푝 h i 1 2 2/ ( − − )( − ) (A5) at the real space position 푗, there can be as many identity 푅퐼 퐼 →푋 푋 = 푟 푁 훽 푀 푛 1 푀 푛 , operators 퐼 푗 as we want but no 푋 푗 operators, otherwise these 푅퐼 푋→푋 퐼 = 푅푋 퐼 →퐼 푋 = 푟, updates would lead to nonvalid configurations. Note the “periodic boundary condition” along the 푝 axis, as shown in Fig.3. with 푛 the number of nonidentity operators before the update. This update can be performed as follows: The 푁2 factor comes from the fact that we label the identity AAAB8HicbVBNS8NAEJ34WetX1aOXYBE8laQI6kEoevFYwX5IG8pmu2mX7m7C7kQoob/CiwdFvPpzvPlv3LY5aOuDgcd7M8zMCxPBDXret7Oyura+sVnYKm7v7O7tlw4OmyZONWUNGotYt0NimOCKNZCjYO1EMyJDwVrh6Hbqt56YNjxWDzhOWCDJQPGIU4JWehRMDXB4XZW9UtmreDO4y8TPSRly1Hulr24/pqlkCqkgxnR8L8EgIxo5FWxS7KaGJYSOyIB1LFVEMhNks4Mn7qlV+m4Ua1sK3Zn6eyIj0pixDG2nJDg0i95U/M/rpBhdBhlXSYpM0fmiKBUuxu70e7fPNaMoxpYQqrm91aVDoglFm1HRhuAvvrxMmtWKf165uq+Wazd5HAU4hhM4Ax8uoAZ3UIcGUJDwDK/w5mjnxXl3PuatK04+cwR/4Hz+AIQxkDw=

AAAB63icdVDLSsNAFJ34rPVVdelmsAiuQhJNH7uiG5cV+oI2lMl00g6dmYSZiVBCf8GNC0Xc+kPu/BsnbQUVPXDhcM693HtPmDCqtON8WGvrG5tb24Wd4u7e/sFh6ei4o+JUYtLGMYtlL0SKMCpIW1PNSC+RBPGQkW44vcn97j2RisaipWcJCTgaCxpRjHQutRBlw1LZsesV3/Vr0LEvK36t6hri1T3fc6BrOwuUwQrNYel9MIpxyonQmCGl+q6T6CBDUlPMyLw4SBVJEJ6iMekbKhAnKsgWt87huVFGMIqlKaHhQv0+kSGu1IyHppMjPVG/vVz8y+unOqoFGRVJqonAy0VRyqCOYf44HFFJsGYzQxCW1NwK8QRJhLWJp2hC+PoU/k86nu1e2fU7r9y4XsVRAKfgDFwAF1RBA9yCJmgDDCbgATyBZ4tbj9aL9bpsXbNWMyfgB6y3T3Ptjoo=

AAAB63icdVDLSsNAFJ34rPVVdelmsAiuQhJNH7uimy4r2Ae0oUwmk3bozCTMTIQS+gtuXCji1h9y5984fQgqeuDC4Zx7ufeeMGVUacf5sNbWNza3tgs7xd29/YPD0tFxRyWZxKSNE5bIXogUYVSQtqaakV4qCeIhI91wcjP3u/dEKpqIOz1NScDRSNCYYqTnUpOgaFgqO3a94rt+DTr2ZcWvVV1DvLrnew50bWeBMlihNSy9D6IEZ5wIjRlSqu86qQ5yJDXFjMyKg0yRFOEJGpG+oQJxooJ8cesMnhslgnEiTQkNF+r3iRxxpaY8NJ0c6bH67c3Fv7x+puNakFORZpoIvFwUZwzqBM4fhxGVBGs2NQRhSc2tEI+RRFibeIomhK9P4f+k49nulV2/9cqN61UcBXAKzsAFcEEVNEATtEAbYDAGD+AJPFvcerRerNdl65q1mjkBP2C9fQJPaY5y

AAAB9XicdVDLSsNAFJ3UV62vqks3g0VwFZLa9LErunFZwT6giWUynaRDJ5MwM1FK7H+4caGIW//FnX/j9CGo6IELh3Pu5d57/IRRqSzrw8itrK6tb+Q3C1vbO7t7xf2DjoxTgUkbxywWPR9JwignbUUVI71EEBT5jHT98cXM794SIWnMr9UkIV6EQk4DipHS0s29K2kYIVcgHjIyKJYss1F1bKcOLfOs6tRrtiblRtkpW9A2rTlKYInWoPjuDmOcRoQrzJCUfdtKlJchoShmZFpwU0kShMcoJH1NOYqI9LL51VN4opUhDGKhiys4V79PZCiSchL5ujNCaiR/ezPxL6+fqqDuZZQnqSIcLxYFKYMqhrMI4JAKghWbaIKwoPpWiEdIIKx0UAUdwten8H/SKZt2xWxcVUrN82UceXAEjsEpsEENNMElaIE2wECAB/AEno0749F4MV4XrTljOXMIfsB4+wRo65Mk

AAAB+HicdVDLSsNAFJ3UV62PVl26GSyiq5DEPndFNy4r2FZoQplMJ+nQySTMTIQa+yVuXCji1k9x5984fQgqeuDC4Zx7ufceP2FUKsv6MHIrq2vrG/nNwtb2zm6xtLfflXEqMOngmMXixkeSMMpJR1HFyE0iCIp8Rnr++GLm926JkDTm12qSEC9CIacBxUhpaVAq3ruShhE6cQXiISODUtkym7WqXW1AyzyrVRt1WxOn6VQdC9qmNUcZLNEelN7dYYzTiHCFGZKyb1uJ8jIkFMWMTAtuKkmC8BiFpK8pRxGRXjY/fAqPtTKEQSx0cQXn6veJDEVSTiJfd0ZIjeRvbyb+5fVTFTS8jPIkVYTjxaIgZVDFcJYCHFJBsGITTRAWVN8K8QgJhJXOqqBD+PoU/k+6jmlXzOZVpdw6X8aRB4fgCJwCG9RBC1yCNugADFLwAJ7As3FnPBovxuuiNWcsZw7ADxhvn0TDk4Y=

9 operators with a lattice site. In practice, the ratio of matrix elements (A4) can be eﬃciently computed since the update only involves a change of two operators on the same site 푗 at Head 0 position 푝1 and 푝2. | i

length=2m Tail 2. Initialization and thermalization | i We initially start with a configuration only involving 푀 ≈ 10 FIG. 9. Example of a Monte Carlo configuration in the projective identity operators, randomly positioned in real space. An algorithm: It is a snake of length 2푚 in the graph where the vertices are basis states and the edges are single spin flips. The “head” and the initial spin configuration |휎i is also generated at random. The 0 thermalization process consists of running consecutively the “tail” of the snake are |휎 i and |휎i, respectively. Its “body” consists of all intermediate basis states connecting |휎i to |휎0i by applying the identity and spin-flip updates and increasing the size of the spin flips of the sequence 푆 . sequence of operators 푀 by about 10% (by randomly adding 2푚 identity operators at the end of the sequence) when 푛 > 0.8푀, with 푛 the number of nonidentity operators. This ensures that The configurations have a convenient graphical representa- 푛 < 푀 in the following, when updates are performed in order tion: It represents a “snake” of length 2푚 on a graph where the to get measurements. vertices are basis states and the edges are single spin flips. The “head” and the “tail” of the snake are |휎0i and |휎i, respectively. Its “body” consists of all intermediate basis states connecting Appendix B: Projective quantum Monte Carlo 0 |휎i to |휎 i by applying the spin flips of the sequence 푆2푚. Including |휎i and |휎0i, there are (2푚 + 1) vertices in total, see The basic idea of projective quantum Monte Carlo [78] lies Fig.9. behind the power method,

푚 |휓gsi ∼ lim H |휙i, (B1) 푚→+∞ 2. Monte Carlo update with h휓gs|휙i ≠ 0. This algorithm was used in Ref. 59 to study the Z model at 훼 = 1/2. The update aims at moving the snake of Fig.9 on the graph: 2 It corresponds to translating the whole snake by one vertex at a time, either by its head or its tail. It can be implemented as 1. Configuration space follows, 1. Select the head or the tail of the snake at random with Based on Eq. (B1), we define the following equivalent of probability 1/2. the “partition function” (or normalization) at order 푚, 2. Independently of which one has been selected, there are 푚 푚 0 Z(푚) = 휙 H H 휙 . (B2) 푁 spins in |휎0i and |휎i which can be potentially flipped. One is selected at random with probability 1/푁. The | i | 0i Í | i Choosing for initial state 휙 = 휙 = {| 휎i} 휎 and using new basis state obtained when applying the spin flip the same notation as Eq. (3.1) for the Hamiltonian, one arrives corresponds to the update proposal, where the head or to the tail will move if it is accepted. 2푚 ∑︁ ∑︁ D ∑︁ 0E Z(푚) = 휎 푋 푗 휎 . (B3) The probability to accept such an update involves the ratio of 푗 {| 휎i} {|휎0 i} the weights of the configurations after “a” and before “b” the move, i.e., 푃accept = min(푅b→a, 1), with, Expanding the power as the product of all the possible sequences {푆 } of operators 푋 of length 2푚, one gets Î 0 2푚 푗 h휎a| 푋 ∈푆a 푋 푗 |휎ai 푗 2푚 ≥ 푅b→a = Î 0 0, (B6) ∑︁ ∑︁ ∑︁ D Ö 0E h휎b| 푋 ∈푆b 푋 푗 |휎 i Z( ) 푗 2푚 b 푚 = 휎 푋 푗 휎 . (B4) 0 {|휎i} {|휎 i} {푆2푚 } 푋푗 ∈푆2푚 which can be readily computed since only two operators 푋 푗 b a One can rewrite the partition function (B4) as differ between the two sequences 푆2푚 and 푆2푚. Because a single update is highly local in the configuration space, we ∑︁ ∑︁ ∑︁ 0 perform 2푚 of them consecutively in what we call an actual Z(푚) = 푊 |휎i, |휎 i, 푆2푚 , (B5) 0 update for this algorithm. {| 휎i} {|휎 i} {푆2푚 } At 훼 = ±1/2, some spin flips become impossible (the matrix 0 where 푊 (|휎i, |휎 i, 푆2푚) is the weight of a configuration with element is strictly zero), while the ones which remain possible 0 0 a probability 푃(|휎i, |휎 i, 푆2푚) = 푊 (|휎i, |휎 i, 푆2푚)/Z(푚) ≥ all have the same matrix element 1. In practice, one can take 0 for the parameter range 훼 ∈ [−0.5, +0.5]. advantage of this and slightly adapt the above algorithm by AAACAHicbVDLSsNAFJ3UV62vqAsXbgaLUDclKQUVXBTduKxgH9iEMJlO2qGTSZiZCCXEhb/ixoUibv0Md/6NkzYLrR64cDjnXu69x48ZlcqyvozS0vLK6lp5vbKxubW9Y+7udWWUCEw6OGKR6PtIEkY56SiqGOnHgqDQZ6TnT65yv3dPhKQRv1XTmLghGnEaUIyUljzzwAmRGvt+epd5jYeag1g8RhfWiWdWrbo1A/xL7IJUQYG2Z346wwgnIeEKMyTlwLZi5aZIKIoZySpOIkmM8ASNyEBTjkIi3XT2QAaPtTKEQSR0cQVn6s+JFIVSTkNfd+bnykUvF//zBokKztyU8jhRhOP5oiBhUEUwTwMOqSBYsakmCAuqb4V4jATCSmdW0SHYiy//Jd1G3W7Wz2+a1dZlEUcZHIIjUAM2OAUtcA3aoAMwyMATeAGvxqPxbLwZ7/PWklHM7INfMD6+Ac1xlec=

AAACB3icbVBNS8NAEN3Ur1q/oh4FCRbBQyjZUrQehKIXDx4q2FpoQtlst+3SzSbsboQScvPiX/HiQRGv/gVv/hs3bQ7a+mDg8d4MM/P8iFGpHOfbKCwtr6yuFddLG5tb2zvm7l5bhrHApIVDFoqOjyRhlJOWooqRTiQICnxG7v3xVebfPxAhacjv1CQiXoCGnA4oRkpLPfPQDZAaYcSSm/TCTWr2qV23oWPDqg3rbtozy07FmcJaJDAnZZCj2TO/3H6I44BwhRmSsgudSHkJEopiRtKSG0sSITxGQ9LVlKOASC+Z/pFax1rpW4NQ6OLKmqq/JxIUSDkJfN2ZXS3nvUz8z+vGalD3EsqjWBGOZ4sGMbNUaGWhWH0qCFZsognCgupbLTxCAmGloyvpEOD8y4ukXa3AWuX8tlZuXOZxFMEBOAInAIIz0ADXoAlaAINH8AxewZvxZLwY78bHrLVg5DP74A+Mzx8zUJb1

AAACBXicbVC7SgNBFJ2Nrxhfq5ZaDAbBYgk7IWgshKCNhUUEEwPZJcxOJsmQ2Qczs0JYtrHxV2wsFLH1H+z8G2eTLTTxwIXDOfdy7z1exJlUtv1tFJaWV1bXiuuljc2t7R1zd68tw1gQ2iIhD0XHw5JyFtCWYorTTiQo9j1O773xVebfP1AhWRjcqUlEXR8PAzZgBCst9cxDx8dqRDBPbtILJzm16hayLVS1UN1Je2bZrthTwEWCclIGOZo988vphyT2aaAIx1J2kR0pN8FCMcJpWnJiSSNMxnhIu5oG2KfSTaZfpPBYK304CIWuQMGp+nsiwb6UE9/TndnNct7LxP+8bqwGdTdhQRQrGpDZokHMoQphFgnsM0GJ4hNNMBFM3wrJCAtMlA6upENA8y8vkna1gmqV89tauXGZx1EEB+AInAAEzkADXIMmaAECHsEzeAVvxpPxYrwbH7PWgpHP7IM/MD5/AEa9loE=

AAACA3icbVDLSsNAFJ3UV62vqDvdBIvgopRMKVgXQtGNCxcV7AOaUCbTSTt0MgkzE6GEgBt/xY0LRdz6E+78GydtFtp64MLhnHu59x4vYlQq2/42Ciura+sbxc3S1vbO7p65f9CRYSwwaeOQhaLnIUkY5aStqGKkFwmCAo+Rrje5zvzuAxGShvxeTSPiBmjEqU8xUloamEdOgNQYI5bcppdO0qhAuwJrFdhw0oFZtqv2DNYygTkpgxytgfnlDEMcB4QrzJCUfWhHyk2QUBQzkpacWJII4Qkakb6mHAVEusnsh9Q61crQ8kOhiytrpv6eSFAg5TTwdGd2sVz0MvE/rx8rv+EmlEexIhzPF/kxs1RoZYFYQyoIVmyqCcKC6lstPEYCYaVjK+kQ4OLLy6RTq8J69eKuXm5e5XEUwTE4AWcAgnPQBDegBdoAg0fwDF7Bm/FkvBjvxse8tWDkM4fgD4zPH1fQlgs=

AAACAHicbVDLSsNAFJ3UV62vqAsXbgaL4KKUpFS0C6HoxoWLCvYBTSiT6aQdOpmEmYlQQjb+ihsXirj1M9z5N07aLLT1wIXDOfdy7z1exKhUlvVtFFZW19Y3ipulre2d3T1z/6Ajw1hg0sYhC0XPQ5IwyklbUcVILxIEBR4jXW9yk/ndRyIkDfmDmkbEDdCIU59ipLQ0MI+cAKkxRiy5S6+cpFGxaxX73EkHZtmqWjPAZWLnpAxytAbmlzMMcRwQrjBDUvZtK1JugoSimJG05MSSRAhP0Ij0NeUoINJNZg+k8FQrQ+iHQhdXcKb+nkhQIOU08HRndq5c9DLxP68fK//STSiPYkU4ni/yYwZVCLM04JAKghWbaoKwoPpWiMdIIKx0ZiUdgr348jLp1Kp2vdq4r5eb13kcRXAMTsAZsMEFaIJb0AJtgEEKnsEreDOejBfj3fiYtxaMfOYQ/IHx+QP5IpVe

AAACAnicbVDLSsNAFJ3UV62vqCtxM1gEF6UkpT66EIpuXLioYB/QhDKZTtuhk0mYmQglBDf+ihsXirj1K9z5N07aLLT1wIXDOfdy7z1eyKhUlvVt5JaWV1bX8uuFjc2t7R1zd68lg0hg0sQBC0THQ5IwyklTUcVIJxQE+R4jbW98nfrtByIkDfi9moTE9dGQ0wHFSGmpZx44PlIjjFh8m1w68VmpVrIrJfvUSXpm0SpbU8BFYmekCDI0euaX0w9w5BOuMENSdm0rVG6MhKKYkaTgRJKECI/RkHQ15cgn0o2nLyTwWCt9OAiELq7gVP09ESNfyonv6c70YDnvpeJ/XjdSgws3pjyMFOF4tmgQMagCmOYB+1QQrNhEE4QF1bdCPEICYaVTK+gQ7PmXF0mrUrar5dpdtVi/yuLIg0NwBE6ADc5BHdyABmgCDB7BM3gFb8aT8WK8Gx+z1pyRzeyDPzA+fwDm3ZXU

AAACBHicbVDLSsNAFJ3UV62vqMtuBovgIpSk1kcXQtGNCxcVbC00oUymk3bo5MHMRCghCzf+ihsXirj1I9z5N07aLLT1wIXDOfdy7z1uxKiQpvmtFZaWV1bXiuuljc2t7R19d68jwphj0sYhC3nXRYIwGpC2pJKRbsQJ8l1G7t3xVebfPxAuaBjcyUlEHB8NA+pRjKSS+nrZ9pEcYcSSm/TCTo6NU6NhWDXDOrHTvl4xq+YUcJFYOamAHK2+/mUPQhz7JJCYISF6lhlJJ0FcUsxIWrJjQSKEx2hIeooGyCfCSaZPpPBQKQPohVxVIOFU/T2RIF+Iie+qzuxkMe9l4n9eL5beuZPQIIolCfBskRczKEOYJQIHlBMs2UQRhDlVt0I8QhxhqXIrqRCs+ZcXSadWterVxm290rzM4yiCMjgAR8ACZ6AJrkELtAEGj+AZvII37Ul70d61j1lrQctn9sEfaJ8/0LSWRw==

(a) 3 3 AAACz3icjVHLTsJAFD3UF+ILdemmkZjghhQ00SXRjUtIBEyAmGkZoKGvtFMNIRi3/oBb/SvjH+hfeGccEpUYnabtmXPvOTP3Xjvy3ERY1mvGWFhcWl7JrubW1jc2t/LbO80kTGOHN5zQC+MrmyXccwPeEK7w+FUUc+bbHm/Zo3MZb93wOHHD4FKMI9712SBw+67DBFGdYsdnYmj3J2x6eJ0vWCVLLXMelDUoQK9amH9BBz2EcJDCB0cAQdgDQ0JPG2VYiIjrYkJcTMhVcY4pcqRNKYtTBiN2RN8B7dqaDWgvPROldugUj96YlCYOSBNSXkxYnmaqeKqcJfub90R5yruN6W9rL59YgSGxf+lmmf/VyVoE+jhVNbhUU6QYWZ2jXVLVFXlz80tVghwi4iTuUTwm7CjlrM+m0iSqdtlbpuJvKlOycu/o3BTv8pY04PLPcc6DZqVUPipV6seF6pkedRZ72EeR5nmCKi5QQ4O8IzziCc9G3bg17oz7z1QjozW7+LaMhw+kuJPv (↵ > 0) AAAC43icjVHLSsNAFD2Nr1pfVZeCBFuhbkraLnQlRTcuK9gHWi2TdLSheZFMhFLqyp07cesPuNV/Ef9A/8I7Ywo+EJ2Q5My595yZe68ZOHYkDOMlpU1MTk3PpGczc/MLi0vZ5ZVG5MehxeuW7/hhy2QRd2yP14UtHN4KQs5c0+FNs78v481LHka27x2JQcBPXXbh2ee2xQRRnex6vu0y0TPN4fHorNIp568K+TZzgh7Td438ViebM4qGWvpPUEpADsmq+dlntNGFDwsxXHB4EIQdMET0nKAEAwFxpxgSFxKyVZxjhAxpY8rilMGI7dP3gnYnCevRXnpGSm3RKQ69ISl1bJLGp7yQsDxNV/FYOUv2N++h8pR3G9DfTLxcYgV6xP6lG2f+VydrETjHjqrBppoCxcjqrMQlVl2RN9c/VSXIISBO4i7FQ8KWUo77rCtNpGqXvWUq/qoyJSv3VpIb403ekgZc+j7On6BRLpYqxfJhOVfdS0adxho2UKB5bqOKA9RQJ+9rPOARTxrXbrRb7e4jVUslmlV8Wdr9O3Oxmec= (b) AAACz3icjVHLTsJAFD3UF+ILdemmkZjghhQ00SXRjUtIBEyAmLYM0NBXplMNIRi3/oBb/SvjH+hfeGccEpUYnabtmXPvOTP3Xif2vURY1mvGWFhcWl7JrubW1jc2t/LbO80kSrnLGm7kR/zKsRPmeyFrCE/47CrmzA4cn7Wc0bmMt24YT7wovBTjmHUDexB6fc+1BVGdYiewxdDpT5zp4XW+YJUstcx5UNagAL1qUf4FHfQQwUWKAAwhBGEfNhJ62ijDQkxcFxPiOCFPxRmmyJE2pSxGGTaxI/oOaNfWbEh76ZkotUun+PRyUpo4IE1EeZywPM1U8VQ5S/Y374nylHcb09/RXgGxAkNi/9LNMv+rk7UI9HGqavCoplgxsjpXu6SqK/Lm5peqBDnExEncozgn7CrlrM+m0iSqdtlbW8XfVKZk5d7VuSne5S1pwOWf45wHzUqpfFSq1I8L1TM96iz2sI8izfMEVVyghgZ5x3jEE56NunFr3Bn3n6lGRmt28W0ZDx+nGpPw (c) Z2 (↵ < 0) AAACz3icjVHLTsJAFD3UF+ILdemmkZjghhQ00SXRjUtIBEyAmHYYoKGvtFMNIRi3/oBb/SvjH+hfeGccEpUYnabtmXPvOTP3Xify3ERY1mvGWFhcWl7JrubW1jc2t/LbO80kTGPGGyz0wvjKsRPuuQFvCFd4/CqKue07Hm85o3MZb93wOHHD4FKMI9717UHg9l1mC6I6xY5vi6HTn7Dp4XW+YJUstcx5UNagAL1qYf4FHfQQgiGFD44AgrAHGwk9bZRhISKuiwlxMSFXxTmmyJE2pSxOGTaxI/oOaNfWbEB76ZkoNaNTPHpjUpo4IE1IeTFheZqp4qlyluxv3hPlKe82pr+jvXxiBYbE/qWbZf5XJ2sR6ONU1eBSTZFiZHVMu6SqK/Lm5peqBDlExEnco3hMmCnlrM+m0iSqdtlbW8XfVKZk5Z7p3BTv8pY04PLPcc6DZqVUPipV6seF6pkedRZ72EeR5nmCKi5QQ4O8IzziCc9G3bg17oz7z1QjozW7+LaMhw+pfJPx (↵ < 0) ZAAAC43icjVHLSsNAFD2Nr1pfVZeCBFuhbkraLnThoujGZQX7QKtlko42NC+SiVBKXblzJ279Abf6L+If6F94Z0zBB6ITkpw5954zc+81A8eOhGG8pLSJyanpmfRsZm5+YXEpu7zSiPw4tHjd8h0/bJks4o7t8bqwhcNbQciZazq8afb3Zbx5ycPI9r0jMQj4qcsuPPvctpggqpNdz7ddJnqmOTwenVU65fxVId9mTtBj+q6R3+pkc0bRUEv/CUoJyCFZNT/7jDa68GEhhgsOD4KwA4aInhOUYCAg7hRD4kJCtopzjJAhbUxZnDIYsX36XtDuJGE92kvPSKktOsWhNySljk3S+JQXEpan6SoeK2fJ/uY9VJ7ybgP6m4mXS6xAj9i/dOPM/+pkLQLn2FE12FRToBhZnZW4xKor8ub6p6oEOQTESdyleEjYUspxn3WliVTtsrdMxV9VpmTl3kpyY7zJW9KAS9/H+RM0ysVSpVg+LOeqe8mo01jDBgo0z21UcYAa6uR9jQc84knj2o12q919pGqpRLOKL0u7fwdu65nl 2 Z2 0.040 0.040 0.040 = 3, 6, 9, 12, 15 = 9, 12, 15 L { } L { } = 6, 9, 12, 15 = 4, 6, 8, 10, 12, 18 0.035 = 6, 9, 12, 15 0.035 L { } 0.035 L { } L { } = 9, 12, 15 L { } = 6, 8, 10, 12, 18 0.030 0.030 0.030 L { } = 8, 10, 12, 18 0.025 0.025 0.025 L { }

0.020 0.020 0.020

0.015 0.015 0.015

0.010 0.010 0.010

0.005 0.005 0.005

0.000 0.000 0.000 0.8 1.0 1.2 0.36 0.38 0.40 0.8 1.0 1.2 0.466 0.464 0.8 1.0 1.2 0.38 0.37 2/⌫ ↵c 2/⌫ ↵c 2/⌫ ↵c

FIG. 10. (a, b, c) Minimum spreading of the crossing points defined in Eq. (C1) between all pairs of linear system sizes in L, according to Eq. (C2). A minimum is observed in all cases with a drift of its position observed as one removes the smallest system sizes from the set L. The values of 2훽/휈 and 훼c reported in TableI correspond to the minimum considering the largest system sizes only (violet curves). The error Z3 bars that we give reflect the difference with respect to the data set containing all sizes (see text). For the 2 model at 훼 < 0, we did not take into account 푁 = 3 × 3 as its crossing point with other system sizes was always leading to outliers.

only suggesting moves of the head/tail to configurations where Appendix C: Locating the transitions the spin is flippable (note that the probability 푃accept needs to be modified accordingly to satisfy detailed balance). We assume a second order phase transition for which the following finite-size critical scaling is expected at criticality 훼 = 훼c [2, 81] [see also Eq. (4.1)], 3. Initialization and thermalization 2훽/휈 퐶 풒order, 0 × 퐿 ∼ constant, (C1) The initialization and thermalization parts of the algorithms increase the length ℓ of the snake until it reaches the desired with 훽 the order parameter critical exponent, 휈 the correlation length critical exponent, and 퐿 the linear system size. value 2푚. We typically start with a snake of length ℓ = 2, generated at random on the graph and perform a number of 퐶 풒order, 0 is measured in quantum Monte Carlo for different updates of the order of the number of lattice sites (as described values of 훼 and linear system sizes 퐿. Both 2훽/휈 and 훼c are above). When this is done, we symmetrically (with respect to unknown. In order to determine them, we set 2훽/휈 as a param- the tail and the head) increase the length of the snake ℓ → ℓ+2. eter and find the value for which the crossing of the different The positions of the new head and tail are selected at random. system sizes is as close as possible to a single point, which We then repeat this whole process until ℓ ≡ 2푚. Although the gives 훼c. the position of the initial snake and the position of the new head In practice, we have a set of data corresponding to different system sizes L = {퐿 , 퐿 , . . . , 퐿 }. We compute the and tail are random, we have to ensure that the configuration 1 2 푁sizes L is valid by making sure that the corresponding operators 푋 푗 2 possible combinations between them. For a given value introduced in the sequence do not lead to zero matrix elements. of the exponent, and for each pair, we get the coordinates of Increasing the length of the snake on the fly allows one to the crossing point between the two curves (푥c, 푦c) (we do a check on whether or not its current size is sufficiently long to linear interpolation between the different data points). From probe the ground state or not (by regularly performing mea- the resulting list of coordinates {(푥c, 푦c)}, we quantify their surements, of the energy for instance), and adjust 푚 accord- spreading by computing the standard deviation of the euclidean √︁ 2 2 ingly. distances 푑c = 푥c + 푦c, √︃ 휎2훽/휈 = 푑2 − 푑 for a given L. (C2) 4. Measurements c c We estimate the best exponent 2훽/휈 from the minimum of 휎 With the projective algorithm, the measurement of an ob- and estimate 훼c as the average over the 푥c coordinates for the servable 푂 takes the form, corresponding best exponent. This method puts all the system sizes on the same level (the smallest and the largest), but we 푂 = 휙 H 푚푂H 푚 휙0 Z(푚). (B7) know that the crossings can exhibit some drifts with the system From the snake configuration perspective, if 푂 is a diagonal sizes. To that end, we repeat the procedure by removing from observable in the computational basis, it is measured on the the set L the smallest system size and the two smallest ones spin configuration positioned in the middle of the snake. If (we are limited on how far we can go by the total size of L). one wants to measure the energy hHi, it can be achieved Results are plotted in Fig. 10 for the three models considered, Í by averaging − 푗 1 − |훼| + 훼휗푗 over the head or tail spin with 휎 showing a well-defined minimum in all cases. The 0 configurations |휎 i and |휎i. values of 2훽/휈 and 훼c reported in TableI correspond to the AAAB63icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0WPRi8cK9gPaUDbbTbt0dxN2J0IJ/QtePCji1T/kzX9j0uagrQ8GHu/NMDMviKWw6LrfTmltfWNzq7xd2dnd2z+oHh61bZQYxlsskpHpBtRyKTRvoUDJu7HhVAWSd4LJXe53nrixItKPOI25r+hIi1AwirnUR5oMqjW37s5BVolXkBoUaA6qX/1hxBLFNTJJre15box+Sg0KJvms0k8sjymb0BHvZVRTxa2fzm+dkbNMGZIwMllpJHP190RKlbVTFWSdiuLYLnu5+J/XSzC88VOh4wS5ZotFYSIJRiR/nAyF4QzlNCOUGZHdStiYGsowi6eSheAtv7xK2hd177J+9XBZa9wWcZThBE7hHDy4hgbcQxNawGAMz/AKb45yXpx352PRWnKKmWP4A+fzByPFjlA= AAAB63icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0WPRi8cK9gPaUDbbTbt0dxN2J0IJ/QtePCji1T/kzX9j0uagrQ8GHu/NMDMviKWw6LrfTmltfWNzq7xd2dnd2z+oHh61bZQYxlsskpHpBtRyKTRvoUDJu7HhVAWSd4LJXe53nrixItKPOI25r+hIi1AwirnUR5oMqjW37s5BVolXkBoUaA6qX/1hxBLFNTJJre15box+Sg0KJvms0k8sjymb0BHvZVRTxa2fzm+dkbNMGZIwMllpJHP190RKlbVTFWSdiuLYLnu5+J/XSzC88VOh4wS5ZotFYSIJRiR/nAyF4QzlNCOUGZHdStiYGsowi6eSheAtv7xK2hd177J+9XBZa9wWcZThBE7hHDy4hgbcQxNawGAMz/AKb45yXpx352PRWnKKmWP4A+fzByPFjlA= AAAB63icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0WPRi8cK9gPaUDbbTbt0dxN2J0IJ/QtePCji1T/kzX9j0uagrQ8GHu/NMDMviKWw6LrfTmltfWNzq7xd2dnd2z+oHh61bZQYxlsskpHpBtRyKTRvoUDJu7HhVAWSd4LJXe53nrixItKPOI25r+hIi1AwirnUR5oMqjW37s5BVolXkBoUaA6qX/1hxBLFNTJJre15box+Sg0KJvms0k8sjymb0BHvZVRTxa2fzm+dkbNMGZIwMllpJHP190RKlbVTFWSdiuLYLnu5+J/XSzC88VOh4wS5ZotFYSIJRiR/nAyF4QzlNCOUGZHdStiYGsowi6eSheAtv7xK2hd177J+9XBZa9wWcZThBE7hHDy4hgbcQxNawGAMz/AKb45yXpx352PRWnKKmWP4A+fzByPFjlA=

AAACBnicbVDLSsNAFJ3UV62vqEsRBotQF5akLVQXQtGNywr2gU0Ik+mkHTp5MDMRSqgbN/6KGxeKuPUb3Pk3TtostPXAhcM593LvPW7EqJCG8a3llpZXVtfy64WNza3tHX13ry3CmGPSwiELeddFgjAakJakkpFuxAnyXUY67ugq9Tv3hAsaBrdyHBHbR4OAehQjqSRHP7R8JIeum9xNnMpDyUIsGqKLU6NcrZvVE0cvGmVjCrhIzIwUQYamo39Z/RDHPgkkZkiInmlE0k4QlxQzMilYsSARwiM0ID1FA+QTYSfTNybwWCl96IVcVSDhVP09kSBfiLHvqs70aDHvpeJ/Xi+W3pmd0CCKJQnwbJEXMyhDmGYC+5QTLNlYEYQ5VbdCPEQcYamSK6gQzPmXF0m7UjZr5fObWrFxmcWRBwfgCJSACeqgAa5BE7QABo/gGbyCN+1Je9HetY9Za07LZvbBH2ifP6VDl00=

AAACB3icbVDLSsNAFJ34rPUVdSlIsAh1E5I+qBuh6MZlBfvAJobJdNoOnUzCzEQooa7c+CtuXCji1l9w5984abPQ1gMXDufcy733+BElQlrWt7a0vLK6tp7byG9ube/s6nv7LRHGHOEmCmnIOz4UmBKGm5JIijsRxzDwKW77o8vUb99jLkjIbuQ4wm4AB4z0CYJSSZ5+5ARQDn0/uZ14pbvyQ9GBNBrCc8ss18rVU08vWKY1hbFI7IwUQIaGp385vRDFAWYSUShE17Yi6SaQS4IonuSdWOAIohEc4K6iDAZYuMn0j4lxopSe0Q+5KiaNqfp7IoGBEOPAV53p1WLeS8X/vG4s+2duQlgUS8zQbFE/poYMjTQUo0c4RpKOFYGIE3WrgYaQQyRVdHkVgj3/8iJplUy7YlavK4X6RRZHDhyCY1AENqiBOrgCDdAECDyCZ/AK3rQn7UV71z5mrUtaNnMA/kD7/AFwhpe7 AAACB3icbVDLSsNAFJ3UV62vqEtBgkWoC0NS62MjFN24rGAf2MQwmU7boZNJmJkIJdSVG3/FjQtF3PoL7vwbJ20W2nrgwuGce7n3Hj+iREjL+tZyc/MLi0v55cLK6tr6hr651RBhzBGuo5CGvOVDgSlhuC6JpLgVcQwDn+KmP7hM/eY95oKE7EYOI+wGsMdIlyAoleTpu04AZd/3k9uRV747eig5kEZ9eH5omZWT8oGnFy3TGsOYJXZGiiBDzdO/nE6I4gAziSgUom1bkXQTyCVBFI8KTixwBNEA9nBbUQYDLNxk/MfI2FdKx+iGXBWTxlj9PZHAQIhh4KvO9Gox7aXif147lt0zNyEsiiVmaLKoG1NDhkYaitEhHCNJh4pAxIm61UB9yCGSKrqCCsGefnmWNMqmXTGPryvF6kUWRx7sgD1QAjY4BVVwBWqgDhB4BM/gFbxpT9qL9q59TFpzWjazDf5A+/wBYsCXsg==

AAACBnicbVDLSsNAFJ3UV62vqEsRBovgqiRS0GXRjcsK9gFNCJPJpB06mYkzE6GErtz4K25cKOLWb3Dn3zhps9DWC8MczrmXe+4JU0aVdpxvq7Kyura+Ud2sbW3v7O7Z+wddJTKJSQcLJmQ/RIowyklHU81IP5UEJSEjvXB8Xei9ByIVFfxOT1LiJ2jIaUwx0oYK7GMvFCxSk8R8+f008BKkR2GcCxkROQ3sutNwZgWXgVuCOiirHdhfXiRwlhCuMUNKDVwn1X6OpKaYkWnNyxRJER6jIRkYyFFClJ/PzpjCU8NEMBbSPK7hjP09kaNEFU5NZ+FSLWoF+Z82yHR86eeUp5kmHM8XxRmDWsAiExhRSbBmEwMQltR4hXiEJMLaJFczIbiLJy+D7nnDbTaat81666qMowqOwAk4Ay64AC1wA9qgAzB4BM/gFbxZT9aL9W59zFsrVjlzCP6U9fkD73KaEQ==

3 3 Z2 (↵ =0.3735) Z2 (↵ = 0.462) Z2 (↵ = 0.3713) q q 0.1 0.1 order 3 3 0.1 order 3 3 4 4 10 10 ⇥ (b) ⇥ ⇥ ⇥ (a) AAACz3icjVHLTsJAFD3UF+ILdemmkZjghhQ00SXRjUtIBEyAmLYM0NBXplMNIRi3/oBb/SvjH+hfeGccEpUYnabtmXPvOTP3Xif2vURY1mvGWFhcWl7JrubW1jc2t/LbO80kSrnLGm7kR/zKsRPmeyFrCE/47CrmzA4cn7Wc0bmMt24YT7wovBTjmHUDexB6fc+1BVGdYiewxdDpT5zp4XW+YJUstcx5UNagAL1qUf4FHfQQwUWKAAwhBGEfNhJ62ijDQkxcFxPiOCFPxRmmyJE2pSxGGTaxI/oOaNfWbEh76ZkotUun+PRyUpo4IE1EeZywPM1U8VQ5S/Y374nylHcb09/RXgGxAkNi/9LNMv+rk7UI9HGqavCoplgxsjpXu6SqK/Lm5peqBDnExEncozgn7CrlrM+m0iSqdtlbW8XfVKZk5d7VuSne5S1pwOWf45wHzUqpfFSq1I8L1TM96iz2sI8izfMEVVyghgZ5x3jEE56NunFr3Bn3n6lGRmt28W0ZDx+nGpPw p AAACz3icjVHLTsJAFD3UF+ILdemmkZjghhQ00SXRjUtIBEyAmGkZoKGvtFMNIRi3/oBb/SvjH+hfeGccEpUYnabtmXPvOTP3Xjvy3ERY1mvGWFhcWl7JrubW1jc2t/LbO80kTGOHN5zQC+MrmyXccwPeEK7w+FUUc+bbHm/Zo3MZb93wOHHD4FKMI9712SBw+67DBFGdYsdnYmj3J2x6eJ0vWCVLLXMelDUoQK9amH9BBz2EcJDCB0cAQdgDQ0JPG2VYiIjrYkJcTMhVcY4pcqRNKYtTBiN2RN8B7dqaDWgvPROldugUj96YlCYOSBNSXkxYnmaqeKqcJfub90R5yruN6W9rL59YgSGxf+lmmf/VyVoE+jhVNbhUU6QYWZ2jXVLVFXlz80tVghwi4iTuUTwm7CjlrM+m0iSqdtlbpuJvKlOycu/o3BTv8pY04PLPcc6DZqVUPipV6seF6pkedRZ72EeR5nmCKi5QQ4O8IzziCc9G3bg17oz7z1QjozW7+LaMhw+kuJPv (c) 2⇡/3, 2⇡/3 3 AAACz3icjVHLTsJAFD3UF+ILdemmkZjghhQ00SXRjUtIBEyAmHYYoKGvtFMNIRi3/oBb/SvjH+hfeGccEpUYnabtmXPvOTP3Xify3ERY1mvGWFhcWl7JrubW1jc2t/LbO80kTGPGGyz0wvjKsRPuuQFvCFd4/CqKue07Hm85o3MZb93wOHHD4FKMI9717UHg9l1mC6I6xY5vi6HTn7Dp4XW+YJUstcx5UNagAL1qYf4FHfQQgiGFD44AgrAHGwk9bZRhISKuiwlxMSFXxTmmyJE2pSxOGTaxI/oOaNfWbEB76ZkoNaNTPHpjUpo4IE1IeTFheZqp4qlyluxv3hPlKe82pr+jvXxiBYbE/qWbZf5XJ2sR6ONU1eBSTZFiZHVMu6SqK/Lm5peqBDlExEnco3hMmCnlrM+m0iSqdtlbW8XfVKZk5Z7p3BTv8pY04PLPcc6DZqVUPipV6seF6pkedRZ72EeR5nmCKi5QQ4O8IzziCc9G3bg17oz7z1QjozW7+LaMhw+pfJPx 0, 0 6 6 6 6 6 6 12 12 ⇥ ± ± ⇥ ⇥ ⇥ ⌧ 9 9 p 9 9 8 8 18 18 , 4⇡/3, 0 ⇥ 2⇡/3, 2⇡/3 3 ⇥ ⇥ ⇥ 12 12 ± ⌥ 12 12 2⇡/3, 2⇡/p3 ⇥ 0, 4⇡/3p3 ⇥ 15 15 ± 15 15 ⇥ ⇥ qorder 0.01 order 0.01 0, 2⇡/p3

q 0.01

⇡, ⇡/p3 ± AAACGHicbVDLSsNAFJ3UV62vqEs3g0WoIDWRii6L3bisYB/QhDCZTNqhk0ycmQgl9DPc+CtuXCjitjv/xkmbhVYPDHM491zuvcdPGJXKsr6M0srq2vpGebOytb2zu2fuH3QlTwUmHcwZF30fScJoTDqKKkb6iSAo8hnp+eNWXu89EiEpj+/VJCFuhIYxDSlGSkueed5yfDpkNcfnLJCTSH/Zw9RzIqRGfphxERAxPXMUSnOfOPXMqlW35oB/iV2QKijQ9syZE3CcRiRWmCEpB7aVKDdDQlHMyLTipJIkCI/RkAw0jVFEpJvND5vCE60EMORCv1jBufqzI0ORzJfWznxhuVzLxf9qg1SF125G4yRVJMaLQWHKoOIwTwkGVBCs2EQThAXVu0I8QgJhpbOs6BDs5ZP/ku5F3W7UL+8a1eZNEUcZHIFjUAM2uAJNcAvaoAMweAIv4A28G8/Gq/FhfC6sJaPoOQS/YMy+ATr4oSY= C 0.001 0.001

0 10 20 30 0 10 20 30 0 10 20 30 40 ⌧ ⌧ ⌧

Z3 Z ' FIG. 11. (a, b, c) Fourier transform of the imaginary-time two-point correlation function (3.4) for the 2 and 2 models at 훼 훼c, computed at 풒order, see TableI. An exponential fit of the data versus imaginary time gives access to the gap Δ between the two Z2 Ising symmetry sectors for a given finite-size system. The value of the gap is reported in Fig.5 with a finite-size critical scaling of the form (4.2).

Z3 minimum considering the largest system sizes only (violet Appendix E: Binder ratio with individual sublattices for the 2 curves). The error bars that we give reflect the difference with model with 훼 > 0 respect to the data set containing all sizes. In that sense, this is more of an upper bound since we see that the difference 3 For the Z2 model with 훼 > 0, in order to check whether or between the position of the minima decreases when removing not ferromagnetic order settles in one, two, or the three sub- the smallest system sizes. lattices A, B, and C (whereas the others remain paramagnetic), see Fig.2 (c), we compute the following Binder ratio,

푚 , 푚4  A 퐵(푚) = with 푚 ≡  (E1) 2 푚A푚B, 푚2  푚A푚B푚C,  which involves the magnetization of each sublattice independently,

Appendix D: Extracting the gap Δ between the two Ising 3 ∑︁ 푧 푚A = 휎 . (E2) symmetry sectors 푁 푖∈A 푖 It is plotted in Fig. 12 (a, b, c) versus 훼, with a crossing of the largest system sizes observed in all cases, meaning that the The gap Δ between the two Z2 Ising symmetry sectors, re- three sublattices experience long-range ordering. For 퐵(푚A) ported in Fig.5 for different system sizes 퐿 × 퐿, is indirectly and 퐵(푚A푚B), the crossing point seems to happen at a larger accessed in quantum Monte Carlo by performing an expo- 훼 value than for the structure factor 퐶(풒order, 0) of the main nential fit ∝ exp(−Δ휏) of the imaginary-time correlation data text. However, this is attributed to the “effectively” smaller displayed in Fig. 11. The imaginary time 휏 is defined within system size when considering the sublattices independently the range [0, 훽/2] with inverse temperature 훽 = 4퐿 consid- (each account for 푁/3 lattice sites only), as compared to the ered. The fit to extract the gap is only performed over the range other Binder ratio 퐵(푚A푚B푚C) or the structure factor since 2 ≡ ( ) showing a genuine exponential decay. we expect that 푚A 퐶 0, 0 .

[1] Lev Davidovich Landau, “On the theory of phase transitions,” [5] Xie Chen, Zheng-Xin Liu, and Xiao-Gang Wen, “Two- Zh. Éksp. Teor. Fiz. 7, 19 (1937). dimensional symmetry-protected topological orders and their [2] S. Sachdev, Quantum Phase Transitions (Cambridge Univer- protected gapless edge excitations,” Phys. Rev. B 84, 235141 sity Press, 2001). (2011). [3] Xie Chen, Zheng-Cheng Gu, Zheng-Xin Liu, and Xiao-Gang [6] Yuan-Ming Lu and Ashvin Vishwanath, “Theory and classiﬁ- Wen, “Symmetry-protected topological orders in interacting cation of interacting integer topological phases in two dimen- bosonic systems,” Science 338, 1604–1606 (2012). sions: A chern-simons approach,” Phys. Rev. B 86, 125119 [4] Xie Chen, Zheng-Cheng Gu, Zheng-Xin Liu, and Xiao- (2012). Gang Wen, “Symmetry protected topological orders and the [7] Zhen Bi, Alex Rasmussen, Kevin Slagle, and Cenke Xu, “Clas- group cohomology of their symmetry group,” Phys. Rev. B 87, siﬁcation and description of bosonic symmetry protected topo- 155114 (2013). logical phases with semiclassical nonlinear sigma models,” AAAB7XicbVBNS8NAEJ3Ur1q/qh69LBbBU0mkoN6KXjxWsB/QhjLZbtq1m03Y3Qgl9D948aCIV/+PN/+N2zYHbX0w8Hhvhpl5QSK4Nq777RTW1jc2t4rbpZ3dvf2D8uFRS8epoqxJYxGrToCaCS5Z03AjWCdRDKNAsHYwvp357SemNI/lg5kkzI9wKHnIKRortXookhH2yxW36s5BVomXkwrkaPTLX71BTNOISUMFat313MT4GSrDqWDTUi/VLEE6xiHrWioxYtrP5tdOyZlVBiSMlS1pyFz9PZFhpPUkCmxnhGakl72Z+J/XTU145WdcJqlhki4WhakgJiaz18mAK0aNmFiCVHF7K6EjVEiNDahkQ/CWX14lrYuqV6te39cq9Zs8jiKcwCmcgweXUIc7aEATKDzCM7zCmxM7L86787FoLTj5zDH8gfP5A49hjyQ= AAAB7XicbVBNS8NAEJ3Ur1q/qh69LBbBU0mkoN6KXjxWsB/QhjLZbtq1m03Y3Qgl9D948aCIV/+PN/+N2zYHbX0w8Hhvhpl5QSK4Nq777RTW1jc2t4rbpZ3dvf2D8uFRS8epoqxJYxGrToCaCS5Z03AjWCdRDKNAsHYwvp357SemNI/lg5kkzI9wKHnIKRortXookhH2yxW36s5BVomXkwrkaPTLX71BTNOISUMFat313MT4GSrDqWDTUi/VLEE6xiHrWioxYtrP5tdOyZlVBiSMlS1pyFz9PZFhpPUkCmxnhGakl72Z+J/XTU145WdcJqlhki4WhakgJiaz18mAK0aNmFiCVHF7K6EjVEiNDahkQ/CWX14lrYuqV6te39cq9Zs8jiKcwCmcgweXUIc7aEATKDzCM7zCmxM7L86787FoLTj5zDH8gfP5A49hjyQ= AAAB7XicbVBNS8NAEJ3Ur1q/qh69LBbBU0mkoN6KXjxWsB/QhjLZbtq1m03Y3Qgl9D948aCIV/+PN/+N2zYHbX0w8Hhvhpl5QSK4Nq777RTW1jc2t4rbpZ3dvf2D8uFRS8epoqxJYxGrToCaCS5Z03AjWCdRDKNAsHYwvp357SemNI/lg5kkzI9wKHnIKRortXookhH2yxW36s5BVomXkwrkaPTLX71BTNOISUMFat313MT4GSrDqWDTUi/VLEE6xiHrWioxYtrP5tdOyZlVBiSMlS1pyFz9PZFhpPUkCmxnhGakl72Z+J/XTU145WdcJqlhki4WhakgJiaz18mAK0aNmFiCVHF7K6EjVEiNDahkQ/CWX14lrYuqV6te39cq9Zs8jiKcwCmcgweXUIc7aEATKDzCM7zCmxM7L86787FoLTj5zDH8gfP5A49hjyQ=

3 3 (a) ⇥ 7 3 3 3 3 AAACz3icjVHLTsJAFD3UF+ILdemmkZjghhQ00SXRjUtIBEyAmGkZoKGvtFMNIRi3/oBb/SvjH+hfeGccEpUYnabtmXPvOTP3Xjvy3ERY1mvGWFhcWl7JrubW1jc2t/LbO80kTGOHN5zQC+MrmyXccwPeEK7w+FUUc+bbHm/Zo3MZb93wOHHD4FKMI9712SBw+67DBFGdYsdnYmj3J2x6eJ0vWCVLLXMelDUoQK9amH9BBz2EcJDCB0cAQdgDQ0JPG2VYiIjrYkJcTMhVcY4pcqRNKYtTBiN2RN8B7dqaDWgvPROldugUj96YlCYOSBNSXkxYnmaqeKqcJfub90R5yruN6W9rL59YgSGxf+lmmf/VyVoE+jhVNbhUU6QYWZ2jXVLVFXlz80tVghwi4iTuUTwm7CjlrM+m0iSqdtlbpuJvKlOycu/o3BTv8pY04PLPcc6DZqVUPipV6seF6pkedRZ72EeR5nmCKi5QQ4O8IzziCc9G3bg17oz7z1QjozW7+LaMhw+kuJPv ⇥ ⇥

6 6 (b) (c) AAACz3icjVHLTsJAFD3UF+ILdemmkZjghhQ00SXRjUtIBEyAmLYM0NBXplMNIRi3/oBb/SvjH+hfeGccEpUYnabtmXPvOTP3Xif2vURY1mvGWFhcWl7JrubW1jc2t/LbO80kSrnLGm7kR/zKsRPmeyFrCE/47CrmzA4cn7Wc0bmMt24YT7wovBTjmHUDexB6fc+1BVGdYiewxdDpT5zp4XW+YJUstcx5UNagAL1qUf4FHfQQwUWKAAwhBGEfNhJ62ijDQkxcFxPiOCFPxRmmyJE2pSxGGTaxI/oOaNfWbEh76ZkotUun+PRyUpo4IE1EeZywPM1U8VQ5S/Y374nylHcb09/RXgGxAkNi/9LNMv+rk7UI9HGqavCoplgxsjpXu6SqK/Lm5peqBDnExEncozgn7CrlrM+m0iSqdtlbW8XfVKZk5d7VuSne5S1pwOWf45wHzUqpfFSq1I8L1TM96iz2sI8izfMEVVyghgZ5x3jEE56NunFr3Bn3n6lGRmt28W0ZDx+nGpPw 6 6 AAACz3icjVHLTsJAFD3UF+ILdemmkZjghhQ00SXRjUtIBEyAmHYYoKGvtFMNIRi3/oBb/SvjH+hfeGccEpUYnabtmXPvOTP3Xify3ERY1mvGWFhcWl7JrubW1jc2t/LbO80kTGPGGyz0wvjKsRPuuQFvCFd4/CqKue07Hm85o3MZb93wOHHD4FKMI9717UHg9l1mC6I6xY5vi6HTn7Dp4XW+YJUstcx5UNagAL1qYf4FHfQQgiGFD44AgrAHGwk9bZRhISKuiwlxMSFXxTmmyJE2pSxOGTaxI/oOaNfWbEB76ZkoNaNTPHpjUpo4IE1IeTFheZqp4qlyluxv3hPlKe82pr+jvXxiBYbE/qWbZf5XJ2sR6ONU1eBSTZFiZHVMu6SqK/Lm5peqBDlExEnco3hMmCnlrM+m0iSqdtlbW8XfVKZk5Z7p3BTv8pY04PLPcc6DZqVUPipV6seF6pkedRZ72EeR5nmCKi5QQ4O8IzziCc9G3bg17oz7z1QjozW7+LaMhw+pfJPx 6 6 ⇥ ⇥ ⇥ C 9 9 9 9 15 9 9

2.5 ⇥ B 6 ⇥ ⇥ 12 12 m 12 12

A 12 12 ⇥ ⇥ ⇥ B 15 15 m 15 15 15 15 m

⇥ A ⇥ ⇥ m 5 A m AAACAnicbVDLSsNAFJ3UV62vqCtxM1iEuimJFNRdrRuXFewDmhAm00k7dCYJMxOhhOLGX3HjQhG3foU7/8ZJmoW2HrhwOOde7r3HjxmVyrK+jdLK6tr6RnmzsrW9s7tn7h90ZZQITDo4YpHo+0gSRkPSUVQx0o8FQdxnpOdPbjK/90CEpFF4r6YxcTkahTSgGCkteeZRy/HpiNW453CkxjJIr2eZIs48s2rVrRxwmdgFqYICbc/8coYRTjgJFWZIyoFtxcpNkVAUMzKrOIkkMcITNCIDTUPEiXTT/IUZPNXKEAaR0BUqmKu/J1LEpZxyX3fmdy56mfifN0hUcOmmNIwTRUI8XxQkDKoIZnnAIRUEKzbVBGFB9a0Qj5FAWOnUKjoEe/HlZdI9r9uN+tVdo9psFXGUwTE4ATVggwvQBLegDToAg0fwDF7Bm/FkvBjvxse8tWQUM4fgD4zPH+UnlyI=

B 10 m AAACDnicbZDLSsNAFIYnXmu9RV26GSyFuimJFNRdrRuXFewFmhAm00k7dDIJMxOhhDyBG1/FjQtF3Lp259s4TQNq64GBj/8/hznn92NGpbKsL2NldW19Y7O0Vd7e2d3bNw8OuzJKBCYdHLFI9H0kCaOcdBRVjPRjQVDoM9LzJ9czv3dPhKQRv1PTmLghGnEaUIyUljyz2nJ8OmK10HNCpMYySK+yH25lM1ecembFqlt5wWWwC6iAotqe+ekMI5yEhCvMkJQD24qVmyKhKGYkKzuJJDHCEzQiA40chUS6aX5OBqtaGcIgEvpxBXP190SKQimnoa878z0XvZn4nzdIVHDhppTHiSIczz8KEgZVBGfZwCEVBCs21YCwoHpXiMdIIKx0gmUdgr148jJ0z+p2o35526g0W0UcJXAMTkAN2OAcNMENaIMOwOABPIEX8Go8Gs/Gm/E+b10xipkj8KeMj2+fHZx/ B 2.0 4 AAACGnicbZDLSsNAFIYn9VbrLerSzWAR6qYkUlB3td24rGAv0IQwmU7aoZMLMxOhhDyHG1/FjQtF3Ikb38ZJGrC2/jDw8Z9zmHN+N2JUSMP41kpr6xubW+Xtys7u3v6BfnjUE2HMMenikIV84CJBGA1IV1LJyCDiBPkuI3132s7q/QfCBQ2DezmLiO2jcUA9ipFUlqObLculY1bzHctHciK85Cb95dYCt9Osk587etWoG7ngKpgFVEGhjqN/WqMQxz4JJGZIiKFpRNJOEJcUM5JWrFiQCOEpGpOhwgD5RNhJfloKz5Qzgl7I1QskzN3FiQT5Qsx8V3Xmey7XMvO/2jCW3pWd0CCKJQnw/CMvZlCGMMsJjignWLKZAoQ5VbtCPEEcYanSrKgQzOWTV6F3UTcb9eu7RrXZKuIogxNwCmrABJegCW5BB3QBBo/gGbyCN+1Je9HetY95a0krZo7BH2lfP5Z6od0= B 3 5 0.30 0.35 0.40 0.45 0.30 0.35 0.40 0.45 0.30 0.35 0.40 0.45 ↵ ↵ ↵

Z3 ∈ [ ] FIG. 12. (a, b, c) Binder ratio 퐵 of Eq. (E1) for the 2 model with 훼 < 0 involving the magnetization 푚X of each sublattice X A, B, C as deﬁned in Eq. (E2).

Phys. Rev. B 91, 134404 (2015). [23] T. Senthil and Michael Levin, “Integer quantum hall effect for [8] Frank Pollmann, Ari M. Turner, Erez Berg, and Masaki Os- bosons,” Phys. Rev. Lett. 110, 046801 (2013). hikawa, “Entanglement spectrum of a topological phase in one [24] Shunsuke Furukawa and Masahito Ueda, “Integer quantum dimension,” Phys. Rev. B 81, 064439 (2010). hall state in two-component bose gases in a synthetic magnetic [9] Frank Pollmann, Erez Berg, Ari M. Turner, and Masaki Os- field,” Phys. Rev. Lett. 111, 090401 (2013). hikawa, “Symmetry protection of topological phases in one- [25] Ying-Hai Wu and Jainendra K. Jain, “Quantum hall effect of dimensional quantum spin systems,” Phys. Rev. B 85, 075125 two-component bosons at fractional and integral fillings,” Phys. (2012). Rev. B 87, 245123 (2013). [10] Yi-Zhuang You, Zhen Bi, Alex Rasmussen, Kevin Slagle, and [26] N. Regnault and T. Senthil, “Microscopic model for the boson Cenke Xu, “Wave function and strange correlator of short- integer quantum hall effect,” Phys. Rev. B 88, 161106 (2013). range entangled states,” Phys. Rev. Lett. 112, 247202 (2014). [27] Tarun Grover and Ashvin Vishwanath, “Quantum phase tran- [11] Zohar Ringel and Steven H. Simon, “Hidden order and flux sition between integer quantum hall states of bosons,” Phys. attachment in symmetry-protected topological phases: A Rev. B 87, 045129 (2013). laughlin-like approach,” Phys. Rev. B 91, 195117 (2015). [28] Scott D. Geraedts and Olexei I. Motrunich, “Exact realization [12] Thomas Scaffidi and Zohar Ringel, “Wave functions of of integer and fractional quantum hall phases in u(1)×u(1) symmetry-protected topological phases from conformal field models in (2+1)d,” Ann. Phys. (N. Y.) 334, 288 – 315 (2013). theories,” Phys. Rev. B 93, 115105 (2016). [29] Yin-Chen He, Subhro Bhattacharjee, R. Moessner, and Frank [13] Juven C. Wang, Luiz H. Santos, and Xiao-Gang Wen, “Bosonic Pollmann, “Bosonic integer quantum hall effect in an interact- anomalies, induced fractional quantum numbers, and degen- ing lattice model,” Phys. Rev. Lett. 115, 116803 (2015). erate zero modes: The anomalous edge physics of symmetry- [30] Scott Geraedts and Olexei I. Motrunich, “Lattice realization protected topological states,” Phys. Rev. B 91, 195134 (2015). of a bosonic integer quantum hall state–trivial insulator transi- [14] Luiz H. Santos and Juven Wang, “Symmetry-protected many- tion and relation to the self-dual line in the easy-plane nccp1 body aharonov-bohm effect,” Phys. Rev. B 89, 195122 (2014). model,” Phys. Rev. B 96, 115137 (2017). [15] Juven C. Wang, Zheng-Cheng Gu, and Xiao-Gang Wen, [31] J. P. Kestner, Bin Wang, Jay D. Sau, and S. Das Sarma, “Field-theory representation of gauge-gravity symmetry- “Prediction of a gapless topological haldane liquid phase in protected topological invariants, group cohomology, and be- a one-dimensional cold polar molecular lattice,” Phys. Rev. B yond,” Phys. Rev. Lett. 114, 031601 (2015). 83, 174409 (2011). [16] F. D. M. Haldane, “Nonlinear field theory of large-spin heisen- [32] Tarun Grover and Ashvin Vishwanath, “Quantum critical- berg antiferromagnets: Semiclassically quantized solitons of ity in topological insulators and superconductors: Emer- the one-dimensional easy-axis néel state,” Phys. Rev. Lett. 50, gence of strongly coupled majoranas and supersymmetry,” 1153–1156 (1983). arXiv:1206.1332 (2012). [17] F.D.M. Haldane, “Continuum dynamics of the 1-d heisenberg [33] Anna Keselman and Erez Berg, “Gapless symmetry-protected antiferromagnet: Identification with the o(3) nonlinear sigma topological phase of fermions in one dimension,” Phys. Rev. B model,” Phys. Lett. A 93, 464 – 468 (1983). 91, 235309 (2015). [18] Ian Affleck and F. D. M. Haldane, “Critical theory of quantum [34] Thomas Scaffidi, Daniel E. Parker, and Romain Vasseur, “Gap- spin chains,” Phys. Rev. B 36, 5291–5300 (1987). less symmetry-protected topological order,” Phys. Rev. X 7, [19] Ian Affleck, Tom Kennedy, Elliott H. Lieb, and Hal Tasaki, 041048 (2017). “Rigorous results on valence-bond ground states in antiferro- [35] Daniel E. Parker, Thomas Scaffidi, and Romain Vasseur, magnets,” Phys. Rev. Lett. 59, 799–802 (1987). “Topological luttinger liquids from decorated domain walls,” [20] Ian Affleck, Tom Kennedy, Elliott H. Lieb, and Hal Tasaki, Phys. Rev. B 97, 165114 (2018). “Valence bond ground states in isotropic quantum antiferro- [36] Daniel E. Parker, Romain Vasseur, and Thomas Scaf- magnets,” Comm. Math. Phys. 115, 477–528 (1988). fidi, “Topologically protected long edge coherence times [21] I Affleck, “Quantum spin chains and the haldane gap,” J. Phys. in symmetry-broken phases,” Phys. Rev. Lett. 122, 240605 Condens. Matter 1, 3047–3072 (1989). (2019). [22] Zheng-Cheng Gu and Xiao-Gang Wen, “Tensor-entanglement- [37] Ruben Verresen, Roderich Moessner, and Frank Pollmann, filtering renormalization approach and symmetry-protected “One-dimensional symmetry protected topological phases and topological order,” Phys. Rev. B 80, 155131 (2009). their transitions,” Phys. Rev. B 96, 165124 (2017). 13

[38] Ruben Verresen, Nick G. Jones, and Frank Pollmann, “Topol- [56] Ching-Yu Huang and Tzu-Chieh Wei, “Detecting and ogy and edge modes in quantum critical chains,” Phys. Rev. identifying two-dimensional symmetry-protected topological, Lett. 120, 057001 (2018). symmetry-breaking, and intrinsic topological phases with [39] Ruben Verresen, Ryan Thorngren, Nick G. Jones, and modular matrices via tensor-network methods,” Phys. Rev. B Frank Pollmann, “Gapless topological phases and symmetry- 93, 155163 (2016). enriched quantum criticality,” arXiv:1905.06969 (2019). [57] Mohsin Iqbal, Kasper Duivenvoorden, and Norbert Schuch, [40] Carlos M. Duque, Hong-Ye Hu, Yi-Zhuang You, Vedika Khe- “Study of anyon condensation and topological phase transitions mani, Ruben Verresen, and Romain Vasseur, “Topological from a z4 topological phase using the projected entangled pair and symmetry-enriched random quantum critical points,” Phys. states approach,” Phys. Rev. B 97, 195124 (2018). Rev. B 103, L100207 (2021). [58] Wen-Tao Xu and Guang-Ming Zhang, “Tensor network state [41] Yuan-Ming Lu and Dung-Hai Lee, “Quantum phase transi- approach to quantum topological phase transitions and their tions between bosonic symmetry-protected topological phases criticalities of z2 topologically ordered states,” Phys. Rev. B in two dimensions: Emergent qed3 and anyon superfluid,” 98, 165115 (2018). Phys. Rev. B 89, 195143 (2014). [59] Maxime Dupont, Snir Gazit, and Thomas Scaffidi, “Evidence [42] Siddhardh C. Morampudi, Curt von Keyserlingk, and Frank for deconfined 푢(1) gauge theory at the transition between toric Pollmann, “Numerical study of a transition between z2 topo- code and double semion,” Phys. Rev. B 103, L140412 (2021). logically ordered phases,” Phys. Rev. B 90, 035117 (2014). [60] Steven R. White, “Density matrix formulation for quantum [43] Lokman Tsui, Fa Wang, and Dung-Hai Lee, “Topological ver- renormalization groups,” Phys. Rev. Lett. 69, 2863–2866 sus landau-like phase transitions,” arXiv:1511.07460 (2015). (1992). [44] Lokman Tsui, Hong-Chen Jiang, Yuan-Ming Lu, and Dung- [61] Steven R. White, “Density-matrix algorithms for quantum Hai Lee, “Quantum phase transitions between a class of sym- renormalization groups,” Phys. Rev. B 48, 10345–10356 metry protected topological states,” Nucl. Phys. B 896, 330 – (1993). 359 (2015). [62] Matthias Troyer and Uwe-Jens Wiese, “Computational com- [45] Yi-Zhuang You, Zhen Bi, Dan Mao, and Cenke Xu, “Quantum plexity and fundamental limitations to fermionic quantum phase transitions between bosonic symmetry-protected topo- monte carlo simulations,” Phys. Rev. Lett. 94, 170201 (2005). logical states without sign problem: Nonlinear sigma model [63] Michael Levin and Zheng-Cheng Gu, “Braiding statistics ap- with a topological term,” Phys. Rev. B 93, 125101 (2016). proach to symmetry-protected topological phases,” Phys. Rev. [46] Yuan-Yao He, Han-Qing Wu, Yi-Zhuang You, Cenke Xu, B 86, 115109 (2012). Zi Yang Meng, and Zhong-Yi Lu, “Bona fide interaction- [64] Beni Yoshida, “Gapped boundaries, group cohomology and driven topological phase transition in correlated symmetry- fault-tolerant logical gates,” Ann. Phys. (N. Y.) 377, 387 – 413 protected topological states,” Phys. Rev. B 93, 115150 (2016). (2017). [47] Yi-Zhuang You, Yin-Chen He, Ashvin Vishwanath, and Cenke [65] Beni Yoshida, “Topological phases with generalized global Xu, “From bosonic topological transition to symmetric fermion symmetries,” Phys. Rev. B 93, 155131 (2016). mass generation,” Phys. Rev. B 97, 125112 (2018). [66] A system in 퐷 spatial dimensions can have nontrivial “weak” [48] Lokman Tsui, Yen-Ta Huang, Hong-Chen Jiang, and Dung- 퐷0-dimensional SPT order, for 퐷0 < 퐷 [4]. The protection of Hai Lee, “The phase transitions between zn×zn bosonic topo- this weak SPT order requires translation symmetry to be re- logical phases in 1+1d, and a constraint on the central charge spected, on top of the protecting on-site symmetry. For 퐷0 = 0 for the critical points between bosonic symmetry protected and a Z2 symmetry, the classification is Z2. For a symme- topological phases,” Nucl. Phys. B 919, 470 – 503 (2017). try generated by 휎푥, parent Hamiltonians of the two zero- [49] Zhen Bi and T. Senthil, “Adventure in topological phase transi- dimensional SPT classes are given by −휎푥 and 휎푥, and the tions in 3 + 1-d: Non-abelian deconfined quantum criticalities twisting unitary which relates the two is given by 휎푧. and a possible duality,” Phys. Rev. X 9, 021034 (2019). [67] M. de Wild Propitius, Topological interactions in broken gauge [50] Nick Bultinck, “Uv perspective on mixed anomalies at criti- theories, Ph.D. thesis, PhD Thesis, 1995 (1995). cal points between bosonic symmetry-protected phases,” Phys. [68] Zheng-Cheng Gu, Juven C. Wang, and Xiao-Gang Wen, “Mul- Rev. B 100, 165132 (2019). tikink topological terms and charge-binding domain-wall con- [51] Samuel Gozel, Didier Poilblanc, Ian Affleck, and Frédéric densation induced symmetry-protected topological states: Be- Mila, “Novel families of su(n) aklt states with arbitrary self- yond chern-simons/bf field theories,” Phys. Rev. B 93, 115136 conjugate edge states,” Nucl. Phys. B 945, 114663 (2019). (2016). [52] Tian-Sheng Zeng, D. N. Sheng, and W. Zhu, “Continuous [69] Huan He, Yunqin Zheng, and Curt von Keyserlingk, “Field phase transition between bosonic integer quantum hall liquid theories for gauged symmetry-protected topological phases: Ë and a trivial insulator: Evidence for deconfined quantum crit- z 3 Non-abelian anyons with abelian gauge group 2 ,” Phys. icality,” Phys. Rev. B 101, 035138 (2020). Rev. B 95, 035131 (2017). [53] T. Senthil, Ashvin Vishwanath, Leon Balents, Subir Sachdev, [70] Pavel Putrov, Juven Wang, and Shing-Tung Yau, “Braiding and Matthew P. A. Fisher, “Deconfined quantum critical statistics and link invariants of bosonic/fermionic topological points,” Science 303, 1490–1494 (2004). quantum matter in 2+1 and 3+1 dimensions,” Annals of Physics [54] Chong Wang, Adam Nahum, Max A. Metlitski, Cenke Xu, and 384, 254 – 287 (2017). T. Senthil, “Deconfined quantum critical points: Symmetries [71] Anders W. Sandvik, “Computational studies of quantum spin and dualities,” Phys. Rev. X 7, 031051 (2017). systems,” AIP Conf. Proc. 1297, 135–338 (2010). [55] Yan Qi Qin, Yuan-Yao He, Yi-Zhuang You, Zhong-Yi Lu, [72] Anders W. Sandvik, Many-Body Methods for Real Materi- Arnab Sen, Anders W. Sandvik, Cenke Xu, and Zi Yang als, Modeling and Simulation, edited by Eva Pavarini, Erik Meng, “Duality between the deconfined quantum-critical point Koch, and Shiwei Zhang, Verlag des Forschungszentrum and the bosonic topological transition,” Phys. Rev. X 7, 031052 Julich, Vol. 9 (Modeling and Simulation, 2019). (2017). 14

[73] Anders W. Sandvik, “Stochastic series expansion method [91] Andrew Smerald, Sergey Korshunov, and Frédéric Mila, with operator-loop update,” Phys. Rev. B 59, R14157–R14160 “Topological aspects of symmetry breaking in triangular- (1999). lattice ising antiferromagnets,” Phys. Rev. Lett. 116, 197201 [74] Olav F. Syljuåsen and Anders W. Sandvik, “Quantum monte (2016). carlo with directed loops,” Phys. Rev. E 66, 046701 (2002). [92] L. S. Levitov, “Equivalence of the dimer resonating-valence- [75] Anders W. Sandvik, “Stochastic series expansion method for bond problem to the quantum roughening problem,” Phys. Rev. quantum ising models with arbitrary interactions,” Phys. Rev. Lett. 64, 92–94 (1990). E 68, 056701 (2003). [93] Eduardo Fradkin, David A. Huse, R. Moessner, V. Oganesyan, [76] H. G. Evertz, “The loop algorithm,” Adv. Phys. 52, 1–66 and S. L. Sondhi, “Bipartite rokhsar–kivelson points and cantor (2003). deconfinement,” Phys. Rev. B 69, 224415 (2004). [77] A W Sandvik, “A generalization of handscomb’s quantum [94] Stefanos Papanikolaou, Kumar S. Raman, and Eduardo Frad- monte carlo scheme-application to the 1d hubbard model,” kin, “Devil’s staircases, quantum dimer models, and stripe J. Phys. A 25, 3667–3682 (1992). formation in strong coupling models of quantum frustration,” [78] F. Becca and S. Sorella, Quantum Monte Carlo Approaches for Phys. Rev. B 75, 094406 (2007). Correlated Systems (Cambridge University Press, Cambridge, [95] Thiago Schlittler, Thomas Barthel, Grégoire Misguich, Julien UK, 2017). Vidal, and Rémy Mosseri, “Phase diagram of an extended [79] Anders W. Sandvik and Juhani Kurkĳärvi, “Quantum monte quantum dimer model on the hexagonal lattice,” Phys. Rev. carlo simulation method for spin systems,” Phys. Rev. B 43, Lett. 115, 217202 (2015). 5950–5961 (1991). [96] Bowen Zhao, Jun Takahashi, and Anders W. Sandvik, “Mul- [80] Anders W. Sandvik, “Finite-size scaling of the ground-state ticritical deconfined quantum criticality and lifshitz point of parameters of the two-dimensional heisenberg model,” Phys. a helical valence-bond phase,” Phys. Rev. Lett. 125, 257204 Rev. B 56, 11678–11690 (1997). (2020). [81] Matthias Vojta, “Quantum phase transitions,” Rep. Prog. Phys. [97] Zheng Yan, Yongzheng Wu, Chenrong Liu, Olav F. Syljuåsen, 66, 2069–2110 (2003). Jie Lou, and Yan Chen, “Sweeping cluster algorithm for quan- [82] J. Ashkin and E. Teller, “Statistics of two-dimensional lattices tum spin systems with strong geometric restrictions,” Phys. with four components,” Phys. Rev. 64, 178–184 (1943). Rev. B 99, 165135 (2019). [83] Amnon Aharony, “Critical behavior of anisotropic cubic sys- [98] Daniel S. Rokhsar and Steven A. Kivelson, “Superconductivity tems,” Phys. Rev. B 8, 4270–4273 (1973). and the quantum hard-core dimer gas,” Phys. Rev. Lett. 61, [84] José Manuel Carmona, Andrea Pelissetto, and Ettore Vi- 2376–2379 (1988). cari, “푛-component ginzburg-landau hamiltonian with cubic [99] Ashvin Vishwanath, L. Balents, and T. Senthil, “Quantum crit- anisotropy: A six-loop study,” Phys. Rev. B 61, 15136–15151 icality and deconfinement in phase transitions between valence (2000). bond solids,” Phys. Rev. B 69, 224416 (2004). [85] Loran Ts. Adzhemyan, Ella V. Ivanova, Mikhail V. Kompani- [100] Eddy Ardonne, Paul Fendley, and Eduardo Fradkin, “Topolog- ets, Andrey Kudlis, and Aleksandr I. Sokolov, “Six-loop 휀 ical order and conformal quantum critical points,” Ann. Phys. expansion study of three-dimensional n-vector model with cu- (N. Y.) 310, 493 – 551 (2004). bic anisotropy,” Nucl. Phys. B. 940, 332 – 350 (2019). [101] R. Moessner, S. L. Sondhi, and Eduardo Fradkin, “Short- [86] B. Mihura and D. P. Landau, “New type of multicritical be- ranged resonating valence bond physics, quantum dimer mod- havior in a triangular lattice gas model,” Phys. Rev. Lett. 38, els, and ising gauge theories,” Phys. Rev. B 65, 024504 (2001). 977–980 (1977). [102] E. Fradkin, Field Theories of Condensed Matter Physics, Field [87] Eytan Domany, M. Schick, J. S. Walker, and R. B. Griffiths, Theories of Condensed Matter Physics (Cambridge University “Classification of continuous order-disorder transitions in ad- Press, 2013). sorbed monolayers,” Phys. Rev. B 18, 2209–2217 (1978). [103] S. V. Isakov, P. Fendley, A. W. W. Ludwig, S. Trebst, and [88] One could also imagine a scenario in which this order only M. Troyer, “Dynamics at and near conformal quantum critical appears on one or two sublattices, while the other ones remain points,” Phys. Rev. B 83, 125114 (2011). paramagnetic, but this scenario seems unlikely. [104] Roderich Moessner and Kumar S Raman, “Quantum dimer [89] Jie Lou, Anders W. Sandvik, and Leon Balents, “Emergence models,” in Introduction to Frustrated Magnetism (Springer, of u(1) symmetry in the 3d 푥푦 model with 푍푞 anisotropy,” 2011) pp. 437–479. Phys. Rev. Lett. 99, 207203 (2007). [105] Maissam Barkeshli, “Transitions Between Chiral Spin Liquids [90] S. E. Korshunov, “Nature of phase transitions in the striped and Z2 Spin Liquids,” arXiv:1307.8194 (2013). phase of a triangular-lattice ising antiferromagnet,” Phys. Rev. [106] Long Zhang and Fa Wang, “Unconventional surface critical B 72, 144417 (2005). behavior induced by a quantum phase transition from the two-dimensional affleck-kennedy-lieb-tasaki phase to a néel- ordered phase,” Phys. Rev. Lett. 118, 087201 (2017).