due to an exchange of soft Coulomb photons. Such contributions, however, were already taken into account within the non-relativistic calculations on a pointlike deuteron. Furthermore, the finite size effects were accounted for, as well, and have to be subtracted from the full result of Eq. (\[elastic\]) to avoid double-counting. This subtraction leads to &&|E\^[el]{}\_[n0]{}=\^2\_[n0]{}(0) \_0\^\[elastic\_subt\]\
&&{ G\_M\^2 (1+\_d)( -) .\
&&- (-)\
&&+16M\_d\^2G\_C’(0) }.We evaluate Eq. (\[elastic\_subt\]) with the most recent deuteron form factors’ parametrization from [@abbott]. We use the parametrization I and II of that Ref. to estimate the uncertainty, and list the result with the uncertainty in Table \[tab2\].\
The inelastic contributions contain two parts, E\^[inel]{}\_[n0]{}=E\^[QE]{}\_[n0]{}+E\^[hadr]{}\_[n0]{}, the quasielastic nucleon knock-out ([*QE*]{}) and hadronic excitation spectrum ([*hadr*]{}) that we will treat separately.\

Hadronic contribution
---------------------

The part of the deuteron excitation spectrum above the pion production threshold can be dealt with very similarly as it was done in Ref. [@Carlson:2011zd] for the proton case. We use the modern deuteron virtual photoabsorption data that were parametrized in terms of resonances plus non-resonant background by Bosted and Christy in [@bostedchristy]. Since the integration over the energy extend beyond the validity of the fit of Ref. [@bostedchristy], we supplement the correct high-energy 