From the MD simulations, we estimated the average x-displacement of the CIT ions during the runs (Figure : ⟨x⟩). The data show how in all cases the systems reach a dynamic equilibrium; i.e., the diffusion along x direction reaching a constant rate during the MD (see Methods section for further details). While these data are extracted from an approximated CG model (where, e.g., the E is constant in intensity and direction, polariza- tion effects on the solvent molecules or on the Au NPs are not explicitly treated, etc.), and as such should be considered as qualitative, they provide interesting insights useful for understanding the global behavior of these superlattices. After a short initial transient phase, in all cases where CIT diffusion is activated by the applied E, the systems reach a steady state characterized by a constant ionic current. The plots in Figure -c also show how the extent of this ionic drift is proportional to the E intensity. From the ⟨x⟩ data over MD time it is possible to estimate the ionic current j = N qv d , where N is the number density of the CIT ions, q = -3e is the CIT charge, and v d is the average drift velocity (in x direction) measured in our simulations. Figure shows the absolute value of the j as a function of the intensity of the applied E. The j plot shows two different regimes: for E lower than ∼ 0.05 -0.07 V nm -1 (at T = 333 -300 K, respectively), the ionic current is relatively negligible and the superlattice is non-conductive, while for higher E intensities, the CIT current density increases quasi-linearly with the amplitude of E. These results indicate a non-ohmic behavior of such superlattice in terms of CIT conduction, in sharp contrast to pure metals, which are typically ohmic conductors. In Figure we also note that the ionic current j increases with increasing temperature, as thermal agitation favors CIT mobility.