Next Article in Journal
Research on Dynamic Response of Cold Rolling Mill with Dynamic Stiffness Compensation
Next Article in Special Issue
Memristors Based on Many-Layer Non-Stoichiometric Germanosilicate Glass Films
Previous Article in Journal
Efficient Neural Network DPD Architecture for Hybrid Beamforming mMIMO
Previous Article in Special Issue
Modulating the Filamentary-Based Resistive Switching Properties of HfO2 Memristive Devices by Adding Al2O3 Layers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Charge Transport Mechanism in the Forming-Free Memristor Based on PECVD Silicon Oxynitride

by
Andrei A. Gismatulin
1,*,
Gennadiy N. Kamaev
1,
Vladimir A. Volodin
1,2,* and
Vladimir A. Gritsenko
1,3
1
Rzhanov Institute of Semiconductor Physics SB RAS, Lavrentiev Ave. 13, 630090 Novosibirsk, Russia
2
Physical Department, Novosibirsk State University, Pirogova Street 2, 630090 Novosibirsk, Russia
3
Faculty of Automation and Computer Engineering, Novosibirsk State Technical University, 20 Prospekt K. Marksa, 630073 Novosibirsk, Russia
*
Authors to whom correspondence should be addressed.
Electronics 2023, 12(3), 598; https://doi.org/10.3390/electronics12030598
Submission received: 29 December 2022 / Revised: 23 January 2023 / Accepted: 23 January 2023 / Published: 25 January 2023
(This article belongs to the Special Issue RRAM Devices: Multilevel State Control and Applications)

Abstract

:
A memristor is a new generation memory that merges dynamic random access memory and flash properties. In addition, it can be used in neuromorphic electronics. The advantage of silicon oxynitride, as an active memristor layer, over other dielectrics it is compatibility with silicon technology. It is expected that SiNxOy-based memristors will combine the advantages of memristors based on nonstoichiometric silicon oxides and silicon nitrides. In the present work, the plasma-enhanced chemical vapor deposition (PECVD) method was used to fabricate a silicon oxynitride-based memristor. The memristor leakage currents determine its power consumption. To minimize the power consumption, it is required to study the charge transport mechanism in the memristor in the high-resistance state and low-resistance state. The charge transport mechanism in the PECVD silicon oxynitride-based memristor in high and low resistance states cannot be described by the Schottky effect, thermally assisted tunneling model, Frenkel effect model of Coulomb isolated trap ionization, Hill–Adachi model of overlapping Coulomb potentials, Makram–Ebeid and Lannoo model of multiphonon isolated trap ionization, Nasyrov–Gritsenko model of phonon-assisted tunneling between traps, or the Shklovskii–Efros percolation model. The charge transport in the forming-free PECVD SiO0.9N0.6-based memristor in high and low resistance states is described by the space charge limited current model. The trap parameters responsible for the charge transport in various memristor states are determined. For the high-resistance state, the trap ionization energy W is 0.35 eV, and the trap concentration Nd is 1.7 × 1019 cm−3; for the low-resistance state, the trap ionization energy W is 0.01 eV, and the trap concentration Nt is 4.6 × 1017 cm−3.

1. Introduction

Silicon oxynitride (SiNxOy) is a promising material, and it has very wide applications in different fields [1]. However, we will focus only on one of the silicon oxynitride applications as a memristor element [2]. A memristor is an element that can switch between high resistance state (HRS) and low resistance state (LRS) and save those states for 10 years at 85 °C [3,4]. SiNxOy is the mixture between silicon oxide (SiO2) and silicon nitride (Si3N4) [5,6]. SiO2 and Si3N4 are still the main materials used in silicon technology. Hence, the undoubted advantage of SiNxOy is that it is compatibility with standard technological processes of modern silicon microelectronics. SiNxOy is already used in logic devices as a thin layer (0.5–0.7 nm) at the interfaces with silicon when using high-k dielectrics to create a transistor gate. It is expected that memristors based on SiNxOy will combine the advantages of memristors based on nonstoichiometric silicon oxide (SiOx) and silicon nitride (SiNx), which are currently studied on a wide front [7,8,9,10]. According to our previous research data, the performance of SiOx-based memristors is higher than that of SiNx-based memristors [9,10]; the switching voltage in the SiOx-based memristors (+10 V and −10 V) is lower than in the SiNx-based memristors (+10 V and −16 V) with similar leakage current in the LRS and HRS. Therefore, the power consumption in the SiOx-based memristors is less than in the SiNx-based memristors. On the other hand, SiNx-based memristors are characterized by a longer storage time [9,10]. The SiNxOy-based memristors have a wide range of control of physical properties through the ability to vary the chemical composition, which opens up the possibility of controlling the memory properties. In addition, it can be expected that SiNxOy-based memristors will have a higher information storage time than SiOx-based memristors. This hypothesis is based on the fact that the information storage time is determined by the defect diffusion involved in the conductive filament forming (filament consists of nitrogen vacancies for SiNx-based memristors and of oxygen vacancies for SiOx-based memristors), and the nitrogen self-diffusion coefficient value (0.212 cm−2/s at 25 °C [11]) is lower than that of oxygen (0.232 cm−2/s [11]). The diffusion coefficient value of nitrogen in SiNx [12] and the diffusion coefficient value of oxygen in SiOx [13] are more difficult to compare due to their very strong diffusion coefficient dependence on temperature and the layer-forming technique [13]. Finally, SiNxOy-based memristors should have a higher performance than SiNx-based memristors due to the lower nitrogen vacancy formation energy compared to the silicon vacancy.
The memristor is very sensitive to the obtaining by the active layer technology and contacts to the active area. As an example, the GeO[SiO2]-based memristor has a better memory window (LRS/HRS ratio) than the GeO[SiO]-based memristor with the same contacts [14]. Silicon oxynitride films can be obtained by various methods: magnetron deposition [15], low-pressure chemical vapor deposition [16] and plasma-enhanced chemical vapor deposition (PECVD) [17]. The advantage of the PECVD method for the silicon oxynitride film is low deposition temperature that can be used to easier integrate the PECVD silicon oxynitride film in the formation of crossbar structures.
Usually, the charge transport in silicon oxynitride layers is described by the Frenkel effect [18,19]. Recently, it has been established that the model of phonon-assisted tunneling between neighboring traps describes the charge transport in the chemical vapor deposition of SiNxOy [20]. We studied the charge transport mechanism in detail in forming-free PECVD SiOx [9] and forming-free PECVD SiNx [10]. We also studied in detail the electronic structure and optical properties of different composition PECVD SiNxOy films for memristor applications [21]. Few papers study the charge transport mechanism in depth. In [16], SiNxOy film was obtained by the oxygen plasma treatment of LPCVD silicon nitride. As a result, the authors obtained the memristor with a SiNxOy layer and a silicon nitride layer. In this work, they explained the charge transport mechanism by the space-charge limited current (SCLC) model, but the modelling of experimental data was made only for the Si3N4-based memristor in the HRS. Thus, the aim of this work is to determine the charge transport mechanism in the PECVD SiNxOy memristor and estimate the trap concentration and trap ionization energy.

2. Materials and Methods

SiNxOy:H films were obtained by PECVD with a gas mixture of SiH4, N2, and O2 on p++-type Si (100) wafers after the RCA-1 cleaning procedure followed by oxide removal in diluted hydrofluoric acid, a final rinsing in deionized water, and drying in a nitrogen stream. However, after all the procedures, the loading of the substrate into the PECVD chamber was carried out through air; accordingly, the silicon was covered with a fresh natural oxide layer. The PECVD unit has a wide-aperture inductive plasma source with frequency 13.56 MHz and a turbomolecular pump. The residual pressure in the reactor was less than 10−6 Torr. The substrate temperature was set at 200 °C and maintained during the deposition process automatically by a heater. The SiH4 flow was constantly supplied to the reaction zone (gas mixture of 10% SiH4 diluted with Ar) and its supply rate was 10 sccm. The films were obtained at the nitrogen flow rate value of 6 sccm. This nitrogen flow value range is optimal (as it was found earlier when studying the PECVD SiNxOy films). If we applied a lower nitrogen flow, the obtained samples were not suitable for the memristor application. If we applied a higher nitrogen flow, the nitrogen content in plasma increased very slowly. Nitrogen gas was diluted with oxygen so that the oxygen concentration in the mixture was less than 1%. The low oxygen content in the gas mixture is needed to obtain the films highly enriched with silicon because oxygen has a very high chemical activity with silicon. In addition, the high oxidation rate is expected because of the presence of Ar [22]. The film was formed at the 6 sccm N2 flow and plasma generator power 150 W. The silicon oxynitride films with a thickness of about 33 nm were obtained with this PECVD method. The film composition was studied using the method of X-ray photoelectron spectroscopy (XPS). Based on the intensity of the peaks of silicon, oxygen, and nitrogen in the XPS spectra, the composition was determined as SiO0.9N0.6. To obtain the metal–dielectric–semiconductor memristor structure, the nickel electrodes were deposited through a metal mask using magnetron sputtering in the Ar atmosphere. The obtained nickel electrodes had a 200 nm thickness and a 0.5 mm2 area. A continuous 200 nm thick nickel layer was deposited on the p++-type back substrate side.
The current–voltage measurements were carried out in the temperature range of 300 to 400 K in the LTS420E Linkam cell by using the Keithley 2400 and temperature controller Linkam T95 devices. The voltage ramp rate used for the current–voltage measurements was 0.9 V/s.
To clarify the charge transport mechanism, the temperature current–voltage (IV) characteristics were taken for the PECVD SiO0.9N0.6-based memristor in the HRS and LRS. These characteristics were compared with two contact-limited models, namely, the Schottky effect model [23] and the thermally assisted tunneling (TAT) mechanisms [24,25], and six bulk-limited models, namely, the Frenkel effect model of Coulomb isolated trap ionization [26,27], the Hill–Adachi (H-A) model of overlapping Coulomb potentials [28,29], the Makram–Ebeid and Lannoo (ME-L) model of multiphonon isolated trap ionization [30], the Nasyrov–Gritsenko (N-G) model of phonon-assisted tunneling between traps [31], the Shklovskii–Efros (S-E) percolation model [32,33], and the SCLC model [34,35,36,37].
In the Schottky effect, the current is limited by the thermal electron emission through the barrier, which is reduced due to image forces (Figure 1a) [23]. The Schottky effect current is described by the following equation:
I = A S T 2 exp ( W 0 ( e 3 4 π ε ε 0 ) 1 / 2 U d k T ) ,
A = 4 π m * k 2 e m e h 3 = 120.4 m * m e [ A c m 2 K 2 ]
where I—current, A—Richardson–Deshman constant, S—contact area, T—temperature, W0—potential barrier height at the metal/SiOxNy interface, U—voltage, d—dielectric thickness, k—Boltzmann constant, e—electron charge, ε = n2—high-frequency dielectric permittivity, n—refractive index, ε0—dielectric constant, m*—electron effective mass, meelectron mass, h—Plank constant.
The TAT model describes the charge transport mechanism, where an electron, after the phonon absorption, is excited to a certain energy value and then tunnels through a triangle barrier (Figure 1b) [24,25]. The current for the TAT model has the following formula [24,25]:
I = 2 π m * k T ( e h ) 2 S U d exp ( W 0 1 6 ( h e × ( U d ) 4 π k T m * ) 2 k T ) ,
The current I through the material-containing traps can be described by the following equation [20]:
I = e N 2 / 3 S P
where N = a−3—trap concentration, a—average distance between traps, S—contact area, and P—trap ionization probability.
The Frenkel effect is due to a lower Coulomb positively charged trap potential in a strong electric field (Figure 2a). The Coulomb trap ionization probability has the following form [26,27]:
P = ν exp ( W ( e 3 π ε ε 0 ) 1 / 2 U d k T ) ,
where ν = W/h—attempt to escape factor, W—trap ionization energy.
The H-A model describes the charge transport mechanism though overlapping Coulomb traps. In the H-A model, the trap ionization probability is described by the following equation (Figure 2b) [28,29]:
P = 2 ν exp ( W e 2 π ε ε 0 s k T ) sinh ( e U d s 2 k T )
According to the ME-L model, the charge transport in a dielectric is controlled by the multiphonon isolated trap ionization (Figure 2c) [30]. The current density is described by Equation (3), and the trap ionization probability in the ME-L model is described by Equations (6) and (7) [30]:
P = exp ( n W p h 2 k T W o p t W t W p h coth n W p h 2 k T ) I n ( W o p t W t W p h sinh ( W p h / 2 k T ) ) P i ( W t + n W p h )
P i = e × ( U d ) 2 2 m * W exp ( 4 3 2 m * e × ( U d ) W 3 / 2 )
where Wt—thermal trap energy, Wopt—optical trap energy, Wph—phonon energy, Inmodified Bessel function, Pitunneling probability through a triangular barrier, and ћ = h/2π—reduced Plank constant.
The electron is excited from the ground state with the multiphonon absorption, and it tunnels to a neighboring trap due to a large integral overlapping value in the N-G model (Figure 2d). In the N-G model, the trap ionization probability is given by the following equation [31]:
P = 2 π W t m * s 2 2 k T ( W o p t W t ) exp ( W o p t W t k T ) exp ( 2 s 2 m * W t ) sinh ( e U d s 2 k T )
The charge transport in a non-stoichiometric dielectric can be described according to the S-E percolation model [32,33]. This model assumes that excited electrons with an energy higher than the percolation energy We are delocalized, and, moving round a random potential, they transfer the charge. The current–voltage characteristics have the following form [32,33]:
I = I 0 exp ( W e ( C e U d a V 0 γ ) 1 1 + γ k T ) ,
where I0 is the pre-exponential factor, We is the percolation energy, a is the space scale of fluctuations, V0 is the energy fluctuation amplitude, C ≈ 0.25 is a numeric constant, and ɣ is the critical index and is equal to 0.9.
The theory of SCLC between flat parallel electrodes is presented in detail in [37]. The current in the SCLC model consists of several sections with different current–voltage dependences. In our case, there are Ohmic and quadratic parts, which are given by the following expressions:
I = I Ohm + I SCLC = S e μ n U d + S 9 8 μ ε ε 0 θ U 2 d 3 ,
n = 2 N d 1 + 1 + 4 g N d N c exp ( E a k T ) ; N c = 2 ( 2 π m * k T h 2 ) 3 2 ; θ = N c N t exp ( W t k T ) ;
where S is the area involved in the charge transport, μ is electron mobility, n is free electron concentration in the dielectric, ε is static dielectric constant, θ is the fraction of free electrons from all injected ones (trapped and free), Nd is donor concentration, g is degeneracy factor, Nc is effective density of states, Ea is donor activation energy, Nt is trap concentration, Wt is trap energy, m* is electron effective mass, and h is the Planck constant.
If we have the third part in the SCLC model not proportional to the square, we will use the following SCLC formula [37]:
I = S e μ n U d + S 9 8 μ ε ε 0 θ U 2 d 3 + S N c μ e 1 l ( ε l N c ( l + 1 ) ) l ( 2 l + 1 l + 1 ) l + 1 U l + 1 d 2 l + 1
where l = Tt/T, and Tt is the temperature parameter that characterizes the exponential trap distribution.
To compare the theoretical charge transport models and the temperature IV characteristics, we used the method of least absolute deviations [38]. Most theoretical models are described in the logarithmic form (Log J~U or J~eU); in addition, most IV characteristics are also described in the form of Log J~U. Therefore, we used the absolute deviation as per the following formula:
Δ max = | Log J theory J experiment | × 100 %
The method of least absolute deviations was applied until the absolute deviation of the theoretical maximum current value from the experiment was 20% or less.
To model the experimental data, we used only the analytical calculation using Mathematica software. To conduct a more accurate simulation, we needed to account for the local field distribution in the sample with the Poisson equation and Shockley–Read–Hall statistics for the trap occupation. However, that required more specialized and complex software for numerical calculation.

3. Results and Discussion

The IV characteristics of the obtained PECVD SiO0.9N0.6-based memristor structures and their switching were investigated. The switching cycles for a PECVD SiO0.9N0.6-based memristor are shown in Figure 3. There were only two states: HRS and LRS. The PECVD SiO0.9N0.6-based memristor in the initial state (before forming) was in the LRS. Hence, the PECVD SiNxOy-based memristor did not require forming, which is an advantage, compared to memristors based on PECVD SiOx:H [9] and PECVD SiNx:H [10]. The memristors based on PECVD SiOx:H [9] and PECVD SiNx:H [10] have a unique forming process: the forming process does not require high voltages for the first switch, and the initial state in those memristors has a smaller resistance than that of LRS; however, after the first cycle, we cannot return to the initial state. Therefore, the memristors based on PECVD SiOx:H and PECVD SiNx:H need a forming stage, but they do not require high voltages for the forming cycle, which avoids the memory cell breakdown.

3.1. Contact-Limited Charge Transport Mechanism in the PECVD SiO0.9N0.6-Based Memristor

The IV characteristics of the PECVD SiO0.9N0.6-based memristor in the HRS and LRS were simulated by contact-limited charge transport models. The Schottky effect and thermally assisted tunneling (TAT) model qualitatively were used to describe the IV characteristics of the PECVD SiNxOy-based memristor in the LRS (Figure 4a) and HRS (Figure 4c). However, the Schottky effect predicted an abnormally small effective mass value for the PECVD SiNxOy-based memristor (for the LRS m* = 2.3 × 10−5 me and for the HRS m* = 1 × 10−12 me). The Schottky effect also gave an abnormally larger high-frequency dielectric constant value for the PECVD SiNxOy-based memristor (for the LRS ε = 18 and for the HRS ε = 9). The TAT model did not describe the charge transport mechanism because when modelling the experimental data, it obtained abnormally small values for the transport area for the PECVD SiNxOy-based memristor (for the LRS S = 4 × 10−16 m2 and for the HRS S = 1.0 × 10−14 m2). In addition, the barrier energies for the Schottky effect and TAT model were 0.15 eV for the LRS and HRS. Those energies represent the barrier energies at the Si/SiO0.9N0.6 interface. In [20], we estimated the energy gap as 4.1 eV for this memristor structure. Thus, the barrier energies at the Si/SiO0.9N0.6 interface for an electron and a hole, at least, will be around 1 eV or more in the HRS.

3.2. Bulk-Limited Charge Transport Mechanism in the PECVD SiO0.9N0.6-Based Memristor

In the memristor, the HRS conduction filament is partially or fully disrupted. In that case, conduction in the HRS is similar to dielectric conduction. Thin dielectric film conductivity can be divided into two processes: charge-carrier injection from the contact into the dielectric, and charge transport in the dielectric bulk. Because current depends on temperature, we excluded the Fowler–Nordheim model [39]. If there is a high trap concentration in the dielectric, then the space charge on them limits the injection from the contact, and the leakage current is limited by the trap ionization rate. The trap ionization can be described by six analytical bulk-limited current models.
The experimental charge transport simulation data with the Frenkel effect for the PECVD SiO0.9N0.6-based memristor in the LRS and in HRS are shown in Figure 5a and Figure 6a. From the current temperature dependence, we found that the trap ionization energy was 0.13 eV. The attempt to escape factor controls the vertical position of the theoretical IV characteristic, and it can be found from the trap ionization energy with formula ν = W/h = 3.1 × 1013 s−1. The trap concentration parameter in the Frenkel model also governs the vertical theoretical curve position, and the high-frequency dielectric constant governs the slope of the theoretical curve. To describe the experimental data by the Frenkel effect, we obtained the abnormally low trap concentration value for the PECVD SiO0.9N0.6-based memristor (for the LRS N = 8 × 106 cm−3 and for the HRS N = 100 cm−3) and the abnormally large high-frequency dielectric constant value for the PECVD SiO0.9N0.6-based memristor in both states (for the LRS ε = 65 and for the HRS ε = 28).
The IV characteristics of the PECVD SiO0.9N0.6-based memristor and the results of the simulation by the H-A model are shown in Figure 5b and Figure 6b. The H-A model describes the charge transport when Coulomb traps overlap; the overlapping happens when the distance between traps is small. Therefore, the trap concentration is large. Hence, we could use the H-A model to obtain an accurate estimate of the trap concentration value (for the LRS—N = 3.5 × 1020 cm−3 and for the HRS—N = 3.5 × 1020 cm−3). The attempt to escape factor could be found with the same formula as in the Frenkel model. Moreover, with the trap ionization energies W = 0.22 eV for the LRS and W = 0.35 eV for the HRS, the theoretical attempt to escape factor values were ν = 5.3 × 1013 s−1 for the LRS and ν = 8.5 × 1013 s−1 for the HRS. However, to describe the experimental data, the H-A model predicts the abnormally small attempt to escape factor values (for the LRS—ν = 8 × 103 s−1 and for the HRS—ν = 1.2 × 101 s−1) for the PECVD SiO0.9N0.6-based memristor in both states. In addition, the abnormally large high-frequency dielectric constant value for the PECVD SiO0.9N0.6-based memristor in both states is predicted by the H-A model (for the LRS—ε = 60 and for the HRS—ε = 20). Thus, the H-A model does not describe the charge transport mechanism in the PECVD SiO0.9N0.6-based memristor.
The ME-L model describes the multiphonon ionization of isolated traps. When we modeled the experimental data with this model, we obtained reasonable thermal trap ionization energy, optical trap ionization energy, and effective mass values (Figure 5c and Figure 6c). However, in the ME-L model, the trap concentration value was predicted to be abnormally small (for the LRS—N = 1.2 × 109 cm−3 and for the HRS—N = 1.0 × 103 cm−3). In addition, for the LRS, even with those reasonable parameters, the ME-L model described the experimental data only at high voltages. Therefore, the ME-L model did not describe the charge transport mechanism in the PECVD SiO0.9N0.6-based memristor in any states.
The IV characteristics of the PECVD SiO0.9N0.6-based memristor and results of simulation by the N-G model are shown in Figure 5d and Figure 6d. The N-G model is similar to H-A model and describes the charge transport mechanism at a high trap concentration, and the N-G model predicts the same trap concentration as in the H-A model (for the LRS—N = 3.5 × 1020 cm−3 and for the HRS—N = 3.5 × 1020 cm−3). It is the same value because the N-G model and the H-A model govern the same theoretical slope and vertical position as in another model. However, the N-G model did not describe the charge transport mechanism, since quantitative agreement of the experimental data and the theoretical models was reached with the abnormally large effective mass value for the PECVD SiO0.9N0.6-based memristor (for the LRS—m*= 47 me and for the HRS—m*= 9.2 me).
For the S-E percolation model, the resulting fluctuation space scale was about 0.6 nm for the LRS and 1.4 nm for the HRS for the PECVD SiO0.9N0.6-based memristor (Figure 5e and Figure 6e). Since the tunneling probability at such spatial scales was high, the classical approximation criterion of the S-E model applicability was not fulfilled. Thus, the S-E percolation model could not be used to describe the charge transport mechanism in the PECVD SiO0.9N0.6-based memristor in any states.
To apply the SCLC model, we did not account for the switching process. In the classical SCLC, an increase in current is related to filling traps, but in the memristor switching process, it is very complex and is usually related to creating and destroying a conduction channel. Thus, the SCLC cannot be applied to bipolar memristors. The SCLC can only be applied individually to the LRS or HRS without the switching process. To account for the differences in memristor states, we introduced the filament size changes from state to state. We did not account for the conductive filament thickness changes, although in the HRS the conductive channel can dissolve partially or completely. In addition, we did not take into account the possibility of structural changes in the conductive channel during the resistance switching. With those assumptions, the current in the LRS and in HRS for the PECVD SiO0.9N0.6-based memristor was described by the SCLC model (Figure 5f and Figure 6f). There are many variable parameters in the SCLC model. The effective mass was taken as m* = 0.5 me, and the static dielectric constant value was taken as 6 for SiO0.9N0.6. The simulation of the experimental IV characteristics in the LRS for the PECVD SiO0.9N0.6-based memristor showed that the effective conducting channel diameter of was 820 nm, the donor defect ionization energy was 0.11 eV, the donor defect concentration was 5.5 × 1019 cm−3, the trap ionization energy was 0.01 eV, and the trap concentration was 4.6 × 1017 cm−3. The simulation of the experimental IV characteristics in the HRS for the PECVD SiO0.9N0.6-based memristor gave the following results: the conductivity spanned the entire contact area (S = 0.005 cm2), the donor defect ionization energy was 0.19 eV, the donor defect concentration was 4.6 × 1017 cm−3, the trap ionization energy was 0.35 eV, and the trap concentration was 1.7 × 1019 cm−3. In the HRS, there was a third section of SCLC with a degree of 2.3. The Table A1 and Table A2 in Appendix A show all fitted parameters in all applied models.
The SCLC characteristics were significantly different from those of the experimentally observed in the HRS at 400 K from the 3 V to the 4 V part of the curve. This difference was due to the pre-switching portion of the IV curve, which was similar to the IV curve in Figure 3 from 6 V to 8 V. At a higher temperature, the switching began at a low voltage, as shown in Figure 6f. We simulated the IV curves only at a fixed resistance. However, at this pre-switching portion, the resistance began to change. Thus, we did not include it in the simulation.

4. Conclusions

The IV characteristics of the obtained PECVD SiO0.9N0.6-based memristor structures and their switching revealed that the PECVD SiO0.9N0.6-based memristor is forming-free. To clarify the charge transport mechanism in the PECVD SiO0.9N0.6-based forming-free memristor, the current at different temperature values in the high resistance state and low resistance states was simulated by the Schottky effect, the thermally assisted tunneling model, the Frenkel effect model of Coulomb isolated trap ionization, the Hill–Adachi model of overlapping Coulomb potentials, the Makram–Ebeid and Lannoo multiphonon isolated trap ionization model, the Nasyrov–Gritsenko model of phonon-assisted tunneling between traps, the Shklovskii–Efros percolation model, and the space-charge limited current model. After filtering out the models by various fitted parameters, the charge transport mechanism in the forming-free PECVD SiOxNy-based memristor in the high and low resistance states could be described by the space-charge-limited current model. The trap ionization energy Wt = 0.35 eV and the trap concentration Nd = 1.7 × 1019 cm−3 were obtained from simulation using the SCLC model in the high resistance state. For the low resistance state, the trap ionization energy was Wt = 0.01 eV, and the trap concentration was Nt = 4.6 × 1017 cm−3.

Author Contributions

Conceptualization, V.A.G. and V.A.V.; methodology, V.A.G., A.A.G. and G.N.K.; validation, V.A.G., V.A.V. and A.A.G.; formal analysis, A.A.G.; investigation, A.A.G.; resources, G.N.K.; writing—original draft preparation, A.A.G. and V.A.V.; writing—review and editing, A.A.G. and V.A.V.; visualization, A.A.G.; supervision, V.A.G. and V.A.V.; project administration, V.A.G. and V.A.V.; funding acquisition, V.A.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian Science Foundation, grant number 22-19-00369.

Acknowledgments

The work was supported by the Russian Science Foundation, project No. 22-19-00369. The electrophysical measurements were conducted using equipment from the Center of Collective Usage “VTAN” NSU supported by the Ministry of Education and Science of Russia by agreement #075-12-2021-697.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A

Table A1. Constant and fitted parameters of the Schottky effect, TAT model, Frenkel effect, and H-A model. The bold and underline parameters are parameters due to which the model does not fit experimental data.
Table A1. Constant and fitted parameters of the Schottky effect, TAT model, Frenkel effect, and H-A model. The bold and underline parameters are parameters due to which the model does not fit experimental data.
StateSchottky EffectTATFrenkel EffectH-A Model
LRSConstant:
S = 0.005 cm2 *
d = 33 nm
Fitted:
W0 = 0.07 eV
ε = 18
m* = 2.3 × 10−5 me
Constant:
d = 33 nm
Fitted:
W0 = 0.04 eV
S = 4 × 10−16 m2
m* = 6.9 me
Constant:
S = 0.005 cm2 *
d = 33 nm
Fitted:
W = 0.13 eV
ν = 3.1 × 1013 s−1
ε = 65
N = 8.0 × 106 cm−3
Constant:
S = 0.005 cm2 *
d = 33 nm
Fitted:
N = 3.5 × 1020 cm−3
W = 0.22 eV
ε = 60
ν = 8.0 × 103 s−1
HRSConstant:
S = 0.005 cm2
d = 33 nm
Fitted:
W0 = 0.15 eV
ε = 9
m* = 1.0 × 10−12 me
Constant:
S = 0.005 cm2
d = 33 nm
Fitted:
W0 = 0.15 eV
S = 1.0 × 10−14 m2
m* = 1.4 me
Constant:
S = 0.005 cm2
d = 33 nm
Fitted:
W = 0.24 eV
ν = 5.8 × 1013 sec−1
ε = 28
N = 1.0 × 102 cm−3
Constant:
S = 0.005 cm2
d = 33 nm
Fitted:
N = 3.5 × 1020 cm−3
W = 0.35 eV
ε = 20
ν = 1.2 × 101 s−1
* Area is a fitted parameter in the LRS, but the modeling had been stopped due to other parameters before the area parameter could be determined.
Table A2. Constant and fitted parameters of the ME-L model, N-G model, S-E model, and SCLC model. The bold and underline parameters are parameters due to which the model does not fit experimental data.
Table A2. Constant and fitted parameters of the ME-L model, N-G model, S-E model, and SCLC model. The bold and underline parameters are parameters due to which the model does not fit experimental data.
StateME-L ModelN-G ModelS-E ModelSCLC
LRSConstant:
S = 0.005 cm2 *
d = 33 nm
Fitted:
N = 1.2 × 109 cm−3
Wt = 0.1 eV
Wopt = 0.2 eV
Wph = 60 meV
m* = 6 me
Constant:
S = 0.005 cm2 *
d = 33 nm
Fitted:
N = 3.5 × 1020 cm−3
Wt = 0.18 eV
Wopt = 0.36 eV
m* = 53 me
Constant:
S = 0.005 cm2 *
d = 33 nm
Fitted:
I0 = 2 × 103 A
W = 0.11 eV
V0 = 0.5 eV
a = 0.6 nm
Constant:
ε = 6
d = 33 nm
m* = 0.5 me
μ = 1 cm2/(V×s)
Fitted:
Ea = 0.11 eV
Nd = 5.5 × 1019 cm−3
Wt = 0.01 eV
Nt = 4.6 × 1017 cm−3
r = 410 nm
HRSConstant:
S = 0.005 cm2
d = 33 nm
Fitted:
N = 1.0 × 103 cm−3
Wt = 0.15 eV
Wopt = 0.3 eV
Wph = 60 meV
m* = 1.1 me
Constant:
S = 0.005 cm2
d = 33 nm
Fitted:
N = 3.5 × 1020 cm−3
Wt = 0.35 eV
Wopt = 0.7 eV
m* = 47 me
Constant:
S = 0.005 cm2
d = 33 nm
Fitted:
I0 = 0.15 A
W = 0.24 eV
V0 = 0.5 eV
a = 1.4 nm
Constant:
S = 0.005 cm2
d = 33 nm
ε = 6
m* = 0.5 me
Fitted:
Ea = 0.19 eV
Nd = 4.6 × 1017 cm−3
Wt = 0.35 eV
Nt = 1.7 × 1019 cm−3
l = 1.3
μ = 9.6 × 10−11 cm2/(V×s)
* Area is a fitted parameter in the LRS, but the modeling had been stopped due to other parameters before the area parameter could be determined.

References

  1. Shi, Y.; He, L.; Guang, F.; Li, L.; Xin, Z.; Liu, R. A Review: Preparation, Performance, and Applications of Silicon Oxynitride Film. Micromachines 2019, 10, 552. [Google Scholar] [CrossRef] [Green Version]
  2. Strukov, D.B.; Snider, G.S.; Stewart, D.R.; Williams, R.S. The missing memristor found. Nature 2008, 453, 80–83. [Google Scholar] [CrossRef] [PubMed]
  3. Lee, M.J.; Lee, C.; Lee, D.; Chang, M.; Hur, J.H.; Kim, Y.-B.; Kim, C.-J.; Seo, D.H.; Seo, S.; Chung, U.-I.; et al. A fast, high-endurance and scalable non-volatile memory device made from asymmetric Ta2O5−x/TaO2−x bilayer structures. Nature Mater. 2011, 10, 625–630. [Google Scholar] [CrossRef]
  4. Prakash, A.; Jana, D.; Maikap, S. TaOx-based resistive switching memories: Prospective and challenges. Nanoscale Res. Lett. 2013, 8, 418. [Google Scholar] [CrossRef] [Green Version]
  5. Meziani, S.; Moussi, A.; Mahiou, L.; Outemzabet, R. Compositional analysis of silicon oxide/silicon nitride thin films. Mater. Sci. 2016, 34, 315–321. [Google Scholar] [CrossRef] [Green Version]
  6. Iwase, Y.; Horie, Y.; Daiko, Y.; Honda, S.; Iwamoto, Y. Synthesis of a Novel Polyethoxysilsesquiazane and Thermal Conversion into Ternary Silicon Oxynitride Ceramics with Enhanced Thermal Stability. Materials 2017, 10, 1391. [Google Scholar] [CrossRef] [Green Version]
  7. Mehonic, A.; Shluger, A.L.; Gao, D.; Valov, I.; Miranda, E.; Ielmini, D.; Bricalli, A.; Ambrosi, E.; Li, C.; Yang, J.J.; et al. Silicon oxide (SiOx): A promising material for resistance switching? Adv. Mater. 2018, 30, 1801187. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Yao, J.; Zhong, L.; Natelson, D.; Tour, J.M. Silicon oxide: A non-innocent surface for molecular electronics and nanoelectronics studies. J. Am. Chem. Soc. 2011, 133, 941. [Google Scholar] [CrossRef] [Green Version]
  9. Gismatulin, A.A.; Voronkovskii, V.A.; Kamaev, G.N.; Novikov, Y.N.; Kruchinin, V.N.; Krivyakin, G.K.; Gritsenko, V.A.; Prosvirin, I.P.; Chin, A. Electronic structure and charge transport mechanism in a forming-free SiOx-based memristor. Nanotechnology 2020, 31, 505704. [Google Scholar] [CrossRef]
  10. Gismatulin, A.A.; Kamaev, G.N.; Kruchinin, V.N.; Gritsenko, V.A.; Orlov, O.M.; Chin, A. Charge transport mechanism in the forming-free memristor based on silicon nitride. Sci. Rep. 2021, 11, 2417. [Google Scholar]
  11. Winn, E.B. The Temperature Dependence of the Self-Diffusion Coefficients of Argon, Neon, Nitrogen, Oxygen, Carbon Dioxide, and Methane. Phys. Rev. 1950, 80, 1024–1027. [Google Scholar] [CrossRef]
  12. Kijima, K.; Shirasaki, S. Nitrogen self-diffusion in silicon nitride. J. Chem. Phys. 1976, 65, 2668–2671. [Google Scholar] [CrossRef]
  13. Roussel, M.; Talbot, E.; Pareige, P.; Gourbilleau, F. Influence of the supersaturation on Si diffusion and growth of Si nanoparticles in silicon-rich silica. J. Appl. Phys. 2013, 113, 063519. [Google Scholar] [CrossRef] [Green Version]
  14. Volodin, V.A.; Geydt, P.; Kamaev, G.N.; Gismatulin, A.A.; Krivyakin, G.K.; Prosvirin, I.P.; Azarov, I.A.; Fan, Z.; Vergnat, M. Resistive Switching in Non-Stoichiometric Germanosilicate Glass Films Containing Ge Nanoclusters. Electronics 2020, 9, 2103. [Google Scholar] [CrossRef]
  15. Chen, D.; Huang, S.; He, L. Effect of oxygen concentration on resistive switching behavior in silicon oxynitride film. J. Semicond. 2017, 38, 043002. [Google Scholar] [CrossRef]
  16. Vasileiadis, N.; Karakolis, P.; Mandylas, P.; Ioannou-Sougleridis, V.; Normand, P.; Perego, M.; Komninou, P.; Ntinas, V.; Fyrigos, I.-A.; Karafyllidis, I.; et al. Understanding the role of defects in silicon nitride-based resistive switching memories through oxygen doping. IEEE Trans. Nanotechnol. 2021, 20, 365. [Google Scholar] [CrossRef]
  17. Leng, K.; Zhu, X.; Ma, Z.; Yu, X.; Xu, J.; Xu, L.; Li, W.; Chen, K. Artificial neurons and synapses based on Al/a-SiNxOy:H/P+-Si device with tunable resistive switching from threshold to memory. Nanomaterials 2022, 12, 311. [Google Scholar] [CrossRef]
  18. Yu, Z.; Aceves, M.; Carrillo, J.; Lopez-Estopier, R. Charge trapping and carrier transport mechanism in silicon-rich silicon oxynitride. Thin Solid Films 2006, 515, 2366–2372. [Google Scholar] [CrossRef]
  19. Yonamoto, Y.; Inaba, Y.; Akamatsu, N. Compositional dependence of trap density and origin in thin silicon oxynitride film investigated using spin dependent Poole–Frenkel current. Appl. Phys. Lett. 2011, 98, 232905. [Google Scholar] [CrossRef]
  20. Novikov, Y.N.; Gismatulin, A.A.; Prosvirin, I.P.; Bobovnikov, P.G.; Krasnikov, G.Y.; Gritsenko, V.A. Short-range order and charge transport in silicon-rich pyrolytic silicon oxynitride. J. Non Cryst. Solids 2023, 599, 121984. [Google Scholar] [CrossRef]
  21. Perevalov, T.V.; Volodin, V.A.; Kamaev, G.N.; Gismatulin, A.A.; Cherkova, S.G.; Prosvirin, I.P.; Astankova, K.N.; Gritsenko, V.A. Electronic structure of silicon oxynitride films grown by plasma-enhanced chemical vapor deposition for memristor application. J. Non Cryst. Solids 2022, 598, 121925. [Google Scholar] [CrossRef]
  22. Antonenko, A.K.; Volodin, V.A.; Efremov, M.D.; Zazulya, P.S.; Kamaev, G.N.; Marin, D.V. Silicon nitride films deposited at substrate temperatures <100 ◦C in a permanent magnet electron cyclotron resonance plasma. Optoelectron. Instrum. Data Process. 2011, 47, 459. [Google Scholar]
  23. Jensen, K.L. General formulation of thermal, field and photoinduced electron emission. J. Appl. Phys. 2007, 102, 024911. [Google Scholar] [CrossRef]
  24. Roberts, G.G.; Polango, J.I. Thermally assisted tunneling in dielectric films. Phys. Stat. Sol. a 1970, 1, 409–420. [Google Scholar] [CrossRef]
  25. Gritsenko, V.A.; Meerson, E.E.; Morokov, Y.N. Thermally assisted hole tunneling at the Au-Si3N4 interface and the energy-band diagram of metal-nitride-oxide-semiconductor structures. Phys. Rev. B 1998, 57, R2081–R2083. [Google Scholar] [CrossRef] [Green Version]
  26. Frenkel, J. On the theory of electric breakdown of dielectrics and electronic semiconductors. Tech. Phys. USSR 1938, 5, 685–695. [Google Scholar]
  27. Frenkel, J. On pre-breakdown phenomena in insulators and electronic semiconductors. Phys. Rev. B 1938, 54, 647. [Google Scholar] [CrossRef]
  28. Hill, R.M. Poole-Frenkel conduction in amorphous solids. Philos. Mag. 1971, 23, 59–86. [Google Scholar] [CrossRef]
  29. Adachi, H.; Shibata, Y.; Ono, S. On electronic conduction through evaporated silicon oxide films. J. Phys. D Appl. Phys. 1971, 4, 988–994. [Google Scholar] [CrossRef]
  30. Makram-Ebeid, S.S.; Lannoo, M. Quantum model for phonon-assisted tunnel ionization of deep levels in a semiconductor. Phys. Rev. B 1982, 25, 6406. [Google Scholar] [CrossRef]
  31. Nasyrov, K.A.; Gritsenko, V.A. Charge transport in dielectrics via tunneling between traps. J. Appl. Phys. 2011, 109, 093705. [Google Scholar] [CrossRef]
  32. Shklovskii, B.I. Percolation Mechanism of Electrical-Conduction in Strong Electric-Fields. Sov. Phys. Semicond. 1979, 13, 53–56. [Google Scholar]
  33. Shklovskii, B.I. & Efros, A.L. Percolation Theory and Conductivity of Highly Inhomogeneous-Media. Usp. Fiz. Nauk. 1975, 117, 401–435. [Google Scholar]
  34. Mott, N.F.; Gurney, R.W. Electronic Processes in Ionic Crystals. J. Phys. Chem. 1941, 45, 1142. [Google Scholar]
  35. Lampert, M.A. Simplified Theory of Space-Charge-Limited Currents in an Insulator with Traps. Phys. Rev. 1956, 103, 1648. [Google Scholar] [CrossRef]
  36. Murgatroyd, P.N. Theory of space-charge-limited current enhanced by Frenkel effect. J. Phys. D: Appl. Phys. 1970, 3, 151–156. [Google Scholar] [CrossRef]
  37. Lampert, M.A.; Mark, P. Current Injection in Solids; Academic Press: New York, USA, 1970; p. 351. [Google Scholar]
  38. Bloomfield, P.; Steiger, W. Least Absolute Deviations, Theory, Applications and Algorithms; Huber, P., Rosenblatt, M., Eds.; Birkhäuser: Boston, MA, USA, 1983; p. 349. [Google Scholar]
  39. Fowler, R.H.; Nordheim, L. Electron emission in intense electric fields. Proc. R. Soc. Lond. A 1928, 119, 173–181. [Google Scholar]
Figure 1. Schematic diagram of contact-limited models: (a) Schottky effect, (b) thermally assisted tunneling model.
Figure 1. Schematic diagram of contact-limited models: (a) Schottky effect, (b) thermally assisted tunneling model.
Electronics 12 00598 g001
Figure 2. Schematic diagram of bulk-limited models: (a) Frenkel effect, (b) Hill–Adachi model of overlapping Coulomb potentials, (c) Makram–Ebeid and Lannoo multiphonon isolated trap ionization model, (d) Nasyrov–Gritsenko phonon-assisted tunneling via neighboring traps model.
Figure 2. Schematic diagram of bulk-limited models: (a) Frenkel effect, (b) Hill–Adachi model of overlapping Coulomb potentials, (c) Makram–Ebeid and Lannoo multiphonon isolated trap ionization model, (d) Nasyrov–Gritsenko phonon-assisted tunneling via neighboring traps model.
Electronics 12 00598 g002
Figure 3. IV characteristics of the PECVD SiO0.9N0.6-based memristor.
Figure 3. IV characteristics of the PECVD SiO0.9N0.6-based memristor.
Electronics 12 00598 g003
Figure 4. The IV characteristics of the PECVD SiO0.9N0.6-based memristor and simulation results by the contact-limited models: (a) Schottky effect in LRS; (b) the TAT model in LRS; (c) the Schottky effect in HRS; and (d) the TAT model in HRS.
Figure 4. The IV characteristics of the PECVD SiO0.9N0.6-based memristor and simulation results by the contact-limited models: (a) Schottky effect in LRS; (b) the TAT model in LRS; (c) the Schottky effect in HRS; and (d) the TAT model in HRS.
Electronics 12 00598 g004
Figure 5. IV characteristics of the PECVD SiO0.9N0.6-based memristor and the results of simulation by the bulk-limited models in the HRS: (a) Frenkel effect; (b) Hill–Adachi model of overlapping Coulomb potentials; (c) Makram–Ebeid and Lannoo model of multiphonon isolated trap ionization; (d) Nasyrov–Gritsenko model of phonon-assisted tunneling between traps; (e) Shklovskii–Efros percolation model; and (f) SCLC model.
Figure 5. IV characteristics of the PECVD SiO0.9N0.6-based memristor and the results of simulation by the bulk-limited models in the HRS: (a) Frenkel effect; (b) Hill–Adachi model of overlapping Coulomb potentials; (c) Makram–Ebeid and Lannoo model of multiphonon isolated trap ionization; (d) Nasyrov–Gritsenko model of phonon-assisted tunneling between traps; (e) Shklovskii–Efros percolation model; and (f) SCLC model.
Electronics 12 00598 g005
Figure 6. IV characteristics of the PECVD SiO0.9N0.6-based memristor, and the results of simulation by the bulk-limited models in HRS: (a) Frenkel effect; (b) Hill–Adachi model of overlapping Coulomb potentials; (c) Makram–Ebeid and Lannoo multiphonon isolated trap ionization model; (d) Nasyrov–Gritsenko model of phonon-assisted tunneling between traps, (e) Shklovskii–Efros percolation model; and (f) the SCLC model.
Figure 6. IV characteristics of the PECVD SiO0.9N0.6-based memristor, and the results of simulation by the bulk-limited models in HRS: (a) Frenkel effect; (b) Hill–Adachi model of overlapping Coulomb potentials; (c) Makram–Ebeid and Lannoo multiphonon isolated trap ionization model; (d) Nasyrov–Gritsenko model of phonon-assisted tunneling between traps, (e) Shklovskii–Efros percolation model; and (f) the SCLC model.
Electronics 12 00598 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gismatulin, A.A.; Kamaev, G.N.; Volodin, V.A.; Gritsenko, V.A. Charge Transport Mechanism in the Forming-Free Memristor Based on PECVD Silicon Oxynitride. Electronics 2023, 12, 598. https://doi.org/10.3390/electronics12030598

AMA Style

Gismatulin AA, Kamaev GN, Volodin VA, Gritsenko VA. Charge Transport Mechanism in the Forming-Free Memristor Based on PECVD Silicon Oxynitride. Electronics. 2023; 12(3):598. https://doi.org/10.3390/electronics12030598

Chicago/Turabian Style

Gismatulin, Andrei A., Gennadiy N. Kamaev, Vladimir A. Volodin, and Vladimir A. Gritsenko. 2023. "Charge Transport Mechanism in the Forming-Free Memristor Based on PECVD Silicon Oxynitride" Electronics 12, no. 3: 598. https://doi.org/10.3390/electronics12030598

APA Style

Gismatulin, A. A., Kamaev, G. N., Volodin, V. A., & Gritsenko, V. A. (2023). Charge Transport Mechanism in the Forming-Free Memristor Based on PECVD Silicon Oxynitride. Electronics, 12(3), 598. https://doi.org/10.3390/electronics12030598

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop