<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Neurosci.</journal-id>
<journal-title>Frontiers in Neuroscience</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Neurosci.</abbrev-journal-title>
<issn pub-type="epub">1662-453X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fnins.2022.868671</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Neuroscience</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Voltage&#x2013;Time Transformation Model for Threshold Switching Spiking Neuron Based on Nucleation Theory</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Yap</surname> <given-names>Suk-Min</given-names></name>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Wang</surname> <given-names>I-Ting</given-names></name>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1745149/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Wu</surname> <given-names>Ming-Hung</given-names></name>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Hou</surname> <given-names>Tuo-Hung</given-names></name>
<xref ref-type="corresp" rid="c002"><sup>&#x002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/163677/overview"/>
</contrib>
</contrib-group>
<aff><institution>Department of Electrical Engineering and Institute of Electronics, National Yang Ming Chiao Tung University</institution>, <addr-line>Hsinchu</addr-line>, <country>Taiwan</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Shimeng Yu, Georgia Institute of Technology, United States</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Can Li, The University of Hong Kong, Hong Kong SAR, China; Jiyong Woo, Kyungpook National University, South Korea</p></fn>
<corresp id="c001">&#x002A;Correspondence: I-Ting Wang, <email>itwang@nycu.edu.tw</email></corresp>
<corresp id="c002">Tuo-Hung Hou, <email>thhou@mail.nctu.edu.tw</email></corresp>
<fn fn-type="other" id="fn004"><p>This article was submitted to Neuromorphic Engineering, a section of the journal Frontiers in Neuroscience</p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>13</day>
<month>04</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>16</volume>
<elocation-id>868671</elocation-id>
<history>
<date date-type="received">
<day>03</day>
<month>02</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>09</day>
<month>03</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2022 Yap, Wang, Wu and Hou.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Yap, Wang, Wu and Hou</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract>
<p>In this study, we constructed a voltage&#x2013;time transformation model (V&#x2013;t Model) to predict and simulate the spiking behavior of threshold-switching selector-based neurons (TS neurons). The V&#x2013;t Model combines the physical nucleation theory and the resistor&#x2013;capacitor (RC) equivalent circuit and successfully depicts the history-dependent threshold voltage of TS selectors, which has not yet been modeled in TS neurons. Moreover, based on our model, we analyzed the currently reported TS devices, including ovonic threshold switching (OTS), insulator-metal transition, and silver- (Ag-) based selectors, and compared the behaviors of the predicted neurons. The results suggest that the OTS neuron is the most promising and potentially achieves the highest spike frequency of GHz and the lowest operating voltage and area overhead. The proposed V&#x2013;t Model provides an engineering pathway toward the future development of TS neurons for neuromorphic computing applications.</p>
</abstract>
<kwd-group>
<kwd>threshold switching selector</kwd>
<kwd>spiking neuron</kwd>
<kwd>nucleation theory</kwd>
<kwd>history-dependent</kwd>
<kwd>neuromorphic computing</kwd>
</kwd-group>
<contract-num rid="cn001">109-2221-E-009-020-MY3</contract-num>
<contract-num rid="cn001">110-2221-E-A49-088-MY3</contract-num>
<contract-num rid="cn001">110-2634-F-009-017</contract-num>
<contract-num rid="cn001">110-2622-8-009-018-SB</contract-num>
<contract-sponsor id="cn001">Ministry of Science and Technology, Taiwan<named-content content-type="fundref-id">10.13039/501100004663</named-content></contract-sponsor>
<counts>
<fig-count count="8"/>
<table-count count="1"/>
<equation-count count="5"/>
<ref-count count="22"/>
<page-count count="7"/>
<word-count count="4719"/>
</counts>
</article-meta>
</front>
<body>
<sec id="S1" sec-type="intro">
<title>Introduction</title>
<p>With the increasing demand for massive data storage and processing, conventional computing systems based on the von-Neumann architecture have encountered their limitations. Frequent data transition between the separated processor and memory units makes conventional computation less efficient. Recently, emerging neuromorphic computing is regarded as the next-generation computing paradigm. Unlike the conventional von-Neumann-based computing system, brain-inspired neuromorphic computing not only provides energy-efficient computation with high parallelism but also shortens the latency of data transmission by realizing in-memory computing within crossbar memory arrays (<xref ref-type="bibr" rid="B6">Ielmini and Wong, 2018</xref>; <xref ref-type="bibr" rid="B5">Hua et al., 2019</xref>; <xref ref-type="bibr" rid="B18">Woo et al., 2019</xref>). In a neuromorphic computing system, an artificial synapse provides an adjustable and long-lasting weight value. In addition, an artificial neuron integrates and processes signals from synapses and then transmits the processed signals to the next neural layer as inputs. Both synapses and neurons have been extensively studied based on solid-state devices for neuromorphic hardware implementation (<xref ref-type="bibr" rid="B8">Lee et al., 2019a</xref>; <xref ref-type="bibr" rid="B18">Woo et al., 2019</xref>; <xref ref-type="bibr" rid="B22">Zhang et al., 2020</xref>). However, the conventional complementary metal-oxide semiconductor- (CMOS-)based neuron circuit occupies large chip areas because it requires a large number of transistors and capacitors for generating spike signals. In contrast, the neuron circuit area can be 10 times smaller by using novel devices, such as magnetoresistance memory (MRAM) (<xref ref-type="bibr" rid="B19">Wu et al., 2019</xref>, <xref ref-type="bibr" rid="B20">2020</xref>; <xref ref-type="bibr" rid="B11">Liang et al., 2020</xref>), phase-change memory (PCM) (<xref ref-type="bibr" rid="B16">Tuma et al., 2016</xref>), and threshold switching (TS) selector (<xref ref-type="bibr" rid="B14">Park et al., 2016</xref>; <xref ref-type="bibr" rid="B15">Song et al., 2018</xref>; <xref ref-type="bibr" rid="B3">Grisafe et al., 2019</xref>; <xref ref-type="bibr" rid="B4">Hatem et al., 2019</xref>; <xref ref-type="bibr" rid="B5">Hua et al., 2019</xref>), which is beneficial for ultrahigh density neuromorphic computing applications (<xref ref-type="bibr" rid="B12">Liang et al., 2021</xref>).</p>
<p>Among several novel device-based neurons, threshold-switching selector-based neurons (TS neurons) are especially promising for ultra-high density neuromorphic architectures due to their simpler and smaller neuronal circuits (<xref ref-type="bibr" rid="B12">Liang et al., 2021</xref>). A circuit-level model solving Kirchhoff&#x2019;s Law based on the resistor&#x2013;capacitor (RC) equivalent circuit has been proposed to describe the behavior of TS neurons (RC Model) (<xref ref-type="bibr" rid="B2">Chen et al., 2016</xref>; <xref ref-type="bibr" rid="B17">Wang et al., 2020</xref>). However, the RC Model oversimplified the TS neuron by assuming constant switching behavior of the TS selector. Indeed, the switching dynamics of the real TS selector is affected by the external electric field, which can be explained using the nucleation theory (<xref ref-type="bibr" rid="B7">Karpov et al., 2008</xref>; <xref ref-type="bibr" rid="B10">Lee et al., 2020</xref>). Specifically, the way the external electric field is previously accumulated determines the device behavior, and we regard this time-dependent phenomenon as history dependence. Consequently, the TS voltage (<italic>V</italic><sub>th</sub>) in the TS selector is history dependent rather than constant. In this study, aiming for constructing a more comprehensive and accurate neuron model, we proposed an improved voltage&#x2013;time transformation model (V&#x2013;t Model) on top of the original RC Model by considering the TS behavior both experimentally and theoretically.</p>
<p>In the following sections, we will first verify the spiking behavior of the TS neuron according to different synaptic weights. A silver- (Ag-)based TS selector was chosen to observe the switching dynamics and the history-dependent <italic>V</italic><sub>th</sub> of the device. Additionally, based on the nucleation theory, we will introduce a V&#x2013;t transformation (V&#x2013;t) equation to describe the variant <italic>V</italic><sub>th</sub> of the TS selector, and a V&#x2013;t Model will be constructed. Furthermore, several types of TS neurons based on the reported TS selectors, including ovonic threshold switching (OTS) (<xref ref-type="bibr" rid="B15">Song et al., 2018</xref>; <xref ref-type="bibr" rid="B4">Hatem et al., 2019</xref>), insulator&#x2013;metal transition (IMT) (<xref ref-type="bibr" rid="B14">Park et al., 2016</xref>), and Ag-based selectors (<xref ref-type="bibr" rid="B3">Grisafe et al., 2019</xref>; <xref ref-type="bibr" rid="B5">Hua et al., 2019</xref>), will be evaluated. The results suggest that the OTS neuron has the fastest spike frequency and a lower history-dependent <italic>V</italic><sub>th</sub>. The V&#x2013;t Model not only successfully depicts and predicts the characteristics of TS neurons, but it also provides a useful engineering guideline for future high-performance neuron circuits for neuromorphic computing applications.</p>
</sec>
<sec id="S2">
<title>Experimental Details and Measurement Setup</title>
<sec id="S2.SS1">
<title>Ag-Based Threshold Switching Selector</title>
<p>In this study, an Ag/hafnium oxide (HfO<sub>x</sub>)/Pt TS selector was fabricated and investigated. The schematic illustration of the Ag-based TS selector is shown in <xref ref-type="fig" rid="F1">Figure 1A</xref>. The Pt bottom layer was first deposited on a Ti/Si substrate using electron beam evaporation, followed by the silicon dioxide (SiO<sub>2</sub>) layer deposition using plasma-enhanced chemical vapor deposition. After the photolithography process, the reactive ion etching of SiO<sub>2</sub> was applied to form a <italic>via</italic> contact with a diameter of 1 &#x03BC;m, which defines the effective device area. Then, the 4.5-nm-insulating HfO<sub>x</sub> layer was deposited using atomic layer deposition. After that, 2-nm-thick Ag was deposited on the HfO<sub>x</sub> layer using electron beam evaporation followed by rapid thermal annealing (RTA) at 500&#x00B0;C for 5 min to form Ag nanoparticles (NPs) as the active electrode. Finally, a 60-nm-thick Ni capping layer was deposited using electron beam deposition to prevent the oxidation of Ag NPs. Electrical measurements were performed using an Agilent B1500A and B1530A waveform generation/fast measurement unit at room temperature. <xref ref-type="fig" rid="F1">Figure 1B</xref> shows the scanning electron microscope (SEM) image of Ag NPs. The size distribution of NPs is shown in the inset. <xref ref-type="fig" rid="F1">Figure 1C</xref> shows the DC current-voltage (I-V) characteristics of the Ag/HfO<sub>x</sub>/Pt TS selector with 500 DC cycles of TS and a compliance current (<italic>I</italic><sub>cc</sub>) of 0.1 mA. The device provides an extremely high on/off ratio (&#x223C;10<sup>9</sup>) and small <italic>V</italic><sub>th</sub> and hold voltage (<italic>V</italic><sub>hold</sub>) for both positive and negative bias, showing typical behaviors of Ag-based TS selectors as reported in the literature (<xref ref-type="bibr" rid="B21">Yoo et al., 2017</xref>).</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption><p>Silver- (Ag-) based threshold switching (TS) selector studied in this study. <bold>(A)</bold> The schematic structure showing the cross-section view of the Ag/hafnium oxide (HfO<sub>x</sub>)/Pt device. <bold>(B)</bold> The top-view scanning electron microscope (SEM) image of Ag nanoparticles (NPs) with an average diameter of 30 nm formed by 500&#x00B0;C rapid thermal annealing (RTA) for 5 min. The inset shows the size distribution plot of the Ag NPs. <bold>(C)</bold> DC current&#x2013;voltage (I&#x2013;V) characteristics of the device with a current compliance of 0.1 mA, showing the success of TS for 500 DC cycles.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-16-868671-g001.tif"/>
</fig>
</sec>
<sec id="S2.SS2">
<title>Threshold Switching Neuron Circuit</title>
<p>To emulate neuromorphic hardware in which synapses and neurons are connected in the neural network (<xref ref-type="fig" rid="F2">Figure 2A</xref>), the measurement setup adopted in this study is illustrated in <xref ref-type="fig" rid="F2">Figure 2B</xref>. The effective resistor connected in series (<italic>R</italic><sub>series</sub>) represents the total resistance of multiple synaptic devices in the synaptic array connecting in parallel to the same TS neuron. A parasitic capacitor (<italic>C</italic><sub>parasitic</sub>) of the TS selector is exploited; therefore, no extra capacitor is needed for signal integration. The evolution of the total current (<italic>I</italic><sub>total</sub>) flowing through <italic>R</italic><sub>series</sub> and the voltage across TS selector (<italic>V</italic><sub>selector</sub>) is described in <xref ref-type="fig" rid="F2">Figure 2C</xref>: when a constant input voltage (<italic>V</italic><sub>input</sub>) is applied to the neuron circuit, most of the voltage initially drops across the TS selector in the off-state. Then, <italic>V</italic><sub>selector</sub> is gradually increased by charging <italic>C</italic><sub>parasitic</sub>. Once <italic>V</italic><sub>selector</sub> reaches <italic>V</italic><sub>th</sub>, the TS selector is switched to the on-state due to the formation of a volatile conducting filament, and an increase in <italic>I</italic><sub>total</sub> can be observed. However, <italic>V</italic><sub>selector</sub> drops right after the TS selector is switched to the on-state due to the discharge of <italic>C</italic><sub>parasitic</sub>, and <italic>I</italic><sub>total</sub> starts to decrease. The TS selector returns to the off-state when <italic>V</italic><sub>selector</sub> reduces to <italic>V</italic><sub>hold</sub> because of the rupture of the volatile conducting filament. <italic>t</italic><sub>on</sub> and <italic>t</italic><sub>off</sub> define the required period of time for the selector to be turned on (<italic>V</italic><sub>selector</sub> to increase from <italic>V</italic><sub>hold</sub> to <italic>V</italic><sub>th</sub>) and off (<italic>V</italic><sub>selector</sub> to decrease from <italic>V</italic><sub>th</sub> to <italic>V</italic><sub>hold</sub>) in the neuron circuit, respectively. When the circuit is biased, a series of continuous current and voltage spikes are generated, and the spike frequency can be calculated as the number of spikes per second (Hz) accordingly. To fulfill the requirement of neural network applications, artificial neurons should be capable of generating different spike frequencies according to the weights of connected synapses, i.e., <italic>R</italic><sub>series</sub>. In the RC Model (<xref ref-type="bibr" rid="B2">Chen et al., 2016</xref>; <xref ref-type="bibr" rid="B17">Wang et al., 2020</xref>), the <italic>t</italic><sub>on</sub> in the TS neuron circuit is obtained by</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption><p><bold>(A)</bold> Inspired by the biological neural network, neurons are connected with synapses through bit lines and word lines in the crossbar memory array. <italic>V</italic><sub>input</sub> received from the pre-neurons is applied on the word lines, and the post-neurons connected to bit lines generate output spikes according to the weight of synapses. <bold>(B)</bold> Equivalent circuit of the measurement setup where the synaptic devices along the same bit line (with a total resistance of <italic>R</italic><sub>series</sub>) are connected in series with a TS selector as a neuron. A parasitic capacitor <italic>C</italic><sub>parasitic</sub> is considered and is connected in parallel with the TS neuron. <bold>(C)</bold> The evolution of <italic>I</italic><sub>total</sub> and <italic>V</italic><sub>selector</sub> shows the spiking behavior of the TS neuron when <italic>V</italic><sub>input</sub> is applied due to the charging and discharging of <italic>C</italic><sub>parasitic</sub>. <italic>t</italic><sub>on</sub> and <italic>t</italic><sub>off</sub> define the required period of time for the TS selector to turned on and off in the neuron circuit, respectively.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-16-868671-g002.tif"/>
</fig>
<disp-formula id="S2.E1"><label>(1)</label><mml:math display="block" id="M1"><mml:mrow><mml:mpadded width="+3.3pt"><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mtext>on</mml:mtext></mml:mrow></mml:msub></mml:mpadded><mml:mo rspace="5.8pt">=</mml:mo><mml:mrow><mml:mo>-</mml:mo><mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mtext>series</mml:mtext></mml:mrow></mml:msub><mml:mo>&#x2062;</mml:mo><mml:mpadded width="+3.3pt"><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mtext>series</mml:mtext></mml:mrow></mml:msub></mml:mpadded></mml:mrow><mml:mo rspace="5.8pt">&#x00D7;</mml:mo><mml:mtext>ln</mml:mtext></mml:mrow><mml:mo>&#x2062;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mo>|</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mtext>input</mml:mtext></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mtext>th</mml:mtext></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mtext>input</mml:mtext></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mtext>hold</mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>|</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mrow></mml:mrow></mml:math></disp-formula>
<p>In this study, we assume that the IR voltage drop on <italic>R</italic><sub>series</sub> is negligible when the TS selector is at its off-state due to the low leakage current (below pA in our case). <xref ref-type="fig" rid="F3">Figure 3A</xref> illustrates the statistically measured <italic>t</italic><sub>on</sub> of the TS neuron when connecting to different <italic>R</italic><sub>series</sub>, and the inset is an example of experimentally obtained current spikes (<italic>I</italic><sub>total</sub>) when <italic>R</italic><sub>series</sub> is 3,300 k&#x03A9;. <xref ref-type="fig" rid="F3">Figure 3B</xref> presents the calculated spike frequency, as shown in <xref ref-type="fig" rid="F3">Figure 3A</xref>. The results indicate that, with the decrease of <italic>R</italic><sub>series</sub>, <italic>t</italic><sub>on</sub> is decreased and the spike frequency is increased accordingly. However, the spike frequency cannot be further increased when <italic>R</italic><sub>series</sub> is &#x003C; 100 k&#x03A9;. It is worth mentioning that <italic>R</italic><sub>series</sub> in the neuron circuit also acts as current compliance, where it controls the morphology and the size of conducting filaments in the TS selector (<xref ref-type="bibr" rid="B1">Chae et al., 2017</xref>). If <italic>R</italic><sub>series</sub> is too small, the filaments of extremely large size become non-volatile and cannot be ruptured even at <italic>V</italic><sub>hold</sub> = 0 V. Consequently, <italic>t</italic><sub>off</sub> increases and limits the spike frequency due to the difficult dissolution of large-size filaments in the TS selector. As a result, the resistance range of <italic>R</italic><sub>series</sub> requires careful adjustment (&#x003E; 100 k&#x03A9; in our case) to prevent the dysfunction of neuron circuits. With a suitable range of <italic>R</italic><sub>series</sub> and with <italic>t</italic><sub>off</sub> being much smaller than <italic>t</italic><sub>on</sub>, the spike frequency is the inverse of <italic>t</italic><sub>on</sub>, thus proportional to the inverse of <italic>R</italic><sub>series</sub>, i.e., the effective total conductance of the synaptic array.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption><p><bold>(A)</bold> The measured <italic>t</italic><sub>on</sub> increases with increasing <italic>R</italic><sub>series</sub> while the calculated spike frequency shown in panel <bold>(B)</bold> decreases with increasing <italic>R</italic><sub>series</sub> in the neuron circuit. The inset in panel <bold>(A)</bold> shows an example of the experimentally obtained spike current (<italic>I</italic><sub>total</sub>) when <italic>R</italic><sub>series</sub> is 3,300 k&#x03A9;. The spike frequency is defined by the number of spikes per second (Hz). The spike frequency is approximately equal to the inverse of <italic>t</italic><sub>on</sub> when <italic>R</italic><sub>series</sub> is greater than 100 k&#x03A9; and <italic>t</italic><sub>on</sub> is much larger than <italic>t</italic><sub>off</sub>.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-16-868671-g003.tif"/>
</fig>
</sec>
</sec>
<sec id="S3" sec-type="results|discussion">
<title>Results and Discussion</title>
<sec id="S3.SS1">
<title>History-Dependent <italic>V</italic><sub>th</sub> of the Threshold Switching Selector in Nueron Circuit</title>
<p>An important assumption of the RC Model in Equation 1 is that the <italic>V</italic><sub>th</sub> of the TS selector is constant. <xref ref-type="fig" rid="F4">Figure 4A</xref> compares the measured <italic>V</italic><sub>th</sub> captured by an oscilloscope with <italic>R</italic><sub>series</sub> of 150 and 470 k&#x03A9;, and the statistical results are indicated in <xref ref-type="fig" rid="F4">Figure 4B</xref>. Instead of remaining constant, the <italic>V</italic><sub>th</sub> of the TS selector varies with <italic>R</italic><sub>series</sub>. Different <italic>R</italic><sub>series</sub> modulate the charging rate of <italic>V</italic><sub>selector</sub> and give rise to the history-dependent <italic>V</italic><sub>th</sub>. The na&#x00EF;ve RC model does not consider the history-dependent <italic>V</italic><sub>th</sub> of the TS selector, thus reducing the accuracy and prediction capability of the neuron model.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption><p><bold>(A)</bold> The oscilloscope waveform of <italic><underline>V</underline></italic><sub>selector</sub> when the TS selector is connected with <italic>R</italic><sub>series</sub> of 150 and 470 k&#x03A9; and <italic>V</italic><sub>input</sub> = 2 V. Corresponding <italic>V</italic><sub>th</sub> is also indicated. <bold>(B)</bold> Statistically measured <italic>V</italic><sub>th</sub> increases with decreasing <italic>R</italic><sub>series</sub>. Instead of remaining constant, the history-dependent <italic>V</italic><sub>th</sub> needs to be carefully considered in the TS neuron model.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-16-868671-g004.tif"/>
</fig>
</sec>
<sec id="S3.SS2">
<title>Voltage&#x2013;Time Transformation Model</title>
<p>To include the characteristic of history-dependent <italic>V</italic><sub>th</sub> into the neuron model, the V&#x2013;t Model is proposed. Starting from considering the switching dynamics of TS selectors when constant voltage stress (V<sub>CVS</sub>) is applied directly on the device, i.e., <italic>V</italic><sub>selector</sub> equals to <italic>V</italic><sub>CVS</sub>. This is the case similar to <xref ref-type="fig" rid="F2">Figure 2B</xref> but without the external <italic>R</italic><sub>series</sub>. The time delay before turning on the selector (<italic>t</italic><sub>on_CVS</sub>) is determined by the nucleation theory (<xref ref-type="bibr" rid="B7">Karpov et al., 2008</xref>; <xref ref-type="bibr" rid="B10">Lee et al., 2020</xref>):</p>
<disp-formula id="S3.E2"><label>(2)</label><mml:math display="block" id="M2"><mml:mrow><mml:mpadded width="+3.3pt"><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mtext>on_CVS</mml:mtext></mml:mrow></mml:msub></mml:mpadded><mml:mo rspace="5.8pt">=</mml:mo><mml:mrow><mml:msub><mml:mi mathvariant="normal">&#x03C4;</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>&#x2062;</mml:mo><mml:mrow><mml:mi>exp</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>W</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>&#x2062;</mml:mo><mml:msup><mml:mi mathvariant="normal">&#x03B1;</mml:mi><mml:mfrac><mml:mn>3</mml:mn><mml:mn>2</mml:mn></mml:mfrac></mml:msup><mml:mo>&#x2062;</mml:mo><mml:msub><mml:mi>E</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>&#x2062;</mml:mo><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mi>k</mml:mi><mml:mo>&#x2062;</mml:mo><mml:mi>T</mml:mi><mml:mo>&#x2062;</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mtext>CVS</mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mrow></mml:mrow></mml:math></disp-formula>
<p>where &#x03C4;<sub>0</sub> is the intrinsic time constant of the device, <italic>W</italic><sub>0</sub> is the nucleation barrier energy without electric field, &#x03B1; is a geometric factor of a nucleus, <italic>E</italic><sub>0</sub> is the voltage acceleration factor, <italic>d</italic> is the effective thickness of the insulating layer, <italic>k</italic> is Boltmann&#x2019;s constant, and <italic>T</italic> is the ambient temperature. We define &#x1D538; as a material-related constant at a fixed <italic>T</italic>, and (2) can be rewritten as</p>
<disp-formula id="S3.E3"><label>(3)</label><mml:math display="block" id="M3"><mml:mrow><mml:mpadded width="+3.3pt"><mml:mi>&#x1D538;</mml:mi></mml:mpadded><mml:mo rspace="5.8pt">=</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mtext>CVS</mml:mtext></mml:mrow></mml:msub><mml:mo>&#x22C5;</mml:mo><mml:mtext>ln</mml:mtext></mml:mrow><mml:mo>&#x2062;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mtext>on_CVS</mml:mtext></mml:mrow></mml:msub><mml:msub><mml:mi mathvariant="normal">&#x03C4;</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mfrac><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>W</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>&#x2062;</mml:mo><mml:msup><mml:mi mathvariant="normal">&#x03B1;</mml:mi><mml:mfrac><mml:mn>3</mml:mn><mml:mn>2</mml:mn></mml:mfrac></mml:msup><mml:mo>&#x2062;</mml:mo><mml:msub><mml:mi>E</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>&#x2062;</mml:mo><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mi>k</mml:mi><mml:mo>&#x2062;</mml:mo><mml:mi>T</mml:mi></mml:mrow></mml:mfrac></mml:mrow></mml:math></disp-formula>
<p>where &#x1D538; and &#x03C4;<sub>0</sub> can be obtained from fitting the measured <italic>V</italic><sub>CVS</sub> and <italic>t</italic><sub>on_CVS</sub>. We assume &#x1D538; remains constant when measuring the same device. Therefore, the V&#x2013;t equation can be used to describe the transformation relation between any two arbitrary CVS voltages, <italic>V</italic><sub>CVS1</sub> and <italic>V</italic><sub>CVS2</sub>, and their corresponding turn-on times, <italic>t</italic><sub>on_CVS1</sub> and <italic>t</italic><sub>on_CVS2</sub> as</p>
<disp-formula id="S3.E4"><label>(4)</label><mml:math display="block" id="M4"><mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi>CVS1</mml:mi></mml:msub><mml:mo>&#x22C5;</mml:mo><mml:mtext>ln</mml:mtext></mml:mrow><mml:mo>&#x2062;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mrow><mml:mtext>on</mml:mtext></mml:mrow><mml:mo>&#x2062;</mml:mo><mml:mi mathvariant="normal">_</mml:mi><mml:mo>&#x2062;</mml:mo><mml:mi>CVS1</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mi mathvariant="normal">&#x03C4;</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mfrac><mml:mo rspace="5.8pt">)</mml:mo></mml:mrow></mml:mrow><mml:mo rspace="5.8pt">=</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mrow><mml:mtext>CVS</mml:mtext></mml:mrow><mml:mo>&#x2062;</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>&#x22C5;</mml:mo><mml:mtext>ln</mml:mtext></mml:mrow><mml:mo>&#x2062;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mrow><mml:mtext>on</mml:mtext></mml:mrow><mml:mo>&#x2062;</mml:mo><mml:mi mathvariant="normal">_</mml:mi><mml:mo>&#x2062;</mml:mo><mml:mi>CVS2</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mi mathvariant="normal">&#x03C4;</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mfrac><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mrow></mml:math></disp-formula>
<p>When connecting the TS selector with <italic>R</italic><sub>series</sub> to form a complete TS neuron circuit, as shown in <xref ref-type="fig" rid="F2">Figure 2B</xref>, <italic>V</italic><sub>selector</sub> becomes time-varying according to the RC equivalent circuit. The time-varying <italic>V</italic><sub>selector</sub> could be approximated using a finite number of CVS steps as depicted in <xref ref-type="fig" rid="F5">Figure 5</xref>, which increase from (<italic>t</italic><sub>1</sub>, <italic>V</italic><sub>1</sub>) to (<italic>t</italic><sub>2</sub>, <italic>V</italic><sub>2</sub>) and eventually to (<italic>t</italic><sub>on</sub>, <italic>V</italic><sub>th</sub>) indicated by the blue line. &#x0394;<italic>V</italic> and &#x0394;<italic>t</italic> are the voltage and time intervals, respectively. A similar conversion between CVS and ramp voltage stress has been reported and validated in resistive switching memory devices (<xref ref-type="bibr" rid="B13">Luo et al., 2013</xref>). As indicated in <xref ref-type="fig" rid="F5">Figure 5</xref>, the stress effect of the (<italic>t</italic><sub>1</sub>, <italic>V</italic><sub>1</sub>) step indicated by the blue-filled rectangle on the device is transformed to that of an equivalent (<italic>t</italic>&#x2032;<sub>1</sub>, <italic>V</italic><sub>2</sub>) step indicated by the red dashed rectangle based on (4), <italic>t</italic>&#x2032;<sub>1</sub> is therefore expressed as:</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption><p>Time-varying <italic>V</italic><sub>selector</sub> in the neuron circuit is approximated using a finite number of constant voltage stress (CVS) steps from (<italic>t</italic><sub>1</sub>, <italic>V</italic><sub>1</sub>) to (<italic>t</italic><sub>2</sub>, <italic>V</italic><sub>2</sub>) and eventually to (<italic>t</italic><sub>on</sub>, <italic>V</italic><sub>th</sub>). &#x0394;<italic>V</italic> and &#x0394;<italic>t</italic> determine the voltage and time intervals, respectively. Based on the proposed V&#x2013;t Model, the transformed (<italic>t</italic><sub>1</sub>, <italic>V</italic><sub>2</sub>) step indicated by the red-dashed rectangle is equivalent to the (<italic>t</italic><sub>1</sub>, <italic>V</italic><sub>1</sub>) step indicated by the blue-filled rectangle. The new <inline-formula><mml:math id="INEQ2"><mml:mover accent="true"><mml:msub><mml:mtext mathvariant="bold">t</mml:mtext><mml:mn mathvariant="bold">2</mml:mn></mml:msub><mml:mo>&#x00AF;</mml:mo></mml:mover></mml:math></inline-formula> of <italic>V</italic><sub>2</sub> is now <italic>t</italic>&#x2032;<sub>1</sub>+&#x0394;<italic>t</italic>, which includes the history effect of the previously accumulated (<italic>t</italic><sub>1</sub>, <italic>V</italic><sub>1</sub>) step.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-16-868671-g005.tif"/>
</fig>
<disp-formula id="S3.E5"><label>(5)</label><mml:math display="block" id="M5"><mml:mrow><mml:mpadded width="+3.3pt"><mml:msubsup><mml:mi>t</mml:mi><mml:mn>1</mml:mn><mml:mo>&#x2032;</mml:mo></mml:msubsup></mml:mpadded><mml:mo rspace="5.8pt">=</mml:mo><mml:mrow><mml:mi>exp</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mrow><mml:mrow><mml:mfrac><mml:msub><mml:mi>V</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:msub><mml:mi>V</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mfrac><mml:mo>&#x22C5;</mml:mo><mml:mtext>ln</mml:mtext></mml:mrow><mml:mo>&#x2062;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:msub><mml:mi mathvariant="normal">&#x03C4;</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mfrac><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:mtext>ln</mml:mtext><mml:mo>&#x2062;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi mathvariant="normal">&#x03C4;</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow></mml:mrow></mml:math></disp-formula>
<p>A new equivalent CVS step of (<inline-formula><mml:math id="INEQ10"><mml:mover accent="true"><mml:msub><mml:mtext mathvariant="bold">t</mml:mtext><mml:mn mathvariant="bold">2</mml:mn></mml:msub><mml:mo>&#x00AF;</mml:mo></mml:mover></mml:math></inline-formula>, <italic>V</italic><sub>2</sub>) with an equivalent stress time of <inline-formula><mml:math id="INEQ11"><mml:mover accent="true"><mml:msub><mml:mtext mathvariant="bold">t</mml:mtext><mml:mn mathvariant="bold">2</mml:mn></mml:msub><mml:mo>&#x00AF;</mml:mo></mml:mover></mml:math></inline-formula> = <italic>t</italic>&#x2032;<sub>1</sub> + &#x0394;<italic>t</italic> at <italic>V</italic><sub>2</sub> includes the history effect of the previous (<italic>V</italic><sub>1</sub>, <italic>t</italic><sub>1</sub>) step. This equivalent stress time is accumulated until it reaches the <italic>t</italic><sub>on_CVS</sub> at the stop voltage, i.e., <italic>V</italic><sub>th</sub>, which could be calculated by Equation 2. Under these circumstances, <italic>V</italic><sub>th</sub> becomes history-dependent and is affected by the RC charging process and <italic>R</italic><sub>series</sub>. The larger <italic>R</italic><sub>series</sub>, the lower <italic>V</italic><sub>th</sub>, and longer <italic>t</italic><sub>on</sub>.</p>
<p>To confirm the feasibility of the V&#x2013;t Model on the prediction of the TS neuron behavior, the simulation results obtained from the RC Model (<xref ref-type="bibr" rid="B2">Chen et al., 2016</xref>; <xref ref-type="bibr" rid="B17">Wang et al., 2020</xref>) and the proposed V&#x2013;t Model are compared in <xref ref-type="fig" rid="F6">Figures 6A,B</xref> with the measurement. The RC Model only describes the RC behavior of the neuron circuit with a constant <italic>V</italic><sub>th</sub> of the TS selector, therefore it not only underestimates <italic>t</italic><sub>on</sub> but also fails to depict the history-dependent <italic>V</italic><sub>th</sub> of the TS selector. In contrast, the proposed V&#x2013;t Model predicted well <italic>t</italic><sub>on</sub> and <italic>V</italic><sub>th</sub> of the TS selector under different <italic>R</italic><sub>series</sub>.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption><p>Measurement data and simulated results of <bold>(A)</bold> <italic>t</italic><sub>on</sub> and <bold>(B)</bold> <italic>V</italic><sub>th</sub> of the TS neuron with different <italic>R</italic><sub>series</sub>. The results suggest a more accurate prediction based on the V&#x2013;t Model than the RC Model (<xref ref-type="bibr" rid="B2">Chen et al., 2016</xref>; <xref ref-type="bibr" rid="B17">Wang et al., 2020</xref>) by considering the history-dependent <italic>V</italic><sub>th</sub> of the TS selector.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-16-868671-g006.tif"/>
</fig>
</sec>
<sec id="S3.SS3">
<title>Prediction of Threshold Switching Neuron Performance Based on V&#x2013;t Model</title>
<p>In this section, we explored the impact of the TS selector on TS neurons, and the effect of <italic>t</italic><sub>on</sub> and &#x1D538; on <italic>V</italic><sub>th</sub> can be predicted based on Equation 2. As shown in <xref ref-type="fig" rid="F7">Figure 7</xref>, when <italic>t</italic><sub>on</sub> approaches &#x03C4;<sub>0</sub>, the voltage required for nucleation (<italic>V</italic><sub>th</sub>) approaches infinity. In addition, the TS selector with larger &#x1D538; requires a higher <italic>V</italic><sub>th</sub> to be turned on. These results indicated that, under the same <italic>t</italic><sub>on</sub>, the TS selector with larger &#x03C4;<sub>0</sub> and &#x1D538; needs a higher applied voltage than the one with smaller &#x03C4;<sub>0</sub> and &#x1D538;. However, the required high applied voltage is not favorable because it not only may result in an irreversible breakdown of the device but also may increase the difficulty of circuit integration. Therefore, the TS selector with larger &#x03C4;<sub>0</sub> and &#x1D538; may require an additional external integration capacitor to maintain a reasonable <italic>V</italic><sub>th</sub>, which on the other hand increases the circuit footprint and <italic>t</italic><sub>on</sub> and decreases the spike frequency. The energy consumption per spike of the neuron circuit could also increase due to slow spiking (<xref ref-type="bibr" rid="B12">Liang et al., 2021</xref>).</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption><p>Relation between <italic>t</italic><sub>on</sub>, &#x03C4;<sub>0</sub>, &#x1D538;, and <italic>V</italic><sub>th</sub> of TS selectors based on the nucleation theory in Equation 2.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-16-868671-g007.tif"/>
</fig>
<p><xref ref-type="table" rid="T1">Table 1</xref> lists the reported parameters of &#x03C4;<sub>0</sub> and &#x1D538; of different TS devices (<xref ref-type="bibr" rid="B14">Park et al., 2016</xref>; <xref ref-type="bibr" rid="B21">Yoo et al., 2017</xref>; <xref ref-type="bibr" rid="B9">Lee et al., 2019b</xref>,<xref ref-type="bibr" rid="B10">2020</xref>), and the simulated <italic>t</italic><sub>on</sub> and <italic>V</italic><sub>th</sub> of the neuron circuit based on the V&#x2013;t Model are indicated in <xref ref-type="fig" rid="F8">Figure 8</xref>. In this study, the <italic>R</italic><sub>series</sub> is given from 10 to 1,000 k&#x03A9;. The value of the integration capacitor in the neuron circuit is adjusted to keep the maximum <italic>V</italic><sub>th</sub> below 1.2 V at <italic>R</italic><sub>series</sub> = 10 k&#x03A9;, and the adopted capacitance corresponding to each TS device is also given. Among IMT, OTS, and Ag-based TS devices, the OTS neuron matched with the lowest capacitance is the most favorable for reducing the neuron circuit area. It is noted that the minimal integration capacitor is limited by the parasitic capacitor of the TS selector itself. As a result, the device area scaling would be necessary to achieve a low enough capacitance value. Moreover, the simulated <italic>t</italic><sub>on</sub> in <xref ref-type="fig" rid="F8">Figure 8A</xref> shows that the OTS neuron is capable of achieving GHz-level spike frequency due to its extremely small &#x03C4;<sub>0</sub> (10 <sup>&#x2013; 21</sup>s), even though its &#x1D538; is larger. Furthermore, in <xref ref-type="fig" rid="F8">Figure 8B</xref>, the OTS selector with the smallest &#x03C4;<sub>0</sub> results in a large <italic>t</italic><sub>on</sub>/&#x03C4;<sub>0</sub>; therefore, the <italic>V</italic><sub>th</sub> is less history dependent. The OTS selector shows promising potential not only in generating high spike frequency but also consuming less area and energy in the neuron circuit.</p>
<table-wrap position="float" id="T1">
<label>TABLE 1</label>
<caption><p>Key parameters extracted from the reported threshold switching (TS) selectors.</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<thead>
<tr>
<td valign="top" align="left">TS selector type</td>
<td valign="top" align="center">&#x1D538; (V&#x22C5;s)</td>
<td valign="top" align="center">&#x03C4; <sub>0</sub> (s)</td>
<td valign="top" align="center">Capacitor (F)<xref ref-type="table-fn" rid="t1fns1">&#x002A;</xref></td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">IMT (<xref ref-type="bibr" rid="B10">Lee et al., 2020</xref>)</td>
<td valign="top" align="center">1.29</td>
<td valign="top" align="center">10<sup>&#x2013;8</sup></td>
<td valign="top" align="center">6 &#x00D7; 10<sup>&#x2013;13</sup></td>
</tr>
<tr>
<td valign="top" align="left">IMT (<xref ref-type="bibr" rid="B14">Park et al., 2016</xref>)</td>
<td valign="top" align="center">1.602</td>
<td valign="top" align="center">10<sup>&#x2013;8</sup></td>
<td valign="top" align="center">7 &#x00D7; 10<sup>&#x2013;13</sup></td>
</tr>
<tr>
<td valign="top" align="left">OTS (<xref ref-type="bibr" rid="B10">Lee et al., 2020</xref>)</td>
<td valign="top" align="center">30.28</td>
<td valign="top" align="center">10<sup>&#x2013;21</sup></td>
<td valign="top" align="center">2 &#x00D7; 10<sup>&#x2013;15</sup></td>
</tr>
<tr>
<td valign="top" align="left">OTS (<xref ref-type="bibr" rid="B9">Lee et al., 2019b</xref>)</td>
<td valign="top" align="center">45.78</td>
<td valign="top" align="center">10<sup>&#x2013;24</sup></td>
<td valign="top" align="center">6 &#x00D7; 10<sup>&#x2013;16</sup></td>
</tr>
<tr>
<td valign="top" align="left">Ag-based (<xref ref-type="bibr" rid="B10">Lee et al., 2020</xref>)</td>
<td valign="top" align="center">3.09</td>
<td valign="top" align="center">10<sup>&#x2013;6</sup></td>
<td valign="top" align="center">3.25 &#x00D7; 10<sup>&#x2013;10</sup></td>
</tr>
<tr>
<td valign="top" align="left">Ag-based (<xref ref-type="bibr" rid="B21">Yoo et al., 2017</xref>)</td>
<td valign="top" align="center">2.92</td>
<td valign="top" align="center">10<sup>&#x2013;6</sup></td>
<td valign="top" align="center">2.95 &#x00D7; 10<sup>&#x2013;10</sup></td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn id="t1fns1"><p><italic>&#x002A;The value of an integrated capacitor in the neuron circuit is adjusted to keep the maximum V<sub>th</sub> below 1.2 V at R<sub>series</sub> = 10 k&#x03A9;.</italic></p></fn>
</table-wrap-foot>
</table-wrap>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption><p>V<bold>&#x2013;</bold>t Model prediction on <bold>(A)</bold> <italic>t</italic><sub>on</sub> and <bold>(B)</bold> <italic>V</italic><sub>th</sub> of the neuron circuit based on insulator&#x2013;metal transition (IMT) (<xref ref-type="bibr" rid="B14">Park et al., 2016</xref>; <xref ref-type="bibr" rid="B10">Lee et al., 2020</xref>), ovonic threshold switching (OTS) (<xref ref-type="bibr" rid="B9">Lee et al., 2019b</xref>,<xref ref-type="bibr" rid="B10">2020</xref>), and Ag-based (<xref ref-type="bibr" rid="B21">Yoo et al., 2017</xref>; <xref ref-type="bibr" rid="B10">Lee et al., 2020</xref>) selectors. <italic>R</italic><sub>series</sub> is assumed to range from 10 k&#x03A9; to 1 M&#x03A9;. The value of an integrated capacitor in each condition is listed in <xref ref-type="table" rid="T1">Table 1</xref>.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-16-868671-g008.tif"/>
</fig>
</sec>
</sec>
<sec id="S4" sec-type="conclusion">
<title>Conclusion</title>
<p>In this study, a V&#x2013;t Model is successfully constructed to simulate the spiking behavior of TS neurons according to the synaptic weight of connected synapses. By considering the history-dependent <italic>V</italic><sub>th</sub> of the TS selector based on the nucleation theory, the proposed V&#x2013;t Model is in good agreement with the measurement results and provides more accurate prediction compared to the conventional RC Model. Moreover, the behavior of TS neurons based on different TS devices, including IMT, OTS, and Ag-based selectors, are simulated and compared using the proposed V&#x2013;t Model. The results show that the OTS selector matched with the lowest capacitance that is the most favorable for reducing the circuit area overhead. Moreover, the OTS selector with the lowest &#x03C4;<sub>0</sub> and <italic>t</italic><sub>on</sub> not only achieves less history-dependent <italic>V</italic><sub>th</sub> but also realizes a high-speed neuron with GHz-level spike frequency. The proposed V&#x2013;t model provides a useful engineering pathway toward the future development of TS neurons.</p>
</sec>
<sec id="S5" sec-type="data-availability">
<title>Data Availability Statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="S6">
<title>Author Contributions</title>
<p>S-MY fabricated the device. S-MY and M-HW performed data analysis. S-MY, I-TW, M-HW, and T-HH contributed to the conception and discussion of the study. S-MY, I-TW, and T-HH drafted manuscript. All authors contributed to the article and approved the submitted version.</p>
</sec>
<sec id="conf1" sec-type="COI-statement">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="pudiscl1" sec-type="disclaimer">
<title>Publisher&#x2019;s Note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
</body>
<back>
<sec id="S7" sec-type="funding-information">
<title>Funding</title>
<p>This work was supported in part by the Ministry of Science and Technology, Taiwan, under Grant Nos. 109-2221-E-009-020-MY3, 110-2221-E-A49-088-MY3, 110-2634-F-009-017, and 110-2622-8-009-018-SB.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chae</surname> <given-names>B.-G.</given-names></name> <name><surname>Seol</surname> <given-names>J.-B.</given-names></name> <name><surname>Song</surname> <given-names>J.-H.</given-names></name> <name><surname>Baek</surname> <given-names>K.</given-names></name> <name><surname>Oh</surname> <given-names>S.-H.</given-names></name> <name><surname>Hwang</surname> <given-names>H.</given-names></name><etal/></person-group> (<year>2017</year>). <article-title>Nanometer-scale phase transformation determins threshold and memory switching mechanism.</article-title> <source><italic>Adv. Mater.</italic></source> <volume>29</volume>:<fpage>1701725</fpage>. <pub-id pub-id-type="doi">10.1002/adma.201701752</pub-id> <pub-id pub-id-type="pmid">28605067</pub-id></citation></ref>
<ref id="B2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname> <given-names>P.-Y.</given-names></name> <name><surname>Seo</surname> <given-names>J.-S.</given-names></name> <name><surname>Cao</surname> <given-names>Y.</given-names></name> <name><surname>Yu</surname> <given-names>S.</given-names></name></person-group> (<year>2016</year>). &#x201C;<article-title>Compact oscillation neuron exploiting metal-insulator-transition for neuromorphic computing</article-title>,&#x201D; in <source><italic>Proceedings of the IEEE/ACM International Conference on Computer-Aided Design (ICCAD)</italic></source>, (<publisher-loc>Austin, TX</publisher-loc>: <publisher-name>IEEE</publisher-name>). <pub-id pub-id-type="doi">10.1145/2966986.2967015</pub-id></citation></ref>
<ref id="B3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Grisafe</surname> <given-names>B.</given-names></name> <name><surname>Jerry</surname> <given-names>M.</given-names></name> <name><surname>Smith</surname> <given-names>J. A.</given-names></name> <name><surname>Datta</surname> <given-names>S.</given-names></name></person-group> (<year>2019</year>). <article-title>Performance enhancement of Ag/HfO2 metal ion threshold switch cross-point selectors.</article-title> <source><italic>IEEE Electron Device Lett.</italic></source> <volume>40</volume> <fpage>1602</fpage>&#x2013;<lpage>1605</lpage>. <pub-id pub-id-type="doi">10.1109/LED.2019.2936104</pub-id></citation></ref>
<ref id="B4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hatem</surname> <given-names>F.</given-names></name> <name><surname>Chai</surname> <given-names>Z.</given-names></name> <name><surname>Zhang</surname> <given-names>W.</given-names></name> <name><surname>Fantini</surname> <given-names>A.</given-names></name> <name><surname>Degraeve</surname> <given-names>R.</given-names></name> <name><surname>Clima</surname> <given-names>S.</given-names></name><etal/></person-group> (<year>2019</year>). &#x201C;<article-title>Endurance improvement of more than five orders in GexSe1-x OTS selectors by using a novel refreshing program scheme</article-title>,&#x201D; in <source><italic>Proceedings of the IEEE International Electron Devices Meeting (IEDM)</italic></source>, (<publisher-loc>San Francisco, CA</publisher-loc>: <publisher-name>IEEE</publisher-name>), <fpage>827</fpage>&#x2013;<lpage>830</lpage>. <pub-id pub-id-type="doi">10.1109/IEDM19573.2019.8993448</pub-id></citation></ref>
<ref id="B5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hua</surname> <given-names>Q.</given-names></name> <name><surname>Wu</surname> <given-names>H.</given-names></name> <name><surname>Gao</surname> <given-names>B.</given-names></name> <name><surname>Zhao</surname> <given-names>M.</given-names></name> <name><surname>Li</surname> <given-names>Y.</given-names></name> <name><surname>Li</surname> <given-names>X.</given-names></name><etal/></person-group> (<year>2019</year>). <article-title>A threshold switching selector based on highly ordered Ag nanodots for X-point memory applications.</article-title> <source><italic>Adv. Sci.</italic></source> <volume>6</volume>:<fpage>1900024</fpage>. <pub-id pub-id-type="doi">10.1002/advs.201900024</pub-id> <pub-id pub-id-type="pmid">31131198</pub-id></citation></ref>
<ref id="B6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ielmini</surname> <given-names>D.</given-names></name> <name><surname>Wong</surname> <given-names>H.-S. P.</given-names></name></person-group> (<year>2018</year>). <article-title>In-memory computing with resistive switching devices.</article-title> <source><italic>Nat. Electron.</italic></source> <volume>1</volume> <fpage>333</fpage>&#x2013;<lpage>343</lpage>. <pub-id pub-id-type="doi">10.1038/s41928-018-0092-2</pub-id></citation></ref>
<ref id="B7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Karpov</surname> <given-names>I. V.</given-names></name> <name><surname>Mitra</surname> <given-names>M.</given-names></name> <name><surname>Kau</surname> <given-names>D.</given-names></name> <name><surname>Spadini</surname> <given-names>G.</given-names></name> <name><surname>Kryukov</surname> <given-names>Y. A.</given-names></name> <name><surname>Karpov</surname> <given-names>V. G.</given-names></name></person-group> (<year>2008</year>). <article-title>Evidence of field induced nucleation in phase change memory.</article-title> <source><italic>Appl. Phys. Lett.</italic></source> <volume>92</volume>:<fpage>173501</fpage>. <pub-id pub-id-type="doi">10.1063/1.2917583</pub-id></citation></ref>
<ref id="B8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lee</surname> <given-names>D.</given-names></name> <name><surname>Kwak</surname> <given-names>M.</given-names></name> <name><surname>Moon</surname> <given-names>K.</given-names></name> <name><surname>Choi</surname> <given-names>W.</given-names></name> <name><surname>Park</surname> <given-names>J.</given-names></name> <name><surname>Yoo</surname> <given-names>J.</given-names></name><etal/></person-group> (<year>2019a</year>). <article-title>Integrate and fire neuron based on various threshold switching devices with scalable device area and ultra-low power operation for neuromorphic system applications.</article-title> <source><italic>Adv. Electron. Mater.</italic></source> <volume>5</volume>:<fpage>1800866</fpage>. <pub-id pub-id-type="doi">10.1002/aelm.201800866</pub-id></citation></ref>
<ref id="B9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lee</surname> <given-names>S.</given-names></name> <name><surname>Yoo</surname> <given-names>J.</given-names></name> <name><surname>Park</surname> <given-names>J.</given-names></name> <name><surname>Hwang</surname> <given-names>H.</given-names></name></person-group> (<year>2019b</year>). <article-title>Field-induced nucleation switching in binary ovonic threshold switches.</article-title> <source><italic>Appl. Phys. Lett.</italic></source> <volume>115</volume>:<fpage>233503</fpage>. <pub-id pub-id-type="doi">10.1063/1.5126913</pub-id></citation></ref>
<ref id="B10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lee</surname> <given-names>S.</given-names></name> <name><surname>Yoo</surname> <given-names>J.</given-names></name> <name><surname>Park</surname> <given-names>J.</given-names></name> <name><surname>Hwang</surname> <given-names>H.</given-names></name></person-group> (<year>2020</year>). <article-title>Understanding of the abrupt resistive transition in different types of threshold switching devices from materials perspective.</article-title> <source><italic>IEEE Trans. Electron Devices</italic></source> <volume>67</volume> <fpage>2878</fpage>&#x2013;<lpage>2883</lpage>. <pub-id pub-id-type="doi">10.1109/TED.2020.2997670IEEE</pub-id></citation></ref>
<ref id="B11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liang</surname> <given-names>F.-X.</given-names></name> <name><surname>Sahu</surname> <given-names>P.</given-names></name> <name><surname>Wu</surname> <given-names>M.-H.</given-names></name> <name><surname>Wei</surname> <given-names>J.-H.</given-names></name> <name><surname>Sheu</surname> <given-names>S.-S.</given-names></name> <name><surname>Hou</surname> <given-names>T.-H.</given-names></name></person-group> (<year>2020</year>). &#x201C;<article-title>Stochastic STT-MRAM spiking neuron circuit</article-title>,&#x201D; in <source><italic>Proceedings of International Symposium on VLSI Technology, Systems and Applications (VLSI-TSA)</italic></source>, (<publisher-loc>Hsinchu</publisher-loc>: <publisher-name>IEEE</publisher-name>), <fpage>151</fpage>&#x2013;<lpage>152</lpage>. <pub-id pub-id-type="doi">10.1109/VLSI-TSA48913.2020.9203701</pub-id></citation></ref>
<ref id="B12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liang</surname> <given-names>F.-X.</given-names></name> <name><surname>Wang</surname> <given-names>I.-T.</given-names></name> <name><surname>Hou</surname> <given-names>T.-H.</given-names></name></person-group> (<year>2021</year>). <article-title>Progress and benchmark of spiking neuron devices and circuits.</article-title> <source><italic>Adv. Intell. Syst.</italic></source> <volume>3</volume>:<fpage>2100007</fpage>. <pub-id pub-id-type="doi">10.1002/aisy.202100007</pub-id></citation></ref>
<ref id="B13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Luo</surname> <given-names>W.-C.</given-names></name> <name><surname>Liu</surname> <given-names>J.-C.</given-names></name> <name><surname>Lin</surname> <given-names>Y.-C.</given-names></name> <name><surname>Lo</surname> <given-names>C.-L.</given-names></name> <name><surname>Huang</surname> <given-names>J.-J.</given-names></name> <name><surname>Lin</surname> <given-names>K.-L.</given-names></name><etal/></person-group> (<year>2013</year>). <article-title>Statistical model and rapid prediction of RRAM SET speed&#x2013;disturb dilemma</article-title> <source><italic>IEEE Trans. Electron Devices</italic></source> <volume>60</volume> <fpage>3760</fpage>&#x2013;<lpage>3766</lpage>. <pub-id pub-id-type="doi">10.1109/TED.2013.2281991</pub-id></citation></ref>
<ref id="B14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Park</surname> <given-names>J.</given-names></name> <name><surname>Cha</surname> <given-names>E.</given-names></name> <name><surname>Karpov</surname> <given-names>I.</given-names></name> <name><surname>Hwang</surname> <given-names>H.</given-names></name></person-group> (<year>2016</year>). <article-title>Dynamics of electroforming and electrically driven insulator-metal transition in NbOx selector.</article-title> <source><italic>Appl. Phys. Lett.</italic></source> <volume>108</volume>:<fpage>232101</fpage>. <pub-id pub-id-type="doi">10.1063/1.4953323</pub-id></citation></ref>
<ref id="B15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Song</surname> <given-names>B.</given-names></name> <name><surname>Xu</surname> <given-names>H.</given-names></name> <name><surname>Liu</surname> <given-names>S.</given-names></name> <name><surname>Liu</surname> <given-names>H.</given-names></name> <name><surname>Li</surname> <given-names>Q.</given-names></name></person-group> (<year>2018</year>). <article-title>Threshold switching behavior of Ag-SiTe-based selector device and annealing effect on its characteristics.</article-title> <source><italic>IEEE J. Electron Devices Soc.</italic></source> <volume>6</volume> <fpage>674</fpage>&#x2013;<lpage>679</lpage>. <pub-id pub-id-type="doi">10.1109/JEDS.2018.2836400</pub-id></citation></ref>
<ref id="B16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tuma</surname> <given-names>T.</given-names></name> <name><surname>Pantazi</surname> <given-names>A.</given-names></name> <name><surname>Le Gallo</surname> <given-names>M.</given-names></name> <name><surname>Sebastian</surname> <given-names>A.</given-names></name> <name><surname>Eleftheriou</surname> <given-names>E.</given-names></name></person-group> (<year>2016</year>). <article-title>Stochastic phase-change neurons.</article-title> <source><italic>Nat. Nanotechnol.</italic></source> <volume>11</volume> <fpage>693</fpage>&#x2013;<lpage>699</lpage>. <pub-id pub-id-type="doi">10.1038/nnano.2016.70</pub-id> <pub-id pub-id-type="pmid">27183057</pub-id></citation></ref>
<ref id="B17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>P.</given-names></name> <name><surname>Khan</surname> <given-names>A. I.</given-names></name> <name><surname>Yu</surname> <given-names>S.</given-names></name></person-group> (<year>2020</year>). <article-title>Cryogenic behavior of NbO2 based threshold switching devices as oscillation neurons.</article-title> <source><italic>Appl. Phys. Lett.</italic></source> <volume>116</volume>:<fpage>162108</fpage>. <pub-id pub-id-type="doi">10.1063/5.0006467</pub-id></citation></ref>
<ref id="B18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Woo</surname> <given-names>J.</given-names></name> <name><surname>Wang</surname> <given-names>P.</given-names></name> <name><surname>Yu</surname> <given-names>S.</given-names></name></person-group> (<year>2019</year>). <article-title>Integrated crossbar array with resistive synapses and oscillation neurons.</article-title> <source><italic>IEEE Electron Device Lett.</italic></source> <volume>40</volume> <fpage>1313</fpage>&#x2013;<lpage>1316</lpage>. <pub-id pub-id-type="doi">10.1109/LED.2019.2921656</pub-id></citation></ref>
<ref id="B19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wu</surname> <given-names>M.-H.</given-names></name> <name><surname>Hong</surname> <given-names>M.-C.</given-names></name> <name><surname>Chang</surname> <given-names>C.-C.</given-names></name> <name><surname>Sahu</surname> <given-names>P.</given-names></name> <name><surname>Wei</surname> <given-names>J.-H.</given-names></name> <name><surname>Lee</surname> <given-names>H.-Y.</given-names></name></person-group> (<year>2019</year>). &#x201C;<article-title>Extremely compact integrate-and-fire STT-MRAM neuron: a pathway toward all-spin artificial deep neural network</article-title>,&#x201D; in <source><italic>Proceedings of IEEE Symp. on VLSI Technol. (VLSI-T), T34-T35</italic></source>, (<publisher-loc>Kyoto</publisher-loc>: <publisher-name>IEEE</publisher-name>). <pub-id pub-id-type="doi">10.23919/VLSIT.2019.8776569</pub-id></citation></ref>
<ref id="B20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wu</surname> <given-names>M.-H.</given-names></name> <name><surname>Huang</surname> <given-names>M.-S.</given-names></name> <name><surname>Zhu</surname> <given-names>Z.</given-names></name> <name><surname>Liang</surname> <given-names>F.-X.</given-names></name> <name><surname>Hong</surname> <given-names>M.-C.</given-names></name> <name><surname>Deng</surname> <given-names>J.</given-names></name></person-group> (<year>2020</year>). &#x201C;<article-title>Compact probabilistic Poisson neuron based on back-hopping oscillation in STT-MRAM for all-spin deep spiking neural network</article-title>,&#x201D; in <source><italic>Proceedings of IEEE Symp. on VLSI Technol. (VLSI-T), JFS4.2.1-JFS4.2.2</italic></source>, (<publisher-loc>Honolulu, HI</publisher-loc>: <publisher-name>IEEE</publisher-name>). <pub-id pub-id-type="doi">10.1109/VLSITechnology18217.2020.9265033</pub-id></citation></ref>
<ref id="B21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yoo</surname> <given-names>J.</given-names></name> <name><surname>Park</surname> <given-names>J.</given-names></name> <name><surname>Song</surname> <given-names>J.</given-names></name> <name><surname>Lim</surname> <given-names>S.</given-names></name> <name><surname>Hwang</surname> <given-names>H.</given-names></name></person-group> (<year>2017</year>). <article-title>Field-induced nucleation in threshold switching characteristics of electrochemical metallization devices.</article-title> <source><italic>Appl. Phys. Lett.</italic></source> <volume>111</volume>:<fpage>163109</fpage>. <pub-id pub-id-type="doi">10.1063/1.4985165</pub-id></citation></ref>
<ref id="B22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>X.</given-names></name> <name><surname>Wu</surname> <given-names>Z.</given-names></name> <name><surname>Lu</surname> <given-names>J.</given-names></name> <name><surname>Wei</surname> <given-names>J.</given-names></name> <name><surname>Lu</surname> <given-names>J.</given-names></name> <name><surname>Zhu</surname> <given-names>J.</given-names></name><etal/></person-group> (<year>2020</year>). &#x201C;<article-title>Fully memristive SNNs with temporal coding for fast and low-power edge computing</article-title>,&#x201D; in <source><italic>Proceedings of IEEE International Electron Devices Meeting (IEDM)</italic></source>, (<publisher-loc>San Francisco, CA</publisher-loc>), <fpage>649</fpage>&#x2013;<lpage>652</lpage>. <pub-id pub-id-type="doi">10.1109/IEDM13553.2020.9371937</pub-id></citation></ref>
</ref-list>
</back>
</article>