Quantum Monte Carlo Methods for Fermion-Boson problems
Fakher F. Assaad (Autumn School on Correlated Electrons, Jülich 20th September 2018)
Outline: Ø Fermion-Boson problems Electron-Phonon, Su-Schrieffer-Heeger (SSH) Unconstrained Gauge theories “De-signer” Hamiltonians
Ø Auxiliary field Quantum Monte Carlo (QMC), Generalities
Ø Application of the Auxiliary field QMC to the SSH model Sign-free formulation Sampling, Hybrid Monte Carlo
Ø Conclusions
SFB1170 ToCoTronics Quantum Monte Carlo Methods for Fermion-Boson problems
Many thanks to
M. Hohenadler J. Hofmann S. Beyl M. Rackowski Z. Wang T. Sato
M. Ulybyshev F. Parisen Toldin J. Schwab E. Huffman Fermion-Boson problems: The Su-Schrieffer-Heeger (SSH) model Su-Schrieffer-Heeger model (SSH)
ˆ2 P i,j k ˆ ˆ h i ˆ 2 H = t cˆi,† cˆj, +ˆcj,† cˆi, + g X i,j cˆi,† cˆj, +ˆcj,† cˆi, + + X i,j h i 2m 2 h i i,j , i,j , i,j " # AAADz3icvVNNa9swGFbsfXTZV7oddxELg452wQ6D9VIo9NLDDiksbSByjazIjhrJNpLcEjSVXff3dtsP2P+YYiclawJjl70gePS8z/uJlJScKR0EP1ue/+Dho8c7T9pPnz1/8bKz++pcFZUkdEgKXshRghXlLKdDzTSno1JSLBJOL5LZycJ/cU2lYkX+Rc9LGgmc5SxlBGtHxbutX2iKtTm18Ah+0EhVIjaI01QjjvOMU8gOrpBk2VQj2RAHSLFMYFur9upoYmPDVvylQROcZVTaO9/Vyre/SW2R36VqCr+H+zD7l85gnWlk/6b/vyOsJlhrpoZNF2OUSkxMnWBgt+ouTd9a0xeuhVo7cxe7NurWgKZ+FHe6QS+oDW6CcAm6YGmDuPMDTQpSCZprwrFS4zAodWSw1Iy41G1UKVpiMsMZHTuYY0FVZOr3aOE7x0xgWkh3cg1rdj3CYKHUXCROKbCeqvu+BbnNN650ehgZlpeVpjlpCqUVh7qAi8cNJ0xSovncAUwkc71CMsVuWdp9gbZbQnh/5E1w3u+FQS88+9g97i/XsQPegLdgD4TgEzgGp2AAhoB4nz3pGe+rf+bf+Lf+t0bqtZYxr8Ef5n//DTqQSdo= hXi ⇣ ⌘ hXi ⇣ ⌘ hXi
cˆ†
AAAC5nichVFNSxxBEH1ONK5G42qOelhcAjnIMiuCehM0wYug4KrgmqVnth2bnS96eoV12Ev+QE7iNcdc9d/ob/Hgm3YUVMQeeqr61avXVV1eGqrMuO7tiPNpdOzzeGVi8svU9NeZ6uzcQZb0tS9bfhIm+sgTmQxVLFtGmVAepVqKyAvlodfbLOKH51JnKon3zSCVJ5EIYnWqfGEIdaoL7TNhcn/4O293RRBIPezkaqmdqSASw0617jZcu2pvnWbp1FGu3aR6hza6SOCjjwgSMQz9EAIZv2M04SIldoKcmKanbFxiiEnm9smSZAiiPf4Dno5LNOa50Mxsts9bQm7NzBq+c/+yih7Zxa2SfkZ7z31hseDdG3KrXFQ4oPWoOGEVd4gbnJHxUWZUMp9q+Tiz6MrgFGu2G8X6UosUffrPOluMaGI9G6nhp2UG1PDs+ZwvENO2WEHxyk8KNdtxl1ZYK61KXCoK6mna4vVZD8fcfD3Ut05rubHeaO6t1DeWy3lXMI9F/OBQV7GBbeyyDB9/8B/XuHGU89e5dK4eqc5ImfMNL5bz7wHWc5r9 i,
ˆ k ˆ ˆ X i,j !0 = Xb, Pb0 = i~ b,b0 AAACB3icbZDLSsNAFIYn9VbrLepSkMEiuJCSFEGXBTcuK9gLNCFMpift2MkkzEyEErpz46u4caGIW1/BnW/jtM1CW38Y+PjPOZw5f5hyprTjfFulldW19Y3yZmVre2d3z94/aKskkxRaNOGJ7IZEAWcCWpppDt1UAolDDp1wdD2tdx5AKpaIOz1OwY/JQLCIUaKNFdjH3pDovDsJco8TMeCA2fk99uScJ4FddWrOTHgZ3AKqqFAzsL+8fkKzGISmnCjVc51U+zmRmlEOk4qXKUgJHZEB9AwKEoPy89kdE3xqnD6OEmme0Hjm/p7ISazUOA5NZ0z0UC3WpuZ/tV6moys/ZyLNNAg6XxRlHOsET0PBfSaBaj42QKhk5q+YDokkVJvoKiYEd/HkZWjXa65Tc28vqo16EUcZHaETdIZcdIka6AY1UQtR9Iie0St6s56sF+vd+pi3lqxi5hD9kfX5A5H5mQo= h i AAAC6HichVFNa9tAEH1W08Z2P+I2x0IRNYGcjBwKTQ4BQ5KSS8GFOAnEwayUtSr0mdU64Ard+gd6Cr3mmGvyZ9rf0kOfNnKhDSUrVjP75s3bmR03i4JcO86PhvVo6fGT5War/fTZ8xcrnZevDvN0pjw58tIoVceuyGUUJHKkAx3J40xJEbuRPHLDnSp+dCFVHqTJgZ5n8jQWfhJMA09oQpPOm3EaS19MCqfctsf5udLFeKqEV4RlEZflpNN1eo5Z9n2nXztd1GuYdn5ijDOk8DBDDIkEmn4EgZzfCfpwkBE7RUFM0QtMXKJEm7kzsiQZgmjIv8/TSY0mPFeaucn2eEvErZhpY437g1F0ya5ulfRz2l/cXwzm//eGwihXFc5pXSq2jOJH4hqfyXgoM66Zi1oezqy60phi03QTsL7MIFWf3h+dXUYUsdBEbOwZpk8N15wv+AIJ7YgVVK+8ULBNx2e0wlhpVJJaUVBP0Vavz3o45v6/Q73vjDZ6W73+p3fdwUY97yZe4y3WOdT3GGAfQ5bh4SuucYNbK7S+WZfW9zuq1ahzVvHXsq5+A/dpm88= rm AAADCHichVHLahRBFD3p+Ejia9Slm8JB4iIM3UFQkEDAB26EEZxkYHoYqnpqeorp6W6qawJxmB/Ib/gD7iTb7NzqSr/FhafLjqBBUk31PXXuvafuravKzFQuDL+vBetXrl67vrG5dePmrdt3WnfvHVTFwia6lxRZYftKVjozue454zLdL62Wc5XpQzV7UfsPj7StTJG/d8elHs5lmpuJSaQjNWrtxZmeuIGIp9It+6uR2vGoS7QdW5NO3VDsCRNPlbQifi7isc6cHC3VjtpeiVGrHXZCv8RFEDWgjWZ1i9YPxBijQIIF5tDI4YgzSFT8BogQoiQ3xJKcJTLer7HCFnMXjNKMkGRn/Kc8DRo257nWrHx2wlsybstMgUfcr72iYnR9qyauaH9yf/Bc+t8bll65rvCYVlFx0yu+Je8wZcRlmfMm8ryWyzPrrhwmeOa7Mayv9EzdZ/JH5yU9ltzMewRe+ciUGsqfj/gCOW2PFdSvfK4gfMdjWumt9ip5oyipZ2nr12c9HHP071AvgoPdThR2ondP2vthM/ANPMBDPOZUn2Ifb9BlHQk+4gu+4ltwEnwKPgenv0ODtSbnPv5awdkv7p6mqg== h i
Limiting cases
AAAB83icbVDLSgNBEOyNrxhfUY9eBoPgKcyKoBch4MVjBBMD2SXMTmaTIfNYZmaFEPIbXjwo4tWf8ebfOEn2oIkFDUVVN91dSSa4dRh/B6W19Y3NrfJ2ZWd3b/+genjUtjo3lLWoFtp0EmKZ4Iq1HHeCdTLDiEwEe0xGtzP/8YkZy7V6cOOMxZIMFE85Jc5LUaQlG5AeRjcI96o1XMdzoFUSFqQGBZq96lfU1zSXTDkqiLXdEGcunhDjOBVsWolyyzJCR2TAup4qIpmNJ/Obp+jMK32UauNLOTRXf09MiLR2LBPfKYkb2mVvJv7ndXOXXscTrrLcMUUXi9JcIKfRLADU54ZRJ8aeEGq4vxXRITGEOh9TxYcQLr+8StoX9RDXw/vLWgMXcZThBE7hHEK4ggbcQRNaQCGDZ3iFtyAPXoL34GPRWgqKmWP4g+DzB1mJkIE= ! =0 !AAAC6HichVFNa9tAEH1W+mG7H3HTYyCImkJPRjaFprdAmtJLwIU4DsTGrOS1vFhfrNYJrvGtfyCnkmuOvTZ/Jv0tPfRpIxfaULJiNbNv3ryd2fGzSOXG824qzsaDh48eV2v1J0+fPd9svNg6ztO5DmQvSKNUn/gil5FKZM8oE8mTTEsR+5Hs+7P9It4/kzpXaXJkFpkcxiJM1EQFwhAaNXYGaSxDMVp6K9cdaBVOjdA6PXcHKpmYxajR9FqeXe5dp106TZSrmzZ+YoAxUgSYI4ZEAkM/gkDO7xRteMiIDbEkpukpG5dYoc7cOVmSDEF0xn/I02mJJjwXmrnNDnhLxK2Z6eI190er6JNd3Crp57S/uL9YLPzvDUurXFS4oPWpWLOKh8QNpmTclxmXzHUt92cWXRlMsGu7Uawvs0jRZ/BH5wMjmtjMRlwcWGZIDd+ez/gCCW2PFRSvvFZwbcdjWmGttCpJqSiop2mL12c9HHP736HedXqd1vtW+/Pb5l6nnHcV23iFNxzqO+zhE7osI8BXfMcPXDsz58L55lzeUp1KmfMSfy3n6jdMYZsf 0 0 !1 Rev. Mod. Phys. 60, 781 (1988) Perfect nesting and Van-Hove singularity13.03.18 HQMC simulations of the SSH model 5 ˆ W H = t cˆi,† cˆj, +ˆcj,† cˆi, (QQ, i⌦ = 0) ln2 m ordered valence/ bondk statesT (VBS) i,j , AAACQXicbZBPSxtBGMZntVZNbZvq0cvQUIhQwm4o6KUgevGmQmIC2bjMTt5NhszMLjPvCmHZr+al36A37148KOLVi5M/gjV9YeDheZ6XmfnFmRQWff/GW1n9sPZxfWOz8mnr85ev1W/bFzbNDYc2T2VqujGzIIWGNgqU0M0MMBVL6MTj42neuQJjRapbOMmgr9hQi0Rwhs6Kqt2QjwQNJSRYp2Gm4uK8/Emdc6pgyCJFf1OfhkYMR7jncpNmmFLX15fN16XEMF50ymIcHdFW+VqOqjW/4c+GLotgIWpkMWdR9W84SHmuQCOXzNpe4GfYL5hBwSWUlTC3kDE+ZkPoOamZAtsvZgRK+sM5A5qkxh2NdOa+3SiYsnaiYtdUDEf2fTY1/5f1ckwO+oXQWY6g+fyiJJfUUZjipANhgKOcOMG4Ee6tlI+YI4IOesVBCN5/eVlcNBuB3wjOf9UOmwscG2SXfCd1EpB9ckhOyBlpE06uyS25Jw/eH+/Oe/Se5tUVb7GzQ/4Z7/kFJzauOw== sha1_base64="hP+6LrUf2d3tZaldqaQQvEKMXyw=">AAAB2XicbZDNSgMxFIXv1L86Vq1rN8EiuCozbnQpuHFZwbZCO5RM5k4bmskMyR2hDH0BF25EfC93vo3pz0JbDwQ+zknIvSculLQUBN9ebWd3b/+gfugfNfzjk9Nmo2fz0gjsilzl5jnmFpXU2CVJCp8LgzyLFfbj6f0i77+gsTLXTzQrMMr4WMtUCk7O6oyaraAdLMW2IVxDC9YaNb+GSS7KDDUJxa0dhEFBUcUNSaFw7g9LiwUXUz7GgUPNM7RRtRxzzi6dk7A0N+5oYkv394uKZ9bOstjdzDhN7Ga2MP/LBiWlt1EldVESarH6KC0Vo5wtdmaJNChIzRxwYaSblYkJN1yQa8Z3HYSbG29D77odBu3wMYA6nMMFXEEIN3AHD9CBLghI4BXevYn35n2suqp569LO4I+8zx84xIo4 sha1_base64="E469Am/V+esbdjPBFF4sNuTXhCI=">AAACNnicbZBPSxtBGMbftVVj/Bd79TJUCgoSdr3opSD10psKxgSycZmdvJsMmZldZt4thGW/Wi/9Br317sWDIn4CJzFCq31g4OF53peZ+aWFko7C8E+w9OHj8spqY625vrG5td3a2bh2eWkFdkSucttLuUMlDXZIksJeYZHrVGE3nZzN+u4PtE7m5oqmBQ40HxmZScHJR0mrF4uxZLHCjPZZXOi0uqwPmU/ONY54otlXFrLYytGYDnxv84Jy5ufNzdHrUma5qLp1NUm+sav6dThp7YXtcC723kQLswcLXSSt3/EwF6VGQ0Jx5/pRWNCg4pakUFg349JhwcWEj7DvreEa3aCaE6jZF58MWZZbfwyxefr3RsW1c1Od+knNaezedrPwf12/pOxkUElTlIRGvFyUlYp5CjOcbCgtClJTb7iw0r+ViTH3RMhDb3oI0dsvvzfXR+0obEeXITRgFz7DPkRwDKfwHS6gAwJ+wi3cw0PwK7gLHl9wLQULbp/gHwVPz7jbrJ8= sha1_base64="a22j3uBbepvqkCSkzfJ5mvCp9vw=">AAACQXicbZBLaxsxFIU1Sdok7iPTdNmNiCm4UMyMN8kmEJpNdonBL/C4g0a+YwtLmkG6UzDD/LVu8g+yyz6bLFpKt91UfgQSJxcEh3PORdKX5FJYDIJbb2t759Xr3b392pu3794f+B8OezYrDIcuz2RmBgmzIIWGLgqUMMgNMJVI6Cez80Xe/wHGikx3cJ7DSLGJFqngDJ0V+4OITwWNJKTYoFGukrJdfaXOuVQwYbGipzSgkRGTKX5xuclyzKjr6++th6XUMF72q3IWf6Od6qEc+/WgGSyHPhfhWtTJeq5i/yYaZ7xQoJFLZu0wDHIclcyg4BKqWlRYyBmfsQkMndRMgR2VSwIV/eycMU0z445GunQfb5RMWTtXiWsqhlO7mS3Ml7JhgenJqBQ6LxA0X12UFpI6CgucdCwMcJRzJxg3wr2V8ilzRNBBrzkI4eaXn4teqxkGzbAd1M9aaxx75BM5Ig0SkmNyRi7IFekSTn6SO/KL/PauvXvvj/d3Vd3y1jsfyZPx/v0HJfauNw== B hXi ⇣ ⌘ ✓ ◆ 2 g2 cˆ† cˆ +ˆc† cˆ 4k i, j, j, i, strong bonds weak bonds i,j ! hXi X g2 ˆ ˆ i,j S iS j +⌘ˆi⌘ˆj
AAAD/nicvVPLbtNAFJ3aPEp4pcCOzYgIVGgb2RYSbJAqsemyCNJWyiTWeHztTGKPzcwYKRpZ4lfYsAAhtnwHO/6GieNWIemiK65l+fjcM/eeeUVlxpX2vD9bjnvt+o2b27c6t+/cvXe/u/PgRBWVZDBgRVbIs4gqyLiAgeY6g7NSAs2jDE6j2dtF/vQTSMUL8UHPSxjlNBU84YxqS4U7ziMSQcqFgY+CSknnL2oyodoc1RjjNxg/sw8+0ERVeWhIBokmGRVpBpjvT4nk6UQTuST2ieJpTutGtdtUYXVo+Dk/NiSmaQqyvshNz3N7m9Ql8otSy8bPrUVMyNJkx8JKxCAjSRmYA5LYr0nHQW1ezuor+G9941baNvoPsxgH2LIrfq9mFzcFSZlH5n1dh3ztf4r3VhgCmq6JWmqKFztNQMQrRyDs9ry+1wTeBH4LeqiN47D7m8QFq3IQmmVUqaHvlXpkqNScWa8dUikoKZvRFIYWCpqDGpnm+Nb4qWVinBTSvkLjhl0dYWiu1DyPrDKneqLWcwvystyw0snrkeGirDQItmyUVBnWBV7cBRxzCUxncwsok9x6xWxC7T5oe2M6dhH89SlvgpOg73t9/13QOwza5dhGj9ETtIt89AodoiN0jAaIOcb54nxzvruf3a/uD/fnUupstWMeon/C/fUXLnBZHw== k Valence Bond Solid h i | Quantum antiferromagnet/P {z SSC + CDW}
Tang & Hirsch, Phys. Rev. B 37, 9546 (1988) 13.03.18 HQMC simulations of the SSH model 9 Fermion-Boson problems: The Su-Schrieffer-Heeger (SSH) model Su-Schrieffer-Heeger model (SSH)
ˆ2 P i,j k ˆ ˆ h i ˆ 2 H = t cˆi,† cˆj, +ˆcj,† cˆi, + g X i,j cˆi,† cˆj, +ˆcj,† cˆi, + + X i,j h i 2m 2 h i i,j , i,j , i,j " # AAADz3icvVNNa9swGFbsfXTZV7oddxELg452wQ6D9VIo9NLDDiksbSByjazIjhrJNpLcEjSVXff3dtsP2P+YYiclawJjl70gePS8z/uJlJScKR0EP1ue/+Dho8c7T9pPnz1/8bKz++pcFZUkdEgKXshRghXlLKdDzTSno1JSLBJOL5LZycJ/cU2lYkX+Rc9LGgmc5SxlBGtHxbutX2iKtTm18Ah+0EhVIjaI01QjjvOMU8gOrpBk2VQj2RAHSLFMYFur9upoYmPDVvylQROcZVTaO9/Vyre/SW2R36VqCr+H+zD7l85gnWlk/6b/vyOsJlhrpoZNF2OUSkxMnWBgt+ouTd9a0xeuhVo7cxe7NurWgKZ+FHe6QS+oDW6CcAm6YGmDuPMDTQpSCZprwrFS4zAodWSw1Iy41G1UKVpiMsMZHTuYY0FVZOr3aOE7x0xgWkh3cg1rdj3CYKHUXCROKbCeqvu+BbnNN650ehgZlpeVpjlpCqUVh7qAi8cNJ0xSovncAUwkc71CMsVuWdp9gbZbQnh/5E1w3u+FQS88+9g97i/XsQPegLdgD4TgEzgGp2AAhoB4nz3pGe+rf+bf+Lf+t0bqtZYxr8Ef5n//DTqQSdo= hXi ⇣ ⌘ hXi ⇣ ⌘ hXi
cˆ†
AAAC5nichVFNSxxBEH1ONK5G42qOelhcAjnIMiuCehM0wYug4KrgmqVnth2bnS96eoV12Ev+QE7iNcdc9d/ob/Hgm3YUVMQeeqr61avXVV1eGqrMuO7tiPNpdOzzeGVi8svU9NeZ6uzcQZb0tS9bfhIm+sgTmQxVLFtGmVAepVqKyAvlodfbLOKH51JnKon3zSCVJ5EIYnWqfGEIdaoL7TNhcn/4O293RRBIPezkaqmdqSASw0617jZcu2pvnWbp1FGu3aR6hza6SOCjjwgSMQz9EAIZv2M04SIldoKcmKanbFxiiEnm9smSZAiiPf4Dno5LNOa50Mxsts9bQm7NzBq+c/+yih7Zxa2SfkZ7z31hseDdG3KrXFQ4oPWoOGEVd4gbnJHxUWZUMp9q+Tiz6MrgFGu2G8X6UosUffrPOluMaGI9G6nhp2UG1PDs+ZwvENO2WEHxyk8KNdtxl1ZYK61KXCoK6mna4vVZD8fcfD3Ut05rubHeaO6t1DeWy3lXMI9F/OBQV7GBbeyyDB9/8B/XuHGU89e5dK4eqc5ImfMNL5bz7wHWc5r9 i,
ˆ k ˆ ˆ X i,j !0 = Xb, Pb0 = i~ b,b0 AAACB3icbZDLSsNAFIYn9VbrLepSkMEiuJCSFEGXBTcuK9gLNCFMpift2MkkzEyEErpz46u4caGIW1/BnW/jtM1CW38Y+PjPOZw5f5hyprTjfFulldW19Y3yZmVre2d3z94/aKskkxRaNOGJ7IZEAWcCWpppDt1UAolDDp1wdD2tdx5AKpaIOz1OwY/JQLCIUaKNFdjH3pDovDsJco8TMeCA2fk99uScJ4FddWrOTHgZ3AKqqFAzsL+8fkKzGISmnCjVc51U+zmRmlEOk4qXKUgJHZEB9AwKEoPy89kdE3xqnD6OEmme0Hjm/p7ISazUOA5NZ0z0UC3WpuZ/tV6moys/ZyLNNAg6XxRlHOsET0PBfSaBaj42QKhk5q+YDokkVJvoKiYEd/HkZWjXa65Tc28vqo16EUcZHaETdIZcdIka6AY1UQtR9Iie0St6s56sF+vd+pi3lqxi5hD9kfX5A5H5mQo= h i AAAC6HichVFNa9tAEH1W08Z2P+I2x0IRNYGcjBwKTQ4BQ5KSS8GFOAnEwayUtSr0mdU64Ard+gd6Cr3mmGvyZ9rf0kOfNnKhDSUrVjP75s3bmR03i4JcO86PhvVo6fGT5War/fTZ8xcrnZevDvN0pjw58tIoVceuyGUUJHKkAx3J40xJEbuRPHLDnSp+dCFVHqTJgZ5n8jQWfhJMA09oQpPOm3EaS19MCqfctsf5udLFeKqEV4RlEZflpNN1eo5Z9n2nXztd1GuYdn5ijDOk8DBDDIkEmn4EgZzfCfpwkBE7RUFM0QtMXKJEm7kzsiQZgmjIv8/TSY0mPFeaucn2eEvErZhpY437g1F0ya5ulfRz2l/cXwzm//eGwihXFc5pXSq2jOJH4hqfyXgoM66Zi1oezqy60phi03QTsL7MIFWf3h+dXUYUsdBEbOwZpk8N15wv+AIJ7YgVVK+8ULBNx2e0wlhpVJJaUVBP0Vavz3o45v6/Q73vjDZ6W73+p3fdwUY97yZe4y3WOdT3GGAfQ5bh4SuucYNbK7S+WZfW9zuq1ahzVvHXsq5+A/dpm88= rm AAADCHichVHLahRBFD3p+Ejia9Slm8JB4iIM3UFQkEDAB26EEZxkYHoYqnpqeorp6W6qawJxmB/Ib/gD7iTb7NzqSr/FhafLjqBBUk31PXXuvafuravKzFQuDL+vBetXrl67vrG5dePmrdt3WnfvHVTFwia6lxRZYftKVjozue454zLdL62Wc5XpQzV7UfsPj7StTJG/d8elHs5lmpuJSaQjNWrtxZmeuIGIp9It+6uR2vGoS7QdW5NO3VDsCRNPlbQifi7isc6cHC3VjtpeiVGrHXZCv8RFEDWgjWZ1i9YPxBijQIIF5tDI4YgzSFT8BogQoiQ3xJKcJTLer7HCFnMXjNKMkGRn/Kc8DRo257nWrHx2wlsybstMgUfcr72iYnR9qyauaH9yf/Bc+t8bll65rvCYVlFx0yu+Je8wZcRlmfMm8ryWyzPrrhwmeOa7Mayv9EzdZ/JH5yU9ltzMewRe+ciUGsqfj/gCOW2PFdSvfK4gfMdjWumt9ip5oyipZ2nr12c9HHP071AvgoPdThR2ondP2vthM/ANPMBDPOZUn2Ifb9BlHQk+4gu+4ltwEnwKPgenv0ODtSbnPv5awdkv7p6mqg== h i
Limiting cases
AAAB83icbVDLSgNBEOyNrxhfUY9eBoPgKcyKoBch4MVjBBMD2SXMTmaTIfNYZmaFEPIbXjwo4tWf8ebfOEn2oIkFDUVVN91dSSa4dRh/B6W19Y3NrfJ2ZWd3b/+genjUtjo3lLWoFtp0EmKZ4Iq1HHeCdTLDiEwEe0xGtzP/8YkZy7V6cOOMxZIMFE85Jc5LUaQlG5AeRjcI96o1XMdzoFUSFqQGBZq96lfU1zSXTDkqiLXdEGcunhDjOBVsWolyyzJCR2TAup4qIpmNJ/Obp+jMK32UauNLOTRXf09MiLR2LBPfKYkb2mVvJv7ndXOXXscTrrLcMUUXi9JcIKfRLADU54ZRJ8aeEGq4vxXRITGEOh9TxYcQLr+8StoX9RDXw/vLWgMXcZThBE7hHEK4ggbcQRNaQCGDZ3iFtyAPXoL34GPRWgqKmWP4g+DzB1mJkIE= ! =0 !AAAC6HichVFNa9tAEH1W+mG7H3HTYyCImkJPRjaFprdAmtJLwIU4DsTGrOS1vFhfrNYJrvGtfyCnkmuOvTZ/Jv0tPfRpIxfaULJiNbNv3ryd2fGzSOXG824qzsaDh48eV2v1J0+fPd9svNg6ztO5DmQvSKNUn/gil5FKZM8oE8mTTEsR+5Hs+7P9It4/kzpXaXJkFpkcxiJM1EQFwhAaNXYGaSxDMVp6K9cdaBVOjdA6PXcHKpmYxajR9FqeXe5dp106TZSrmzZ+YoAxUgSYI4ZEAkM/gkDO7xRteMiIDbEkpukpG5dYoc7cOVmSDEF0xn/I02mJJjwXmrnNDnhLxK2Z6eI190er6JNd3Crp57S/uL9YLPzvDUurXFS4oPWpWLOKh8QNpmTclxmXzHUt92cWXRlMsGu7Uawvs0jRZ/BH5wMjmtjMRlwcWGZIDd+ez/gCCW2PFRSvvFZwbcdjWmGttCpJqSiop2mL12c9HHP736HedXqd1vtW+/Pb5l6nnHcV23iFNxzqO+zhE7osI8BXfMcPXDsz58L55lzeUp1KmfMSfy3n6jdMYZsf 0 0 !1 Rev. Mod. Phys. 60, 781 (1988) Perfect nesting and Van-Hove singularity13.03.18 HQMC simulations of the SSH model 5 ˆ W H = t cˆi,† cˆj, +ˆcj,† cˆi, (QQ, i⌦ = 0) ln2 m ordered valence/ bondk statesT (VBS) i,j , AAACQXicbZBPSxtBGMZntVZNbZvq0cvQUIhQwm4o6KUgevGmQmIC2bjMTt5NhszMLjPvCmHZr+al36A37148KOLVi5M/gjV9YeDheZ6XmfnFmRQWff/GW1n9sPZxfWOz8mnr85ev1W/bFzbNDYc2T2VqujGzIIWGNgqU0M0MMBVL6MTj42neuQJjRapbOMmgr9hQi0Rwhs6Kqt2QjwQNJSRYp2Gm4uK8/Emdc6pgyCJFf1OfhkYMR7jncpNmmFLX15fN16XEMF50ymIcHdFW+VqOqjW/4c+GLotgIWpkMWdR9W84SHmuQCOXzNpe4GfYL5hBwSWUlTC3kDE+ZkPoOamZAtsvZgRK+sM5A5qkxh2NdOa+3SiYsnaiYtdUDEf2fTY1/5f1ckwO+oXQWY6g+fyiJJfUUZjipANhgKOcOMG4Ee6tlI+YI4IOesVBCN5/eVlcNBuB3wjOf9UOmwscG2SXfCd1EpB9ckhOyBlpE06uyS25Jw/eH+/Oe/Se5tUVb7GzQ/4Z7/kFJzauOw== sha1_base64="hP+6LrUf2d3tZaldqaQQvEKMXyw=">AAAB2XicbZDNSgMxFIXv1L86Vq1rN8EiuCozbnQpuHFZwbZCO5RM5k4bmskMyR2hDH0BF25EfC93vo3pz0JbDwQ+zknIvSculLQUBN9ebWd3b/+gfugfNfzjk9Nmo2fz0gjsilzl5jnmFpXU2CVJCp8LgzyLFfbj6f0i77+gsTLXTzQrMMr4WMtUCk7O6oyaraAdLMW2IVxDC9YaNb+GSS7KDDUJxa0dhEFBUcUNSaFw7g9LiwUXUz7GgUPNM7RRtRxzzi6dk7A0N+5oYkv394uKZ9bOstjdzDhN7Ga2MP/LBiWlt1EldVESarH6KC0Vo5wtdmaJNChIzRxwYaSblYkJN1yQa8Z3HYSbG29D77odBu3wMYA6nMMFXEEIN3AHD9CBLghI4BXevYn35n2suqp569LO4I+8zx84xIo4 sha1_base64="E469Am/V+esbdjPBFF4sNuTXhCI=">AAACNnicbZBPSxtBGMbftVVj/Bd79TJUCgoSdr3opSD10psKxgSycZmdvJsMmZldZt4thGW/Wi/9Br317sWDIn4CJzFCq31g4OF53peZ+aWFko7C8E+w9OHj8spqY625vrG5td3a2bh2eWkFdkSucttLuUMlDXZIksJeYZHrVGE3nZzN+u4PtE7m5oqmBQ40HxmZScHJR0mrF4uxZLHCjPZZXOi0uqwPmU/ONY54otlXFrLYytGYDnxv84Jy5ufNzdHrUma5qLp1NUm+sav6dThp7YXtcC723kQLswcLXSSt3/EwF6VGQ0Jx5/pRWNCg4pakUFg349JhwcWEj7DvreEa3aCaE6jZF58MWZZbfwyxefr3RsW1c1Od+knNaezedrPwf12/pOxkUElTlIRGvFyUlYp5CjOcbCgtClJTb7iw0r+ViTH3RMhDb3oI0dsvvzfXR+0obEeXITRgFz7DPkRwDKfwHS6gAwJ+wi3cw0PwK7gLHl9wLQULbp/gHwVPz7jbrJ8= sha1_base64="a22j3uBbepvqkCSkzfJ5mvCp9vw=">AAACQXicbZBLaxsxFIU1Sdok7iPTdNmNiCm4UMyMN8kmEJpNdonBL/C4g0a+YwtLmkG6UzDD/LVu8g+yyz6bLFpKt91UfgQSJxcEh3PORdKX5FJYDIJbb2t759Xr3b392pu3794f+B8OezYrDIcuz2RmBgmzIIWGLgqUMMgNMJVI6Cez80Xe/wHGikx3cJ7DSLGJFqngDJ0V+4OITwWNJKTYoFGukrJdfaXOuVQwYbGipzSgkRGTKX5xuclyzKjr6++th6XUMF72q3IWf6Od6qEc+/WgGSyHPhfhWtTJeq5i/yYaZ7xQoJFLZu0wDHIclcyg4BKqWlRYyBmfsQkMndRMgR2VSwIV/eycMU0z445GunQfb5RMWTtXiWsqhlO7mS3Ml7JhgenJqBQ6LxA0X12UFpI6CgucdCwMcJRzJxg3wr2V8ilzRNBBrzkI4eaXn4teqxkGzbAd1M9aaxx75BM5Ig0SkmNyRi7IFekSTn6SO/KL/PauvXvvj/d3Vd3y1jsfyZPx/v0HJfauNw== B hXi ⇣ ⌘ ✓ ◆ 2 g2 cˆ† cˆ +ˆc† cˆ 4k i, j, j, i, strong bonds weak bonds i,j ! hXi X g2 ˆ ˆ i,j S iS j +⌘ˆi⌘ˆj
AAAD/nicvVPLbtNAFJ3aPEp4pcCOzYgIVGgb2RYSbJAqsemyCNJWyiTWeHztTGKPzcwYKRpZ4lfYsAAhtnwHO/6GieNWIemiK65l+fjcM/eeeUVlxpX2vD9bjnvt+o2b27c6t+/cvXe/u/PgRBWVZDBgRVbIs4gqyLiAgeY6g7NSAs2jDE6j2dtF/vQTSMUL8UHPSxjlNBU84YxqS4U7ziMSQcqFgY+CSknnL2oyodoc1RjjNxg/sw8+0ERVeWhIBokmGRVpBpjvT4nk6UQTuST2ieJpTutGtdtUYXVo+Dk/NiSmaQqyvshNz3N7m9Ql8otSy8bPrUVMyNJkx8JKxCAjSRmYA5LYr0nHQW1ezuor+G9941baNvoPsxgH2LIrfq9mFzcFSZlH5n1dh3ztf4r3VhgCmq6JWmqKFztNQMQrRyDs9ry+1wTeBH4LeqiN47D7m8QFq3IQmmVUqaHvlXpkqNScWa8dUikoKZvRFIYWCpqDGpnm+Nb4qWVinBTSvkLjhl0dYWiu1DyPrDKneqLWcwvystyw0snrkeGirDQItmyUVBnWBV7cBRxzCUxncwsok9x6xWxC7T5oe2M6dhH89SlvgpOg73t9/13QOwza5dhGj9ETtIt89AodoiN0jAaIOcb54nxzvruf3a/uD/fnUupstWMeon/C/fUXLnBZHw== k Valence Bond Solid h i | Quantum antiferromagnet/P {z SSC + CDW}
Tang & Hirsch, Phys. Rev. B 37, 9546 (1988) 13.03.18 HQMC simulations of the SSH model 9 Fermion-Boson problems: The Su-Schrieffer-Heeger (SSH) model
Symmetries ˆ2 P i,j k ˆ ˆ h i ˆ 2 H = t cˆi,† cˆj, +ˆcj,† cˆi, + g X i,j cˆi,† cˆj, +ˆcj,† cˆi, + + X i,j h i 2m 2 h i i,j , i,j , i,j " # AAADz3icvVNNa9swGFbsfXTZV7oddxELg452wQ6D9VIo9NLDDiksbSByjazIjhrJNpLcEjSVXff3dtsP2P+YYiclawJjl70gePS8z/uJlJScKR0EP1ue/+Dho8c7T9pPnz1/8bKz++pcFZUkdEgKXshRghXlLKdDzTSno1JSLBJOL5LZycJ/cU2lYkX+Rc9LGgmc5SxlBGtHxbutX2iKtTm18Ah+0EhVIjaI01QjjvOMU8gOrpBk2VQj2RAHSLFMYFur9upoYmPDVvylQROcZVTaO9/Vyre/SW2R36VqCr+H+zD7l85gnWlk/6b/vyOsJlhrpoZNF2OUSkxMnWBgt+ouTd9a0xeuhVo7cxe7NurWgKZ+FHe6QS+oDW6CcAm6YGmDuPMDTQpSCZprwrFS4zAodWSw1Iy41G1UKVpiMsMZHTuYY0FVZOr3aOE7x0xgWkh3cg1rdj3CYKHUXCROKbCeqvu+BbnNN650ehgZlpeVpjlpCqUVh7qAi8cNJ0xSovncAUwkc71CMsVuWdp9gbZbQnh/5E1w3u+FQS88+9g97i/XsQPegLdgD4TgEzgGp2AAhoB4nz3pGe+rf+bf+Lf+t0bqtZYxr8Ef5n//DTqQSdo= hXi ⇣ ⌘ hXi ⇣ ⌘ hXi
Partial particle-hole symmetry ˆ 1 ˆ iQQ i P↵ cˆii,† P↵ = ↵, e · cˆii, +(1 ↵, )ˆcii,† AAAC03icjVJLaxRBEO4ZX3F9rXr0UrgIEZNlJgh6EQJePG7ATQLbm6Wmp2e2Sc+D7prA0gyKePXPefMv+CvsmUxANwoWNHz9VX1V1VWd1FpZiqIfQXjj5q3bd3buju7df/Dw0fjxk2NbNUbIuah0ZU4TtFKrUs5JkZantZFYJFqeJOfvO//JhTRWVeVH2tRyWWBeqkwJJE+txj/5GsnN2pXjqOs1tmduP24Belr4G08xz6XpAuoicardA55Iwha2lfAOeCo14RWzNwTKM6egVx+1XKQVwZDK17kqtJ0eOt8r4FpmtBvDPvwjNzcqX9PL/+l3NZ5E06g3uA7iAUzYYLPV+DtPK9EUsiSh0dpFHNW0dGhICS3bEW+srFGcYy4XHpZYSLt0/U5aeOGZFLLK+FMS9OzvCoeFtZsi8ZEF0tpu+zryb75FQ9nbpVNl3ZAsxWWhrNFAFXQLhlQZKUhvPEBhlO8VxBoNCvLfYOSHEG8/+To4PpjG0TQ+ej05PBjGscOesedsl8XsDTtkH9iMzZkIZsFF8Cn4HM5DF34Jv16GhsGgecr+sPDbL0Kz43E= sha1_base64="hP+6LrUf2d3tZaldqaQQvEKMXyw=">AAAB2XicbZDNSgMxFIXv1L86Vq1rN8EiuCozbnQpuHFZwbZCO5RM5k4bmskMyR2hDH0BF25EfC93vo3pz0JbDwQ+zknIvSculLQUBN9ebWd3b/+gfugfNfzjk9Nmo2fz0gjsilzl5jnmFpXU2CVJCp8LgzyLFfbj6f0i77+gsTLXTzQrMMr4WMtUCk7O6oyaraAdLMW2IVxDC9YaNb+GSS7KDDUJxa0dhEFBUcUNSaFw7g9LiwUXUz7GgUPNM7RRtRxzzi6dk7A0N+5oYkv394uKZ9bOstjdzDhN7Ga2MP/LBiWlt1EldVESarH6KC0Vo5wtdmaJNChIzRxwYaSblYkJN1yQa8Z3HYSbG29D77odBu3wMYA6nMMFXEEIN3AHD9CBLghI4BXevYn35n2suqp569LO4I+8zx84xIo4 sha1_base64="rQ20K2/kGjcawcQHYn398eupReg=">AAACyHicjVJba9RAFJ7ES+tadfXVl4NFqNguiS/6UhB88XELblvY2S6TyUkydHJh5qSwDEERf6Fv/gV/hZM0Bd0qeGDgm+87tzlnkkYrS1H0Iwjv3L13f2f3weTh3qPHT6ZP905t3RqJC1nr2pwnwqJWFS5IkcbzxqAoE41nyeWHXj+7QmNVXX2iTYOrUuSVypQU5Kn19CcvBLl5t3Zc6KYQ3YU7ijuAgZb+xlOR52h6h6ZMnOoOgSdIooPtSDgGnqImccMcjo544RQM0Scdl2lNMKbydW4KbaeHXnsNXGNGBzEcwT9yc6Pygl79T7/r6X40iwaD2yAewT4bbb6efudpLdsSK5JaWLuMo4ZWThhSUmM34a3FRshLkePSw0qUaFdu2EkHLz2TQlYbfyqCgf09wonS2k2ZeM9SUGG3tZ78m7ZsKXu3cqpqWsJKXhfKWg1UQ79gSJVBSXrjgZBG+V5BFsIISf4bTPwQ4u0n3wanb2ZxNItPIrbLnrMX7IDF7C17zz6yOVswGcyDq+Bz8CVchC78ej2uMBjn9oz9YeG3X3FX4kI= sha1_base64="NiaOTykuxlLW3Jby+kQLz42VB8Y=">AAAC03icjVJLaxRBEO4ZX3F9ZNWjl8JFiJgsM7noRQh48bgBNwlsb5aanp7ZJj0PumsCSzMo4tU/582/4K+wZzIB3ShY0PD1V/VVVVd1UmtlKYp+BOGt23fu3tu5P3rw8NHj3fGTpye2aoyQc1HpypwlaKVWpZyTIi3PaiOxSLQ8TS7ed/7TS2msqsqPtKnlssC8VJkSSJ5ajX/yNZKbtSvHUddrbM/dQdwC9LTwN55inkvTBdRF4lS7DzyRhC1sK+Ed8FRqwmtmfwiU505Brz5uuUgrgiGVr3NdaDs9dL7XwLXMaC+GA/hHbm5UvqZX/9PvajyJplFvcBPEA5iwwWar8XeeVqIpZElCo7WLOKpp6dCQElq2I95YWaO4wFwuPCyxkHbp+p208NIzKWSV8ack6NnfFQ4LazdF4iMLpLXd9nXk33yLhrK3S6fKuiFZiqtCWaOBKugWDKkyUpDeeIDCKN8riDUaFOS/wcgPId5+8k1wcjiNo2l8HE2ODodx7LDn7AXbYzF7w47YBzZjcyaCWXAZfAo+h/PQhV/Cr1ehYTBonrE/LPz2C0Fz420= ˆ ˆ ˆ 1 ˆ ˆ P↵, H =0 P S iP = ⌘ˆi AAACOnicbVC7SgNBFJ31bXxFLW0Gg2Bj2BVBGyFgYxnRRCG7Lncnk2TI7IOZu0pY9rts/Ao7CxsLRWz9ACebLdR4YOBwzrncuSdIpNBo28/WzOzc/MLi0nJlZXVtfaO6udXWcaoYb7FYxuomAM2liHgLBUp+kygOYSD5dTA8G/vXd1xpEUdXOEq4F0I/Ej3BAI3kVy/cAWDWzP3MTRNQKr7Pb7MDJ6eF7iZhkF3muS8onQ7SU/oj5nKEcdKv1uy6XYBOE6ckNVKi6Vef3G7M0pBHyCRo3XHsBL0MFAomeV5xU80TYEPo846hEYRce1lxek73jNKlvViZFyEt1J8TGYRaj8LAJEPAgf7rjcX/vE6KvRMvE1GSIo/YZFEvlRRjOu6RdoXiDOXIEGBKmL9SNgAFDE3bFVOC8/fkadI+rDt23bk4qjUOyzqWyA7ZJfvEIcekQc5Jk7QIIw/khbyRd+vRerU+rM9JdMYqZ7bJL1hf32CLr38= AAACF3icbVBNS8NAEN34bf2qevSyWAQPUhIR9CIUvHisYFuhCWGy3TSLmw92J0IJ+Rde/CtePCjiVW/+G7dtDtr6YODx3gwz84JMCo22/W0tLC4tr6yurdc2Nre2d+q7e12d5orxDktlqu4C0FyKhHdQoOR3meIQB5L3gvursd974EqLNLnFUca9GIaJCAUDNJJfb1JX8hD71I0Ai3bpFy7ILILyZKpcl9RVYhihR+kltf16w27aE9B54lSkQSq0/fqXO0hZHvMEmQSt+46doVeAQsEkL2turnkG7B6GvG9oAjHXXjH5q6RHRhnQMFWmEqQT9fdEAbHWozgwnTFgpGe9sfif188xvPAKkWQ58oRNF4W5pJjScUh0IBRnKEeGAFPC3EpZBAoYmihrJgRn9uV50j1tOnbTuTlrtE6rONbIATkkx8Qh56RFrkmbdAgjj+SZvJI368l6sd6tj2nrglXN7JM/sD5/ACE0now= " " h i
Time reversal symmetry (Survives finite doping and and hence allows simulations at finite chemical potential)
cˆ cˆ ˆ 1 ˆ with ˆ2 T ↵ " T = ↵ # AAAB+nicbVBNS8NAEN3Ur1q/Uj16WSyCF0tSBL0IBS8eK/QL2lg22027dLMJuxOlxPwULx4U8eov8ea/cdvmoK0PBh7vzTAzz48F1+A431ZhbX1jc6u4XdrZ3ds/sMuHbR0lirIWjUSkuj7RTHDJWsBBsG6sGAl9wTr+5Gbmdx6Y0jySTZjGzAvJSPKAUwJGGtjl/phA2szu01qGr/G5iwd2xak6c+BV4uakgnI0BvZXfxjRJGQSqCBa91wnBi8lCjgVLCv1E81iQidkxHqGShIy7aXz0zN8apQhDiJlSgKeq78nUhJqPQ190xkSGOtlbyb+5/USCK68lMs4ASbpYlGQCAwRnuWAh1wxCmJqCKGKm1sxHRNFKJi0SiYEd/nlVdKuVV2n6t5dVOq1PI4iOkYn6Ay56BLV0S1qoBai6BE9o1f0Zj1ZL9a79bFoLVj5zBH6A+vzB1gukq0= cˆ cˆ T = 1
AAACtnicbVFNTxsxEPVuS0u3pQQ49mI1qtQL0S4gFamqhMSFI5USoMqGaNaZJBZe27K9pZG1P5FLb/wbvJscIGFkS09vPt7zuNCCW5emj1H85u3Wu/fbH5KPn3Y+73b29q+sqgzDAVNCmZsCLAouceC4E3ijDUJZCLwu7s6b/PVfNJYr2XcLjaMSZpJPOQMXqHHnIZ+D8/361h9mNaU5CD0HmuQFzrj0ugRn+L86aatYPfZ5pcEYdV/TPH/GTtS9XPJJjnKy1tgPk3+FoSo4aYz6pUy9qXP42shW6idtzoaNl3LjTjftpW3QTZCtQJes4nLc+R90WFWidEyAtcMs1W7kwTjOBIbplUUN7A5mOAxQQol25Nu11/RbYCZ0qky40tGWfd7hobR2URahMvib2/VcQ76WG1ZuejryXOrKoWRLoWklqFO0+UM64QaZE4sAgBkevFI2BwPMhZ9OwhKy9SdvgqujXpb2st8n3bOj1Tq2yRfylXwnGflBzsgFuSQDwqLj6E9URCw+jW9jjGfL0jha9RyQFxHrJ+Tu2co= ✓ #◆ ✓ " ◆ Fermion-Boson problems: The Su-Schrieffer-Heeger (SSH) model
Symmetries ˆ2 P i,j k ˆ ˆ h i ˆ 2 H = t cˆi,† cˆj, +ˆcj,† cˆi, + g X i,j cˆi,† cˆj, +ˆcj,† cˆi, + + X i,j h i 2m 2 h i i,j , i,j , i,j " # AAADz3icvVNNa9swGFbsfXTZV7oddxELg452wQ6D9VIo9NLDDiksbSByjazIjhrJNpLcEjSVXff3dtsP2P+YYiclawJjl70gePS8z/uJlJScKR0EP1ue/+Dho8c7T9pPnz1/8bKz++pcFZUkdEgKXshRghXlLKdDzTSno1JSLBJOL5LZycJ/cU2lYkX+Rc9LGgmc5SxlBGtHxbutX2iKtTm18Ah+0EhVIjaI01QjjvOMU8gOrpBk2VQj2RAHSLFMYFur9upoYmPDVvylQROcZVTaO9/Vyre/SW2R36VqCr+H+zD7l85gnWlk/6b/vyOsJlhrpoZNF2OUSkxMnWBgt+ouTd9a0xeuhVo7cxe7NurWgKZ+FHe6QS+oDW6CcAm6YGmDuPMDTQpSCZprwrFS4zAodWSw1Iy41G1UKVpiMsMZHTuYY0FVZOr3aOE7x0xgWkh3cg1rdj3CYKHUXCROKbCeqvu+BbnNN650ehgZlpeVpjlpCqUVh7qAi8cNJ0xSovncAUwkc71CMsVuWdp9gbZbQnh/5E1w3u+FQS88+9g97i/XsQPegLdgD4TgEzgGp2AAhoB4nz3pGe+rf+bf+Lf+t0bqtZYxr8Ef5n//DTqQSdo= hXi ⇣ ⌘ hXi ⇣ ⌘ hXi
U(N) symmetry is manifest cˆi,1 cˆi,1 . . . U . ,U†U =1 0 . 1 ! 0 . 1 cˆi,N cˆi,N BAAACnXicrVFda9RAFJ3EqjV+rfZNHzq4CD6UJSmChSIsiKAipQXTLWzW5WZyNzt0MhNmbqpLyL/qL/HNf+PsdqFfPvYyA4dz7pl75968VtJRHP8Nwnsb9x883HwUPX7y9Nnz3ouXx840VmAqjDL2JAeHSmpMSZLCk9oiVLnCUX76aamPztA6afQPWtQ4qaDUciYFkKemvfMsx1Lqtq6ArPzdRdkcqBXdtJU7ScezLMrOCkNuhS6lA5+Iurhis7KcE1hrfvEo5Xfz6g7P9pcnSn+2WQFlibZLP/Jk2uvHg3gV/DZI1qDP1nE47f3JCiOaCjUJBc6Nk7imSQuWpFDoqzYOaxCnUOLYQw0Vukm7mm7H33qm4DNj/dXEV+xVRwuVc4sq95m+77m7qS3J/2njhmZ7k1bquiHU4qLQrFGcDF+uihfSoiC18ACElb5XLuZgQZBfaOSHkNz88m1wvDtI4kFy9L4/3F2PY5O9Zm/YO5awD2zIvrBDljIRvAqGwdfgW7gdfg6/hwcXqWGw9myxaxGO/gGSesw0 sha1_base64="hP+6LrUf2d3tZaldqaQQvEKMXyw=">AAAB2XicbZDNSgMxFIXv1L86Vq1rN8EiuCozbnQpuHFZwbZCO5RM5k4bmskMyR2hDH0BF25EfC93vo3pz0JbDwQ+zknIvSculLQUBN9ebWd3b/+gfugfNfzjk9Nmo2fz0gjsilzl5jnmFpXU2CVJCp8LgzyLFfbj6f0i77+gsTLXTzQrMMr4WMtUCk7O6oyaraAdLMW2IVxDC9YaNb+GSS7KDDUJxa0dhEFBUcUNSaFw7g9LiwUXUz7GgUPNM7RRtRxzzi6dk7A0N+5oYkv394uKZ9bOstjdzDhN7Ga2MP/LBiWlt1EldVESarH6KC0Vo5wtdmaJNChIzRxwYaSblYkJN1yQa8Z3HYSbG29D77odBu3wMYA6nMMFXEEIN3AHD9CBLghI4BXevYn35n2suqp569LO4I+8zx84xIo4 sha1_base64="j/guVYSPPNYdE2vfGJxZ4L/h430=">AAACknicrVFRS9xAEN6kVm209uyjfeiiFPogR9KXFkQQpNBKEQWjwuV6TDZzucXNbtidaI+Qf9Vf0jf/jXvnQa320WEXPr5vZmf2m7xW0lEc3wbhi6WXyyurr6K19dcbb3qb6+fONFZgKowy9jIHh0pqTEmSwsvaIlS5wov86nCmX1yjddLoM5rWOKyg1HIsBZCnRr3fWY6l1G1dAVn5q4uyCVArulErd5OOZ1mUXReG3Bz9lY59IuriQZmV5YTAWnPDo5Q/z6u7PNubnSj92WYFlCXaLt3nyai3E/fjefCnIFmAHbaIk1HvT1YY0VSoSShwbpDENQ1bsCSFQt+1cViDuIISBx5qqNAN27m7Hf/gmYKPjfVXE5+zDytaqJybVrnP9HNP3GNtRv5PGzQ0/jJspa4bQi3uG40bxcnw2ap4IS0KUlMPQFjpZ+ViAhYE+YVG3oTk8ZefgvNP/STuJ6cxW2Xv2Db7yBL2mR2wb+yEpUwEW8FB8D04Ct+HX8Mf93aFwcK3t+yfCI/vALcYyzs= sha1_base64="du6ngX7CvuMdrE6W3SsrjC/wsjQ=">AAACnXicrVHfa9RAEN7Eamv80dO+6UMXD8GHciR9aaEIByKoSGmh6RUu5zHZzOWWbnbD7qR6hPxX/iV9879x73pgbX102IWP75tvZ3Ymr5V0FMe/gvDBxsNHm1uPoydPnz3f7r14ee5MYwWmwihjL3JwqKTGlCQpvKgtQpUrHOWXH5b66Aqtk0af0aLGSQWlljMpgDw17f3MciylbusKyMofXZTNgVrRTVu5l3Q8y6LsqjDkVuiPdOwTURe3bFaWcwJrzXcepfz/vLrHs6PlidJvbVZAWaLt0vc8mfb68SBeBb8PkjXos3WcTHvXWWFEU6EmocC5cRLXNGnBkhQKfdXGYQ3iEkoce6ihQjdpV9Pt+FvPFHxmrL+a+Iq97Wihcm5R5T7T9z13d7Ul+S9t3NDscNJKXTeEWtwUmjWKk+HLVfFCWhSkFh6AsNL3ysUcLAjyC438EJK7X74PzvcHSTxITuP+cH89ji32mr1h71jCDtiQfWInLGUieBUMg8/Bl3A3/Bh+DY9vUsNg7dlhf0U4+g2ROsww C B C @ A @ A Majorana fermions
i A : ii,1, = cˆii, +ˆcii,† , ii,2, = i cˆii, cˆii,† AAADt3icjVHbbhMxEJ3tcmnDLYVHXiwipCLSaLdCooCQylW8IBVE2kp1ibwbx7GyN7xORVnlAf4SvoFP4IGz7qaCRih45Z3xmTNnxp6oSHRpg+C7t+JfuHjp8upa68rVa9dvtNdv7pX51MSyH+dJbg4iUcpEZ7JvtU3kQWGkSKNE7keTF3V8/1iaUufZB3tSyKNUqEyPdCwsoEH7Jy/SqNIzxnXG2DP2mDHGlUhTMaiaUDfs8lKrVMzYU8afMJ7Ikd1gfCxsFc9OaUyfce7PIx8rPhRKSbPA4Uarsb3XhZr7ztXb+qPepl5Wb/O/67UG7U7QC9xii07YOB1q1m7e/kGchpRTTFNKSVJGFn5Cgkp8hxRSQAWwI6qAGXjaxSXNqIXcKVgSDAF0gr/C6bBBM5xrzdJlx6iSYBtkMrqL/dopRmDXVSX8EvYX9heHqX9WqJxy3eEJbATFNaf4FrilMRjLMtOGOe9leWZ9K0sj2na30eivcEh9z/hM5yUiBtjERRi9ckwFjcidj/ECGWwfHdSvPFdg7sZDWOGsdCpZoyigZ2Dr10c/GHN4fqiLzt5WLwx64bsHnZ3nzcBX6TbdoQ1M9SHt0BvaRR+x99777H31vvmP/IE/8sen1BWvyblFfy3/02+bv+Nq 2 ˆ , ˆ =2
AAADXnichVFNa9tAEB1ZbZq4TeOkl0IPFTWlKRgjmUJ6KQT6QS+FFOokkA1mJa/ljVcfSOtAIvQr+uvaU/9GIYe83UiGJpSsWM3smzdvZ3bCXMlS+/4vp+M+eLj2aH2j+/jJ5tOt3vbOYZkti0iMo0xlxXHIS6FkKsZaaiWO80LwJFTiKFx8NPGjc1GUMkt/6ItcnCY8TuVMRlwDmvR+MiVmmlUem3NdsZgnCa8n1dmAcZXP+YCVMgYyuBWXAxYK3Ybf1J7nsULGc83gfvBGOE6F0twwz+rWbzVNar1i3GispLrdSa/vD327vLtO0Dh9atZB1vtNjKaUUURLSkhQShq+Ik4lvhMKyKcc2ClVwAp40sYF1dRF7hIsAQYHusA/xumkQVOcjWZpsyPcorALZHr0GvuLVQzBNrcK+CXsFfalxeL/3lBZZVPhBWwIxQ2r+A24pjkY92UmDbOt5f5M05WmGb233UjUl1vE9BmtdD4hUgBb2IhHny0zhkZoz+d4gRR2jArMK7cKnu14CsutFVYlbRQ59ApY8/qoB2MObg/1rnM4Ggb+MPj+rr8/aga+Ti/oFe1iqnu0T1/pAHVE9Nd56ew6bzt/3DV30926oXacJucZ/bPc59eVwMMC j,↵, i, , i,j ↵, , ⇣ ⌘ ⇣ ⌘ { 0 } 0 i B : ii,2, = cˆii, +ˆcii,† , ii,1, = i cˆii, cˆii,† AAADs3icjVHbbhMxEJ3tcmnDLcAjLxYRUhFptFuBQCCkqlzEC6hIpK0Ul8jeOBsre5PXqVRWy3/CF/ARPHDW2VbQqCpeeWd85syZsUcWiS5tEPzw1vwrV69dX9/o3Lh56/ad7t17+2W+MJEaRnmSm0MpSpXoTA2ttok6LIwSqUzUgZy/aeIHx8qUOs++2JNCHaUizvRUR8ICGnd/8SKVla4Z1xljuy8ZYzwWaSrGVRvpb/d5qeNU1Ow12+KJmtpNPhO2iuolh+kzwhPWRr5WfCLiWJkVDjc6ntnHfcZfue9csfDvYpoty7GL6m39d73OuNsLBoFbbNUJW6dH7drLuz+J04RyimhBKSnKyMJPSFCJb0QhBVQAO6IKmIGnXVxRTR3kLsBSYAigc/xjnEYtmuHcaJYuO0KVBNsgk9Ej7PdOUYLdVFXwS9jf2N8cFl9YoXLKTYcnsBKKG07xI3BLMzAuy0xb5mkvl2c2t7I0pRfuNhr9FQ5p7hmd6bxFxACbuwijd44ZQ0O68zFeIIMdooPmlU8VmLvxBFY4q5xK1ioK6BnY5vXRD8Ycnh/qqrO/PQiDQfj5aW9ntx34Oj2gh7SJqT6nHfpAe+gj8j551qu97/4zf+RLf7Kkrnltzn36Z/npH07b4nk= 2 N,⇣2 ⌘ ⇣ ⌘Pˆ2 ˆ 1 ˆ i,j k ˆ 2 H = t + gX i,j i ˆi,↵, ˆj,↵, + h i + X i,j 2 h i 2m 2 h i i,j ,↵=1, i,j " # AAAEEHichVFNb9NAEB3XfLTho2k5clkRIRU1RHaEBJdKlSioF1CQSBspTqK1u3G28ZfsTaXW8p/gb/AHuCGu3DhS/goHntcuEmlL11rP7MybN2923SSQmbKsn8aKeev2nbura4179x88XG9ubB5k8SL1RN+LgzgduDwTgYxEX0kViEGSCh66gTh056/L/OGJSDMZRx/VaSJGIfcjOZUeVwhNNow9Z8ZVvl+wHeZki3CSOwGP/EAw2T52Uu22mcODZMZ37LaTST/krBjn79vdgjFnmnIvt4scBycQU7XFniu2zXymaQfFVXyAptKfqWdMVjDH52HIgZXtqlPdp1hKHy+lxzkkoNu1ygutaVip1Fy9qxWNMQCGCIvtCjsvJ/rfCGVBNcWINRqTZsvqWHqxy45dOy2qVy9unpNDRxSTRwsKSVBECn5AnDJ8Q7LJogSxEeWIpfCkzgsqqIHaBVACCI7oHH8fp2EdjXAuOTNd7aFLgJ2iktFT7Lea0QW67CrgZ7C/sc90zL+2Q66ZS4WnsC4Y1zTjO8QVzYC4qTKskRdabq4sp1I0pVd6Ggl9iY6Uc3p/efaQSRGb6wyjNxrpg8PV5xPcQATbh4Lyli8YmJ74CJZrKzRLVDNy8KWw5e1DD57ZXn7Uy85Bt2NbHfvDi9Zut37wVXpMT2gLr/qSdmmfetDhGZ+NH8a58cv8ZH4xv5rfKuiKUdc8on+W+f0Pg+EI1A== h iX ⇣ ⌘ hXi ˆi,1,1 ˆi,1,1 . . O(2N) symmetry is now manifest 0 . 1 0 . 1 ˆ ˆ B i,i,N C O B i,i,N C ,OT O =1 B ˆi,2,1 C ! B ˆi,2,1 C B C B C B . C B . C B . C B . C B C B C B ˆi,2,N C B ˆi,2,N C BAAADKniczVJNaxsxENVuv9Ltl9MeexE1hR6MWZlCCyVg6KWnOoU4CXhdM6uV1yL6WKTZpGbZ39NL/0ovObSEXvtDKjs+pE6g1w4SPN7M441Gk1dKekzTiyi+dfvO3Xs795MHDx89ftLZfXrobe24GHOrrDvOwQsljRijRCWOKydA50oc5SfvV/mjU+G8tOYAl5WYaiiNnEsOGKjZbjTMclFK01Qa0MkvbZItAJusBK2hnTWyx3qspVmWZKeFRb9GWxWy97G9iR/8UzkIyiQTprhi72S5QHDOntFkRJP/rL0ezd6tTjL63By0oz3KZp1u2k/XQa8DtgFdson9Wec8KyyvtTDIFXg/YWmF0wYcSq5EMKy9qICfQCkmARrQwk+b9Ve39GVgCjq3LlyDdM1eVTSgvV/qPFSGlhd+O7cib8pNapy/nTbSVDUKwy+N5rWiaOlqb2ghneColgEAdzL0SvkCHHAM25WEIbDtJ18Hh4M+S/vs0+vucLAZxw55Tl6QV4SRN2RIPpB9MiY8+hp9j35EP+Nv8Xl8Ef+6LI2jjeYZ+Svi338AFFsAdA== C B C @ A @ A Fermion-Boson problems: Unconstrained gauge theories
ˆ2 P i,j k ˆ ˆ h i ˆ 2 H = t cˆi,† cˆj, +ˆcj,† cˆi, + g X i,j cˆi,† cˆj, +ˆcj,† cˆi, + + X i,j h i 2m 2 h i i,j , i,j , i,j " # AAADz3icvVNNa9swGFbsfXTZV7oddxELg452wQ6D9VIo9NLDDiksbSByjazIjhrJNpLcEjSVXff3dtsP2P+YYiclawJjl70gePS8z/uJlJScKR0EP1ue/+Dho8c7T9pPnz1/8bKz++pcFZUkdEgKXshRghXlLKdDzTSno1JSLBJOL5LZycJ/cU2lYkX+Rc9LGgmc5SxlBGtHxbutX2iKtTm18Ah+0EhVIjaI01QjjvOMU8gOrpBk2VQj2RAHSLFMYFur9upoYmPDVvylQROcZVTaO9/Vyre/SW2R36VqCr+H+zD7l85gnWlk/6b/vyOsJlhrpoZNF2OUSkxMnWBgt+ouTd9a0xeuhVo7cxe7NurWgKZ+FHe6QS+oDW6CcAm6YGmDuPMDTQpSCZprwrFS4zAodWSw1Iy41G1UKVpiMsMZHTuYY0FVZOr3aOE7x0xgWkh3cg1rdj3CYKHUXCROKbCeqvu+BbnNN650ehgZlpeVpjlpCqUVh7qAi8cNJ0xSovncAUwkc71CMsVuWdp9gbZbQnh/5E1w3u+FQS88+9g97i/XsQPegLdgD4TgEzgGp2AAhoB4nz3pGe+rf+bf+Lf+t0bqtZYxr8Ef5n//DTqQSdo= hXi ⇣ ⌘ hXi ⇣ ⌘ hXi
ˆ ˆ k !mX i,j iP i,j ! = ˆ† h i h i 0 Bosonic variables. b = AAAC6HichVFNa9tAEH1W08Z2P+I2x0IRNYGcjBwKTQ4BQ5KSS8GFOAnEwayUtSr0mdU64Ard+gd6Cr3mmGvyZ9rf0kOfNnKhDSUrVjP75s3bmR03i4JcO86PhvVo6fGT5War/fTZ8xcrnZevDvN0pjw58tIoVceuyGUUJHKkAx3J40xJEbuRPHLDnSp+dCFVHqTJgZ5n8jQWfhJMA09oQpPOm3EaS19MCqfctsf5udLFeKqEV4RlEZflpNN1eo5Z9n2nXztd1GuYdn5ijDOk8DBDDIkEmn4EgZzfCfpwkBE7RUFM0QtMXKJEm7kzsiQZgmjIv8/TSY0mPFeaucn2eEvErZhpY437g1F0ya5ulfRz2l/cXwzm//eGwihXFc5pXSq2jOJH4hqfyXgoM66Zi1oezqy60phi03QTsL7MIFWf3h+dXUYUsdBEbOwZpk8N15wv+AIJ7YgVVK+8ULBNx2e0wlhpVJJaUVBP0Vavz3o45v6/Q73vjDZ6W73+p3fdwUY97yZe4y3WOdT3GGAfQ5bh4SuucYNbK7S+WZfW9zuq1ahzVvHXsq5+A/dpm88= i,j p rm AAACf3icdVFdixMxFM2MH7vWr6qPvly2iIJrmSnC6oOw4IuPFexuoanlTnpnGjfJjElGKGH+hj/MN/+LD6YfoO66FwKHc8+5Sc4tGiWdz7KfSXrj5q3bB4d3enfv3X/wsP/o8ZmrWytoImpV22mBjpQ0NPHSK5o2llAXis6Li/eb/vk3sk7W5pNfNzTXWBlZSoE+Uov+d75CH4ruc+BLrCqy3SJwhaZSBPL4C3C7wx3AO+ClRRGA15oqBA1b7/RaxysAudOMr9N0ELj7an0YwZ+xHXSL/iAbZtuCqyDfgwHb13jR/8GXtWg1GS8UOjfLs8bPA1ovhaKux1tHDYoLrGgWoUFNbh62+XXwLDJLKGsbj/GwZf92BNTOrXURlRr9yl3ubcj/9WatL9/MgzRN68mI3UVlq8DXsFkGLKUl4dU6AhRWxreCWGEM2ceV9WII+eUvXwVno2GeDfOPrweno30ch+wpO2IvWM5O2Cn7wMZswgT7lRwlL5PjNEmfp8M020nTZO95wv6p9O1vlY7AFQ== h i 2!m
ˆ 1 ˆ ˆ ˆ ˆ H = g b†i,j + b i,j cˆi,† cˆj, +ˆcj,† cˆi, + !0 b†i,j b i,j 2!0m h i h i h i h i i,j b AAADkHiclVJdb9MwFHUTPkZh0I1HXiyqSUOdqqQgbQ9MDPEy8TQkuk2qS+S4TuLNdjLbQaqs/B1+EG/8G5wP1rEOAVeydHTOvefea9244EybIPjR8/x79x883HjUf/xk8+mzwdb2qc5LReiU5DxX5zHWlDNJp4YZTs8LRbGIOT2LLz/U+tlXqjTL5WezLOhc4FSyhBFsHBVt9b6hDBt7XEEIDyHSpYgs4limnEK2dwGRarHTUydfKWNRojCxYWUnEOWCpjgKoIBV7YA4TcwubCzj6otFC5ymVFV/9BzdyC0yLE0u/lqDFEsz86pu146rWSpw1fZu7IgrZnsdv5riWrv4pY3WqTvSr61WnesYrdZvB4mrf9/8f9eOBsNgHDQB10HYgSHo4iQafEeLnJSCSkM41noWBoWZW6wMI5xWfVRqWmByiVM6c1BiQfXcNgdVwR3HLGCSK/ekgQ17s8JiofVSuIV3BDaZvq3V5F3arDTJwdwyWZSGStI2SkoOTQ7r64QLpigxfOkAJoq5WSHJsLs342647z4hvL3yOjidjMNgHH56MzyadN+xAV6Al2AXhGAfHIFjcAKmgHib3mvvrXfob/sH/jv/fZvq9bqa5+C38D/+BAGKJmU= sha1_base64="hP+6LrUf2d3tZaldqaQQvEKMXyw=">AAAB2XicbZDNSgMxFIXv1L86Vq1rN8EiuCozbnQpuHFZwbZCO5RM5k4bmskMyR2hDH0BF25EfC93vo3pz0JbDwQ+zknIvSculLQUBN9ebWd3b/+gfugfNfzjk9Nmo2fz0gjsilzl5jnmFpXU2CVJCp8LgzyLFfbj6f0i77+gsTLXTzQrMMr4WMtUCk7O6oyaraAdLMW2IVxDC9YaNb+GSS7KDDUJxa0dhEFBUcUNSaFw7g9LiwUXUz7GgUPNM7RRtRxzzi6dk7A0N+5oYkv394uKZ9bOstjdzDhN7Ga2MP/LBiWlt1EldVESarH6KC0Vo5wtdmaJNChIzRxwYaSblYkJN1yQa8Z3HYSbG29D77odBu3wMYA6nMMFXEEIN3AHD9CBLghI4BXevYn35n2suqp569LO4I+8zx84xIo4 sha1_base64="gurAvKVpNEBy9ziNqFumI0QMxv8=">AAADhXiclVJdb9MwFHUTPkbZoNsrLxbTpKFOVbI9sAcmQLzscUh0m1SX6MZ1Um+2E2wHqbLyd/hBvPFvcJKyjnUIuJKlo3N8z/3QTUvBjY2iH70gfPDw0eONJ/2nm1vPng+2N89NUWnKxrQQhb5MwTDBFRtbbgW7LDUDmQp2kV5/aPSLr0wbXqhPdlGyqYRc8YxTsJ5KtnvfyBysO60xxieYmEomjghQuWCYH1xhojvs9dzLX7R1JNNAXVy7Q0wKyXJIIixx3TgQwTK7j1vLtP7syAzynOn6j57DW3/LOShbyL/mEM3zuX3VlOvaNTyXUHe1Wzvqk/nBkl91caNd/dKG69Q932+sVpWbGK7G7xpJ63+f/H/HTga70ShqA6+DeAl20TLOksF3MitoJZmyVIAxkzgq7dSBtpwKVvdJZVgJ9BpyNvFQgWRm6tqDqvGeZ2Y4K7R/yuKWvZ3hQBqzkH7gPQl2bu5qDXmfNqlsdjx1XJWVZYp2hbJKYFvg5jrxjGtGrVh4AFRz3yumc/D3Zv0N9/0S4rsjr4Pzw1EcjeKPEdpAL9BLtI9i9Bq9Q6foDI0RDbaCo+BNcBLuhMfh225dQW+5tx30W4TvfwI+UyV4 sha1_base64="4Yswar4j5FyT7VQucq24bhmK+jY=">AAADkHiclVJNb9QwEPUmQMtC6bYcuVisKhVttUqWAz1QUcSl4lQktq20XiLH6yRubSe1HaSVlb/DD+LGv8H5oFu6RcBIlp7em3kzY01ccKZNEPzoef6Dh482Nh/3nzzderY92Nk903mpCJ2SnOfqIsaacibp1DDD6UWhKBYxp+fx1YdaP/9KlWa5/GyWBZ0LnEqWMIKNo6Kd3jeUYWNPKgjhEUS6FJFFHMuUU8gOLiFSLXZ66uRrZSxKFCY2rOwEolzQFEcBFLCqHRCnidmHjWVcfbFogdOUquqPnqNbuUWGpcnFX2uQYmlmXtXt2nE1SwWu2t6NHXHF7KDjV1PcaJe/tNE6dU/6jdWqcx2j1frtIHH175v/79rRYBiMgybgOgg7MARdnEaD72iRk1JQaQjHWs/CoDBzi5VhhNOqj0pNC0yucEpnDkosqJ7b5qAquOeYBUxy5Z40sGFvV1gstF4Kt/CewCbTd7WavE+blSY5nFsmi9JQSdpGScmhyWF9nXDBFCWGLx3ARDE3KyQZdvdm3A333SeEd1deB2eTcRiMw0/B8HjSfccmeAFegn0QgjfgGJyAUzAFxNvyXntvvSN/1z/03/nv21Sv19U8B7+F//EnAEomYQ== r hXi ⇣ ⌘ X ⇣ ⌘ X
Local Z2 conservation law
nˆb +ˆnb +ˆnb +ˆnb c i,i+a i,i a i,i+a i,i a nˆi Qˆ =( 1) h xi h xi h y i h y i ( 1) AAAC13ichVJNSyNBEO0Zd103rhrd416KDYKihhlZ0IsQ8LLHCMYPkjjUdDpJk56eobtHHIfBg8vi1b+2t/0T/gY7ySxk/SxoeLz3qrq6qsNEcG0876/jzn34OP9p4XNl8cvS8kp1de1Ex6mirEVjEauzEDUTXLKW4Uaws0QxjELBTsPR4Vg/vWRK81gemyxh3QgHkvc5RWOpoPoA0BmiyY+KgAMcQBkbO/7mRV5qsrjIwyLIOwLlQDDg2xy2AIMrq6spVYBlXnfvTNz/zNZtzW/Xzmbd79bOZhopXuie2hxuhSKo1ry6Nwl4DvwS1EgZzaD6p9OLaRoxaahArdu+l5hujspwKlhR6aSaJUhHOGBtCyVGTHfzyV4KWLdMD/qxskcamLCzGTlGWmdRaJ0RmqF+qo3Jl7R2avr73ZzLJDVM0ulF/VSAiWG8ZOhxxagRmQVIFbe9Ah2iQmrsV6jYIfhPn/wcnOzWfa/uH/2oNXbLcSyQb+Q72SA+2SMN8pM0SYtQ59i5dm6dX+65e+P+du+mVtcpc76S/8K9fwRxYdss i ˆ ˆ H,Qi =0 AAAC8nichVFNSxxBEH2OJn7kw9UcvTQuITmEZUYEvQhCTPAiKLgq7C5Lz9g72+zszNDTu6DL/oT8AW8h19y8mh8Sf4sH37RjIJFgDz31+lXV66quME90YX3/94w3O/fi5fzC4tKr12/eLtdWVk+KbGQi1YyyJDNnoSxUolPVtNom6iw3Sg7DRJ2Gg8+l/3SsTKGz9Nhe5KozlHGqezqSllS39qGdqJ5tiXZf2sn+9NMDOJp2tWgbHfdtR4gd4Yture43fLfEUxBUoI5qHWa1W7RxjgwRRhhCIYUlTiBR8GshgI+cXAcTcoZIO7/CFEvMHTFKMUKSHfAf89Sq2JTnUrNw2RFvSbgNMwXec391iiGjy1sVcUF7x33puPi/N0ycclnhBW1IxUWneEDeos+I5zKHVeRjLc9nll1Z9LDtutGsL3dM2Wf0R2ePHkNu4DwCX1xkTI3Qncd8gZS2yQrKV35UEK7jc1rprHIqaaUoqWdoy9dnPRxz8O9Qn4KTjUbgN4KjzfquXw18AWtYx0dOdQu72Mch64jwDde4wS/Pelfed+/HQ6g3U+W8w1/L+3kP3DGdYQ== h i Boson parity on star Fermion parity on site F. F. ASSAAD and TARUN GROVER PHYS. REV. X 6, 041049 (2016) and velocity fluctuations in Dirac metals [26]. It consists of N species of fermions on the vertices of a square lattice interacting with quantum Ising spins that live on the links (d) (e) of the same square lattice. However, in contrast to the Hamiltonian in Ref. [25], our Hamiltonian has a local Z2 (a) invariance, which leads to LxLy conserved operators on a Lx × Ly lattice. Furthermore, we restrict ourselves to half- filling; thus, at low energies, instead of a Fermi surface, we obtain either Fermi points or no fermions (e.g., a Mott insulator). These differences completely alter the phase diagram of our model compared to Ref. [25]. The phase diagram for various values of N is summa- rized in Fig. 1. For N>1, we find that when the kinetic energy of the aforementioned Ising spins is small compared to the kinetic energy of the fermions, the ground state resembles the deconfined phase of the Z2 gauge theory coupled to 2N Dirac fermions (the N 1 case is special in that the deconfined fermions are unstable¼ to infinitesimal interactions). However, as we will discuss in detail, there is an important distinction between this phase and the (b) (c) deconfined phase of a conventional Z2 gauge matter theory (Refs. [27–29]). Namely, the local Z2 invariance in our FIG. 1. Schematic zero-temperature phase diagram of the model is an actual symmetry of the Hamiltonian, in model in Eq. (1) (a), as well as cartoons (b)–(e) of a selected contrast to the so-called “gauge symmetry” in gauge-matter number of phases. (a) We observe a Z2 Dirac deconfined phase theories, which just corresponds to redundancy in the (Z2D), a Néel antiferromagnet phase (AFM) [or a superconductor description of the physical states. This is because we do (SC), depending on the pattern of particle-hole symmetry break- not project the Hilbert space of our Hamiltonian to gauge- ing], a charge density wave (CDW) phase, and a valence-bond solid (VBS). For N 1, we do not find evidence for a Z D phase invariant states; that is, we do not impose Gauss’s law. ¼ 2 Instead, in our model, Gauss’s law in the ground state beyond h 0, consistent with the arguments in the main text. The phase¼ transitions from the Z D to AFM/SC (N 2) and emerges because of spontaneous symmetry breaking as the 2 VBS (N 3) are seemingly continuous. At N 3, we observe¼ a locally conserved operators order in a certain definite deconfined¼ quantum critical point (DCQP) between¼ the VBS and pattern. One physical significance of this result is that, AFM phases. (b)–(e) Cartoons of the corresponding phases. in contrast to regular Z2 gauge theory, in ourFermion model,-Boson equal- problems:Circles correspondUnconstrained to fermions, gauge and theories the color code corresponds space, unequal-time, non-gauge-invariant correlation func- to the flavor index. Circles with two colors represent a pair of tions do not vanish identically. Furthermore, there is a fermions on a site with corresponding flavors. The low-energy finite-temperature Ising transitionˆ corresponding1 ˆ to theˆ properties of the Z2D phase resembleˆ SUˆN fermions propa- H = g b†i,j + b i,j gating freelycˆi,† cˆj, in +ˆ space-timecj,† cˆi, and+ ! connected0 b†i,j byðb i,jÞ a Z gauge string restoration of the aforementioned symmetry.2!0m h i h i h i h i 2 i,j b AAADkHiclVJdb9MwFHUTPkZh0I1HXiyqSUOdqqQgbQ9MDPEy8TQkuk2qS+S4TuLNdjLbQaqs/B1+EG/8G5wP1rEOAVeydHTOvefea9244EybIPjR8/x79x883HjUf/xk8+mzwdb2qc5LReiU5DxX5zHWlDNJp4YZTs8LRbGIOT2LLz/U+tlXqjTL5WezLOhc4FSyhBFsHBVt9b6hDBt7XEEIDyHSpYgs4limnEK2dwGRarHTUydfKWNRojCxYWUnEOWCpjgKoIBV7YA4TcwubCzj6otFC5ymVFV/9BzdyC0yLE0u/lqDFEsz86pu146rWSpw1fZu7IgrZnsdv5riWrv4pY3WqTvSr61WnesYrdZvB4mrf9/8f9eOBsNgHDQB10HYgSHo4iQafEeLnJSCSkM41noWBoWZW6wMI5xWfVRqWmByiVM6c1BiQfXcNgdVwR3HLGCSK/ekgQ17s8JiofVSuIV3BDaZvq3V5F3arDTJwdwyWZSGStI2SkoOTQ7r64QLpigxfOkAJoq5WSHJsLs342647z4hvL3yOjidjMNgHH56MzyadN+xAV6Al2AXhGAfHIFjcAKmgHib3mvvrXfob/sH/jv/fZvq9bqa5+C38D/+BAGKJmU= sha1_base64="hP+6LrUf2d3tZaldqaQQvEKMXyw=">AAAB2XicbZDNSgMxFIXv1L86Vq1rN8EiuCozbnQpuHFZwbZCO5RM5k4bmskMyR2hDH0BF25EfC93vo3pz0JbDwQ+zknIvSculLQUBN9ebWd3b/+gfugfNfzjk9Nmo2fz0gjsilzl5jnmFpXU2CVJCp8LgzyLFfbj6f0i77+gsTLXTzQrMMr4WMtUCk7O6oyaraAdLMW2IVxDC9YaNb+GSS7KDDUJxa0dhEFBUcUNSaFw7g9LiwUXUz7GgUPNM7RRtRxzzi6dk7A0N+5oYkv394uKZ9bOstjdzDhN7Ga2MP/LBiWlt1EldVESarH6KC0Vo5wtdmaJNChIzRxwYaSblYkJN1yQa8Z3HYSbG29D77odBu3wMYA6nMMFXEEIN3AHD9CBLghI4BXevYn35n2suqp569LO4I+8zx84xIo4 sha1_base64="gurAvKVpNEBy9ziNqFumI0QMxv8=">AAADhXiclVJdb9MwFHUTPkbZoNsrLxbTpKFOVbI9sAcmQLzscUh0m1SX6MZ1Um+2E2wHqbLyd/hBvPFvcJKyjnUIuJKlo3N8z/3QTUvBjY2iH70gfPDw0eONJ/2nm1vPng+2N89NUWnKxrQQhb5MwTDBFRtbbgW7LDUDmQp2kV5/aPSLr0wbXqhPdlGyqYRc8YxTsJ5KtnvfyBysO60xxieYmEomjghQuWCYH1xhojvs9dzLX7R1JNNAXVy7Q0wKyXJIIixx3TgQwTK7j1vLtP7syAzynOn6j57DW3/LOShbyL/mEM3zuX3VlOvaNTyXUHe1Wzvqk/nBkl91caNd/dKG69Q932+sVpWbGK7G7xpJ63+f/H/HTga70ShqA6+DeAl20TLOksF3MitoJZmyVIAxkzgq7dSBtpwKVvdJZVgJ9BpyNvFQgWRm6tqDqvGeZ2Y4K7R/yuKWvZ3hQBqzkH7gPQl2bu5qDXmfNqlsdjx1XJWVZYp2hbJKYFvg5jrxjGtGrVh4AFRz3yumc/D3Zv0N9/0S4rsjr4Pzw1EcjeKPEdpAL9BLtI9i9Bq9Q6foDI0RDbaCo+BNcBLuhMfh225dQW+5tx30W4TvfwI+UyV4 sha1_base64="4Yswar4j5FyT7VQucq24bhmK+jY=">AAADkHiclVJNb9QwEPUmQMtC6bYcuVisKhVttUqWAz1QUcSl4lQktq20XiLH6yRubSe1HaSVlb/DD+LGv8H5oFu6RcBIlp7em3kzY01ccKZNEPzoef6Dh482Nh/3nzzderY92Nk903mpCJ2SnOfqIsaacibp1DDD6UWhKBYxp+fx1YdaP/9KlWa5/GyWBZ0LnEqWMIKNo6Kd3jeUYWNPKgjhEUS6FJFFHMuUU8gOLiFSLXZ66uRrZSxKFCY2rOwEolzQFEcBFLCqHRCnidmHjWVcfbFogdOUquqPnqNbuUWGpcnFX2uQYmlmXtXt2nE1SwWu2t6NHXHF7KDjV1PcaJe/tNE6dU/6jdWqcx2j1frtIHH175v/79rRYBiMgybgOgg7MARdnEaD72iRk1JQaQjHWs/CoDBzi5VhhNOqj0pNC0yucEpnDkosqJ7b5qAquOeYBUxy5Z40sGFvV1gstF4Kt/CewCbTd7WavE+blSY5nFsmi9JQSdpGScmhyWF9nXDBFCWGLx3ARDE3KyQZdvdm3A333SeEd1deB2eTcRiMw0/B8HjSfccmeAFegn0QgjfgGJyAUzAFxNvyXntvvSN/1z/03/nv21Sv19U8B7+F//EnAEomYQ== r As the kinetic energy of the IsinghXi degrees of⇣ freedom ⌘(b).X The⇣ symmetry-broken⌘ phasesX correspond to the confined is increased, we find that the aforementioned gapless phases of the model. At N 3, the AFM phase (e) has the fundamental (conjugate) representation¼ of SU(3) on sublattice Dirac deconfined phase enters a symmetry-broken confined A (B). The corresponding Young tableaux are included. The VBS phase. The nature of symmetry breaking depends on the phase (d) corresponds to a pattern of intersite SU(3) singlets. value of N. For N 1, it corresponds to a charge density Retain only one phonon excita/on per bond, float AAAC3nichVHJSgNBEK2Me9zicvPSGARPYUYEl5PgghdFwRjBiPRM2rHJbHR3AiqePYlXj171j/RbPPimnQgqYg09Vf2q6nUtfhZJbVz3teT09Q8MDg2PlEfHxicmK1PTxzrtqEDUgzRK1YnPtYhkIupGmkicZErw2I9Ew29v5v5GVygt0+TIXGXiLOZhIi9kwA2g88psU8sw5mydMY81g1ZqNNs/r1TdmmuF/Ta8wqhSIQdp5Y2a1KKUAupQTIISMrAj4qTxnZJHLmXAzugGmIIlrV/QLZWR20GUQAQH2sY/xO20QBPcc05tswO8EuEoZDJawNmxjD6i81cFbA39jnNtsfDPF24sc17hFbQPxhHLuAfc0CUi/suMi8heLf9n5l0ZuqBV241EfZlF8j6DL54teBSwtvUw2raRITh8e+9iAgl0HRXkU+4xMNtxC5pbLSxLUjBy8CnofPqoB2v2fi71t1Ffqq3VvMPl6sZSse9hmqN5WsRSV2iDdukAZQQY9hM904vDnTvn3nn4DHVKRc4MfRPn8QPpj5WD :1 N wave, while for N ¼ 2, because of an enlarged symmetry, ··· ¼ II. MODEL AND SYMMETRIES k the symmetry-broken phase can correspond to a super- !0 =
AAAC6HichVFNa9tAEH1W08Z2P+I2x0IRNYGcjBwKTQ4BQ5KSS8GFOAnEwayUtSr0mdU64Ard+gd6Cr3mmGvyZ9rf0kOfNnKhDSUrVjP75s3bmR03i4JcO86PhvVo6fGT5War/fTZ8xcrnZevDvN0pjw58tIoVceuyGUUJHKkAx3J40xJEbuRPHLDnSp+dCFVHqTJgZ5n8jQWfhJMA09oQpPOm3EaS19MCqfctsf5udLFeKqEV4RlEZflpNN1eo5Z9n2nXztd1GuYdn5ijDOk8DBDDIkEmn4EgZzfCfpwkBE7RUFM0QtMXKJEm7kzsiQZgmjIv8/TSY0mPFeaucn2eEvErZhpY437g1F0ya5ulfRz2l/cXwzm//eGwihXFc5pXSq2jOJH4hqfyXgoM66Zi1oezqy60phi03QTsL7MIFWf3h+dXUYUsdBEbOwZpk8N15wv+AIJ7YgVVK+8ULBNx2e0wlhpVJJaUVBP0Vavz3o45v6/Q73vjDZ6W73+p3fdwUY97yZe4y3WOdT3GGAfQ5bh4SuucYNbK7S+WZfW9zuq1ahzVvHXsq5+A/dpm88= m conductor or a Néel AFM, depending on the sign of the Our model reads r infinitesimal symmetry-breaking field. The N 3 case is even richer. We find that after exiting the deconfined¼ phase, N the model enters a VBS phase, and as the tuning parameter ˆ ˆ † ˆ H Z i;j cˆi;αcˆj;α H:c: − Nh X i;j : 1 controlling the kinetic energy of Ising spins is tuned further, ¼ h i þ h i ð Þ i;j α 1 i;j the valence-bond-solid quantum melts and yields a Néel Xh i X¼ Xh i AFM phase. The phase transition between the VBS and † k AFM phase, if second order, corresponds to a deconfined Here, cˆi;α creates a fermion of flavor α inh a Wannier!0 = state AAAC6HichVFNa9tAEH1W08Z2P+I2x0IRNYGcjBwKTQ4BQ5KSS8GFOAnEwayUtSr0mdU64Ard+gd6Cr3mmGvyZ9rf0kOfNnKhDSUrVjP75s3bmR03i4JcO86PhvVo6fGT5War/fTZ8xcrnZevDvN0pjw58tIoVceuyGUUJHKkAx3J40xJEbuRPHLDnSp+dCFVHqTJgZ5n8jQWfhJMA09oQpPOm3EaS19MCqfctsf5udLFeKqEV4RlEZflpNN1eo5Z9n2nXztd1GuYdn5ijDOk8DBDDIkEmn4EgZzfCfpwkBE7RUFM0QtMXKJEm7kzsiQZgmjIv8/TSY0mPFeaucn2eEvErZhpY437g1F0ya5ulfRz2l/cXwzm//eGwihXFc5pXSq2jOJH4hqfyXgoM66Zi1oezqy60phi03QTsL7MIFWf3h+dXUYUsdBEbOwZpk8N15wv+AIJ7YgVVK+8ULBNx2e0wlhpVJJaUVBP0Vavz3o45v6/Q73vjDZ6W73+p3fdwUY97yZe4y3WOdT3GGAfQ5bh4SuucYNbK7S+WZfW9zuq1ahzVvHXsq5+A/dpm88= rm quantum critical point [30], reported earlier in the context centered aroundPHYSICAL lattice REVIEW site Xi 6,of041049 a square (2016) lattice, and the sum of quantum magnets [31,32]. Our results are consistent runs over bonds of nearest neighbors. Each bond, b i; j , Simple Fermionic Model of Deconfined Phases and Phase Transitions ¼ h i with a second-order transition. accommodates a two-levelF. F. Assaad1 and Tarun system Grover2,3 spanned by the Hilbert 1Institut für Theoretische Physik und Astrophysik, Universität Würzburg, 97074 Würzburg, Germany 2Department of Physics, University of California at San Diego, La Jolla, California 92093, USA 3Kavli Institute for Theoretical Physics, University of California, Santa Barbara, California 93106, USA (Received 7 August 2016; revised manuscript received 2 November 2016; published 12 December 2016) Using quantum Monte Carlo simulations, we study a series of models of fermions coupled to quantum Ising spins on a square lattice with N flavors of fermions per site for N 1, 2, and 3. The models have an ¼ extensive number of conserved quantities but are not integrable, and they have rather rich phase diagrams consisting of several exotic phases and phase transitions that lie beyond the Landau-Ginzburg paradigm. 041049-2 In particular, one of the prominent phases for N>1 corresponds to 2N gapless Dirac fermions coupled to
an emergent Z2 gauge field in its deconfined phase. However, unlike a conventional Z2 gauge theory, we do not impose “Gauss’s Law” by hand; instead, it emerges because of spontaneous symmetry breaking. Correspondingly, unlike a conventional Z2 gauge theory in two spatial dimensions, our models have a finite-temperature phase transition associated with the melting of the order parameter that dynamically imposes the Gauss’s law constraint at zero temperature. By tuning a parameter, the deconfined phase undergoes a transition into a short-range entangled phase, which corresponds to Néel antiferromagnet or superconductor for N 2 and a valence-bond solid for N 3. Furthermore, for N 3, the valence-bond ¼ ¼ ¼ solid further undergoes a transition to a Néel phase consistent with the deconfined quantum critical phenomenon studied earlier in the context of quantum magnets.
DOI: 10.1103/PhysRevX.6.041049 Subject Areas: Condensed Matter Physics, Strongly Correlated Materials
I. INTRODUCTION structure in their wave functions, leading to the infamous Monte Carlo fermion sign problem [15 17]. Similarly, Ground states of strongly interacting electronic systems – fractional quantum Hall phases again possess a nontrivial can exhibit an extremely rich variety of phases. One way sign structure to their wave functions [18]; therefore, they are to organize our understanding is to classify the zero- temperature phases by their entanglement structure at so far amenable only via techniques such as exact diago- nalization [19,20] and the density matrix renormalization long distances and low energies [1–9]. Gapped phases that possess a local order parameter are characterized by short- group [21], which are restricted to only small two- range entanglement; that is, the reduced density matrix of a dimensional systems because of exponential scaling of large subsystem A can be understood simply by patching numerical cost with system size. However, long-range entangled phases exist that do not suffer from the together density matrices of smaller subsystems Ai whose union is A [10]. This is no longer true for gapless phases Monte Carlo sign problem. Two prominent examples are (1) interacting Dirac fermions with an even number of such as Fermi liquids [11–13] or gapped topological phases such as a fractional quantum Hall liquid, and such flavors [22] and (2) a gapped Z2 topological ordered system phases are thus said to possess “long-range entanglement” such as a toric code Hamiltonian [14] in a magnetic field [7,8,14]. [23]. The absence of a sign problem for the former is related Experience as well as heuristic arguments suggest that to the positive fermion determinant, while in the latter Hamiltonians whose ground states possess long-range case, the Hamiltonian has nonpositive off-diagonal elements entanglement are relatively difficult to simulate on a in a local basis, allowing one to sample the corresponding classical computer. For example, even a phase as ubiqui- Boltzmann weight efficiently [24]. In this paper, we study a tous and as well understood as a Fermi liquid is rather hard sign-problem-free model that has a ground state with to simulate numerically because fermions at finite density features of both of the aforementioned long-range entangled with repulsive interactions tend to have an intricate sign phases together, namely, Dirac fermions coupled to a fluctuating Z2 gauge field. By tuning parameters in the Hamiltonian, we also study competition with symmetry- breaking short-range entangled phases, which leads to novel, Published by the American Physical Society under the terms of the Creative Commons Attribution 3.0 License. Further distri- strongly interacting, quantum critical phenomena. bution of this work must maintain attribution to the author(s) and Our Hamiltonian [Eq. (1)] is partially motivated by a the published article’s title, journal citation, and DOI. recent study of nematic instability of a Fermi surface [25]
2160-3308=16=6(4)=041049(19) 041049-1 Published by the American Physical Society F. F. ASSAAD and TARUN GROVER PHYS. REV. X 6, 041049 (2016) and velocity fluctuations in Dirac metals [26]. It consists of N species of fermions on the vertices of a square lattice interacting with quantum Ising spins that live on the links (d) (e) of the same square lattice. However, in contrast to the Hamiltonian in Ref. [25], our Hamiltonian has a local Z2 (a) invariance, which leads to LxLy conserved operators on a Lx × Ly lattice. Furthermore, we restrict ourselves to half- filling; thus, at low energies, instead of a Fermi surface, we obtain either Fermi points or no fermions (e.g., a Mott insulator). These differences completely alter the phase diagram of our model compared to Ref. [25]. The phase diagram for various values of N is summa- rized in Fig. 1. For N>1, we find that when the kinetic energy of the aforementioned Ising spins is small compared to the kinetic energy of the fermions, the ground state resembles the deconfined phase of the Z2 gauge theory coupled to 2N Dirac fermions (the N 1 case is special in that the deconfined fermions are unstable¼ to infinitesimal interactions). However, as we will discuss in detail, there is an important distinction between this phase and the (b) (c) deconfined phase of a conventional Z2 gauge matter theory (Refs. [27–29]). Namely, the local Z2 invariance in our FIG. 1. Schematic zero-temperature phase diagram of the model is an actual symmetry of the Hamiltonian, in model in Eq. (1) (a), as well as cartoons (b)–(e) of a selected contrast to the so-called “gauge symmetry” in gauge-matter number of phases. (a) We observe a Z2 Dirac deconfined phase theories, which just corresponds to redundancy in the (Z2D), a Néel antiferromagnet phase (AFM) [or a superconductor description of the physical states. This is because we do (SC), depending on the pattern of particle-hole symmetry break- not project the Hilbert space of our Hamiltonian to gauge- ing], a charge density wave (CDW) phase, and a valence-bond solid (VBS). For N 1, we do not find evidence for a Z D phase invariant states; that is, we do not impose Gauss’s law. ¼ 2 Instead, in our model, Gauss’s law in the ground state beyond h 0, consistent with the arguments in the main text. The phase¼ transitions from the Z D to AFM/SC (N 2) and emerges because of spontaneous symmetry breaking as the 2 VBS (N 3) are seemingly continuous. At N 3, we observe¼ a locally conserved operators order in a certain definite deconfined¼ quantum critical point (DCQP) between¼ the VBS and pattern. One physical significance of this result is that, AFM phases. (b)–(e) Cartoons of the corresponding phases. in contrast to regular Z2 gauge theory, in our model, equal- Circles correspond to fermions, and the color code corresponds space, unequal-time, non-gauge-invariant correlation func- to the flavor index. Circles with two colors represent a pair of tions do not vanish identically. Furthermore, there is a fermions on a site with corresponding flavors. The low-energy finite-temperature Ising transition corresponding to the properties of the Z2D phase resemble SU N fermions propa- ð Þ restoration of the aforementioned symmetry. gating freely in space-time and connected by a Z2 gauge string As the kinetic energy of the Ising degrees of freedom (b). The symmetry-broken phases correspond to the confined is increased, we find that the aforementioned gapless phases of the model. At N 3, the AFM phase (e) has the fundamental (conjugate) representation¼ of SU(3) on sublattice Dirac deconfined phase enters a symmetry-broken confined A (B). The corresponding Young tableaux are included. The VBS phase. The nature of symmetry breaking depends on the phase (d) corresponds to a pattern of intersite SU(3) singlets. value of N. For N 1, it corresponds to a charge density wave, while for N ¼ 2, because of an enlarged symmetry, the symmetry-broken¼ phase can correspond to a super- II. MODEL AND SYMMETRIES conductor or a Néel AFM, depending on the sign of the Our model reads infinitesimal symmetry-breaking field. The N 3 case is Fermion-Boson problems: Unconstrained gauge theories even richer. We find that after exiting the deconfined¼ phase, N the model enters a VBS phase, and as the tuning parameter ˆ ˆ † ˆ H Z i;j cˆi;αcˆj;α H:c: − Nh X i;j : 1 controlling the kinetic energy of Ising spins is tuned further, ¼ h i þ h i ð Þ i;j α 1 i;j the valence-bond-solid quantum melts and yields a Néel Xh i X¼ Xh i AFM phase. The phase transition between the VBSF. F. ASSAAD and and TARUN GROVER PHYS. REV. X 6, 041049 (2016) Sign free for allHere, values cofˆ† N.creates a fermion of flavor α in a Wannier state AFM phase, if second order, corresponds to a deconfinedand velocity fluctuations ini;α Dirac metals [26]. It consists of quantum critical point [30], reported earlier in theN=2nN contextspecies Antiunitary of fermionscenteredsymmetry on aroundthe vertices squaring lattice of ato square site -1 i of lattice a square lattice, and the sum of quantum magnets [31,32]. Our results are consistentinteracting withruns quantum over Ising bonds spins of that nearest live on neighbors. the links Each(d) bond, b i; j , (e) with a second-order transition. N=2n+1of the same Majonara squareaccommodatesrepresentation lattice. However, a à two-level inO(2N) contrast system to the spanned by the¼ Hilberth i Hamiltonian in Ref. [25], our Hamiltonian has a local Z2 (a) invariance, which leads to LxLy conserved operators on a Lx × Ly lattice. Furthermore, we restrict ourselves to half- filling; thus, at low energies, instead of a Fermi surface, we obtain041049-2 either Fermi points or no fermions (e.g., a Mott insulator). These differences completely alter the phase diagram of our model compared to Ref. [25]. The phase diagram for various values of N is summa- rized in Fig. 1. For N>1, we find that when the kinetic energy of the aforementioned Ising spins is small compared to the kinetic energy of the fermions, the ground state resembles the deconfined phase of the Z2 gauge theory coupled to 2N Dirac fermions (the N 1 case is special in that the deconfined fermions are unstable¼ to infinitesimal interactions). However, as we will discuss in detail, there F. F. Assaad and Tarun Grover, Phys. Rev. X 6 (2016), 041049. is an important distinction between this phase and the (b) (c) S.deconfined Gazit, F. F. Assaad, phase ofS. Sachdev, a conventional A. Vishwanath,Z gauge and matter C. Wang, theoryPNAS, (2018) 2 (Refs. [27–29]). Namely, the local Z2 invariance in our FIG. 1. Schematic zero-temperature phase diagram of the model is an actual symmetry of the Hamiltonian, in model in Eq. (1) (a), as well as cartoons (b)–(e) of a selected contrast to the so-called “gauge symmetry” in gauge-matter number of phases. (a) We observe a Z2 Dirac deconfined phase theories, which just corresponds to redundancy in the (Z2D), a Néel antiferromagnet phase (AFM) [or a superconductor description of the physical states. This is because we do (SC), depending on the pattern of particle-hole symmetry break- not project the Hilbert space of our Hamiltonian to gauge- ing], a charge density wave (CDW) phase, and a valence-bond solid (VBS). For N 1, we do not find evidence for a Z D phase invariant states; that is, we do not impose Gauss’s law. ¼ 2 Instead, in our model, Gauss’s law in the ground state beyond h 0, consistent with the arguments in the main text. The phase¼ transitions from the Z D to AFM/SC (N 2) and emerges because of spontaneous symmetry breaking as the 2 VBS (N 3) are seemingly continuous. At N 3, we observe¼ a locally conserved operators order in a certain definite deconfined¼ quantum critical point (DCQP) between¼ the VBS and pattern. One physical significance of this result is that, AFM phases. (b)–(e) Cartoons of the corresponding phases. in contrast to regular Z2 gauge theory, in our model, equal- Circles correspond to fermions, and the color code corresponds space, unequal-time, non-gauge-invariant correlation func- to the flavor index. Circles with two colors represent a pair of tions do not vanish identically. Furthermore, there is a fermions on a site with corresponding flavors. The low-energy finite-temperature Ising transition corresponding to the properties of the Z2D phase resemble SU N fermions propa- ð Þ restoration of the aforementioned symmetry. gating freely in space-time and connected by a Z2 gauge string As the kinetic energy of the Ising degrees of freedom (b). The symmetry-broken phases correspond to the confined is increased, we find that the aforementioned gapless phases of the model. At N 3, the AFM phase (e) has the fundamental (conjugate) representation¼ of SU(3) on sublattice Dirac deconfined phase enters a symmetry-broken confined A (B). The corresponding Young tableaux are included. The VBS phase. The nature of symmetry breaking depends on the phase (d) corresponds to a pattern of intersite SU(3) singlets. value of N. For N 1, it corresponds to a charge density wave, while for N ¼ 2, because of an enlarged symmetry, the symmetry-broken¼ phase can correspond to a super- II. MODEL AND SYMMETRIES conductor or a Néel AFM, depending on the sign of the Our model reads infinitesimal symmetry-breaking field. The N 3 case is even richer. We find that after exiting the deconfined¼ phase, N the model enters a VBS phase, and as the tuning parameter ˆ ˆ † ˆ H Z i;j cˆi;αcˆj;α H:c: − Nh X i;j : 1 controlling the kinetic energy of Ising spins is tuned further, ¼ h i þ h i ð Þ i;j α 1 i;j the valence-bond-solid quantum melts and yields a Néel Xh i X¼ Xh i AFM phase. The phase transition between the VBS and † AFM phase, if second order, corresponds to a deconfined Here, cˆi;α creates a fermion of flavor α in a Wannier state quantum critical point [30], reported earlier in the context centered around lattice site i of a square lattice, and the sum of quantum magnets [31,32]. Our results are consistent runs over bonds of nearest neighbors. Each bond, b i; j , with a second-order transition. accommodates a two-level system spanned by the¼ Hilberth i
041049-2 RESEARCH NEWS & VIEWS conditions. Both groups agree that, under Judith Berman is in the Department of 4. Hede, K. J. Natl Cancer Inst. 97, 87–89 (2005). conditions optimized by geneticists for growth Genetics, Cell Biology, and Development, and 5. Rancati, G. et al. Cell 135, 879–893 (2008). 6. Selmecki, A., Gerami-Nejad, M., Paulson, C., Forche, in conventional laboratories, aneuploid cells the Department of Microbiology, University of A. & Berman, J. Mol. Microbiol. 68, 624–641 (2008). usually divide less rapidly than cells with the Minnesota, Minneapolis, Minnesota 55455, 7. Torres, E. M. et al. Science 317, 916–924 (2007). normal chromosomal complement. USA. 8. Torres, E. M. et al. Cell 143, 71–83 (2010). An important distinction in the methods e-mail: [email protected] 9. Geiger, T., Cox, J. & Mann, M. PLoS Genet. 6, used to generate the aneuploid strains might e1001090 (2010). 2,7,8 1. http://8e.devbio.com/article.php?ch=4&id=24 10. Springer, M., Weissman, J. S. & Kirschner, M. W. explain the differences in the findings . 2. Pavelka, N. et al. Nature 468, 321–325 (2010). Mol. Syst. Biol. 368 (2010). 7,8 6, Torres et al. engineered yeast cells with a 3. Duesberg, P., Li, R., Fabarius, A. & Hehlmann, R. 11. Weaver, B. A. A. & Cleveland, D. W. Cancer Res. 67, haploid (single) set of chromosomes to carry Contrib. Microbiol. 13, 16–44 (2006). 10103–10105 (2007). one extra chromosome and then selected for faster growth using conventional lab condi- tions for 9–14 days — a time frame during HIGH-TEMPERATURE SUPERCONDUCTIVITY which mutations are expected to accumulate. By contrast, Pavelka et al. analysed strains that often carried multiple aneuploid chromosomes and, importantly, minimized the number of Mind the pseudogap generations before analysis. This illustrates a crucial truism of experimental genetics: you The discovery of predicted collective electronic behaviour in copper-oxide get what you select for. superconductors in the non-superconducting state provides clues to unlocking the Pavelka and co-workers also directly 24-year-old mystery of high-temperature superconductivity. SEE LETTER P.283 address a controversy concerning the role of excess proteins in aneuploid cells. Previously, Torres et al.7 proposed that there is a specific CHANDRA VARMA in the dome-shaped superconductivity region set of genes and proteins that are regulated in of the phase diagram (Fig. 1). This point would response to aneuploidy in general. In their he phenomenon of high-temperature occur at zero temperature and would involve more recent study8, they showed that some superconductivity is a beautiful and a change in the symmetry of the mater ials’ 20% of proteins exhibit levels that do not well-posed scientific problem with electronic structure. Because Tc is determined track with gene copy number, and that a large Tmany facets. On page 283 of this issue, Li et al.1 by the materials’ collective electronic excita- proportion of these proteins are members of report observing a special kind of intense col- tions in the non-superconducting state, it is macromolecular complexes. By contrast, other lective electronic fluctuation in the most un arguable that the coupling of electrons to groups9,10 have found that the levels of most mysterious phase of matter exhibited by high- such excitations in the strange-metal region proteins generally reflect changes in chromo- temperature superconducting copper-oxide and their modifications in the pseudogap some copy number and that less than 5% of materials (cuprates). Taken together with pre- region lead to high-FermionTc superconductivity.-Boson problems: ”De-signer” Hamiltonians the proteins exhibit ‘dosage compensation’ vious experimental2–6 and theoretical7 work, If it exists, a quantum critical point in the — whereby the relative protein level is inde- this observation significantly narrows the pendent of gene-copy number. Pavelka et al. range of directions likely to be fruitful in the La2-xSrxCuO4, HgBa2CuO4+δ specifically test this hypothesis by quantitative quest to understand high-temperature super- Insulator and mass spectrometry of about 2,000 proteins in conductivity. The authors performed their antiferromagnetic each of five aneuploid strains and do not find experiments in two samples of HgBa2CuO4+δ, Hˆ = t cˆ† cˆ +ˆc† cˆ compelling evidence for specific dosage com- which has one of the simplest crystal structures i, j, j, i, ? Strange Temperature i,j , pensation of protein-complex components. of any of the cuprate families and is ideal for T* metal Crossover hXi ⇣ ⌘ Overall, these studies2,7–10 are consistent with such studies. 1 + U (ˆni, 1/2) (ˆni, 1/2) µ nˆi, the idea that aneuploidy is not a single, unique Li and colleagues’ experiments pertain to " # i i, Fermi AAAELnichVFdb9MwFL1Z+NjKxzp45MWiYhqsLUlf4AVpEh/aC9KQ6DZpHpWTuqnXxAmOs6lE+U/8Df4AEkiIV954BYkHbtxkW+nQHDn3+txzjq9tLwlFqh3ni7VkX7l67frySuPGzVu3V5trd3bTOFM+7/txGKt9j6U8FJL3tdAh308UZ5EX8j1v8rys7x1zlYpYvtXThB9GLJBiJHymERqsWR71eCBkzt9LphSbPioadMx0vl2QZ2SdkI7GP02zaJDTkMkg5ES0j6gyaRsrIogYKQghNOQjvUGM2i8GuWjPisW7nA5ZEHBVnBWPzooF2bwQXxSdcyyoEsFYPyz3pQRb3CTr/apPUcz1ImfKLMHzxScF6biPe6fqC4jD+ETWVIJcUpM7hEZZvUfdypx0BjUol8NzFzpotpyuYwZZTNwqaUE1duLmV6AwhBh8yCACDhI05iEwSPE7ABccSBA7hBwxhZkwdQ4FNFCbIYsjgyE6wX+Aq4MKlbguPVOj9nGXEKdCJYEHOF8ZRw/Z5a4c8xTjH5wfDBb8d4fcOJcdTjF66LhiHF8jrmGMjMuUUcWse7lcWZ5KwwiemtMI7C8xSHlO/9TnBVYUYhNTIfDSMAP08Mz6GG9AYuxjB+Ut1w7EnHiIkZnIjYusHBn6KYzl7WM/+Mzuv4+6mOz2uq7Tdd/0Wltu9eDLcA/uwwa+6hPYgm3YwT5865P10/pl/bY/2p/tb/b3GXXJqjR3YW7YP/4Cd+4SMQ== state and that all aneuploid strains do not share the pseudogap region of the phase diagram of Pseudo- liquid X X a single, common phenotype or protein profile. the cuprates (Fig. 1), a sort of precursor state gap phase Tc Rather, different aneuploid strains use differ- to the superconducting phase that most con- ent mechanisms for optimal growth under densed-matter physicists regard as the Rosetta Doping different conditions. This conclusion may be Stone for discovering the physical principles Superconductivity QCP less satisfying than a single, simple answer, that underlie the cuprates’ behaviour. On C. Varma, Nature 468 (2010). especially given the crucial implications for entering the pseudogap region, at a temper- Figure 1 | Phase diagram of the cuprates. At very cancer cells: it remains unclear whether ature below T* but above the temperature à Frac%onaliza%onlow levels of electron–hole doping,Emergent cuprates are gauge theories cancer cells divide uncontrollably because below which superconductivity emerges (T ), c insulating and antiferromagneticSpin liquids (the (Z materials’2 gauge theory in deconfined phase) they are aneuploid and/or because they have all cuprates’ thermodynamic and electronic- neighbouring spins point in opposite directions). accumulated mutations that allow them to transport properties change by a large amount At increased doping levels, they become tolerate aneuploidy. But it should be remem- owing to the materials’ loss of low-energy Intertwinedconducting, and theorders exact temperature à Superconduc%vity and doping versus stripes bered that work on cancer cells themselves11 electronic excitations. level determine whichSuperconduc%vity phase of matter they will versus an%ferromagne%sm suggests that not all aneuploidies are equal: The pseudogap region is bounded on one be in. At temperatures below Tc, they become aneuploidy can either promote or inhibit side by a region of remarkably simple but superconducting, and at temperatures above Tc but tumorigenesis, depending on the context. unusual properties, which do not fit into the Quantumbelow T* they phasefall into the transi%ons pseudogap phase. à TheNema%c transi%ons Pavelka and colleagues’ work2 therefore Fermi-liquid-type model that has been used to boundary of the pseudogap region at low doping levels is unknown. The transition(See also between heavy the fermions) supports the idea that, whereas mutations describe metals at low temperatures for about Fermi-liquid phase and the strange-metal phase can facilitate the proliferation of aneuploid a hundred years. Some researchers got to grips occurs gradually (by crossover). QCP denotes the cells, aneuploidy itself can be sufficient to with understanding this ‘strange-metal’ region quantum critical point at which the temperature à Designer models: are numerically1 tractable and capture aspects of the above physics provide a growth advantage under a broad early in the history of high-Tc superconduct- T* goes to absolute zero. Li and colleagues’ study range of stress conditions. ivity, by hypothesizing a quantum critical point pertains to the pseudogap phase.
184 | NATURE | VOL 468 | 11 NOVEMBER 2010 © 2010 Macmillan Pu blishers Limited. All rights reserved 3
Sec. I we define the lattice Hamiltonian and describe its dimensionless coupling strength between the pseudospin phase diagram; in the next section, we provide evidence and the fermion bond density. for the statementsIsing nematic above regarding quantum bosonic critical correlations point in a metal: a Monte Carlo study (III A), superconductivity (III C), and single-fermion cor- The pseudospins are not related to the physical spins relations (III D); finally, in Sec. 1,IV, we discuss various2, of the electrons; their2 ordering corresponds1 to a nematic Yoni Schattner, ⇤ Samuel Lederer, ⇤ Steven A. Kivelson, and Erez Berg caveats concerning1 the interpretation of our results, and transition. The model (3) should be viewed as an e↵ec- Department of Condensed Matter Physics, The Weizmanntive Institute lattice model, of Science, designed Rehovot, to give 76100, a nematic Israel QCP. Mi- their bearing on both2 priorDepartment theoretical of Physics, work and Stanford the in- University, Stanford, CA 94305, USA croscopically, the nematic degrees of freedom could rep- terpretation of experiments. (Dated: November 12, 2015) resent (via Hubbard-Stratonovich transformation) inter- 3 The Ising nematic quantum critical point (QCP) associatedactions withinvolving the zero the temperature same electrons transition that form the Fermi II.from MODEL a symmetric AND to PHASE a nematic DIAGRAMmetal is an exemplar ofsurface, metallic quantum or another, criticality. independent We have degree carried of freedom (such Sec. I we defineout a minus the lattice sign-free Hamiltonian quantum Monte and Carlo describe study its of thisasdimensionless a QCP phonon for a mode two coupling dimensional that becomes strength lattice soft betweenmodel at a structural the pseudospin tran- phase diagram;with sizes in the up nextto 24 section,24 sites. we The provide system evidence remains non-superconductingsition).and the So fermion long asdown bond the to properties density. the lowest of acces- the QCP are uni- Our modelsible is temperatures. illustrated in⇥ The Fig. results3, and exhibit is defined scaling on behavior over the accessible ranges of temperature, fora two-dimensional the statements square above lattice. regarding Every bosonic lattice correlations site hosts versal, the low-energy behavior does not depend on its (imaginary) time, and distance. This scaling behavior hasThe remarkable pseudospins similarities are not with related recently to the physical spins (aIII single A), superconductivity (spinful)measured fermionic properties ( degreeIII of C the), andof Fe-based freedom. single-fermion superconductors Each link cor- proximatemicroscopic to theirorigin. putative nematic QCP. relationshas a pseudospin-1 (III D); finally,/2 degree in Sec. of freedomIV, we that discuss couples various to of the electrons; their ordering corresponds to a nematic caveatsthe fermion concerning bond-density. the interpretation The system of is our described results, by and transition.For ↵ = 0, The the system model ( is3) composed should be of viewed two decoupled as an e↵ec- theirthe following bearing Hamiltonian: onI. both INTRODUCTION prior theoretical work and the in-maticsetstive order of latticedegrees [54– model,59 of] freedom: and designed nematic free quantum fermions, to give a critical and nematic pseudospins fluctu- QCP. Mi- terpretation of experiments. ationsgovernedcroscopically, [60–65 by] haveHb, the which been nematic hasobserved the degrees form in many of of a d freedom of= the 2 transverse mate- could2 rep- H = Hf + Hb + Hint, (3) rialsfieldresent mentioned Ising (via model. Hubbard-Stratonovichabove. At zero temperature, transformation) the pseudospins inter- A hallmark of strongly correlatedthe Matsubara electron frequency, systems and isA, b, and are positive dimensionful constants. We find that theundergo exponents in a second-orderTABLE I: Critical quantum exponents for a phase 2d quantum transition Ising QCP from Ouractions simulations involvingwith z =1(classical3 are the limited samed Ising), to electrons finite and for the system Ising that nematic sizes form QCP – the up Fermi thewhere competition of ground statesthis with expression di↵erent take values kinds =1 of.0 0.1 and =1.0 0.1. II. MODEL AND PHASE DIAGRAM ± toa 24surface, paramagnet±24 lattice orin another, a to metal sites. an from “antiferromagnet” Withinindependent our DQMC our simulations. numerical degree In thethat limit of accuracy, breaks freedom of van- 90 (such order.[1] In this context, a central(See set Fig. of8.) unsettled theo- ⇥ ishing coupling to the fermions, exponents from our DQMC This implies, among other things, thatwerotational find theas uniform a evidence phonon symmetry for mode a continuous thatat h = becomeshc0 nematic. This soft transition quantum at a structural phase is in the tran- retical issuesHf concerns= t the characterci† cj µ of quantumci† ci , critical simulations are consistent with 3d classical Ising values.z ˜ is susceptibility D(h, T, 0, 0) has a Curie-Weiss3 dimensional form the classical “apparent dynamical Ising universality exponent” defined belowclass. Eq. At (1). T = 0, Our model is illustratedi,j , Fermion in Fig. -Bosoni,3 , and⌘ isproblems: defined on transition”De-signer”sition). with So Hamiltonians critical long as exponents the properties that are of the significantly QCP are uni- points (QCPs) in metalshX [i2–4]. Suchwith metallic anX e↵ective WeissQCPs temperature have whichthe varies pseudospin linearly degrees of freedom are gapped in both the with (hc h), as can be seen in Fig. di1.Eq.↵erentversal,1 can from be the those low-energy of a nematic behavior QCP does in an not insulator. depend on its beena two-dimensional identified in several square heavy lattice.z fermionz Every compounds[ latticex site5 hosts]; /⌫ Critical Exponents ⌫ ⌘ z˜ Hb = V ⌧i,jviewed⌧k,l ash a scaling⌧i,j relation,, D(h, T, q, 0) =nematic⇠ (x, y), and isotropic phases. At finite temperatures, ⌫ z˜ Specifically,microscopic in the2d Ising origin. disordered QCP phase0.63 near1.24 the0.04 QCP,1 where.0 evidencea single for (spinful) quantum fermionic critical degree behaviorwhere of ⇠ freedom. hash hc also , x Each= beenT ⇠ , linky = q ⇠.InEq.1, Isingi,j ; k,l nematic quantum⇠ | i,j | critical point| | in a metal:2d metallic a MonteIsing ne- Carlo study is the conventionally defined susceptibilitythea exponent, dimensionless line of the second-order quantum classical control0.5 0.1 parameterd1.0=0. 21 0 Ising.0 h0.3 2 transitions.0hc0,the.3 has a pseudospin-1hh/2X degreei h ii of freedomhX thati couples to matic QCP found in the ruthenate Sr3Ru2O7 correlation[6, 7], and length the exponent cuprate⌫ /(2 ⌘), where ⌘ is ± ± ± ± z 1, 2,thermodynamicextendsFor from↵ = the0, (zero the2 QCP frequency) system in the ish nematic1 composedT plane. two-point of two decoupledcor- andthe iron-based fermionHint bond-density.= superconductors↵t ⌧Yonii,j c The† Schattner, [cthe8]. system anomalous. Quantum⇤ is dimensionSamuel described critical- of Lederer, the⌘(4) nematic by ⇤ Steven field, and A. we Kivelson, and Erez Berg i j relator (defined in Eq. 5) is consistent with the following Phys. Rev. X 6,1Department 031028i,j , (2016) of Condensedhave introduced Matter an “apparent Physics, dynamical The Weizmann exponent,”sets Institute ofz ˜ = degrees of Science, of freedom: Rehovot, free 76100, fermions, Israel and pseudospins itythe is followingalso often Hamiltonian: invokedh asi a possible1 explanation of the Forn non-zero ↵, the phase diagram is similar, but ex- X 2Department(⌫ ) . The values of Physics, of these Stanford exponents University,functional derivedgoverned from Stanford, our description: by H CA, 94305, which USA has the form of a d = 2 transverse “strange” or “bad” metal behaviorMonte seen Carlo in data a are variety given in of Table I. hibits quantitativeb and qualitative modifications. Fig. 1 Here, c† createsH a fermion= H + onH site+ToH illustratej with, the spin significance(Dated: = of, November(3)these, exponents,field 12, 2015) note Ising model. At zero temperature, the pseudospins such materialsj [9–16]. f b int shows the phase diagram, obtainedA by DQMC, for both The Ising nematicthat working quantum backwards critical from point the" # exponents (QCP) associated correspond- with the zero temperature transition i, j denotes a pair of nearest-neighboring to the d + z dimensional sites on classical the IsingD modelundergo(h, T, withq,i! an second-order= 0) quantum phase transition2 (1) from wherehTo date,i there existsfrom no satisfactory a symmetric to theory a nematic ofmetal metallicis an exemplar↵ of= metallic 0 and quantum↵ = criticality. 0.⇡ TheT We↵+ have=b(h 0 carried transitionhc)+ q between the square lattice, t and µ are thedynamical hopping exponent strengthz = 1 (also and listed in TableaI paramagnet) would 6 to an “antiferromagnet”6 | that| breaks 90 QCPs in d = 2 or 3 spatialout a dimensions minus sign-freeyielda although=1 quantum.96, =1 a Monte. number24, and Carlo 2 ⌘ study=1.96.nematic of Indeed,this QCP we and for a isotropic two dimensional phases lattice remains model second order, and chemical potential, respectively,with sizes† up tofind⌧ 24i,j critical(24a = sites.exponents†x, y, The z consistent) system denote remains withwhere this non-superconducting expectationrotationalT is the symmetry temperature, down to at thehq lowest=ish thec acces-0. wave-vector, This transition! isis in the of field theoreticHf = approachest ci have cj beenµ⇥ attemptedci ci , [17– extrapolates to a new QCP, shifted relative to hc0.n More Pauli matrices that actsible on temperatures. theby pseudospin DQMC The when results we that set exhibit to lives zero scaling the on coupling behavior between3 dimensional over the the accessible classical ranges Isingof temperature, universality class. At T = 0, 34]. Among the unsettledi,j , issues are:i, a) the values of striking is the change in the slope with which the phase the bond connectingh(imaginary) theXi neighboring time,critical and fluctuations sitesX distance.i and and Thisj the, V> fermions. scaling0 behaviorthe has pseudospin remarkable degrees similarities of with freedom recently are gapped in both the critical exponents, andmeasured which properties propertiesz Includingz of can the results be Fe-based expressed atx finite superconductors (Matsubara)boundary, frequencies, proximate2 T toN their(h), putative approaches nematic the QCP. QCP. For ↵ = 0, the is the nearest-neighborHb = V Ising interaction⌧i,jwe⌧k,l find ah wide between intermediate⌧i,j , neigh- regime wherenematic the full and isotropic phases. At finite temperatures, as scaling functions involving thesedynamical exponents, fluctuation b) the spectrum ex- is wellslope described diverges by at low temperature, consistent with the ex- boring pseudospins (here,i,j ; k,l i, j ; k, l denotesi,j a pair of a line of second-order classical d = 2 Ising transitions tent to which metallic QCPs arethe preempted simple “functional by approximant”super- D(h,pectation T, q,i!n) for the transverse field Ising transition, where hh XiI.h hh INTRODUCTIONii i h ii hXi maticT order⇡ [54–59] and nematic quantum critical fluctu- nearest-neighbor bonds), h V is( theh, T, q strength,i!n), where of a trans- extends from⌫z the QCP in the h T plane. conducting [35–43] or other formsz A [38, 44, 45] of auxil- TN 1.5 h hc0 with ⌫ =0.63 and z = 1. In contrast, Hint = ↵t ⌧i,j· ci† cj . (4) ations [60–65] have been observed in many of the mate- verse field that acts on the pseudospins, and ↵ is theA / | | iary order, that gap out the Fermi surface,(h, T, q,i! ) c) whether forrials↵ mentioned. >(2)0, we above. find that TN (hc h). On the disor- A hallmarki,j of, strongly correlated electronn systems is 2 Quantum there exists a “non-FermihXi liquid” metalA in the⌘ quantumT + b(h hc)+ q + c !OurnForn simulations non-zero are↵ limited, the phase to finite/ diagram system sizes is similar, – up but ex- the competition of ground states with di↵erent kinds of| | dered| | side of the transition,critical we define a crossover line by critical regime. Controlled theoreticalHowever, methods, in contrast with that the can thermodynamictohibits correlation 24 1 24 quantitative lattice sites. Within and qualitative our numerical modifications. accuracy, Fig. 1 Here, c† order.[creates1] In a this fermion context, on a site centralj with set of spin unsettled = theo-, , identifying⇥ hcross(T ) as the value of h at which the static be usedj both to benchmark thecj† function, field theories as shown in and Fig. 10 for, the" # dynamicalweshows find response evidence the phase for a continuous diagram, nematic obtained quantum by DQMC, phase for both i, j denotesretical a issues pair concerns of nearest-neighbor theexhibits character noticeable of quantum deviations sites on critical from the this formsusceptibility for the few FIG. at 2: fixed FermionT Green’shas function, fallenG(k to, ⌧ = half /2), as of a func-its value lowest values of q ; q 2⇡/L where L istransition the system withtion critical of momentum exponentsk across thatthe Brillouin are significantly zone for three dif- comparisonh i withpoints experiments, (QCPs) in metals are [2 greatly–4]. Such needed. metallic Deter- QCPs have at↵h== 0hc and.NematicTcross↵ =(h), 0. the The inverse↵ = of 0 transitionhcross, also between vanishes the square lattice, t and µ are thewidth. hopping It is unclearx,y,z| | strength whether| | ⇠ these and deviationsdi↵ shoulderent0.5 be fromferent those values6 of of a the nematic transverse QCP6 field: inh =3 an.35 insulator. (in the nematic minant quantumbeen identified Monte Carlo in several (DQMC)[ heavy fermion46⌧ , 47]–in compounds[ par- 5]; nematic andphase), isotropich =3.525 (near phases the QCP), remains and h =3 second.7(inthedis- order, and chemical potential, respectively,dismissed⌧ a ( asa finitei,j= x, size y, e z↵ects,) denote or taken aslinearlySpecifically, indicators ofwith in theh upon disordered approaching phase near the theSymmetric QCP, QCP, where although its evidence for quantum criticaladi↵i,jerent behavior form of dynamical has also scaling been closer to criticality, ordered phase). In the nematically ordered phase, the data ticular, in cases whereV the fermion sign problem can slopetheextrapolates dimensionless is steeper to quantum than a new that control QCP, of T parameterN shifted(h). relativeh hc,the to hc0. More Pauli matricesfound in that the actruthenate on the Sri.e.Ru pseudospin atO lower[6, temperatures7], and that thelives and cuprate larger on system sizes than is averaged over both orientations of the order parameter. 3 2 7 G(k, ⌧ = /2) is proportional to the integral over the spec- be circumvented [48–50]–may serveare presently this purpose. accessible. Taking Such Eq. (2) at facethermodynamicstriking value,0 one is the (zero change frequency) in the nematic slope two-point with which cor- the phase the bond connectingand iron-based the superconductors neighboring sites[8]. Quantumi and j, critical-V>0 2 tral2.5 function A(k,3!) over an3.5 energy window4 of width T 4.5[see methods have been successfully appliedfinds the dynamical in several critical prob- exponent z relator2,The as would (defined linear behavior in Eq. 5) is of consistentT (h) for with small the followingh h is also ity is also often invoked as a possible explanation of the ⇡ boundary,Eq.TN (11()].h), The approaches temperatureN is T =1 the/8, QCP. and the systemc For size↵ = 0, the is the nearest-neighbor Ising interactionbe obtained from between a naive dynamical neigh- scalingfunctional assumption description: hc0 hc lems in which“strange” critical or bosonic “bad” metal fluctuationst behavior are seen coupled in a variety to of seenslope for diverges otheris L model= at20. Thelow parameters. data temperature, were taken from In systems consistentFig. with4hwe either show with pe- the ex- boring pseudospins (here, i, jz˜ =; zk,. This l lastdenotes result is at a odds pair with of theoretical expec- riodic or anti-periodic boundary conditions. The Fermi sur- fermions withsuch a materialsDirac-like [9– dispersion16hh]. tations,i h [51 [102ii–53] and]. so requires special analysis,phasepectation especially diagramfaces for for are the twoseen transverse as values peaks in ofG(Ak field the, ⌧ = fermion Ising/2). The transition, lowerdensity, right con- where given that there are reasons to question the validity of panel shows the Fermi surface of the bare band structure. nearest-neighbor bonds), h V is the strength of a trans-FIG.D 1:(h, Phase T, q,i! diagramn = 0)⌫z of the model obtained2 by(1) DQMC, arXiv:1511.03282v1 [cond-mat.supr-con] 10 Nov 2015 To date, there exists no satisfactory theory of metallic trolled by the chemical potential µ. As the fermion den- In this paper, we report a· DQMCnaive dynamical study scaling of a [3, two-16]. TN h hc0 ⇡withT +⌫b=0(h .h63c)+ and qz = 1. In contrast, At µAAAC0HichVFLS8NAEJ7GV1tfVY9egkXwVBIR9CIUfOBFqGIf0BZJ0m1cmhebpFiLiFdvXvWX6W/x4LdrKmgRN2xm9ptvvp3ZsSOPx4lhvOW0mdm5+YV8obi4tLyyWlpbb8RhKhxWd0IvFC3bipnHA1ZPeOKxViSY5dsea9qDIxlvDpmIeRhcJaOIdX3LDXifO1YCqNnxU/1QN65LZaNiqKVPO2bmlClbtbD0Th3qUUgOpeQTo4AS+B5ZFONrk0kGRcC6NAYm4HEVZ3RPReSmYDEwLKAD/F2c2hka4Cw1Y5Xt4BYPWyBTp23sU6Vogy1vZfBj2A/sO4W5f94wVsqywhGsDcWCUjwHntANGP9l+hlzUsv/mbKrhPp0oLrhqC9SiOzT+dY5RkQAG6iITieK6ULDVuchXiCAraMC+coTBV113IO1lGVKJcgULegJWPn6qAdjNn8Pddpp7FZMo2Je7JWrRjbwPG3SFu1gqvtUpTOqoQ7Z4TO90Kt2qd1qD9rjF1XLZTkb9GNpT59rgJCm sha1_base64="I3qzNHs0O153gIiFHHCMQ4HedNo=">AAACu3ichVLLSsNAFD2Nr7ZWrWs3wSK4KokbXQo+EEGoYB9QiyTpNIamSZhMK7X4A279Ov0WF56MqaBFOmFy75x777mPGTcJg1RZ1nvBWFldW98olsqblfLW9k610krjsfRE04vDWHZcJxVhEImmClQoOokUzsgNRdsdnmX29kTINIijOzVNRG/k+FEwCDxHEWo8VGtW3dLLXFTsXKkhX3H1A/foI4aHMUYQiKCoh3CQ8uvChoWEWA8zYpJaoO0CLygzdkwvQQ+H6JB/n6dujkY8Z5ypjvaYJeSWjDRxwH2pGV16Z1kF9ZTyk/tZY/6/GWaaOatwSumSsaQZb4grPNJjWeQo95zXsjwy60phgBPdTcD6Eo1kfXo/POe0SGJDbTFxoT19crj6POEEIsomK8imPGcwdcd9SkdLoVminNEhn6TMps96eMv23ztdVFpHdduq27cWitjDPg55mcc4xRUaTO8xySvejGsjMZ6+X4NRyJ/FLn4tY/oFhOiMYg== sha1_base64="ZGB5CT+PsS1f6/CwOpGBa+beq2c=">AAACxXichVLLLgRBFD3aa4zBsLXpEInVpNuGjUTiERsJYswkiFT3lFGZfqW6RzAR8Qm2fBnfYuFUaRJEVKf63jr33nMfVUEWqbzwvJchZ3hkdGy8MlGdrE1Nz9Rna8d52tehbIZplOp2IHIZqUQ2C1VEsp1pKeIgkq2gt2nsrSupc5UmR8VNJs9i0U3UhQpFQah1Gvfdddc7ry96Dc8u97fil8oiyrWf1l9xig5ShOgjhkSCgnoEgZzfCXx4yIidYUBMU1PWLnGHKmP79JL0EER7/Hd5OinRhGfDmdvokFkibs1IF0vcO5YxoLfJKqnnlG/ctxbr/plhYJlNhTeUARknLOMe8QKX9PgvMi49P2v5P9J0VeACa7Ybxfoyi5g+wy+eLVo0sZ61uNi2nl1yBPZ8xQkklE1WYKb8yeDajjuUwkppWZKSUZBPU5rpsx5es//zUn8rxysN32v4Bx4qmMcClnmZq9jALvaZ3jT2iCc8O4fOtXP/8SCcofJlzOHbch7eAS0nj3w= sha1_base64="gMuuCJSeRn+qX2ljhfq7jG2igVI=">AAAC0HichVFLS8NAEJ7GV1tfVY9egkXwVDZe9CIUfOBFqGIf0BZJ0m0MzYtNUqxFxKs3r/rL9Ld48Ns1FbRIN2xm9ptvvp3ZsSLPjRPG3nPa3PzC4lK+UFxeWV1bL21sNuIwFTav26EXipZlxtxzA15P3MTjrUhw07c83rQGxzLeHHIRu2FwnYwi3vVNJ3D7rm0mgJodP9WPdHZTKrMKU0ufdozMKVO2amHpgzrUo5BsSsknTgEl8D0yKcbXJoMYRcC6NAYm4LkqzumBishNweJgmEAH+Ds4tTM0wFlqxirbxi0etkCmTrvYZ0rRAlveyuHHsJ/Y9wpz/r1hrJRlhSNYC4oFpXgBPKFbMGZl+hlzUsvsTNlVQn06VN24qC9SiOzT/tE5QUQAG6iITqeK6UDDUuchXiCAraMC+coTBV113IM1leVKJcgUTegJWPn6qAdjNv4Oddpp7FcMVjEuWbnKsoHnaZt2aA9TPaAqnVMNdcgOX+iV3rQr7U571J6+qVouy9miX0t7/gJqQJCi =0 as , the a function model has of the the transversesame symmetries field h and as |the temperature| SSH modelT . verse fieldQCPs that in actsd = 2 on or 3 the spatial pseudospins, dimensionsIn the fermionic although and sector we↵ a number findis results the consistentsity is/ with reduced,| both| h h and the range over which dimensional sign problem-free lattice model that exhibits Here,forV ↵=>t, 0,↵ we=0 find.5, µ thatc= Tc00.N5t. The(hc solidh). line On marks the disor- of field theoretic approachesstrongly have renormalized been attempted FermiAt liquid, µ [ 17=0 “marginal – , theTwhere Fermi( theh)T is liq-O(2N)is linear thea smallsymmetry temperature, become power of smaller, Tis .q Fig.brokenis2 theshows indicating/ wave-vector, down maps of to the thatU(N) low-energy!n is the e↵ect an “Ising nematic” QCP in a metal at finite fermionAAAC03ichVFLS8NAEJ7GV1tfVY9egkXwVDYi6LHgAy9CBdMW2yKbdFtD83KTFGrxIl69edX/pb/Fg1/WVNAi3bCZ2W9mvnlZoetEMWPvOW1ufmFxKV8oLq+srq2XNjbrUZBIW5h24AayafFIuI4vzNiJXdEMpeCe5YqGNThO7Y2hkJET+FfxKBQdj/d9p+fYPAZ03fYSve2LO53dlMqswtTRpxUjU8qUnVpQ+qA2dSkgmxLySJBPMXSXOEX4WmQQoxBYh8bAJDRH2QU9UBGxCbwEPDjQAf59vFoZ6uOdckYq2kYWF1ciUqdd3DPFaME7zSqgR5CfuPcK6/+bYayY0wpHkBYYC4rxAnhMt/CYFellnpNaZkemXcXUoyPVjYP6QoWkfdo/PCewSGADZdHpVHn2wWGp9xAT8CFNVJBOecKgq467kFxJoVj8jJGDT0Km00c9WLPxd6nTSn2/YrCKcXlQrrJs4Xnaph3aw1YPqUrnVEMdNrK80Cu9aaY21h61p29XLZfFbNGvoz1/ATnfkic= Ndered side of the transition, we define a crossover line by 34]. Among the unsettleduid” issues [66], or are: weakly a) non-Fermi the values liquid6 of behaviorthe transition down to temperature,spectral weightT asN a,betweenthenematicandthe function of momentum in the dis- our lowest temperatures (T 0.02E ,whereof theE couplingis ordered between phase, at electrons the QCP, and near in the the ordered Fermi phase. surface density; incritical the nematic exponents, phase, and which the properties discrete can lattice be expressedmin ro-⇡ symmetricF identifyingF phases.hcross The(T line) as extrapolates the value oftoh theatT which=0QCP the static FIG. 3: Illustration of the latticec modelj†the Fermi (3 energy).). Spinful In particular, electrons at the QCP,and the the e↵2ec- nematicThe existence[ modes67] becomes of “cold-spots” weaker. along the The zone diago- fact that tational symmetryas scaling is functions spontaneously involvingtive quasiparticle thesebroken exponents, from weight,CZ b)4 ( thetoT ), defined ex- at inh termscsusceptibility. The of the dashednals, line where at marksfixed the quasiparticlesTTcrosshas, where are fallen relatively the to nematicweakly half scat- of sus- its value reside on the sites, while the Ising pseudospins livek onF the C . This istent a particularly to which metallic relevant QCPssingle-fermion QCP are preempted given Green’s that function by ne- super- in Eq. (ceptibility14),T remainsN (h) sub- appears reachestered, 50% linear is apparent. of atits low magnitude Finally, temperature no superconducting at h = forhc transition both3.525 values at 2 x,y,z atTh = hc. Tcross(h), the inverse of hcross⇡, also vanishes bonds. Theconducting pseudospins [35– interact43] or other withstantial. forms their However, [neighbors⌧38, 44 it monotonically, 45] antifer- of auxil- decreasestheof same on fermion cooling,1.5 temperature. densityis found down is The consistent to pale our lowest (grey) with temperatures, lines this show being although the a corre-theuniver- with downwardi,j curvature. In a Fermi liquid,linearlyZ (T ) withsuperconductingh upon approaching susceptibility in the the s-wave QCP, channel although is its romagnetically,iary order, and are that coupled gap out to the the Fermi fermion surface, bond c) density. whether kF would approach a positive limit as T spondingsal0, whileproperty temperaturesin a peaked of the about for metallich theQuantumh case. (See QCP.↵ Fig.= 013 (where.) the fermions V slope is steeper than⇡ thatc of TN (h). there exists a “non-Fermiweak liquid” non-Fermi metal liquid, in the it would quantum vanishand! in proportion pseudospins to areThe decoupled). remainder ofcritical the In paper this is case, organized the as QCP follows: occurs in critical regime. Controlled theoretical methods, that can 1 These authors have contributed equally to this work. at hc0 3.06. ⇤ be used both to benchmark the field theories and for ⇡The linear behavior of TN (h) for small hc h is also comparison with experiments, are greatly needed. Deter- t seen0.5 for otherNematic model parameters. In Fig. 4 we show the minant quantum Monte Carlo (DQMC)[46, 47]–in par- phase diagram for two values ofSymmetric the fermion density, con- ticular, in cases where the fermion sign problem can trolled by the chemical potential µ. As the fermion den- be circumvented [48–50]–may serve this purpose. Such 0 methods have been successfully applied in several prob- sity is2 reduced,2.5 both3 hc 3.5hc0 and4 the range4.5 over which hc0 hc lems in which critical bosonic fluctuations are coupled to TN (h) is linear become smaller, indicatingh that the e↵ect fermions with a Dirac-like dispersion [51–53]. of the coupling between electrons near the Fermi surface FIG. 1: Phase diagram of the model obtained by DQMC, FIG. 3:arXiv:1511.03282v1 [cond-mat.supr-con] 10 Nov 2015 IllustrationIn this ofpaper, the lattice we report model a DQMC (3). Spinful study ofelectrons a two- asand a function the nematic of the transverse modes field becomesh and temperature weaker. TheT . fact that reside on thedimensional sites, while sign problem-free the Ising pseudospins lattice model that live exhibits on the Here,T (Vh)= appearst, ↵ =0 linear.5, µ = at0 low.5t. temperature The solid line marks for both values an “Ising nematic” QCP in a metal at finite fermion N bonds. The pseudospins interact with their neighbors antifer- theof transition fermion temperature, density isTN consistent,betweenthenematicandthe with this being a univer- romagnetically,density; and in are the coupled nematic to phase, the fermion the discrete bond lattice density. ro- symmetric phases. The line extrapolates to the T =0QCP tational symmetry is spontaneously broken from C4 to atsalhc.property The dashed of line the marks metallicTcross, where QCP. the nematic sus- C2. This is a particularly relevant QCP given that ne- ceptibility reaches 50% of its magnitude at h = hc 3.525 at the same temperature. The pale (grey) lines show⇡ the corre- sponding temperatures for the case ↵ = 0 (where the fermions and pseudospins are decoupled). In this case, the QCP occurs at hc0 3.06. ⇤These authors have contributed equally to this work. ⇡ Quantum Monte Carlo Methods for Fermion-Boson problems
Fakher F. Assaad (Autumn School on Correlated Electrons, Jülich 20th September 2018)
Outline: Ø Fermion-Boson problems Electron-Phonon, Su-Schrieffer-Heeger (SSH) Unconstrained Gauge theories “De-signer” Hamiltonians
Ø Auxiliary field Quantum Monte Carlo (QMC), Generalities
Ø Application of the Auxiliary field QMC to the SSH model Sign-free formulation Sampling, Hybrid Monte Carlo
Ø Conclusions
SFB1170 ToCoTronics Quantum Monte Carlo: Formulation