Cohesive property of magnetized neutron star surfaces: Computations and implications