Results of numerical simulations of fusion rate d(d,p)t, for low-energy deuteron beam, colliding with deuterated metallic matrix (Raiola [1,2]) confirm analytical estimates given in Ref.[3] (M. Coraddu et al., this issue), taking into account quantum tails in the momentum distribution function of target particles, and predict an enhanced astrophysical factor in the 1 keV region in qualitative agreement with experiments.