## Abstract

We present a classification and description of the principal resistive switching and transport mechanisms in chalcogonides materials. We classify the model according to how many material dimensions are involved in the resistive switching mechanism. In this way, we describe the phase change model (3D), the interface modulation model (2D) and models where the switching mechanism depends on the formation of a conduction filament (1D). Among the conduction filament models, we include the thermochemical oxygen diffusion mechanism, the oxidation/reduction mechanism and the quantum point effect.

### Keywords

- chalcogonides
- resistive switching
- physical models

## 1. Introduction

Typically, a resistive switching material changes its resistance between two states: High Resistive State (HRS or OFF‐state) and Low Resistive State (LRS or ON‐state). The most common structure for a resistive switching devices is an insulator between two metals or metal/insulator/metal (MIM) structure. One of the most important applications of these kind of devices is for non‐volatile memories or Resistive RAMs (ReRAMs).

The metal elements of the MIM structure are called top and bottom electrodes. An electrical stimulus is necessary to apply between these electrodes to change the resistive state of the insulator material. In order to determine the resistive state (HRS or LRS), a low voltage is applied on the electrodes, and the current, which flows through the insulator, is measured (*I*_{HRS} or *I*_{LRS}). There are several orders of magnitudes of difference between *I*_{HRS} and *I*_{LRS} currents. The change from the HRS to the LRS is called SET process and the change from the LRS to the HRS is called RESET process.

Depending on the voltage polarity applied on the electrodes, there are two schemes to change the resistive state: unipolar and bipolar. On the one hand, in the unipolar scheme, the resistive state change does not depend on the voltage polarity and there are two threshold voltages: one for the RESET process (*V*_{RESET}) and one for the SET process (*V*_{SET}) with the same polarity as we can see in Figure 1. On the other hand, in the bipolar scheme, *V*_{RESET} and *V*_{SET} have different polarities (Figure 2).

There is no universal theory or model which explains the electron conduction in the two resistive states and the SET and RESET processes because there are many factors that affect the switching behaviour, such as the type of the insulator material, fabrication process, nature of the dielectric breakdown, among others. However, Waser and Wutting proposed a classification based on the type of the resistive switching mechanisms, such as nanochemical materials, ferroelectric tunnelling, electrostatics effects, phase change mechanism, thermochemical mechanism, redox‐based effect, electrochemical effect, molecular switching effect and magnetoresistive effect [1].

Chalcogonides are one of the most used materials in the fabrication of resistive memory devices. The switching mechanisms related to these materials are: phase change memory effect, thermochemical memory effect, redox‐based memory effect and interface defect modulation. The phase change mechanism affects the complete volume of the insulator material and it is considered as a 3D mechanism. A 2D resistive switching mechanism is the modulation of the defect density at the metal/insulator interface. Finally, when the resistivity material depends on the formation of a conduction filament (CF), 1D mechanisms are involved in the resistive switching process. Thermochemical diffusion of oxygen, reduction/oxidation of the CF and quantum point contact effects are typical 1D mechanisms.

## 2. Phase change model

Heavier chalcogonides, such as tellurides and selenides, show different electrical and optical properties in their amorphous and crystalline phases (Figure 3). The resistive switching of these materials is unipolar. On the one hand, the amorphous phase of these materials has high resistivity and low reflectivity and, on the other hand, crystalline phase has low resistivity and high reflectivity [2]. These properties are being exploited in the development of optical storage products, such as compact disks (CD), digital versatile disks (DVD), high‐definition digital versatile disks (HD‐DVD) and Blu‐ray disks (BR). The first kind of materials used for optical storage were good glasses, such as Te‐based alloys like Te_{85}Ge_{13} doped with Sb, S and P. These materials show good electrical switching properties in the amorphous and crystalline phases but the crystallization time was in the order of microseconds, too high to be considered for optical storage. The second generation of these materials shows shorter crystallization time and good optical properties. Among these materials, we have GeTe, Ge_{11}Te_{80}Sn_{4}Au_{23}, GeTe‐Sb_{2}Te_{3}, GeBiTe and GeInSbTe. After, a third family of materials was later discovered which includes alloys of Sb_{2}Te doped with Ag, In and Ge [2].

There are two important temperature thresholds in these materials: the melting temperature (*Tm*) and the glass‐transition temperature (*Tg*) with *T*_{m} > *T*_{g}. The process to change the material phase from crystalline to amorphous and vice versa is as follows [2]:

To write a bit, a short high‐pulse laser or current is applied on the crystalline material to reach

*Tm*temperature.The material is cooling down rapidly with a rate higher than 10

^{9}K/s. In a very short time, the material reaches the amorphous phase without passing through the crystalline one.To erase the bit, a long short‐pulse laser or current is applied on the amorphous material. The material temperature increases over

*Tg*. There is an increment of the electron mobility and the material changes to the crystalline phase.

For electrical storage devices, the resistive switching property is fundamental and all these materials have it. However, not all these materials have the reflectivity switching property. To have this property, a very short time (few tens of nanoseconds) is needed for cooling down the material from the liquid phase to amorphous phase (step 2) [2].

When a low voltage is applied to the material in the amorphous phase, a very low current is measured due to the high resistance. When the voltage reaches a value around 0.7 V, the resistivity decreases and the material reach the so‐called ON‐amorphous phase [1, 2]. In these phases, the current increases significantly and enables enough heat to recrystallize the material. During the phase change, material defects play an important role. In the amorphous phase, the current is controlled by the Pool‐Frenkel conduction, where carriers are trapped in defect sites according to the following equation (electron hopping mechanism) [1]:

where *A* is the contact area, *V* the applied voltage, *N*_{T} the integral of the trap distribution, *Δz* the intertrap distance, *τ0* the scape time for a trapping electron, *E*_{F} the Fermi energy, *E*_{C} the conduction band energy, *q* the elementary electron charge and *t*_{h} the thickness of the material. Because of the total conduction is presented in the complete material volume, this mechanism is called 3D. Experimental results show that the defect density in the material is very high and most of the defects are negative U‐centres [1]. On the other hand, computational simulations show that in the crystalline phase, vacancy defects predominate with a concentration of 25% [1].

## 3. Metal/insulator interface modulation model

In this model, the resistive switching is presented at the metal/insulator interface. In other words, there is a contact resistance switching behaviour. This interface depending mechanism is presented in Peroskite oxides in which the material resistivity strongly depends on the interface area and the switching mechanism is always bipolar. A typical material that shows this behaviour is the Nb‐doped SrTiO_{3} [3].

The origin of this resistivity change can be understood by examining the metal/insulator interface band diagram as shown in Figure 4. The insulator oxide is usually doped with different metals. Depending of doped metal and its density, the insulator behaves as semiconductor at the interface. This provokes a Schottky barrier contact instead of a pure ohmic contact [3].

An electric field applied on the metal electrodes can electrochemically modify the oxygen vacancy density at the interface. For an n‐type semiconductor, an increment of the oxygen vacancies density reduces the depletion layer, *W*_{d}, in the energy band diagram provoking an increment of the tunnel electron conduction and, therefore, a decrement of the contact resistance (LRS). If an electric field is applied on the opposite direction, the number of oxygen vacancies in the n‐type semiconductor decreases and provokes an increment of *W*_{d} and, therefore, the contact resistance increases (HRS). On the other hand, for a p‐type semiconductor, the increment of oxygen vacancies increases *W*_{d} and the contact resistance (HRS), and a diminution of the oxygen vacancies decreases *W*_{d} and the contact resistance (LRS). In this model, the metal work function plays a very important role because the band bending strongly depends on this parameter.

Sawa showed very good resistive switching results for Ti/Pr_{07}Ca_{0.3}MnO_{3}/SrRuO_{3} (Ti/PCMO/SRO) cells, where SRO has metal properties and PCMO acts as p‐type semiconductor. As well, the SRO/SrTi_{0.99}Nb_{0.01}O_{3}/Ag (SRO/Nb:STO/Ag) cell showed a good resistive switching behaviour where Nb:STO acts as an n‐type semiconductor [3].

Another way to change the contact resistance is by adding a thin semiconductor layer of an oxide material between the metal and the insulator materials [3]. The semiconductor layer transforms the contact resistance from ohmic to a Schottky barrier. Without this layer, there is no switching resistance for some insulators. Sawa shows experimental results for Ti/Sm_{0.7}Ca_{0.3}MnO_{3}(n unit cells)/La_{0.7}Sr_{0.3}MnO_{3}/SRO (Ti/SCMO(n)/LSMO/SRO) where SCMO is a p‐type semiconductor [3]. It was demonstrated that for n = 5 unit cells, there was a very good hysteresis in the I‐V curve. As well, the SRO/SrTiO_{3}/SrTi_{0.99}Nb_{0.01}O_{3}/Ag (SRO/STO/Nb:STO) cell showed a resistivity changes but not a good hysteresis in the I‐V curve [3].

## 4. One dimensional models for resistive switching materials

In several transition metal oxides, when a voltage is applied on the electrodes of a pristine MIM cell, the current measured is very low. When the voltage increases up to a threshold value, the electric field applied provokes a dielectric breakdown. When this occurs, a conduction filament (CF) is formed in the insulator as shown in Figure 5. The necessary potential to form this filament is called forming voltage (*VF*). The CF is formed due to the Joule‐heating effect, which leads a temperature increment in the insulator. The dielectric breakdown is driven by a thermal runway. When a voltage is applied on a transition metal oxide, the resistance starts to decrease but not in a permanent way. In this process, the material temperature increases up to certain value (related to a threshold voltage) when a local redox reaction begins and causes structural defects (grain boundaries or dislocations) [4].

The voltage polarity to change states in cells with CF can be unipolar or bipolar. For unipolar behaviour, the electron conduction mechanism is related to thermochemical changes in the filament due to the Joule‐heating effect. On the other hand, for bipolar behaviour, the electron conduction through the CF depends on redox effects or quantum point contact effects. In both cases, the switching mechanisms (SET and RESET processes) are related to the thermochemical oxygen diffusion in the CF. We can consider the CF as a 1D parameter because the material resistivity is area independent. It is very important to note that in CF programmable devices, there exist a variability problem because not all filaments are equal or similar [5].

### 4.1. Thermochemical oxygen diffusion model

Ielmeni, Nardi and Cagli have been developed a physical model for NiO cells with very good concordance with experimental results [4]. In this material, the RESET process happens in small steps, whereas the change of states occurs suddenly during the SET process. The *V*_{RESET}and *I*_{RESET} parameters depend on the resistance material in the HRS. *I*_{RESET} always decreases when the resistance increases. On the other hand, the *V*_{RESET} in function of *R* curve has an *U* behaviour. That is, for low values of *R*, *V*_{RESET} decreases when *R* increases, and for high values of *R*, *V*_{RESET} increases when *R* increases. For the RESET process, the CF temperature depends on the square of the applied voltage as follows [4]:

where *T*_{0} is the room temperature and *R*_{th} is the effective thermal resistance. By using Eq. (2), we can obtain the *V*_{RESET} voltage:

where Δ*T*_{RESET} is the critical temperature increment for the onset oxidation. The ratio *R*/*R*_{th} is almost constant according to the Wiedemann‐Franz law for metals [4]. This means that *V*_{RESET} is almost constant and *I*_{RESET} decreases with respect to *R*. This is not true in experimental results. For explaining the increment of *V*_{RESET} with respect to high values of *R*, we have to study the size‐dependent Joule‐heating effect. The parameter *R*_{th} is the parallel of two resistances: *R*_{th}′ and *R*_{th}′′, where *R*_{th}′ only depends on the CF and *R*_{th}′′ depends on the rest of the material (bulk oxide). *R*_{th}′ can be computed by considering the thermal nanofilament conductivity, *k*_{th}, as follows [4]:

where *t*_{h} is the oxide thickness and *A*_{CF} is CF area. As *R*_{th} is inverse proportional to *A*_{CF}, for high values of *A*_{CF} (low resistance) *R*_{th} is approximately *R*_{th}′. On the other hand, for low values of ACF, we have that *R*_{th} ≈ *R*_{th}′′. As well, when *k*_{th} increases, *R*_{th}′ predominates over *R*_{th}′′.

For explaining the behaviour of the *V*_{RESET}‐*R* curve for low values of *R*, we have to consider that the ratio *R*/*R*_{th} is almost constant in Eq. (3) and Δ*T*_{RESET} must increase in order to obtain a metal diffusion in the filament, which is a filament area dependent process. Hence, *V*_{RESET} increases for low values of *R*. This size‐dependent diffusion is considering in following Arrhenius expression developed in Ref. [4]:

where *E*_{A} is the activation energy, *k* the Boltzmann constant, *t*_{0} and *ϕ*_{0} are constants, *ϕ* is the CF diameter and *t*_{RESET}= 1s is the reset time. Ielmini et al. showed very good results between their model and experimental data of NiO devices [4].

The electron conduction mechanism in the LRS strongly depends on the activation energy. For low values of *E*_{A}, the filament has a metallic behaviour and the resistance is given by:

where *T*_{0} is the room temperature, *R*_{0m} is the metallic resistant at *T*_{0} and *α* is the temperature coefficient. On the other hand, for high values of *E*_{A}, electron conduction is driven by the Pool‐Frenkel model in semiconductors and the resistance follows the following equation:

where *R*_{0s} is the extrapolated resistance at infinite *T*, *k* is the Boltzmann constant and *E*_{AC} is the activation energy for conduction. Both conduction behaviours are related to position of Fermi level (*E*_{F}). Inside the CF filament, there are oxygen vacancies, whereas the bulk oxide is doped by oxygen. An insulator doped by oxygen behaves as a p‐type semiconductor and, on the other hand, oxygen vacancies provoke an n‐type behaviour. Therefore, the conduction filament in LRS behaves as an n‐type semiconductor and the electron conduction is modulated by the concentration of oxygen vacancies, which is directly related to the *E*_{F} position. When the oxygen vacancies concentration is too high, the CF behaves as a degenerately doped semiconductor and *E*_{F} is very close or above the conduction band [4].

As was mentioned before, the SET process happens suddenly and it strongly depends on the resistance of the HRS (before the SET process). The HRS resistance can be described by the Pool‐Frenkel model and is given by Ielmini et al. [4]:

where τ*0* is the attempt‐to‐escape characteristic time for a carrier from a specific state, *t*_{h} is the material thickness (or filament length), *NT* is the dopant density, *A*_{CF} is the filament area and Δ*z* is the distance between positive charged defects. Ielmini et al. showed, from experimental data, that necessary power for the setting process (*P*_{SET}) is directly proportional to *R*^{−0.5}, which means that [4]:

The tendencies for *V*_{SET} and *I*_{SET} parameters in function of *R* were experimentally confirmed the Ielmini group in Ref. [4].

### 4.2. Trap‐assisted‐tunnelling model

In bipolar cells, the transition between the LRS and the HRS is commonly related to the formation and rupture of the CF. A typical material which presents a CF with bipolar behaviour is the HfO*2*. As well as in unipolar cells, oxygen vacancies play an important role. The most accepted theory for forming the CF in a virgin cell is that the oxygen atoms migrate from the CF to the insulator/metal interface due to the Joule‐heating effect. When the CF is already formed, to change from the LRS to the HRS, the CF is oxidized (oxygen atoms migrate from the electrode to the CF), whereas to change from the HRS to the LRS, the CF is reduced leaving oxygen vacancies and forming percolation paths (oxygen atoms migrate from the CF to the electrode).

Guan, Yu and Wong have developed a model for explaining the carrier conduction through the CF in bipolar cells where the principal transport mechanism is the trap‐assisted‐tunnelling (TAT). The continuity transport equation in the oxide region is given by Guan et al. [6]:

where *f*_{n} is the electron occupation probability of the *n*th trap, *R*_{mn} is the electron hopping rate from trap *m* to trap *n*, *R*_{n}*oL*/*R*_{n}*oR* are the electron hopping rate from trap *n* to the right/left electrode and *R*_{n}*iL*/*R*_{n}*iR* the electron hopping rate from the right/left electrode to trap *n*. It is well known that oxygen vacancies contribute to the TAT. In quasi‐steady state, Eq. (11) transform to [6]:

The current can be computed by evaluating the electron flow near the electrode:

The hopping rate can be computed by the Mott hopping model as [7]:

where *R*_{0} ≈ 10^{12}Hz is the vibration electron frequency, *r*_{mn} = |*r*_{m}−*r*_{n}| is the distance between vacancies *n* and *m*, *a*_{0} is the attenuation length wave function, *V*_{mn}*H*≈ −*FH*(*r*_{m}−*r*_{n}) ≈ *VH*(*r*_{m})−*VH*(*r*_{n}) is the barrier change due to an external electric field and *VH*(*r*_{n}) is the homogeneous component of the potential solution of the Poisson equation. The hopping rates from a trap to an electrode are [6]:

where e *R0*_{tunnel} is the tunnel coupling strength between a trap and an electrode, *NL,R* is the number of states at a given energy in an electrode and *E*_{v}*+*/*E*_{v}*•* are the energy of an empty/filled trap given by Guan et al. [6]:

*TnL,R,+•* is the tunnel probability from the left/right electrode into a trap given by the Wentzel‐Kramers‐Brilloin approximation [6]:

where *xn* is the *x*th component of *rn*, *L* is the oxide thickness and *m** is the tunnelling effective mass in the oxide. *F*_{in}*L,R* is the Fermi integral which represents the filled states in an electrode above *E*_{v}*+* and takes into account the inject electrons from the electrode into the trap *n*:

where *E*_{F}*L,R* is the Fermi level of the right/left electrode and *VL,R* is the applied voltage on the left/right electrode. On the other hand, *F*_{out}*L,R* is the Fermi integral which takes into account the number of empty states in an electrode bellow *E*_{v}*−* which can accept electrons from the trap *n*:

According to Guan et al. model, the generation oxygen vacancies are given by [7]:

where the time, *t*, is within the interval [*τ*, *τ* + *t*], *Feq* is the local electric field of an ion, 1/*t*_{0}≈ 10^{13}Hz is the oxygen vibration frequency, *E*_{a} ≈ 1eV is a parameter related to the height of the potential barrier and *γ* is a coefficient which represents the local enhancement due to the electric field. This rate dominates the SET process.

On the other hand, during the equilibrium state (absence of *Feq*), the oxygen vacancy recombination rate is given by Guan et al. [6]:

Therefore, the recombination rate for a non‐equilibrium state is [7]:

where *β* is a parameter related to concentration of oxygen ions which can be computed by the following approximation [6]:

where *L*_{p} is decaying length of ion concentration, *u*(*x*,*t*) is a function related to the oxygen diffusion and can be approximated by the complementary error function. *v* is velocity of the oxygen ions waveform given by Yu et al. [7]:

where *a* is the lattice constant, *E*_{m} is the migration barrier, *γdrift* is the enhancement coefficient related to the dielectric material and *F* is the electric field left by an oxygen ion. Eq. (13) is coupled with the solution of the Poisson equation to obtain the potential distribution in the cell:

where *ρ* is the volumetric charge density and *ε* is the material permittivity. The border conditions for Eq. (24) are: *V*(*x*=0)=*VL* and *V*(*x*=*L*)=*VR.* Guan et al. showed very god results of their model for experimental data of HfO*x* devices [7].

### 4.3. Quantum point contact model

The complete quantum point contact (QPC) model was developed by Miranda and Suñe [8]. Originally, the model was developed for explaining the soft and hard‐dielectric‐breakdown in SiO*2*. If the dimension of the narrowest point of the CF is in the order of the Fermi wavelength, *λ*_{F}, quantum point contact effects are presented. There are some experimental works, where the QPC model could explain well the transport conduction in the HRS and LRS for HfO*2* devices [9, 10, 11].

According to the QPC model, the first quantized sub‐band behaves as a potential barrier for the incoming electrons as shown in Figure 6. We used a parabolic potential as potential barrier with the following physical parameters: *Φ* being the potential barrier height measured at the Fermi level, *t*_{B} is the potential thickness at the Fermi level, *R* is a series resistance external to the constriction, *V* is the applied voltage on the electrodes, *q* is the elementary electron charge and *I* is the filament current that flows in the *x* direction.

The potential barrier height is defined by the cross‐sectional area of the constriction and determines two conduction states. For the HRS, the top of the potential barrier is above or inside the energy window and the dominant conduction mechanism is tunnelling (this description is valid only for low‐voltages). On the other hand, if the top of the potential barrier is below the energy window, the cell is in the LRS and the conduction mechanism is essentially ballistic (transmission probability close to 1). The conduction in the LRS is independent of the potential barrier. By assuming a parabolic potential barrier in the narrow constriction and by using the Landauer formalism for 1D quantum conductors and the zero‐temperature limit for the parabolic potential barrier, the current through the filament is [9, 10]:

where *T*(*E*) is the electron transmission probability, *N* is the number of active channels in the filament and *h* is the Planck’s constant. For a parabolic potential barrier, there is an analytical expression for *T*(*E*):

where *α* is a shape parameter related to *t*_{B}. By integrating over the total energy window, we have that the filament current is [9, 10]:

If we consider *V* ≫ *IR* and only one active filament for the HRS (not multiple filaments), we obtain [9, 10]:

For the LRS, we suppose an ideally ballistic transport that (*T*(*E*) ≈ 1) and the current is [9, 10]:

where *G*_{0}= 2*q*^{2}/*h* is the quantum conductance unit.

We show in Figure 7 experimental results of I‐V curves for the HRS and LRS and its fitting with the QPC model for HfO*2* cells. We have found in Ref. [10] that *α*∝*Φn* with *n* = −0.35. There is a QPC expression which relates *Φ* and α by considering *t*_{B} constant [9, 10]:

where *m** is the electron effective mass in the constriction. Moreover, the constriction radius, *r*_{B}, can be extracted by using another QPC equation [9, 10]:

where *z*_{0}= 2.404 is the first zero of the Bessel function *J*_{0}. In Ref. [10], we have experimentally found that *r*_{B}= 1.14 nm with a standard deviation of 0.06 nm for bipolar HfO_{2} cells, which agrees with the quantum approach of this model. Miranda et al. showed that *Φ* has a linear dependence with respect to the temperature given by Avellán et al. [12]:

where *Φ0* is a potential height at a *T*_{0} temperature (a given temperature) and *γ* is a temperature coefficient. For extracting *γ*, we can use the following expression [12]:

This dependence has been probed in experimental results of HfO*2* cells in Ref. [10].

## 5. Conclusions

We have presented a classification of physical models for explaining the resistive switching mechanisms in chalcogonides materials. In the literature, there are many physical models proposed for explaining the electron conduction and switching mechanism in specific materials and fabrication process conditions. In the present work, we divide the models according to the number of material dimensions involved on the resistive switching mechanism. The phase change mechanism (PCM) is presented in some Te‐alloys used an optical storage devices. In this switching mechanism, the material changes from the amorphous phase to the crystalline phase. Because of the resistivity affects the complete cell volume, the phase change mechanism is considered a 3D model. On the other hand, the modulation of the resistive contact is a 2D model because the defect concentration only affects the metal/insulator interface. This mechanism is presented in some perovskite materials.

The principal component of 1D models is the presence of a conduction filament. The filament is formed in a virgin cell by applying a certain threshold voltage. Depending on the cell polarity, the transport mechanisms can be: thermochemical diffusion of oxygen, filament oxidation‐reduction or quantum point contact. For unipolar cells (like NiO cells), the resistive switching and the carrier conduction are controlled by the thermochemical diffusion of oxygen in the CF due to the Joule‐heating effect. For bipolar cells (like HfO_{2}), the switching mechanism is related to the oxidation/reduction of CF. If the CF radius at the narrowest part of the CF is in the order of the Fermi wavelength, the transport is driven by the quantum point contact effect, otherwise the transport depends on the oxidation/reduction of CF.