In this article, we calculate the form factors and the coupling constant of the vertex D * s D s φ using the threepoint QCD sum rules. We consider the contributions of the vacuum condensates up to dimension 7 in the operator product expansion. And all possible off-shell cases are considered, φ, D s and D * s , resulting in three different form factors. Then we fit the form factors into analytical functions and extrapolate them into time-like regions, which giving the coupling constant for the process.