supported in part by the Grants-in-Aid for Scientific Research from the Ministry of Education, Science and Culture (Nos. 2373, 2375). T. I., Y. K. and Y. T. are supported by Japan Society for Promotion of Science.

Appendix A: Mean-field improvement {#appendix-a-mean-field-improvement .unnumbered}
==================================

The mean-field improvement[@MF] in our paper uses $$u= 1 - \frac{g^2 C_F}{2}T$$ with $T= 0.15493$, which is the value for the link in Feynman gauge. It may be better to use $u$ from $K_c$ or plaquette in DWQCD. In that case $$u = 1 - \frac{g^2 C_F}{2}(T + \delta T),$$ where $\delta T =0.00793$ for $K_c$, or $\delta T = -0.02993$ for plaquette. Accordingly we have to modify renormalization factors as follows: $$\begin{aligned}
z_w^{MF}(T+\delta T) &=& z_w^{MF}(T)+\frac{2w_0}{1-w_0^2}16\pi^2\times
2\delta T, \\
z_2^{MF}(T+\delta T) &=&z_2^{MF}(T)+16\pi^2\times \delta T /2, \\
z_{\Gamma,\rm bilinear}^{MF}(T+\delta T) &=& z_{\Gamma,\rm bilinear}^{MF}(T)
+16\pi^2 \times \delta T/2, \\
z_{\Gamma,\rm 4-quark}^{MF}(T+\delta T) &=& z_{\Gamma,\rm 4-quark}^{MF}(T)
+16\pi^2 \times\delta T/2 \times 2 C_F, \\
z_{\Gamma,\rm 3-quark}^{MF}(T+\delta T) &=& z_{\Gamma,\rm 3-quark}^{MF}(T)
+16\pi^2 \times\delta T/2 \times \frac{3}{2} C_F .\end{aligned}$$

Appendix B: Naive Dimensional Regularization(NDR) {#appendix-b-naive-dimensional-regularizationndr .unnumbered}
=================================================

In this appendix we compile the finite part of the renormalization constant in the ${{\overline {\rm MS}}}$ subtraction scheme with the Naive Dimensional Regularization: $$\begin{aligned}
