Differential, total, and cumulative cross section calculations for neutral current neutrino scattering on 128,130 Te isotopes are performed in the context of the quasiparticle random phase approximation by utilizing realistic twonucleon forces. These isotopes are the main contents of detectors of ongoing experiments with multiple neutrino physics goals (COBRA and CUORE at Gran Sasso), including potential low-energy astrophysical neutrino (solar, supernova, geoneutrinos) detection. The incoming neutrino energy range adopted in our calculations (ε ν 100 MeV) covers the low-energy β-beam neutrinos and the pion-muon stopped neutrino beams existing or planned to be conducted at future neutron spallation sources. The aim of these facilities is to measure neutrino-nucleus cross sections at low and intermediate neutrino energies with the hope of shedding light on open problems in neutrino-induced reactions on nuclei and neutrino astrophysics. Such probes motivate theoretical studies on weak responses of various nuclear systems; thus the evaluated cross sections may be useful in this direction. V. TSAKSTARA AND T. S. KOSMAS PHYSICAL REVIEW C 83, 054612 (2011) mag(5) J (q) J i | 2The definitions of the eight multipole operators M (5)Obviously, for normal parity transitions, C V A contains contributions of T el J and T mag5 J operators while for abnormal parity ones it contains matrix elements of T mag J and T el5 J . The reduced one-body matrix elements of the eight multipole operators entering Eqs. (9) and (10) are evaluated through the use of the analytic formulas derived in Ref. [8] (see the Appendix).The magnitude of the three-momentum transfer q ≡ |q| and the square of the four-momentum transfer q 2 µ entering Eqs. (4)-(10) are written in terms of the kinematical parameters (laboratory scattering angle θ and lepton energies ε i and 054612-3 V. TSAKSTARA AND T. S. KOSMAS PHYSICAL REVIEW C 83, 054612 (2011)