In this work, we calculate the sub-leading power contributions to the radiative leptonic D → γ ν decay. For the first time, we provide the analytic expressions of next-toleading power contributions and the error estimation associated with the power expansion of O(Λ QCD /m c ). In our calculation, we adopt two different models of the D-meson distribution amplitudes φ + D,I and φ + D,II . Within the framework of the QCD factorization as well as the dispersion relation, we evaluate the soft contribution up to the next-toleading logarithmic accuracy, and the higher-twist contribution from the two-particle and three-particle distribution amplitudes is also considered. Finally, we find that all the sub-leading power contributions are significant at λ D (µ 0 ) = 354 MeV, and the nextto-leading power contributions will lead to 143% in φ + D,I and 120% in φ + D,II corrections to leading power vector form factors with E γ = 0.5 GeV. As the corrections from the higher-twist and local sub-leading power contributions will be enhanced with the growing of the inverse moment, it is difficult to extract an appropriate inverse moment of the Dmeson distribution amplitude. The predicted branching fractions are (1.88 +0.36 −0.29 ) × 10 −5 for φ + D,I and (2.31 +0.65 −0.54 ) × 10 −5 for φ + D,II , respectively.