Is anyone using MFT to calculate J , then which code are you using...is it possible by DFT .. IF yes then what will be the syntax for it , can u explain by giving a simple example.
Without knowing any info on the system it is hard to help regarding MFT, in particular.
But for DFT usually people try to employ, directly, a broken-symmetry approach on systems with two coupled spins. If the system has more than two magnetic atoms you'd need to solve a particular set of equations that can be obtained considering the electronic system, as you can read in some papers like Article A Few Comments on the Application of Density Functional Theo...
This kind of modelling is possible using Gaussian, Orca, well, any of the major packages that we use for DFT.