A new formula for the evaluation of the modal radiation Q factor is derived. The total Q of selected structures is to be calculated from the set of eigenmodes with associated eigen-energies and eigen-powers. Thanks to the analytical expression of these quantities, the procedure is highly accurate, respecting arbitrary current densities flowing along the radiating device. The electric field integral equation, Delaunay triangulation, method of moments, Rao-Wilton-Glisson basis function and the theory of characteristic modes constitute the underlying theoretical background. In terms of the modal radiation Q, all necessary relations are presented and the essential points of implementation are discussed. Calculation of the modal energies and Q factors enable us to study the effect of the radiating shape separately to the feeding. This approach can be very helpful in antenna design. A few examples are given, including a thin-strip dipole, two coupled dipoles a bowtie antenna and an electrically small meander folded dipole. Results are compared with prior estimates and some observations are discussed. Good agreement is observed for different methods.