In the present paper, we study the properties of ϕ-meson longitudinal leading-twist light-cone distribution amplitude ϕ2;ϕ∥(x,μ) by starting from a light-cone harmonic oscillator model for its wave function. To fix the input parameters, we derive the first ten order ξ moments of ϕ2;ϕ∥(x,μ) by using the QCD sum rules approach under the background field theory. The curves of ϕ2;ϕ∥(x,μ=2 GeV) tend to be a single-peak behavior, which is consistent with the latest lattice QCD result. To show how the twist-3 light-cone distribution amplitudes (LCDAs) affect the results, we consider two scenarios for the ϕ-meson chiral twist-3 LCDAs ϕ3;ϕ⊥(x) and ψ3;ϕ⊥(x), i.e., the ones using the Wandzura-Wilczek approximation with ϕ2;ϕ∥(x,μ) (S1) and the ones using self-consistent conformal expansion with second-order Gegenbauer moments a2;ϕ2 in this work (S2). As an application, we derive the Ds+→ϕ transition form factors (TFFs) by using the QCD light-cone sum rules. The TFFs at large recoil point for those two scenarios are given separately. As for the two TFF ratios γV and γ2, we obtain γV(S1)=1.755−0.005+0.008, γ2(S1)=0.852−0.133+0.135, γV(S2)=1.723−0.021+0.023, and γ2(S2)=0.785−0.104+0.100. After extrapolating those TFFs to the physically allowable region, we then obtain the transverse, longitudinal, and total decay widths for semileptonic decay Ds+→ϕℓ+νℓ. Then the branching fractions are B(S1)(Ds+→ϕe+νe)=(2.347−0.191+0.342)×10−3, B(S1)(Ds+→ϕμ+νμ)=(2.330−0.190+0.341)×10−3, B(S2)(Ds+→ϕe+νe)=(2.367−0.132+0.256)×10−3, and B(S2)(Ds+→ϕμ+νμ)=(2.349−0.132+0.255)×10−3, which show good agreement with the data issued by the BESIII, CLEO, and Collaborations. We finally calculate Ds+→ϕℓ+νℓ polarization and asymmetry parameters, which can be measured and tested in future experiments. Published by the American Physical Society 2024