A generic method to control hysteresis and memory effect in Van der Waals hybrids

The diverse properties of two-dimensional materials have been utilized in a variety of architecture to fabricate high quality electronic circuit elements. Here we demonstrate a generic method to control hysteresis and stable memory effect in Van der Waals hybrids with a floating gate as the base layer. The floating gate can be charged with a global back gate-voltage, which it can retain in a stable manner. Such devices can provide a very high, leakage-free effective gate-voltage on the field-effect transistors due to effective capacitance amplification, which also leads to reduced input power requirements on electronic devices. The capacitance amplification factor of ∼10 can be further enhanced by increasing the area of the floating gate. We have exploited this method to achieve highly durable memory action multiple genre of ultra-thin 2D channels, including graphene, MoS2, and topological insulators at room temperature.

One major area of thrust with Van der Waals hybrids has been to achieve non-volatile, durable memory action, which provide several advantages such as reduced power consumption, higher integration density, and enhanced life-times [43,44]. In this work, we have demonstrated a generic device architecture to achieve nonvolatile memory effect in stacks made using 2D materials in a floating gate geometry [45][46][47][48]. While the floating gate geometry has been exploited in case of transition metal dichalcogenide-based in 2D atomic/molecular layers, it's applicability for a broader range of 2D materials, and an optimization of the architecture itself, remains unexplored. The floating gate can be charged by the application of a back-gate voltage(V G ). The charge remains stored even after the gate-voltage is switched off due to lower back transfer probability, and gives rise to higher effective gate-voltage compared to that of the standard 285nm SiO 2 . The cyclic charge trapping due to V G -dependent band-bending leads to a large hysteresis window, which can be accurately tuned by controlling the area of the floating gate. This effect can be utilized in a variety of devices such as extremely robust low power, Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. non-volatile memory applications using MoS 2 , which is difficult to achieve with graphene due to the absence of any band-gap. The high effective V G can also be implemented to achieve ambipolar transport in topological insulators, and highly repeatable large hysteresis window, even at room temperature. Attaining ambipolar transport usually requires a large gate-voltage, and low temperature due to high residual bulk doping. Such devices can be utilized for exploring several properties of Dirac Fermions in TIs when access to the surface states is required.
The durable memory action in our devices is made possible due to the presence of the floating gate. Floating gate MOSFETs (FGMOS) have been extensively used as in memory applications such as EPROMs (Erasable Programmable Read Only Memories) and EEPROMSs (Electrically Erasable Programmable Read Only Memories), and can be used in neural networks as well [46,49,50]. With substantial improvement of the fabrication techniques and the emergence of two-dimensional materials, such circuits can be implemented using entirely Van der Waals materials, leading to further miniaturization. The advantage of using a floating gate rests in the amplification of the equivalent capacitance (C eqv ), and hence the effective V G on the channel, which is our case, is graphene, MoS 2 , and TIs. A typical schematic to understand the capacitance amplification in the present scenario is illustrated in figure 1(a) . Let C 1 be the capacitance between the graphene channel and the 2DM underlayer. The thickness of this capacitor is equal to the thickness of hBN(d hBN ), which also acts as the tunnel barrier in these devices. The area of the capacitor is equal to the area (A) of the channel. Let C 2 be the capacitance between the Si back gate and the 2DM underlayer, whose thickness is equal to the thickness of SiO 2 (d S ) , and the area is equal to the area of the underlayer(A 2D ) . Let us now assume that the area of the underlayer is equal to the area of the channel (A 2DM =A), which is illustrated in figure 1(a). Since C 1 and C 2 both have the same area and the dielectric constant of hBN and SiO 2 are also similar(ò r ∼4), we can write = C nC 1 2 , where n=d S /d hBN . If V G is applied on the Si gate, and V 0 is the voltage at the underlayer, since the charge induced on the underlayer and the channel is same, we can write In real devices, however, the area of the underlayer is always larger, and can be expressed as A 2D =A+xA. Here A is the part of the underlayer that is lying underneath the channel, and xA is the extended part of the underlayer located outside the channel area. In such a scenario, the capacitance between the Si gate and the underlayer ( ¢ C 2 ) can be expressed as ¢ = + C C xC 2 2 2 , which is schematically represented in figure 1(b). This implies that the voltage at the underlayer V A can be expressed as, ⎛ The capacitance amplification (A V ) in our devices is defined as the ratio of V A and V 0 . It can be expressed as Using the relation = C kC 1 2 and considering the fact that both > > k x , 1 , equation (3) can be rewritten as clearly enunciates that the effective V G can be further enhanced by the increasing the area of the underlayer, a property which has been exploited in our present work. The devices studied in this manuscript have been fabricated from several two dimensional materials such as single layer graphene(SLG), molybdenum disulphide(MoS 2 ), Bi-based chalcogenide Bi 1.6 Sb 0.4 Te 2 Se topological insulators (TIs), and boron nitride. While the top channel is formed by graphene, MoS 2 , and TIs, the floating gate is primarily made of single or multi-layer graphene, and MoS 2 . In certain devices, the area of the floating gate has been extended by depositing 20nm Au, which is connected only to the 2D underlayer only. In all our devices, hBN forms the tunnel barrier between the channel and the floating gate. In order to process these devices, graphene, MoS 2 , Bi 1.6 Sb 0.4 Te 2 Se, and hBN were separately exfoliated on standard 285nm SiO 2 /Si substrate, followed by the fabrication of the required hetero-structure in a home-made transfer assembly. The contacts were patterned by standard electron beam lithography followed by thermal evaporation of 5/50 nm Cr/Au for source, and drain electrodes. All resistivity measurements were performed using a lock-in amplifier with a carrier frequency of 227Hz, in a home-made variable temperature cryostat.
The preliminary transfer characteristics of a typical device using the floating gate geometry with single layer graphene as the top channel(figure 2(a)) is shown in figure 2(b), which exhibits two distinct features. Firstly, the transport in the graphene channel with 2DM underlayer displays a strong hysteresis between the forward and backward V G sweeps, while the standard hysteresis-free R vs V G is restored when the 2D underlayer is absent. Secondly, the sharpness of the transfer characteristics is significantly enhanced in the presence of the 2D layer, suggesting an increase in either the carrier mobility in the channel or the effective capacitance. Since hBN is the common substrate to the entire channel, both effects can only be due to the presence of the 2D underlayer. Both these effects, however, are absent when the 2D underlayer is electrically grounded, and the transfer characteristics is determined by the extent of V G -dependent screening by the doped 2D layer(see supplementary section A available online at stacks.iop.org/MRX/7/014004/mmedia). The extent of (anti-)hysteresis can be reduced, and eventually completely removed, by decreasing the sweep range below a critical value V CG (±12V for the device in figure 2(c)), which suggests that the total trapped charge in the 2DM layer remains fixed as long as the sweep range is below V CG around the Dirac point. The enhanced sharpness of the transfer characteristics in the presence of the 2DM underlayer, however, has negligible effect on the range of the V G sweep.
To understand it's origin, we have extracted the carrier density n in the channel directly from Hall measurements. The measured n for forward and backward sweeps over a sweep range of 60V is shown in figure 2(d). The dependence of n on V G is also (anti-)hysteretic, with the Dirac points coinciding with the reversal in the sign of the carriers. Around the Dirac points n varies linearly with V G , and the total capacitance(C G ) sis obtained from the slope. Evidently, C G is considerably larger than C 0 (slope of the solid line in figure 2 where d BN and d SiO 2 are the thicknesses of the hBN and SiO 2 layers, respectively. C 2DM is usually 10-100 times larger that C BN and C SiO 2 , and can be neglected. This confirms that the sharpness in the R−V G characteristics is indeed is due to enhanced C G , rather than an improbable increase in the carrier mobility in the channel region just above the 2DM underlayer. By restricting the sweep range below V CG , the variation of n with V G can be made completely non-hysteretic, irrespective of the location of the Dirac point, i.e. the extent of trapped charge in the 2DM underlayer. Figure 2(e) illustrates this at three locations of the Dirac point.
The clockwise hysteresis (i.e. anti-hysteresis) of R−V G in figure 2(f) is the effect of cyclic charge trapping by the 2DM due to V G -dependent bending of the hBN bands, and has been recently observed in memory devices based on trilayer Van der Waals hybrid of similar architecture [45,47]. When V G increases in positive or negative directions from the Dirac point, the transfer of electrons or holes from the graphene channel to the 2DM causes an effective gate voltage shift, that can be erased only by reversing the polarity of V G . The enhancement in V A is limited by the band gap Δ BN of hBN, and for voltage | | > D V D B N , electric field-induced band bending within hBN will cause charge to flow across the dielectric, thereby pinning the Fermi level of the channel. This allows us to estimate both D~D V C C G C BN BN G , as well as the saturation carrier density n max ∼ε 0 ε BN Δ BN /ed BN , both of which could be directly verified from their dependence on d BN . Further amplification of the effective capacitance can be achieved by enhancing the area of the floating gate, which can be done by connecting the underlayer with Au. Using this, A V has been found to increase by an order of magnitude (See Supplementary section B).
One of the shortcomings of graphene is the lack of a band-gap, rendering it unsuitable for non-volatile memory based applications. To circumvent this, we have utilized the floating gate technique in an MoS 2 FET to evince the possibility of high performance and power efficient non-volatile memory devices. The device structure is shown in (figure 3(a)) where the graphene and molybdenum disulphide layer have been interchanged from ( figure 2(a)). Similar hysteresis is also observed when MoS 2 forms the top current carrying layer as demonstrated in figure 3(b) due to the screening of back gate electric field by the trapped charges in the floating gate. Such structures have been previously investigated in [45,47,51]. One of the major parameters that determine electronic device performances is the sub-threshold swing. It is defined as the amount of gate voltage bias required to change the source-drain current(I sd ) by one decade, and is given by the following expression, 10 10 Here S is the surface potential in the channel. The second term in the equation (5) has Boltzman theoretical limit of 60mV/decade at room temperature while the first term, known as the body factor is given by where C S is the quantum capacitance of the channel. The equivalent capacitance for our extended floating gate structure is the capacitance of the hBN dielectric, which is much larger than that of the 285nm SiO 2 . From equation (6) it is clear that higher the equivalent capacitance, lower is the value of sub-threshold swing (S) leading to faster switching, and reduced power consumption. The higher capacitance per unit area of hBN allows us to reduce the body factor significantly leading to almost ideal subthreshold slope in all the measured devices for over four decades of I sd as shown in the inset of figure 3(b). The large memory window and the near ideal subthreshold swing ( figure 3(b)) make the extended floating gate MoS 2 devices perfect candidates for non-volatile memory storage. The binary memory action is depicted in figure 3(c) with a ON/OFF ratio of ≈10 3 at a drain bias of 50mV, with a pulse time of 1 section. The pulse rise time has a significant effect on the on/off ratio, which is expected due to increased storage of charges on the floating gate, and is depicted for pulse widths of 0.1, 1, and 10msec (Supplementary section C.1). The devices also show excellent retentivity, with negligible change of I SD upto ∼10000 secs for all three pulse times (Supplementary section C.1). To obtain the switching transition, a program and erase pulse of +9 and −9 V respectively was applied to the silicon back gate. The memory action is robust and repeatability has been checked for 100 program and erase cycles ( figure 3(d)). We find a very low cycle to cycle variability in the memory ratio, defined as the ratio of the on and off state current, demonstrating the reliability of these non-volatile memory devices (inset of figure 3(d)). Furthermore, the memory ratio in our devices, which can be as high as 10 8 , is comparable to state of the art sensors fabricated with two-dimensional materials, that have reported in literature (Supplementary section C.2) [52][53][54][55][56][57][58][59]. Additionally, the low bias requirements and program state current value makes the device power efficient with an observed power dissipation of ≈10pW in the program state in our MoS 2 based memory devices [51]. A lower stress on the gate dielectric due to reduced program and magnitude improves the reliability of these devices making them ideal for low power portable memory applications.
We have further exploited this technique by using topological insulators(TIs) as the conducting channel. High bulk conductivity due to doping from defects in the bulk makes it difficult to tune the Fermi level across the Dirac point, and achieve ambipolar transport. This severely limits the T range for practical applications of TIs since the surface transport regime remains inaccessible. This problem can be circumnavigated either by going to very low temperatures, where the bulk doping is suppressed, which is not a feasible solution for several practical applications, or by using an ionic liquid top gate, can degrade the device characteristics due to intercalation of ions [60][61][62][63][64][65]. To address this, we have utilized capacitance amplification to access the Dirac Fermion dominated transport regime, even at high T. To fabricate the devices, the TI Bi 1.6 Sb 0.4 Te 2 Se was used, which offers reduced bulk number density due to compensation doping, and the transport below T<50 K for samples with thickness d<100 nm is primarily through the surface [66][67][68]. The scanning electron micrograph of a typical device with MoS 2 as the floating gate is shown in figure 4(a). The R−V G as shown in figure 4(b) shows ambipolar transport at T=22 K, and remarkably at T=300 K as well. Similar transfer characteristics were observed in devices with multi-layer graphene as the floating gate(Supplementary section D.1). Such large and controlled hysteresis in TI-FETs clearly demonstrate that they can be used in non-volatile memory operations at room temperature. The percentage change in R, however, is ∼10%-20%, which implies that only the bottom surface is affected by the V G , and the top surface remains unaffected due to screening of V G by the bottom surface and the bulk. The measured R is possibly limited by the top surface, and any residual bulk conducting channels. To verify the ambipolar nature of transport, we have performed Hall measurements in the sample at different V G -s across the , and is equal to 7.5×10 −4 F, which is ∼6 times larger than the 285nm SiO 2 dielectric. The R-V G data for several sweep rates clearly indicates that the hysteresis is independent of the sweep rate(Supplementary section D.2). From R−V G at different temperatures, the impurity density as well as mobility of the TIs can be extracted from T=22-300K (Supplementary section D.3). The mobility of the sample increases as the T is reduced, as expected. However, the value of mobility is low(< 50 cm 2 V −1 S −1 ) at low T, even though the samples are on BN. This is in stark contrast to a medium like graphene, since μ increases manifold when the SiO 2 substrate is replaced by the atomically flat BN, where the effect of trapping-detrapping and dangling bonds on the electrical transport in the sample is significantly reduced. This clearly signifies that the low mobility in TI samples can be attributed to the contribution of bulk defects which act as the dominant scattering centers [68][69][70][71], and provide an intrinsic limit to the mobility.
In summary, we have utilized capacitance amplification due to a floating gate in Van der Waals hybrids fabricated out of graphene, MoS 2 , and topological insulators. The increase in effective capacitance leads to an increase in the effective gate-voltage on the conducting channel, which shows a large hysteresis due to cyclic trapping and detrapping of charges. This effect has been used in a variety of ultra-thin materials like graphene, MoS 2 , and topological insulators to demonstrate highly stable, non-volatile memory response at room temperature. Our work provides a framework for employing hybrids of atomically thin quantum materials for memory-based applications requiring high durability, and low power .