Semiempirical Two-Dimensional Model of the Bipolar Resistive Switching Process in Si-NCs/SiO₂ Multilayers

In this work, the SET and RESET processes of bipolar resistive switching memories with silicon nanocrystals (Si-NCs) embedded in an oxide matrix is simulated by a stochastic model. This model is based on the estimation of two-dimensional oxygen vacancy configurations and their relationship with the resistive state. The simulation data are compared with the experimental current-voltage data of Si-NCs/SiO2 multilayer-based memristor devices. Devices with 1 and 3 Si-NCs/SiO2 bilayers were analyzed. The Si-NCs are assumed as agglomerates of fixed oxygen vacancies, which promote the formation of conductive filaments (CFs) through the multilayer according to the simulations. In fact, an intermediate resistive state was observed in the forming process (experimental and simulated) of the 3-BL device, which is explained by the preferential generation of oxygen vacancies in the sites that form the complete CFs, through Si-NCs.