We calculate the nuclear structure corrections to the Lamb shift in muonic deuterium by using state-of-the-art nucleon-nucleon potentials derived from chiral effective field theory. Our calculations complement previous theoretical work obtained from phenomenological potentials and the zero range approximation. The study of the chiral convergence order-by-order and the dependence on cutoff variations allows us to improve the estimates on the nuclear structure corrections and the theoretical uncertainty coming from nuclear potentials. This will enter the determination of the nuclear radius from ongoing muonic deuterium experiments at PSI.

Stimulated by the proton radius conundrum, measurements of the Lamb shift in various light muonic atoms are planned at PSI. The aim is to extract the rms charge radius with high precision, limited by the uncertainty in the nuclear polarization corrections. We present an ab initio calculation of the nuclear polarization for μ(4)He(+) leading to an energy correction in the 2S-2P transitions of δ(pol)(A)=-2.47 meV ±6%. We use two different state-of-the-art nuclear Hamiltonians and utilize the Lorentz integral transform with hyperspherical harmonics expansion as few-body methods. We take into account the leading multipole contributions, plus Coulomb, relativistic, and finite-nucleon-size corrections. Our main source of uncertainty is the nuclear Hamiltonian, which currently limits the attainable accuracy. Our predictions considerably reduce the uncertainty with respect to previous estimates and should be instrumental to the μ(4)He(+) experiment planned for 2013.

The measurement of the Lamb shift in muonic hydrogen and the subsequent emergence of the proton-radius puzzle have motivated an experimental campaign devoted to measuring the Lamb shift in other light muonic atoms, such as muonic deuterium and helium. For these systems it has been shown that twophoton exchange nuclear structure corrections are the largest source of uncertainty and consequently the bottleneck for exploiting the experimental precision to extract the nuclear charge radius. Utilizing techniques and methods developed to study electromagnetic reactions in light nuclei, recent calculations of nuclear structure corrections to the muonic Lamb shift have reached unprecedented precision, reducing the uncertainty with respect to previous estimates by a factor of 5 in certain cases. These results will be useful for shedding light on the nature of the proton-radius puzzle and other open questions pertaining to it. Here, we review and update calculations for muonic deuterium and tritium atoms, and for muonic helium-3 and helium-4 ions. We present a thorough derivation of the formalism and discuss the results in relation to other approaches where available. We also describe how to assess theoretical uncertainties, for which the language of chiral effective field theory furnishes a systematic approach that could be further exploited in the future.Keywords: two-photon exchange, muonic atoms, few-nucleon dynamics arXiv:1806.03101v1 [nucl-th] 8 Jun 2018 ‡ Using perturbation theory up to third order, Friar calculated in [26] the nuclear finite-size effect through order (Zα) 6 .

To shed light on the deuteron radius puzzle we analyze the theoretical uncertainties of the nuclear structure corrections to the Lamb shift in muonic deuterium. We find that the discrepancy between the calculated two-photon exchange correction and the corresponding experimentally inferred value by Pohl et al.[1] remain. The present result is consistent with our previous estimate, although the discrepancy is reduced from 2.6 σ to about 2 σ. The error analysis includes statistic as well as systematic uncertainties stemming from the use of nucleon-nucleon interactions derived from chiral effective field theory at various orders. We therefore conclude that nuclear theory uncertainty is more likely not the source of the discrepancy.

Energy-dependent sum rules are useful tools in many fields of physics. In nuclear physics, they typically involve an integration of the response function over the nuclear spectrum with a weight function composed of integer powers of the energy. More complicated weight functions are also encountered, e.g., in nuclear polarization corrections of atomic spectra. Using the Lorentz integral transform method and the Lanczos algorithm, we derive a computationally efficient technique for evaluating such sum rules that avoids the explicit calculation of both the continuum states and the response function itself. Our numerical results for electric dipole sum rules of the 4 He nucleus with various energy-dependent weights show rapid convergence with respect to the number of Lanczos steps. This demonstrates the usefulness of the method in a variety of electroweak reactions.