Suppose, we bipartite the system in the middle in two equal half. In that case, is there any simple way to calculate EE like 1d system?
Hi, so yes you can calculate this kind of entanglement (I assume you mean using DMRG to first compute a ground-state MPS of a 2D system?).
The way to do it is to arrange the numbering of your sites to first visit every site in the left half of your system, then every site in the right half. If you do the "zig zag" path of an MPS through a 2D cylinder or rectangle this property will already be fulfilled.
Then simply compute the entanglement entropy the same way you would for a 1D system, cutting on the bond that separates the two halves of your system.
Here are some related questions and answers: