Hadronic decay modes D0→(Kˉ0,Kˉ∗0)η,η′ and
(D+,Ds+)→(π+,ρ+)η,η′ are studied in the generalized
factorization approach. Form factors for (D,Ds+)→(η,η′) transitions
are carefully evaluated by taking into account the wave function normalization
of the eta and eta'. The predicted branching ratios are generally in agreement
with experiment except for D0→Kˉ0η′,D+→π+η and
Ds+→ρ+η′; the calculated decay rates for the first two decay modes
are too small by an order of magnitude. We show that the weak decays D0→K−π+ and D+→K+Kˉ0 followed by resonance-induced final-state
interactions (FSI), which are amenable technically, are able to enhance the
branching ratios of D0→Kˉ0η′ and D+→π+η dramatically
without affecting the agreement between theory and experiment for D0→Kˉ0η and D+→π+η′. We argue that it is difficult to understand
the observed large decay rates of Ds+→ρ+η′ and ρ+η
simultaneously; FSI, W-annihilation and the production of excess eta' from
gluons are not helpful in this regard. The large discrepancy between the
factorization hypothesis and experiment for the ratio of Ds+→ρ+η′
and Ds+→η′e+ν remains as an enigma.Comment: 15 pages, 1 figure, to appear in Phys. Rev. D. Form factors for D to
eta and eta' transitions are slightly change