1 Introduction
Artificial neurons are crucial for constructing neuromorphic systems. Although CMOS-based circuits have been the most common approach for implementing artificial neurons [
1−
3], they suffer from limitations such as low integration density and high power consumption. The development of alternative scalable neurons, including those based on memristors, is therefore essential for building large-scale neuromorphic computing systems. Biological neurons transmit and process information by spatiotemporally integrating input signals and firing action potentials [
4,
5]. Several models have been developed to mimic these behaviors, including the McCulloch-Pitts (MP) [
6], Hodgkin-Huxley (H-H) [
7], Integrate-and-Fire (IF) [
8], and Leaky Integrate-and-Fire (LIF) models [
9]. The IF model, pioneered by Louis Lapicque in 1907, offers a simplified mathematical framework for neuronal dynamics [
10]. It describes how a neuron integrates synaptic inputs until its membrane potential reaches a threshold, triggering a spike followed by reset [
11]. While computationally efficient, the IF model lacks biological realism by omitting membrane potential leakage. The LIF model addresses this by incorporating a leakage term that mimics the natural decay of membrane potential, providing a more physiologically plausible representation [
12]. The LIF model is particularly valuable for spiking neural networks due to its ability to simulate leakage, spatiotemporal integration, and firing functions inherent to biological neurons. In biological systems, these short-term dynamics arise from ion concentration differences across the cell membrane, which generate capacitive-like effects. Interestingly, similar capacitive effects occur in many memristive devices [
13], enabling the implementation of LIF behaviors with suitably designed circuits. Thus, memory devices exhibiting capacitive effects are promising for constructing neurons with leaky integration and firing capabilities.
Both memristors and capacitors possess similar two-terminal structures, with a dielectric layer sandwiched between two metal electrodes. Therefore, memristors can exhibit capacitive characteristics during the operating stage under certain conditions. Up to now, the influence of the capacitive-coupling memory effect has been observed in many memristors [
14−
16]. In addition, it has been found that numerous ferroelectric materials, such as Ba
0.6Sr
0.4TiO
3 (BST) [
17], BaTiO
3 (BTO) [
18], Pb(Zr
0.4Ti
0.6)O
3 (PZT) [
19] and so on, can be used not only as parallel-plate capacitors but also as the dielectric layers of memristors. However, due to the high dielectric constant of ferroelectric materials, their tunability is limited. Introducing another substance to form a composite structure with ferroelectrics is a promising method for developing new materials with special coupling properties suitable for exploitation. Among the above materials, the Ba
0.6Sr
0.4TiO
3 (BST) compound has been widely used in tunable microwave devices and dynamic memory devices due to its material properties of high dielectric constant and low dielectric loss [
17]. Additionally, CeO
2 has been widely applied in memristor devices because of its high oxygen ion conductivity, valence states of copper ions, inherent resistive switching (RS) characteristics, large on/off ratio, remarkable plasticity and compatibility with current CMOS technology [
20−
22]. Recent studies have further highlighted the potential of advanced materials in optimizing memristor performance, providing valuable references for the design of composite structures [
23−
25]. In oxide-based memristors, the capacitive effect is mainly attributed to the migration of oxygen ions. Hence, the composite thin-film devices combining BST and CeO
2 are expected to achieve the characteristics of capacitive-coupling memory.
Notably, existing studies have explored memristor-based implementations of neuromorphic systems with functionalities relevant to LIF neurons, such as the brain-inspired in-memory computing system leveraging memristive circuits for neuronal communication [
26], the multimodal local-global neuromorphic computing system (MLG-NCS) for affective video content analysis [
27], and the time-frequency hybrid neuromorphic architecture for battery state-of-health estimation [
28]. Unlike these works that focus on system-level integration for specific tasks (neuronal communication, affective analysis, battery estimation), this study specifically emphasizes optimizing memristive device-level capacitive effects to enhance the biological realism and scalability of LIF neuron implementations, addressing the trade-off between computational efficiency and physiological plausibility that is less highlighted in the aforementioned task-oriented systems.
Herein, this study proposes and fabricates a Pt/La0.5Sr0.5CoO3(LSCO)/BST:CeO2/LSCO device with capacitive-coupled memristive effects. High-quality BST:CeO2 composite films were deposited on single-crystalline SrTiO3 (STO) substrates through a series of techniques such as magnetron sputtering and pulsed laser deposition. Initially, the structural and elemental characteristics of the device were thoroughly investigated. Under electrical modulation, the device’s transition from capacitive behavior to pure memory behavior was revealed by varying the voltage magnitude. Furthermore, the device’s capacity to mimic a variety of biological synaptic functions in the pure memristor state was explored. Through integrating the device into a custom-designed circuit, the LIF function was simulated. Especially, based on spiking neural network (SNN), this device is expected to be applied to monitor heartbeats and could early diagnose cardiovascular diseases. Thereby, this research lays a solid foundation for the development of advanced neuromorphic computing systems with augmented performance and functionality.
2 Result and discussion
Figure 1(a) presents the structure of the Pt/LSCO/BST:CeO2/LSCO device. The cross-sectional TEM image of the device is shown in Fig. S1. Meanwhile, the EDS elemental mapping results have effectively verified the elemental distribution within the Pt/LSCO/BST:CeO2/LSCO device, as depicted in Fig. S2. Figure S2(b) shows the element distribution of Pt (top electrode). In Figs. S2(c)−(i), the elements La, Ti, Co, Ce, Sr, O, and Ba were uniformly distributed in the dielectric layer, indicating that the chemical composition of the Pt/LSCO/BST:CeO2/LSCO was homogeneous. The X-ray diffraction (XRD) spectra (Fig. 1(b)) shows the BST thin films display a (001) orientation with respect to the STO substrate, indicating the successful formation of pure perovskite thin films. The Φ-scanning of the BST (001) planes, in view of the tetragonal BST crystal structure, the four peaks are uniformly distributed, which reveals the fourfold symmetry of the epitaxial BST film, as shown in Fig. 1(c). To further prove this conclusion, the highresolution transmission electron microscopy (HRTEM) images of the initial Pt/LSCO/BST:CeO2/LSCO device are shown in Figs. 1(d)−(f), in which clear lattice fringes can be observed. The crystal plane indices, such as (001), (100), and (110), indicate a consistent crystal orientation relationship among the BST:CeO2, LSCO, and STO layers at the interfaces. As revealed in Fig. 1(e), the atom positions near the interface can be clearly resolved, indicating that the BST:CeO2 film grew perfectly and coherently on the LSCO (bottom electrode) lattice because they have the same type of crystal structure and almost same lattice constants. Moreover, the interfaces between the material layers are highly flat without any noticeable lattice distortion or defects. The flatness of the interfaces is beneficial for reducing interface scattering, thereby enhancing the electrical performance of the device. Additionally, the lattice fringes in the images are extremely sharp, suggesting that the crystal structures of the materials are nearly perfect with minimal crystal defects. These observations collectively confirm the superior epitaxial properties of the device, indicating its potential for excellent performance in high-performance electronic device applications.
Subsequently, in order to investigate the electrical properties of the device,
I−V measurements were carried out on the Pt/LSCO/BST:CeO
2/LSCO device. A voltage was applied to the tested devices and scans were executed between
Vmax (maximum applied voltage) at a constant scanning rate of 0.01 V/s. The acquired
I−V curves with
Vmax values of 0.5, 1.5, and 3.0 V are illustrated in Figs. 2(a)−(c) (linear scale) and Figs. 2(d)−(f) (logarithmic scale), respectively. At
Vmax of 0.5 V, Figs. 2(a) and 2(d) show that the
I−V curve has a non-zero crossing without an intersection point. It is well known that such a non-zero-crossing
I−V curve indicates the existence of a capacitive effect in the memristor [
29−
33]. Upon increasing
Vmax to 1.5 V, a crossing behavior manifested within the
I−V curve was discernible [Fig. 2(b)], with the intersection point situated in the first quadrant [Fig. 2(e)]. Nevertheless, this
I−V curve remained non-zero-crossing. When the test voltage was further increased to
Vmax = 3.0 V, a zero-crossing
I−V hysteresis curve indicating the standard memory behavior of the device was observed [Figs. 2(c) and (f)]. These
I−V characteristics are reproducible and can be controlled by selecting
Vmax. The above investigations reveal that by augmenting
Vmax, the transformation of the device from capacitive behavior to capacitively-coupled memory behavior and ultimately to pure memory behavior can be accomplished. Based on the above analysis, further research on the pure memristive performance of the device was conducted. The device exhibits stable memristive behavior with over 100 cycles of operation without significant degradation, and the high and low resistance states are well-defined and concentrated, as shown in Fig. S3. Furthermore, its intrinsic memristive properties enable the emulation of biological synapse functions (Fig. S4). Figure 2(g) illustrates the carrier transport mechanism in the functional layer of the device, showing the transitions between the HRS and the LRS in both the low voltage district (LVD) and the high voltage district (HVD) of the device. It can be seen that when a positive voltage is applied to the Pt electrode of the device, the positively charged oxygen vacancies (V
o2+) in the BST:CeO
2 functional layer will move towards the bottom electrode along the direction of the electric field, and shield the negatively charged region at the bottom interface. The positive and negative ions (V
o2+ and O
2-) moving in the electric field carry displacement current, resulting in a capacitive effect in the functional layer at low voltages. The displacement current presents a non-zero-crossing
I−V curve [
34]. Moreover, as the voltage increases, the ion concentration increases correspondingly, which further alters the memory behavior caused by the internal electromotive force [
35], so the capacitive-coupled memory effect will be exhibited in the device. Under the application of a high voltage, local Joule heating is generated in the device due to the high current, improving the mobility of oxygen vacancies, and the high electric field can effectively accelerate the migration of electrons and ions [
36−
40]. After reaching a certain threshold, conductive filaments (oxygen vacancies) will eventually form within the BST:CeO
2 composite oxide layer. Once the conductive filaments are formed, they provide a path for the conductive current to move between the two electrodes, surpassing any displacement current. At this time, the capacitive effect is largely suppressed completely, and a pure memory effect is exhibited. After the voltage polarity is reversed, due to the Coulomb repulsion effect, the oxygen vacancies are pushed back to the top electrode, causing the conductive filaments to partially dissolve and form a relatively large gap in the functional layer. The carriers of the conduction current are expected to pass through this gap [
41,
42]. Therefore, after such a reset process, the device returns to the HRS, and this process under high voltage is similar to that reported in standard memristive devices.
Based on the capacitive-coupled memristive characteristic of the device, it can be incorporated into a circuit to demonstrate an artificial neuron with leaky, integrate and fire functions. The physical diagram is shown in Fig. S5. In biological neurons, the difference of ion concentration inside and outside cells leads to leaky integrate-and-fire potential (LGP) with short-term kinetics. The difference of ion concentration can induce the cell with concentration change. This effect has been widely used in many memristor devices [
30]. Ion migration in the switching process of memristor will produce an internal electric field, which can realize short-term dynamic changes and can be used to simulate LGP in biological neurons. More interestingly, the capacitive effect can simulate the hyper potential of biological spikes. As shown in Fig. 3(a), a simple circuit can be designed by using the capacitive-coupled memristive devices to realize LIF function, so as to simulate the process of neuronal transmission.
Therefore, the memristive device with capacitive effect is very important to construct neurons with LIF function. The circuit of the artificial neuron with LIF functions is shown in Fig. 3(b). The circuit contains two devices M1 and M2 of Pt/LSCO/BST:CeO2/LSCO structure. An input signal is applied on M1 to simulate LGP. When the LGP exceeds the threshold, the comparator generates a rising edge and emits a pulse to stimulate M2. The voltage on R3 (1 kΩ) is emitted as a spike similar to biological excitation. In a biological neuronal system, each neuron transmits neural signals through multiple synapses. When a number of dendrites of a neuron receive signals simultaneously, these signals will be summed together in a process called spatial integration. Temporal integration is the sum of the signals arriving at different frequencies of the neuron. Therefore, by controlling the amplitude and frequency of the input pulses applied to M1, the conductance of M1 can be changed to simulate the spatial and temporal integration of the LGP. A spike output will be generated when LGP is above the threshold, and the conductance of M1 will slowly recover to the high resistance state when the input stimulus is small.
To investigate the response of the neuron device to input signals, we first applied input pulses with a fixed pulse width of 0.1 ms and a fixed pulse interval of 0.5 ms, while adjusting the pulse amplitudes to four specific values: −0.5 V, −0.8 V, −1 V, −1.3 V, −1.5 V, and −1.8 V, as illustrated in Fig. 4(a). When the sum of these input signals (with the aforementioned pulse parameters) arriving simultaneously at the neuron reaches a sufficiently high level, the LGP will rise above the preset threshold (250 mV, consistent with the threshold setting in the LIF circuit described earlier), at which point a spike output will be generated. Figure 4(b) simulates temporal integration by varying the frequency of the input pulse signal, with the input pulse amplitude fixed at −1.5 V and frequencies set to 1 kHz, 1.25 kHz, 1.67 kHz, 2.5 kHz, and 5 kHz respectively. When the frequency of the input signal is low (1 kHz) the neuron will not produce output spikes. When the frequency of the input signal is high enough, output spikes will be generated and the frequency of the spike output will increase as the frequency of the input signal increases. To further verify the firing conditions of the device, six continuous voltage pulses with different frequency ranges (from 1 kHz to 3 kHz) and different amplitude ranges (from −0.5 V to −2 V) were applied to the neuronal device in various combinations, and the discharge conditions were recorded [Figs. 4(b)−(e)]. As depicted in Fig. 4(d), according to the differences in the applied conditions, it can be divided into two regions: the non-firing region (pink) and the firing region (green). When the frequency and amplitude are relatively small, the generated LIF potential (LGP) cannot reach Vth [Fig. 4(c)], which corresponds to the non-firing condition. As the frequency and amplitude increase, the LIF neuron gradually transitions to the firing region and finally reaches the firing state [Fig. 4(e)]. The statistical results are instructive for the selection of experimental parameters in subsequent work.
Based on the properties of the LIF neuron circuit validated by the above tests, a SNN for arrhythmia classification was constructed to verify the effectiveness of this LIF neuron circuit in biomedical signal processing tasks. The network structure is shown in Fig. 5(a), employing a two-layer fully connected architecture. The input layer node count matches the feature dimension of the ECG signals, set to 235, corresponding to the number of features in a single ECG signal from the MIT-BIH arrhythmia database. The output layer contains four neurons, each corresponding to one of the four rhythm categories. In the data processing stage, a rate-coding strategy was used to convert continuous ECG signals into spike trains. Specifically, each sample was divided into 10 time steps along the temporal dimension, and a binarized spike sequence was generated by comparing randomly generated thresholds with the signal amplitude. This encoding method effectively preserves the signal’s amplitude characteristics and provides SNNs with impulse signals compatible with the input characteristics of biological neurons. The experimental data were collected from the MIT-BIH arrhythmia database, which covers diverse arrhythmia types, such as normal sinus rhythm, premature ventricular beats, and premature atrial beats. Four specific categories were selected for training and classification: Normal beat (N), Supra ventricular escape beat (SVEB), Ventricular escape beat (E), and Fusion of ventricular and normal beat (F). After training, the SNN achieved 91% classification accuracy on the test set. The test accuracy over 50 training epochs is shown in Fig. 5(b). The confusion matrix, visualizing the model’s classification performance for different arrhythmia categories [Fig. 5(c)], demonstrates effective discrimination across all four ECG data classes. Figure 5(d) displays typical ECG signal samples from each category, along with the outputs of their corresponding output neurons.
These results indicate that the network possesses strong discriminative ability for major arrhythmia types, extracting key features in ECG signals through its spiking patterns. This outcome validates the effectiveness of the LIF neuronal circuit based on BST:CeO2 memristors in real-world biomedical signal processing tasks, providing a theoretical and experimental foundation for applying memristor-based neuromorphic computing systems in medical diagnosis.
3 Experimental
The device with a structure of Pt/LSCO/BST:CeO2/LSCO was fabricated on a single-crystalline STO substrate. The specific steps are as follows: the single-crystalline STO was ultrasonically cleaned in acetone and alcohol for 10 minutes respectively to remove surface impurities, and then the cleaned substrate was placed into a magnetron sputtering chamber. Firstly, the chamber was evacuated until reaching a vacuum of 2 × 10−4 Pa, and then the total gas flow rate was configured to 100 sccm with the gas flow ratio of argon to oxygen set at 3:1. After that, the gate valve was precisely adjusted to maintain the pressure at 1.4 Pa, the radio frequency power was set to 50 W and the deposition temperature was set to 550 °C. After pre-sputtering for 20 minutes, formal sputtering was carried out for 30 minutes. The annealing temperature was 600 °C, and the thickness of the grown LSCO thin film was 50 nm. The BST thin film was fabricated via pulsed laser deposition. Initially, the chamber was evacuated to attain a background vacuum of 2 × 10−4 Pa. Subsequently, the gas flow ratio of argon to oxygen was adjusted to 75:25. Concurrently, the deposition pressure within the chamber was stabilized at 5 Pa and the deposition temperature was set at 885 °C. The laser energy was maintained within the range of 314 to 316 mJ, while simultaneously the laser frequency was fixed at 1 Hz. Meanwhile, set the frequency of magnetron sputtering at 50 W for doping CeO2. The LSCO thin film with a thickness of 20 nm was repeatedly grown on the BST thin film. Finally, the Pt circular top electrode with the thickness of 50 nm (diameter : 100 μm) was fabricated by means of magnetron sputtering with a mask.
The cross-sectional image was characterized by HRTEM. The crystal structure of BST was characterized by XRD. The DC electrical performance of the device is tested by Keithley 4200SCS source measurement unit. Pulse measurements for the Pt/LSCO/BST:CeO2/LSCO devices were performed with the Rigol DG5072 function/arbitrary waveform generator and the Rigol DS4022 oscilloscope.
4 Conclusions
Artificial neurons with functions like LIF, and spike output are of utmost importance for highly efficient brain-inspired computing. In this study, a functional layer was formed by integrating Ba0.6Sr0.4TiO3 and CeO2 within the memory device, which verified its capacitive effect and non-volatile storage capacity. When the voltage was ramped up from zero, the device initially entered a capacitively-coupled memory state at low voltages and subsequently transitioned to a pure memory state at higher voltages. This work not only simulated the biological synaptic plasticity function of the device in the pure memristor state, but also constructed a neuron model that concisely and effectively replicated the LIF functions of biological neurons based on the unique capacitive coupling memristive behavior of the device, thereby demonstrating the reliability of the simulation of neuron transmission. Finally, a SNN based on this device was constructed to recognize the ECG dataset, achieving an identification accuracy of 91%. The above research collectively demonstrated the potential of this capacitive coupling memristive device in neural network applications. Future efforts will focus on enhancing the device’s operational stability across varying temperatures and further elucidating the underlying ion migration and capacitive-coupled switching mechanisms. At the circuit level, developing more bio-plausible and energy-efficient neuron models is essential for enabling large-scale, low-power spiking neural networks. On the system level, advanced learning algorithms and robustness to low SNR signals will be prioritized to facilitate real-world biomedical applications, such as real-time health monitoring and edge-computing diagnostics.