where t is time, J is an electric current, E is an electric field, and κ is conductivity, i.e., a kinetic coefficient. After solving this equation, we conduct Gaussian fitting of ρ(r) by an algorithm that Guo et al. recently proposed (Fig. ). That is, we calculate (a, b, c) that minimizes an error function S as shown by