Subsurface deformations in nematic liquid crystals: The hexagonal lattice approach