This chapter discusses the methods of solving PEPS or other two-dimensional tensor network states, including variational optimization and the annealing simulation. The variation optimization determines the local terms by minimizing the ground-state energy. The annealing simulation takes the full or simple update strategy to filter out the ground state through the imaginary time evolution. The nonlinear effect arises in evaluating the derivative of uniform PEPS and is avoided by utilizing automatic differentiation. Both variational optimization and the annealing simulation involve a contraction of double-layer tensor network states. This contraction is the primary technical barrier in the study of PEPS. A nested tensor network approach is introduced to combat this difficulty.

This chapter introduces two commonly used methods of determining the local tensors of an MPS. The first is the variational optimization method, which determines an MPS by minimizing the energy expectation value. This method is equivalent to solving a generalized eigenequation around the extreme point of the ground-state energy. The second is an update method based on an imaginary time evolution, which cools down a quantum state from finite to zero temperature. We discuss three update approaches: update via canonicalization, full update, and simple update. For an MPS, the canonicalization approach is accurate and easy to implement. However, the full and simple update can be generalized to higher dimensions and applied to, for example, PEPS. The full update is a global minimization approach. It is accurate but has a higher computational cost than the simple update. The simple update is a local optimization approach based on an entanglement mean-field approximation and is easy to implement. Finally, we discuss the purification technique and apply it to evaluate the thermal density matrix or solve a quenched disorder problem in the framework of MPS.