How can we find out the force of interaction between different ions in the silicon supercell from DFT outputs? For example, if we dope silicon bulk with some other element say A, with different ratios and different positions in the supercell, then what will be the force of interaction between two different ions and what will be net force of interaction between all silicon atoms and dopant (A) atoms?
Thank you.