(Invited) Simulation and Modeling of the Switching Dynamics in Resistive Switching Devices
ECM cells consist of an ion-conducting layer sandwiched between a Ag/Cu active electrode and an inert counter electrode. The switching mechanism is based on the growth and dissolution of Ag/Cu filament within the ion-conducting host matrix involving the physical and electrochemical processes illustrated in Fig. 1. The VCM stack typically consists of a mixed ionic-electronic conducting layer, e.g. a metal oxide, a high work function metal electrode and an oxidizable counter electrode. The switching process involves the movement of oxygen vacancies within the MIEC and subsequent local redox reactions in the cation sublattice.
Fig.1: Illustration of the electrochemical processes occurring during SET switching in an ECM cell: charge-transfer reactions, nucleation, ion migration, filamentary growth. From  – Reproduced by permission of the PCCP Owner Societies
In order to fulfill the requirements of a nonvolatile memory a resistive switching device has to be switchable with a voltage amplitude of a few volts in the ns-regime and to offer a read-disturb immunity of up to ten years. This is equivalent to the driving physical processes exhibiting a nonlinearity of more than 15 orders of magnitude.To understand the underlying electrochemical and physical processes is hence of fundamental importance.
In the present paper, simulation models for ECM and VCM cells will be presented with special emphasis on the underlying physical and electrochemical processes. Based on simulation results the SET and RESET switching kinetics will be discussed in detail.
ECM compact models, continuum models and kinetic Monte Carlo (kMC) models[5, 6] will be presented. All ECM models include the electrochemical and physical processes illustrated in Fig. 1: nucleation, electron-transfer reactions at the boundaries, ion migration and filamentary growth. The models are capable of capturing the observed experimental phenomena: nonlinear switching kinetics, multilevel switching, current-voltage characteristics and switching variability.The multilevel switching is attributed to different tunneling gaps between the growing filament and the counter electrode at low current compliance and a filament widening if a galvanic contact is achieved. By comparison of the simulation results with experimental data the processes limiting the switching speed are discussed. Furthermore, the influence of mechanical stress on the filamentary growth will be addressed using kMC simulations.
The developed VCM simulation model is based on the solution of the drift-diffusion equation for oxygen vacancies and electrons coupled to Poisson’s equation and the heat transfer equation. Electron transport across the metal/oxide interface via thermionic emission and tunneling is also included within the model. Using this simulation model the nature of the gradual RESET and abrupt SET transition will be discussed. The gradual RESET can be explained as a result of a drift-diffusion balance. Moreover, the simulation results reveal that the switching kinetics in VCM devices are dominated by temperature-accelerated ion movement. Furthermore, the transition between bipolar and complementary switching can be explained by the variation of the metal/oxide interface properties.
This work was supported in part by the Deutsche Forschungsgemeinschaft (SFB 917).
 S. Menzel, S. Tappertzhofen, R. Waser, I. Valov, PCCP 2013, 15, 6945.
 R. Waser, R. Dittmann, G. Staikov, K. Szot, Adv. Mater. 2009, 21, 2632.
 S. Menzel, M. Salinga, U. Böttger, M. Wimmer, Adv. Func. Mater. 2015 (submitted)
 S. Menzel, N. Adler, J. van den Hurk, S. Tappertzhofen, I. Valov, R. Waser, 2013 5th IEEE International Memory Workshop (IMW) 2013, 92.
 S. Menzel, B. Wolf, S. Tappertzhofen, I. Valov, U. Böttger, R. Waser, 2014 6th IEEE International Memory Workshop (IMW), 2014
 S. Menzel, U. Böttger, R. Waser, J. Appl. Phys. 2012, 111, 014501/1.