We present a lattice calculation of the mass difference between neutron and proton, for which we find Mn−Mp=1.73(69)MeV. This is obtained at 1st order in the QED coupling αEM and in the mass difference between u and d quarks ΛQCDmd−mu. We adopt a purely hadronic scheme to renormalize the theory and provide a prescription to separate the QED and strong IB contributions. The simulation is carried out using the ETMC gauge configurations with Nf=2+1+1 dynamical quarks. We extrapolate among 3 values of the lattice spacing and pion masses in the range Mπ≃200−450 MeV.