For the magnetic moment of the neutron in units of nuclear magnetons, we thus find µ conn n ( m$_{π}$ = 390 MeV ) = − 1 . 63(10)(4)(5) [ µ$_{N}$ ] . We have appended a superscript to reflect that our computation includes only connected contributions. The three uncertainties quoted are from: (i) statistics and fitting, (ii) the systematic due to the fit window, and (iii) conversion to units of physical nuclear magnetons. For (i), we take the largerst value of the uncertainty from the two fits to the field-strength dependence (which both gave the same value for µ $^{latt}$). There are additional sources of systematic uncertainty that we have not unaccounted for, namely the effects of finite lattice spacing and finite lattice volume. To convert the lattice electric polarizability to physical units, we compare the fit func- tion in lattice units, Eq. (21), to the energy in physical units, Eq. (3), and find where the uncertainty arises from scale setting, and is specifically three times the uncertainty in the lattice spacing. For the neutron electric polarizability, we thus find α$_{E}$ = e 2 2 π a$_{t}$a 2 s α latt E = 0 . 0776(58) × α latt E × 10 − 4 fm $^{3}$, (26) α n E $^{conn}$( m$_{π}$ = 390 MeV ) = 3 . 3(1 . 5)(2)(3) × 10 − 4 fm $^{3}$, taking the central value and uncertainties from fit II. The notation and sources of uncertainty on the electric polarizability are as for the magnetic moment. B. Proton For the proton, we perform a similar analysis. Fits to proton correlation functions are carried out using a two-state fit function. This allows us to remove excited state con- tamination. We perform simultaneous time-correlated fits to both boost projected proton correlation functions using the fit function G$_{±}$ ( t, n ) = Z ( n ) [ 1 ± ˜ µ $^{latt}$( n ) E latt ξ ] D ( t, E ( n ) 2 [ 1 − ( ˜ µ $^{latt}$( n ) E latt ξ ) $^{2}$] ∓ E latt ξ , E latt ξ ) + Z $^{′}$( n ) [ 1 ± ˜ µ ′ $^{latt}$( n ) E latt ξ ] D ( t, E $^{′}$( n ) 2 [ 1 − ( ˜ µ ′ $^{latt}$( n ) E latt ξ ) $^{2}$] ∓ E latt ξ , E latt ξ ) , (27) with D ( x, E $^{2}$, E ) as the relativistic propagator function given in Eq. (12). As the overall amplitudes Z ( n ) and Z $^{′}$( n ) enter the fit function linearly, we utilize variable projection to eliminate them from the simultaneous fits. For each value of the electric field E (or equivalently the integer n ), there are then four parameters in the fit: the ground state rest energy, E ( n ), the ground state anomalous magnetic coupling, ˜ µ ( n ), as well the rest energy netons, µ latt N = e 2 $_{M}$, with M as the lattice value of the nucleon mass. With these units, there is no uncertainty from scale setting, however, they introduce additional pion mass dependence of the extracted moment. 13