nanomaterials

Article Experimental Studies on the Dynamic Memcapacitance Modulation of the ReO3@ReS2 Composite Material-Based Diode

Joanna Borowiec 1,2,* , Mengren Liu 3 , Weizheng Liang 4, Theo Kreouzis 2 , Adrian J. Bevan 2 , Yi He 1, Yao Ma 1 and William P. Gillin 1,2 1 College of Physics, Sichuan University, Chengdu 610064, China; [email protected] (Y.H.); [email protected] (Y.M.); [email protected] (W.P.G.) 2 Materials Research Institute and School of Physics and Astronomy, Queen Mary University of London, Mile End Road, London E1 4NS, UK; [email protected] (T.K.); [email protected] (A.J.B.) 3 Sichuan University—Pittsburgh Institute, Chengdu 610207, China; [email protected] 4 The Peac Institute of Multiscale Sciences, Chengdu 610031, China; [email protected] * Correspondence: [email protected]; Tel.: +86-028-854-12323

 Received: 4 October 2020; Accepted: 19 October 2020; Published: 23 October 2020 

Abstract: In this study, both memcapacitive and memristive characteristics in the composite material based on the disulfide (ReS2) rich in rhenium (VI) (ReO3) surface overlayer (ReO3@ReS2) and in the oxide (ITO)/ReO3@ReS2/aluminum (Al) device configuration is presented. Comprehensive experimental analysis of the ReO3@ReS2 material properties’ dependence on the memcapacitor electrical characteristics was carried out by standard as well as frequency-dependent current–voltage, capacitance–voltage, and conductance–voltage studies. Furthermore, determination of the charge carrier conduction model, charge carrier mobility, of the trap states, density of the available charge carrier, free-carrier concentration, effective density of states in the conduction band, activation energy of the carrier transport, as well as ion hopping was successfully conducted for the ReO3@ReS2 based on the experimental data. The ITO/ReO3@ReS2/Al charge carrier conduction was found to rely on the mixed electronic–ionic processes, involving electrochemical metallization and lattice atoms migration in response to the externally modulated electric field strength. The chemical potential generated by the electronic–ionic ITO/ReO3@ReS2/Al resistive memory cell non-equlibrium processes to the occurrence of the nanobattery effect. This finding supports the possibility of a nonvolatile memory cell with a new operation principle based on the potential read function.

Keywords: rhenium trioxide; rhenium disulfide; transition dichalcogenide; memristor; memcapacitor; nanobattery

1. Introduction In 2009, Di Ventra et al. [1] extended the memristor definition, which was first introduced in 1971 by Chua [2], to memcapacitor and meminductor as two newly generalized energy storage components, which are also referred to as memelements. A memcapacitor is a multi-dimensional electronic device with memory and capacitance capability, where its state can be controlled via external stimuli such as electric field (EF) or voltage (V). Memcapacitor application as a resistive random access memory (RRAM) is one of the most promising in the field of next-generation nonvolatile memory (NVM), owing to its superior scalability, low power consumption, and high speed. In an ideal memcapacitor, the capacitance depends only on the history of the charge stored on the plates or the voltage across them. The memcapacitor main characteristic is a hysteretic loop, which may or may not pass through

Nanomaterials 2020, 10, 2103; doi:10.3390/nano10112103 www.mdpi.com/journal/nanomaterials Nanomaterials 2020, 10, 2103 2 of 13 the origin [1] when driven by a periodic input. The complex mechanisms of the ionic transport, among other identified mechanisms, inherited in memristive systems introduce new challenges from the perspectives of modeling, characterization, and system architectures. The unique memory characteristics and nanoscale dimensions of the memcapacitor have a wide application prospect in the next-generation nonvolatile memory computing, with multibit high-density data storage without a power source [3]. Moreover, it represents an important research direction aiming to solve the von Neumann bottleneck in modern computers [4]. Memristive devices provide a promising platform to perform biologically inspired computing [5,6], as well as the potential to realize functionalities of neurons and synapses in order to build artificial neural networks [7–9]. The memcapacitor application interest, apart from aforementioned, is extended to other areas such as low-pass filters and circuit design [10]. Furthermore, it is suitable for building chaotic circuits for secure communication systems [11,12], which have attracted vast worldwide interest from both academia and industry. In addition, the memcapacitor stimulated considerable interest in the field of nanoionics, since it showed promise for integration with the complementary metal oxide semiconductor (CMOS) technology [13,14]. Apart from information storage and processing, memcapacitors, unlike memristors, can store energy in the form of the EF, therefore opening new possibilities in the concept of the energy storage and redistribution [15,16]. Despite the prospects of innovative application of memcapacitor, it is not yet available commercially in any form. Nonetheless, progress in the field has been mostly achieved in the theoretical work concerning simulation of the memcapacitive circuit models [17,18], with only several reports exhibiting memcapacitive characteristics realized in practical solid-state devices [19–23]. Among them, Chu et al. have reported memcapacitive and memristive properties of a device composed of an Ni-DNA nanowire sandwiched between Au electrodes, where the dual properties were attributed to the change of the Ni [19]. Nonvolatile memcapacitance, induced by oxygen ions migration between indium tin oxide (ITO) and the HfOx layer was reported by Park et al. in the developed ITO/HfOx/Si metal-oxide semiconductor device [20], and by Yang et al. in a reactive electrode (Mo, Al)/HfOx/n-Si structure [21]. In another work of Noh et al. memristive + and memcapacitive attributes were demonstrated in the Ti/Pt-Fe2O3 core–shell nanoparticles/p -Si system as an effect of the nanoparticles potential variation under charging processes [22]. In addition, the memcapacitive property was also realized in a work of Zhao et al. in a new type of photoelectric memcapacitor with a planar Au/La1.875Sr0.125NiO4/Au structure [23]. With reference to the resistance switching phenomenon, it is worth noting that apart from oxygen ions migration, for many systems, it may also involve more complex redox reactions or thermally activated processes [24]. Apart from the aforementioned examples, the limited (at present) number of the practically realized resistance switching memory cells with capacitive properties is caused by the lack of the functional materials, whose capacitance or resistance can be controlled by external variables [25–27]. Despite the promising feature of information storage with very low energy loss, the memcapacitive systems still remain the least studied in the class of memelements. In this study, the material presented is composed of rhenium disulfide (ReS2) rich in the surface overlayer formed of rhenium (VI) oxide (ReO3)[28]. The ReS2 belongs to the class of 2D transition metal dichalcogenides; meanwhile, ReO3 is a member of the broad family of transition metal . Recently, they have attracted significant research interest due to their rich physics, which is desirable from the perspective of their successful integration into a number of electronic and optoelectronic devices, such as thin film transistors, solar and fuel cells, photodetectors, and digital logic devices [29–35]. Herein, the memcapacitive ITO/ReO3@ReS2/Al device characteristics are presented in detail, as follows. The determination of the charge carrier conduction model, quantitative analysis of its electrical properties, as well as an explanation of the resistance switching mechanism, in addition to the frequency modulated capacitance and conductance profile, were successfully carried out based on the experimental data fitting. It was found that resistance switching relies on the electrochemical metallization (ECM) effect and valence change of the sublattice rhenium (Re) atoms, which are 3+ 2 induced by migration of the aluminum cations (Al ) and the lattice oxygen (O −) atoms, respectively. Nanomaterials 2020, 10, 2103 3 of 13

Most importantly, the ITO/ReO3@ReS2/Al electronic–ionic resistive memory cell was controlled by non-equilibrium states (chemical potential generation) evidenced in the occurrence of the nanobattery effect. The establishment of the electromotive force (emf ) is clearly a function of the chemical and the charge transport processes taking place in the electrochemically active system.

2. Materials and Methods

Fabrication of the ITO/ReO3@ReS2/Al memcapacitor: First, the substrates of indium tin oxide (ITO, 185 5 nm) on a glass (Saide Glass Co. Ltd., Dongguan, China) were patterned using the ≈ ± photolithography method. Prior to the device fabrication, the ITO/glass substrates were cleaned three times in acetone, ethyl, and isopropyl alcohol by treatment in an ultrasonication bath for 10 min and dried in a flow of high-purity N2 gas. The water:ethyl alcohol suspension of the ReO3@ReS2 (v/v of 1:0.1, 1 concentration 1.73 mg mL− ) was prepared by ultrasonication (1 h) and used to form the active layer by the spin-coating process (0.1 mL, 500 rpm, 30 s) on the cleaned ITO/glass substrate. The active material layer thickness was 85.55 3.2 nm. In the final step, the aluminum (Al, Alfa Aesar, Lancashire, UK) ± top electrode ( 120 4 nm) was deposited using shadow mask by thermal evaporation in a physical ≈ ± vapor deposition (PVD) chamber (Chi-Vac Research and Development Co. Ltd., Dalian, China), under 4 vacuum of 10− Pa. The acetone, ethyl, and isopropyl alcohol were provided by Keshi Chemical Co. Ltd., Chengdu, China. Characterization of the ITO/ReO3@ReS2/Al memcapacitor: Electrical characterization measurements of the fabricated devices were performed with a Keysight B2901A Precision Source Measure Unit (Keysight, Shanghai, China) in air and at ambient temperature (22 1 C, relative ± ◦ humidity 30%), by applying voltage (V) and measuring current (I) between the top (Al) and the ≈ 2 bottom (ITO) electrode. A compliance current (Icc) of 1 10 A was applied for SET (device switch into × − ON state) and RESET (device switch into OFF state) processes, respectively, unless otherwise stated. The memcapacitive characteristics of the ITO/ReO3@ReS2/Al device were examined by recording capacitance–voltage (C-V) measurements, by sweeping DC voltage with a superimposed AC signal of 30 mA, on the TTPX Probe Station (Lake Shore Cryotronics, Inc., Westerville, OH, USA) controlled with 4200-SCS Semiconductor Characterization System (Keithley Instruments, Inc., Cleveland, OH, USA). For all the measurements, the voltage was applied to the top electrode, while the bottom electrode was grounded. Impedance spectroscopy measurements were performed on the Metrohm Autolab with PGSTAT30 instrument (Eco Chemie B.V., Utrecht, The Netherlands).

3. Results and Discussion

3.1. Physicochemical Characterization of the ReO3@ReS2 6+ The ReO3@ReS2 solid electrolyte is composed of the rhenium (VI) (Re (4f7/2, 44.98 eV), oxygen 2 4+ 2 (O − 1s, 530.59 eV), rhenium (IV) (Re 4f7/2, 41.85 eV), as well as (S − 2p3/2, 162.96 eV), which is confirmed from the high-resolution X-ray photoelectron spectroscopy (XPS, AXIS UltraDLD, Kratos Analytical Ltd., Manchester, UK) analysis results. Based on the data fitting, the ratio of the rhenium (VI) oxide (ReO3) to the rhenium disulfide (ReS2) was estimated to be 0.8:1.0 (for details, see Figure S1 in the electronic supplementary material (ESM)). It should be noted that in the ReO3 6+ 7+ material, the Re atoms are present in the mixed Re and Re oxidation state; therefore, the ReO3 should not be considered as a highly conducting metal oxide, as it is expected from its doped crystal form [36]. Moreover, the transition metal oxides frequently exhibit lattice disorder associated with structural nonstoichiometry [37]. The results from the energy-dispersive X-ray (EDX, FEI Quanta FEG 250, FEI, Hillsboro, OR, USA) spectroscopy confirmed the presence of the Re, S, and O elements in the ReO3@ReS2 composite material (for details, see Figure S2 in the ESM). Nanomaterials 2020, 10, 2103 4 of 13

3.2. Electronic Characterization of the ITO/ReO3@ReS2/Al Diode

The initial current–voltage (I-V) sweep of the ITO/ReO3@ReS2/Al device presented in Figure1 exhibitsNanomaterials the formation 2020, 10, x FOR of conductive PEER REVIEW nanofilaments (CNFs) taking place at low bias voltage, VSET, form4 of 13 2.57 0.06 V. The abrupt bipolar resistance switching between the high (HRS) and low (LRS) resistance ± state,resistance in the forwardstate, in (FB)the forward as well as(FB) in theas well reverse as in bias the (RB) reverse sweep bias direction, (RB) sweep is typically direction, observed is typically in theobserved case of inducedin the case memristance of induced or memristance redox reaction or driven redox devicesreaction [38 driven]. The devices significant [38] variation. The significant in the resistancevariation and in the the operating resistance voltages and the with operating a consecutive voltages voltage with sweep a consecutive can be ascribed voltage to sweep the random can be 3 2 redistributionascribed to theand random incorporation redistribution of the mobile and ionsincorporation in the ReO of3@ReS the 2 mobilecomposite ions material in the ReO network@ReS undercomposite an external material electric network field under (EF) (see an Figureexternal1). electric field (EF) (see Figure 1).

-4 10 2 4

-5 10

-6 10 1

-7 10

3 -8 HRS 10 LRS

-3 -2 -1 0 1 2 3 V (V)

FigureFigure 1. 1.Initial Initial sweep sweep of of the the ITO ITO/ReO/ReO3@ReS3@ReS2/Al2/Al memcapacitor, memcapacitor, measured measured at at a frequencya frequency (f ): (f): 0.6 0.6 Hz. Hz. 2 SweepingSweeping voltage voltage range: range: ±33 V V,, and compliance current current (I (ccI):cc ):1 × 1 10−210 A. ArrowsA. Arrows indicate indicate voltage voltage sweep ± × − sweepsequence. sequence.

TheThe space space charge charge limited limited current current (SCLC) (SCLC) conduction conduction model model of of the the charge charge carrier, carrier, shown shown in in FigureFigure2a,b 2a (initial,b (initial scan) scan) as identified as identified from from the appearance the appearance of the ohmicof the regionsohmic (regionsIαV, FB (I sweep),αV, FB pointssweep), onpoints the formation on the formation of the CNFs of the within CNFs the within ReO3 the@ReS ReO2 electrolyte3@ReS2 electrolyte layer [39 layer]. [39]. TheThe ohmic ohmic conduction conduction (LRS (LRS under under FB, FB,V V< 0.62 voltages V in the ( HRSV > 0.62/LRS V under in the FBHRS/LRS and RB under sweep) FB is and attributed RB sweep) to the is dominatedattributed to bulkthe ReOdominated3@ReS2 bulkconductivity. ReO3@ReS2 The conductivity. hysteresis The curve hysteresis of the fifth curve consecutive of the fifthI–V consecutivesweep is I– equallyV sweep representedis equally represented by the SCLC by conduction the SCLC conduction model (see model Figure 2(seec,d). Figure Under 2c consecutive,d). Under consecutive voltage sweeps, voltage fromsweeps, the initial from to the the fifth initial scan, to the the memcapacitor fifth scan, the electrical memcapacitor properties electrical underwent properties observable underwent change. Theobservable dynamic processeschange. The taking dynamic place processes in the device taking resulted place in in the th decreasee device resulted of the charge in the carrier decrease mobility of the 11 11 2 1 1 (µcharge, from 6.95carrier 10mobility− to 5.11(μ, from10 −6.95cm × 10V−1−1 tos− 5.11) carrier × 10−1 concentration1 cm2V−1s−1) carrier in thermal concentration equilibrium in thermal (no, ×20 19× 3 20 20 3 fromequilibrium 1.24 10 (no,to from 4.28 1.2410 × 10cm20 to− ),4.28 number × 1019 ofcm trap−3), number states (N oft, trap from states 7.93 (N10t, fromto 3.79 7.93 × 101020 tocm 3.79− ), × × × 19 19 3 × × and1020 free-carrier cm−3), and concentration free-carrier (concentrationn, from 9.11 (10n, fromto 4.08 9.11 10× 10cm19 to− ), 4.08 as well× 10 as19 thecm−3 e),ff ective as well density as the × 14 × 14 3 of states in the conduction band (Nc, from 4.84 10 to 1.78 10 cm ). Corresponding14 results14 from−3 effective density of states in the conduction× band (N×c, from 4.84− × 10 to 1.78 × 10 cm ). theCorresponding evaluation of theresults electric from properties the evaluation of the of device the electric attained properties from the of experimental the device attained data fitting from are the listedexperimental in Table S1 data in thefitting ESM. are listed in Table S1 in the ESM.

Nanomaterials 2020, 10, x FOR PEER REVIEW 5 of 13

Nanomaterials 2020, 10, 2103 5 of 13

-2 (a) -3 (b) -4 -4

-6 -5 -2 -3 -4 -4 -6 -5 -2 LRS LRS HRS -3 HRS -4 -4

-6 -5

-1.0 -0.5 0.0 0.5 -1.0 -0.5 0.0 0.5 log V (V) log V (V)

-2 (c) -3 (d) -4 -4 -5 -6 -6 -2 -3 -4 -4 -5 -6 -6 -2 LRS LRS HRS -3 HRS -4 -4 -5 -6 -6 -1.5 -1.0 -0.5 0.0 0.5 -1.5 -1.0 -0.5 0.0 0.5

log V (V) log V (V)

FigureFigure 2. 2.Plot Plot of of log logJ Jvs. vs. loglog VV of the ITO/ReO ITO/ReO3@ReS@ReS22/Al/Al for for the the initial initial scan scan under under (a ()a negative) negative and and (b) (bpositive) positive bias; bias; and and for for the the fifth fifth consecutive consecutive scan underunder ((cc)) negativenegative andand ( d(d)) positive positive bias. bias. Sweeping Sweeping 2 voltagevoltage range: range: 3±3 V, V,Icc Icc:: 1 1 × 1010 −2 AA,, and and frequencyf : 0.6 Hz. (f): 0.6 Hz. ± × − 11 2 1 1 TheThe low low carrier carrier mobility mobility (µ ()μ of) of the the order order of of10 − 10−1cm1 cmV2−V−1s−s−1 isis an an indication indication of of the the mixed mixed electronic–ionicelectronic–ionic conduction conduction of the of thecharge charge carrier carrier in the in ReO the3 @ReS ReO32@ReSelectrolyte2 electrolyte (for details (for of details the charge of the carriercharge conduction, carrier conduction please see, Figureplease2 seeand Figure Table S1). 2 and The bipolarTable S1). switching The bipolar mode is switching attributed mode to the is 3+ 2 electrochemicalattributed to the redox electrochemical reactions and redox the highreactions solubility and the of thehigh Al solubilityand O −ofions the inAl3+ the and ReO O23−@ReS ions 2in materialthe ReO network.3@ReS2 material The inversely network. proportional The inversely surface area proportional dependence surface on the area memcapacitor dependence resistance on the (Rmemcapacitor, see Figure S3 resistance in the ESM) (R, implies see Figure that the S3 CNF-type in the ESM) resistance implies switching that the is CNF not-type the dominant resistance effswitchingect in the is ITO not/ReO the3 @ReSdominant2/Al, effect but it in rather the ITO/ReO occurs over3@ReS the2/Al whole, but area it rather of the occurs cell interface over the [27 w,41hole]. 2- Underarea of FB the (positively cell interface biased [27,41 top].electrode Under FB (TE)), (positively the EF-induced biased top O electrodemigration (TE)), (drift the andEF-induced diffusion) O2- inmigration the lattice (drift sites leadsand diffusion) to the consecutive in the lattice growth sites of leads the oxygen to the vacancy consecutive (Vo ) growth channels of expanding the oxygen 6+ towardvacancy the (V anode.o) channels Simultaneously, expanding the toward oxidation the anode process. Simultaneously, of the transition the metal oxidation cations process (Re ) in of the the Altransition electrode metalvicinity cations prevents (Re the6+) V inO channels the Al from electrode short-circuiting vicinity prevents the electrodes the V [O42 channels]. In the HRS, from 3+ theshort oxygen-deficient-circuiting the channelselectrodes serve [42]. asIn thethe activeHRS, the sites oxygen for the-deficient Al ions channels precipitation serve andas the subsequent active sites growthfor the of Al the3+ highlyions precipitation conductive metal and subsequent NFs, propagating growth in of the the direction highly from conductive the anodic metal to theNFs, cathodicpropagating interface in the (see direction Equation from (S2) in the the anodic ESM and to the Scheme cathodic1a) [ 43 interface]. The ITO (see/ReO Equation3@ReS2 /S2Al in resistance the ESM 3+ switchingand Scheme is driven 1a) [43] by. the The ECM ITO/ReO process3@ReS [24].2/Al The resistance mobile Al switchi, driftingng is in driven the ion by conducting the ECM ReOproces3@ReSs [24]2 . layer,The mobile is formed Al3+ under, drifting anodic in the dissolution ion conducting of the ReO active3@ReS Al2 TE,layer, which is formed is experimentally under anodic verified dissolution by

Nanomaterials 2020, 10, 2103 6 of 13 Nanomaterials 2020, 10, x FOR PEER REVIEW 6 of 13

of the active Al TE, which is experimentally verified by the temperature-depressed conductivity of the temperature-depressed conductivity of the ITO/ReO3@ReS2/Al (as evidenced by the derived Al the ITO/ReO3@ReS2/Al (as evidenced by the derived3 Al1 resistance temperature coefficient (α) of 3.9 × resistance temperature coefficient (α) of 3.9 10− K− being in agreement with the theoretical value, −3 −1 × see10 FigureK being S4 in agreement the ESM) [ 44with]. Thethe theoretical less mobile v Realue,6+ cationssee Figure compensate S4 in the ESM) the oxygen [44]. The deficiency less mobile by 6+ trappingRe cations electrons compensate released the from oxygen the cathodedeficiency (see by Equation trapping (S3) electrons in the ESM), released where from the the electrons cathode fill (see up theEquation empty (S3) states in the of theESM), Re 5whered band. the Furthermore, electrons fill theup the drift empty of the states O2- toward of the Re the 5d cathode band. Furthermore, gives rise to the drift of the O2- toward the cathode gives rise to the evolution of molecular oxygen (O2(g), see the evolution of molecular oxygen (O2(g), see Equation (S4) in the ESM), which is revealed through the Equation (S4) in the ESM), which is revealed through the formation2 of surface bubbles (see Figure S5 in formation of surface bubbles (see Figure S5 in the ESM). The O − depletion (ReO3-x, where x is the the ESM). The O22− depletion (ReO3-x, where x is the number of the O2−) generates a high number of VO, number of the O −) generates a high number of VO, which are considered to act as the effective donors inwhich the n-type are considered oxide semiconductors to act as the [ effective45]. This donors process in is portrayedthe n-type through oxide semiconductors an initial decrease [45] of. This the electrolyteprocess is portrayed/BE contact through resistance an initial (gradual decrease increase of the of electrolyte/BE the HRS current, contact see Figureresistance S6 in(gradual the ESM) increase with of the HRS current, see Figure S6 in the ESM) with consecutive7 1voltage sweeps (maximum EF of 5.8 × 107 consecutive voltage sweeps (maximum EF of 5.8 10 Vm− ), due to the depletion of the available Vm−1), due to2 the depletion of the available interfacial× O2− at the electron injecting electrode interface. interfacial O − at the electron injecting electrode interface. Finally, the reset of the cell into the OFF stateFinally, is completedthe reset of underthe cell RB into polarity the OFF sweep, state whichis completed induces under the electrochemical RB polarity sweep, dissolution which induce of thes topthe partselectrochemical of the Al CNFs dissolution (see Scheme of the top1b). parts of the Al CNFs (see Scheme 1b).

Scheme 1. Schematic representation of the ion redistribution within the ITO/ReO3@ReS2/Al under Scheme 1. Schematic representation of the ion redistribution within the ITO/ReO3@ReS2/Al under (a) (a) forward (V > 0) and (b) reverse (V < 0) bias voltage. forward (V > 0) and (b) reverse (V < 0) bias voltage. The experimentally validated EF-driven ion transport was found to be in agreement with the thermallyThe experimentally induced ion hopping validated in theEF- FBdriven regime ion (Motttransport and was Gurney found model, to be see in Equationagreement (S5) with in the ESM)thermally [46]. induced The activation ion hopping energy in (W the0, at FB 295 regime K) for (Mott the hopping and Gurney in the model, absence see of Equation EF was estimated (S5) in the to a ESM) [46]. The activation energy ( , at 295 K) for the hopping in the absence of EF was estimated to be 0.34 (y(1)) and 0.27 eV (y(2)) in HRS and LRS, respectively (ln JαE regime, see Figure S7 in the ESM), whichbe 0.34 is (y comparable(1)) and 0.27 toeV the(y(2) average) in HRSW and0 for LRS, the respectively electron hopping. (ln JαE Furthermore, regime, see Figure the electrons S7 in the emitted ESM), a fromwhich the is trap comparable states in the to thermally the average activated Poole–Frenkelfor the electron (P-F) hopping. process were Furthermore, also found the to contribute electrons emitted from the trap states in the thermally1 1activated/2 Poole–Frenkel (P-F) process were also found to the charge carrier conduction (ln JE− vs. E plot, see Figure S8 in the ESM). The barrier height of to contribute to the charge carrier conduction (ln JE−1 vs. E1/2 plot, see Figure S8 in the ESM). The the traps (qϕT), derived from the fitting of the P-F plot was equal to 0.69 0.01 eV. The temperature barrier height of the traps (qφT), derived from the fitting of the P-F plot was± equal to 0.69 ± 0.011 eV.1 -dependent activation energy (Ea) for the carrier transport of the ITO/ReO3@ReS2/Al (ln J vs. kb− T− , The temperature-dependent activation energy (Ea) for the carrier transport of the ITO/ReO3@ReS2/Al see Figure S9 in the ESM) exhibited two regions: the charge-hopping conduction in the ReO3@ReS2 −1 −1 framework(ln J vs. kb T at, lower see Figure temperatures S9 in the (T ESM)< 34.2 exhibited K), comprising two regions: the di fftheusion charge of carrier-hopping between conduction existing in the ReO3@ReS2 framework at lower temperatures (T < 34.2 K), comprising the diffusion of carrier vacancy and/or interstitial sites (Ea = 0.51 eV, extrinsic region), and, as expected for ionic conductors, between existing vacancy and/or interstitial− sites (Ea = −0.51 eV, extrinsic region), and, as expected carrier diffusion and simultaneous processess of formation of new vacancy sites (Ea = 0.13 eV, intrinsic region)for ionic taking conductors, place at carrier elevated diffusion temperatures and simultaneous (T > 34.2 K). Theprocessess dynamic of controlformation over of the new mobile vacancy ions redistributionsites (Ea = 0.13 within eV, intrinsic the electrolyte region) network,taking place achieved at elevated through temperatures adequate selection (T > 34.2 of theK). EFThe strength dynamic as wellcontrol as the over number the mobile of the voltage ions redistribution sweep cycles (see within Figure the1, and electrolyte Figure S6 network, in the ESM), achieved can e ff throughectively impactadequate the selection memcapacitor of the EF electronic strength state as well [27]. as The the memristor number of resistance, the voltage after sweep five cycles consecutive (see Figure voltage 1, sweepsand S6 in (∆ V:the ESM),3 V), was can demonstratedeffectively impact to be the in agreementmemcapacitor with electronic the characteristic state [27] shrinking. The memristor of the ± deviceresistance, hysteresis after five loop consecutive with rise of voltage the scan sweeps frequency (ΔV: ( f±3), asV), a was consequence demonstrated of a slow to be response in agreement of the chargewith the carrier characteristic to the high-frequency shrinking of EFthe oscillationdevice hysteresis (FB, see loop Figure with3a) [rise2]. of However, the scan the frequency principle (f of), as the a consequence of a slow response of the charge carrier to the high-frequency EF oscillation (FB, see zero-crossing hysteresis of the ITO/ReO3@ReS2/Al is violated due to the establishment of the chemical Figure 3a) [2]. However, the principle of the zero-crossing hysteresis of the ITO/ReO3@ReS2/Al is violated due to the establishment of the chemical potential gradient (nanobattery) within the device

Nanomaterials 2020, 10, 2103 7 of 13 Nanomaterials 2020, 10, x FOR PEER REVIEW 7 of 13 Nanomaterials 2020, 10, x FOR PEER REVIEW 7 of 13 potentialactive layer gradient [47]. At (nanobattery) frequencies withinof 0.6 and the 6Hz, device the active generated layer electromotive [47]. At frequencies force (emf of 0.6, see and Figure 6Hz, 4) theactive generated layer [ electromotive47]. At frequencies force (ofemf 0.6, see and Figure 6Hz,4 )the resulted generated in the electromotive cell voltage ( forceV )(emf reaching, see Figure values 4) resulted in the cell voltage (Vemf) reaching values of 0.88 and 1.55 V, respectively.emf ofresulted 0.88 and in 1.55 theV, cell respectively. voltage (Vemf) reaching values of 0.88 and 1.55 V, respectively.

(a) (b) -6 (a) -6 (b) 4 x10-6 4 x10-6 4 x10 4 x10 0.06 Hz 60 Hz 60 Hz 0.060.6 Hz Hz 600 Hz 600 Hz 0.66 Hz Hz 6000 Hz 6000 Hz 2 6 Hz 2 ) 2 ) 2 A A ) ) ( (

A A I I ( (

I I

0 0 0 0

-2 -2 -2 -2 -3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3 V (V) V (V) V (V) V (V) Figure 3. Plot of I vs. V as a function of scan frequency in the range of (a) 0.06–6 Hz and (b) 60–6000 FigureFigure 3. 3.Plot Plot of ofI vs. I vs.V asV as a functiona function of scanof scan frequency frequency in the in rangethe range of (a of) 0.06–6 (a) 0.06 Hz–6 and Hz (band) 60–6000 (b) 60– Hz.6000 Hz. Sweeping voltage range: ±3 V, and Icc: 1 × 10−2 A. 2 −2 SweepingHz. Sweeping voltage voltage range: range:3 V, and±3 V,Icc and: 1 Icc10: 1 × 10A. A. ± × −

FigureFigure 4. 4. TheThe II vs. VV dependencedependence on on the the scan scan frequency frequency at at 0.6 0.6 and and 6 Hz 6 Hz of ofthe the ITO/ReO ITO/ReO3@ReS3@ReS2/Al.2/Al. The Figure 4. The I vs. V dependence on the scan frequency at 0.6 and 6 Hz of the ITO/ReO3@ReS2/Al. The −2 2 Thesweeping sweeping voltage voltage range: range: ±3 V3,V, and and IccI:cc 1: × 1 10 10 A.− A. sweeping voltage range: ±3± V, and Icc: 1 × 10×−2 A. TheThe appearance appearance ofof the emfemf, ,which which has has been been reported reported in seve inral several devices devices [48,49] [,48 justifies,49], justifies the necessity the The appearance of the emf, which has been reported in several devices [48,49], justifies the necessity necessityof the expansion of the expansion of the memristive of the memristivetheory to account theory for tothe account extended for memcapacitive the extended and memcapacitive meminductive of the expansion of the memristive theory to account for the extended memcapacitive and meminductive andelements meminductive observed in elements practical observed memristive in devices practical [1,48,50] memristive. The presence devices of [1 ,the48, 50active]. The element presence is clearly of elements observed in practical memristive devices [1,48,50]. The presence of the active element is clearly thea active function element of the is clearlydevice a non function-equilibrium of the state,device which non-equilibrium is illustrated state, by whichinhomogeneous is illustrated charge by a function of the device non-equilibrium state, which is illustrated by inhomogeneous charge inhomogeneousredistribution within charge the redistribution active layer, within and its the dynamic active layer, response and itsto thedynamic varied response frequency to themodulation. varied redistribution within the active layer, and its dynamic response to the varied frequency modulation. frequencyThese processes modulation. generating These the processes emf in the generating ITO/ReO3@ReS the emf2/Alin memcapacitor the ITO/ReO can3@ReS be associated2/Al memcapacitor either with These processes generating the emf in the ITO/ReO3@ReS2/Al memcapacitor can be associated either with canthe be concentration associated either gradient with the betw concentrationeen Al3+ and gradient the majority between charge Al3+ carrierand the and/or majority the charge Al3+ potential carrier the concentration gradient between Al3+ and the majority charge carrier and/or the Al3+ potential anddifference/or the Al between3+ potential the TE diff anderence CNF between [50]. However, the TE and at frequencies CNF [50]. However, above 6 Hz at, frequenciesthe current abovetrend follows 6 Hz, difference between the TE and CNF [50]. However, at frequencies above 6 Hz, the current trend follows thethe current open circuit trend voltage follows (OCV) the open hysteresis, circuit which voltage is typic (OCV)al for hysteresis, capacitors which (see Figure is typical 3b). for capacitors the open circuit voltage (OCV) hysteresis, which is typical for capacitors (see Figure 3b). (see FigureMoreover,3b). the dynamic extraction of the O2- from the ReO3@ReS2 toward the Al electrode under Moreover, the dynamic extraction of the O2- from the ReO3@ReS2 toward the Al electrode under FB sweep [51,52] decreased the electrolyte relative permittivity (εr, initial ≈220 to εr, final ≈36) and in effect FB sweep [51,52] decreased the electrolyte relative permittivity (εr, initial ≈220 to εr, final ≈36) and in effect caused a variance in the device capacitance (C, decrease of 83.5% where Cinitial = 45.5 nF), which is caused a variance in the device capacitance (C, decrease of 83.5% where Cinitial = 45.5 nF), which is restrained by the active material dielectric constant. The capacitance modulation is strictly a function restrained by the active material dielectric constant. The capacitance modulation is strictly a function

Nanomaterials 2020, 10, 2103 8 of 13

2- Moreover, the dynamic extraction of the O from the ReO3@ReS2 toward the Al electrode under FB sweep [51,52] decreased the electrolyte relative permittivity (ε 220 to ε 36) and in r, initial ≈ r, final ≈ effNanomaterialsect caused 20 a variance20, 10, x FOR in PEER the device REVIEW capacitance (C, decrease of 83.5% where Cinitial = 45.5 nF), which8 of is 13 restrained by the active material dielectric constant. The capacitance modulation is strictly a function d ofof the the mean mean Schottky Schottky barrier barrier width width (W (Wd)) as as well well as as the the chemical chemical potential potential gradient; gradient; however, however it, it is is independentindependent from from the the formation formation and and rupture rupture of of the the CNFs CNFs [53 [53]]. The. The applied applied BV BV level level not not only only modulates modulates thethe structural structural properties properties of of the the material material but but also also the the corresponding corresponding active active layer layer/electrode/electrode interfaces, interfaces, whichwhich is is portrayed portrayed through through the the dynamic dynamic variation variation of of the the current current with with consecutive consecutive number number of of sweeping sweeping cyclescycles (see (see Figures Figure S6S6 and and S10 S10 in in the the ESM). ESM). The The Schottky Schottky barrier barrier (q ϕ(qbφ)b profile) profile of of the the electrolyte electrolyte/electrode/electrode junctionjunction can can be be electrically electrically modulated modulated either either by by movement movement of of the the vacancies vacancies away away/toward/toward the the cell cell interfaceinterface by by applying applying reverselyreversely polarizedpolarized bias bias stress stress [54] [54,], or or by by carrier carrier trapping/detrapping trapping/detrapping at defect at defect sites sitesnear near the junction the junction interface interface [23]. [The23]. process The process of filling of fillingoff the o interfaceff the interface trap states trap with states injected with injected electrons F electrons(HR and (HR LR and state LR in state FB, in see FB, Figure see Figure S10 S10 in inthe the ESM) ESM) increase increasess the the Fermi level level (E (EF)) and and e ff effectivelyectively b diminishesdiminishes the theq ϕqφb for for electron electron injection. injection. The TheC-V C-Vcurves curves at at the the electrode electrode interfaces interfaces show show hysteretic hysteretic d behaviorbehavior (see (see Figure Figure S11 S11 in in the the ESM), ESM), where where the theW Wd is is narrower narrower in in the the LRS, LRS, facilitating facilitating the the electrons electrons tunnelingtunneling through through the the thin thin Schottky Schottky barrier barrier (this (this is is in in agreement agreement with with the the positive positive shift shift of of the the flat-band flat-band FB voltage,voltage,∆ ΔVVFB,, asas discusseddiscussed below)below) [ 55[55]].. ThisThis mechanismmechanism hashas beenbeen proven proven to to be be also also valid valid for for devices devices based on materials [56]. Significantly increased ITO/ReO3@ReS2/Al rectification (FB range) based on perovskite materials [56]. Significantly increased ITO/ReO3@ReS2/Al rectification (FB range) waswas observed observed under under enhanced enhancedEF EFstrength strength upup toto 8.28.2 × 101077 VmVm−1 (ΔV:1 (∆V: ±7 V,7 V,see see Figure Figure S12 S12 in the in the ESM). ESM). The × − ± elevated EF induced motion of the O2- (better2- described by the migration of VO), which in turn triggers The elevated EF induced motion of the O (better described by the migration of VO), which in turn triggersthe electrolyte the electrolyte reconfiguration reconfiguration and extension and extension of the conducting of the conducting channels, channels, finally yields finally device yields electrode device electrodebridging bridging (R = 6.35 ( R× 10= 6.35−3 Ω in 10FB,3 seeΩ inFigure FB, see S12 Figure in the S12ESM). inthe ESM). × − With reference to the aforementioned mechanisms, the ITO/ReO3@ReS2/Al memristor C-V With reference to the aforementioned mechanisms, the ITO/ReO3@ReS2/Al memristor C-V dependence on the scan frequency (fAC) relies strongly on the device interface profile and bulk dependence on the scan frequency (f AC) relies strongly on the device interface profile and bulk distributiondistribution of ofthe the charge charge carrier carrier and trap and states trap states (see Figure (see 5 Figure)[ 57]. 5)Under [57]. exertion Under ofexertion the large of negative the large BVnegative ( 20V

-8 -8 1.5 x10 1.5 x10 (a) (b)

1.0

1.0 0.5 ) ) F F ( (

C C

0.5 0.0 1kHz 10kHz 5MHz 5kHz 1MHz 1kHz 1MHz -0.5 5kHz 5MHz 0.0 10kHz

-4 -2 0 2 4 -4 -2 0 2 4 V (V) V (V)

FigureFigure 5. 5.The TheC Cvs. vs.V Vcurves curves ofof thethe ITOITO/ReO/ReO33@ReS22/Al/Al in in the the frequency frequency range range of of 11 kHz kHz–5–5 MHz MHz for for (a) (astandard) standard and and (b ()b short) short-circuited-circuited devices. devices. In In(b) (, bthe), the C-VC-V curvecurve at aat sca an scan frequency frequency of 1 ofkHz 1 kHzwas wasfitted fittedwith with a polynomial a polynomial function function (original (original data datacan be can found be found in Figure in Figure S13). S13).Higher Higher resolution resolution C-V curvesC-V curvesacquired acquired at frequencies at frequencies from from1 kHz 1 kHzto 5 MHz to 5 MHz (a), and (a), andat 1 atand 1 and5 MHz 5 MHz (b) (areb) areshown shown in Figure in Figure S14. 2 S14.Sweeping Sweeping voltage voltage range: range: ±5 V,5 and V, and Icc: 1Icc ×: 10 1 −2 A.10 A. ± × −

A large capacitance is formed due to the accumulation of the incremental charge at the inversion layer of the semiconductor (C = 1.20 ± 0.02 × 10−8 F), which remains in the equilibrium state irrespective of the level of the applied voltage. A subsequent increase of the potential (V > VT) induces the depletion of the minority carrier from the inversion region, leading to the generation of the depletion layer (Ei crosses below EF, see Scheme S1 and Equation (S6) in the ESM). As expected, the highest capacitance dependence on the scan frequency occurs in the strong inversion regime. The inversion layer is not formed at high frequencies (f ≥ 5 kHz, non-equilibrium state), where the

Nanomaterials 2020, 10, 2103 9 of 13

A large capacitance is formed due to the accumulation of the incremental charge at the inversion layer of the semiconductor (C = 1.20 0.02 10 8 F), which remains in the equilibrium state irrespective ± × − of the level of the applied voltage. A subsequent increase of the potential (V > VT) induces the depletion of the minority carrier from the inversion region, leading to the generation of the depletion layer (Ei crosses below EF, see Scheme S1 and Equation (S6) in the ESM). As expected, the highest capacitance dependenceNanomaterials 20 on20 the, 10, scanx FORfrequency PEER REVIEW occurs in the strong inversion regime. The inversion layer is9 not of 13 formed at high frequencies (f 5 kHz, non-equilibrium state), where the minority carriers are not ≥ ableminority to respond carriers to theare appliednot able externalto respond high-frequency to the applied AC external signal, high and- thefrequency growing AC charge signal, number and the isgrowing kept at thecharge edge number of the depletionis kept at the region. edge Underof the depletion these circumstances, region. Under the these depletion circumstances, layer width the (Wdepletiondl) keeps layer increasing width beyond (Wdl) ke itseps maximum increasing thermal beyond equilibrium its maximum value, thermal leading equilibrium to a capacitance value, thatleading further to decreases a capacitance with voltagethat further (see Equation decreases (S7) with in the voltage ESM). (see Importantly, Equation the (S7) strong in the positive ESM). shiftImportantly, of the flat-band the strong voltage positive (∆VFB shift~2.85 of V)the observed flat-band for voltage the C-V (ΔVsweepFB ~2.85 acquired V) observed at frequency for the C of-V 1sweep kHz (see acquired Figure 5ata) frequency is a clear indication of 1 kHz (see of the Figure process 5a) ofis trappinga clear indication of the electrons of the process injected of from trapping the electrodeof the electr [58].ons injected from the electrode [58]. Interestingly,Interestingly, the the appearance appearance of of the the negative negative capacitance, capacitance, in in both both the the FB FB and and RB RB region, region, was was observedobserved after after short-circuiting short-circuitingthe the devicedevice electrodeselectrodes (see(see Figure5 5bb).). The increase of of the the frequency frequency as aswell well as as the the BV BV lead to to a a capacitance drop in thethe FBFB region.region. TheThe riserise ofof the the voltage voltage generated generated an an instantinstant injection injection of both of both majority majority and minority and minority charge carriers charge into carriers the deviceinto theactive device layer (influenced active layer by(influenced the properties by the of theproperties injecting of contact the injecting and the contact solid electrolyte and the solid charging electrolyte processes). charging The followingprocesses). slowThe redistributionfollowing slow of redistribution the minority of carrier the minority (slower carrierVO diff (slowerusion in VO comparison diffusion in to comparison the drift of to the the electrons)drift of the within electrons) the ReO w3ithin@ReS the2 brought ReO3@ReS on the2 brought decrease on of the the decrease current and,of the in ecurrentffect, influenced and, in effect, the changeinfluenced of the the di ffchangeusion capacitanceof the diffusion of the capacitance system, as of well the assystem, caused as thewell appearance as caused the ofthe appearance inductive of efftheect inductive at low frequencies effect at low [59]. frequencies This mechanism [59]. This was mechanism identified to was be inidentified agreement to be with in theagreement drop of with the 1 −1 ITOthe/ReO drop3@ReS of the2/Al ITO/ReO conductance3@ReS (normalized2/Al conductance with angular(normalized frequency, with Gangularω− ) with frequency, a subsequent Gω ) rise with of a thesubsequent BV (see Figure rise of6). the BV (see Figure 6).

-7 2.5 x10 (a) (b) -6 10 x10 2.0

1kHz 8 1kHz 5kHz 5kHz 1.5 ) ) 10kHz 10kHz F F ( (

6 1 1 - - ω ω G G 1.0 4

0.5 2

0.0 0

-4 -2 0 2 4 -4 -2 0 2 4 V (V) V (V)

−11 FigureFigure 6. 6.The TheG Gωω− vs. VV sweepsweep of of the the ITO/ReO ITO/ReO3@ReS3@ReS2/Al2/Al memcapacitor memcapacitor for for (a) ( a standard) standard and and (b) −2 2 (bshort) short-circuited-circuited device. device. Voltage Voltage sweep sweep direction direction from from −V Vto to+V,+V, IccI: cc1 :× 1 10 10 A. A. − × − TheThe aforementioned aforementioned mechanism mechanism is supported is supported by the by low the series low resistance series resistance of the device of (seethe Figuredevice S15(see inFigure the ESM) S15. in Moreover, the ESM). at Moreover, low voltages, at low the voltagesG of the, the memristor G of the is memristor low, and theis low, transient and the current transient is dominatedcurrent is by dominated the displacement by the component displacement related component to the release related and to escape the release of the electrons and escape from of the the traps,electrons recharging from the processes, traps, recharging and injection processes properties, and of injection the emitter properties barrier of [57 the]. The emitter decrease barrier of the[57]. inversionThe decrease layer of capacitance the inversion (RB layer region) capacitance can be understood (RB region) from can thebe understood perspective from of widening the perspective of the Wofdl with widening rise of of the the scan Wdl frequency,with rise of which the scanis caused frequency, by the which minority is caused charge depletion. by the minority In addition, charge 6 thedepletion. constant In value addition, of the the capacitance constant value at 0 Vof biasthe capacitance (C = 6.02 at0.05 0 V bias10− (CF), = 6.02 irrespective ± 0.05 × of10−6 the F), ± × irrespective of the frequency, is more likely to be associated with the device geometric capacitance than the charging of the interface states [60]. Nonetheless, the heating effects causing a drop of the device resistance or charge carrier generation from the CNF dissolution might be considered as other probable explanations of the phenomenon.

4. Conclusions In conclusion, reliable and controllable memcapacitive behavior was investigated in the two-terminal device structure utilizing ReO3@ReS2 as a solid electrolyte. The electronic–ionic conduction, with the charge carrier mobility of 10−11 cm2V−1s−1, was identified in the memcapacitor

Nanomaterials 2020, 10, 2103 10 of 13 frequency, is more likely to be associated with the device geometric capacitance than the charging of the interface states [60]. Nonetheless, the heating effects causing a drop of the device resistance or charge carrier generation from the CNF dissolution might be considered as other probable explanations of the phenomenon.

4. Conclusions In conclusion, reliable and controllable memcapacitive behavior was investigated in the two-terminal device structure utilizing ReO3@ReS2 as a solid electrolyte. The electronic–ionic conduction, with the 11 2 1 1 charge carrier mobility of 10− cm V− s− , was identified in the memcapacitor from the SCLC model fitting. It was shown that a moderate rise of the BV across the active layer can generate a large enough EF to achieve the control over the device electronic state. The origin of the occurrence of the emf, as high as 1.55 V, was associated with the frequency-modulated device charge carrier dynamic state variation. Furthermore, the experimental work presented here indicates a strong contribution to the concept of practical memcapacitor modeling and most importantly device structure engineering, from the perspective of understanding the composition of the ionic materials utilized for practical devices fabrication. The chemical potential generated by the electronic–ionic ITO/ReO3@ReS2/Al memcapacitor non-equlibrium processes opens up new possibilities toward memory cells whose operation relies on the potential read function.

Supplementary Materials: The following are available online at http://www.mdpi.com/2079-4991/10/11/2103/s1, Figure S1: High resolution XPS spectra, Figure S2: SEM and EDX spectra, Table S1: Electronic properties of the 1 ITO/ReO3@ReS2/Al, Figure S3: Device area dependent I vs. V plot, Figure S4: Normalized resistance (RRo− ) o vs. temperature (T-T ) plot, Figure S5: SEM image of the ITO/ReO3@ReS2/Al device, Figure S6: The I vs. V plot 1 1/2 in the voltage range 5 V, Figure S7: The ln J vs. E plot, Figure S8: The ln JE− vs. E plot, Figure S9: The ln J 1 1 ± 2 1/2 vs. kb− T− Arrhenius plot, Figure S10: The ln IT− vs. V plot, Figure S11: The C vs. V plot, Figure S12: The I vs. V plot in the voltage range 7 V, Scheme S1: Energy band and charge density diagram, Figure S13: The C vs. V plot at scan frequency of 1 kHz,± Figure S14: The C vs. V plot, Figure S15: Nyquist plot of the impedance spectra. Author Contributions: Formal analysis, J.B., W.L., A.J.B. and W.P.G.; Funding acquisition, J.B.; Investigation, J.B., M.L. and Y.H.; Methodology, J.B.; Resources, J.B. and Y.M.; Supervision, W.P.G.; Validation, T.K. and A.J.B.; Writing—original draft, J.B.; Writing—review & editing, J.B. All authors have read and agreed to the published version of the manuscript. Funding: This research was funded by the National Natural Science Foundation of China (Grant No. 61574095) China, and Sichuan University Full-time Post-doctoral Research and Development Fund (Grant No. 2019SCU12070) Chengdu, China. Conflicts of Interest: The authors declare no conflict of interest.

References

1. Di Ventra, M.; Pershin, Y.V.; Chua, L.O. Circuit elements with memory: Memristors, memcapacitors and meminductors. Proc. IEEE 2009, 97, 1717. [CrossRef] 2. Chua, L.O. Meristor—The missing circuit element. IEEE Trans. Circuit Theory 1971, 18, 507. [CrossRef] 3. Zidan, M.A.; Strachan, J.P.; Lu, W.D. The future of electronics based on memristive systems. Nat. Electron. 2018, 1, 22. [CrossRef] 4. Backus, J. Can programming be liberated from the von Neumann style? A functional style and its algebra of programs. Commun. ACM 1978, 21, 21613. [CrossRef] 5. Rahmani, M.K.; Kim, M.H.; Hussain, F.; Abbas, Y.; Ismail, M.; Hong, K.; Mahata, C.; Choi, C.; Park, B.G.; Kim, S. Memristive and synaptic characteristics of nitride-based heterostructures on Si substrate. Nanomaterials 2020, 10, 994. [CrossRef] 6. Van de Burgt, Y.; Lubberman, E.; Fuller, E.J.; Keene, S.T.; Faria, G.C.; Agarwal, S.; Marinella, M.J.; Talin, A.A.; Salleo, A. A non-volatile organic chemical devices as a low-voltage artificial synapse for neuromorphic computing. Nat. Mater. 2017, 16, 414. [CrossRef][PubMed] 7. Wang, Z.; Rao, M.; Han, J.W.; Zhang, J.; Lin, P.; Li, Y.; Li, C.; Song, W.; Asapu, S.; Midya, R.; et al. Capacitive neural network with neuro-transistors. Nat. Commun. 2018, 9, 3208. [CrossRef] Nanomaterials 2020, 10, 2103 11 of 13

8. Ryu, H.; Kim, S. Pseudo-interface switching of a two-terminal TaOx/HfO2 synaptic device for neuromorphic applications. Nanomaterials 2020, 10, 1550. [CrossRef] 9. Wang, Z.; Joshi, S.; Savelev, S.; Song, W.; Midya, R.; Li, Y.; Rao, M.; Asapu, S.; Zhuo, Y.; Jiang, H.; et al. Fully memristive neural networks for pattern classification with unsupervised learning. Nat. Electron. 2018, 1, 137. [CrossRef] 10. Li, Y.; Yang, C.; Yu, Y.; Fernández Díez, F.J. Research on low pass filter based on memristor and memcapacitor. In Proceedings of the 36th Chinese Control Conference, Dalian, China, 26–27 July 2017; Volume 7. 11. Bao, B.C.; Wang, N.; Xu, Q.; Wu, H.G.; Hu, Y.H. A simple third-order memristive band pass filter chaotic circuit. IEEE Trans. Circuits Syst. II Exp. Briefs 2016, 64, 977. [CrossRef] 12. Wang, C.H.; Zhou, L.; Wu, R.P. The design and realization of a hyper-chaotic circuit based on a flux-controlled memristor with linear memductance. J. Circuit Syst. Comp. 2018, 27, 1850038. [CrossRef] 13. Jeong, D.S.; Thomas, R.; Katiyar, R.S.; Scott, J.F.; Kohlstedt, H.; Petraru, A.; Hwang, C.S. Emerging memories: Resistive switching mechanisms and current status. Rep. Prog. Phys. 2012, 75, 076502. [CrossRef][PubMed] 14. Zheng, Y.; Fischer, A.; Sawatzki, M.; Hai, D.; Matthias, D.; Glitzky, L.A.; Reineke, S.; Mannsfeld, S.C.B. Introducing pinMOS memory: A novel, nonvolatile organic memory device. Adv. Funct. Mater. 2020, 30, 1907119. [CrossRef] 15. Cohen, G.Z.; Pershin, Y.V.; Di Ventra, M. Lagrange formalism of memory circuit elements: Classical and quantum formulations. Phys. Rev. B 2012, 85, 165428. [CrossRef] 16. Dubi, Y.; Di Ventra, M. Colloquium: Heat flow and thermoelectricity in atomic and molecular junctions. Rev. Mod. Phys. 2011, 83, 131. [CrossRef]

17. Khan, A.K.; Lee, B.H. Monolayer MoS2 metal insulator transition based memcapacitor modeling with extension to a ternary device. AIP Adv. 2016, 6, 095022. [CrossRef] 18. Zhao, Q.; Wang, C.H.; Zhang, X. A universal emulator for memristor, memcapacitor, and meminductor and its chaotic circuit. Chaos 2019, 29, 013141. [CrossRef] 19. Chu, H.L.; Lai, J.J.; Wu, L.Y.; Chang, S.L.; Liu, C.M.; Jian, W.B.; Chen, Y.C.; Yuan, C.J.; Wu, T.S.; Soo, Y.L.; et al. Exploration and characterization of the memcapacitor and memristor properties of Ni–DNA nanowire devices. NPG Asia Mater. 2017, 9, e430. [CrossRef] 20. Park, D.; Yang, P.; Kim, H.J.; Beom, K.; Lee, H.H.; Kang, C.J.; Yoon, T.S. Analog reversible nonvolatile memcapacitance in metal-oxide-semiconductor memcapacitor with ITO/HfOx/Si structure. Appl. Phys. Lett. 2018, 113, 162102. [CrossRef] 21. Yang, P.; Noh, Y.J.; Baek, Y.J.; Zheng, H.; Kang, C.J.; Lee, H.H.; Yoon, T.S. Memcapacitive characteristics in reactive-metal (Mo, Al) HfOx/n-Si structures through migration of oxygen by applied voltage. Appl. Phys. Lett. 2016, 108, 052108. [CrossRef] 22. Noh, Y.J.; Baek, Y.J.; Hu, Q.; Kang, C.J.; Choi, Y.J.; Lee, H.H.; Yoon, T.S. Analog memristive and memcapacitive

characteristics of Pt-Re2O3 core-shell nanoparticles assembly on p+-Si substrate. IEEE Trans. Nanotechnol. 2015, 14, 798. [CrossRef] 23. Zhao, L.; Fan, Z.; Cheng, S.; Hong, L.; Li, Y.; Tian, G.; Chen, D.; Hou, Z.; Qin, M.; Zeng, M.; et al. An artificial optoelectronic synapse based on a photoelectric memcapacitor. Adv. Electron. Mater. 2020, 6, 1900858. [CrossRef] 24. Waser, R.; Dittmann, R.; Staikov, G.; Szot, K. Redox-based resistive switching memories-nanoionic mechanisms, prospects, and challenges. Adv. Mater. 2009, 21, 2632. [CrossRef] 25. Akinaga, H.; Shima, H. Resistive random access memory (ReRAM) based on metal oxides. Proc. IEEE 2010, 98, 2237. [CrossRef] 26. Sawa, A. Resistive switching in transition metal oxides. Mater. Today 2008, 11, 28. [CrossRef] 27. Waser, R.; Aono, M. Nanoionics-based resistive switching memories. Nat. Mater. 2007, 6, 833. [CrossRef] 28. Borowiec, J.; Liang, W.; Boi, F.S.; He, Y.; Wang, S.L.; Gillin, W.P. promoted sulfidation of

ammonium perrhenate: Presence of nanobattery in the ReS2 composite material based memcapacitor. Chem. Eng. J. 2019, 392, 123745. [CrossRef] 29. Corbet, C.M.; McClellan, C.; Rai, A.; Sonde, S.S.; Tutuc, E.; Banerjee, S.K. Field effect transistors with current

saturation and voltage gain in ultrathin ReS2. ACS Nano 2015, 9, 363. [CrossRef] 30. Zhang, E.; Jin, Y.; Yuan, X.; Wang, W.; Zhang, C.; Tang, L.; Liu, S.; Zhou, P.; Hu, W.; Xiu, F. ReS2-based field-effect transistors and photodetectors. Adv. Funct. Mater. 2015, 25, 4076. [CrossRef] Nanomaterials 2020, 10, 2103 12 of 13

31. Liu, F.; Zheng, S.; He, X.; Chaturvedi, A.; He, J.; Chow, W.L.; Mion, T.R.; Wang, X.; Zhou, J.; Fu, Q.; et al.

Highly sensitive detection of polarized light using anisotropic 2D ReS2. Adv. Funct. Mater. 2016, 26, 1169. [CrossRef] 32. Zhang, Q.; Tan, S.; Mendes, R.G.; Sun, Z.; Chen, Y.; Kong, X.; Xue, Y.; Rummeli, M.H.; Wu, X.; Chen, S.; et al.

Extremely weak van der Waals coupling in vertical ReS2 nanowalls for high-current-density -ion batteries. Adv. Mater. 2016, 28, 2616. [CrossRef][PubMed] 33. Shim, J.; Oh, S.; Kang, D.H.; Jo, S.H.; Ali, M.H.; Choi, W.Y.; Heo, K.; Jeon, J.; Lee, S.; Kim, M.; et al. Phosphorene/rhenium disulfide heterojunction-based negative differential resistance device for multi-valued logic. Nat. Commun. 2016, 7, 13413. [CrossRef][PubMed] 34. Cazzanelli, E.; Castriota, M.; Marino, S.; Scaramuzza, N.; Purans, J.; Kuzmin, A.; Kalendarev, R.; Mariotto, G.; Das, G. Characterization of rhenium oxide films and their application to liquid crystal cells. J. Appl. Phys. 2009, 105, 114904. [CrossRef] 35. Xu, H.; Yang, L.Y.; Tian, H.; Yin, S.G.; Zhang, F.L. Rhenium oxide as the interfacial buffer layer for polymer photovoltaic cells. Optoelectron. Lett. 2010, 6, 176. [CrossRef]

36. Pearsall, T.P.; Lee, C.A. Electronic transport in ReO3: Dc conductivity and Hall effect. Phys. Rev. B 1974, 10, 2190. [CrossRef] 37. Kroeger,F.A. The Chemistry of Imperfect Crystals; North-Holland Publishing Company: Amsterdam, The Netherland, 1973. 38. Kozicki, M.N.; Mitkova, M.; Valov, I. Resistive Switching: From fundamentals of nanoionic redox processes to memristive device applications, Electrochemical metallization memories. In Resistive Switching: From Fundamentals of Nanoionic Redox Processes to Memristive Device Applications; Wiley-VCH: Weinheim, Germany, 2016; Chapter 17. 39. Sze, S.M. Physics of Semiconductor Devices; John Wiley & Sons: Hoboken, NJ, USA, 1981. 40. Silinsh, E.A. On the physical nature of traps in molecular crystals. Phys. Stat. Sol. A 1970, 3, 817. [CrossRef] 41. Sim, H.; Choi, H.; Lee, D.; Chang, M.; Choi, D.; Son, Y.; Lee, E.H.; Kim, W.; Park, Y.D.; Yoo, I.K.; et al.

Resistance-switching characteristics of polycrystalline Nb2O5 for nonvolatile memory application. Electron Devices Meeting IEDM Tech. Dig. 2005, 26, 292.

42. Jeong, D.S.; Schroeder, H.; Breuer, U.; Waser, R. Characteristic electroforming behavior in Pt/TiO2/Pt resistive switching cells depending on atmosphere. J. Appl. Phys. 2008, 104, 123716. [CrossRef] 43. Yang, Y.; Zhang, X.; Qin, L.; Zheng, Q.; Qiu, X.; Huang, R. Probing nanoscale oxygen ion motion in memristive systems. Nat. Commun. 2017, 8, 15173. [CrossRef][PubMed] 44. Serway, R.A. Principles of Physics; Saunders College Publishing: Texas, TX, USA, 1998. 45. Liu, M.; Borowiec, J.; Sun, L.J.J.; Konop, M.; Rahman, M.M.; Taallah, A.; Boi, F.S.; Gillin, W.P. Experimental studies on the conduction mechanism and electrical properties of the inverted Ba doped ZnO nanoparticles based memristor. Appl. Phys. Lett. 2019, 115, 073505. [CrossRef] 46. O’Dwyer, J.J. The Theory of Electrical Conduction and Breakdown in Solid Dielectrics; Clarendon: Oxford, UK, 1973. 47. Chua, L.O. Resistance switching memories are memristors. Appl. Phys. A Mater. Sci. Process 2011, 102, 765. [CrossRef] 48. Tsuruoka, T.; Terabe, K.; Hasegawa, T.; Valov, I.; Waser, R.; Aono, M. Effects of moisture on the switching characteristics of oxide-based, gapless-type atomic switches. Adv. Funct. Mater. 2012, 22, 70. [CrossRef] 49. Corinto, F.; Civalleri, P.P.; Chua, L.O. A Theoretical Approach to Memristor Devices. IEEE J. Emerg. Sel. Topics Circuits Syst. 2015, 5, 123. [CrossRef] 50. Valov, I.; Linn, E.; Tappertzhofen, S.; Schmelzer, S.; van den Hurk, J.; Lentz, F.; Wasser, R. Nanobatteries in redox-based resistive switches require extension of memristor theory. Nat. Commun. 2013, 4, 1771. [CrossRef] [PubMed] 51. Kasap, S.O. Principles of Electronic Materials and Devices. McGraw-Hill, Co.: Boston, MA, USA, 2006; Chapter 7. 52. Lai, Q.; Zhang, L.; Li, Z.; Stickle, W.F.; Williams, R.S.; Chen, Y. Analog memory capacitor based on field-configurable ion-doped polymers. Appl. Phys. Lett. 2009, 95, 213503. [CrossRef] 53. Lee, E.; Gwon, M.; Kim, D.W.; Kim, H. Resistance state-dependent barrier inhomogeneity and transport

mechanism in resistive-switching Pt/SrTiO3 junctions. Appl. Phys. Lett. 2011, 98, 132905. [CrossRef] 54. Wu, S.X.; Xu, L.M.; Xing, X.J.; Chen, S.M.; Yuan, Y.B.; Liu, Y.J.; Yu, Y.P.; Li, X.Y.; Li, S.W. Reverse-bias-induced

bipolar resistance switching in Pt/TiO2/SrTi0.99Nb0.01O3/Pt devices. Appl. Phys. Lett. 2008, 93, 043502. [CrossRef] Nanomaterials 2020, 10, 2103 13 of 13

55. Sawa, A.; Fujii, T.; Kawasaki, M.; Tokura, Y. Interface transport properties and resistance switching in perovskite-oxide heterojunctions, Strongly correlated electron materials: Physics and nanoengineering. Proc. SPIE 2005, 5932, 59322C. 56. Sawa, A.; Fujii, T.; Kawasaki, T. Interface transport properties and resistance switching in perovskite-oxide heterojunctions. Appl. Phys. Lett. 2006, 88, 232112. [CrossRef] 57. Nicollian, E.H.; Brews, J.R. Metal Oxide Semiconductor (MOS) Physics and Technology; John Wiley and Sons: Hoboken, NJ, USA, 1982. 58. Bonafos, C.; Spiegel, Y.; Normand, P.;Ben-Assayag, G.; Groenen, J.; Carrada, M.; Dimitrakis, P.;Kapetanakis, E.;

Sahu, B.S.; Slaoui, A.; et al. Controlled fabrication of Si nanocrystal delta-layers in thin SiO2 layers by plasma immersion ion implantation for nonvolatile memories. Appl. Phys. Lett. 2013, 103, 253118. [CrossRef] 59. Misawa, T. Impedance of bulk semiconductor in junction diode. J. Phys. Soc. Jpn. 1957, 12, 882. [CrossRef] 60. Wang, Q.; Chen, J.; Tang, H.; Li, X. Anomalous capacitance in temperature and frequency characteristics of a TiW/P-InP Schottky barrier. Semicond. Sci. Technol. 2016, 31, 065023. [CrossRef]

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).