Abstract: The recent development in the field of two-dimensional magnetic materials urges reliable theoretical methodology for determination of magnetic properties. Among the available methods, ab initio four-state energy mapping based on density functional theory stands out as a powerful technique to calculate the magnetic exchange interaction in the Heisenberg spin model. Although the required formulas were explained in earlier works, the considered Hamiltonian in those studies always corresponded to the specific case that the off-diagonal part of J matrix is antisymmetric, which may be misleading in other cases. Therefore, using the most general form of the Heisenberg spin Hamiltonian, we here derive the generic formulas. With a proper choice of four different magnetic states, a single formula governs all elements of the exchange interaction matrix for any considered pair of spin sites.