The effect of doping Sb on the electronic structure and the device characteristics of Ovonic Threshold Switches based on Ge-Se.

The Ovonic Threshold Switch (OTS) based on an amorphous chalcogenide material has attracted much interest as a promising candidate for a high-performance thin-film switching device enabling 3D-stacking of memory devices. In this work, we studied on the electronic structure of amorphous Sb-doped Ge(0.6)Se(0.4) (in atomic mole fraction) film and its characteristics as to OTS devices. From the optical absorption spectroscopy measurement, the band gap (Eg) was found to decrease with increasing Sb content. In addition, as Sb content increased, the activation energy (Ea) for electrical conduction was found to decrease down to about one third of Eg from a half. As to the device characteristics, we found that the threshold switching voltage (Vth) drastically decreased with the Sb content. These results, being accountable in terms of the changes in the bonding configuration of constituent atoms as well as in the electronic structure such as the energy gap and trap states, advance an effective method of compositional adjustment to modulate Vth of an OTS device for various applications.

Electronic switches including transistors and diodes are the most elementary building blocks of numerous electronic devices such as logic, memory, and display devices. As a result of incessant scale-down, switches based on single crystalline silicon are now considered to confront their limits in both performances and scalability. Scientists and engineers are struggling to find solutions which often exploit new materials, new device structures, and new operation principles. Among them, increasing attention has been paid to the Ovonic Threshold Switching (OTS) device owing to its high-current-driving capability and simple device structure1. S. R. Ovshinsky firstly reported reversible electrical switching phenomenon using a chalcogenide composed of Te, As, Si, and Ge in 1968, featured by a sudden huge drop in its resistance at a certain bias and the recovery of its initial high resistance state upon removal of the bias23. About 40 years after the discovery, a successful integration into a phase-change memory was demonstrated and a switching speed being as fast as several nsec was also shown promising45. However, there still exist a few issues to deal with for achieving a reliable OTS device. For one issue, most of the chalcogenide materials used in the OTS device consist of environment-unfriendly compounds containing As. To resolve this issue, we investigated the feasiblity of more environment-friendly and simple binary chalcogenide materials, GexSe1-x (x = 0.4~0.7), and demonstrated their reliable characteristics as materials for OTS devices6. As another issue, one should be able to tailor the operational characteristics of such OTS in effective and controllable manners for these devices to be integrated in various functional devices. For that matter, threshold switching voltage (Vth), a fundamental parameter characterizing an OTS, is to be singled out. According to a recent model developed by D. Ielmini, Vth depends on the energy gap (Eg), the density (ntrap) and the energy level (Etrap) of trap states78. Doping of group-V elements as well as varying the kind of chalcogen atoms (S, Se, Te) could modify the electronic structure of the amorphous chalcogenide with concomitant modulation of the switching characteristics, as shown with N- and Bi-doped Ge0.6Se0.4 systems910. Herein, we have extended the approach to Sb and systematically studied the effect of its doping in Ge0.6Se0.4 on the electronic structures and the characteristics of OTS devices. In the previous reports111213, Sb was shown to result in the reduction of Eg when doped in various amorphous GexSey materials.

Results and Discussion

Figure 1(a) shows the results of the optical absorption spectroscopy measurements for various Sbx(Ge0.6Se0.4)1-x films. In many amorphous semiconductors, the dependence of the absorption coefficient (α) on the photon energy (hν) is well expressed by Tauc's relation, for hν > Eg, with a broad tail structure due to trap states for hν < Eg14. In Figure 1(a), it is found that depends linearly on hν for each film in the high photon energy regions. From the slope of each curve, Eg of a respective film was obtained and the resultant Eg values are plotted as a function of Sb content in Figure 1(c). It is apparent that Eg decreases with increasing Sb content.

Figure 1

(a) Tauc plot ((αhν)1/2 vs. hν) of Sb(Ge0.6Se0.4)1- films. The dashed lines are the linear fits to the data. (b) Semilog-plot of the resistance (R) of the Sb(Ge0.6Se0.4)1- films as a function of the reciprocal of the temperature (T) (Arrhenius plot). The solid lines are the linear fits to the data. The curves are vertically shifted for clarity (multiplications factors are indicated on the left of the curve). (c) The energy gap (Eg, solid red circle) and the activation energy for electrical conduction (Ea: solid black square) as a function of Sb contents. Eg/3 is also plotted for comparison with Ea.

In Figure 1(b), we plotted the electrical resistance (R) as a function of the reciprocal temperature (1/T) for each film. T was increased from room temperature to 420°C at a rate of 5°C/min. For all the samples, variation in resistance with temperature was found to be well described by Arrhenius' law (), where Ea and k are the activation energy for the electrical conduction and Boltzmann constant, respectively. It should be mentioned that the observed data were also attempted to fit to the variable-range hopping model with an R ~ exp[(T0/T)1/4] dependence but in vain (see Supplementary Information, Figure S1). In Figure 1(b), linear fitting curves are also shown (solid lines) while their slopes representing Ea's of the films. Interestingly, Ea was found to approach Eg/3 with increasing Sb content from Eg/2 in the Sb-free case as shown in Figure 1(c). In most of the amorphous chalcogenide materials examined, Ea was usually found to be about one half of Eg and this was understood as resulting from the Fermi level pinning at the midgap due to a huge amount of charged defect centers13. Such a relationship (Ea ~ Eg/3) was originally observed in BiGe0.2Se0.8-15 and interpreted in terms of the incorporation of positively charged impurities16. In addition, such a relationship was found in recent studies on Ge0.4Te0.6-Sb17 and Bi(Ge0.5Se0.5)1-10 Considering the electronegativity of Sb (χSb = 2.05 eV) not much different from that of Bi (χBi = 2.02 eV), but much lower than that of Se (χSb = 2.55 eV), the amount of positively charged impurities is expected to be similarly large in Sb(Ge0.6Se0.4)1- with a high x, leading to the observation.

Meanwhile, the crystallization temperature (Tx) values can be also determined from the measured R vs. T curves as R decreases drastically upon crystallization (see Figure S2 in Supplementary Information). Tx was found to decrease with increasing Sb content but was higher than 350°C even at the highest Sb content tested i.e. x = 0.22. This guarantees that the Sb-doped GeSe is not liable to crystallize during device fabrication around at 350°C, a prerequisite for a material to be adopted in the OTS.

Figure 2(a) and 2(b) show schematic illustrations of an OTS device using Sb(Ge0.6Se0.4)1- and of the electrical measurement setup for device characterization, respectively. To find Vth of each device, triangular voltage pulses were applied instead of DC voltage which forces enormous thermal and electrical stresses resulting in irreversible changes in device characteristics. A pulse composed of 1 μsec (rising time), 100 nsec (width), and 1 μsec (falling time) was applied while both the voltage across and the current through the device were measured with two oscilloscopes, OSC1 and OSC2, respectively. The size of a sample was varied from 2*2 μm2 to 20*20 μm2 and the data presented in this paper were obtained with devices of 5*5 μm2 in size. Figure 2(c) shows typical current (I) – voltage (V) curves of OTS devices of each composition. Vth is defined as the voltage at which a drastic change in device resistance comes about, initiating a snap-back in the I–V curve. In Figure 2(d), Vth is shown to decrease remarkably as Sb content increases. The observed reduction in Vth agrees qualitatively with the reduction of Eg, which is expected from the aforementioned Ielmini's model78. But quantitatively, such a remarkable reduction of Vth by about 33.3% (from 3.24 V to 2.16 V as Sb content increases from 0 to 22 at. %) may not be solely ascribed to the reduction of Eg by about 17% (Eg = 0.94 and 0.78 eV for Ge0.6Se0.4 and Sb0.22(Ge0.6Se0.4)0.78, respectively). According to our calculation using the Ielmini's model (See Figure S3 in Supplementary Information), Vth was shown to linearly depend on Eg with a slope of about one. Apart from Eg, two material parameters have significant effects on the characteristics of an OTS device. One is the density of trap states (nT) and the other is the energy difference (ΔE) between the trap level and the edge of the neighboring extended band (either valence band (VB) or conduction band (CB)). In our previous work18, we derived an expression describing the dependence of the delay time (tdel) of an OTS on the bias voltage based on a simplified model assuming thermally-activated hopping transport and trap states of a single kind, from which nT and ΔE and can be estimated. In Figure 3(a), a set of panels is shown with each representing a contour map of the current waveforms at varying Va vlaues, flowing through a typical OTS device of the respective film composition. The color scale represents the current levels measured with OSC2 shown in Figure 2(b). For these measurements, we applied a 200 ns-long voltage pulse of varying amplitude while registering the applied voltage (Va) with OSC1 and the current flowing through the device with OSC2. A representative waveform at Va = 3.68 V is shown in the inset for an undoped sample (leftmost panel). The black dashed lines indicate the beginning of the applied voltage pulse. It is apparent that Vth decreases as Sb content increases, consistent with Figure 2(c). Furthermore, it is observed that a threshold switching event occurs at tdel after the voltage pulse is applied and that tdel decreases with increasing Va. In ref. 18, the dependence of tdel on the bias voltage was shown to have the following expression.

Figure 2

(a) and (b) Schematic illustrations of a cross-point OTS device and the measurement setup. (P/G: pulse generator, SMU: source-measure unit, Rload: load resistor). (c) Characteristic current (I)-voltage (V) curves for typical OTS devices using Sb(Ge0.6Se0.4)1-. It was measured by applying a 2 μs-long triangular pulse and recording the voltage- and current-pulse waveforms by osc1 and osc2, respectively. The curves are vertically shifted in the order of Sb contents. (d) The dependence of Vth on the Sb contents, where Vth is defined in (c).

Figure 3

(a) Contour maps of current waveforms recorded in osc. 2 (in figure 2(b)) with varying x. Color scale represents the amplitude of current. 200 ns-long rectangular pulses were applied with increasing amplitude and the beginning of the voltage pulse is indicated by the dashed lines. For x = 0 case, a current waveform at Va = 3.68 V is displayed in the inset. (b) The symbols represent tdel as a function of Va for each composition. The solid lines are fitting curves by Eq. (1) in the text. Error bars means the standard deviation. (c) and (d) ntrap and ΔE( = EC − ET) normalized by that of the x = 0 sample as a function of Sb contents, respectively.

In Eq. (1), τ0, E, ET, Va, and k are the escape-trial time, conduction band minimum, energy level of trap states, applied voltage, and the Boltzmann constant, respectively. And q* is the effective charge defined by eΔz/ua, where Δz and ua are the inter-trap distance and the thickness of the amorphous chalcogenides layer, respectively. In Figure 3(b), the curves of measured tdel(Va) values are shown for OTS devices with Sb(Ge0.6Se0.4)1- together with the fitting curves due to Eq. (1), where a data point with an error bar represents an average and a standard deviation taken over about 20 measurements using more than five devices. It is clear that, at low Va region, the rate of change in tdel with Va gets larger as Sb content increases. Conceivably, this is associated with the increase of the trap density (ntrap) with Sb content. With increasing trap density, traps come closer to each other and the barrier height between them gets lower due to the Poole-Frenkel mechanism19. Such lowering of the barrier height leads to an exponential decrease in the escape-time from a trap hence in tdel. The calculated ntrap and ΔE values are plotted in Figure 3(c) and 3(d), respectively. In Figure 3(c), the ntrap is found to increase as Sb content increases, consistent with the above explanation. In Figure 3(d), ΔE is shown to decrease as Sb content increases, indicating that the trap states come closer to VB or CB.

The observed variations in ntrap and ΔE are thought to be associated with the change in the electronic structrure accompanying the one in bond configuration by Sb-doping in Ge0.6Se0.4. In amorphous materials, the atomic structure is determined by the average coordination number () given by the average of the coordination numbers of constituent atomic species weighted by the respective composition of each species. In our Sb(Ge0.6Se0.4)1- system, ranges from 2.8 (for x = 0) to 3.156 (for x = 0.22), which are higher than Tanaka's threshold ( = 2.67)20, indicating that the investigated materials in this work have 3D-network structures. The bond configurations constituting 3D-network structures may be conjectured according to the chemical ordered network model (CONM) proposed by Bicerano and Ovshinsky21. For this purpose, the bonding energies between constituent atoms should be known and they were calculated by using Pauling's relation in a previous paper as listed in Table 122. From the listed bonding energies, the followings may be inferred as to bonding configurations. Firstly, the higher bonding energy of a Ge–Se pair than those of the other pairs i.e. Ge–Ge and Se–Se would lead to the preferential formation of Ge(Se1/2)4 tetrahedron in Ge–Se binary materials. Secondly, in the particular case of Ge0.6Se0.4 material characterized by a relative Se-deficient composition, all Se atoms are bonded with Ge atoms to form Ge(Se1/2)4 tetrahedra and the remaining Ge atoms form homopolar Ge–Ge bonds. When Sb atoms are incorporated in Ge0.6Se0.4, Ge(Se1/2)4 remain intact because of the highest bonding energy whereas a portion of Ge homopolar bonds would be replaced with Ge–Sb bonds instead because of the preference of the heteropolar bonds to the homopolar bonds.

Table 1

Energy of possible bonds in GexSeySbz compounds

Bonds	Energy (kcal/mol)
Ge–Se	49.42
Sb–Se	43.96
Se–Se	44.00
Ge–Ge	37.60
Ge–Sb	33.76
Sb–Sb	30.22

This change in bond configuration can result in significant changes in the electronic structure. The suggested change is schematized in Figure 4. In amorphous chalcogenides, it is assumed that a chalcogen atom has a two-fold coordination with its neighbors and the lone-pair (LP) electrons are configured spontaneously to form valence-alternation pairs (VAPs)323. VAPs are composed of two oppositely charged chalcogen atoms, leading respectively to donor-like trap (C3+) and acceptor-like traps (C1−) inside the band gap. We assume that this electronic structure may apply essentially to undoped amorphous Ge0.6Se0.4. For Sb-doped Ge0.6Se0.4, it should be noted that the preferred bond between Ge and Sb has the bonding energy lower than Ge homopolar bonds by about 3.9 kcal/mol. Weak bonds are expected to create localized anti-bonding states in the band gap leading to the Mott delocalization23. In addition, since an Sb atom has one more electron than a Ge atom, an Sb atom is expected to be ionized on forming a bond with Ge. These positively-charged impurities might create extra defect states broadening the trap states near CB. In this picture, it should be addressed that the assumed coordination number of an Sb atom is larger than the expeced ( = 3) from the “8-N” rule, where N represents the number of valence electrons of a constituent atomic species (5 for Sb). The existence of over-coordinated Sb atoms is supported by a previous EXAFS study24 and a density functional calculations25, where the coordination number of Sb in Ge2Sb2Te5 was found to be higher than three. This model can naturally explain the above observations and the resultant drastic reduction of Vth by Sb-doping.

Figure 4

Schematic illustration of the change in bond configuration and the electronic structure by doping Sb in Ge0.6Se0.4.

In Figure S4, we plotted the off-state resistance (Roff), the on-state resistance (Ron), and the ratio between them (Roff/Ron) as a function of Sb contents. Roff was measured by biasing DC 0.2 V and (Ron) was measured by dividing the voltage amplitude of an applied pulse (200 ns-long) by the current amplitude. For x ≤ 0.10, Roff/Ron is found to be higher than 105.

Considering that Vth changed mainly in the range of 0 < x < 0.1 (see Figure 2(d)), this work is believed to provide not only a method to modulate the threshold voltage of an OTS device for various applications by doping group-V elements, but also a promising candidate chalcogenide material to be employed for the development of a high-performance switching device beyond the Si era.

Methods

Sb(Ge0.6Se0.4)1- thin films were fabricated by co-sputtering technique using Sb, Ge, and Ge0.4Se0.6 targets, while the film composition was varied with changing RF (radio-frequency) power imposed on each target. The Sb contents were determined to be 0, 8.3, 10.4, 17.4 and 22.5 at. % respectively by Electron Probe Micro Analysis (EPMA) and X-Ray Fluorescence (XRF) analysis. UV/VIS spectroscopy was performed with a 100 nm-thick film deposited on a corning glass substrate (t = 0.5 mm) in the wavelength range of 300~1600 nm using UV/VIS/NIR spectrometer (UV 3600, Shimatzu). The electrical resistance (R) was also measured as a function of temperature (T) under the Ar atmosphere of 1 Torr using a standard four-point probe method. From the resulting R vs.T plot the activation energy (Ea) for electrical conduction as well as the crystallization temperature (Tx) was determined for each composition. The former is known to be related to Eg in amorphous chalcogenide due to the pinning of the Fermi level. The latter sets a limit in thermal loading below which the local structure of the used amorphous chalcogenide material can be secured. OTS devices were fabricated with each device having a cross-point structure composed of metal(bottom electrode, Pt 50 nm)/Sb(Ge0.6Se0.4)1-(100 nm)/metal(top electrode, TiN 50 nm) using photolithography and lift-off techniques. The size of a cross-point device was defined with the overlapping area of the top and bottom electrode metals. The device characteristics were investigated using a measurement setup composed of a source-measurement unit (Keithley 236), a multi-channel oscilloscope (Tektronix TDS 5104), and a pulse generator (agilent 81110A) (Figure 2(b)).

Author Contributions

S.L. planned and managed the project. S.-Y.S. performed the experiments with the help of J.M.C., J.S. and H.-W.A. B.-k.C., Y.G.C. and S.L. took care of the data analyses and wrote the manuscript.