1. Introduce quantum phase transition in the Bose-Hubbard model. 2. Emphasize vanishing energy scales and universality at the quantum critical point. 3. Introduce microscopic theories to describe the quantum phase transition, and explicitly show the vanishing of energy scales and the critical exponents. 4. Illustrate the emergent Lorentz symmetry and the Higgs mode at the quantum critical point of the Bose-Hubbard model with the particle-hole symmetry. 5. Discuss an experimental probe of the superfluid to the Mott insulator phase transition. 6. Show that the repulsive and the attractive Fermi-Hubbard models are related by the particle-hole symmetry. 7. Discuss the origin of antiferromagnetic order in the repulsive Fermi-Hubbard model at half-filling. 8. Introduce the enlarged
SO(4) symmetry of the Fermi-Hubbard model at half-filling and zero spin imbalance and its physical concequence. 9. Introduce important unsolved challenge issues in the Fermi-Hubbard model. 10. Introduce the concept of eigenstate thermalization hypothesis and many-body localization as opposite to thermalization. 11. Introduce a few metrics to characterize the many-body localization. 12. Emphasize the role of entanglement entropy in characterizing quantum thermalization, and discuss how to measure entanglement entropy in ultracold atom systems.
We adapt and optimize the projected-pair-entangled-state (PEPS) algorithm on finite lattices (fPEPS) for two-dimensional Hubbard models and apply the algorithm to the Hubbard model with nearest-neighbor hopping on a square lattice. In particular, we formulate the PEPS algorithm using projected entangled pair operators, incorporate SU(2) symmetry in all tensor indices, and optimize the PEPS using both iterative-diagonalization-based local bond optimization and gradient-based optimization of the PEPS. We discuss the performance and convergence of the algorithm for the Hubbard model on lattice sizes of up to $8\ifmmode\times\else\texttimes\fi{}8$ for PEPS states with U(1) symmetric bond dimensions of up to
D=8 and SU(2) symmetric bond dimensions of up to
D=6. Finally, we comment on the relative and overall efficiency of schemes for optimizing fPEPS.
Intertwining of spin, charge, and pairing correlations in the repulsive two-dimensional Hubbard model is shown through unrestricted variational calculations, with projected wave functions free of symmetry breaking. A crossover from incommensurate antiferromagnetism to stripe order naturally emerges in the hole-doped region when increasing the on-site coupling. Although effective attractive pairing interactions are identified, they are strongly fragmented in several modes including
d-wave pairing and more exotic channels related to an underlying stripe. We demonstrate that the entanglement of a mean-field wave function by symmetry restoration can largely account for interaction effects, and that our approach is exact for a two-site cluster.