Lattice QCD calculation of the two-photon exchange contribution to the muonic-hydrogen Lamb shift

We develop a method for lattice QCD calculation of the two-photon exchange (TPE) contribution to the muonic-hydrogen Lamb shift. To demonstrate the feasibility of this method, we also present an exploratory study with a gauge ensemble at $m_\pi = 142$ MeV. By adopting the infinite-volume reconstruction (IVR) method along with an optimized subtraction scheme, we obtain a preliminary result of the TPE contribution which agrees well with previous calculation using other methods and one magnitude smaller compare to the large $\sim300~\mu$eV discrepancy for the proton radius puzzle.