A main target of high energy physics is to test the standard model (SM)
and to search for new physics (NP) beyond the SM [1]}, [2]}, [3]} after the
discovery of the Higgs boson [4]}, [5]}.
Before the collision energy of the center of mass of the running collider
reaches the threshold energy of new physics, the SM theoretical
predictions of various physical observations should be accurately calculated [6]}.
In order to predict the observables precisely, one should evaluate Feynman integrals
exactly in the time-space dimension \(D=4-2\varepsilon \)  [7]} at first.
