
==== Front
Nat Electron
Nat Electron
Nature Electronics
2520-1131
Nature Publishing Group UK London

38681725
1137
10.1038/s41928-024-01137-9
Article
Nanofluidic logic with mechano–ionic memristive switches
http://orcid.org/0000-0001-5923-4918
Emmerich Theo theo.emmerich@epfl.ch

1
http://orcid.org/0000-0003-2372-2340
Teng Yunfei 12
http://orcid.org/0000-0001-8069-2411
Ronceray Nathan nathan.ronceray@epfl.ch

1
http://orcid.org/0000-0001-8975-413X
Lopriore Edoardo 3
Chiesa Riccardo 3
Chernev Andrey 1
http://orcid.org/0000-0002-6855-5606
Artemov Vasily 1
http://orcid.org/0000-0001-9416-189X
Di Ventra Massimiliano 4
http://orcid.org/0000-0002-3426-7702
Kis Andras 3
http://orcid.org/0000-0001-8194-2785
Radenovic Aleksandra aleksandra.radenovic@epfl.ch

12
1 https://ror.org/02s376052 grid.5333.6 0000 0001 2183 9049 Laboratory of Nanoscale Biology, Institute of Bioengineering, Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne, Switzerland
2 https://ror.org/02s376052 grid.5333.6 0000 0001 2183 9049 NCCR Bio-Inspired Materials, Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne, Switzerland
3 https://ror.org/02s376052 grid.5333.6 0000 0001 2183 9049 Laboratory of Nanoscale Electronics and Structures, Institute of Electrical and Microengineering & Institute of Materials Science and Engineering, Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne, Switzerland
4 grid.266100.3 0000 0001 2107 4242 Department of Physics, University of California, San Diego, La Jolla, CA USA
19 3 2024
19 3 2024
2024
7 4 271278
9 7 2023
21 2 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Neuromorphic systems are typically based on nanoscale electronic devices, but nature relies on ions for energy-efficient information processing. Nanofluidic memristive devices could thus potentially be used to construct electrolytic computers that mimic the brain down to its basic principles of operation. Here we report a nanofluidic device that is designed for circuit-scale in-memory processing. The device, which is fabricated using a scalable process, combines single-digit nanometric confinement and large entrance asymmetry and operates on the second timescale with a conductance ratio in the range of 9 to 60. In operando optical microscopy shows that the memory capabilities are due to the reversible formation of liquid blisters that modulate the conductance of the device. We use these mechano–ionic memristive switches to assemble logic circuits composed of two interactive devices and an ohmic resistor.

Nanofluidic devices with a large entrance asymmetry can function as memristive switches—operating on the second timescale and with a conductance ratio of up to sixty—and can be assembled into basic logic circuits.

Subject terms

Nanopores
Fluidics
https://doi.org/10.13039/100010661 EC | Horizon 2020 Framework Programme (EU Framework Programme for Research and Innovation H2020) 101020445 894369 Kis Andras Radenovic Aleksandra https://doi.org/10.13039/501100001711 Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (Swiss National Science Foundation) TMPFP-2217134 157739 Emmerich Theo Kis Andras https://doi.org/10.13039/100010665 EC | EU Framework Programme for Research and Innovation H2020 | H2020 Priority Excellent Science | H2020 Marie Sklodowska-Curie Actions (H2020 Excellent Science - Marie Sklodowska-Curie Actions) 881603 Kis Andras issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcMain

A fundamental difference between artificial computers and biological brains is the nature of their information carriers. Computers rely on electrons and holes, while brains employ a range of different ions to process data1, using nanoscale synaptic ion channels to perform information processing and storage at energy costs orders of magnitude lower than solid-state digital circuits2,3. Understanding and harnessing ionic transport under nanometric confinement is a goal of nanofluidics4. This has been achieved in different geometries, including zero-dimensional nanopores5,6, one-dimensional nanotubes7–9 and two-dimensional (2D) slits10–13. Such platforms have been used to emulate biological ionic transport8,9,14–25 and could potentially be used to perform computational operations. Recent work has shown, in particular, that nanofluidic channels filled with aqueous electrolytes can store information by exhibiting a memristive effect2,26–31.

A memristive device, which can also be called a memristor32, is a passive two-terminal electrical component with a programmable conductivity that depends on its previous history of operation. Memristors can act as the artificial equivalent of biological synapses due to their ability to store information as a conductance value33, thus enabling data storage and processing with a single device. They can be used in electronic neural networks to adjust the connection strengths at nodes between crossbars, acting as the basic unit for brain-inspired, or neuromorphic, computing34.

Unlike previous fluidic memristors35–37, micropipettes coated with polyimidazolium brushes and bidimensional slits can exhibit memory effects in simple electrolytes within the electrochemical window of water, specifically at voltages below 1.23 V (refs. 26,27). For both devices, this behaviour is understood using phenomenological modelling by considering the combined effects of asymmetric entrances and slow interfacial diffusion, resulting in cycles of accumulation and depletion of ions inside the channel. Such a phenomenology recalls the induction of long-term potentiation through local changes in calcium ion concentration occurring in biological synapses38,39. Yet, direct experimental confirmations regarding the mechanisms of nanoscale fluidic memristors are still lacking, which restricts efforts to improve the performances of these devices.

The large volume occupied by micropipettes and the complex fabrication procedure of 2D slits, as well as the slow speed or small hysteresis displayed by these devices, also hinders their use as synaptic components on the circuit scale. To enable practical computing applications—such as logic operations, pattern recognition and image processing—it is essential to connect multiple memristive devices28,40–44. In this Article, we report devices—termed ‘highly asymmetric channels’ (HACs)—that are designed for circuit-scale neuromorphic nanofluidics. By combining electrokinetic measurements with optical observations, we show that the memory effect in the HACs arises from a combined mechano–ionic effect. With the HACs, we implement a nanofluidic logic operation by connecting two fluidic cells.

Highly asymmetric channels as nanofluidic memristors

To fabricate HACs (Fig. 1a and Supplementary Information Section 1), we start with silicon nitride (SiN) windows presenting a single circular aperture in their centre (step 0). A thin discontinuous layer of palladium is then evaporated (step 1). Finally, a graphite crystal (with thickness between 50 and 150 nm) is deposited above the aperture (step 2). Since step 0 and step 1 are achieved at wafer scale, this current production approach is only limited by step 2 achieved at the single-chip scale. This step is, however, a simple dry-transfer of a large pristine van der Waals crystal not requiring high positioning accuracy. The main fabrication improvement lies in eliminating the intricate etching process necessary for creating 2D channels10,13. With further developments, the complete nanofabrication of HACs could be realized at wafer scale. Nevertheless, the current process already enables us to produce HACs in batches of several tens of units by avoiding delicate fabrication steps at the single-device level. This scalability is required for building nanofluidic circuits, where several devices have to function simultaneously. In this study, we present results for the 32 devices referenced in Supplementary Table 1 with corresponding experiments.Fig. 1 Highly asymmetric channels.

a, Nanofabrication process flow. Step 0: the starting point is silicon chips (in green) with an SiN membrane (in yellow) of around 20 × 20 microns, a thickness of 20 nm and an aperture with a diameter of approximately 100 nm. Step 1: palladium islands (in grey) are placed by evaporative deposition. As a result of careful tuning of the deposition parameters, palladium islands form, leaving space for ions to flow around them. Step 2: dry transfer emplaces a graphite crystal (in blue) with a lateral dimension of 20–50 µm. b, TEM images. Top, the top view of the SiN aperture after step 1; bottom, a cross-sectional view of a completed HAC device. c, Device and setup. Left, sketch of a finished device, with colours as in a. Black lines indicate typical field lines and fluidic paths. For clarity, the graphite crystal, aperture and islands are not at scale. Right, experimental setup for nanofluidic measurements. d, I–V characteristics (Device A, 50 mHz, 1 M KCl). The I–V curve is composed of loops self-crossing at the origin, the signature of a memristive effect. Arrows indicate the direction of the sweep. The dots show individual data points highlighting the abrupt switching and quasi-discrete conductance states. Inset, G–V curve extracted from d. The conductance is defined as the instantaneous current/voltage ratio, G(t) = I(t)/V(t). Here, the conductance is between 6 and 120 nS, yielding a conductance ratio of 20.

Source data

Characterization using transmission electron microscopy (TEM) provides information regarding the geometry of HACs (Fig. 1b and Supplementary Figs. 2 and 3). The Pd islands have a characteristic lateral dimension of 5–10 nm and a spacing and a height of a few nanometres. Ions and water molecules experience single-digit nanometric confinement by flowing around palladium islands, converging to or diverging from the SiN aperture depending on the sign of the applied potential. HACs present a large asymmetry between their two entrances as highlighted in Fig. 1c. The inner entrance, through the membrane aperture, has a lower area Ain = πDh where h is the island’s height and D is the aperture diameter. The outer entrance through the crystal edge has a larger area Aout ≈ πLh, where L is the crystal characteristic dimension of between 20 and 50 μm. Thus, the ratio of entrance areas LD is on the order of several hundred, depending on the size of the top layer crystal. We expect this ratio will promote ionic accumulation when counter-ions enter the channels through the crystal edge and depletion when they enter through the aperture, according to the theoretical framework presented by ref. 27.

To perform nanofluidic measurements, HACs are placed into a fluidic cell separating two reservoirs filled with potassium chloride aqueous solution at pH 5.5 (Fig. 1c). Upon the application of a sinusoidal potential, HACs exhibit a clear bipolar memristive signature27 (the conductance increases at positive voltage and decreases at negative ones with memory retention at low voltage as shown in Fig. 1d), favourable for computing operations with programming voltage pulses. At 1 M KCl, they operate at frequencies in the 30–300 mHz range, with a conductance ratio between 9 and 60 depending on the device (Fig. 1d and Supplementary Table 1). The memory effect in HACs is therefore typically one order of magnitude stronger and two orders of magnitude faster compared to slits exhibiting similar bipolar dynamics, while presenting a larger hysteresis than micropipettes26,27. Our improved performance allows for complete setting and resetting HACs with 2 seconds of voltage pulses (Supplementary Fig. 6). We additionally verified that HACs can be set to intermediate levels by applying short voltage pulses, showing potential for implementing neuromorphic functions such as spike-timing-dependent plasticity (Supplementary Fig. 7). Besides exhibiting a combination of fast speed and large conductance ratio, the response of HACs also qualitatively differs from previous devices by exhibiting delayed switching45. When the applied voltage becomes positive, the response is initially linear until a given threshold is reached where the conductance dramatically increases (Fig. 1d). This abrupt transition from the OFF state to the ON state justifies referring to this phenomenon as ‘switching’. Such dynamics are reminiscent of solid-state electrochemical metallization memory cells46 and have not been reported in nanofluidics. The switching threshold (ON/OFF behaviour) enables reading the memory state without disturbing the programmed state and is beneficial for several neuromorphic computing applications, such as bio-realistic Hebbian learning or conditional logic42,47. These initial findings demonstrate the successful fabrication of scalable nanofluidic switches that exhibit excellent performance and operate effectively using simple monovalent salt solutions while staying within the electrochemical window of water.

The nature of the threshold behaviour is studied by tuning the applied sinusoidal bias V = A sin(ωt). The voltage at which the switching occurs depends on the bias amplitude (A) (Fig. 2a and Supplementary Fig. 8 with two additional devices). Consequently, switching is not triggered at a specific voltage threshold value, regardless of the chosen voltage waveform (Supplementary Fig. 9). Current traces recorded at different frequencies highlight the switching nature of ionic transport in HACs (Fig. 2b and Supplementary Fig. 11), occurring at timescale faster than 100 ms. The delay time τ before switching increases with decreasing frequency of the applied bias signal. Regardless of the applied frequency, the switching occurs when a given amount of charge Q has flowed out of the device (Fig. 2b–c). We recover such conservation of the charge threshold with three additional devices, as shown in Supplementary Fig. 11. The complete I–V cycles at various frequencies, leading to linearity at higher frequencies as expected for a memristor, are displayed in Supplementary Fig. 10 for four devices. In Supplementary Fig. 10, we also show frequency sweeps down to 1 mHz, showing a diode-like behaviour corresponding to the low-frequency limit of a charge-threshold bipolar memristor.Fig. 2 Nature of the switching threshold (Device A-1M KCl).

a, I–V characteristics for different applied sinusoidal potential amplitudes (A). The applied bias frequency is 50 mHz. b, Device switching at different frequencies of the applied positive sinusoidal bias voltage with an amplitude of 1 V. The dots are individual data points. The shaded area represents the charge threshold Q=∫0τI(t)dt, where I(t) is the measured current. This quantity is conserved across frequencies. c, Dependence of the conductance on the cumulative amount of charge, q(t)=∫0tI(t)dt flowing out of the device, extracted from b. The conductance is the instantaneous current/voltage ratio. The conductance abruptly increases when the charge threshold is reached. Here, the charge threshold is approximately equal to 50 nC. Inset, frequency dependence of the amount of charge Q flowing out of the device before reaching the threshold.

Source data

Memory has a mechano–ionic origin

To understand the origin of such memristive dynamics, we observe HACs in operando through wide-field light reflection imaging, as illustrated in Fig. 3a. Such direct observation of nanofluidic processes and, in particular, memory, has been identified as a major challenge to be reached for the field28,43. Briefly, the device chip is mounted in a custom fluidic cell similar to Fig. 1c, but with a top reservoir large enough to accommodate a water-dipping objective. The reflected light intensity from the multi-layer structure of the device is set by thin film interference, which was previously shown to enable monitoring deformations of 2D devices at the 10 nm scale48. See Supplementary Information Section 4.1 for details regarding in our in operando optical setup.Fig. 3 In operando microscopy reveals the memory mechanism.

a, Setup schematic for correlative electrokinetic and optical measurements. The reflected signal from the chip backside is collected with a water-dipping objective while measuring ionic current. b, I–V characteristic (Device 2, 1 V, 30 mHz) with optical images of the SiNx window at different time points marked by white dots and numbers. White arrows illustrate blister motion. c, Optical signal from the pore region denoted by the yellow square in the images in b and its correlation with the device conductance during the I–V acquisition. d, Threshold mechanism observed in b and c: a liquid-filled blister forms at positive bias (branch 0 to 1 in b) between the graphite and Pd–SiNx walls, and the threshold occurs when the blister front crosses the nanopore (at 1 in b).

Source data

Our combined electrokinetic and optical results are presented in Fig. 3b. Concomitantly with the I–V curve acquisition, we observe the emergence of interference fringes at positive voltage and their disappearance at negative ones. These patterns arise from the formation of a (sub-)micrometric liquid blister between the Pd–SiN surface and the graphitic surface. Since monochromatic light is used for imaging (λ = 635 nm), a contrast inversion corresponds to a blister height difference of half a wavelength. The blister modulates ion conduction dramatically by uncovering the SiN aperture, which suppresses the HAC’s resistance at its bottleneck, explaining the larger conductance at positive voltage and overall bistability. As the thickness of the graphite crystal is substantially larger than the thickness of the SiN, we infer that most of the deformation arises from the SiN membrane.

For the devices in Supplementary Videos 1 and 2, the switching occurs when the blister edge moves across the SiN aperture (Fig. 3b–d). In this mechanism, the conductance state of the device correlates perfectly with the reflected intensity signal collected from the pore region (integrated in a 3 x 3 µm window around the pore), as shown in Fig. 3c. As observed in Fig. 3b and sketched in Fig. 3d, when a small off-centred blister is present, the conductance is low (OFF state, at 0), but after the threshold corresponding to the enlarged blister crossing the pore (at 1) the conductance reaches a higher value (ON state, at 2). Within the simplified radial model introduced in Supplementary Section 5, the conductance ratio can be expressed as: GON/GOFF=ln(rout/rin)/ln(rout/rblister). With rin = 50 nm, rout = 50 μm and rblister = 30 μm, this formula predicts the right order of magnitude of observed conductance ratios.

In other cases, the blister forms directly above the hole and we thus attribute the switching to a sudden increase of blister volume occurring outside the SiN window (Supplementary Videos 3 and 4). This dichotomy points to the importance of optimizing HACs to control the blister motion and thereby the performance of the memristor. This will enable the building of neuromorphic nanofluidic circuits with a large number of components. All videos containing matching optical signals and I–V plots can be found in Supplementary Information.

In operando imaging shows that blisters are expanding at positive voltages and shrinking at negative voltages. We identified two plausible mechanisms to explain this phenomenology. First, the focusing of excess counter-ions can result in coulombic repulsion, overcoming the palladium–graphite adhesion near the edge of the blister. Second, as the blister morphological changes align with the expected electro-osmotic flow for negatively charged surfaces, the electro-osmotic flow could modulate the size of small-scale pre-existing blisters, by creating an effective electro-osmotic pressure that can overcome the adhesion between palladium spacers and graphite. In both cases, there is a competition between surface charge (inducing coulombic repulsion or electro-osmotic pressure) and palladium–graphite adhesion. We show in Supplementary Information Sections 5.3, 5.4 and 5.5 that electrostatic repulsion induced by charge focusing, as well as the repulsion induced by electro-osmotic pressure, could both overcome palladium–graphite adhesion. Residual strain in the SiN membrane may also play a role in the blister formation.

We verify our key assumptions with a set of control experiments. Using programming pulses, we notice HACs remain in a high conductance state after setting, similar to synaptic long-term potentiation as the blister requires a negative voltage to vanish (Supplementary Figs. 12 and 21), confirming that conductivity is directly correlated to the blister state. Pressure-driven streaming currents, as well as osmotic currents induced by the concentration gradient, show that HACs exhibit K+ ionic selectivity arising from a net negative surface charge (Supplementary Fig. 16). The conductance of HACs drops by two orders of magnitude between 1 M and 1 mM KCl and the memristive effect can only be observed for concentrations above 100 mM (Supplementary Fig. 14). This points to the importance of the surface charge, which—for mechanically exfoliated pristine graphite as used in this work—increases with the salt concentration to reach an absolute value above 0.1 C m−2 at 1 M (ref. 13). Molecular dynamics simulations also showed that the interface of pristine graphite and hexagonal boron nitride gets electrified in water49. While the effect can still be observed with a hexagonal boron nitride cover layer, this is not the case with mica (Supplementary Fig. 15). We also performed in operando optical measurements with mica controls, directly showing considerably reduced blister formation occurring at negative voltages, as opposed to the large deformations observed for graphite at positive voltages only. The emergence of the effect with the graphite cap at high salt concentration (equivalent to high surface charge) and the highly reduced blister dynamics with mica indicates that the memory effect requires a high charge of the top layer interface. For these reasons, as well as because their spacing is larger than twice the Debye length (0.3 nm) at 1 M, we neglect the contribution of palladium islands on the ionic transport and consider them as neutral spacers. However, they play a role on the memristive effect through their adhesion energy.

The origin of memory in HACs is thus related to reversible mechanical deformation induced by converging counter-ion fluxes, explaining the observed charge threshold for switching. These ingredients have already been used to describe out-of-the-equilibrium and nonlinear ion transport phenomena at the nanoscale26,50. However, the large electrostatic pressure that can build up in nanoscale systems was never reported to trigger mechanical deformations at micron scale. HACs represent the optimal geometry to favour such effects by combining large asymmetry with single-digit confinement in the presence of a highly charged interface, resulting in dramatic mechano–ionic memristive dynamics.

Building logic circuits with interactive devices

Having established scalable nanofluidic device fabrication and obtained direct experimental proof regarding the underlying switching mechanism, we can utilize this knowledge to harness the potential of HACs in the realm of ionic computing. Building a logic circuit with two aqueous memristors influencing each other represents a paradigm shift for nanofluidics, where devices have until now been measured independently. We can now wire them for neuromorphic computing applications by taking advantage of the performances and reliability of HACs. This confers a new purpose to nanofluidic devices, in addition to their role as a technological platform for uncovering fundamentals of molecular transport at the smallest scales. Using two parallel HACs connected to a resistor as shown in Fig. 4a, we implement the material implication (IMP) logic gate as demonstrated by ref. 42 with solid-state memristors. The first HAC is called the P-switch, and its conductance state remains unchanged during the logic operation. The second HAC is the Q-switch; switching it from the OFF state (low conductance, defined as 0) to the ON state (high conductance, defined as 1) is only possible when P is in a low conductance state. Such conditional switching enables the implementation of the first two rows of the IMP truth table (Fig. 4b), which are the non-trivial cases. The two last rows are vacuous truth (or trivial cases) as Q is already switched ON at the start of the logic operation. In this demonstration of ionic computing, both nanofluidic memristors interact through the resistor to realize a conditional logic operation. The IMP gate represents a milestone for nanofluidic memristive action, as it can be used to derive any other classical logic gate commonly employed in digital computing42.Fig. 4 Nanofluidic logic.

a, Circuit schematic. Two HACs are connected in parallel with a variable resistor set to 6 MΩ. The working electrodes of each device are connected to the channels of a source-measurement unit. The shared ground electrode of the two cells is connected to the resistor. b, IMP truth table. The first two columns (P and Q) represent input states, and the neighbouring P' and Q' columns are the corresponding outputs of the IMP gate. A Greek letter is allocated to each one of the logic cases. c, Illustration of the working principle of conditional switching implementing the non-trivial cases of the IMP truth table: left, the α case and right, the β case. The effective voltage VQ applied to the Q-switch is sufficient to reach the charge threshold of the Q device within the pulse duration when P is in the 0 state (α case) and insufficient when P is in the 1 state (β case). d, Applied voltage and measured current for each logic case on both Q and P lines (Devices B and C, 1 M KCl). The measured pulses before the operation (with grey background) give the IMP table inputs and the ones after the operation provide the outputs. The Q-line current during operation reaches the threshold in the α case, indicated by an abrupt increase. e, Read current relative variations of the Q-line for the α case (rose background, P = 0) and β case, (blue background, P = 1). f, Read currents normalized by their respective minimum value in the α case both devices. The transition between the range corresponding to state 0 (blue background) and state 1 (rose background) occurs when the conductance of the corresponding device increases by at least 40% relative to its minimum value.

Source data

As HACs can be set and reset at least ten times (Supplementary Fig. 17), they are suitable for building logic circuits. We also checked the endurance of HACs under sinusoidal operation (Supplementary Figs. 18 and 23). We observe, through electrokinetic and optical measurement, perfect stability up to 17 cycles. We could also measure 300 cycles on another device, showing HACs’ robustness and setting an upper bound for nanofluidic memristors. Individual HACs can still operate when connected in series with the aforementioned resistor, even though this reduces the effective conductance ratio (Supplementary Fig. 28). To connect the switches, we use a single electrode with two chlorinated tips grounded through a variable resistor (Supplementary Fig. 27). The P and Q switches are both first placed into their targeted initial conductance state by applying a 2 s voltage pulse of high amplitude (+1 V for state 1 and −1 V for state 0). Input states of the memristors are then probed by applying read pulses to the working electrodes of each switch. The IMP gate is operated by applying voltage pulses simultaneously to each switch (Fig. 4c). In particular, a conditional voltage pulse of amplitude Vcond, which is not high enough to reach the switching threshold, is applied to the working electrode of the P-switch, while a larger voltage pulse, Vset is applied to the working electrode of the Q-switch. When the P-switch is in the OFF state, the current IPα flowing across it is of a lower value (α case shown in the rose (left) portion of Fig. 4c). Thus, the voltage drop across the resistor VRα=R(IPα+IQα) is sufficiently low, resulting in a large enough voltage drop VQα across Q to reach the charge threshold required to switch Q. Conversely, when P is in the ON state (β case shown in the right panel of Fig. 4c), the IPβ and VRβ are larger, making the voltage drop across Q (VQβ) insufficient for switching Q. This is why a switching threshold is essential for a successful implementation of IMP logic as illustrated in Supplementary Fig. 17, where we observe in single device experiments that setting HACs is only possible with 2 s programming pulses having an amplitude high enough for the blister to reach the ON-state regime. The value of the resistor R is chosen such that the potential drop given by the low conductance state current of P is low enough to enable the pulsed setting of the Q-switch. On the other hand, the value of R needs to provide a sufficient potential drop at the high-conductance state of P to avoid the setting of Q. R is therefore chosen such that it has a value between the devices resistances in the set and conditional states of the pulsed operation, RSET < R < RCOND (ref. 42). Finally, the output states are recorded similarly to the input states.

Figure 4d displays the experimental results obtained for the IMP gate implementation with memristive HACs. The Q-line current trace during the IMP gate operation reaches the switching threshold when P is set to 0 (α case) but not when P is set to 1 (β case). The conductance of Q increases by 50% in the α case and only 10% in the β case (Fig. 4e). We are thus able to condition the switching of a nanofluidic memristor by the conductance state of another device. We recover the IMP truth table by normalizing the read pulses of each memristor by their current minimum, while defining a posteriori an arbitrary borderline delimiting range corresponding to states 1 and 0 (Fig. 4f). Its existence serves as a proof of concept, demonstrating the potential to construct logic circuits using nanofluidic memristors as building blocks.

For robustness, we carry out another conditional switching experiment with two HACs and show switching with a HAC, as Q-switch and a variable resistor, as P-switch that emulates another HAC (Supplementary Figs. 29 and 30). We have also performed numerical simulations of an equivalent electrical circuit consisting of two charge-threshold memristive switches with typical conductance values and ratio, as in our pulse programming experiments. Qualitatively, simulation reproduces experimental results, showing that IMP logic is indeed achievable with HACs (Supplementary Fig. 31).

Logic operations with HACs are made possible by the presence of a charge switching threshold, a fast speed enabling operations with 2 s programming pulses, as well as a large conductance ratio compensating the performance loss caused by the presence of the resistor. The combination of these features is exclusive to HACs, making logic gating unattainable with previously reported nanofluidic memristors.

Conclusions

Our HACs are scalable and compact nanofluidic memristors that can be set or reset within a few seconds, with a conductance ratio reaching sixty. In operando optical observations show that their large entrance asymmetry and single-digit confinement result in a switching behaviour related to the reversible formation of liquid blisters. This suggests several parameters that can be tuned for further optimization. In particular, choosing the membrane size and stiffness, the surface charge and adhesion properties of materials, and developing guidelines to control the blister position and dynamics will help increase performance and reliability, as well as reduce the device-to-device variation.

The strongly nonlinear dynamics of the HACs allow us to implement a Boolean operation with two interacting devices, providing the fundamental building block for future aqueous computing machines. Our HACs cannot currently match the performance of solid-state Pt/TiO2/Pt memristive switches, which operate at the microsecond timescale42, but further developments—such as optimizing their design and connecting them with water channels to fabricate fully liquid circuits—should lead to improvements. By taking inspiration from electronic crossbar arrays, as well as the brains of living organisms, our approach could potentially lead to the creation of nanofluidic neural networks.

Methods

Device fabrication

To fabricate HACs, we start with homemade SiN membranes having a single hole of approximately 100 nm width at their centre (Fig. 1a, step 0 and Supplementary Fig. 1). Using electron-beam evaporation, we then deposit a discontinuous layer of palladium that self-organizes into clusters or ‘islands’, leaving space for ions to flow around them (Fig. 1a, step 1 and Supplementary Fig. 2). Finally, we transfer a bidimensional crystal above the hole to close the system using a droplet-shaped polydimethylsiloxane stamp covered with polypropylene carbonate (Fig. 1a, step 2 and Supplementary Fig. 3). Here, we mostly used graphite but also tested the effect of using hexagonal boron nitride and mica as the top layer material. Further details are given in Supplementary Information Section 1.

Nanofluidic measurements

Ag/AgCl electrodes are used to apply the potential and measure the resulting current. Data acquisition setup for single device measurements consists of a low noise current amplifier (FEMTO DLPCA-200, Electro Optical Components) and a digital-to-analogue converter (NI 63 series, NI). The measurements are sampled at 100 kHz with an acquisition time of 100 ms. We used KCl, CaCl2 and AlCl3 solutions with concentrations ranging from 1 mM to 1 M. For logic experiments with two channels, the working electrodes of each cell are connected to a dual channel source-measurement unit (2636B, Keithley) using a sampling rate of 10 Hz.

Supplementary information

Supplementary Information Supplementary Figs. 1–30, Discussion and Tables 1 and 2.

Supplementary Video 1 In operando measure presented in Fig. 3 (sinusoidal).

Supplementary Video 2 In operando measurement presented in Supplementary Fig. 20 (sinusoidal).

Supplementary Video 3 In operando measurement supplementary (sinusoidal).

Supplementary Video 4 In operando measurement supplementary (sinusoidal).

Supplementary Video 5 In operando measurement presented in Supplementary Fig. 21 (pulses–retention).

Supplementary Video 6 In operando measurement presented in Supplementary Fig. 22 (mica control 1).

Supplementary Video 7 In operando measurement presented in Supplementary Fig. 22 (mica control 2).

Supplementary Data Statistical source data for Supplementary figures.

Source data

Source Data Fig. 1 Statistical source data.

Source Data Fig. 2 Statistical source data.

Source Data Fig. 3 Statistical source data.

Source Data Fig. 4 Statistical source data.

Supplementary information

The online version contains supplementary material available at 10.1038/s41928-024-01137-9.

Acknowledgements

We thank the EPFL Center of MicroNanoTechnology (CMi) for their help with chip fabrication and the EPFL Center for Electron Microscopy (CIME) for access to electron microscopes, in particular L. Navratilova and V. Boureau, for assistance with TEM cross-sectional imaging. We thank E. Matioli and G. M. Marega and B. Coquinot for fruitful discussions. Y.T. thanks NCCR Bio-Inspired Materials. Funding: European Union’s H2020 Framework Programme/ERC Advanced Grant agreement number 101020445 (T.E., N.R. and A.R). Swiss National Science Foundation grant No. TMPFP2-217134 (T.E.). European Union’s Horizon 2020 research and innovation programme under grant agreements No 894369 (Marie Curie Sklodowska ITN network ‘2-Exciting’) and No. 881603 (Graphene Flagship Core 3 Phase) as well as the Swiss National Science Foundation (grants no. 157739 and 205114) (E.L., R.C. and A.K). NSF grant No. ECCS-2229880 (M.D.).

Author contributions

A.R. and T.E. conceived the project and designed the experiments. A.R. supervised the Laboratory of Nanoscale Biology team and A.K. supervised the Laboratory of Nanoscale Electronics and Structures team. N.R. built the microscope and T.E built the fluidic cell for in operando optical measurements. E.L. and R.C. built the ionic logic setup with inputs from Y.T. and T.E. T.E., Y.T. and N.R. performed the measurements with input from E.L. and R.C. Y.T. made the SiN windows and performed palladium deposition. T.E. transferred the 2D crystals. A.C. did the TEM cross-sectional image. N.R. developed the theoretical modelling with inputs from T.E. and V.A. T.E and N.R. analysed the experimental data. Y.T. performed the numerical simulations with inputs from A.R. T.E. wrote the paper with inputs from Y.T. and N.R. All authors contributed to editing of the paper.

Peer review

Peer review information

Nature Electronics thanks Boya Radha and Ping Yu for their contribution to the peer review of this work.

Funding

Open access funding provided by EPFL Lausanne.

Data availability

Source data are provided with this paper. Additional data related to this work are available from the corresponding authors upon request.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Theo Emmerich, Yunfei Teng, Nathan Ronceray.
==== Refs
References

1. Voglis G Tavernarakis N The role of synaptic ion channels in synaptic plasticity EMBO Rep. 2006 7 1104 1110 10.1038/sj.embor.7400830 17077866
Voglis, G. & Tavernarakis, N. The role of synaptic ion channels in synaptic plasticity. EMBO Rep. 7, 1104–1110 (2006).17077866 10.1038/sj.embor.7400830
2. Noy A Darling SB Nanofluidic computing makes a splash Science 2023 379 143 144 10.1126/science.adf6400 36634195
Noy, A. & Darling, S. B. Nanofluidic computing makes a splash. Science 379, 143–144 (2023).36634195 10.1126/science.adf6400
3. Howard N Energy paradox of the brain Brain Sci. 2012 1 35
Howard, N. Energy paradox of the brain. Brain Sci. 1, 35 (2012).
4. Kavokine N Netz RR Bocquet L Fluids at the nanoscale: from continuum to subcontinuum transport Annu. Rev. Fluid Mech. 2021 53 377 410 10.1146/annurev-fluid-071320-095958
Kavokine, N., Netz, R. R. & Bocquet, L. Fluids at the nanoscale: from continuum to subcontinuum transport. Annu. Rev. Fluid Mech. 53, 377–410 (2021).10.1146/annurev-fluid-071320-095958
5. Garaj S Graphene as a subnanometre trans-electrode membrane Nature 2010 467 190 193 10.1038/nature09379 20720538
Garaj, S. et al. Graphene as a subnanometre trans-electrode membrane. Nature 467, 190–193 (2010).20720538 10.1038/nature09379
6. Feng J Single-layer MoS2 nanopores as nanopower generators Nature 2016 536 197 200 10.1038/nature18593 27409806
Feng, J. et al. Single-layer MoS2 nanopores as nanopower generators. Nature 536, 197–200 (2016).27409806 10.1038/nature18593
7. Li Z Breakdown of the Nernst-Einstein relation in carbon nanotube porins Nat. Nanotechnol. 2023 18 177 183 10.1038/s41565-022-01276-0 36585518
Li, Z. et al. Breakdown of the Nernst-Einstein relation in carbon nanotube porins. Nat. Nanotechnol. 18, 177–183 (2023).36585518 10.1038/s41565-022-01276-0
8. Tunuguntla RH Enhanced water permeability and tunable ion selectivity in subnanometer carbon nanotube porins Science 2017 357 792 796 10.1126/science.aan2438 28839070
Tunuguntla, R. H. et al. Enhanced water permeability and tunable ion selectivity in subnanometer carbon nanotube porins. Science 357, 792–796 (2017).28839070 10.1126/science.aan2438
9. Secchi E Massive radius-dependent flow slippage in carbon nanotubes Nature 2016 537 210 213 10.1038/nature19315 27604947
Secchi, E. et al. Massive radius-dependent flow slippage in carbon nanotubes. Nature 537, 210–213 (2016).27604947 10.1038/nature19315
10. Radha B Molecular transport through capillaries made with atomic-scale precision Nature 2016 538 222 225 10.1038/nature19363 27602512
Radha, B. et al. Molecular transport through capillaries made with atomic-scale precision. Nature 538, 222–225 (2016).27602512 10.1038/nature19363
11. Esfandiar A Size effect in ion transport through angstrom-scale slits Science 2017 358 511 513 10.1126/science.aan5275 29074772
Esfandiar, A. et al. Size effect in ion transport through angstrom-scale slits. Science 358, 511–513 (2017).29074772 10.1126/science.aan5275
12. Gopinadhan K Complete steric exclusion of ions and proton transport through confined monolayer water Science 2019 363 145 148 10.1126/science.aau6771 30630924
Gopinadhan, K. et al. Complete steric exclusion of ions and proton transport through confined monolayer water. Science 363, 145–148 (2019).30630924 10.1126/science.aau6771
13. Emmerich T Enhanced nanofluidic transport in activated carbon nanoconduits Nat. Mater. 2022 21 696 702 10.1038/s41563-022-01229-x 35422506
Emmerich, T. et al. Enhanced nanofluidic transport in activated carbon nanoconduits. Nat. Mater. 21, 696–702 (2022).35422506 10.1038/s41563-022-01229-x
14. Agre P Aquaporin water channels–from atomic structure to clinical medicine J. Physiol. 2002 542 3 16 10.1113/jphysiol.2002.020818 12096044
Agre, P. et al. Aquaporin water channels–from atomic structure to clinical medicine. J. Physiol. 542, 3–16 (2002).12096044 10.1113/jphysiol.2002.020818
15. Gouaux E MacKinnon R Principles of selective ion transport in channels and pumps Science 2005 310 1461 1465 10.1126/science.1113666 16322449
Gouaux, E. & MacKinnon, R. Principles of selective ion transport in channels and pumps. Science 310, 1461–1465 (2005).16322449 10.1126/science.1113666
16. Coetzee W Molecular diversity of K+ channels Ann. N. Y. Acad. Sci. 1999 868 233 255 10.1111/j.1749-6632.1999.tb11293.x 10414301
Coetzee, W. et al. Molecular diversity of K+ channels. Ann. N. Y. Acad. Sci. 868, 233–255 (1999).10414301 10.1111/j.1749-6632.1999.tb11293.x
17. Ajo-Franklin CM Noy A Crossing over: nanostructures that move electrons and ions across cellular membranes Adv. Mater. 2015 27 5797 5804 10.1002/adma.201500344 25914282
Ajo-Franklin, C. M. & Noy, A. Crossing over: nanostructures that move electrons and ions across cellular membranes. Adv. Mater. 27, 5797–5804 (2015).25914282 10.1002/adma.201500344
18. Lanyi JK Bacteriorhodopsin Annu. Rev. Physiol. 2004 66 665 688 10.1146/annurev.physiol.66.032102.150049 14977418
Lanyi, J. K. Bacteriorhodopsin. Annu. Rev. Physiol. 66, 665–688 (2004).14977418 10.1146/annurev.physiol.66.032102.150049
19. Busath D Szabo G Gramicidin forms multi-state rectifying channels Nature 1981 294 371 373 10.1038/294371a0 6171731
Busath, D. & Szabo, G. Gramicidin forms multi-state rectifying channels. Nature 294, 371–373 (1981).6171731 10.1038/294371a0
20. Holt JK Fast mass transport through sub-2-nanometer carbon nanotubes Science 2006 312 1034 1037 10.1126/science.1126298 16709781
Holt, J. K. et al. Fast mass transport through sub-2-nanometer carbon nanotubes. Science 312, 1034–1037 (2006).16709781 10.1126/science.1126298
21. Xie Q Fast water transport in graphene nanofluidic channels Nat. Nanotechnol. 2018 13 238 245 10.1038/s41565-017-0031-9 29292381
Xie, Q. et al. Fast water transport in graphene nanofluidic channels. Nat. Nanotechnol. 13, 238–245 (2018).29292381 10.1038/s41565-017-0031-9
22. Xiao K Artificial light-driven ion pump for photoelectric energy conversion Nat. Commun. 2019 10 74 10.1038/s41467-018-08029-5 30622279
Xiao, K. et al. Artificial light-driven ion pump for photoelectric energy conversion. Nat. Commun. 10, 74 (2019).30622279 10.1038/s41467-018-08029-5
23. Zhang, H. et al. Bioinspired artificial single ion pump. J. Am. Chem. Soc. 135, 16102–16110 (2013).
24. Xu Y Yazbeck R Duan C Anomalous mechanosensitive ion transport in nanoparticle-blocked nanopores J. Chem. Phys. 2021 154 224702 10.1063/5.0046086 34241226
Xu, Y., Yazbeck, R. & Duan, C. Anomalous mechanosensitive ion transport in nanoparticle-blocked nanopores. J. Chem. Phys. 154, 224702 (2021).34241226 10.1063/5.0046086
25. Davis SJ Pressure-induced enlargement and ionic current rectification in symmetric nanopores Nano Lett. 2020 20 8089 8095 10.1021/acs.nanolett.0c03083 33048551
Davis, S. J. et al. Pressure-induced enlargement and ionic current rectification in symmetric nanopores. Nano Lett. 20, 8089–8095 (2020).33048551 10.1021/acs.nanolett.0c03083
26. Xiong T Neuromorphic functions with a polyelectrolyte-confined fluidic memristor Science 2023 379 156 161 10.1126/science.adc9150 36634194
Xiong, T. et al. Neuromorphic functions with a polyelectrolyte-confined fluidic memristor. Science 379, 156–161 (2023).36634194 10.1126/science.adc9150
27. Robin P Long-term memory and synapse-like dynamics in two-dimensional nanofluidic channels Science 2023 379 161 167 10.1126/science.adc9931 36634187
Robin, P. et al. Long-term memory and synapse-like dynamics in two-dimensional nanofluidic channels. Science 379, 161–167 (2023).36634187 10.1126/science.adc9931
28. Robin P Bocquet L Nanofluidics at the crossroads J. Chem. Phys. 2023 158 160901 10.1063/5.0143222 37093148
Robin, P. & Bocquet, L. Nanofluidics at the crossroads. J. Chem. Phys. 158, 160901 (2023).37093148 10.1063/5.0143222
29. Hou Y Learning from the brain: bioinspired nanofluidics J. Phys. Chem. Lett. 2023 14 2891 2900 10.1021/acs.jpclett.2c03930 36927003
Hou, Y. et al. Learning from the brain: bioinspired nanofluidics. J. Phys. Chem. Lett. 14, 2891–2900 (2023).36927003 10.1021/acs.jpclett.2c03930
30. Li C Synaptic iontronic devices for brain-mimicking functions: fundamentals and applications ACS Appl. Bio Mater. 2020 4 71 84 10.1021/acsabm.0c00806
Li, C. et al. Synaptic iontronic devices for brain-mimicking functions: fundamentals and applications. ACS Appl. Bio Mater. 4, 71–84 (2020).10.1021/acsabm.0c00806
31. Xie B Perspective on nanofluidic memristors: from mechanism to application Chem. Asian J. 2022 17 e202200682 10.1002/asia.202200682 35994236
Xie, B. et al. Perspective on nanofluidic memristors: from mechanism to application. Chem. Asian J. 17, e202200682 (2022).35994236 10.1002/asia.202200682
32. Pershin YV Kim J Datta T Di Ventra M An experimental demonstration of the memristor test Physica E 2022 142 115290 10.1016/j.physe.2022.115290
Pershin, Y. V., Kim, J., Datta, T. & Di Ventra, M. An experimental demonstration of the memristor test. Physica E 142, 115290 (2022).10.1016/j.physe.2022.115290
33. Pershin YV Di Ventra M Experimental demonstration of associative memory with memristive neural networks Neural Netw. 2010 23 881 886 10.1016/j.neunet.2010.05.001 20605401
Pershin, Y. V. & Di Ventra, M. Experimental demonstration of associative memory with memristive neural networks. Neural Netw. 23, 881–886 (2010).20605401 10.1016/j.neunet.2010.05.001
34. Yao P Fully hardware-implemented memristor convolutional neural network Nature 2020 577 641 646 10.1038/s41586-020-1942-4 31996818
Yao, P. et al. Fully hardware-implemented memristor convolutional neural network. Nature 577, 641–646 (2020).31996818 10.1038/s41586-020-1942-4
35. Sheng Q Xie Y Li J Wang X Xue J Transporting an ionic-liquid/water mixture in a conical nanochannel: a nanofluidic memristor Chem. Commun. 2017 53 6125 6127 10.1039/C7CC01047H
Sheng, Q., Xie, Y., Li, J., Wang, X. & Xue, J. Transporting an ionic-liquid/water mixture in a conical nanochannel: a nanofluidic memristor. Chem. Commun. 53, 6125–6127 (2017).10.1039/C7CC01047H
36. Bu Y Ahmed Z Yobas L A nanofluidic memristor based on ion concentration polarization Analyst 2019 144 7168 7172 10.1039/C9AN01561B 31724657
Bu, Y., Ahmed, Z. & Yobas, L. A nanofluidic memristor based on ion concentration polarization. Analyst 144, 7168–7172 (2019).31724657 10.1039/C9AN01561B
37. Zhang P Nanochannel-based transport in an interfacial memristor can emulate the analog weight modulation of synapses Nano Lett. 2019 19 4279 4286 10.1021/acs.nanolett.9b00525 31150262
Zhang, P. et al. Nanochannel-based transport in an interfacial memristor can emulate the analog weight modulation of synapses. Nano Lett. 19, 4279–4286 (2019).31150262 10.1021/acs.nanolett.9b00525
38. Evans RC Blackwell KT Calcium: amplitude, duration, or location? Biol. Bull. 2015 228 75 83 10.1086/BBLv228n1p75 25745102
Evans, R. C. & Blackwell, K. T. Calcium: amplitude, duration, or location? Biol. Bull. 228, 75–83 (2015).25745102 10.1086/BBLv228n1p75
39. Kukushkin NV Carew TJ Memory takes time Neuron 2017 95 259 279 10.1016/j.neuron.2017.05.029 28728021
Kukushkin, N. V. & Carew, T. J. Memory takes time. Neuron 95, 259–279 (2017).28728021 10.1016/j.neuron.2017.05.029
40. Li C Analogue signal and image processing with large memristor crossbars Nat. Electron. 2018 1 52 59 10.1038/s41928-017-0002-z
Li, C. et al. Analogue signal and image processing with large memristor crossbars. Nat. Electron. 1, 52–59 (2018).10.1038/s41928-017-0002-z
41. Chu M Neuromorphic hardware system for visual pattern recognition with memristor array and CMOS neuron IEEE Trans. Ind. Electron. 2014 62 2410 2419 10.1109/TIE.2014.2356439
Chu, M. et al. Neuromorphic hardware system for visual pattern recognition with memristor array and CMOS neuron. IEEE Trans. Ind. Electron. 62, 2410–2419 (2014).10.1109/TIE.2014.2356439
42. Borghetti J ‘Memristive’ switches enable ‘stateful’ logic operations via material implication Nature 2010 464 873 876 10.1038/nature08940 20376145
Borghetti, J. et al. ‘Memristive’ switches enable ‘stateful’ logic operations via material implication. Nature 464, 873–876 (2010).20376145 10.1038/nature08940
43. Bocquet L Nanofluidics coming of age Nat. Mater. 2020 19 254 256 10.1038/s41563-020-0625-8 32099111
Bocquet, L. Nanofluidics coming of age. Nat. Mater. 19, 254–256 (2020).32099111 10.1038/s41563-020-0625-8
44. Xiong T Li W Yu P Mao L Fluidic memristor: bringing chemistry to neuromorphic devices Innovation 2023 4 100435 37215530
Xiong, T., Li, W., Yu, P. & Mao, L. Fluidic memristor: bringing chemistry to neuromorphic devices. Innovation 4, 100435 (2023).37215530
45. Wang FZ Delayed switching in memristors and memristive systems IEEE Electron Device Lett. 2010 31 755 757 10.1109/LED.2010.2049560
Wang, F. Z. et al. Delayed switching in memristors and memristive systems. IEEE Electron Device Lett. 31, 755–757 (2010).10.1109/LED.2010.2049560
46. Valov I Waser R Jameson JR Kozicki MN Electrochemical metallization memories—fundamentals, applications, prospects Nanotechnology 2011 22 254003 10.1088/0957-4484/22/25/254003 21572191
Valov, I., Waser, R., Jameson, J. R. & Kozicki, M. N. Electrochemical metallization memories—fundamentals, applications, prospects. Nanotechnology 22, 254003 (2011).21572191 10.1088/0957-4484/22/25/254003
47. Boyn S Learning through ferroelectric domain dynamics in solid-state synapses Nat. Commun. 2017 8 14736 10.1038/ncomms14736 28368007
Boyn, S. et al. Learning through ferroelectric domain dynamics in solid-state synapses. Nat. Commun. 8, 14736 (2017).28368007 10.1038/ncomms14736
48. Ronceray, N. et al. Elastocapillarity-driven 2D nano-switches enable zeptoliter-scale liquid encapsulation. Nature Commun. 15, 185 (2024).
49. Grosjean B Bocquet ML Vuilleumier R Versatile electrification of two-dimensional nanomaterials in water Nat. Commun. 2019 10 1656 10.1038/s41467-019-09708-7 30971700
Grosjean, B., Bocquet, M. L. & Vuilleumier, R. Versatile electrification of two-dimensional nanomaterials in water. Nat. Commun. 10, 1656 (2019).30971700 10.1038/s41467-019-09708-7
50. Jubin L Poggioli A Siria A Bocquet L Dramatic pressure-sensitive ion conduction in conical nanopores Proc. Natl Acad. Sci. USA 2018 115 4063 4068 10.1073/pnas.1721987115 29610303
Jubin, L., Poggioli, A., Siria, A. & Bocquet, L. Dramatic pressure-sensitive ion conduction in conical nanopores. Proc. Natl Acad. Sci. USA 115, 4063–4068 (2018).29610303 10.1073/pnas.1721987115
