We present the self-consistent, non-perturbative analysis of isospin mixing using the nuclear density functional approach and the rediagonalization of the Coulomb interaction in the good-isospin basis. The largest isospin-breaking effects are predicted for N =Z nuclei and they quickly fall with the neutron excess. The unphysical isospin violation on the mean-field level, caused by the neutron excess, is eliminated by the proposed method. We find a significant dependence of the magnitude of isospin breaking on the parametrization of the nuclear interaction term. A rough correlation has been found between the isospin mixing parameter and the difference of proton and neutron rms radii. The theoretical framework described in this study is well suited to describe a variety of phenomena associated with isospin violation in nuclei, in particular the isospin symmetry-breaking corrections to superallowed Fermi beta decays. [2], is largely preserved by strong interactions; a small violation of isospin on the hadronic level is due to the difference in the masses of the up and down quarks [3]. In atomic nuclei, the main source of isospin breaking is the electromagnetic interaction [4,5]. Since the isovector and isotensor parts of electromagnetic force are much weaker than the strong interaction between nucleons, many effects associated with isospin breaking in nuclei be can treated in a perturbative way. With this caveat, the formalism of isotopic spin is a very powerful concept in nuclear structure and reactions [6,7], where many spectacular examples of isospin symmetry can be found.The main effect of Coulomb force in nuclei is to exert a long-range overall polarization effect on nuclear states whose detailed structure is dictated by the short-ranged strong force. The net effect of such a polarization is a result of two competing trends: the nuclear force is strongly attractive in the isoscalar neutron-proton channel, while the Coulomb force acts against this attraction by making neutron and proton states different. In order to explain this interplay, self-consistent feedback between strong and electromagnetic fields must be considered to best locate the point of the nuclear equilibrium.An excellent example of this interplay is the systematic behavior of nuclear binding energies: with increasing mass number, the stability line bends away from the N =Z line towards the neutron-rich nuclei. The effect of electromagnetic force on nuclear binding is clearly nonperturbative. Even in medium-mass nuclei, which are of principal interest in this study, energy balance between strong and Coulomb forces is not tremendously favorable, e.g., 342 MeV versus 72 MeV in 40 Ca. The situation becomes dramatic in superheavy nuclei and in the neutron star crust, where not only the binding but also spectra are strongly impacted by the Coulomb frustration effects resulting from a self-consistent, non-perturbative feedback between strong and electromagnetic parts of the nuclear Hamiltonian [8,9].The strong motivator for studies of isospin breakin...