The Lakshmanan-Porsezian-Daniel equation (LPD) is a nonlinear Schrödinger type (NLS) equation which describe the effect on the integrable properties of Heisenberg bilinear spin chains under the classical limit by the biquadratic interactions. In this paper, the integrability of LPD equation is illustrated by obtaining the Lax pair via the AKNS scheme. The exact soliton solutions via the inverse scattering transformation (IST) method are obtained and illustrated in 3D graphical representation.