Search

Fluidic iontronics offer a unique capability for emulating the chemical processes found in neurons. We extract multiple distinct chemically regulated synaptic features from a single conical microfluidic channel carrying functionalized surface groups, using finite-element calculations of continuum transport equations. Such channels have long been employed for fluidic sensing and are therefore experimentally well established. By modeling a Langmuir-type surface reaction on the channel wall we couple fast voltage-induced volumetric salt accumulation with a long-term channel surface charge modulation by means of fast charging and slow discharging. These nonlinear charging dynamics are understood through an analytic approximation rooted in first-principles. We show how short-and long-term potentiation and depression, frequency-dependent plasticity, and chemical-electrical signal coincidence detection (acting like a chemical-electrical AND logic gate), akin to the NMDA mechanism for Hebbian learning in biological synapses, can all be emulated with a single channel.

In this paper, building upon the discovery of asymmetric rectified electric fields (AREF) in recent experiments [S.H. Hashemi et al., Physical Review Letters 121, 185504 (2018)], we explore the generation of AREF by applying a sawtooth-like voltage to 1:1 electrolytes with equal diffusion coefficients confined between two planar blocking electrodes. This differs from an earlier approach based on a sinusoidal AC voltage applied to 1:1 electrolytes with unequal diffusion coefficients. By numerically solving the full Poisson-Nernst-Planck equations, we demonstrate that AREF can be generated by a slow rise and a fast drop of the potential (or vice versa), even for electrolytes with equal diffusion coefficients of the cations and anions. We employ an analytically constructed equivalent electric circuit to explain the underlying physical mechanism. Importantly, we find that the strength of AREF can be effectively tuned from zero to its maximal value by only manipulating the time-dependence of the driving voltage, eliminating the necessity to modify the electrolyte composition between experiments. This provides valuable insights to control the manipulation of AREF, which facilitates enhanced applications in diverse electrochemical systems.

A hitherto unexploited characteristic feature of emerging iontronic devices for information processing is the intrinsic mobility of the medium (water) of dissolved ions in aqueous electrolytes, which therefore not only respond to voltage but also to pressure. Here we study a microfluidic memristor, in the form of a conical channel, exposed to simultaneously applied time-dependent voltage and pressure drops, through numerical solutions of the Poisson-Nernst-Planck-Stokes equations for ion and fluid transport. We show that the channel's memristive properties can be enhanced, reduced or instantaneously reset by a suitable pressure, and we leverage this finding with two examples of time series processing of simultaneously applied voltage and pressure pulses. We not only show that the distinction between different voltage time series can be improved by enhancing the conductance response with corresponding pressure pulses, but also that the bandwidth of information transfer through the channel can be doubled by letting the pressure pulses represent a second independent time series.

In pursuit of neuromorphic (brain-inspired) devices, memristors (memory-resistors) have emerged as effective components for emulating neuronal circuitry. Here we formally define a class of Simple Volatile Memristors (SVMs) based on a simple conductance equation of motion from which we build a simple mathematical theory on the dynamics of isolated SVMs and SVM-based spiking circuits. Notably, SVMs include various fluidic iontronic devices that have recently garnered significant interest due to their unique quality of operating within the same medium as the brain. Specifically we show that symmetric SVMs produce non self-crossing current-voltage hysteresis loops, while asymmetric SVMs produce self-crossing loops. Additionally, we derive a general expression for the enclosed area in a loop, providing a relation between the voltage frequency and the SVM memory timescale. These general results are shown to materialise in physical finite-element calculations of microfluidic memristors. An SVM-based circuit has been proposed that exhibits all-or-none and tonic neuronal spiking. We generalise and analyse this spiking circuit, characterising it as a two-dimensional dynamical system. Moreover, we demonstrate that stochastic effects can induce novel neuronal firing modes absent in the deterministic case. Through our analysis, the circuit dynamics are well understood, while retaining its explicit link with the physically plausible underlying system.

Fluidic iontronics is emerging as a distinctive platform for implementing neuromorphic circuits, characterized by its reliance on the same aqueous medium and ionic signal carriers as the brain. Drawing upon recent theoretical advancements in both iontronic spiking circuits and in dynamic conductance of conical ion channels, which form fluidic memristors, we expand the repertoire of proposed neuronal spiking dynamics in iontronic circuits. Through a modelled circuit containing channels that carry a bipolar surface charge, we extract phasic bursting, mixed-mode spiking, tonic bursting, and threshold variability, all with spike voltages and frequencies within the typical range for mammalian neurons. These features are possible due to the strong dependence of the typical conductance memory retention time on the channel length, enabling timescales varying from individual spikes to bursts of multiple spikes within a single circuit. These advanced forms of neuronal-like spiking support the exploration of aqueous iontronics as an interesting platform for neuromorphic circuits.

The brain's remarkable and efficient information processing capability is driving research into brain-inspired (neuromorphic) computing paradigms. Artificial aqueous ion channels are emerging as an exciting platform for neuromorphic computing, representing a departure from conventional solid-state devices by directly mimicking the brain's fluidic ion transport. Supported by a quantitative theoretical model, we present easy to fabricate tapered microchannels that embed a conducting network of fluidic nanochannels between a colloidal structure. Due to transient salt concentration polarisation our devices are volatile memristors (memory resistors) that are remarkably stable. The voltage-driven net salt flux and accumulation, that underpin the concentration polarisation, surprisingly combine into a diffusionlike quadratic dependence of the memory retention time on the channel length, allowing channel design for a specific timescale. We implement our device as a synaptic element for neuromorphic reservoir computing. Individual channels distinguish various time series, that together represent (handwritten) numbers, for subsequent in-silico classification with a simple readout function. Our results represent a significant step towards realising the promise of fluidic ion channels as a platform to emulate the rich aqueous dynamics of the brain.

Conical channels filled with an aqueous electrolyte have been proposed as promising candidates for iontronic neuromorphic circuits. This is facilitated by a novel analytical model for the internal channel dynamics [Kamsma et al., arXiv:2301.06158, 2023], the relative ease of fabrication of conical channels, and the wide range of achievable memory retention times by varying the channel lengths. In this work, we demonstrate that the analytical model for conical channels can be generalized to channels with an inhomogeneous surface charge distribution, which we predict to exhibit significantly stronger current rectification and more pronounced memristive properties in the case of bipolar channels, i.e. channels where the tip and base carry a surface charge of opposite sign. Additionally, we show that the use of bipolar conical channels in a previously proposed iontronic circuit features hallmarks of neuronal communication, such as all-or-none action potentials and spike train generation. Bipolar channels allow, however, for circuit parameters in the range of their biological analogues, and exhibit membrane potentials that match well with biological mammalian action potentials, further supporting its potential for bio-compatibility.

Experiments have shown that the conductance of conical channels, filled with an aqueous electrolyte, can strongly depend on the history of the applied voltage. These channels hence have a memory and are promising elements in brain-inspired (iontronic) circuits. We show here that the memory of such channels stems from transient concentration polarization over the ionic diffusion time. We derive an analytic approximation for these dynamics which shows good agreement with full finite-element calculations. Using our analytic approximation, we propose an experimentally realisable Hodgkin-Huxley iontronic circuit where micrometer cones take on the role of sodium and potassium channels. Our proposed circuit exhibits key features of neuronal communication such as all-or-none action potentials upon a pulse stimulus and a spike train upon a sustained stimulus.

Inspired by recent experimental observations of anomalously large decay lengths in concentrated electrolytes, we revisit the Restricted Primitive Model (RPM) for an aqueous electrolyte. We investigate the asymptotic decay lengths of the one-body ionic density profiles for the RPM in contact with a planar electrode using classical Density Functional Theory (DFT), and compare these with the decay lengths of the corresponding two-body correlation functions in bulk systems, obtained in previous Integral Equation Theory (IET) studies. Extensive Molecular Dynamics (MD) simulations are employed to complement the DFT and IET predictions. Our DFT calculations incorporate electrostatic interactions between the ions using three different (existing) approaches: one based on the simplest mean field treatment of Coulomb interactions (MFC), whilst the other two employ the Mean Spherical approximation (MSA). The MSAc invokes only the MSA bulk direct correlation function whereas the MSAu also incorporates the MSA bulk internal energy. Although MSAu yields profiles that agree best with MD simulations in the near field, in the far field we observe that the decay lengths are consistent between IET, MSAc, and MD simulations, whereas those from MFC and MSAu deviate significantly. Using DFT we calculated the solvation force, which relates directly to surface force experiments. We find that its decay length is neither qualitatively nor quantitatively close to the large decay lengths measured in experiments and conclude that the latter cannot be accounted for by the primitive model. The anomalously large decay lengths found in surface force measurements require an explanation that lies beyond primitive models.

We theoretically study the electrokinetic transport properties of nano-fluidic devices under the influence of a pressure, voltage or salinity gradient. On a microscopic level the behaviour of the device is quantified by the Onsager matrix ${\bf L}$, a generalised conductivity matrix relating the local driving forces and the induced volume, charge and salt flux. Extending ${\bf L}$ from a local to a global linear-response relation is trivial for homogeneous electrokinetic systems, but in this manuscript we derive a generalised conductivity matrix ${\bf G}$ from ${\bf L}$ that applies also to heterogeneous electrokinetic systems. This extension is especially important in the case of an imposed salinity gradient, which gives necessarily rise to heterogeneous devices. Within this formalism we can also incorporate a heterogeneous surface charge due to, for instance, a charge regulating boundary condition, which we show to have a significant impact on the resulting fluxes. The predictions of the Poisson-Nernst-Planck-Stokes theory show good agreement with exact solutions of the governing equations determined using the Finite Element Method under a wide variety of parameters. Having established the validity of the theory, it provides an accessible method to analyse electrokinetic systems in general without the need of extensive numerical methods. As an example, we analyse a Reverse Electrodialysis "blue energy" system, and analyse how the many parameters that characterise such a system affect the generated electrical power and efficiency.

We theoretically study the electrokinetic problem of a pressure-induced liquid flow through a narrow long channel with charged walls, going beyond the classical Helmholtz-Schmolukowski picture by considering the surprisingly strong combined effect of (i) Stern layer conductance and (ii) dynamic charge-regulating rather than fixed surface charges. We find that the water flow induces, apart from the well-known streaming potential, also a strongly heterogeneous surface charge and zeta-potential on chemically homogeneous channel walls. Moreover, we identify a novel steady state with a nontrivial 3D electric flux with 2D surface charges acting as sources and sinks. For a pulsed pressure drop our findings also provide a first-principles explanation for ill-understood experiments on the effect of flow on interfacial chemistry [D. Lis et al., Science 344, 1138 (2014)]., Comment: Accepted to be published in Physical Review Letters

We revisit the classical problem of flow of electrolyte solutions through charged capillary nanopores or nanotubes as described by the capillary pore model (also called "space charge" theory). This theory assumes very long and thin pores and uses a one-dimensional flux-force formalism which relates fluxes (electrical current, salt flux, fluid velocity) and driving forces (difference in electric potential, salt concentration, pressure). We analyze the general case with overlapping electric double layers in the pore and a nonzero axial salt concentration gradient. The 3$\times$3 matrix relating these quantities exhibits Onsager symmetry and we report a significant new simplification for the diagonal element relating axial salt flux to the gradient in chemical potential. We prove that Onsager symmetry is preserved under changes of variables, which we illustrate by transformation to a different flux-force matrix given by Gross and Osterle (1968). The capillary pore model is well-suited to describe the nonlinear response of charged membranes or nanofluidic devices for electrokinetic energy conversion and water desalination, as long as the transverse ion profiles remain in local quasi-equilibrium. As an example, we evaluate electrical power production from a salt concentration difference by reverse electrodialysis, using an efficiency vs. power diagram. We show that since the capillary pore model allows for axial gradients in salt concentration, partial loops in current, salt flux or fluid flow can develop in the pore. Predictions for macroscopic transport properties using a reduced model where the potential and concentration are assumed to be invariant with radial coordinate ("uniform potential" or "fine capillary pore" model), are close to results of the full model.

Even though chiral nematic phases were the first liquid crystals experimentally observed more than a century ago, the origin of the thermodynamic stability of cholesteric states is still unclear. In this Letter we address the problem by means of a novel density functional theory for the equilibrium pitch of chiral particles. When applied to right-handed hard helices, our theory predicts an entropy-driven cholesteric phase, which can be either right- or left-handed, depending not only on the particle shape but also on the thermodynamic state. We explain the origin of the chiral ordering as an interplay between local nematic alignment and excluded-volume differences between left- and right-handed particle pairs.

Due to remarkable advances in colloid synthesis techniques, systems of squares and cubes, once an academic abstraction for theorists and simulators, are nowadays an experimental reality. By means of a free minimization of the free-energy functional, we apply Fundamental Measure Theory to analyze the phase behavior of parallel hard squares and hard cubes. We compare our results with those obtained by the traditional approach based on the Gaussian parameterization, finding small deviations and good overall agreement between the two methods. For hard squares our predictions feature at intermediate packing fraction a smectic phase, which is however expected to be unstable due to thermal fluctuations. This implies that for hard squares the theory predicts either a vacancy-rich second-order transition or a vacancy-poor weakly first-order phase transition at higher density. In accordance with previous studies, a second-order transition with a high vacancy concentration is predicted for hard cubes.

Inspired by the observations of a remarkably stable biaxial nematic phase [E.v.d. Pol et al., Phys. Rev. Lett. 103, 258301 (2009)], we investigate the effect of size polydispersity on the phase behavior of a suspension of boardlike particles. By means of Onsager theory within the restricted orientation (Zwanzig) model we show that polydispersity induces a novel topology in the phase diagram, with two Landau tetracritical points in between which oblate uniaxial nematic order is favored over the expected prolate order. Additionally, this phenomenon causes the opening of a huge stable biaxiality regime in between uniaxial nematic and smectic states.

A huge amount of entropy is produced at places where fresh water and seawater mix, for example at river mouths. This mixing process is a potentially enormous source of sustainable energy, provided it is harnessed properly, for instance by a cyclic charging and discharging process of porous electrodes immersed in salt and fresh water, respectively [D. Brogioli, Phys. Rev. Lett. 103, 058501 (2009)]. Here we employ a modified Poisson-Boltzmann free-energy density functional to calculate the ionic adsorption and desorption onto and from the charged electrodes, from which the electric work of a cycle is deduced. We propose optimal (most efficient) cycles for two given salt baths involving two canonical and two grand-canonical (dis)charging paths, in analogy to the well-known Carnot cycle for heat-to-work conversion from two heat baths involving two isothermal and two adiabatic paths. We also suggest a slightly modified cycle which can be applied in cases that the stream of fresh water is limited., Comment: 7 Figures

We apply the recently developed triangular tessellation technique as presented in [J. de Graaf et al., Phys. Rev. E 80, 051405 (2009)] to calculate the free energy associated with the adsorption of anisotropic colloidal particles at a flat interface. From the free-energy landscape, we analyze the adsorption process, using a simplified version of Langevin dynamics. The present result is a first step to understand the time-dependent behavior of colloids near interfaces. This study shows a wide range of adsorption trajectories, where the emphasis lies on a strong dependence of the dynamics on the orientation of the colloid at initial contact with the interface. We believe that the observed orientational dependence in our simple model can be recovered in suitable experimental systems., Comment: 14 pages, 8 figures, 1 table

Using simulations of hard rods in smectic-A states, we find non-gaussian diffusion and heterogeneous dynamics due to the equilibrium periodic smectic density profiles, which give rise to permanent barriers for layer-to-layer diffusion. This relaxation behavior is surprisingly similar to that of non-equilibrium supercooled liquids, although there the particles are trapped in transient (instead of permanent) cages. Interestingly, we also find stringlike clusters of up to 10 inter-layer rods exhibiting dynamic cooperativity in this equilibrium state., Comment: 10 pages, 4 figures

Nonlinear ionic screening theory for heterogeneously charged spheres is developed in terms of a mode-decomposition of the surface charge. A far-field analysis of the resulting electrostatic potential leads to a natural generalization of charge renormalization from purely monopolar to dipolar, quadropolar, etc., including mode-couplings. Our novel scheme is generally applicable to large classes of surface heterogeneities, and is explicitly applied here to Janus spheres with differently charged upper and lower hemispheres, revealing strong renormalization effects for all multipoles., Comment: 2 figures

We present a study of the isotropic-nematic interface in a system of hard spherocylinders. First we compare results from Monte Carlo simulations and Onsager density functional theory for the interfacial profiles of the orientational order parameter and the density. Those interfacial properties that are not affected by capillary waves are in good agreement, despite the fact that Onsager theory overestimates the coexistence densities. Then we show results of a Monte Carlo study of the capillary waves of the interface. In agreement with recent theoretical investigations (Eur.Phys.J. E {\bf 18} 407 (2005)) we find a strongly anistropic capillary wave spectrum. For the wave-numbers accessed in our simulations, the spectrum is quadratic, i.e.elasticity does not play a role. We conjecture that this effect is due to the strong bending rigidity of the director field in suspensions of spherocylinders., Comment: 8 pages, 10 figures

We theoretically study the nematic ordering transition of rods that are able to elastically adjust their mutually excluded volumes. The model rods, which consist of a hard core surrounded by a deformable shell, mimic the structure of polymer-coated, rod-like fd virus particles that have recently been the object of experimental study [K.~Purdy et al., Phys. Rev. Lett. \textbf{94}, 057801 (2005)]. We find that fluids of such soft rods exhibit an isotropic-nematic phase transition at a density higher than that of the corresponding hard-rod system of identical diameter, and that at coexistence the order parameter of the nematic phase depends non monotonically on the elastic properties of the polymer coating. For binary mixtures of hard and soft rods, the topology of the phase diagram turns out to depend sensitively on the elasticity of shell. The lower nematic-nematic critical point, discovered in mixtures of bare and polymer-coated fd virus particles, is not reproduced by the theory., Comment: 8 pages, 5 figures

Within an Onsager-like density functional theory we explore the thermodynamic and structural properties of an isotropic and nematic fluid of hard needle-like colloids in contact with a hard substrate coated with a soft short-ranged attractive or repulsive layer. As a function of the range and the strength of the soft interactions we find wetting and drying transitions, a pre-drying line, and a symmetry-breaking transition from uniaxial to biaxial in the wetting and drying film., Comment: 7 pages, 2 figures

We study sedimentation equilibrium of low-salt suspensions of binary mixtures of charged colloids, both by Monte Carlo simulations of an effective colloids-only system and by Poisson-Boltzmann theory of a colloid-ion mixture. We show that the theoretically predicted lifting and layering effect, which involves the entropy of the screening ions and a spontaneous macroscopic electric field [J. Zwanikken and R. van Roij, Europhys. Lett. {\bf 71}, 480 (2005)], can also be understood on the basis of an effective colloid-only system with pairwise screened-Coulomb interactions. We consider, by theory and by simulation, both repelling like-charged colloids and attracting oppositely charged colloids, and we find a re-entrant lifting and layering phenomenon when the charge ratio of the colloids varies from large positive through zero to large negative values.

We have measured equilibrium sedimentation profiles in a colloidal model system with confocal microscopy. By tuning the interactions, we have determined the gravitational length in the limit of hard-sphere-like interactions, and using the same particles, tested a recent theory [R.van Roij, J. Phys. Cond. Mat. 15, S3569, (2003)], which predicts a significantly extended sedimentation profile in the case of charged colloids with long-ranged repulsions, due to a spontaneously formed macroscopic electric field. For the hard-sphere-like system we find that the gravitational length matches that expected. By tuning the buoyancy of the colloidal particles we have shown that a mean field hydrostatic equilibrium description even appears to hold in the case that the colloid volume fraction changes significantly on the length scale of the particle size. The extended sedimentation profiles of the colloids with long-ranged repulsions are well-described by theory. Surprisingly, the theory even seems to hold at concentrations where interactions between the colloids, which are not modeled explicitly, play a considerable role.

Within a second virial theory, we study bulk phase diagrams as well as the free planar isotropic-nematic interface of binary mixtures of nonadditive thin and thick hard rods. For species of the same type the excluded volume is determined only by the dimensions of the particles, whereas for dissimilar ones it is taken to be larger or smaller than that, giving rise to a nonadditivity that can be positive or negative. We argue that such a nonadditivity can result from modelling of soft interactions as effective hard-core interactions. The nonadditivity enhances or reduces the fractionation at isotropic-nematic ($IN$) coexistence and may induce or suppress a demixing of the high-density nematic phase into two nematic phases of different composition ($N_1$ and $N_2$), depending on whether the nonadditivity is positive or negative. The interfacial tension between co-existing isotropic and nematic phases show an increase with increasing fractionation at the $IN$ interface, and complete wetting of the $IN_2$ interface by the $N_1$ phase upon approach of the triple point coexistence. In all explored cases bulk and interfacial properties of the nonadditive mixtures exhibit a striking and quite unexpected similarity with the properties of additive mixtures of different diameter ratio., Comment: 12 pages, revised version, submitted to JCP

Quantitative predictions are presented of a depletion-induced torque and force acting on a single colloidal hard rod immersed in a solvent of hard spheres close to a planar hard wall. This torque and force, which are entirely of entropic origin, may play an important role for the key-lock principle, where a biological macromolecule (the key) is only functional in a particular orientation with respect to a cavity (the lock).