“…For the boundary parameters out of this region, stable boundary bound states exist in the ground state [41,42,43,44]. However, the energy is indeed a smooth function about the boundary parameters as demonstrated in the diagonal boundary field case [44,45]. (3)An interesting fact is that the contributions of a + , a − , β to the energy are completely separated and the surface energy does not depend on θ ± at all (same effect was also obtained in [46] where the surface energy and the finite size correction were derived for some constraint boundary parameters), which indicate that the two boundary fields behave independently in the thermodynamic limit.…”