We have investigated unsubstituted and methyl substituted polyradical chains of meta-xylylene by using density functional theory-broken symmetry methodology (DFT-BS). Optimization of geometry in the high-spin and low-spin states have been done at B3LYP/6-31G(d,p) and M06-2X/6-31G(d,p) levels in unrestricted methodology. Single-point calculations on the high-spin optimized geometries have been done by using the 6-311G(d,p) basis set. Each polyradical has been found to be nonplanar with a high-spin ground state. Each has a coupling constant larger than thermal energy. For each group of polyradicals, the coupling constant has been found to exponentially decrease with increase in the number of phenylene groups. The B3LYP infinite chain limit has been estimated for both the unsubstituted and substituted polyradicals. The individual inter-radical-site coupling constants have been estimated for the triradicals and tetraradicals using HDVV Hamiltonian in ORCA 3.0.1 code. These are also generally large and positive, revealing a strong intersite ferromagnetic interaction. The intersite coupling constant too decreases with increasing distance between the radical centers. Finally, we have used CRYSTAL09 package for calculations on the infinitely long one-dimensional and periodic polyradical chains. The coupling constants estimated from the periodic calculations are quite large at about 500 cm(-1) and somewhat greater than the limiting values calculated for the polyradicals with an increasing number of phenylene groups. This happens as the individual polyradicals of increasing size progressively deviate from periodicity, and thereby lessens the strength of through-bond spin-spin coupling. The calculated band gap of similar to 4.5 eV indicates that the infinitely long one-dimensional chains must be ferromagnetic and electron insulators.