In this work we introduce a new family of ten-step linear multistep methods for the integration of orbital problems. The new methods are constructed by adopting a new methodology which improves the phase lag characteristics by vanishing both the phase lag function and its first derivatives at a specific frequency. The efficiency of the new family of methods is proved via error analysis and numerical applications.