The decay mode ψ ( 2 S ) → J / ψ + γ γ is proposed in order to experimentally identify the effects of the coupling of charmonium states to the continuum D ¯¯¯ D states. To have a better understanding of such a two-photon decay process, in this work we restrict ourselves to investigate the contribution of the discrete part, in which the photons are mainly produced via the intermediate states χ c J ( n P ) . Besides calculating the resonance contributions of χ c J ( 1 P ) ( J = 0 , 1 , 2 ) , we also take into account the contributions of the higher excited states χ c J ( 2 P ) and the interference effect among the 1 P and 2 P states. We find that the contribution of the 2 P states and the interference terms to the total decay width is very tiny. However, for specific regions of the Dalitz plot, off the resonance peaks, we find that these contributions are sizable and should also be accounted for. We also provide the photon spectrum and study the polarization of J / ψ .HadronPhysics2 (Grant Agreement No. 227431) (EU)