Since the magnitude of the current, especially at low bias, see Eq 3, is mainly determined by coupling, number of molecules and coupling are thus convoluted. Therefore, in modeling, the accuracy of coupling largely influenced by the number of molecules used in the fitting. When explaining charge transport in molecular junction, coupling strength and number of molecules (N) in the junctions are like twin sisters that are hard to separate, especially when using Eq 3, where N and 𝛤 are all pre-factors of the expression. However, for Eq 1 and Eq 2, 𝛤 also influences the shape of the I-V in addition to the magnitude, to be more precisely the non-linearity of the I-V. Therefore, number of molecules N can be obtained, in principle, as an independent parameter, using Eq 1 and Eq 2 at large enough bias that reached the nonlinear part of the I-V. In this study, our I-V plots (figure ) did not consider the effect of number of molecules or channel i.e. number of molecules/channel assumed or considered to be unit. Therefore, the misinterpretation of fitting results in this study can originate from mis-estimation of the number of molecules. However, our map in Figure can still be used to examine whether proper parameters were obtained by using correct modeling methods since it was developed under idea single molecule condition.