Artificial Synaptic Emulators Based on MoS2 Flash Memory Devices

Aug 14, 2018 - We fabricated MoS2-based flash memory devices by stacking MoS2 and hexagonal boron nitride (hBN) layers on an hBN/Au substrate and ...
0 downloads 0 Views 1MB Size
Subscriber access provided by University of South Dakota

Functional Inorganic Materials and Devices

Artificial synaptic emulators based on MoS2 flash memory devices with double floating gates Sum-Gyun Yi, Myung Uk Park, Sung Hyun Kim, Chang Jun Lee, Junyoung Kwon, Gwan-Hyoung Lee, and Kyung-Hwa Yoo ACS Appl. Mater. Interfaces, Just Accepted Manuscript • DOI: 10.1021/acsami.8b10203 • Publication Date (Web): 14 Aug 2018 Downloaded from http://pubs.acs.org on August 15, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 27 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Artificial synaptic emulators based on MoS2 flash memory devices with double floating gates Sum-Gyun Yi1, Myung Uk Park1, Sung Hyun Kim1, Chang Jun Lee1, Junyoung Kwon2, GwanHyoung Lee2 and Kyung-Hwa Yoo1*. 1

Department of Physics, Yonsei University, 50 Yonsei-ro, Seoul, 03722, Republic of Korea.

2

Department of Materials Science and Engineering, Yonsei University, 50 Yonsei-ro, Seoul,

03722, Republic of Korea. KEYWORDS. synaptic emulator, neuromorphic, MoS2, double floating gate, flash memory device

Abstract

We fabricated MoS2-based flash memory devices by stacking MoS2 and hBN layers on an hBN/Au substrate, and demonstrated that these devices can emulate various biological synaptic functions, including potentiation and depression processes, spike-rate-dependent plasticity, and spike-timing dependent plasticity. In particular, compared to a flash memory device prepared on an hBN substrate, the device fabricated on the hBN/Au exhibited considerably more symmetric and linear bidirectional gradual conductance change curves, 1

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

which may be attributed to the device structure incorporating double floating gate. For the device on the hBN/Au, electron transfers may occur between the floating gate MoS2 and Au, as well as between the floating gate MoS2 and the channel MoS2, allowing for more control over electron tunneling and injection. To test our hypothesis, we also fabricated a MoS2based flash memory device on an hBN/Pd substrate and found similar behavior to the device fabricated on hBN/Au. Our results demonstrate that flexible synaptic electronics may be implemented using MoS2-based flash memory devices with double floating gates.

2

ACS Paragon Plus Environment

Page 2 of 27

Page 3 of 27 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Introduction Significant efforts have been made to develop biologically inspired neuromorphic networks with the capacity to mimic the massive parallelism and low-power operation found in the human brain, as the human brain can outperform supercomputers in many respects.1-4 To date, a large number of reports have investigated neural networks using complementary metal-oxide-semiconductor (CMOS) technology.5-7 However, implementing the functions of biological synapses in CMOS-based circuits remains challenging because synaptic weights are non-volatile and adjustable, which is uncommon in conventional CMOS devices.8,9 Previously, a variety of two-terminal memristors have been explored for the implemention of synaptic devices, and synaptic learning rules such as spike-rate-dependent plasticity (SRDP) and spike-timing dependent plasticity (STDP) have been demonstrated in various material systems including AgS2,10,11 WOx,12 Ag-doped TiO2,13 and Ge2Sb2Te5.14 More recently, three-terminal devices have been proposed as artificial synapses for emulating biological synaptic functions since three-terminal devices can realize learning and memory functions concurrently.15-19 In recent years, atomically thin two dimensional (2D) materials such as graphene, semiconducting transition metal dichalcogenides, and insulating hexagonal boron nitride (hBN) have received considerable attention owing to their possible application in ultrathin, flexible, and nearly transparent electronic and optoelectronic devices20-23. For example, MoS2 and hBN have been incorporated into nonvolatile memory structures serving as the floating gate and tunnel barrier, respectively22,23. In this study, we investigate whether a threeterminal memory device fabricated from stacked 2D materials can be applied to emulate biological synaptic functions. To do so, we fabricated a MoS2-based flash memory device on an hBN/Au substrate and compared it to one fabricated on an hBN substrate. The former 3

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

exhibited symmetric characteristics of potentiation and depression, whereas the latter exhibited asymmetric and non-linear directional gradual conductance change curves. The advantage of the synaptic emulator based on the flash memory device compared to the two terminal memristor is that the floating-gate synaptic emulator allows simultaneous memory reading and writing, so additional selection devices, such as transistors and diodes, are not required to construct a synaptic array.24 In addition, the flash memory device is operated by electron injection and tunneling rather than ion movements, so the operating speed of floating-gate synaptic device will be faster than that of the ion-gated synaptic transistors.15-19 To understand the different behavior of MoS2-based flash memory devices on hBN and hBN/Au substrates, we measured the workfunction of MoS2 on hBN and hBN/metal substrates using a scanning Kelvin probe microscope (SKPM),25-27 since all of atoms are exposed to the surface in 2D materials, so their properties, such as workfunction (i.e., the energy required to extract an electron from a material to a vacuum) may be easily modulated by the surrounding environment.28 Indeed, the measurements revealed that the workfunction of MoS2 on the hBN/metal differs from that of MoS2 on hBN, and that it depends on the metal underneath the hBN layer. Based on these findings, we propose energy band diagrams showing how Au acts as another floating gate for flash memory devices on hBN/Au, resulting in more controllable electron tunneling and injection. Because the symmetry of potentiation and depression is the most important requirement among the various synaptic behaviors,29 these results demonstrate that flexible synaptic electronics may be implemented using MoS2based flash memory devices with double floating gates.

Results and discussion

4

ACS Paragon Plus Environment

Page 4 of 27

Page 5 of 27 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

First, we investigated whether the workfunction of few-layer MoS2 is affected by metals with different workfunctions through an hBN layer. A specimen was prepared by depositing Au (42 nm) and Al (42 nm) films on a Si substrate and subsequently transferring a stacked layer consisting of hBN (3 nm) and MoS2 (5 nm) over the Au and Al films (Figures 1a and S1). The contact potential difference (VCPD) image was acquired using the SKPM (Figure 1b), since the VCPD between a SKPM tip and a sample is given by

ܸ஼௉஽ =

ఝ೟೔೛ ିఝೞೌ೘೛೗೐

(1)

ି௘

where ϕtip and ϕsample are the workfunctions of the tip and sample, respectively, and e is the electronic charge. The green line in Figure 1b represents the VCPD measured across uncovered Au and Al films with hBN (Figure 1c). The difference in VCPD values between Au and Al films was about 0.4 V, which is smaller than the reference value of 0.8-1.2 eV.30 However, when the CPD image was measured before transfer, the difference in VCPD values between the Au and Al films was estimated to be about 0.8 V (Figure S3), implying that contamination occurring during the transfer of hBN and MoS2 might decrease the VCPD difference between the Au and Al films. Figure 1d shows the line profiles of VCPD in the y - direction denoted using brown, blue, and pink lines in Figure 1b. The MoS2 flake exhibited different values of VCPD depending on the material layers. Assuming that the workfunction of MoS2 measured on hBN is close to the reported value of 4.3 – 4.6 eV,22 the workfunction of MoS2 increased by ~0.1 eV on hBN/Au, and decreased by ~0.08 eV on hBN/Al. The line profile of VCPD in the x – direction, denoted using a red line in Figure 1c, also showed that the VCPD of MoS2 on hBN/Au was higher than

5

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

that of MoS2 on hBN/Al by ~0.18 eV. These results indicate that the workfunction of MoS2 might be modulated by Au or Al via the hBN layer. To confirm the workfunction modulation of MoS2 through the hBN layer, we fabricated MoS2 field effect transistors (FETs) with hBN as the encapsulating layer on SiO2/Si substrates, followed by depositing Au or Al pattern on the hBN (Figures 2a and 2b). Figure 2c shows the source-drain current (ISD) – gate voltage (VG) transfer curves measured before and after the deposition of the Au or Al pattern, where the highly doped Si substrate was used as a gate electrode. After depositing the Au pattern, the threshold voltage (Vth) shifted to a positive value. By contrast, the Al pattern shifted the Vth to a negative value (Figure 2d). These results corroborate the results of the SKPM analysis that is, the metals located underneath the hBN layer modulate the workfunction of MoS2. Next, we constructed a flash memory device by transferring the stacked layers of channel MoS2/upper hBN/underlying MoS2/lower hBN onto a SiO2 (300 nm)/Si substrate, as shown in Figure S1 (device I, Figure 3a). Figure 3d shows the ISD – VG transfer curve measured for device I by sweeping VG from -40 to +40 V, and back to -40 V. A large memory window (∆V) of ~20 V was observed, indicating that the underlying MoS2 in device I acts as a charge trapping layer. 22,23 To determine whether ∆V was affected by the thickness of the underlying MoS2, we fabricated various devices with different thicknesses of the underlying MoS2. However, there was no clear relationship between ∆V and the thickness of the underlying MoS2 (Figure S4). For comparison, we also fabricated a FET device on the hBN/SiO2/Si substrate. The ISD – VG transfer curve of this FET device showed no hysteresis (Figure 3e), confirming that the memory effects observed with device I were caused by the charge trapped in the underlying MoS2. Moreover, we noted that as the VG was swept, ISD changed more rapidly for device I than for the FET device. In other words, when the ISD was normalized by 6

ACS Paragon Plus Environment

Page 6 of 27

Page 7 of 27 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

ISD at VG=40 V and differentiated with respect to VG, the higher values of dISD/dVG were found for device I than for the FET device. To explore whether the characteristics of the flash memory device are affected by the workfunction modulation of the underlying MoS2 layer, we fabricated device II on an hBN/Au substrate (Figures 3b, 3c, and S1). The ISD of device II depended strongly on VG, implying that the Au layer did not screen significantly the electric field induced by VG (Figures 3f). In addition, when the ISD-VG transfer curves were measured in different ranges of VG, a smaller ∆V was observed for the narrower VG sweep range, which is consistent with conventional flash memory devices (Figure S5).31,32 However, when compared to device I, device II exhibited a larger ∆V (≈ 40 V) and higher maximum values in dISD/dVG. These results suggest that the workfunction of the underlying MoS2 layer was increased by the Au layer, potentially contributing to more rapid changes in the carrier concentration of the MoS2 channel upon application of VG. To investigate the data retention property of device II, we measured ISD at VSD = 0.1 V in the low resistance state (LRS) and high resistance state (HRS) as a function of time (Figure S6). No significant change in resistance magnitudes was observed within 104 s for both the LRS and HRS. In addition to device II, we fabricated device III, which incorporated Al instead of Au, since Al decreased the workfunction of the underlying MoS2 layer (Figure 3g). In contrast to device II, the ∆V of device III was smaller than the ∆V of device I, and the maximum values in dISD/dVG of device III were comparable to those of device I. These measurements indicate that the characteristics of MoS2–based flash memory devices might be dependent on the workfunction of the underlying MoS2 layer. To further investigate this hypothesis, we fabricated another device with Pd (device IV), whose workfunction (5.22-5.6 eV) is close to that of Au.30 Similarly to device II, device IV exhibited a larger ∆V and higher maximum values in dISD/dVG when compared to device I (Figure 3h). These observations further 7

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

corroborate the hypothesis that the characteristics of MoS2 flash memory devices are influenced by the workfunction modulation of the underlying MoS2 layer. Figure 3i shows the values of ∆V estimated for various devices with different thicknesses of the upper hBN layer. The ∆V increased with decreasing thickness of the upper hBN layer. A relationship between ∆V and the thickness of the lower hBN layer was also investigated; however, a clear relationship was not found (Figure S7). Figure 4a illustrates the energy band diagrams of device I. The energy bandgap (Eg) and electron affinity (χ) are approximately 1.2 eV and 4 eV for the multi-layer MoS2, 6 eV and 2 eV for hBN, and 9 eV and 0.9 eV for SiO2, respectively.33-35 Therefore, the barrier height for electron tunneling from the channel MoS2 (MoS2(1)) to the underlying MoS2 (MoS2(2)) through the upper hBN layer (BN(1)) (Φe1) is approximately 2 eV. When a high positive voltage is applied to VG (VG > 0), electrons can tunnel from MoS2(1) to MoS2(2) through BN(1), most likely via the Fowler-Nordheim tunneling,36 but cannot tunnel through BN(2) because of the SiO2. The resultant accumulation of electrons in MoS2(2) acts to screen the bottom-gate electric field that can reach the MoS2 channel, which results in a positive shift in the value of the Vth. However, when a negative voltage is applied to VG (VG < 0), electrons are transferred back from MoS2(2) to MoS2(1); simultaneously, holes can tunnel through BN(1) and become trapped in MoS2(2). As a result, the Vth shifts to a negative value, leading to ∆V. In addition, this electron tunneling between MoS2(1) and MoS2(2) might cause a larger variation in the carrier concentration of the channel MoS2 for device I than for the FET device, leading to a higher dISD/dVG. According to the results of SKPM analysis, the workfunction of the underlying MoS2 is expected to increase from ~4.4 eV to ~4.5 eV for device II and to decrease to ~4.32 eV for device III. Therefore, the barrier height for electron tunneling from MoS2(2) to Au or Al 8

ACS Paragon Plus Environment

Page 8 of 27

Page 9 of 27 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

through the lower hBN layer (BN(2)) (Φ e2) is estimated to be 1.9 eV for device II (Figure 4b) and 2.08 eV for device III, respectively (Figure 4c), although Φ e1 is 2 eV for both devices. In the case of device II, when VG > 0, electrons can tunnel through BN(1) as well as BN(2), and thus electrons are trapped in MoS2(2) as well as the Au film (Figure 4b). Consequently, the electric field that can reach MoS2(1) is more screened, which in turn shifts the Vth to a more positive value. When VG < 0, electrons are transferred from MoS2(2) and Au to MoS2(1), and holes are trapped in MoS2(2). As a result, device II yields larger ∆V and higher dISD/dVG when compared to device I, as shown in Figure 3. By contrast, the Φe2 of device III may be higher than the Φe2 of device II, and may increase further because of the oxide film (AlOx) formed on the Al surface (Figure 4c). Therefore, electron tunneling from MoS2 (2) to Al is most likely not permitted, resulting in a narrower ∆V for device III than for device II. Finally, we investigated whether the MoS2-based flash memory devices can be applied to synaptic devices. To mimic the potentiation and depression processes, the conductance change (∆S/S0) was measured for devices I (Figure 5a) and II (Figure 5b), while a series of 20 negative and positive voltage pulses with a width of 1 ms and a period of 1 s was applied to VG, where S0 is the initial conductance. The applied voltage amplitudes were -28 and +23 V for device I, and -12.5 and +24.2 V for device II, respectively. The application of negative (potentiating) voltage pulses resulted in a stepwise increase in ∆S/S0 for both devices, which was ascribed to a non-volatile element. However, as the voltage pulses were applied repeatedly, the ∆S/S0 of device I increased more slowly, while the ∆S/S0 of device II increased steadily. Furthermore, when the positive (depressing) voltage pulses were applied 20 times, the ∆S/S0 of device II decreased linearly down to the initial value, whereas the

∆S/S0 of device I decreased nonlinearly. Figure 5c illustrates the ∆S/S0 of device II measured 9

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

while a series of potentiating and depressing pulses were applied repeatedly. The measured

∆S/S0 exhibited linear and repeatable behaviors. Compared to the reported ion-gated synaptic transistors (2.5~5 V, 10 ms), 15-18 the higher voltage pulses were applied for the potentiation and depression process. However, the pulse width was shorter, so the power consumption was not higher than that of the reported ion-gated synaptic transistors. In fact, the 10 ~ 25 V requirement for device II was a consequence of the 300 nm-thick SiO2 gate oxide. If the thickness of SiO2 is reduced, the voltage pulse height and width can be significantly reduced, resulting in very low power consumption. The only difference between devices I and II is the presence of the Au layer under the lower hBN layer. It suggests that the presence of Au might contribute to improving the linearity and symmetry of bidirectional gradual conductance curves for device II. In the case of device II, charge transfer might occur between the underlying MoS2 and the Au layer, as well as between the channel MoS2 and the underlying MoS2 (Figure 4). As a result, when the voltage pulses were repeatedly applied, the stored charge on the underlying MoS2 might saturate more slowly for device II than for device I, improving the linearity of device II. To evaluate this possibility, we carried out similar measurements for device IV. Device IV exhibited bidirectional gradual conductance changes that were similar to those observed for to device II (Figure S8), supporting the hypothesis that the improved symmetry and linearity might be attributed to Au or Pd acting as another floating gate. These results demonstrate that the MoS2-based flash memory devices with double floating gates possibly simulate the potentiation and depression processes. To determine whether spike-rate-dependent plasticity (SRDP) can also be realized using device II, ∆S/S0 was measured after the application of a series of 10 identical negative voltage pulses with different delay times (∆t) to VG, in which the voltage amplitude was -29 V and 10

ACS Paragon Plus Environment

Page 10 of 27

Page 11 of 27 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

the pulse width was 1 ms (Figure 5d). The ∆S/S0 remained nearly unchanged for ∆t > 15 ms. However, for ∆t < 15 ms, the ∆S/S0 increased with decreasing ∆t, thereby indicating that short-term plasticity and long-term potentiation might be mimicked using device II.37 The effects of symmetric presynaptic and postsynaptic pulses do not cancel each other and the net effect is dominated by the timing difference between the presynaptic and postsynaptic effects, leading to the learning rule STDP.13,16 To realize STDP functionality, voltage pulses with an exponentially decaying function of time difference, td = tpre - tpost, was applied to VG, where tpre and tpost are the time when the presynaptic and postsynaptic pulses spike, respectively. Figure 5e shows the ∆S/S0 measured as a function of td. When the presynaptic spike occurred before the postsynaptic spike (td > 0, potentiation), the conductance (synaptic weight) decreased with increasing ∆t. When the postsynaptic spike occurred prior to the presynaptic spike (td < 0, depression), reverse modulation was observed. In addition, the plot of ∆S/S0 versus td could be well fitted with exponential decay functions, demonstrating that neural synaptic STDP behavior may be achieved using device II.13

Conclusion SKPM analysis showed that the workfunction of MoS2 placed on hBN/metal depends on the metal. In addition, this workfunction modulation of MoS2 through the hBN layer was supported by ISD - VG transfer curves measured for MoS2-based FET devices. To investigate whether the characteristics of MoS2-based flash memory device are affected by the workfunction modulation of MoS2 acting as a floating gate, we fabricated MoS2-based flash memory devices on hBN, hBN/Au, and hBN/Al substrates. The device fabricated on hBN/Au yielded larger memory windows than the device prepared on hBN, whereas the device fabricated on hBN/Al showed comparable memory 11

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

windows to the device on hBN. Furthermore, the symmetry and linearity of the bidirectional gradual conductance change curves improved significantly with the device fabricated on the hBN/Au substrate compared to the devices prepared on hBN or hBN/Al. To explain this, we proposed energy band diagrams that show how Au acts as another floating gate, while Al does not, owing to the oxide film that forms on the Al surface. For flash memory devices with double floating gates, electron transfers occur between the underlying MoS2 and Au, as well as between the underlying MoS2 and the channel MoS2. Thus, electron tunneling and injection may be more controllable, leading to the symmetry of potentiation and depression processes. This hypothesis was investigated by fabricating a device on an hBN/Pd substrate. This device exhibited behavior similar to that of the device fabricated on hBN/Au. In addition, various synaptic functions, including potentiation and depression processes, SRDP, and STDP, could be emulated using flash memory devices fabricated on hBN/Au, demonstrating that MoS2-based flash devices with double floating gates may be applied to synaptic emulators. To show proof- of-concept, we fabricated devices using the pick-up method. However, if large scale MoS2 and hBN layers are available, the synaptic array may be constructed using conventional lithography techniques.

Methods Device fabrication The flash memory devices were fabricated as shown in Figure S1. First, MoS2 and hBN flakes were mechanically exfoliated from MoS2 (SPI supply) and hBN (HQ graphene) single crystals onto a SiO2/Si substrate. Subsequently, they were transferred to another substrate using a hot pick-up method with a polypropylene carbonate film (PPC, Sigma-Aldrich, CAS 25511-85-7).38 Briefly, 1-µm thick PPC film was spin-coated onto a SiO2/Si substrate and 12

ACS Paragon Plus Environment

Page 12 of 27

Page 13 of 27 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

peeled off from the substrate using a scotch-tape with a center hole. The PPC film was transferred onto a polydimethylsiloxane (PDMS) block on a slide glass, followed by heating at 120 °C for 1 min to enhance the adhesion of PPC to PDMS. The hBN flake prepared on the SiO2/Si substrate was picked up using the PPC film on the PDMS block, and the MoS2 flake was picked up by placing the hBN layer on the MoS2 flake. In the same way, the hBN flake was picked up, affording an hBN/MoS2/hBN stacked layer. This layer was released onto a metal/SiO2/Si or SiO2/Si substrate at 180°C by dissolving the PPC film in chloroform, where metals such as Au (30 nm), Al (30 nm), and Pd (30 nm) were deposited using an electron-beam evaporator. Next, the channel MoS2 flake was transferred onto the stacked sheets (hBN/MoS2/hBN) using the PDMS stamping method to avoid contamination. Finally, Cr (3 nm)/Au (30 nm) electrodes were made using electron-beam lithography and lift-off techniques. The thicknesses of the MoS2 and hBN flakes measured using AFM were approximately 10 ~14 nm and 10 ~ 34 nm, respectively, as summarized in Table S1. Contact potential difference (CPD) and electrical measurements The CPD was measured using a scanning Kelvin probe microscope (NX10, Parks System, Korea) at ambient conditions. The tip was a conductive Si tip coated with Au (PPP-NCHAu AFM, Nanosensors). Electrical measurements were carried out using a semiconductor characterization system (4200-SCS, Keithley). To measure the synaptic functions, voltage pulses were applied using an ultra-fast I-V module for 4200-SCS (4225-PMU, Keithley).

13

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1. Workfunction of MoS2 modulated by Au and Al. (a) Schematic diagram of the specimen prepared by stacking metals (Au and Al), hBN, and MoS2 layers for SKPM measurements. (b) SKPM image of the specimen prepared by stacking metals (Au and Al), hBN, and MoS2 layers. The topography image measured simultaneously with SKPM is shown in Figure S2. The hBN layer is denoted with white dotted curves. (c) The x-directional line profiles of VCPD obtained for uncovered Au and Al with hBN and MoS2 layers (green line in b) and for the MoS2 layer on hBN layer and Au (or Al) (red line in b). (d) The ydirectional line profiles of VCPD obtained for the MoS2 layer on hBN layer and Au (brown line in b), for the MoS2 layer on hBN layer (blue line in b), and for the MoS2 layer on hBN layer and Al (pink line in b).

14

ACS Paragon Plus Environment

Page 14 of 27

Page 15 of 27 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Figure 2. Effects of metals on the MoS2 FET devices through the hBN layer. (a) Schematic diagram of the MoS2 FET device with hBN layer and Au (or Al). (b) Optical image of the fabricated device. (c) The ISD-VG transfer curves measured at VSD=0.1 V before and after depositing Au on the hBN layer. (d) The ISD-VG transfer curves measured at VSD=0.1 V before and after depositing Al on the hBN layer.

15

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3. Memory effects of various MoS2-based flash memory devices. (a) Schematic diagram of a MoS2-based flash memory device fabricated on SiO2/Si substrate (device I). (b) Schematic diagram of a MoS2-based flash memory device fabricated on metal/SiO2/Si substrate. Devices II, III, and IV are prepared on Au, Al, and Pd/SiO2/Si substrates, respectively. (c) An optical image of device II. (d-h) The ISD-VG transfer curves (black symbols) and dISD/dVG (red curves) for (d) device I, (e) FET, (f) device II, (g) device III, and (h) device IV. The ISD-VG transfer curves were obtained by sweeping VG from -40 to +40 V and back to -40 V, and ISD was normalized by ISD at VG=40 V for comparison among different devices. (i) Memory windows (∆V) versus the thickness of the upper hBN layer estimated for various MoS2-based flash memory devices. The filled and open symbols represent ∆V estimated from the ISD-VG transfer curves in the VG sweep range of -40 ~ 40 V and -60 ~ 60 V, respectively. 16

ACS Paragon Plus Environment

Page 16 of 27

Page 17 of 27 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Figure 4. The energy band diagrams of devices I, II, and III. (a) The energy band diagrams in a flat band state, equilibrium, VG > 0, and VG < 0 for device I. EC, EV, Eg, χ, ϕ, and Φ1 are the conduction band energy, valence band energy, energy bandgap, electron 17

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

affinity, workfunction, and tunneling barrier height from the MoS2(1) to MoS2(2) through BN(1). For VG > 0, electrons tunnel from MoS2(1) to MoS2(2) through BN(1) and are subsequently trapped in MoS2(2). For VG < 0, electrons tunnel from MoS2(2) to MoS2(1) through BN(1), while holes tunnel from MoS2(1) to MoS2(2). (b) The band diagrams in a flat band state, equilibrium, VG > 0, and VG < 0 for device II. Φ2 is the tunneling barrier height from MoS2(2) to Au through BN(2). For VG > 0, electrons tunnel from MoS2(1) to MoS2(2) through BN(1), as well as from MoS2(2) to Au through BN(2). For VG < 0, electrons tunnel from Au to MoS2(2) through BN(2), as well as from MoS2(2) to MoS2(1) through BN(1). (c) The band diagrams in a flat band state, equilibrium, VG > 0, and VG < 0 for device III. AlOx forms on the Al surface naturally, leading to an increase in Φ2 and prohibiting electron transfer between Al and MoS2.

18

ACS Paragon Plus Environment

Page 18 of 27

Page 19 of 27 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Figure 5. Simulation of synaptic functions. (a) Conductance change measured at ISD=0.1 V for device I while repeated potentiation (-28 V for 1 ms, spaced 1 s apart) and depression (23 V for 1 ms, spaced 1 s apart) pulses are applied to VG. The conductance change was normalized by the initial conductance (∆S/S0). (b) ∆S/S0 measured at ISD=0.1 V for device II while repeated potentiation (-12.5 V for 1 ms, spaced 1 s apart) and depression (24.2 V for 1 ms, spaced 1 s apart) pulses are applied to VG. (c) ∆S/S0 measured at ISD=0.1 V for device II while potentiation and depression gate voltages are applied repeatedly. (d) ∆S/S0 measured at ISD=0.1 V for device II after gate pulses (-29 V for 1 ms) with various delay times (∆t) are applied 10 times, demonstrating spike-rate dependent plasticity (SRDP). (e) ∆S/S0 measured at ISD=0.1 V for device II after presynaptic (Vpeak = 35 V) and postsynaptic pulses (Vpeak = -35 V) are applied with a time difference td, as shown in the left inset, demonstrating spike-timing dependent plasticity (STDP). The source and drain are analogues of the preneuron and postneuron terminals, respectively. The red curves are an exponential decay function fitted to the data. 19

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

20

ACS Paragon Plus Environment

Page 20 of 27

Page 21 of 27 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Additional information Supplementary information is available in the online version of the paper. Reprints and permissions information is available online at http://pubs.acs.org. Correspondence and requests for materials should be addressed to K.Y. *Corresponding author e-mail: K-H Yoo, [email protected] Author contributions K.Y., S.Y., and G.L. conceived and designed the study. S.Y. J.K. M.P. and C.L. fabricated devices and performed the electrical measurements. S.K. performed the SKPM measurements. K.Y. and S.Y. wrote the manuscript. All authors discussed the results and commented on the manuscript. Competing financial interests The authors declare that they have no competing interests. Acknowledgements This work has been financially supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT and Future Planning (Grant Nos. 2016R1A2B3011980).

21

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

References 1. Cochocki, A., Unbehauen, R. Neural Networks for Optimization and Signal Processing; Wiley: New York, 1993. 2. Brette, R.; Rudolph, M.; Carnevale, T.; Hines, M.; Beeman, D.; Bower, J. M.; Diesmann, M.; Morrison, A.; Goodman, P. H.; Harris, F. C.; Zirpe, M.; Natschläger, T.; Pecevski, D.; Ermentrout, B.; Djurfeldt, M.; Lansner, A.; Rochel, O.; Vieville, T.; Muller, E.; Davison, A.P.; Boustani, S.E.; Destexhe, A. Simulation of Networks of Spiking Neurons: a Review of Tools and Strategies. J. Comput. Neurosci. 2007, 23, 349-398. 3. Prezioso, M.; Merrikh-Bayat, F.; Hoskins, B. D.; Adam, G. C.; Likharev, K. K.; Strukov, D. B. Training and Operation of an Integrated Neuromorphic Network Based on Metal-Oxide Memristors. Nature 2015, 521, 61-64. 4. Yang, J. J.; Strukov, D. B.; Stewart, D. R. Memristive Devices for Computing. Nat. Nanotechnol. 2013, 8, 13-24. 5. Carve, M. Analog VLSI and Neural Systems; Addison Wesley: Boston, 1989. 6. Mahowald, M.; Douglas, R. A Silicon Neuron. Nature 1991, 354, 515-518. 7. Morie, T.; Fujita, O.; Uchimura, K. Self-Learning Analog Neural Network LSI with HighResolution Non-Volatile Analog Memory and a Partially-Serial Weight-Update Architecture. IEICE T. Electron. 1997, E80c, 990-995. 8. Sun, J. CMOS and Memristor Technologies for Neuromorphic Computing Applications; http://www.eecs.berkeley.edu/Pubs/TechRpts/2015/EECS-2015-219.html 2015.

22

ACS Paragon Plus Environment

Page 22 of 27

Page 23 of 27 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

9. Calhoun, B.H.; Wang, A.; Chandrakasan, A. Modeling and Sizing for Minimum Energy Operation in Subthreshold Circuits. IEEE J. Solid-St. Circ. 2005, 40, 1778-1786. 10. Ohno, T.; Hasegawa, T.; Tsuruoka, T.; Terabe, K.; Gimzewski, J. K.; Aono, M. ShortTerm Plasticity and Long-Term Potentiation Mimicked in Single Inorganic Synapses. Nat. Mater. 2011, 10, 591-595. 11. Jo, S.H.; Chang, T.; Ebong, I.; Bhadviya, B.B.; Mazumder, P.; Lu, W. Nanoscale Memristor Device as Synapse in Neuromorphic Systems. Nano Lett. 2010, 10, 1297-1301. 12. Chang, T.; Jo, S. H.; Lu, W. Short-Term Memory to Long-Term Memory Transition in a Nanoscale Memristor. ACS Nano 2011, 5, 7669-7676. 13. Yan, X. B.; Zhao, J. H.; Liu, S.; Zhou, Z. Y.; Liu, Q.; Chen, J. S.; Liu, X. Y. Memristor with Ag-Cluster-Doped TiO2 Films as Artificial Synapse for Neuroinspired Computing. Adv. Funct. Mater. 2018, 28, 1705320. 14. Li, Y.; Zhong, Y. P.; Xu, L.; Zhang, J. J.; Xu, X. H.; Sun, H. J.; Miao, X. S. Ultrafast Synaptic Events in a Chalcogenide Memristor. Sci. Rep. 2013, 3, 1619. 15. Nishitani, Y.; Kaneko, Y.; Ueda, M.; Morie, T.; Fujii, E. Three-Terminal Ferroelectric Synapse Device with Concurrent Learning Function for Artificial Neural Networks. J. Appl. Phys. 2012, 111, 124108. 15. Shi, J.; Ha, S. D.; Zhou, Y.; Schoofs, F.; Ramanathan, S. A Correlated Nickelate Synaptic Transistor. Nat. Commun. 2013, 4, 2676.

23

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

16. Jiang, J.; Guo, J.; Wan, X.; Yang, Y.; Xie, H.; Niu, D.; Yang, J.; He, J.; Gao, Y.; Wan, Q. 2D MoS2 Neuromorphic Devices for Brain-Like Computational Systems. Small 2017, 13, 1700933. 17. Shen, A. M.; Chen, C. L.; Kim, K.; Cho, B.; Tudor, A.; Chen, Y. Analog Neuromorphic Module Based on Carbon Nanotube Synapses. ACS Nano 2013, 7, 6117-6122. 18. Zhu, J.; Yang, Y.; Jia, R.; Liang, Z.; Zhu, W.; Rehman, Z. U.; Bao, L.; Zhang, X.; Cai, Y.; Song, L.; Huang, R. Ion Gated Synaptic Transistors Based on 2D van der Waals Crystals with Tunable Diffusive Dynamics. Adv. Mater. 2018, 30, 1800195.

19. Dai, S.; Wu, X.; Liu, D.; Chu, Y.; Wang, K.; Yang, B.; Huang, J. Light-Stimulated Synaptic Devices Utilizing Interfacial Effect of Organic Field-Effect Transistors. ACS Appl. Mater. Inter. 2018, 10, 21472−21480

20. Novoselov, K. S.; Mishchenko, A.; Carvalho, A.; Castro Neto, A. H. 2D Materials and van der Waals Heterostructures. Science 2016, 353, aac9439. 21. Fiori, G.; Bonaccorso, F.; Iannaccone, G.; Palacios, T.; Neumaier, D.; Seabaugh, A.; Banerjee, S. K.; Colombo, L. Electronics Based on Two-Dimensional Materials. Nat. Nanotechnol. 2014, 9, 768-779. 22. Choi, M. S.; Lee, G. H.; Yu, Y. J.; Lee, D. Y.; Lee, S. H.; Kim, P.; Hone, J.; Yoo, W. J. Controlled Charge Trapping by Molybdenum Disulphide and Graphene in Ultrathin Heterostructured Memory Devices. Nat. Commun. 2013, 4, 1624. 23. Vu, Q. A.; Shin, Y. S.; Kim, Y. R.; Nguyen, V. L.; Kang, W. T.; Kim, H.; Luong, D. H.; Lee, I. M.; Lee, K.; Ko, D. S.; Heo, J.; Park, S.; Lee, Y. H.; Yu, W. J. Two-Terminal 24

ACS Paragon Plus Environment

Page 24 of 27

Page 25 of 27 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Floating-Gate Memory with van der Waals Heterostructures for Ultrahigh On/off Ratio. Nat. Commun. 2016, 7, 12725. 24. Diorio, C; Hasler, P; Minch, B. A.; Mead, C. A. A Single-Transistor Silicon Synapse. IEEE T. Electron. Dev. 1996, 43, 1972-1980. 25. Yu, Y. J.; Zhao, Y.; Ryu, S.; Brus, L. E.; Kim, K. S.; Kim, P. Tuning the Graphene Work Function by Electric Field Effect. Nano Lett. 2009, 9, 3430-3434. 26. Smith, C. E.; Xie, Z. T.; Baldea, I.; Frisbie, C. D. Work Function and Temperature Dependence of Electron Tunneling Through an N-type Perylene Diimide Molecular Junction with Isocyanide Surface Linkers. Nanoscale 2018, 10, 964-975. 27. Shi, Y.; Kim, K. K.; Reina, A.; Hofmann, M.; Li, L. J.; Kong, J. Work Function Engineering of Graphene Electrode via Chemical Doping. ACS Nano 2010, 4, 2689-2694. 28. Lee, S. Y.; Kim, U. J.; Chung, J.; Nam, H.; Jeong, H. Y.; Han, G. H.; Kim, H.; Oh, H. M.; Lee, H.; Kim, H.; Roh, Y. G.; Kim, J.; Hwang, S. W.; Park, Y.; Lee, Y. H. Large Work Function Modulation of Monolayer MoS2 by Ambient Gases. ACS Nano 2016, 10, 61006107. 29. Woo, J.; Moon, K.; Song, J.; Kwak, M.; Park, J.; Hwang, H. Optimized Programming Scheme Enabling Linear Potentiation in Filamentary HfO2 RRAM Synapse for Neuromorphic Systems. IEEE T. Electron. Dev. 2016, 63, 5064-5067. 30. Rumble, J. CRC Handbook of Chemistry and Physics; Taylor & Francis: U.K, 2017.

25

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

31. Lin, C. C.; Wu, Y. H.; Lin, Y. S.; Wu, M. L.; Chen, L. L. Electrical Characteristics for Flash Memory with Germanium Nitride as the Charge-Trapping Layer. IEEE T. Nanotechnol. 2013, 12, 436-441. 32. Joe, S. M.; Lee, J. H. Investigation of Hysteresis Phenomenon in Floating-Gate NAND Flash Memory Cells. IEEE T. Electron. Dev. 2015, 62, 2738-2744. 33. Ganatra, R.; Zhang, Q. Few-Layer MoS2: A Promising Layered Semiconductor. ACS Nano 2014, 8, 4074-4099. 34. Cassabois, G.; Valvin, P.; Gil, B. Hexagonal Boron Nitride is an Indirect Bandgap Semiconductor. Nat. Photonics 2016, 10, 262-267. 35. El-Kareh, Badih. Fundamentals of Semiconductor Processing Technology; Springer Verlag: Germany, 1995. 36. Sze, S.M., Ng, Kwok K. Physics of Semiconductor Devices; Wiley: New York, 2007. 37. Wang, Z. R.; Joshi, S.; Savel'ev, S. E.; Jiang, H.; Midya, R.; Lin, P.; Hu, M.; Ge, N.; Strachan, J. P.; Li, Z. Y.; Wu, Q.; Barne, M.; Li, G. L.; Xin, H. L.; Williams, R. S.; Xia, Q. F.; Yang, J. J. Memristors with Diffusive Dynamics as Synaptic Emulators for Neuromorphic Computing. Nat. Mater. 2017, 16, 101-108. 38. Pizzocchero, F.; Gammelgaard, L.; Jessen, B. S.; Caridad, J. M.; Wang, L.; Hone, J.; Boggild, P.; Booth, T. J. The Hot Pick-Up Technique for Batch Assembly of van der Waals Heterostructures. Nat. Commun. 2016, 7, 11894.

26

ACS Paragon Plus Environment

Page 26 of 27

Page 27 of 27 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

SYNOPSIS TOC

27

ACS Paragon Plus Environment