In resistive random access memory devices (ReRAM), in order to achieve high endurance and good retention characteristics repetitive control over the resistance change between low resistance (LRS) and high resistance (HRS) states is desired. The binary transition metal oxides TiOx, NiOx, HfOx, AlOx, TaOx have been specifically promising to achieve this goal and are the preferential candidates for the next generation of resistance change based devices. To address the critical issues of ReRAM operation, however, atomistic modeling approaches based on quantum mechanical principles are needed to achieve the optimization of material properties and structures of ReRAM cells, as well as of selection devices and memory array configurations.To explain the process of resistance change observed experimentally, various switching mechanisms were proposed recently and discussed in the technical literature. Generally, in transition metal oxide system such as TiO2, HfO2, NiO, Al2O3, Ta2O5, the diffusion of oxygen vacancies to cluster and create filaments [1-8] and the diffusion of oxygen atoms away from the oxide region to form a thin interfacial reduced oxide, are considered the key ionic processes taking place during forming and switching. Filamentary models for transition metal oxides had been proposed theoretically and the formation energy implications of a conductive filament channel formation corresponding to the “ON” state or LRS [4-9] was investigated for TiO2, HfO2 and Al2O3 (Fig. 1). It was found that during the clustering process of oxygen vacancies, conductive channels are formed by electron delocalization trends induced by the vacancies. The accompanying electronic effects of charge trapping under applied electrical field during the switching process were also recently addressed [7-9]. Hole injection into a reduced transition metal oxide containing a formed filament were found to favor the dissolution, since isolated vacancies prefer the 2+ charge state configuration. On the other hand electron injection induces filament formation and a transition of vacancy charge state to 1+or neutral configuration. The energy barriers of oxygen vacancy out-diffusion from a filament in various charge states are shown in Fig. 2 as a function of chemical potential and Fermi level control. In addition, preferential impurity doping [10] in these types of systems can favorably affect the forming of the conductive filaments and the transition process between the “ON” and “OFF” states (Fig. 3). In conclusion, several competing processes may undergo simultaneously at the microscopic level during the switching of ReRAM devices. Atomistic simulations based on density functional theory can provide an accurate description of the role of each component and describe the mechanism of ReRAM switching.[1] B. Magyari-Köpe, S.G. Park, H.-D. Lee and Y. Nishi, J. Mater. Sci. 47, pp 7498, 2012. [2] S.G. Park, B. Magyari-Köpe, and Y. Nishi, Phys. Rev. B, 82, 115109, 2010. [3] H.-D. Lee, B. Magyari-Köpe, and Y. Nishi, Phys. Rev. B 81, 193202, 2010. [4] S.G. Park, B. Magyari-Köpe, and Y. Nishi, Tech. Digest VLSI Tech. Symp., 2011. [5] B. Magyari-Köpe, M. Tendulkar, S.G. Park, H.-D. Lee, and Y. Nishi, Nanotechnology 22, 254029, 2011. [6] S.G. Park, B. Magyari-Köpe, and Y. Nishi, IEEE Electron Device Letters, 32, 197, 2011. [7] K. Kamiya, M.Y. Yang, S.G. Park, B. Magyari-Köpe, Y. Nishi, M. Niwa, and K. Shiraishi, Appl. Phys. Lett.,100, 073502, 2012. [8] K. Kamiya, M.Y. Yang, B. Magyari-Köpe, M. Niwa, Y. Nishi, and K. Shiraishi, IEEE Trans. Electron Devices, vol. 60, pp. 3400, 2013. [9] K. Kamiya, M.Y. Yang, B. Magyari-Köpe, M. Niwa, Y. Nishi, and K. Shiraishi, Proc. IEEE IEDM, pp. 478, 2012. [10] L. Zhao, S.G. Park, B. Magyari-Köpe and Y. Nishi, Tech. Digest VLSI Tech. Symp. , 2013.