The bipartite ground state entanglement in a finite linear harmonic chain of particles is numerically investigated. The particles are subjected to an external on-site periodic potential belonging to a family parametrized by the unit interval encompassing the sine-Gordon potential at both ends of the interval. Strong correspondences between the soliton entanglement entropy and the kink energy distribution profile as functions of the sub-chain length are found.